PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plantsigLink to Publisher's site
 
Plant Signal Behav. 2013 January 1; 8(1): e22534.
Published online 2012 October 26. doi:  10.4161/psb.22534
PMCID: PMC3745560

Identification of an EMS-induced causal mutation in a gene required for boron-mediated root development by low-coverage genome re-sequencing in Arabidopsis

Abstract

Next-generation sequencing (NGS) technologies enable the rapid production of an enormous quantity of sequence data. These powerful new technologies allow the identification of mutations by whole-genome sequencing. However, most reported NGS-based mapping methods, which are based on bulked segregant analysis, are costly and laborious. To address these limitations, we designed a versatile NGS-based mapping method that consists of a combination of low- to medium-coverage multiplex SOLiD (Sequencing by Oligonucleotide Ligation and Detection) and classical genetic rough mapping. Using only low to medium coverage reduces the SOLiD sequencing costs and, since just 10 to 20 mutant F2 plants are required for rough mapping, the operation is simple enough to handle in a laboratory with limited space and funding. As a proof of principle, we successfully applied this method to identify the CTR1, which is involved in boron-mediated root development, from among a population of high boron requiring Arabidopsis thaliana mutants. Our work demonstrates that this NGS-based mapping method is a moderately priced and versatile method that can readily be applied to other model organisms.

Keywords: next-generation sequencing, low- to medium-coverage sequencing, mutant mapping, Arabidopsis thaliana, boron

Introduction

Next-generation sequencing (NGS) technologies have superseded conventional Sanger sequencing by capillary electrophoresis in many instances of large-scale analysis, such as de novo sequencing, transcriptomics, DNA methylation analysis, metagenomics and population genetic studies.1-4 These new technologies also provide a powerful tool for identifying ethyl methanesulfonate (EMS)- and N-ethyl-N-nitrosourea (ENU)-induced the mutations by whole-genome re-sequencing. For example, mutations underlying an aberrant neuronal phenotype in Caenorhabditis elegans were identified by re-sequencing of the mutant genomes.5 In another case, mutants of the ethanol producing yeast Pichia stipitis were sequenced to identify mutations that facilitated efficient fermentation.6 The ability to identify mutations via NGS technologies has greatly reduced the amount of time needed for conventional map-based cloning.

In plant research, as in research in a variety of model organisms, these NGS technologies have been successfully applied to identify the mutations underlying phenotypes of interest. Schneeberger, et al.7 developed a method called SHOREmap that uses an Illumina Genome Analyzer (GA) to identify causative mutations of A. thaliana.7 By sequencing a genomic DNA sample prepared from a pool of 500 mutant F2 plants that were obtained by crossing a mutant with another wild-type accession, they have identified a causative mutation. This method is based on an approach termed bulked segregant analysis.8,9 For example, a mutant in the Columbia (Col-0) background was crossed to the polymorphic Landsberg erecta (Ler) accession, followed by selfing of the F1 progeny to generate an F2 population in which the mutant phenotype segregates according to Mendelian rules. Then, F2 plants exhibiting the mutant phenotype were pooled and the bulk segregant pool was subjected to deep sequencing. The causative region was confined to a small genomic region by analyzing the index of enrichment of Col-0-type single nucleotide polymorphisms (SNPs) and, finally, causal mutations were identified from the sequence data obtained by deep sequencing analysis.7 Recently, Austin, et al.10 developed a Next-Gen Mapping (NGM) method by modifying the above-mentioned SHOREmap technique, and succeeded in identifying three genes involved in cell wall biology.10 Moreover, a spontaneous mutation in a nonreference A. thaliana accession and EMS-induced mutations in a nonreference accession background were successfully identified using deep sequencing.11,12 These modifications of bulked segregant analysis are extremely useful for identifying mutations in A. thaliana and can also be applied in crops and other organisms with fully sequenced genomes. However, deep sequencing remains expensive and laborious, as approximately 100 or more mutant F2 plants are required for this type of bulked segregant analysis.

To address these problems, we designed a versatile NGS-based mapping method that incorporates SOLiD (Sequencing by Oligonucleotide Ligation and Detection). This mapping method is based on a combination of low- to medium-coverage SOLiD13 and classical genetic rough mapping. Sequencing at just low to medium coverage reduced costs. Furthermore, since rough mapping required only 10 to 20 F2 plants with the mutant phenotype, experiments using this strategy do not require a lot of space. Using this method, we rapidly identified CTR1, which is involved in boron-mediated growth. Here, we describe this NGS-based mapping method and discuss its applications.

Results

Screening for high boron requiring mutants

Boron is an essential nutrient for plants and boron deficiency is a major cause of reduced crop production.14 Boron maintains the structure and function of the plant cell wall by cross-linking the pectic polysaccharide rhamnogalacturonan (RG) II.15 Several boron transporters that are upregulated under boron limitation have been identified in A. thaliana16-18 and their polar localization and degradation through trafficking pathways in plant cells have been demonstrated.19,20 Although significant progress has been made in our understanding of boron transport mechanisms, the precise role of boron in plant growth and development remains unclear.

To obtain insight into boron function in plants, the EMS-mutagenized M2 seeds of A. thaliana were screened for mutants that required more boron than the wild type for root elongation. Approximately 20,000 seeds were sown onto normal medium (30 µM B) and short-root plants were transferred to medium containing 1 mM boron after 7 d. After growth on high boron medium for 7 d, plants that exhibited increased root elongation at 1 mM boron were selected. From this screening, we isolated 13 mutants. We named one of these mutants b26–6/ctr1–16, as it is allelic to the ctr1 mutants described later (Fig. 1A).

figure psb-8-e22534-g1
Figure 1. Identification and characterization of the b26–6/ctr1–16 mutants. (A) The seeds were sown on MGRL medium containing 0.3 µM, 30 µM and 1 mM boron and grown for 2 weeks. (B) Identification of the causal ...

Rough mapping

The b26–6/ctr1–16 mutant in the Col-0 background was crossed with Ler wild-type plants for rough mapping. The F2 population segregated into wild type and mutant type at a ratio of 3:1, indicating that the mutant phenotype is caused by a single recessive mutation. Genomic DNA was isolated from 12 F2 plants that exhibited the mutant phenotype and the mutation was assigned to a chromosome using simple sequence length polymorphism (SSLP) markers F15A17 and T32M21. A candidate region with the mutation was rough mapped to between 0.70 Mb and 1.26 Mb on chromosome 5, a region that spanned 175 putative genes annotated in TAIR9 (Fig. 1B and Table 2).

Table thumbnail
Table 2. EMS treatment conditions and SNP filtering in the b26–6/ctr1–16 mutant

SOLiD sequencing

To identify point mutations, we sequenced the genomic DNA of the b26–6/ctr1–16 mutant by SOLiD. We constructed sequence libraries from the b26–6/ctr1–16 mutant and seven other mutants derived from Lehle Seeds using the SOLiD barcoding system to distinguish the eight samples (Fig. 2). The 8-plex libraries were sequenced on a single SOLiD slide. In total, 378.4 M reads were obtained, of which 58.4 M were assigned to the b26–6/ctr1–16 mutant library (see Table 1 for details). Of all the b26–6/ctr1–16 mutant library reads, 73.2% were mapped to the TAIR9 release of the A. thaliana Col-0 genome. The median value of per-base sequence depth was 10 and the genome coverage was 91.8% (Table 1 and Fig. S1).

figure psb-8-e22534-g2
Figure 2. Scheme of the method used to identify mutations described in this manuscript. This method is based on a combination of two approaches: low- (< 5 × per site per individual, on average) or medium- (5–20 × ...
Table thumbnail
Table 1. Mapping of SOLiD sequence tags to the Arabidopsis genome (TAIR9 release)

SNP detection

SNPs refer to sites that differ from the TAIR9 release of the A. thaliana Col-0 genome. We used diBayes for SNP calling with a low stringency setting. Although this setting may increase the risk of reporting false positives, considering that the sequence depth was relatively low (median = 10), we chose this aggressive but sensitive SNP calling strategy to avoid overlooking true SNPs of low sequence coverage. With this strategy, the diBayes program listed 2162 homozygous SNP sites (Table 2).

We assessed the pattern of genome-wide SNPs by plotting the SNP frequencies using a bin size of 1 Mb (Fig. S2A). Almost all SNPs were distributed evenly throughout the genome at an average of 18 SNPs per 1 Mb, with a few regions being enriched in SNPs. This SNP enrichment may be due to variations that were already present in the parental line, rather than to an EMS mutagenesis bias, as a substantial number of SNPs were shared in these regions among eight Lehle Seeds libraries sequenced in this experiment (Fig. S2B). By subtracting shared SNPs, which were regarded as background, the number of SNPs was reduced to 2046. We then extracted SNPs that involved G-to-A (or complementary C-to-T) nucleotide substitutions, which are typical of EMS mutagenesis,21,22 and obtained 462 SNPs. We classified the SNPs according to gene models: 179 occurred in exons and 0 in splice donor/acceptor sites. In total, 177 genes were found to be mutagenized in the genome of the b26–6/ctr1–16 mutant. The b26–6/ctr1–16 mutant was isolated from Lehle Seeds, which are widely used by Arabidopsis researchers. Based on this experiment and additional unpublished ones (Kamiya T and Fujiwara T, unpublished data), an average of 450 mutations was found to exist in the genomes of the Lehle Seeds mutants, approximately 200 mutations of which were exonic. This information can be used in the design of experiments to isolate A. thaliana mutants.

Our re-sequencing method detected three SNPs within the 0.56-Mb candidate genomic region on chromosome 5, which was narrowed down by genetic rough mapping and found to contain 175 genes (Fig. 1B). Two of the three mutations resulted in non-synonymous amino acid changes (Table 3).

Table thumbnail
Table 3. Mutations that result in amino acid substitutions in the mapped region

Validation of candidate SNPs

Genome re-sequencing by SOLiD identified three candidate SNPs for the b26–6/ctr1–16 mutant in the region narrowed down by rough mapping. Two of the three mutations resulted in an amino acid substitution in At5g02990 and At5g03730, respectively. At5g03730 is annotated as CONSTITUTIVE TRIPLE RESPONSE1 (CTR1),23 and the root phenotype of ctr1 mutants, especially an abundance of hairy roots, was quite similar to that of the mutant identified in this study, b26–6/ctr1–16.24,25 b26–6/ctr1–16 has a point mutation that changes the 610th codon from CTT to TTT, which results in an amino acid substitution of Leu to Phe. Leu610 is highly conserved among species and localizes to the kinase domain of CTR1.

To show that CTR1 is the causal gene of b26–6/ctr1–16, the boron response of two ctr1 alleles, ctr1–1 and sis1/ctr1–12, was examined (Fig. 1A).23,25 Both alleles displayed a similar phenotype to b26–6/ctr1–16 under normal growth conditions and, like b26–6/ctr1–16, root elongation was restored to normal rates by treatment with 1 mM boron. These results indicate that the causal gene of b26–6/ctr1–16 is indeed CTR1.

Discussion

In this paper, we demonstrate that a versatile NGS-based mapping method, which combines low- to medium-coverage sequencing by SOLiD with classical genetic rough mapping, may be used to identify causal candidate mutations in A. thaliana. This method is superior to other NGS-based mapping methods in terms of time and cost. Recently, a number of NGS-based mapping methods using bulked segregant analysis have been proposed as means to identify causative mutations.5-7,10-12 Although these methods are undoubtedly effective, they are all laborious and time consuming, as they involve the preparation of approximately 100 or more mutant F2 plants for bulked segregant analysis. This is particularly problematic when dealing with mutants that have a phenotype that is difficult to propagate or score. In addition, when a mutant is crossed to the other accession for mapping, many SNPs are mixed. These background variations occasionally make it difficult to distinguish between the mutant and wild-type phenotype. To avoid this problem, we separated genetic mapping and NGS analysis using SOLiD with low- to medium-coverage. In our method, rough mapping requires only 10 to 20 mutant F2 plants, which can readily be produced in a small laboratory space. In fact, rough mapping using 12 mutant F2 plants was sufficient to identify the causal gene in this study. Based on this experiment and additional unpublished ones (Kamiya T and Fujiwara T, unpublished data), we have identified mutations in intervals of up to 3 Mb that were narrowed down by rough mapping by using approximately ten mutant F2 segregants.

Furthermore, since low- to medium-coverage sequencing can be used to identify a mutation, our method facilitates the simultaneous analysis of multiple mutants (a slide can be shared by 16~20 mutants under the current conditions) by using a multiplexed barcoded samples in a single sequencing run. Thus, we can handle multiple mutants in a single run, and/or multiple laboratories can share slides of barcoded DNA. This method reduces the cost of each sample.

Applicability to other mutants

We applied our NGS-based mapping technique to mutants other than b26–6/ctr1–16, i.e., #6–10, #7–3 and #9–2, to identify the causative mutations, and succeeded in detecting a small number of testable candidate mutations. In A. thaliana, the CLAVATA3 (CLV3) peptide negatively regulates the size of the shoot meristem via the CLV1, CLV2, CRN and RPK2 receptors.26-28 Moreover, overexpression of CLV3 or some other CLV3/ESR-like (CLE) gene, and CLE peptide treatments reduce the size of the root meristem and the length of the roots, suggesting that the CLE signaling pathway also acts in the root.29 To identify new factors that regulate root meristem development, we performed mutant screening as an index of the CLE peptide-resistance phenotype of root formation in A. thaliana, and isolated these #6–10, #7–3 and #9–2 mutants. In these three cases, we detected 9, 6 and 13 SNPs in 3.1-, 2.7- and 5.8-Mb regions, using 17, 15 and 19 F2 mapping populations, respectively (Yamada M, Tabata R and Sawa S, unpublished data; Table. 4).

We also used this method to identify the causal mutations of mutants derived from non-Col-0 backgrounds. In each case, a small number of candidate mutations could be detected. Therefore, our method appears to be suitable for detecting mutations in a non-reference A. thaliana accession background (data not shown).

Moreover, we have already used our method to identify strong mutations in mutagenized rice plants (data not shown). Rice has a complex genome that is three times larger than that of A. thaliana,30 indicating that our method may readily be applied to various other model plants that can be crossed to a mapping line. We will soon provide a platform to identify the mutations responsible for phenotypes of interest in rice and other plants, as we have already done for those in A. thaliana.

Effective size of mapping population

Examining a large number of F2 plants for traditional map-based cloning and bulked segregant analysis is generally time-consuming and labor-intensive. Each NGS-based mapping method using bulked segregant analysis proposed by Uchida, et al.12 and Austin, et al.10 encouraged the use of 80 and at least 50 mutant F2 plants, respectively. Because these NGS-based mapping methods are excellent for identifying causal mutations, they are useful in cases in which a sufficient number of mutant F2 plants are readily collectible. However, in this study, we have shown that only 10 to 20 mutant F2 plants are needed to identify causal mutations using our NGS-based mapping method, as indicated in a number of examples. Our method is particularly effective in situations in which it is difficult to propagate or score mutant F2 plants.

Robustness with respect to location of mutation

The causal mutation of the b26–6/ctr1–16 mutant is located in a low-recombination region near the telomere of chromosome 5. Thus, we have identified a causal mutation even in a region of the genome that is difficult to sequence, i.e., near the telomere. Moreover, the mutations underlying #6–10, #7–3 and #9–2, were located on various chromosomes and regions, demonstrating that the NGS-based mapping method presented here has utility across the A. thaliana genome; however, it remains unclear whether this method is applicable to centromeric regions.

The relationship between CTR1 and boron

The short root phenotype of b26–6/ctr1–16 mutant plants was recovered by supplementation of high boron (Fig. 1A). Considering that CTR1 is a suppressor of ethylene signaling,23 of which loss-of-function mutant shows constitutive active ethylene response, there are two possible reasons for the recovery. One is that high boron inhibits constitutive ethylene signaling in b26–6/ctr1–16 mutant. It has been reported that low boron induces ethylene responsive gene and ethylene responsive reporter gene expression in tobacco BY-2 cell and A. thaliana roots, respectively.31,32 In other words, boron inhibits the induction of ethylene responsive gene. Although there is no direct evidence that high boron inhibits ethylene signaling, it is possible that ethylene signaling downstream of CTR1 is suppressed by high boron leading to root elongation. The other possibility is that high boron rescues the defect in RG-II cross-linking with borate, which is the only reported direct function of boron, in the mutant.33,34 In ctl1/arm mutant, which produces increased amount of ethylene,35 cell wall composition of root is changed; Increase in arabinose and galactose: Decrease in xyloglucan (XyG) and mannose (Man), and methylesterification of homogalacturonans. Those changes in cell wall components may affect the structure of RG-II. For example, Man containing nucleotide sugar, GDP-D-Man, is required for the synthesis of GDP-L-fucose, which is a component of RG-II. The phenotype of mur1 mutant, which has defect in the gene required for GDP-L-fucose from GDP-D-Man, was recovered by high boron like b26–6/ctr1–16.36-38 Furthermore, the phenotype of mur3 and hsr8/mur4, which have defect in RG-II composition, was also recovered by high boron.38 Those result suggest that short root in b26–6/ctr1–16 is caused by the defect in boron cross-linking in RG-II.

Conclusion

In this study, we developed a versatile NGS-based mapping method that combined low- to medium-coverage sequencing by SOLiD with classical genetic rough mapping and successfully applied this method to identify the CTR1, which is involved in boron-mediated root development. In fact, we demonstrated that this is a cost effective and versatile method for identifying EMS-induced causal mutations. We hope that this approach will be a useful tool for many Arabidopsis researchers and will be applied to studies on other model plants.

Materials and Methods

Plant materials and growth conditions

EMS-mutagenized Columbia (Col-0 gl1–1) seeds were purchased from Lehle Seeds. ctr1–1 (CS8057) and sis1/ctr1–14 (CS3874) were obtained from the Arabidopsis Biological Resource Center. The seeds were surface-sterilized and sowed on MGRL growth medium39 solidified with 1.2% gellan gum (Wako) supplemented with 1% sucrose.

Mutant selection and genetic mapping

Approximately 20,000 plants were sowed onto normal MGRL medium (30 µM boron). After incubation for 2 d at 4°C, the plates were placed vertically and the plants were grown at 22°C under a 16 h light/8 h dark photoperiod. After 7 d, the plants with short roots were transferred to medium containing 1 mM boron. After 7 d, the plants with partially recovered roots were selected and the phenotype was confirmed in the M3 generation.

For mapping, the b26–6/ctr1–16 mutant was crossed with wild-type Landsberg erecta (Ler) and F2 seeds were obtained. Genomic DNA was isolated from F2 plants that exhibited the mutant phenotype and the gene was mapped using simple sequence length polymorphism (SSLP) markers and the following primer pairs: F15A17-F (5′-ACCCAAACTTGGCTCACAAC-3′), F15A17-R (5′-CAAAATCATCTCCCCTTGGA-3′) and T32M21-F (5′-CAAACGTAAAACATAAAAGAGAACCA-3′), T32M21-R (5′-TCCGTTGCTTAGAACATTTGC-3′).

Library construction for whole-genome sequencing

Genomic DNA of the M4 generation of the b26–6/ctr1–16 homozygous mutant was isolated using the DNeasy Plant Mini Kit (QIAGEN). The genomic DNA was sheared into 100- to 150-bp fragments using the Covaris S2 system (Covaris, Inc.) in a 120-μL reaction containing 10 mM TE buffer in a Covaris microTube and the following program: 20% duty cycle, 5 intensity and 200 cycles per burst for 60 sec at 5°C. Fragments were arranged into a SOLiD barcoded fragment library using the SOLiD Fragment Library Construction Kit (Life Technologies).

SOLiD sequencing

The fragment libraries were amplified by emulsion PCR (ePCR) at a library concentration of 0.5 pM. After ePCR, the beads were modified at the 3′ terminal and deposited onto a SOLiD sequencing slide, according to the manufacturer’s instructions. The libraries were sequenced to 50 base pairs using an Applied Biosystems SOLiD 3 Plus System.

Data analysis of sequencing reads

Color space reads were mapped to the Arabidopsis thaliana genome reference TAIR9 (TAIR_chr_all.fas) using BioScope 1.3 software (Life Technologies) with default parameters. Single nucleotide polymorphisms (SNPs) were called using a diBayes SNP caller, a component of BioScope 1.3, which implements a Bayesian algorithm that includes color space error detection. diBayes was executed with a parameter setting defined as “low call stringency” using the default parameters (major parameters: dibayes.reads.min.mapping.qv = 8; dibayes.reads.only.uniquely.mapped.allow = no; dibayes.snp.both.strands = no; dibayes.snps.min.base.qv = 26; dibayes.hom.min.nonref.start.pos = 1). An Integrative Genomics Viewer (IGV) http://www.broadinstitute.org/igv/)40 was used to visualize the mapped reads and called SNPs along with gene models. To predict the functional impact of the SNPs, we categorized them into coding (synonymous or non-synonymous), intronic, intergenic and splicing site using custom scripts. Gene models and the annotations were based on TAIR9.

NOTE: Author, please cite Table 4 in text.

Table thumbnail
Table 4. EMS treatment conditions and SNP filtering in the three A. thaliana mutants

Supplementary Material

Additional material

Acknowledgments

We thank Yuko Kawara and Yayoi Inui-Tsujimoto for technical assistance and Dr. N. Uchida for his comments on the manuscript.

Funding

This study was performed under the NIBB cooperative research program (11–103). This work was supported by: Grant-in Aid for Creative Scientific Research; Grant-in-Aid for Young Scientists S to S.S. (19677001) from Japan Society of the Promotion of Science; Grant-in-Aid for Scientific Research for Priority Areas from the Ministry of Education, Culture, Sports, Science and Technology, Japan to S.S. (221S0002, 23119517, 23012034, 24114001 and 24114009), a Grant-in-Aid for Scientific Research S (to T.F.) and a Grant-in-Aid for Scientific Research on Innovative Areas (to T.F.) from the Ministry of Education, Culture, Sports, Science and Technology.

Author Contributions

Conceived and designed the experiments: R.T., T.K., T.F. and S.S. Performed the experiments: R.T., T.K., K.Y. and M.Y. Analyzed the data: S.S. and K.Y. Contributed reagents/materials/analysis tools: M.H., T.F. and S.S. Wrote the paper: R.T., T.K. and S.S.

Disclosure of Potential Conflicts of Interest

No potential conflicts of interest were disclosed.

Footnotes

References

1. Liti G, Carter DM, Moses AM, Warringer J, Parts L, James SA, et al. Population genomics of domestic and wild yeasts. Nature. 2009;458:337–41. doi: 10.1038/nature07743. [PMC free article] [PubMed] [Cross Ref]
2. MacLean D, Jones JD, Studholme DJ. Application of ‘next-generation’ sequencing technologies to microbial genetics. Nat Rev Microbiol. 2009;7:287–96. [PubMed]
3. Li R, Fan W, Tian G, Zhu H, He L, Cai J, et al. The sequence and de novo assembly of the giant panda genome. Nature. 2010;463:311–7. doi: 10.1038/nature08696. [PubMed] [Cross Ref]
4. Metzker ML. Sequencing technologies - the next generation. Nat Rev Genet. 2010;11:31–46. doi: 10.1038/nrg2626. [PubMed] [Cross Ref]
5. Sarin S, Prabhu S, O’Meara MM, Pe’er I, Hobert O. Caenorhabditis elegans mutant allele identification by whole-genome sequencing. Nat Methods. 2008;5:865–7. doi: 10.1038/nmeth.1249. [PMC free article] [PubMed] [Cross Ref]
6. Smith DR, Quinlan AR, Peckham HE, Makowsky K, Tao W, Woolf B, et al. Rapid whole-genome mutational profiling using next-generation sequencing technologies. Genome Res. 2008;18:1638–42. doi: 10.1101/gr.077776.108. [PubMed] [Cross Ref]
7. Schneeberger K, Ossowski S, Lanz C, Juul T, Petersen AH, Nielsen KL, et al. SHOREmap: simultaneous mapping and mutation identification by deep sequencing. Nat Methods. 2009;6:550–1. doi: 10.1038/nmeth0809-550. [PubMed] [Cross Ref]
8. Giovannoni JJ, Wing RA, Ganal MW, Tanksley SD. Isolation of molecular markers from specific chromosomal intervals using DNA pools from existing mapping populations. Nucleic Acids Res. 1991;19:6553–8. doi: 10.1093/nar/19.23.6553. [PMC free article] [PubMed] [Cross Ref]
9. Michelmore RW, Paran I, Kesseli RV. Identification of markers linked to disease-resistance genes by bulked segregant analysis: a rapid method to detect markers in specific genomic regions by using segregating populations. Proc Natl Acad Sci USA. 1991;88:9828–32. doi: 10.1073/pnas.88.21.9828. [PubMed] [Cross Ref]
10. Austin RS, Vidaurre D, Stamatiou G, Breit R, Provart NJ, Bonetta D, et al. Next-generation mapping of Arabidopsis genes. Plant J. 2011;67:715–25. doi: 10.1111/j.1365-313X.2011.04619.x. [PubMed] [Cross Ref]
11. Laitinen RA, Schneeberger K, Jelly NS, Ossowski S, Weigel D. Identification of a spontaneous frame shift mutation in a nonreference Arabidopsis accession using whole genome sequencing. Plant Physiol. 2010;153:652–4. doi: 10.1104/pp.110.156448. [PubMed] [Cross Ref]
12. Uchida N, Sakamoto T, Kurata T, Tasaka M. Identification of EMS-induced causal mutations in a non-reference Arabidopsis thaliana accession by whole genome sequencing. Plant Cell Physiol. 2011;52:716–22. doi: 10.1093/pcp/pcr029. [PubMed] [Cross Ref]
13. Nielsen R, Paul JS, Albrechtsen A, Song YS. Genotype and SNP calling from next-generation sequencing data. Nat Rev Genet. 2011;12:443–51. doi: 10.1038/nrg2986. [PMC free article] [PubMed] [Cross Ref]
14. Marschner H. Mineral Nutrition of Higher Plants. San Diego: Academic Press, 1995.
15. O’Neill MA, Ishii T, Albersheim P, Darvill AG. Rhamnogalacturonan II: structure and function of a borate cross-linked cell wall pectic polysaccharide. Annu Rev Plant Biol. 2004;55:109–39. doi: 10.1146/annurev.arplant.55.031903.141750. [PubMed] [Cross Ref]
16. Takano J, Noguchi K, Yasumori M, Kobayashi M, Gajdos Z, Miwa K, et al. Arabidopsis boron transporter for xylem loading. Nature. 2002;420:337–40. doi: 10.1038/nature01139. [PubMed] [Cross Ref]
17. Takano J, Wada M, Ludewig U, Schaaf G, von Wirén N, Fujiwara T. The Arabidopsis major intrinsic protein NIP5;1 is essential for efficient boron uptake and plant development under boron limitation. Plant Cell. 2006;18:1498–509. doi: 10.1105/tpc.106.041640. [PubMed] [Cross Ref]
18. Takano J, Miwa K, Fujiwara T. Boron transport mechanisms: collaboration of channels and transporters. Trends Plant Sci. 2008;13:451–7. doi: 10.1016/j.tplants.2008.05.007. [PubMed] [Cross Ref]
19. Takano J, Tanaka M, Toyoda A, Miwa K, Kasai K, Fuji K, et al. Polar localization and degradation of Arabidopsis boron transporters through distinct trafficking pathways. Proc Natl Acad Sci USA. 2010;107:5220–5. doi: 10.1073/pnas.0910744107. [PubMed] [Cross Ref]
20. Kasai K, Takano J, Miwa K, Toyoda A, Fujiwara T. High boron-induced ubiquitination regulates vacuolar sorting of the BOR1 borate transporter in Arabidopsis thaliana. J Biol Chem. 2011;286:6175–83. doi: 10.1074/jbc.M110.184929. [PubMed] [Cross Ref]
21. Sega GA. A review of the genetic effects of ethyl methanesulfonate. Mutat Res. 1984;134:113–42. doi: 10.1016/0165-1110(84)90007-1. [PubMed] [Cross Ref]
22. Burns PA, Allen FL, Glickman BW. DNA sequence analysis of mutagenicity and site specificity of ethyl methanesulfonate in Uvr+ and UvrB- strains of Escherichia coli. Genetics. 1986;113:811–9. [PubMed]
23. Kieber JJ, Rothenberg M, Roman G, Feldmann KA, Ecker JR. CTR1, a negative regulator of the ethylene response pathway in Arabidopsis, encodes a member of the raf family of protein kinases. Cell. 1993;72:427–41. doi: 10.1016/0092-8674(93)90119-B. [PubMed] [Cross Ref]
24. Cao XF, Linstead P, Berger F, Kieber J, Dolan L. Differential ethylene sensitivity of epidermal cells is involved in the establishment of cell pattern in the Arabidopsis root. Physiol Plant. 1999;106:311–7. doi: 10.1034/j.1399-3054.1999.106308.x. [PubMed] [Cross Ref]
25. Gibson SI, Laby RJ, Kim D. The sugar-insensitive1 (sis1) mutant of Arabidopsis is allelic to ctr1. Biochem Biophys Res Commun. 2001;280:196–203. doi: 10.1006/bbrc.2000.4062. [PubMed] [Cross Ref]
26. Kinoshita A, Betsuyaku S, Osakabe Y, Mizuno S, Nagawa S, Stahl Y, et al. RPK2 is an essential receptor-like kinase that transmits the CLV3 signal in Arabidopsis. Development. 2010;137:3911–20. doi: 10.1242/dev.048199. [PubMed] [Cross Ref]
27. Kiyohara S, Sawa S. CLE signaling systems during plant development and nematode infection. Plant Cell Physiol. 2012 doi: 10.1093/pcp/pcs136. In press. [PubMed] [Cross Ref]
28. Yamada M, Sawa S. The function of CLE and other plant peptide hormones during plant root development. Curr Opin Plant Biol. In press.
29. Fiers M, Golemiec E, Xu J, van der Geest L, Heidstra R, Stiekema W, et al. The 14-amino acid CLV3, CLE19, and CLE40 peptides trigger consumption of the root meristem in Arabidopsis through a CLAVATA2-dependent pathway. Plant Cell. 2005;17:2542–53. doi: 10.1105/tpc.105.034009. [PubMed] [Cross Ref]
30. International Rice Genome Sequencing Project The map-based sequence of the rice genome. Nature. 2005;436:793–800. doi: 10.1038/nature03895. [PubMed] [Cross Ref]
31. Koshiba T, Kobayashi M, Ishihara A, Matoh T. Boron nutrition of cultured tobacco BY-2 cells. VI. Calcium is involved in early responses to boron deprivation. Plant Cell Physiol. 2010;51:323–7. doi: 10.1093/pcp/pcp179. [PubMed] [Cross Ref]
32. Martín-Rejano EM, Camacho-Cristóbal JJ, Herrera-Rodríguez MB, Rexach J, Navarro-Gochicoa MT, González-Fontes A. Auxin and ethylene are involved in the responses of root system architecture to low boron supply in Arabidopsis seedlings. Physiol Plant. 2011;142:170–8. doi: 10.1111/j.1399-3054.2011.01459.x. [PubMed] [Cross Ref]
33. Kobayashi M, Matoh T, Azuma J. Two chains of rhamnogalacturonan II are cross-linked by borate-diol ester bonds in higher plant cell walls. Plant Physiol. 1996;110:1017–20. [PubMed]
34. O’Neill MA, Warrenfeltz D, Kates K, Pellerin P, Doco T, Darvill AG, et al. Rhamnogalacturonan-II, a pectic polysaccharide in the walls of growing plant cell, forms a dimer that is covalently cross-linked by a borate ester. In vitro conditions for the formation and hydrolysis of the dimer. J Biol Chem. 1996;271:22923–30. [PubMed]
35. Hermans C, Porco S, Vandenbussche F, Gille S, De Pessemier J, Van Der Straeten D, et al. Dissecting the role of CHITINASE-LIKE1 in nitrate-dependent changes in root architecture. Plant Physiol. 2011;157:1313–26. doi: 10.1104/pp.111.181461. [PubMed] [Cross Ref]
36. Bonin CP, Potter I, Vanzin GF, Reiter WD. The MUR1 gene of Arabidopsis thaliana encodes an isoform of GDP-D-mannose-4,6-dehydratase, catalyzing the first step in the de novo synthesis of GDP-L-fucose. Proc Natl Acad Sci USA. 1997;94:2085–90. doi: 10.1073/pnas.94.5.2085. [PubMed] [Cross Ref]
37. O’Neill MA, Eberhard S, Albersheim P, Darvill AG. Requirement of borate cross-linking of cell wall rhamnogalacturonan II for Arabidopsis growth. Science. 2001;294:846–9. doi: 10.1126/science.1062319. [PubMed] [Cross Ref]
38. Li Y, Smith C, Corke F, Zheng L, Merali Z, Ryden P, et al. Signaling from an altered cell wall to the nucleus mediates sugar-responsive growth and development in Arabidopsis thaliana. Plant Cell. 2007;19:2500–15. doi: 10.1105/tpc.106.049965. [PubMed] [Cross Ref]
39. Fujiwara T, Hirai MY, Chino M, Komeda Y, Naito S. Effects of sulfur nutrition on expression of the soybean seed storage protein genes in transgenic petunia. Plant Physiol. 1992;99:263–8. doi: 10.1104/pp.99.1.263. [PubMed] [Cross Ref]
40. Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, et al. Integrative genomics viewer. Nat Biotechnol. 2011;29:24–6. doi: 10.1038/nbt.1754. [PMC free article] [PubMed] [Cross Ref]

Articles from Plant Signaling & Behavior are provided here courtesy of Landes Bioscience