We thank Amy Chen for her help in preparing for the video shoot and for demonstrating the technique. This work conducted by ENIGMA: Ecosystems and Networks Integrated with Genes and Molecular Assemblies (http://enigma.lbl.gov), a Scientific Focus Area Program at Lawrence Berkeley National Laboratory, was supported by the Office of Science, Office of Biological and Environmental Research, of the U. S. Department of Energy under Contract No. DE-AC02-05CH11231.

Note: The protocol below is tailored for determination of gene targets of the RR DVU3023 from the bacterium Desulfovibrio vulgaris Hildenborough. It can be adapted to any other transcriptional regulator of interest.
1. Clone and Purify RR
<ol>
	<li>Clone the RR gene, specifically DVU3023, from D. vulgaris Hildenborough into an Escherichia coli expression vector such that the gene is C-terminally His-tagged and expression is under the control of a T7 promoter. 
	Note: Several cloning methods may be used and is determined by the researcher. Alternate affinity tags may also be used.</li>
	<li>Transform the expression construct into E. coli BL21 (DE3) expression strain.</li>
	<li>Grow up 1 L of the BL21 expression strain at 37 &#176;C. At mid-log phase, add IPTG to 0.5 mM to induce protein expression. Continue growth at room temperature for 24 hr.</li>
	<li>Centrifuge cells at 5,000 x g for 10 min at 4 &#176;C. Resuspend cells in lysis buffer (20 mM sodium phosphate, pH 7.4, 500 mM NaCl, 40 mM imidazole, 1 mg/ml lysozyme, 1x benzonase nuclease). Lyse cells using a French press at 4 &#176;C.</li>
	<li>Centrifuge lysate at 15,000 x g for 30 min at 4 &#176;C. Filter using a 0.45 &#181;m syringe filter.</li>
	<li>Wash a 1 ml Ni Sepharose column on an FPLC instrument using 10 ml of wash buffer (20 mM sodium phosphate, pH 7.4, 500 mM NaCl, 40 mM imidazole).</li>
	<li>Load lysate on the column. Wash column with 20 ml of wash buffer.</li>
	<li>Elute DVU3023 using a gradient of 0-100% elution buffer (20 mM sodium phosphate, pH 7.4, 500 mM NaCl, 500 mM imidazole).</li>
	<li>Load eluted fractions on a desalting column, and wash with a desalting buffer (20 mM sodium phosphate pH 7.4, 100 mM NaCl).</li>
	<li>Centrifuge protein sample in a high molecular weight cutoff centrifuge filter to concentrate the protein. Add DTT to 0.1 mM and glycerol to 50% and store protein at -20 &#176;C. 
	Note: Purification methods will need to be optimized individually for each protein under study.</li>
</ol>
2. Determine Gene Target for RR Using Electrophoretic Mobility Shift Assay (EMSA)
<ol>
	<li>PCR amplify 400 bp region upstream of the candidate target gene DVU3025, using D. vulgaris Hildenborough genomic DNA as the template, and a forward unlabeled primer and a reverse 5&#8217;-biotin-labeled primer. 
	Note: Tips to select a candidate target gene: Often RRs bind the upstream regions of their own gene/operon, or they may regulate proximally encoded genes. If the RR has orthologs in other species, look for genes that are conserved in the neighborhood. Alternate methods include choosing candidate genes based on regulon predictions (e.g.,&#160;RegPrecise5), or predicting sigma54-dependent promoters6 for RRs that are themselves sigma54-dependent.</li>
	<li>Run the PCR product on a 1% agarose gel, cut out the 400 bp sized product, and purify the DNA using a gel extraction clean up kit.</li>
	<li>Mix DVU3023 protein (0.5 pmol) with 100 fmol of biotinylated DNA substrate in 10 mM Tris HCl, pH 7.5, 50 mM KCl, 5 mM MgCl2, 1 mM DTT, 25% glycerol and 1 &#181;g/ml poly dI.dC (non-specific competitor DNA) in a total volume of 20 &#181;l. Also set up a reaction without any protein as a control. Incubate the reactions at room temperature on the bench for 20 min.</li>
	<li>Note: These are standard conditions listed and other reaction components may be added depending on the RR of interest. If activation is desired, add 50 mM acetyl phosphate to the reaction (often RRs will bind DNA in vitro without activation).</li>
	<li>Pre-run a precast mini 6% polyacrylamide-0.5x TBE gel in 0.5x TBE buffer at 100 V for 30 min. Add 5 &#181;l of 5x loading buffer (0.1% bromophenol blue, 0.1% xylene cyanol, 30% glycerol in 1x TBE) to the binding reactions and load 18 &#181;l of the reactions on to the gel. Run at 100 V for 2 hr. 
	Note: In order to avoid overheating, which can lead to disassociation of the protein-DNA complexes, fill the outer buffer chamber in the electrophoresis system with running buffer such that the majority of the gel is insulated by buffer. Alternately, the gel may be run at 4 &#176;C.</li>
	<li>Cut a charged nylon membrane to the size of the gel and soak it in 0.5x TBE for at least ten minutes. Remove the gel from the cassette. Cut any ridges off the gel and sandwich the gel and membrane between two thick filter papers soaked in 0.5x TBE, and place inside a semi-dry blotting apparatus and run at 20 V for 30 min.</li>
	<li>Place the membrane inside a commercial UV light crosslinker instrument and set the time to 3 min. 
	Note: The membrane may now be stored dry at room temperature for several days before proceeding to the next step.</li>
	<li>Use a commercially available chemiluminescent detection kit that employs a streptavidin-horse radish peroxidase conjugate to develop the blot as per manufacturer&#8217;s instructions.</li>
	<li>Image the blot using a computer hooked up to a CCD equipped camera and look for a shift in the DNA substrate mobility in the presence of the RR, which indicates that the RR binds the DNA being tested. 
	Note: The specificity of the protein-DNA shift can be tested by including in the reaction an excess of unlabeled competitor DNA, which should eliminate or decrease the shift seen.</li>
</ol>
3. Verify Target Enrichment after Genomic DNA-protein Binding
<ol>
	<li>Genomic DNA-protein binding reaction
	<ol>
		<li>Shear genomic DNA to an average size of 500 bp by sonicating 100 &#181;l genomic DNA (at concentrations 100-200 ng/&#181;l) in a 1.5 ml microfuge tube using a microtip with 9 low amplitude pulses of 1 sec, with 2 sec gap in between each. 
		Note: If the DNA splatters to the sides of the tube, spin down the contents after every 3 pulses. The exact conditions used will vary with the sonicator instrument used. The sheared DNA fragments may range from 100-1000 bp, but the average size should be within 400-600 bp. Too much shearing may result in fewer intact binding sites, and either too much or too little shearing will affect library preparation during the downstream whole genome amplification steps.</li>
		<li>Mix 2-3 &#181;g of sheared genomic DNA with RR protein (0.5 pmol of DVU3023) in 10 mM Tris-HCl pH 7.5, 1 mM DTT, 50 mM KCl, 5 mM MgCl2, 25% glycerol, and 50 mM acetyl phosphate. Incubate reactions at 25 &#176;C for 30 min. Transfer 10 &#181;l of this reaction to a 1.5 ml tube and label as input DNA. 
		Note: Acetyl phosphate is added to activate the protein by in vitro phosphorylation. Phosphorylation stimulates DNA binding, but many RRs also bind DNA in vitro without activation7. The amount of protein added to the reaction will depend on activity of the protein prep. The EMSA may be used to optimize this amount.</li>
	</ol>
	</li>
	<li>Affinity purify the protein-bound DNA
	<ol>
		<li>Add 30 &#181;l of Ni-NTA agarose resin to a 0.6 ml microfuge tube. Centrifuge at 100 x g for 1 min to collect the resin at the bottom, and remove the supernatant. Add 100 &#181;l of wash buffer (10 mM Tris-HCl, pH 7.5, 5 mM MgCl2, 50 mM KCl, 25% glycerol), flick the tube to mix, and centrifuge at 100 x g for 2 min. Remove supernatant.</li>
		<li>Add the remaining 90 &#181;l of the binding reaction to the washed Ni-NTA resin and incubate in a rotary shaker for 30 min.</li>
		<li>Centrifuge at 100 x g for 2 min, and remove the supernatant (unbound DNA). Add 100 &#181;l of wash buffer, flick the tube to mix contents, and centrifuge at 100 x g for 2 min. Remove supernatant. Repeat the wash step twice more.</li>
		<li>Add 35 &#181;l of elution buffer (20 mM sodium phosphate buffer pH 8, 500 mM NaCl, 500 mM imidazole) to the resin and mix by vortexing. Incubate on the bench at RT for 5 min. Centrifuge at 100 x g for 2 min. Transfer the supernatant to a new 1.5 ml tube and label as the protein-bound DNA fraction.</li>
		<li>Also add 35 &#181;l of the elution buffer to the input DNA.</li>
		<li>Purify the input and the protein-bound DNA fractions using a PCR purification kit.</li>
	</ol>
	</li>
	<li>Whole genome amplify the input and protein-bound DNA samples.
	<ol>
		<li>Add 10 &#181;l of the input and protein-bound DNA to separate PCR tubes. Add 1 &#181;l of 10x Fragmentation buffer, 2 &#181;l of Library Preparation buffer, and 1 &#181;l of Library Stabilization solution to each sample. Mix well by vortexing. Heat in a thermocycler at 95 &#176;C for 2 min, and chill on ice. 
		Note: The starting amount of DNA should be at least 10 ng in order to avoid introducing amplification bias.</li>
		<li>Add 1 &#181;l of Library Preparation Enzyme, mix by pipetting, and incubate in thermal cycler at 16 &#176;C/20 min, 24 &#176;C/20 min, 37 &#176;C for 20 min, 75 &#176;C/5 min, and hold at 4 &#176;C.</li>
		<li>To each tube, add 47.5 &#181;l of water, 7.5 &#181;l of 10x amplification master mix, and 5 &#181;l of polymerase. Mix well and heat at 95 &#176;C/3 min, followed by 20 cycles of 94 &#176;C/15 sec, 65 &#176;C/5 min. Hold reactions at 4 &#176;C.</li>
		<li>Purify the amplified DNA samples using a PCR purification kit and measure DNA concentration spectrophotometrically.</li>
	</ol>
	</li>
	<li>Verify target enrichment in the protein-bound DNA using qPCR
	<ol>
		<li>Design qPCR primers to amplify 200 bp of the upstream region of the EMSA-verified target gene (DVU3025) using any freely available primer design software.</li>
		<li>Set up triplicate qPCR reactions for each DNA template with each primer set. Prepare a master mix for each primer set. 1x master mix contains 10 &#181;l of 2x Sybr Green qPCR mix, 0.5 &#181;M of each primer and water to a total of 18 &#181;l. Aliquot 18 &#181;l of the master mix per well of a 96-well PCR plate.</li>
		<li>Dilute the amplified and purified input and protein-bound DNA samples to 5 ng/&#181;l with water and use as DNA template. Add 2 &#181;l of DNA template to the wells.</li>
		<li>Seal the plate with an ultra clear qPCR sealing film, spin down plate in a centrifuge at 200 x g/1 min. Place plate in a real time qPCR machine. Cycle as follows using the associated qPCR software: 95 &#176;C/1 min, and 40 cycles of 95 &#176;C/10 sec, 59 &#176;C/15 sec, and 70 &#176;C/35 sec.</li>
		<li>Also set up triplicate qPCR reactions with primers to amplify upstream regions of an unrelated gene (negative control).</li>
		<li>Calculate &#916;CT by subtracting CT value of protein-bound DNA from that of the input DNA. Calculate fold enrichment of target gene in the protein-bound DNA as 2&#916;CT. 
		Note: If the target upstream region was enriched in the protein-bound DNA vs. the input DNA, it indicates that the binding reactions and affinity purification were successful. Proceed to step 4. If enrichment of the target upstream region was not observed, repeat the binding reactions (step 3.1) with different amount of protein.</li>
	</ol>
	</li>
</ol>
4. DNA Labeling and Array Hybridization
<ol>
	<li>Label Input DNA with Cy3 and Enriched DNA with Cy5
	 ote: Cy3 and Cy5 dyes are light sensitive and care should be taken to keep light exposure to a minimum.
	<ol>
		<li>Mix 1 &#181;g DNA with 40 &#181;l Cy3/Cy5-labeled 9-mers and adjust volume to 80 &#181;l with water.</li>
		<li>Heat denature at 98 &#176;C for 10 min in the dark (in a thermocycler). Quick chill on ice for 2 min.</li>
		<li>Add 2 &#181;l Klenow polymerase (3&#8217;-5&#8217;-exo-, 50,000 U/ml), 5 mM dNTPs, and 8 &#181;l water to each reaction, mix well, and incubate at 37 &#176;C for 2 hr in the dark (in a thermocycler).</li>
		<li>Add EDTA to 50 mM to stop the reaction and NaCl solution to 0.5 M.</li>
		<li>Transfer samples to 1.5 ml tubes containing 0.9 volume of isopropanol, incubate in dark for 10 min, and centrifuge at 12,000 x g for 10 min. The pellet should be pink for Cy3-labeled DNA and blue for Cy5-labeled DNA.</li>
		<li>Wash pellet with 80% ethanol (500 &#181;l) and centrifuge at 12,000 x g for 2 min. Air-dry pellets for 5-10 min in dark. 
		Note: Pellets may be stored at -20 &#176;C.</li>
		<li>Resuspend pellet in 25 &#181;l water. Measure DNA concentration spectrophotometrically.</li>
		<li>Pool together 6 &#181;g each of the Cy3- and Cy5-labeled DNA in a 1.5 ml tube and vacuum dry in a centrifuge on low heat in the dark (cover the centrifuge lid if it is transparent). 
		Note: Pellets may be stored at -20 &#176;C until ready for hybridization.</li>
	</ol>
	</li>
	<li>Microarray Hybridization
	<ol>
		<li>Turn the hybridization system on 3-4 hr prior to use and set the temperature to 42 &#176;C.</li>
		<li>Prepare a hybridization solution master mix such that 1x solution contains 11.8 &#181;l 2x Hybridization buffer, 4.7 &#181;l Hybridization Component A, and 0.5 &#181;l alignment oligo.</li>
		<li>Resuspend the pellets in 5 &#181;l water. Add 13 &#181;l of this mix to the sample. Vortex for 15 sec, incubate at 95 &#176;C in a dry bath for 5 min. Keep samples at 42 &#176;C in the hybridization system until ready for loading.</li>
		<li>Prepare the microarray slide-mixer assembly, according to manufacturer&#8217;s protocol.</li>
		<li>Place the mixer-slide assembly within the hybridization system. Load 16 &#181;l of sample into the fill port, seal the ports with an adhesive film, turn mixing on in the system, and hybridize for 16-20 hr at 42 &#176;C.</li>
		<li>Prepare 1x wash buffers I (250 ml), II (50 ml) and III (50 ml) by diluting the 10x commercially available buffers in water, and to each add DTT to 1 mM. Warm Buffer I to 42 &#176;C.</li>
		<li>Slide the mixer-slide into a disassembly tool, and place inside a dish containing warm Buffer I. Peel the mixer off, while vigorously shaking the disassembly tool by hand.</li>
		<li>Place the slide into a container with 50 ml of Wash buffer I and shake vigorously for 2 min by hand.</li>
		<li>Transfer slide into a second container with 50 ml Wash buffer II and shake vigorously for 1 min by hand.</li>
		<li>Transfer slide into a third container with 50 ml Wash buffer III, and shake vigorously for 15 sec by hand.</li>
		<li>Quickly blot the edges of the slide on a paper towel, and place in a slide rack. Centrifuge at 200 x g for 2 min to dry the slide. Place within slide case, wrap it with foil, and store in a desiccator.</li>
	</ol>
	</li>
	<li>Scanning the Array
	<ol>
		<li>Place the slide within the microarray scanner according to the instrument&#8217;s instructions.</li>
		<li>Use a scanner software to set the wavelengths as 532 nm = Cy3 and 635 nm = Cy5, and the initial photomultiplier gains between 350 and 400.</li>
		<li>Preview the slide to locate the array on the slide. Select the array region for scanning.</li>
		<li>Scan the array and adjust the photomultiplier settings such that the array features are mostly yellow. The histogram should show the red and green curves superimposed or as close to each other as possible. The curves should end above the 10-5 normalized counts at saturation.</li>
		<li>Save both 532 and 635 nm images separately.</li>
	</ol>
	</li>
	<li>Array Data Analysis
	<ol>
		<li>Import the images into an array analysis software. Using the software, create two pair reports (.pair) for each array (for Cy3 and Cy5 images). Create scaled log2ratio files using the pair files. Use log2ratio files to search for peaks using a sliding window of 500 bp. Map the peak loci to the upstream regions of gene targets.</li>
		<li>Check if the positive target (determined using EMSA and confirmed by qPCR, DVU3025 in this example) appears within the top peaks. 
		Note: If the positive target is not among the top hits, it is possible that the RR under consideration has several gene targets and replicate DAP-chips can be conducted and peaks common to all replicates can be used to generate a target list.</li>
	</ol>
	</li>
</ol>
5. Binding Site Motif Prediction and Validation
<ol>
	<li>Retrieve sequences for upstream regions (400 bp) for the top gene targets as generated by DAP-chip analysis. Apply MEME on these sequences through Microbes Online (meme.nbcr.net) to predict motifs. 
	Note: Enhancer binding proteins such as sigma54-dependent regulators can bind several hundred base pairs upstream of the start site. For other transcriptional regulators, a shorter region such as 200 bp upstream will be sufficient.</li>
	<li>Design top and bottom strand DNA oligomers that contain the predicted binding site motif flanked by 10 bases on either end. Order the top strand 5&#8217; biotinylated.</li>
	<li>Mix the top and bottom strand oligos in 1:1.5 ratio in 10 mM Tris HCl, pH 8.0, 1 mM EDTA, and 50 mM NaCl in a total reaction volume of 20 &#181;l in PCR tubes. Heat to 95 &#176;C for 5 min in a thermal cycler, followed by slow cooling to 25 &#176;C. Dilute ten fold and use as wild type dsDNA substrate for EMSA.</li>
	<li>Prepare modified substrates using the same steps as above, but design oligomers to carry 4-6 substitutions in the conserved bases of the binding motif.</li>
	<li>Set up binding reactions with the RR and the wild type or modified substrates and examine using EMSA as described in step 2.</li>
	<li>Note: The predicted motif is validated if the RR binds the wild type substrate but not the modified one.</li>
</ol>
6. Conservation of Motif in Other Related Bacterial Species
<ol>
	<li>Select genomes of interest that contain orthologs of DVU3023. Obtain sequence and annotation files for the genomes from the NCBI site.</li>
	<li>Use a programming language such as Perl to write scripts that will use the motif sequences from the DAP-chip targets to build a position weighted matrix and use the matrix to score similar motifs present in other sequenced genomes.</li>
</ol>

<ol>
	<li>Galperin, M. Y. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Diversity+of+structure+and+function+of+response+regulator+output+domains.">Diversity of structure and function of response regulator output domains.</a> Curr Opin Microbiol. 13, 150-159 (2010).</li><li>Casino, P., Rubio, V., Marina, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+mechanism+of+signal+transduction+by+two-component+systems.">The mechanism of signal transduction by two-component systems.</a> Curr Opin Struct Biol. 20, 763-771 (2010).</li><li>Aparicio, O., Geisberg, J. V., Struhl, K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Chromatin+immunoprecipitation+for+determining+the+association+of+proteins+with+specific+genomic+sequences+in+vivo.">Chromatin immunoprecipitation for determining the association of proteins with specific genomic sequences in vivo.</a> Curr Protoc Cell Biol. 17 (17), (2004).</li><li>McCleary, W. R., Stock, J. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Acetyl+phosphate+and+the+activation+of+two-component+response+regulators.">Acetyl phosphate and the activation of two-component response regulators.</a> J Biol Chem. 269, 31567-31572 (1994).</li><li>Novichkov, P. S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=RegPrecise+web+services+interface:+programmatic+access+to+the+transcriptional+regulatory+interactions+in+bacteria+reconstructed+by+comparative+genomics.">RegPrecise web services interface: programmatic access to the transcriptional regulatory interactions in bacteria reconstructed by comparative genomics.</a> Nucleic Acids Res. 40, 604-608 (2012).</li><li>Studholme, D. J., Buck, M., Nixon, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Identification+of+potential+sigma(N)-dependent+promoters+in+bacterial+genomes.">Identification of potential sigma(N)-dependent promoters in bacterial genomes.</a> Microbiology. 146 (12), 3021-3023 (2000).</li><li>Rajeev, L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Systematic+mapping+of+two+component+response+regulators+to+gene+targets+in+a+model+sulfate+reducing+bacterium.">Systematic mapping of two component response regulators to gene targets in a model sulfate reducing bacterium.</a> Genome Biol. 12, 99 (2011).</li><li>Crooks, G. E., Hon, G., Chandonia, J. M., Brenner, S. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=WebLogo:+a+sequence+logo+generator.">WebLogo: a sequence logo generator.</a> Genome Res. 14, 1188-1190 (2004).</li><li>Barbieri, C. M., Mack, T. R., Robinson, V. L., Miller, M. T., Stock, A. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Regulation+of+response+regulator+autophosphorylation+through+interdomain+contacts.">Regulation of response regulator autophosphorylation through interdomain contacts.</a> J Biol Chem. 285, 32325-32335 (2010).</li><li>Bourret, R. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Receiver+domain+structure+and+function+in+response+regulator+proteins.">Receiver domain structure and function in response regulator proteins.</a> Curr Opin Microbiol. 13, 142-149 (2010).</li><li>Gao, R., Stock, A. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Molecular+strategies+for+phosphorylation-mediated+regulation+of+response+regulator+activity.">Molecular strategies for phosphorylation-mediated regulation of response regulator activity.</a> Curr Opin Microbiol. 13, 160-167 (2010).</li><li>Hung, D. C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Oligomerization+of+the+response+regulator+ComE+from+Streptococcus+mutans+is+affected+by+phosphorylation.">Oligomerization of the response regulator ComE from Streptococcus mutans is affected by phosphorylation.</a> J Bacteriol. 194, 1127-1135 (2012).</li><li>Ladds, J. C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+response+regulator+Spo0A+from+Bacillus+subtilis+is+efficiently+phosphorylated+in+Escherichia+coli.">The response regulator Spo0A from Bacillus subtilis is efficiently phosphorylated in Escherichia coli.</a> FEMS Microbiol Lett. 223, 153-157 (2003).</li><li>Wen, Y., Feng, J., Scott, D. R., Marcus, E. A., Sachs, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Involvement+of+the+HP0165-HP0166+two-component+system+in+expression+of+some+acidic-pH-upregulated+genes+of+Helicobacter+pylori.">Involvement of the HP0165-HP0166 two-component system in expression of some acidic-pH-upregulated genes of Helicobacter pylori.</a> J Bacteriol. 188, 1750-1761 (2006).</li><li>Liu, X., Noll, D. M., Lieb, J. D., Clarke, N. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=DIP-chip:+rapid+and+accurate+determination+of+DNA-binding+specificity.">DIP-chip: rapid and accurate determination of DNA-binding specificity.</a> Genome Res. 15, 421-427 (2005).</li><li>Gossett, A. J., Lieb, J. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=DNA+Immunoprecipitation+(DIP)+for+the+Determination+of+DNA-Binding+Specificity.">DNA Immunoprecipitation (DIP) for the Determination of DNA-Binding Specificity.</a> CSH Protoc. 2008, (2008).</li><li>Mascher, T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+Streptococcus+pneumoniae+cia+regulon:+CiaR+target+sites+and+transcription+profile+analysis.">The Streptococcus pneumoniae cia regulon: CiaR target sites and transcription profile analysis.</a> J Bacteriol. 185, 60-70 (2003).</li><li>Park, P. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=ChIP-seq:+advantages+and+challenges+of+a+maturing+technology.">ChIP-seq: advantages and challenges of a maturing technology.</a> Nat Rev Genet. 10, 669-680 (2009).</li><li>Lefrancois, P., Zheng, W., Snyder, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=ChIP-Seq+using+high-throughput+DNA+sequencing+for+genome-wide+identification+of+transcription+factor+binding+sites.">ChIP-Seq using high-throughput DNA sequencing for genome-wide identification of transcription factor binding sites.</a> Methods Enzymol. 470, 77-104 (2010).</li><li>Ho, J. W., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=ChIP-chip+versus+ChIP-seq:+lessons+for+experimental+design+and+data+analysis.">ChIP-chip versus ChIP-seq: lessons for experimental design and data analysis.</a> BMC Genomics. 12, 134 (2011).</li><li>Liu, X., Lee, C. K., Granek, J. A., Clarke, N. D., Lieb, J. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Whole-genome+comparison+of+Leu3+binding+in+vitro+and+in+vivo+reveals+the+importance+of+nucleosome+occupancy+in+target+site+selection.">Whole-genome comparison of Leu3 binding in vitro and in vivo reveals the importance of nucleosome occupancy in target site selection.</a> Genome Res. 16, 1517-1528 (2006).</li></ol>

The authors have no conflict of interest to disclose.

The DAP-chip method described here was successfully used to determine the gene targets for several RRs in Desulfovibrio vulgaris Hildenborough7 of which one is shown here as a representative result. For RR DVU3023, choosing a candidate gene target was straightforward. DVU3025 is located immediately downstream of the RR gene, and the RR and target genes are conserved in several Desulfovibrio species, and additionally DVU3025 has a predicted sigma54-dependent promoter. The EMSA provides a simple method to rapidly test the RR for binding to the candidate target gene, and also allows the assessment of the activity of the purified protein sample, as well as determine optimal protein-DNA binding conditions.
The DNA binding activity can also be tested in the presence and absence of acetyl phosphate to see if phosphorylation is necessary for DNA binding. Not all RRs are phosphorylated by small molecule donors9, in these cases the purified cognate sensor kinase, if available, may be used to activate the protein. There are also atypical RRs known that lack key active site residues in the receiver domain and are not activated by phosphorylation10. For the majority of RRs studied, phosphorylation stimulates DNA binding11. However there are examples where phosphorylation does not affect in vitro DNA binding12, and there are examples where phosphorylation is required for binding13, and also cases where the subset of promoters bound increases with phosphorylation14. The DAP-chip method may be performed with and without phosphorylation to determine if there are differences in the set of promoters bound by the RR. However, it should also be noted that some RRs may be purified in a functionally phosphorylated form from E. coli13.
The protocol described here is for a His-tagged protein such that the protein-bound DNA is affinity purified using Ni-NTA agarose resin, but it can easily be adapted for any kind of tagged protein by using the appropriate affinity resin for pull-down. Following the DAP-WGA steps, the qPCR step provides an additional control to ensure that the protein-gDNA binding conditions were appropriate and that the protein-bound DNA was successfully enriched by the affinity purification. Greater than 3-fold enrichment is sufficient to proceed with the hybridization step. Unlike the EMSA reaction, there is no non-specific competitor DNA like poly-dI.dC added to the gDNA binding reaction since it interferes with the WGA step. Due to this if the protein sample has non-specific DNA binding activity then clear enrichment of the confirmed target DNA will not be observed by qPCR. Thus the qPCR control step can be used to optimize the amount of protein used in the gDNA binding reaction. Commercially available kits for whole genome amplification claim to introduce no amplification bias when following their guidelines for minimum starting DNA material and low number of amplification cycles. The amplification bias may be checked using qPCR with various primer sets on amplified versus unamplified DNA. Any effect on downstream array hybridization may also be determined by differentially labeling and hybridizing pooled amplified and unamplified genomic DNA.
The list of peaks that are generated following array data analysis is usually long with several peaks having a false discovery rate value of 0. Therefore a combination of other peak attributes such as high log2R scores and cutoff_p values (as generated by the array analysis software used in this study) are used to cull the list down to the highest confidence gene targets. The presence of the pre-determined gene target among the top five hits greatly strengthens the confidence in the data set. Performing this assay for a number of regulatory proteins from the same organism will also help to identify &#8220;sticky&#8221; DNA sequences that appear in several data sets7. Additionally if a binding site motif can be predicted and validated then the presence of the motif in peaks lower down in the list may also be used to choose a conservative target gene list. Enrichment of DNA sequences other than promoter regions may be artifacts of array hybridization or may indicate the regulation of previously unidentified open reading frames or small RNAs7. For regulators where a gene target could not be predetermined using EMSA, the DAP-chip assay may be performed &#8220;blind&#8221; 7. In such cases also, identification of a binding site motif will improve the selection of gene targets. Determining the protein concentration to be used in the assay will depend on the non-specific binding activity of the protein, which may be assessed by EMSA using randomly chosen DNA substrates. Lower concentrations work better for those proteins with high nonspecific binding activity. The reliability of a blind DAP-chip may be improved by performing replicate assays with different protein concentrations. For RRs with few targets, just two or even a single assay may be sufficient to generate a target list which may be validated using subsequent EMSAs. The DAP-chip data for such RRs usually show a clear jump in the log2R scores or cutoff_p values beyond the initial few peaks. For RRs with several targets, data from three replicates may be analyzed to generate a list of common peaks, some of which may be selected for EMSA validation.
The ability to predict binding site motifs based on the DAP-chip targets adds to the value of this method and vastly increases the information gained. EMSA again provides an efficient protocol to experimentally verify the predictions. Software scripts in Perl or other programming language may be used to rapidly search through available genome sequences. The results will identify both orthologous target genes as well as target genes uniquely regulated in the genomes searched. In the representative result shown here, three of the four upstream target regions identified for DVU3023 have annotated gene functions related to lactate uptake and oxidation. Additionally since the binding site motif for DVU3023 was validated, other genomes could be searched for similar motif sites. Together the results indicate that the DVU3022-3023 TCS is well conserved in the Desulfovibrio and other related sulfate-reducing bacteria, and that it regulates genes for lactate transport and oxidation to acetate. Since lactate is the primary carbon and electron source used by these organisms, DVU3023 likely plays a key role in their physiology. Among the top DAP-chip peaks for RR DVU3023, there was also an intergenic region. Although a binding site motif was not found in this region, the presence of a predicted sigma-54-dependent promoter suggests that there may be an unidentified orf or sRNA encoded, and that it could be a true target for DVU3023. This finding highlights the value of the experimental DAP-chip approach combined with binding site predictions as opposed to determining target sites based on computational predictions alone.
Similar in vitro methods such as the one described here have only been used in a few cases15-17 and very rarely for a previously unstudied regulator. The optimization EMSA and qPCR steps prior to the array hybridization will substantially aid in analyzing the results when the method is to be used for a novel regulator. If array printing becomes an impediment, the DAP steps may be combined with next generation sequencing to obtain binding sites18,19. As more array-based methods become substituted with sequencing based strategies, such future adaptation of the method circumvents the need to design custom tiling microarrays for the organism of interest. ChIP-Seq technologies result in greater sensitivity and specificity of detection of peaks when compared to ChIP-chip methods, and are also rapidly becoming more cost-effective20. Although this article focused on two component response regulators, this method can be used for any prokaryotic or eukaryotic transcription factor15. Nucleosome occupancies often have an effect on targets that are regulated and comparing the target binding sites for eukaryotic transcription factors by in vivo and in vitro analyses can reveal these effects21.

The ability of bacteria to survive and thrive is critically dependent on how well they are able to perceive and respond to perturbations in their environments, and this in turn is dependent on their signal transduction systems. The number of signaling systems a bacterium encodes has been called its &#8220;microbial IQ&#8221; and can be an indication of both variability of its environment and its ability to sense multiple signals and fine tune its response1. Two component signal transduction systems (TCS) are the most prevalent signaling systems used by bacteria, and they consist of a histidine kinase (HK) that senses the external signal and transmits via phosphorylation to an effector response regulator (RR)2. RRs can have a variety of output domains and thus different effector modes, but the most common response is transcriptional regulation via a DNA binding domain1. The signals sensed and the corresponding functions of the vast majority of TCSs remain unknown.
Although in vivo methods such as ChIP-chip are routinely used for determination of genomic binding sites of transcription factors3, they can only be used for bacterial two component system RRs if the activating conditions or signals are known. Often the environmental cues that activate a TCS are harder to determine than their gene targets. The in vitro microarray based assay described here can be used to effectively and rapidly determine the gene targets and predict functions of TCSs. This assay takes advantage of the fact that RRs can be phosphorylated and thus activated in vitro using small molecule donors like acetyl phosphate4.
In this method, named DAP-chip for DNA-affinity-purified-chip (Figure 1), the RR gene of interest is cloned with a His-tag in E. coli, and a subsequently purified tagged protein is allowed to bind to sheared genomic DNA. The protein-bound DNA is then enriched by affinity-purification, the enriched and input DNA are amplified, fluorescently labeled, pooled together and hybridized to a tiling array that is custom made to the organism of interest (Figure 1). Microarray experiments are subject to artifacts and therefore additional steps are employed to optimize the assay. One such step is to attempt to determine one target for the RR under study using electrophoretic mobility shift assays (EMSA) (see workflow in Figure 2). Then, following binding to genomic DNA and the DAP steps, the protein-bound and input DNA are examined by qPCR to see if the positive target is enriched in the protein-bound fraction relative to the input fraction, thus confirming optimal binding conditions for the RR (Figure 2). After array hybridization, the data are analyzed to find peaks of higher intensity signal indicating genomic loci where the protein had bound. Functions may be predicted for the RR based on the gene targets obtained. The target genomic loci are used to predict binding site motifs, which are then experimentally validated using EMSAs (Figure 2). The functional predictions and gene targets for the RR may then be extended to closely related species that encode orthologous RRs by scanning those genomes for similar binding motifs (Figure 2). The DAP-chip method can provide a wealth of information for a TCS where previously there was none. The method can also be used for any transcriptional regulator if the protein can be purified and DNA binding conditions can be determined, and for any organism of interest with a genome sequence available.
<img alt="Figure 1" fo:content-width="5in" src="/files/ftp_upload/51715/51715fig1highres.jpg" width="500" /> 
Figure 1. The DNA-affinity-purified-chip (DAP-chip) strategy7. The RR gene from the organism of interest is cloned with a carboxy-terminal His-tag into an E. coli expression strain. Purified His-tagged protein is activated by phosphorylation with acetyl phosphate, and mixed with sheared genomic DNA. An aliquot of the binding reaction is saved as input DNA, while the rest is subjected to affinity purification using Ni-NTA resin. The input and the RR-bound DNA are whole genome amplified, and labeled with Cy3 and Cy5, respectively. The labeled DNA is pooled together and hybridized to a tiling array, which is then analyzed to determine the gene targets. Figure modified and reprinted using the creative commons license from7.
<img alt="Figure 2" fo:content-width="5in" src="/files/ftp_upload/51715/51715fig2highres.jpg" width="500" /> 
Figure 2. Summary of workflow. For any purified tagged protein, begin by determining a target using EMSA. Allow protein to bind genomic DNA and then DNA-affinity-purify (DAP) and whole genome amplify (WGA) the enriched and input DNA. If a gene target is known, use qPCR to ensure that the known target is enriched in the protein-bound fraction. If no target could be determined, proceed directly to DNA labeling and array hybridization. If enrichment by qPCR could not be observed, then repeat the protein-gDNA binding and DAP-WGA steps using different protein amounts. Use array analysis to find peaks and map them to target genes. Use the upstream regions of target genes to predict binding site motifs. Validate the motifs experimentally using EMSAs. Use the motif to scan the genomes of related species encoding orthologs of the RR under study, and predict genes targeted in those species as well. Based on the gene targets obtained, the physiological function of the RR and its orthologs may be predicted. Figure modified and reprinted using the creative commons license from7.

<table border="0" cellpadding="0" cellspacing="0" width="719">
	<tbody>
		<tr height="17">
			<td height="17" style="height:17px;width:255px;">Name of Material/Equipment</td>
			<td style="width:248px;">Company</td>
			<td style="width:113px;">Catalog number</td>
			<td style="width:104px;">Comments</td>
		</tr>
		<tr height="34">
			<td height="34" style="height:34px;width:255px;">HisTrapFF column (Ni-Sepharose column)</td>
			<td>GE Lfe Sciences, Pittsburgh, PA, USA</td>
			<td>17-5255-01</td>
			<td></td>
		</tr>
		<tr height="17">
			<td height="17" style="height:17px;width:255px;">Akta explorer (FPLC instrument)</td>
			<td colspan="2">GE Lfe Sciences, Pittsburgh, PA, USA</td>
			<td></td>
		</tr>
		<tr height="17">
			<td height="17" style="height:17px;width:255px;">HiPrep 26/10 Desalting column</td>
			<td>GE Lfe Sciences, Pittsburgh, PA, USA</td>
			<td>17-5087-01</td>
			<td></td>
		</tr>
		<tr height="17">
			<td height="17" style="height:17px;width:255px;">Qiaquick Gel extraction kit</td>
			<td>Qiagen Inc, Valencia, CA, USA</td>
			<td align="right">28704</td>
			<td></td>
		</tr>
		<tr height="17">
			<td height="17" style="height:17px;width:255px;">Biotin-labeled oligonucleotides</td>
			<td>Integrated DNA Technologies</td>
			<td>N/A</td>
			<td></td>
		</tr>
		<tr height="34">
			<td height="34" style="height:34px;width:255px;">6% polyacrylamide-0.5X TBE precast mini DNA retardation gel</td>
			<td>Life Technologies, Grand Island, NY, USA</td>
			<td>EC63652BOX</td>
			<td>Alternately, you can pour your own gel.</td>
		</tr>
		<tr height="17">
			<td height="17" style="height:17px;width:255px;">Nylon membrane</td>
			<td>EMD Millipore, Billerica, MA, USA</td>
			<td>INYC00010</td>
			<td></td>
		</tr>
		<tr height="17">
			<td height="17" style="height:17px;width:255px;">Trans-Blot SD Semi-dry electrophoretic transfer cell</td>
			<td>Biorad, Hercules, CA, USA</td>
			<td>170-3940</td>
			<td></td>
		</tr>
		<tr height="17">
			<td height="17" style="height:17px;width:255px;">Extra thick blot paper, 8 x 13.5 cm</td>
			<td>Biorad, Hercules, CA, USA</td>
			<td>170-3967</td>
			<td></td>
		</tr>
		<tr height="17">
			<td height="17" style="height:17px;width:255px;">UV crosslinker Model XL-1000</td>
			<td>Fisher Scientific</td>
			<td style="width:113px;">11-992-89</td>
			<td></td>
		</tr>
		<tr height="34">
			<td height="34" style="height:34px;width:255px;">Nucleic Acid chemiluminescent detection kit (Pierce)</td>
			<td>Thermo fisher Scientific, Rockford, IL, USA</td>
			<td align="right">89880</td>
			<td></td>
		</tr>
		<tr height="17">
			<td height="17" style="height:17px;width:255px;">Ni-NTA agarose resin</td>
			<td>Qiagen Inc, Valencia, CA, USA</td>
			<td align="right">30210</td>
			<td></td>
		</tr>
		<tr height="119">
			<td height="119" style="height:119px;width:255px;">GenomePlex Whole genome amplification kit (Fragmentation buffer, library preparation buffer, library stabilization solution, library preparation enzyme, 10X amplification master mix, WGA polymerase )</td>
			<td>Sigma-Aldrich, St. Louis, MO, USA</td>
			<td>WGA2-50RXN</td>
			<td></td>
		</tr>
		<tr height="17">
			<td height="17" style="height:17px;width:255px;">Nanodrop ND-1000</td>
			<td colspan="2">Thermo Scientific, Wilmington, DE, USA</td>
			<td>For quantitation of DNA</td>
		</tr>
		<tr height="34">
			<td height="34" style="height:34px;width:255px;">Perfecta Sybr Green SuperMix, with ROX</td>
			<td>Quanta biosciences</td>
			<td>95055-500</td>
			<td>Any Sybr Green PCR mix may be used</td>
		</tr>
		<tr height="34">
			<td height="34" style="height:34px;width:255px;">PlateMax Ultra clear heat sealing film for qPCR</td>
			<td>Axygen</td>
			<td></td>
			<td></td>
		</tr>
		<tr height="34">
			<td height="34" style="height:34px;width:255px;">96 well clear low profile PCR microplate</td>
			<td>Life Technologies, Grand Island, NY, USA</td>
			<td colspan="2">PCR-96-LP-AB-C</td>
		</tr>
		<tr height="34">
			<td height="34" style="height:34px;width:255px;">Applied Biosystems StepOne Plus Real time PCR system</td>
			<td>Life Technologies, Grand Island, NY, USA</td>
			<td align="right">4376600</td>
			<td>Any real time PCR system may be used</td>
		</tr>
		<tr height="17">
			<td height="17" style="height:17px;width:255px;">Qiaquick PCR purification kit</td>
			<td>Qiagen Inc, Valencia, CA, USA</td>
			<td align="right">28104</td>
			<td>Any PCR clean up kit may be used</td>
		</tr>
		<tr height="17">
			<td height="17" style="height:17px;width:255px;">Cy3/Cy5-labeled nonamers</td>
			<td>Trilink biotechnologies, San Diego, CA, USA</td>
			<td colspan="2">N46-0001, N46-0002</td>
		</tr>
		<tr height="34">
			<td height="34" style="height:34px;width:255px;">Klenow polymerase 50,000U/ml, 3&#39;-5&#39; exo-</td>
			<td>New England Biolabs, Ipswich, MA</td>
			<td>M0212M</td>
			<td></td>
		</tr>
		<tr height="17">
			<td height="17" style="height:17px;width:255px;">Hybridization system</td>
			<td>Roche-Nimblegen, Madison, WI, USA</td>
			<td>N/A</td>
			<td>This company no longer makes arrays or related items, so alternate sources such as Agilent or Affymetrix will need to be used,</td>
		</tr>
		<tr height="34">
			<td height="34" style="height:34px;width:255px;">Custom printed microarrays and mixers</td>
			<td>Roche-Nimblegen, Madison, WI, USA</td>
			<td>N/A</td>
			<td></td>
		</tr>
		<tr height="51">
			<td height="51" style="height:51px;width:255px;">Hybridization kit (2X Hybridization buffer, Hybridization component A, Alignment oligo)</td>
			<td>Roche-Nimblegen, Madison, WI, USA</td>
			<td>N/A</td>
			<td></td>
		</tr>
		<tr height="34">
			<td height="34" style="height:34px;width:255px;">Wash buffer kit (10X Wash buffer I, II, III, 1 M DTT)</td>
			<td>Roche-Nimblegen, Madison, WI, USA</td>
			<td>N/A</td>
			<td></td>
		</tr>
		<tr height="17">
			<td height="17" style="height:17px;width:255px;">GenePix 4200A microarray scanner</td>
			<td colspan="2">Molecular Devices, Sunnyvale CA, USA</td>
			<td>This model has been replaced by superior ones</td>
		</tr>
		<tr height="17">
			<td height="17" style="height:17px;width:255px;">GenePix Pro microarray software</td>
			<td colspan="2">Molecular Devices, Sunnyvale CA, USA</td>
			<td></td>
		</tr>
		<tr height="34">
			<td height="34" style="height:34px;width:255px;">Nimblescan v.2.4, ChIP-chip analysis software</td>
			<td>Roche-Nimblegen, Madison, WI, USA</td>
			<td>N/A</td>
			<td></td>
		</tr>
	</tbody>
</table>

dna-affinity-purified chip (dap-chip) method to determine gene targets for bacterial two component regulatory systems

In vivo methods such as ChIP-chip are well-established techniques used to determine global gene targets for transcription factors. However, they are of limited use in exploring bacterial two component regulatory systems with uncharacterized activation conditions. Such systems regulate transcription only when activated in the presence of unique signals. Since these signals are often unknown, the in vitro microarray based method described in this video article can be used to determine gene targets and binding sites for response regulators. This DNA-affinity-purified-chip method may be used for any purified regulator in any organism with a sequenced genome. The protocol involves allowing the purified tagged protein to bind to sheared genomic DNA and then affinity purifying the protein-bound DNA, followed by fluorescent labeling of the DNA and hybridization to a custom tiling array. Preceding steps that may be used to optimize the assay for specific regulators are also described. The peaks generated by the array data analysis are used to predict binding site motifs, which are then experimentally validated. The motif predictions can be further used to determine gene targets of orthologous response regulators in closely related species. We demonstrate the applicability of this method by determining the gene targets and binding site motifs and thus predicting the function for a sigma54-dependent response regulator DVU3023 in the environmental bacterium Desulfovibrio vulgaris Hildenborough.

The above method was applied to determine the global gene targets of the RRs in the model sulfate reducing bacterium Desulfovibrio vulgaris Hildenborough7. This organism has a large number of TCSs represented by over 70 RRs, indicating the wide variety of possible signals that it senses and responds to. In vivo analyses on the functions of these signaling systems are hard to perform since their signals and thus their activating conditions are unknown. Here the DAP-chip method was used to determine the gene targets and thus predict possible functions for a representative RR DVU3023.
DVU3023 is a sigma-54-dependent RR encoded in an operon with its cognate HK (Figure 3A). The C-terminal His-tagged gene was cloned into and purified from E. coli. For initial target determination, the purified RR was tested for binding to its downstream operon which is a ten gene-operon (DVU3025-3035) consisting of lactate uptake and oxidation genes. RR DVU3023 shifted the upstream region of DVU3025 (Figure 3B). Next RR DVU3023 was allowed to bind to sheared D. vulgaris Hildenborough genomic DNA. Although phosphorylation was not required for binding to the promoter region of DVU3025, acetyl phosphate was added to the reaction in case it was required for binding to other promoters. Following affinity purification of the protein-bound DNA fraction, qPCR was used to show that the upstream region of DVU3025 is enriched (8.45 fold) in the protein-bound fraction (CT= 6.9) relative to the input DNA (Ct =9.98) (Figure 3C), thus indicating that the binding conditions used were appropriate for RR DVU3023. Lack of enrichment of the promoter region of a randomly chosen gene (DVU0013) was used as a negative control.
The protein-bound and input DNA samples were then labeled fluorescently and hybridized to a D. vulgaris Hildenborough tiling array that had a high probe density in the intergenic regions. The top four peaks were chosen as the most likely targets for DVU3023 (Figure 3D). These four peaks were followed by several others, which were also identified in DAP-chip analyses for several other RRs and hence appear to be sticky DNA (Table 1). Figure 3E is a schematic representation of the genes regulated by DVU3023. The positive target DVU3025 was the first peak obtained with the highest score. Two gene targets are two other singly encoded lactate permeases (DVU2451 and DVU3284). The fourth gene target does not lie in an upstream region, but in the intergenic region between two convergently transcribed genes/operons (DVU0652 and DVU0653). This is a large intergenic region, and additionally also encodes a predicted sigma54-dependent promoter. It is possible that there is an undiscovered sRNA encoded in this region that is regulated by RR DVU3023.
Using the upstream regions of the targets obtained by DAP-chip, MEME was used to predict a binding site motif (Figure 3F, Table 2). EMSA substrates carrying the specific motif upstream of DVU3025 were designed to confirm that RR DVU3023 recognizes and binds the predicted motif. The motif was further validated by making substitutions in the conserved bases within the motif which eliminated the binding shift (Figure 3G). This validated motif was then used in Perl-generated scripts to scan the genome sequences of other closely related sulfate reducing bacteria that had orthologs for DVU3023. Loci were chosen as possible gene targets when the motif was located in upstream regions of open reading frames (Table 3). Using the motif sequences predicted for the orthologous RRs, a consensus binding site motif was generated (Figure 3F) which closely resembled the one obtained for D. vulgaris Hildenborough alone.
<img alt="Figure 3" fo:content-width="5in" src="/files/ftp_upload/51715/51715fig3highres.jpg" width="500" /> 
Figure 3.&#160;Determining genomic targets for D. vulgaris response regulator DVU3023. A. DVU3023 is encoded in an operon with its cognate HK. The downstream operon has a sigma54-dependent promoter (bent black arrow) and was used as a candidate target gene. B. Purified RR DVU3023 bound and shifted the upstream region of DVU30257. C. q-PCR of upstream regions of DVU3025 (positive target) and DVU0013 (chosen as a negative control). E is protein-bound enriched DNA fraction, I is input DNA. D. Top four peaks obtained after DAP-chip analysis. Start and End refer to DNA coordinates at the start and end of the peaks; score refers to the log2R ratio of the fourth highest probe in the peak; Fdr = false discovery rate; cutoff_p is the cutoff percentage at which the peak was identified. E. Schematic representation of the gene targets for DVU3023 based on the DAP-chip peaks. Numbers in boxes indicate the peak number in D. The HK-RR genes are in green, target genes in blue and other genes in grey. Black bent arrows are sigma54-dependent promoters, green filled circles are predicted binding site motifs. Gene names are as follows: por &#8211;pyruvate ferredoxin oxidoreductase; llp &#8211; lactate permease; glcD- glycolate oxidase; glpC- Fe-S cluster binding protein; pta &#8211; phosphotransacetylase; ack- acetate kinase; lldE &#8211; lactate oxidase subunit; lldF/G &#8211; lactate oxidase subunit; MCP &#8211; methyl-accepting chemotaxis protein. F. Weblogo8 images of the predicted binding site motif. Top &#8211; derived from DAP-chip targets; Bottom &#8211; derived from binding sites present in genomes with orthologs of DVU30237. G. Validation of predicted binding site motif using EMSA. DVU3023 shifted the wild type motif (w) but not the modified motif (m). Sequences for the w and m motifs are shown on the right7. Figure modified and reprinted using the creative commons license from Rajeev&#160;et al7.
<img alt="Table 1" fo:content-width="7in" src="/files/ftp_upload/51715/51715table1.jpg" width="700" /> 
Table 1.&#160;Top 20 peaks from the DAP-chip array analysis. The table is divided into three sections. Peak attributes show details of the peaks as generated by array analysis software, where Location refers to whether the peak was found in the genome or the extrachromosomal plasmid, Start and End refer to the start and end loci for the peak, Score refers to the log2R ratios for the fourth highest probe in the peak, Fdr refers to the false discovery rate value, and cutoff_p is the cutoff percentage at which the peak was identified. The other two sections Start coordinate mapping and End coordinate mapping map the start and end loci, respectively, of the peak to the gene. In these sections Gene strand refers to the strand coding the gene, Offset indicates distance from the start of the gene (positive values indicate loci is upstream of gene, while negative values indicate loci is within the gene), Overlap gene value is TRUE if the locus overlaps a gene, DVU refers to the DVU# of the gene that the coordinates map to, and Description indicates the gene annotation. Table modified and reprinted using the creative commons license from Rajeev&#160;et al7. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/51715/51715table1.jpg" target="_blank">Please click here to view a larger version of this table.</a>
<img alt="Table 2" fo:content-width="5in" src="/files/ftp_upload/51715/51715table2.jpg" width="500" /> 
Table 2.&#160;Sequences used to build the consensus DAP-chip target based motif in Figure 3. Table modified and reprinted using the creative commons license from Rajeev&#160;et al7.
<img alt="Table 3" fo:content-width="6in" src="/files/ftp_upload/51715/51715table3.jpg" width="600" /> 
Table 3.&#160;Binding site motifs for DVU3023 orthologs present in other sequenced Desulfovibrio and related species. For each genome scanned, the organism name is indicated, followed by the locus tag for the DVU3023 ortholog and its percent identity to DVU3023 in parentheses. Score indicates the value assigned by the Perl program based on similarity to the input sequences. Description states the gene annotations for the genes in the target operon. Table modified and reprinted using the creative commons license from Rajeev&#160;et al7. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/51715/51715table3.jpg" target="_blank">Please click here to view a larger version of this table.</a>

Watch this Scientific Journal Video about DNA-affinity-purified Chip DAP-chip Method to Determine Gene Targets for Bacterial Two component Regulatory Systems at JoVE.com

DNA-affinity-purified Chip DAP-chip Method to Determine Gene Targets for Bacterial Two component Regulatory Systems

This video article describes an in vitro microarray based method to determine the gene targets and binding sites for two component system response regulators.

DNA-affinity-purified Chip (DAP-chip) Method to Determine Gene Targets for Bacterial Two component Regulatory Systems

In vivo methods such as ChIP-chip are well-established techniques used to determine global gene targets for transcription factors. However, they are of ...

dna-affinity-purified-chip-dap-chip-method-to-determine-gene-targets

Nanyang Technological University

Research

JoVE Journal

Biology

16.7K Views.  Lawrence Berkeley National Laboratory. This video article describes an in vitro microarray based method to determine the gene targets and binding sites for two component system response regulators. 

Video: DNA-affinity-purified Chip DAP-chip Method to Determine Gene Targets for Bacterial Two component Regulatory Systems

Bacterial Gene Expression Analysis Using Microarrays

Hi-C: A Method to Study the Three-dimensional Architecture of Genomes.

The three-dimensional folding of chromosomes compartmentalizes the genome and and can bring distant functional elements, such as promoters and enhancers, into close spatial proximity 2-6. Deciphering the relationship between chromosome organization and genome activity will aid in understanding genomic processes, like transcription and replication. However, little is known about how chromosomes fold. Microscopy is unable to distinguish large numbers of loci simultaneously or at high resolution. To date, the detection of chromosomal interactions using chromosome conformation capture (3C) and its subsequent adaptations required the choice of a set of target loci, making genome-wide studies impossible 7-10.
We developed Hi-C, an extension of 3C that is capable of identifying long range interactions in an unbiased, genome-wide fashion. In Hi-C, cells are fixed with formaldehyde, causing interacting loci to be bound to one another by means of covalent DNA-protein cross-links. When the DNA is subsequently fragmented with a restriction enzyme, these loci remain linked. A biotinylated residue is incorporated as the 5' overhangs are filled in. Next, blunt-end ligation is performed under dilute conditions that favor ligation events between cross-linked DNA fragments. This results in a genome-wide library of ligation products, corresponding to pairs of fragments that were originally in close proximity to each other in the nucleus. Each ligation product is marked with biotin at the site of the junction. The library is sheared, and the junctions are pulled-down with streptavidin beads. The purified junctions can subsequently be analyzed using a high-throughput sequencer, resulting in a catalog of interacting fragments.
Direct analysis of the resulting contact matrix reveals numerous features of genomic organization, such as the presence of chromosome territories and the preferential association of small gene-rich chromosomes. Correlation analysis can be applied to the contact matrix, demonstrating that the human genome is segregated into two compartments: a less densely packed compartment containing open, accessible, and active chromatin and a more dense compartment containing closed, inaccessible, and inactive chromatin regions. Finally, ensemble analysis of the contact matrix, coupled with theoretical derivations and computational simulations, revealed that at the megabase scale Hi-C reveals features consistent with a fractal globule conformation.

The Hi-C method allows unbiased, genome-wide identification of chromatin interactions (1). Hi-C couples proximity ligation and massively parallel sequencing. The resulting data can be used to study genomic architecture at multiple scales: initial results identified features such as chromosome territories, segregation of open and closed chromatin, and chromatin structure at the megabase scale.

The three-dimensional folding of chromosomes compartmentalizes the genome and and can bring distant functional elements, such as promoters and enhancers, ...

Chromatin Immunoprecipitation (ChIP) to Assay Dynamic Histone Modification in Activated Gene Expression in Human Cells

In response to a variety of extracellular ligands, the STAT (signal transducer and activator of transcription) transcription factors are rapidly recruited from their latent state in the cytoplasm to cell surface receptors where they are activated by phosphorylation at a single tyrosine residue1. They then dimerize and translocate to the nucleus to drive the transcription of target genes, affecting growth, differentiation, homeostasis and the immune response. Not surprisingly, given their widespread involvement in normal cell processes, dysregulation of STAT function contributes to human disease, particularly to cancers2 and autoimmune diseases3. It is well established that transcription is regulated by alterations to the chromatin template4,5. These alterations include the activities of ATP-dependent complexes, as well as covalent histone modifications and DNA methylation6. Because STAT activation of gene expression is both rapid and transient, it requires specific mechanisms for modulating the chromatin template at STAT-dependent gene loci. To define these mechanisms, we characterize the histone modifications and the enzymatic activities that generate them at gene loci that respond to STAT signaling. This protocol describes chromatin immunoprecipitation, a method that is valuable for the study of STAT signaling to chromatin in activated gene expression.

Chromatin Immunoprecipitation ChIP to Assay Dynamic Histone Modification in Activated Gene Expression in Human Cells

This protocol describes how chromatin immunoprecipitation (ChIP) is used to study the dynamic alterations to the chromatin template that regulate transcription induced by a signal transduction pathway.

In response to a variety of extracellular ligands, the STAT (signal transducer and activator of transcription) transcription factors are rapidly recruited ...

Genome-wide Analysis using ChIP to Identify Isoform-specific Gene Targets

Recruitment of transcriptional and epigenetic factors to their targets is a key step in their regulation. Prominently featured in recruitment are the protein domains that bind to specific histone modifications. One such domain is the plant homeodomain (PHD), found in several chromatin-binding proteins. The epigenetic factor RBP2 has multiple PHD domains, however, they have different functions (Figure 4). In particular, the C-terminal PHD domain, found in a RBP2 oncogenic fusion in human leukemia, binds to trimethylated lysine 4 in histone H3 (H3K4me3)1. The transcript corresponding to the RBP2 isoform containing the C-terminal PHD accumulates during differentiation of promonocytic, lymphoma-derived, U937 cells into monocytes2. Consistent with both sets of data, genome-wide analysis showed that in differentiated U937 cells, the RBP2 protein gets localized to genomic regions highly enriched for H3K4me33. Localization of RBP2 to its targets correlates with a decrease in H3K4me3 due to RBP2 histone demethylase activity and a decrease in transcriptional activity. In contrast, two other PHDs of RBP2 are unable to bind H3K4me3. Notably, the C-terminal domain PHD of RBP2 is absent in the smaller RBP2 isoform4. It is conceivable that the small isoform of RBP2, which lacks interaction with H3K4me3, differs from the larger isoform in genomic location. The difference in genomic location of RBP2 isoforms may account for the observed diversity in RBP2 function. Specifically, RBP2 is a critical player in cellular differentiation mediated by the retinoblastoma protein (pRB). Consistent with these data, previous genome-wide analysis, without distinction between isoforms, identified two distinct groups of RBP2 target genes: 1) genes bound by RBP2 in a manner that is independent of differentiation; 2) genes bound by RBP2 in a differentiation-dependent manner.
To identify differences in localization between the isoforms we performed genome-wide location analysis by ChIP-Seq. Using antibodies that detect both RBP2 isoforms we have located all RBP2 targets. Additionally we have antibodies that only bind large, and not small RBP2 isoform (Figure 4). After identifying the large isoform targets, one can then subtract them from all RBP2 targets to reveal the targets of small isoform. These data show the contribution of chromatin-interacting domain in protein recruitment to its binding sites in the genome.

Here we are presenting a chromatin immunoprecipitation (ChIP) procedure for genome-wide location analysis of protein isoforms that differ in a histone-binding domain. We are applying it to ChIP-Seq analysis to identify the targets of the KDM5A/JARID1A/RBP2 histone demethylase.

Recruitment of transcriptional and epigenetic factors to their targets is a key step in their regulation. Prominently featured in recruitment are the ...

Using SecM Arrest Sequence as a Tool to Isolate Ribosome Bound Polypeptides

Extensive research has provided ample evidences suggesting that protein folding in the cell is a co-translational process1-5. However, the exact pathway that polypeptide chain follows during co-translational folding to achieve its functional form is still an enigma. In order to understand this process and to determine the exact conformation of the co-translational folding intermediates, it is essential to develop techniques that allow the isolation of RNCs carrying nascent chains of predetermined sizes to allow their further structural analysis.
SecM (secretion monitor) is a 170 amino acid E. coli protein that regulates expression of the downstream SecA (secretion driving) ATPase in the secM-secA operon6. Nakatogawa and Ito originally found that a 17 amino acid long sequence (150-FSTPVWISQAQGIRAGP-166) in the C-terminal region of the SecM protein is sufficient and necessary to cause stalling of SecM elongation at Gly165, thereby producing peptidyl-glycyl-tRNA stably bound to the ribosomal P-site7-9. More importantly, it was found that this 17 amino acid long sequence can be fused to the C-terminus of virtually any full-length and/or truncated protein thus allowing the production of RNCs carrying nascent chains of predetermined sizes7. Thus, when fused or inserted into the target protein, SecM stalling sequence produces arrest of the polypeptide chain elongation and generates stable RNCs both in vivo in E. coli cells and in vitro in a cell-free system. Sucrose gradient centrifugation is further utilized to isolate RNCs.
The isolated RNCs can be used to analyze structural and functional features of the co-translational folding intermediates. Recently, this technique has been successfully used to gain insights into the structure of several ribosome bound nascent chains10,11. Here we describe the isolation of bovine Gamma-B Crystallin RNCs fused to SecM and generated in an in vitro translation system.

We describe here a technique that is now routinely used to isolate stably bound ribosome nascent chain complexes (RNCs). This technique takes advantage of the discovery that a 17 amino acid long SecM "arrest sequence" can halt translation elongation in a prokaryotic (E. coli) system, when inserted into (or fused to the C-terminus) of virtually any protein.

Extensive research has provided ample evidences suggesting that protein folding in the cell is a co-translational process1-5. However, the exact pathway ...

The Use of Ex Vivo Whole-organ Imaging and Quantitative Tissue Histology to Determine the Bio-distribution of Fluorescently Labeled Molecules

Fluorescent labeling is a well-established process for examining the fate of labeled molecules under a variety of experimental conditions both in vitro and in vivo. Fluorescent probes are particularly useful in determining the bio-distribution of administered large molecules, where the addition of a small-molecule fluorescent label is unlikely to affect the kinetics or bio-distribution of the compound. A variety of methods exist to examine bio-distribution that vary significantly in the amount of effort required and whether the resulting measurements are fully quantitative, but using multiple methods in conjunction can provide a rapid and effective system for analyzing bio-distributions.
Ex vivo whole-organ imaging is a method that can be used to quickly compare the relative concentrations of fluorescent molecules within tissues and between multiple types of tissues or treatment groups. Using an imaging platform designed for live-animal or whole-organ imaging, fluorescence within intact tissues can be determined without further processing, saving time and labor while providing an accurate picture of the overall bio-distribution. This process is ideal in experiments attempting to determine the tissue specificity of a compound or for the comparison of multiple different compounds. Quantitative tissue histology on the other hand requires extensive further processing of tissues in order to create a quantitative measure of the labeled compounds. To accurately assess bio-distribution, all tissues of interest must be sliced, scanned, and analyzed relative to standard curves in order to make comparisons between tissues or groups. Quantitative tissue histology is the gold standard for determining absolute compound concentrations within tissues.
Here, we describe how both methods can be used together effectively to assess the ability of different administration methods and compound modifications to target and deliver fluorescently labeled molecules to the central nervous system1.

The Use of Ex Vivo Whole-organ Imaging and Quantitative Tissue Histology to Determine the Bio-distribution of Fluorescently Labeled Molecules

Ex vivo whole-organ imaging is a rapid method for determining the relative concentrations of fluorescently labeled compounds within and between tissues or treatment groups. Conversely, quantitative fluorescence histology, while more labor intensive, allows for the quantification of the absolute tissue levels of labeled molecules.

Fluorescent labeling is a well-established process for examining the fate of labeled molecules under a variety of experimental conditions both in vitro ...

Using a Whole-mount Immunohistochemical Method to Study the Innervation of the Biliary Tract in Suncus murinus

This work describes the whole-mount immunohistochemistry staining method in detail, using neurofilament protein antibody to label the innervation of the biliary tract in Suncus murinus (S. murinus&#160;). First, the specimen was dissected from S. murinus and fixed in 4% paraformaldehyde (PFA). Second, an enzymatic treatment and potential endogenous peroxidase inactivation were performed. The specimen was then exposed to the primary antibody, anti-neurofilament protein antibody, for 3-6 days. It was then incubated with the secondary antibody conjugated with horseradish peroxidase. The color reaction was revealed by reacting the specimen with a 3,3&#39;-diaminobenzidine (DAB) substrate. This method can be applied to analyze the innervation of all visceral organs. Furthermore, this protocol can also be adapted to test other neuronal antibodies, but optimization of the antibodies should be done first. This method was originally introduced by Kuratani and Tanaka1,2,3.

Using a Whole-mount Immunohistochemical Method to Study the Innervation of the Biliary Tract in Suncus murinus

A whole-mount immunohistochemical approach, to visualize neurofilament protein expression in the extrahepatic biliary tract in Suncus murinus. is presented here. This protocol can be used to analyze the innervation of all visceral organs in S. murinus or other species.

This work describes the whole-mount immunohistochemistry staining method in detail, using neurofilament protein antibody to label the innervation of the ...

ODELAY: A Large-scale Method for Multi-parameter Quantification of Yeast Growth

Growth phenotypes of microorganisms are a strong indicator of their underlying genetic fitness and can be segregated into 3 growth regimes: lag-phase, log-phase, and stationary-phase. Each growth phase can reveal different aspects of fitness that are related to various environmental and genetic conditions. High-resolution and quantitative measurements of all 3 phases of growth are generally difficult to obtain. Here we present a detailed method to characterize all 3 growth phases on solid media using an assay called One-cell Doubling Evaluation of Living Arrays of Yeast (ODELAY). ODELAY quantifies growth phenotypes of individual cells growing into colonies on solid media using time-lapse microscopy. This method can directly observe population heterogeneity with each growth parameter in genetically identical cells growing into colonies. This population heterogeneity offers a unique perspective for understanding genetic and epigenetic regulation, and responses to genetic and environmental perturbations. While the ODELAY method is demonstrated using yeast, it can be utilized on any colony forming microorganism that is visible by bright field microscopy.

We present a method for quantifying growth phenotypes of individual yeast cells as they grow into colonies on solid media using time-lapse microscopy termed, One-cell Doubling Evaluation of Living Arrays of Yeast (ODELAY). Population heterogeneity of genetically identical cells growing into colonies can be directly observed and quantified.

Growth phenotypes of microorganisms are a strong indicator of their underlying genetic fitness and can be segregated into 3 growth regimes: lag-phase, ...

The Lactate Dehydrogenase Sequestration Assay &#8212; A Simple and Reliable Method to Determine Bulk Autophagic Sequestration Activity in Mammalian Cells

Bulk autophagy is characterized by the sequestration of large portions of cytoplasm into double/multi-membrane structures termed autophagosomes. Here a simple protocol to monitor this process is described. Moreover, typical results and experimental validation of the method under autophagy-inducing conditions in various types of cultured mammalian cells are provided. During bulk autophagy, autophagosomes sequester cytosol, and thereby also soluble cytosolic proteins, alongside other autophagic cargo. LDH is a stable and highly abundant, soluble cytosolic enzyme that is non-selectively sequestered into autophagosomes. The amount of LDH sequestration therefore reflects the amount of bulk autophagic sequestration. To efficiently and accurately determine LDH sequestration in cells, we employ an electrodisruption-based fractionation protocol that effectively separates sedimentable from cytosolic LDH, followed by measurement of enzymatic activity in sedimentable fractions versus whole-cell samples. Autophagic sequestration is determined by subtracting the proportion of sedimentable LDH in untreated cells from that in treated cells. The advantage of the LDH sequestration assay is that it gives a quantitative measure of the autophagic sequestration of endogenous cargo, as opposed to other methods that either involve ectopic expression of sequestration probes or semi-quantitative protease protection analyses of autophagy markers or receptors.

Here a simple and well-validated protocol for measuring bulk autophagic sequestration activity in mammalian cells is described. The method is based on quantifying the proportion of lactate dehydrogenase (LDH) in sedimentable cell fractions compared to total cellular LDH levels.

Bulk autophagy is characterized by the sequestration of large portions of cytoplasm into double/multi-membrane structures termed autophagosomes. Here a ...

High-Throughput Live Imaging of Microcolonies to Measure Heterogeneity in Growth and Gene Expression

Precise measurements of between- and within-strain heterogeneity in microbial growth rates are essential for understanding genetic and environmental inputs into stress tolerance, pathogenicity, and other key components of fitness. This manuscript describes a microscope-based assay that tracks approximately 105&#160;Saccharomyces cerevisiae microcolonies per experiment. After automated time-lapse imaging of yeast immobilized in a multiwell plate, microcolony growth rates are easily analyzed with custom image-analysis software. For each microcolony, expression and localization of fluorescent proteins and survival of acute stress can also be monitored. This assay allows precise estimation of strains&#39; average growth rates, as well as comprehensive measurement of heterogeneity in growth, gene expression, and stress tolerance within clonal populations.

Yeast growth phenotypes are precisely measured through highly parallel time-lapse imaging of immobilized cells growing into microcolonies. Simultaneously, stress tolerance, protein expression, and protein localization can be monitored, generating integrated datasets to study how environmental and genetic differences, as well as gene-expression heterogeneity among isogenic cells, modulate growth.

Precise measurements of between- and within-strain heterogeneity in microbial growth rates are essential for understanding genetic and environmental ...

DNA-affinity-purified Chip (DAP-chip) Method to Determine Gene Targets for Bacterial Two component Regulatory Systems

Summary

Explore More Videos

Chapters in this video

Bacterial Gene Expression Analysis Using Microarrays

Hi-C: A Method to Study the Three-dimensional Architecture of Genomes.

Chromatin Immunoprecipitation (ChIP) to Assay Dynamic Histone Modification in Activated Gene Expression in Human Cells

Genome-wide Analysis using ChIP to Identify Isoform-specific Gene Targets

Using SecM Arrest Sequence as a Tool to Isolate Ribosome Bound Polypeptides

The Use of Ex Vivo Whole-organ Imaging and Quantitative Tissue Histology to Determine the Bio-distribution of Fluorescently Labeled Molecules

Using a Whole-mount Immunohistochemical Method to Study the Innervation of the Biliary Tract in Suncus murinus

ODELAY: A Large-scale Method for Multi-parameter Quantification of Yeast Growth

The Lactate Dehydrogenase Sequestration Assay — A Simple and Reliable Method to Determine Bulk Autophagic Sequestration Activity in Mammalian Cells

High-Throughput Live Imaging of Microcolonies to Measure Heterogeneity in Growth and Gene Expression