PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of hhmipaabout author manuscriptssubmit a manuscriptHHMI Howard Hughes Medical Institute; Author Manuscript; Accepted for publication in peer reviewed journal
 
Curr Protoc Plant Biol. Author manuscript; available in PMC 2017 August 4.
Published in final edited form as:
PMCID: PMC5544034
HHMIMSID: HHMIMS836730

Profiling of transcription factor binding events by chromatin immunoprecipitation sequencing (ChIP-seq)

Abstract

In multi-cellular organisms, gene expression is orchestrated by thousands of transcription factors (TF). Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is a robust tool to investigate gene expression because this technique profiles in vivo protein-DNA interaction at a genome-wide scale. Eight years after the first ChIP-seq paper, there are limited reports of ChIP-seq experiments in plants, especially for sequence-specific DNA binding TFs This lag greatly prevents our understanding of transcriptional regulation in an entire kingdom. In order to bridge the technical gap, we describe a ChIP-seq procedure that we have successfully applied to dozens of sequence-specific DNA binding TFs. The basic protocol includes procedures to isolate nuclei, sonicate chromatin, immunoprecipitate TF-DNA complex, and recover ChIP-enriched DNA fragments. The support protocol also describes practices to optimize library preparation by a gel-free DNA size selection. Lastly, examples are given to optimize library amplification using real-time PCR.

Keywords: ChIP-seq, transcription factor, Arabidopsis

INTRODUCTION

ChIP-seq is a robust approach to investigate transcriptional regulation in vivo (Ren et al., 2000; Johnson et al., 2007). It is an improvement over ChIP-chip with better signal-to-noise ratio, unambiguous and genome-wide sequence information. ChIP-seq can be applied to any protein that is associated with chromatin, including histones, chromatin remodelers, RNA polymerases, mediators, general and sequence-specific DNA binding transcription factors. Sequence-specific DNA binding TFs are of particular interest because of their large number and functional diversity; yet they are also the most technically challenging to study because of their low protein abundance and rapid turnover in cells. The progress of ChIP-seq in plants has lagged far behind other eukaryotic organisms: a mere 1% of all ChIP-seq samples submitted to Gene Expression Omnibus (GEO) are from plants. Eight years after the publication of the initial ChIP-seq paper, fewer than 30 sequence-specific DNA binding TFs in Arabidopsis have been subjected to ChIP-seq (Heyndrickx et al., 2014; Fan et al., 2014; Pfeiffer et al., 2014). This unit aims to bridge this gap by providing a technical guidance of ChIP-seq experiments in plants.

BASIC PROTOCOL 1

Enrichment of transcription factor bound DNA fragments by chromatin immunoprecipitation

Introductory paragraph

This protocol describes a step-by-step ChIP procedure. It contains experimental modules for cross-linking samples using formaldehyde, preparation of antibody beads, isolation of nuclei, sonication of chromatin, immunoprecipitation of the TF-DNA complex, and recovery of ChIP-enriched DNA fragments. Arabidopsis etiolated seedlings and floral buds are used as examples here. The protocol can also be applied to many other tissues and plant species. A general overview of ChIP-seq procedures is illustrated in Figure 1, and several key steps are shown in Figure 2.

Figure 1
Flowchart of ChIP-seq procedures with a suggested timeline.
Figure 2
Key steps to prepare chromatin for sonication.

Materials

Sample material

Arabidopsis tissues.

Reagents and solutions

  • 10 mM HEPES-NaOH, pH 7.4
  • formaldehyde (37% in H2O, Sigma Aldrich)
  • fixation buffer (see recipe)
  • 2 M glycine
  • PBS-BSA (see recipe)
  • Dynabeads® protein G / Dynabeads protein A (Thermo Fisher Scientific)
  • ChIP-grade antibody (10 μg per ChIP)
  • Liquid nitrogen
  • Extraction buffer I (see recipe)
  • Extraction buffer II (see recipe)
  • Extraction buffer III (see recipe)
  • Nuclei lysis buffer (see recipe)
  • ChIP dilution buffer (see recipe)
  • Low salt wash buffer (see recipe)
  • High salt wash buffer (see recipe)
  • Final wash buffer (see recipe)
  • Elution buffer (see recipe)
  • Proteinase K (20 mg/ml, Thermo Fisher Scientific)
  • Phenol:chloroform:isoamyl alcohol (25:24:1, pH 8.0)
  • 5 M NaCl
  • 100% ethanol
  • glycogen (5 mg/ml)

Equipment

  • Vacuum pump (ME 4 NT diaphragm pump, VACUUBRAND)
  • Desiccator with a stopcock (Bel-Art)
  • DynaMag™-2 Magnet (Thermo Fisher Scientific)
  • Mortar
  • Pestle
  • Spatula
  • 50 ml conical centrifuge tubes
  • Low binding DNase- and RNase-free microcentrifuge tubes and pipette tips
  • Cheesecloth
  • Miracloth (EMD Millipore)
  • Refrigerated centrifuge for 50 ml tubes
  • Refrigerated microcentrifuge
  • Bioruptor® (Diagenode)
  • Tube rotator
  • Thermomixer (Eppendorf)
  • Phase Lock Gel® heavy (5 PRIME) / MaXtract High Density (Qiagen)
  • Qubit® fluorometer (Thermo Fisher Scientific) / TapeStation (Agilent Technologies)
  • Note: use RNase- and DNase-free reagents, deionized distilled water, and low binding plasticware, and barrier pipette tips throughout the experiment.

Cross-link plant samples with 1% formaldehyde

  1. Harvest plant materials into fixation buffer. Use 30 ml fixation buffer per gram of plant material.

    Use fresh formaldehyde to achieve efficient and reproducible fixation. Avoid handling samples with amine-rich solutions such as Tris buffer before fixation. The amount of samples needed depends on the TF expression level and tissue type. In our hands, 3–4 grams of etiolated seedlings or 1–2 grams of flower buds are sufficient to ChIP many sequence-specific DNA binding TFs. More tissue can be used if the abundance of the TF is very low.

  2. Cross-link sample with fixation buffer under vacuum infiltration for 20 minutes. Break and re-apply vacuum after 5 and 15 minutes to facilitate penetration of formaldehyde into the sample (Figure 2A).

    Fixation time can be adjusted according to the thickness of plant tissue. In general, do not cross-link samples for more than 25 minutes for Arabidopsis tissues; otherwise over-fixation may decrease sonication efficiency.

  3. Quench formaldehyde by adding 0.67 ml of 2 M fresh glycine to every 10 ml of fixation buffer. Apply vacuum for another 5 minutes.
  4. Rinse sample with 10 mM HEPES-NaOH, pH7.4. Remove excess liquid from sample.
  5. (optional) Check cross-linked samples by western blot using antibody against TF-of-interest.

    Use uncross-linked sample as a control. A smeared band or bands of higher molecular weight is expected in cross-linked samples.

  6. Snap freeze sample in liquid nitrogen. Store at − 80 °C until use.

    Sample can be stored for at least a few weeks.

Preparation of antibody beads

Start this step 6+ hours before sonication.

  • 7.
    For each ChIP, add 100 μl Dynabeads® protein G/A to a microcentrifuge tube.

    Choose Dynabeads® protein G or protein A according to their affinity for antibody species and subclasses of IgG.

  • 8.
    Wash beads with 1.5 ml pre-chilled PBS-BSA for three times. Collect beads by DynaMag™-2 Magnet and remove supernatant after each round of wash.
  • 9.
    Mix 10 μg of ChIP-grade antibody and 250 μl of PBS-BSA with Dynabeads® protein G/A.

    The antibody is a key factor to a successful ChIP. Based on our experience, GFP tagged TFs ChIPped with a polyclonal GFP antibody (Thermo Fisher Scientific, cat. # A-11122) or FLAG tagged TFs ChIPped with a monoclonal FLAG antibody (Sigma Aldrich, cat. # F1804) usually give a strong signal and clean background. It is also crucial to include “mock IP” controls. If an antibody against a native protein is used in ChIP, mock IP can be carried out by using IgG from the same species in which the antibody was raised. If an antibody against an epitope tag is used, mock IP can either use the same antibody to ChIP wild-type samples, or even better use the same antibody to ChIP transgenic plants expressing the epitope tag alone.

  • 10.
    Bind antibody to Dynabeads® protein G/A on a tube rotator at 4 °C for 6+ hours.

Isolation of nuclei

  • 11.
    Grind cross-linked samples from step 6 in liquid nitrogen to a fine powder.

    Carry out all procedures at 4 °C from step 11 to step 30. Many buffers contain detergent. Avoid foaming in all steps.

  • 12.
    Mix the powder with seven volumes of Extraction Buffer I in a 50 ml conical centrifuge tube.
  • 13.
    Cut cheesecloth and Miracloth to 4 inch squares. Layer three pieces of cheesecloth above one piece of Miracloth. Pre-wet cheesecloth and Miracloth with Extraction Buffer I.
  • 14.
    Filter plant homogenate from step 13 through cheesecloth and Miracloth into a fresh 50 ml conical centrifuge tube (Figure 2B).
  • 15.
    Repeat filtration through one layer of pre-wet Miracloth.
  • 16.
    Centrifuge filtrate at 2,880 g for 20 minutes (Figure 2C). Discard supernatant.
  • 17.
    Resuspend pellet gently and thoroughly with 1.2 ml of Extraction Buffer II. Transfer sample to a 1.5 ml microcentrifuge tube.
  • 18.
    Centrifuge sample at 12,000 g for 10 minutes (Figure 2D). Discard supernatant.

    For light-grown samples, repeat steps 17 and 18 as needed if too many chloroplasts are present as evidenced by a deep green pellet.

  • 19.
    Resuspend pellet gently and thoroughly with 500 μl of Extraction Buffer III. Carefully layer the resuspended pellet on top of another 500 μl of Extraction Buffer III in a fresh microcentrifuge tube (Figure 2E).
  • 20.
    Centrifuge sample at 16,000 g for an hour (Figure 2F). Discard supernatant.

Nuclei lysis and chromatin sonication

  • 21.
    Resuspend pellet from step 20 gently and thoroughly with 300 μl of Nuclei Lysis Buffer.
  • 22.
    Gently rotate sample for 20 minutes on the tube rotator.

    This is to increase sonication efficiency (F. Turck, personal communication).

  • 23.
    Sonicate chromatin in Bioruptor® in cold room for 25 cycles of 30-second ON and 2-minute OFF at HIGH setting. Gently vortex and then centrifuge tubes after every 10 cycles of sonication.

    Sonication efficiency of Bioruptor® may vary greatly by model. Adjust settings accordingly for best performance. Carry out a western blot to ensure sonication setup does not destroy protein or protein-DNA complex. Wear ear protection as sonication generates high-frequency sound waves that may damage hearing.

  • 24.
    Centrifuge sonicated chromatin at maximum speed for 10 minutes (Figure 2G). Transfer nuclear lysate in supernatant to a fresh microcentrifuge tube.
  • 25.
    Repeat step 24.
  • 26.
    Dilute sonicated nuclear lysate with 4 volumes (ca. 1,200 μl) of ChIP dilution buffer.
  • 27.
    Collect antibody beads from step 10 on a magnet. Wash antibody beads with 1.5 ml PBS-BSA for three times. Collect beads on magnet after each round of wash. Discard supernatant.
  • 28.
    Freeze 30 μl diluted nuclear lysate as input sample. Mix the remaining lysate with antibody beads for ChIP.

ChIP and reverse crosslink

  • 29.
    Mix sonicated nuclear lysate with antibody beads overnight using tube rotator. Collect antibody beads on magnet. Discard supernatant **.
  • 30.
    Wash beads for a total of six rounds with the following buffers. Collect beads on magnet and remove supernatant after each round of wash **.
    • One round of quick wash, and one round of 5-minute wash on tube rotator with Low Salt Wash Buffer.
    • One round of quick wash, and one round of 5-minute wash on tube rotator with High Salt Wash Buffer.
    • One round of quick wash, and one round of 5-minute wash on tube rotator with Final Wash Buffer.
  • 31.
    Resuspend antibody beads with 100 μl of Elution Buffer. Incubate tubes at 65 °C for 15 minutes to elute TF-DNA complex. Collect beads on magnet. Transfer supernatant to a fresh tube.
  • 32.
    Repeat step 31. Combine eluates from steps 31 and 32 **.
  • 33.
    (** optional) Use a small aliquot of supernatant from IP, wash, and elution steps for western blot to ensure TF-of-interest is efficiently immunoprecipitated by the antibody.
  • 34.
    Take input sample from step 28 out of freezer. Add 170 μl of Elution Buffer.
  • 35.
    Reverse cross-linked samples from steps 32 and 34 in a Thermomixer at 65 °C for 6+ hours.
  • 36.
    Add 1 μl (20 μg) Proteinase K to each sample and incubate at 65 °C for 2 hours.
  • 37.
    Transfer samples to 1.5 – 2 ml Phase Lock Gel® heavy tubes. Extract DNA with an equal volume of phenol/chloroform/isoamyl alcohol (25:24:1, pH8.0).

    The manufacturer of 5 PRIME has a detailed protocol on how to handle Phase Lock Gel®. Handle phenol/chloroform/isoamy alcohol in chemical hood.

  • 38.
    Centrifuge Phase Lock Gel® heavy tubes at 12,000 g for 5 minutes. Recover the aqueous upper phase to a new microcentrifuge tube.
  • 39.
    Add to samples 4 μl of 5 M NaCl, 2 μl of 5 μg/μl glycogen, and 3 volume (ca. 600 μl) of 100% ethanol. Invert tubes a few times to mix. Incubate samples at − 80 °C for 3+ hours.

    Use NaCl instead of NaOAc to prevent precipitation of SDS.

  • 40.
    Precipitate DNA by centrifuging at maximum speed at 4 °C for 30 minutes.
  • 41.
    Wash pellet twice with 1 ml of 75% ethanol. Dry the pellet.

    Dried pellet may become transparent.

  • 42.
    Resuspend DNA pellet in 50 μl H2O.
  • 43.
    Measure DNA concentration of input sample by Qubit or TapeStation to estimate the amount of DNA used in ChIP.

    Unlike input DNA, concentration of ChIPped DNA is often very low. ChIP of low abundance TFs from 3 grams of etiolated Arabidopsis seedlings may frequently enrich less than 5 ng of fragmented DNA.

SUPPORT PROTOCOL 1 (optional)

Optimization of ChIP-seq library preparation

Many suppliers such as Illumina and New England BioLabs already have detailed protocols to construct ChIP-seq libraries. To avoid redundancy this protocol focuses on how to optimize the size distribution of ChIP-enriched DNA. A tight size distribution of ChIP-enriched DNA is desired. Excluding long fragments improves cluster formation on Illumina sequencer. Tight size distribution of ChIP-enriched DNA also allows more accurate prediction of fragment length and peak calling by analysis packages. A few examples are shown in Figure 3.

Figure 3
Bead-based fractionation improves size distribution of sonicated DNA.

Materials

  • ChIP-enriched DNA, a few nanograms or more.
  • AMPure XP beads (Beckman Coulter)
  • 80% (v/v) ethanol, freshly prepared.

Equipment

  • DynaMag™-2 Magnet (Thermo Fisher Scientific)
  • Thermomixer (Eppendorf)
  • Note: use RNase- and DNase-free reagents, deionized distilled water, and low binding plasticware, and barrier pipette tips throughout the experiment.

Double size selection of ChIP-enriched DNA

  1. Add 37.5 μl (0.75 volume) of AMPure XP beads to 50 μl of ChIP-enriched DNA in a 1.5 ml microcentrifuge tube.

    The volume of AMPure beads can be increased or decreased for a smaller or larger DNA size cutoff, respectively. We recommend empirically testing DNA size fractionation within the range from 0.65 volume to 0.85 volume of beads.

  2. Vigorously vortex tube for 10 – 15 times. Incubate sample at room temperature for 10 minutes.
  3. Collect beads on magnet. Transfer supernatant to a fresh microcentrifuge tube.

    Supernatant contains DNA fragment smaller than 500 bp.

  4. Add 62.5 μl AMPure beads to supernatant.
  5. Vigorously vortex tube for 10 – 15 times. Incubate sample at room temperature for 10 minutes.
  6. Collect beads on magnet. Discard supernatant.
  7. Wash beads twice with freshly prepared 80% ethanol. Remove supernatant after each round of wash.
  8. Dry beads in a Thermomixer at 37 °C for 3-5 minutes.

    The purpose is to remove leftover ethanol from previous wash steps.

  9. Gently resuspend dried beads with 50 μl H2O. Vigorously vortex tube for 10 – 15 times. Incubate sample at room temperature for 5 minutes.
  10. Collect beads on magnet. Transfer supernatant containing the eluted DNA into a fresh microcentrifuge tube.
  11. Check the concentration and size of DNA on TapeStation (Figure 3D).

    If TapeStation is not available or if DNA concentration is too low, size-fractionate input sample from step 42 of the basic protocol and run a gel to infer the size of ChIP-enriched DNA. SYBRgold will interfere with DNA migration in the gel (Figure 3B vs Figure 3C). Therefore, we recommend running samples in ethidium bromide gel, or stain gel with SYBRgold after the run.

REAGENTS AND SOLUTIONS

  • Fixation buffer (1% (v/v) formaldehyde in 10 mM HEPES-NaOH).
  • PBS-BSA (1x PBS, 0.5 (w/v) bovine serum albumin), pre-chill before use.
  • Extraction buffer I (10 mM Tris-HCl pH 7.4, 0.4 M sucrose, β-mercaptoethanol 0.035%*, 1 mM PMSF*, 50 uM MG132*, 1x Roche cOmplete protease inhibitor, EDTA-free*), pre-chill before use.

    * Add fresh before experiments. Same for extraction buffers II and III, nuclei lysis buffer and ChIP dilution buffer. Handle β-mercaptoethanol in chemical hood.

  • Extraction buffer II (10 mM Tris-HCl pH 7.4, 0.25 M sucrose, 1% Triton X-100, 10 mM MgCl2, 1mM PMSF*, 50 uM MG132*, 1x Roche cOmplete protease inhibitor, EDTA-free*), pre-chill before use.
  • Extraction buffer III (10 mM Tris-HCl pH 7.4, 1.7 M sucrose, 0.15% Triton X-100, 2 mM MgCl2, 1 mM PMSF*, 50 uM MG132*, 1x Roche cOmplete protease inhibitor, EDTA-free*), pre-chill before use.
  • Nuclei lysis buffer (50 mM Tris-HCl pH 7.4, 0.5% sarkosyl, 100 mM NaCl, 2 mM EDTA, 1 mM PMSF*, 50 μM MG132*, 1x Roche cOmplete protease inhibitor, EDTA-free*), pre-chill before use.
  • ChIP dilution buffer (50 mM Tris-HCl pH 7.4, 1.25% Triton X-100, 100 mM NaCl, 2 mM EDTA, 1 mM PMSF*, 50 μM MG132*, 1x Roche cOmplete protease inhibitor, EDTA-free*), pre-chill before use.
  • Low salt wash buffer (50 mM Tris-HCl pH 7.4, 150 mM NaCl, 2 mM EDTA, 0.5% Triton X-100), pre-chill before use.
  • High salt wash buffer (50 mM Tris-HCl pH 7.4, 500 mM NaCl, 2 mM EDTA, 0.5% Triton X-100), pre-chill before use.
  • Final wash buffer (50 mM Tris-HCl pH 7.4, 50 mM NaCl, 2 mM EDTA), pre-chill before use.
  • Elution buffer (50 mM Tris-HCl pH 8.0, 10 mM EDTA, 1% SDS)

COMMENTARY

Background Information

Chromatin immunoprecipitation has been used to study protein-DNA interactions for almost 40 years (Jackson, 1978). The technique has evolved multiple times by the advances of DNA detection methods. In the early days, ChIPped DNA were examined by low throughput approaches such as Southern blot or PCR (Solomon et al., 1988; Orlando and Paro, 1993; Hecht et al., 1996). Understanding protein-DNA interactions at a genome-wide scale only became possible after the invention of ChIP-chip by Ren et al. (Ren et al., 2000) and Iyer et al. (Iyer et al., 2001), via the development of high-density DNA microarrays (Schena et al., 1995). ChIP-chip has facilitated numerous breakthroughs in biology, including elucidation of transcriptional regulatory networks in Saccharomyces cerevisiae (Lee et al., 2002) and characterization of DNA regulatory sequences in human (Consortium et al., 2007). On the other hand, ChIP-chip suffers from all the limitations of microarray technology. For instance, the design of array depends on a priori knowledge of genome sequence; high background signal makes it difficult to detect weak protein-DNA interactions; cross-hybridization is problematic, especially for highly homologous or repetitive regions; to obtain micrograms of DNA for hybridization, ChIPped DNA needs to be amplified by several orders of magnitude and the process may introduce biases. In 2007, Johnson et al. addressed these issues by combining ChIP with high-throughput DNA sequencing technologies (Johnson et al., 2007). With the help of sophisticated statistical analysis (Landt et al., 2012), ChIP-seq has become the most prevalent method to confidently identify protein-DNA interaction in vivo. More recently, Rhee et al. developed an improved ChIP-seq protocol called ChIP-exo (Rhee and Pugh, 2011). ChIP-exo utilizes 5'-to-3' exonuclease to trim DNA to the precise protein binding location. Exonuclease also cleans up ChIP background by digesting naked DNA. However, the extra steps associated with exonuclease digestion often result in greater sample loss and subsequently low library complexity. As an answer to this problem, He et al.'s ChIP-nexus protocol improves the ligation efficiency during library preparation, and tracks DNA over-amplification by unique, randomized barcodes (He et al., 2015).

Several factors contribute to the difficulty of performing ChIP-seq in plants. First, unlike animal cells in which nuclei can be extracted by mild detergents, extraction of plant nuclei usually requires vigorous physical disruption because of cell walls. This is arguably the step that causes the greatest sample loss in a ChIP procedure and prevents parallel handling of samples. Secondly, plant tissues often contain high level of phenolic compounds and polysaccharides, which may be problematic for PCR amplification during library preparation. Thirdly, there is limited selection of ChIP-grade antibodies in plants, and as a consequence researchers have to spend months to generate transgenic lines to express epitope-tagged proteins before ChIP-seq experiments can be carried out. Besides this protocol, several other laboratories have also published detailed procedures of ChIP-chip or ChIP-seq for Arabidopsis (Kaufmann et al., 2010; Reimer and Turck, 2010). We recommend readers combine knowledge of all protocols to decide the best practice.

Critical Parameters and Troubleshooting

ChIP-seq is a long procedure. We include multiple quality control steps in the protocols to ensure its success. Formaldehyde fixation is the first key step in the experiment. Both under- and over-fixation will result in inefficient ChIP of transcription factors. Therefore, we suggest using fresh formaldehyde, applying accurate control of cross-link time, and measuring the level of cross-link by western blot (step 5, basic protocol). Ideally, sonication should shear DNA to a consistent, tightly distributed size smaller than 500 bp whereas still preserving the protein-DNA complexes. The size distribution of sheared DNA can be examined by TapeStation (step 43, basic protocol), and can be further tightened by a bead-based size selection (support protocol 1, Figure 3). Antibody quality is crucial for successful ChIP. We suggest monitoring the effectiveness of IP by examining the presence of TF-of-interest after IP, wash, and final elution (step 33, basic protocol). Finally, over-amplification may bias library composition. We suggest determining the number of PCR cycles by real-time PCR to avoid over-amplification of ChIP-seq libraries (Figure 4).

Figure 4
Optimized library amplification determined by real-time PCR.

Anticipated Results

Using the basic protocol, 3 grams of etiolated seedlings is expected to yield more than 10 μg of input DNA. Depending on the TF and antibody, the ChIP-enriched DNA can be 5 ng or less. Light-grown samples especially flowers will usually yield more DNA. In most cases, twelve or fewer cyclers of PCR is sufficient to amplify enough DNA for sequencing.

Time Considerations

Once the samples are harvested and cross-linked, it takes 4 – 6 days to complete ChIP. The most time consuming steps in the procedure are IP (usually overnight, can be shortened to a few hours if the immunoprecipitated protein is very abundant), reverse cross-linking (6 hours to overnight), ethanol precipitation of reverse cross-linked DNA (3 hours to overnight), and adaptor ligation (overnight).

ACKNOWLEDGEMENT

We thank F.Turck, K.N.Chang, S.C.Huang, H.Qiao, U. Padmale, N. Krogan, M. Zander, M. Lewsey, and M. Urich for providing useful discussions on ChIP and DNA size fractionation. L.S. was supported by Salk Pioneer postdoc fellowship. This work was supported by grants from NSF (MCB-1024999 to J.R.E.). J.R.E. is an investigator of the Howard Hughes Medical Institute and the Gordon and Betty Moore Foundation.

LITERATURE CITED