An Array-based Comparative Genomic Hybridization Platform for Efficient Detection of Copy Number Variations in Fast Neutron-induced Medicago truncatula Mutants

Summary

This protocol provides experimental steps and information about reagents, equipment, and analysis tools for researchers who are interested in carrying out whole genome array-based comparative genomic hybridization (CGH) analysis of copy number variations in plants.

Abstract

Mutants are invaluable genetic resources for gene function studies. To generate mutant collections, three types of mutagens can be utilized, including biological such as T-DNA or transposon, chemical such as ethyl methanesulfonate (EMS), or physical such as ionization radiation. The type of mutation observed varies depending on the mutagen used. For ionization radiation induced mutants, mutations include deletion, duplication, or rearrangement. While T-DNA or transposon-based mutagenesis is limited to species that are susceptible to transformation, chemical or physical mutagenesis can be applied to a broad range of species. However, the characterization of mutations derived from chemical or physical mutagenesis traditionally relies on a map-based cloning approach, which is labor intensive and time consuming. Here, we show that a high-density genome array-based comparative genomic hybridization (aCGH) platform can be applied to efficiently detect and characterize copy number variations (CNVs) in mutants derived from fast neutron bombardment (FNB) mutagenesis in Medicago truncatula, a legume species. Whole genome sequence analysis shows that there are more than 50,000 genes or gene models in M. truncatula. At present, FNB-induced mutants in M. truncatula are derived from more than 150,000 M1 lines, representing invaluable genetic resources for functional studies of genes in the genome. The aCGH platform described here is an efficient tool for characterizing FNB-induced mutants in M. truncatula.

Introduction

Legumes (Fabaceae) are the third largest family of flowering plants, with many economically important species such as soybean (Glycine max) and alfalfa (Medicago sativa). Legume plants can interact with nitrogen-fixing soil bacteria, generally called Rhizobia to develop root nodules in which the atmospheric dinitrogen is reduced to ammonia for use by the host plant. As such, cultivation of legume crops requires little input of nitrogen fertilizers and thus contributes to sustainable agriculture. Legume crops produce leaves and seeds with high protein content, serving as excellent forage and grain crops. However, cultivated legume species generally have complex genome structures, making functional studies of genes that play key roles in legume-specific processes cumbersome. Medicago truncatula has been widely adopted as a model species for legume studies primarily because (1) it has a diploid genome with a relatively small haploid genome size (~550 Mbp); (2) plants can be stably transformed for gene functional studies; and (3) it is closely related to alfalfa (M. sativa), the queen of forages, and many other economically important crops for translational studies. Recently, the genome sequence of M. truncatula cv Jemalong A17 has been released^1,2. Annotation of the genome shows that there are more than 50,000 predicted genes or gene models in the genome. To determine the function of most of the genes in the M. truncatula genome is a challenging task. To facilitate functional studies of genes, a comprehensive collection of mutants in the range of over 150,000 M₁ lines has been generated using fast neutron bombardment (FNB) mutagenesis in M. truncatula cv Jemalong A17^3,4. Fast neutron, a high energy ionization mutagen, has been used in generating mutants in many plant species including Arabidopsis^5,6, rice (Oryza sativa)⁷, tomato (Solanum lycopersicum), soybean (Glycine soja; G. max)^8,9, barley (Hordeum vulgare), and Lotus japonicus¹⁰. A large portion of mutations derived from FNB mutagenesis are due to DNA deletions that range in size from a few base pairs to mega base pairs^9,11. Many phenotype-associated genes have been successfully identified and characterized^{4,12,13,14,15,16,17,18,19}. Previously, molecular cloning of the underlying genes from FNB mutants relied on a map-based approach, which is time consuming and limits the number of mutants to be characterized at the molecular level. Recently, several complimentary approaches including transcript-based methods, genome tiling array-based comparative genomic hybridization (CGH) for DNA copy number variation detection, and whole genome sequencing, have been employed to facilitate the characterization of deletion mutants in diverse organisms including animals and plants^{20,21,22,23,24,25,26,27,28,29,30,31}.

To facilitate the characterization of FNB mutants in M. truncatula, a whole-genome array-based comparative genomic hybridization (CGH) platform has been developed and validated. As reported in animal systems, the array-based CGH platform allows detection of copy number variations (CNVs) at the whole genome level in M. truncatula FNB mutants. Furthermore, lesions can be confirmed by PCR and deletion borders can be identified by sequencing. Overall, the array CGH platform is an efficient and effective tool in identifying lesions in M. truncatula FNB mutants. Here, the array CGH procedure and PCR characterization of deletion borders in an M. truncatula FNB mutant are illustrated.

The following protocol provides experimental steps and information about reagents, equipment and analysis tools for researchers who are interested in carrying out whole genome array-based comparative genomic hybridization (CGH) analysis of copy number variations in plants. As an example, Medicago truncatula FN6191 mutant was used to identify deletion regions and candidate genes associated with mutant phenotypes. M. truncatula FN6191 mutant, originally isolated from a fast neutron bombardment-induced deletion mutant collection³² (see Table of Materials), exhibited a hyper-nodulation phenotype after inoculation with the soil bacterium, Sihorhizobium meliloti Sm1021, in contrast to wild type plants.

Protocol

NOTE: Figure 1 shows the five steps for the array CGH protocol. They are: 1) Preparation of plant materials; 2) Isolation of high quality DNA samples; 3) Labeling and purification of DNA samples; 4) Hybridization, washing, and scanning of whole genome arrays; and 5) CGH data analysis. M. truncatula whole genome tiling arrays contain a total of 971,041 unique oligo probes targeting more than 50,000 genes or gene models in the genome (See Table of Materials). The unique probes are spaced approximately every 150 base pairs (bp) in exonic regions and 261 bp in intronic regions of the M. truncatula genome.

1. Preparation of Plant Materials

1. Scarify wild type (WT; M. truncatula cv Jemalong A17) and FN6191 mutant seeds with concentrated sulfuric acid for 8 min (see Table of Materials). Remove and discard sulfuric acid to a liquid waste container.
2. Rinse seeds with autoclaved deionized water three times.
3. Surface sterilize seeds with 20% bleach solution for 10 min.
4. Rinse seeds with autoclaved deionized water three times.
5. Store seeds at 4 °C for 3 days. After vernalization, transfer and germinate seeds in soil for one week in a growth chamber with 16 h/8 h light/dark cycle and 150 µE/m²/s light intensity.
6. Transplant seedlings to 1 gallon pots and grow them in a greenhouse with 16 h/8 h light/dark cycle, 150 µE/m²/s light intensity for three to four weeks. Collect young leaf samples from the plants for DNA isolation.

2. Isolation of High Quality DNA Samples

1. Take 1 g of young leaf tissue from a single plant and freeze it in liquid nitrogen immediately.
2. Grind frozen leaf tissue in liquid nitrogen with a mortar and pestle to fine powder.
3. Isolate genomic DNA with a DNA isolation kit (see Table of Materials).
4. Resuspend purified DNA samples in 50 µL of autoclaved double deionized H₂O (ddH₂O) or 1x TE buffer (10 mM Tris-HCl, 1 mM EDTA, pH 7.0).
5. Measure DNA concentrations with a spectrophotometer (see Table of Materials).
6. Evaluate 260 nm/280 nm and 260 nm/230 nm ratios of DNA samples.
   NOTE: High quality DNA samples should have a 260 nm/280 nm ratio ≥1.8 and a 260 nm/230 nm ratio ≥2.0.
7. Further evaluate DNA samples by running them in 1% agarose gels.
   NOTE: High quality DNA samples should have high molecular weight and should not be contaminated by RNA. Ideally, the final concentration of DNA samples should be ~150 ng/µL.

3. DNA Labeling and Purification

1. Dilute 1 µg of genomic DNA samples with ddH₂O to a final volume of 20 µL in a 500 µL centrifuge tube.
2. Add 5 µL of a random primer (see Table of Materials) to the tube.
3. Vortex the tube briefly and put it into a thermocycler to incubate and denature DNA at 98 °C for 10 min.
4. Take the tube out of the thermocycler and put it immediately on ice water to chill for 5 min.
5. Prepare two labeling mixes as shown in Table 1.
6. Add 25 µL of labeling mixes 1 and 2 to reference DNA and sample DNA tubes, respectively.
7. Mix the labeling mix and DNA by pipetting three times and spin the tubes briefly.
8. Put the tubes into a thermocycler to incubate at 37 °C for 2 h and then at 65 °C for 20 min to inactivate the exo-Klenow enzyme (see Table of Materials).
9. Add 430 µL of 1x TE buffer to each tube.
10. Mix the contents and spin the tubes briefly.
11. Transfer the solution in the tube into a purification column with a 2 mL collection tube (see Table of Materials).
12. Centrifuge the purification column at 14,000 x g for 10 min.
13. Discard the pass through solution.
14. Add 480 µL 1x TE buffer to each column.
15. Centrifuge the column again at 14,000 x g for 10 min.
16. Discard the pass-through solution.
17. Transfer the labeled DNA sample that remains at the bottom portion of the column to a new centrifuge tube.
18. Measure the volume of the labeled DNA sample using a pipette and bring the final volume to 80 µL with 1x TE buffer.
19. Measure the concentration of the labeled DNA sample using a spectrophotometer (see Table of Materials).
20. Mix an equal amount of labeled sample DNA and reference DNA and bring the final volume to 160 µL with 1x TE buffer.
21. Add 256 µL of 2x hybridization buffer, 50 µL of human Cot-1 DNA and 50 µL of 10× aCGH blocking agent (see Table of Materials) and mix them well by pipetting for three times.
22. Spin the tubes briefly and put them into a thermocycler to incubate at 98 °C for 10 min.
23. Incubate the tubes at 37 °C for 20 min.

4. Hybridization, Washing, and Scanning of the Genome Arrays

1. Pre-warm the hybridization oven to 67 °C at least 4 h before hybridization starts.
2. Place a gasket slide onto the base of a hybridization chamber.
3. Load 490 µL of hybridization solution from Step 3.23 onto the gasket slide.
4. Cover the gasket slide with a Medicago truncatula genome microarray chip to form a hybridization chamber.
5. Put on the hybridization chamber covers and tighten the chamber firmly with clamps.
6. Put the assembled hybridization chamber into the hybridization oven and incubate at 67 °C for 40 - 48 h.
7. Prepare three slide washing dishes: two with 250 mL washing buffer I kept at room temperature, and one with washing buffer II kept at 37 °C.
8. Prepare two slide washing jars: one with 70 mL acetonitrile and one with 70 mL stabilization and drying solution; keep both in a fume hood at room temperature (see Table of Materials).
9. Remove the hybridization chamber from the oven and remove the chamber clamps and cover.
10. Move the microarray chip and gasket slide to a slide washing dish with wash buffer I and separate the chip from the cover slide.
11. Move the microarray chip to the second slide washing dish with wash buffer I and gently circulate the solution with a stir bar on a magnetic stir plate at room temperature for 5 min.
12. Transfer the microarray chip to the third slide washing dish with pre-warmed washing buffer II and wash it gently with a stir bar on a magnetic stir plate at 37 °C for 1 min.
13. Transfer the microarray chip to a slide washing jar with acetonitrile in a fume hood and wash it for 30 s.
14. Transfer the microarray chip to a slide washing jar with stabilization and drying solution and wash it for 30 s.
15. Take out the chip slowly and dry it for 1 min.
16. Scan the microarray chip using a scanner (see Table of Materials) under 2 µm resolution. Set parameters to save multiple images and whole slide scans. Set channels 1 and 2 gains to 100% and cancel the auto gain.

5. CGH Data Analysis

1. Extract the Cy3 and Cy5 fluorescence intensities for each probe on the array using scanning software (see Table of Materials).
2. Use segmentation analysis of the software to detect copy number variations (CNVs) between the sample DNA and reference DNA. Use a threshold of the segment mean log₂ sample/reference ratio ± 2.5 standard deviations for determining CNVs.
3. Inspect the distribution of log₂ sample/reference ratios across the whole genome using a signal mapping software (see Table of Materials).

Results

Figure 2 shows the distribution of normalized log₂ ratios of mutant versus WT signals across the whole genome. Analysis of CGH data revealed an approximate 22 kb deletion on chromosome 4 that encompasses the entire SUNN gene³³ and several other annotated genes in FN6191 mutant (Figure 2, Figure 3). The candidate deleted region was covered by 73 consecutive ...

Discussion

We have developed an array-based CGH platform for the detection and characterization of fast neutron bombardment (FNB)-induced mutants in M. truncatula cv. Jemalong A17. To demonstrate the use of the array CGH method in detecting gene mutations, we performed aCGH analysis of the mutant FN6191, which exhibited a hyper-nodulation phenotype in contrast to wild type plants, when inoculated with S. meliloti Sm1021. For segmentation analysis, a segment was deemed significant if the log₂ ratio mean ...

Materials

Name	Company	Catalog Number	Comments
Medicago truncatula genome array, 1 x 1 M	Agilent	G4123A
Medicago truncatula FN6191 (mutant)	In house	FN6191
Medicago truncatula Jemalong A17 (reference)	In house	A17
Sulfuric acid	Sigma-Aldrich	320501
DNeasy Plant Mini Kit	Qiagen	69104
Nanodrop Spectrophotometer	Thermo Scientific	1000D
SureTag DNA Labeling Kit	Agilent	5190-3400
Random primer	Agilent	5190-3399
Acetonitrile	Sigma-Aldrich	271004-1L
Thermocycler	MJ research	PTC-200
Centrifuge	Labnet international Inc	Spectrafuge 24D
Stabilization and Drying Solution	Agilent	5185-5979
Oligo aCGH/ChIP-on-chip Hybridization Kit	Agilent	5188-5380
Hybridization Chamber gasket slides	Agilent	G2505
Human Cot-1 DNA	Agilent	5190-3393
Oligo aCGH/ChIP-on-chip Wash Buffer 1 and 2	Agilent	5188-5221
Hybridization Chamber, stainless	Agilent	G2534A
Hybridization oven	Agilent	G2545A
Purification Columns	Agilent	5190-3391
Laser scanner	Roche	MS200
NimbleScan 2.6	Roche Nimblegen	5225035001
Signal Map 1.9	Roche Nimblegen	Signalmap1.9