Tsukasa Nunome

Simple Sequence Repeat-Based Comparative Genomics Between Brassica rapa and Arabidopsis thaliana: The Genetic Origin of Clubroot Resistance

An SSR-based linkage map was constructed in Brassica rapa. It includes 113 SSR, 87 RFLP, and 62 RAPD markers. It consists of 10 linkage groups with a total distance of 1005.5 cM and an average distance of 3.7 cM. SSRs are distributed throughout the linkage groups at an average of 8.7 cM. Synteny between B. rapa and a model plant, Arabidopsis thaliana, was analyzed. A number of small genomic segments of A. thaliana were scattered throughout an entire B. rapa linkage map. This points out the complex genomic rearrangements during the course of evolution in Cruciferae. A 282.5-cM region in the B. rapa map was in synteny with A. thaliana. Of the three QTL (Crr1, Crr2, and Crr4) for clubroot resistance identified, synteny analysis revealed that two major QTL regions, Crr1 and Crr2, overlapped in a small region of Arabidopsis chromosome 4. This region belongs to one of the disease-resistance gene clusters (MRCs) in the A. thaliana genome. These results suggest that the resistance genes for clubroot originated from a member of the MRCs in a common ancestral genome and subsequently were distributed to the different regions they now inhabit in the process of evolution.

Isolation and characterization of microsatellites in Brassica rapa L.

Abstract.We report here the isolation and characterization of microsatellites, or simple sequence repeats (SSRs), in Brassica rapa. The size-fractionated genomic library was screened with (GA)15 and (GT)15 oligonucleotide probes. A total of 58 clones were identified as having the microsatellite repeats, and specific primer pairs were designed for 38 microsatellite loci. All primer pairs, except two, amplified fragments having the sizes expected from the sequences. Of the 36 primer pairs, 35 amplified polymorphic loci in 19 cultivars of B. rapa, while monomorphism was observed in only one primer pair. A total of 232 alleles was identified by the 36 primer pairs in 19 cultivars of B. rapa, and these primer pairs were examined also in nine Brassicaceae species. Most of the 36 primer pairs amplified the loci in the Brassicaceae species. Segregation of the microsatellites was studied in an F2 population from a cross of doubled-haploid lines DH27 × G309. The microsatellites segregated in a co-dominant manner. These results indicate that the microsatellites isolated in this study were highly informative and could be useful tools for genetic analysis in B. rapa and other related species.

Draft genome sequence of eggplant (Solanum melongena L.): the representative solanum species indigenous to the old world.

Unlike other important Solanaceae crops such as tomato, potato, chili pepper, and tobacco, all of which originated in South America and are cultivated worldwide, eggplant (Solanum melongena L.) is indigenous to the Old World and in this respect it is phylogenetically unique. To broaden our knowledge of the genomic nature of solanaceous plants further, we dissected the eggplant genome and built a draft genome dataset with 33,873 scaffolds termed SME_r2.5.1 that covers 833.1 Mb, ca. 74% of the eggplant genome. Approximately 90% of the gene space was estimated to be covered by SME_r2.5.1 and 85,446 genes were predicted in the genome. Clustering analysis of the predicted genes of eggplant along with the genes of three other solanaceous plants as well as Arabidopsis thaliana revealed that, of the 35,000 clusters generated, 4,018 were exclusively composed of eggplant genes that would perhaps confer eggplant-specific traits. Between eggplant and tomato, 16,573 pairs of genes were deduced to be orthologous, and 9,489 eggplant scaffolds could be mapped onto the tomato genome. Furthermore, 56 conserved synteny blocks were identified between the two species. The detailed comparative analysis of the eggplant and tomato genomes will facilitate our understanding of the genomic architecture of solanaceous plants, which will contribute to cultivation and further utilization of these crops.

Gene expression analysis in cadmium-stressed roots of a low cadmium-accumulating solanaceous plant, Solanum torvum

Solanum torvum Sw. cv. Torubamubiga (TB) is a low cadmium (Cd)-accumulating plant. To elucidate the molecular mechanisms of the Cd acclimation process in TB roots, transcriptional regulation was analysed in response to mild Cd treatment: 0.1 μM CdCl2 in hydroponic solution. A unigene set consisting of 6296 unigene sequences was constructed from 18 816 TB cDNAs. The distribution of functional categories was similar to tomato, while 330 unigenes were suggested to be TB specific. For expression profiling, the SuperSAGE method was adapted for use with Illumina sequencing technology. Expression tag libraries were constructed from Cd-treated (for 3 h, 1 d, and 3 d) and untreated roots, and 34 269 species of independent tags were collected. Moreover, 6237 tags were ascribed to the TB or eggplant (aubergine) unigene sequences. Time-course changes were examined, and 2049 up- and 2022 down-regulated tags were identified. Although no tags annotated to metal transporter genes were significantly regulated, a tag annotated to AtFRD3, a xylem-loading citrate transporter, was down-regulated. In addition to induction of heavy metal chaperone proteins, antioxidative and sulphur-assimilating enzymes were induced, confirming that oxidative stress developed even using a mild Cd concentration. Rapid repression of dehydration-related transcription factors and aquaporin isoforms suggests that dehydration stress is a potential constituent of Cd-induced biochemical impediments. These transcriptional changes were also confirmed by real-time reverse transcription-PCR. Further additions of TB unigene sequences and functional analysis of the regulated tags will reveal the molecular basis of the Cd acclimation process, including the low Cd-accumulating characteristics of TB.

Read2Marker: a data processing tool for microsatellite marker development from a large data set.

Microsatellites, or simple sequence repeats (SSRs), are widely considered to be one of the best sources of DNA polymorphism that can be utilized for developing genetic markers (1). Recently, several methods for making microsatellite-enriched genomic libraries were reported (2,3) and, as a consequence, the random sequencing of a microsatellite-enriched genomic library becomes a possible strategy for the mass development of microsatellite markers. Despite the progress in microsatellite isolation, it remains timeconsuming and laborious to handle the large amounts of sequence data. A set of script programs, read2Marker, was developed as a solution to this problem (Figure 1). In cooperation with phred/ phrap (4), BLAST (5), Primer3 (6), and T-Coffee (7), read2Marker can process an unlimited number of raw output data files from an automatic sequencer. It can run on a MacOS X or Linux platform and generates a report on many aspects of the nonredundant microsatellite-containing clones such as whole sequence, microsatellite core sequences, and PCR primers. A set of trace data obtained by sequencing from both ends of a clone is base-called by phred and assembled by phrap. Based on quality values (QVs), the high QV part of the assembled sequence is used for further processing. When data in plain text format without QVs are to be processed, all of the nucleotides are automatically assigned a QV of 20 to disable the quality threshold aspect of the read2Marker program. Subsequently, a microsatellite core is screened in the clone sequence by a newly developed tool, srchssr2. Initially, srchssr2 searches for repeated sequences of four dinucleotide motifs (GA, GT, AT, GC) and 10 trinucleotide motifs (AAC, AAG, ACC, ACG, ACT, AGC, AGG, ATC, TAT, CGC) independently in the clone sequence, and the core with the highest repeat number is selected. Next, a core structure neighboring the selected longest core is rescreened. If the gap between these cores is shorter than a maximum gap length set by the user in advance, then the two cores and gap sequence are connected to make a new core sequence. These procedures are repeated, and finally the composite microsatellite core in the sequence is identified. Finally, a PCR primer set for amplifying the microsatellite core is designed using the Primer3 program. The algorithm for finding microsatellites with srchssr2 is quite simple, but it is efficient enough to find cores. It can find composite repeats composed of different motifs and extracts 5′and 3′-flanking sequences that will be used for a redundancy check in the next step. Several SSR-finding programs with more sophisticated algorithms such as Msatfinder (M.I. Thurston and D. Field, www.bioinf.ceh.ac.uk/msatfinder), TROLL (8), Sputnik (C. Abajian, espressosoftware.com/pages/sputnik. jsp), and Tandem repeats finder (9) are available, and some of them could be substituted for srchssr2 with some modification, especially in cases where screening for motifs of more than three nucleotides is required. In order to remove redundant sequences, a homology check is performed using the BLAST algorithm. Two outer sequences of the clone flanking the microsatellite core are independently set as query sequences and BLASTed to a database comprising unique clones. If at least one of the two core-flanking sequences does not provide a hit from the database, the clone is judged as unique and added to

A protocol for the construction of microsatellite enriched genomic library

An improved protocol for constructing microsatellite-enriched libraries was developed. The procedure depends on digesting genomic DNA with a restriction enzyme that generates blunt-ends, and on ligating linkers that, when dimerized, create a restriction site for a different blunt-end producing restriction enzyme. Efficient ligation of linkers to the genomic DNA fragments is achieved by including restriction enzymes in the ligation reaction that eliminate unwanted ligation products. After ligation, the reaction mixture is subjected to subtractive hybridization without purification. DNA fragments containing microsatellites are captured by biotin-labeled oligonucleotide repeats and recovered using streptavidin-coated beads. The recovered fragments are amplified by PCR using the linker sequence as primer, and cloned directly into a plasmid vector. The linker has the sequence GTTT on the 5′ end, which promotes efficient adenylation of the 3′ ends of the PCR products. Consequently, the amplified fragments could be cloned into vectors without purification. This procedure enables efficient enrichment and cloning of microsatellite sequences, resulting in a library with a low level of redundancy.

Identification and Characterization of Crr1a, a Gene for Resistance to Clubroot Disease (Plasmodiophora brassicae Woronin) in Brassica rapa L.

Clubroot disease, caused by the obligate biotrophic protist Plasmodiophora brassicae Woronin, is one of the most economically important diseases of Brassica crops in the world. Although many clubroot resistance (CR) loci have been identified through genetic analysis and QTL mapping, the molecular mechanisms of defense responses against P. brassicae remain unknown. Fine mapping of the Crr1 locus, which was originally identified as a single locus, revealed that it comprises two gene loci, Crr1a and Crr1b. Here we report the map-based cloning and characterization of Crr1a, which confers resistance to clubroot in Brassica rapa. Crr1aG004, cloned from the resistant line G004, encodes a Toll-Interleukin-1 receptor/nucleotide-binding site/leucine-rich repeat (TIR-NB-LRR) protein expressed in the stele and cortex of hypocotyl and roots, where secondary infection of the pathogen occurs, but not in root hairs, where primary infection occurs. Gain-of-function analysis proved that Crr1aG004 alone conferred resistance to isolate Ano-01 in susceptible Arabidopsis and B. rapa. In comparison, the susceptible allele Crr1aA9709 encodes a truncated NB-LRR protein, which lacked more than half of the TIR domain on account of the insertion of a solo-long terminal repeat (LTR) in exon 1 and included several substitutions and insertion-deletions in the LRR domain. This study provides a basis for further molecular analysis of defense mechanisms against P. brassicae and will contribute to the breeding of resistant cultivars of Brassica vegetables by marker-assisted selection. Data deposition The sequence reported in this paper has been deposited in the GenBank database (accession no. AB605024).

Genome-wide survey of artificial mutations induced by ethyl methanesulfonate and gamma rays in tomato

Summary Genome‐wide mutations induced by ethyl methanesulfonate (EMS) and gamma irradiation in the tomato Micro‐Tom genome were identified by a whole‐genome shotgun sequencing analysis to estimate the spectrum and distribution of whole‐genome DNA mutations and the frequency of deleterious mutations. A total of ~370 Gb of paired‐end reads for four EMS‐induced mutants and three gamma‐ray‐irradiated lines as well as a wild‐type line were obtained by next‐generation sequencing technology. Using bioinformatics analyses, we identified 5920 induced single nucleotide variations and insertion/deletion (indel) mutations. The predominant mutations in the EMS mutants were C/G to T/A transitions, while in the gamma‐ray mutants, C/G to T/A transitions, A/T to T/A transversions, A/T to G/C transitions and deletion mutations were equally common. Biases in the base composition flanking mutations differed between the mutagenesis types. Regarding the effects of the mutations on gene function, >90% of the mutations were located in intergenic regions, and only 0.2% were deleterious. In addition, we detected 1 140 687 spontaneous single nucleotide polymorphisms and indel polymorphisms in wild‐type Micro‐Tom lines. We also found copy number variation, deletions and insertions of chromosomal segments in both the mutant and wild‐type lines. The results provide helpful information not only for mutation research, but also for mutant screening methodology with reverse‐genetic approaches.

Genetic mapping of AFLP markers in Japanese bunching onion (Allium fistulosum)

SummaryThe first genetic linkage map of Japanese bunching onion (Allium fistulosum) based primarily on AFLP markers was constructed using reciprocally backcrossed progenies. They were 120 plants each of (P1)BC1 and (P2)BC1 populations derived from a cross between single plants of two inbred lines: D1s-15s-22 (P1) and J1s-14s-20 (P2). Based on the (P2)BC1 population, a linkage map of P1 was constructed. It comprises 164 markers – 149 amplified fragment length polymorphisms (AFLPs), 2 cleaved amplified polymorphic sequences (CAPSs), and 12 simple sequence repeats (SSRs) from Japanese bunching onion, and 1 SSR from bulb onion (A. cepa) – on 15 linkage groups covering 947 centiMorgans (cM). The linkage map of P2 was constructed with the (P1)BC1 population and composed of 120 loci – 105 AFLPs, 1 CAPS, and 13 SSRs developed from Japanese bunching onion and 1 SSR from bulb onion – on 14 linkage groups covering 775 cM. Both maps were not saturated but were considered to cover the majority of the genome. Nine linkage groups in P2 map were connected with their counterparts in P1 map using co-dominant anchor markers, 13 SSRs and 1 CAPS.

A simulation-based breeding design that uses whole-genome prediction in tomato

Efficient plant breeding methods must be developed in order to increase yields and feed a growing world population, as well as to meet the demands of consumers with diverse preferences who require high-quality foods. We propose a strategy that integrates breeding simulations and phenotype prediction models using genomic information. The validity of this strategy was evaluated by the simultaneous genetic improvement of the yield and flavour of the tomato (Solanum lycopersicum), as an example. Reliable phenotype prediction models for the simulation were constructed from actual genotype and phenotype data. Our simulation predicted that selection for both yield and flavour would eventually result in morphological changes that would increase the total plant biomass and decrease the light extinction coefficient, an essential requirement for these improvements. This simulation-based genome-assisted approach to breeding will help to optimise plant breeding, not only in the tomato but also in other important agricultural crops.