Skip Navigation


DNA Research Advance Access originally published online on November 2, 2009
DNA Research 2009 16(6):311-323; doi:10.1093/dnares/dsp020
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary Data
Right arrowOA All Versions of this Article:
16/6/311    most recent
dsp020v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Request Permissions
Google Scholar
Right arrow Articles by Li, F.
Right arrow Articles by Nishio, T.
PubMed
Right arrow PubMed Citation
Right arrow Articles by Li, F.
Right arrow Articles by Nishio, T.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2009. Published by Oxford University Press on behalf of Kazusa DNA Research Institute.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/2.5/uk/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

A Brassica rapa Linkage Map of EST-based SNP Markers for Identification of Candidate Genes Controlling Flowering Time and Leaf Morphological Traits

Feng Li, Hiroyasu Kitashiba, Kiyofumi Inaba {dagger} and Takeshi Nishio*

Laboratory of Plant Breeding and Genetics, Graduate School of Agricultural Science, Tohoku University, 1-1 Tsutsumidori-Amamiyamachi, Aoba, Sendai, Miyagi 981-8555, Japan

Received 30 June 2009; accepted 20 September 2009.


   Abstract
Top
 Abstract
1. Introduction
 2. Materials and methods
 3. Results and discussion
 Supplementary data
 Funding
 Acknowledgements
 References
 
For identification of genes responsible for varietal differences in flowering time and leaf morphological traits, we constructed a linkage map of Brassica rapa DNA markers including 170 EST-based markers, 12 SSR markers, and 59 BAC sequence-based markers, of which 151 are single nucleotide polymorphism (SNP) markers. By BLASTN, 223 markers were shown to have homologous regions in Arabidopsis thaliana, and these homologous loci covered nearly the whole genome of A. thaliana. Synteny analysis between B. rapa and A. thaliana revealed 33 large syntenic regions. Three quantitative trait loci (QTLs) for flowering time were detected. BrFLC1 and BrFLC2 were linked to the QTLs for bolting time, budding time, and flowering time. Three SNPs in the promoter, which may be the cause of low expression of BrFLC2 in the early-flowering parental line, were identified. For leaf lobe depth and leaf hairiness, one major QTL corresponding to a syntenic region containing GIBBERELLIN 20 OXIDASE 3 and one major QTL containing BrGL1, respectively, were detected. Analysis of nucleotide sequences and expression of these genes suggested possible involvement of these genes in leaf morphological traits.

Key words: DNA markers; synteny; bolting time; leaf lobe; leaf hairiness


   1. Introduction
Top
 Abstract
 1. Introduction
2. Materials and methods
 3. Results and discussion
 Supplementary data
 Funding
 Acknowledgements
 References
 
The Brassica genus comprises a number of important vegetables and condiment crops and is also a source of oil seed and fodder. Among the Brassica species, B. rapa (AA, 2n = 20), B. nigra (BB, 2n = 16), and B. oleracea (CC, 2n = 18) are diploid species, and B. juncea (AABB, 2n = 36), B. napus (AACC, 2n = 38), and B. carinata (BBCC, 2n = 34) are amphidiploids having combinations of the genomes of these diploids. Genetic maps with molecular markers are necessary to understand the origin and relationship among the genomes of the Brassica species and can be utilized in applied genetics and breeding of Brassica crops. In past decades, genetic maps of Brassica have been constructed with a range of marker types, including RAPD (Random Amplified Polymorphic DNA), AFLP (Amplified Fragment Length Polymorphism), RFLP (Restriction Fragment Length Polymorphism), and SSR (Simple Sequence Repeats) markers, and these maps have contributed to genetic analysis of quantitative traits. RFLP and SSR markers anchored in genetic maps have also been used for comparison of the chromosome organization among Brassica species, even between Brassica and its related model plant Arabidopsis thaliana.1Go–3Go The numbers of currently available mapped RFLP and SSR markers in Brassica are both more than 1000 in total.3Go–6Go

Brassica rapa includes a variety of vegetables such as Chinese cabbage, Pakchoi, and turnip as well as oilseed crops such as turnip rape and sarson. More than 20 genetic linkage maps have been generated in B. rapa, some of which have been applied to quantitative trait locus (QTL) analysis of morphological traits,7Go,8Go and agronomical traits including fatty acid content,9Go glucosinolate accumulation,10Go phytate and phosphate concentrations,11Go and disease resistance.12Go,13Go Because Brassica belongs to the same family as A. thaliana, in which the complete genome sequence has been determined and functions of many genes have been well characterized, a clear synteny map between them can offer a powerful tool to find candidate genes after QTL analysis. By synteny analysis of QTL regions of B. rapa with the A. thaliana genome, Saito et al.12Go have fine mapped the clubroot resistance gene, Crr3, and Zhang et al.14Go have successfully cloned a gene controlling hairiness and seed coat color traits. However, AFLP and RAPD markers are anonymous markers, which make it difficult to compare and combine linkage maps. Some RFLP probes give rise to hybridization with multiple loci because of genomic segment replication,1Go,15Go and SSR markers are usually located in non-coding sequences, which are less conserved between the replication blocks than coding regions.5Go These shortcomings may limit their application in genome analysis. Other marker types such as SCAR (Sequence Characterized Amplified Region) and CAPS (Cleaved Amplified Polymorphic Sequence) have also been produced and used to construct genetic maps,5Go,14Go but the limited number of markers with polymorphism between different lines in a species does not meet the need of fine mapping. At present, in most detailed linkage map of B. rapa, which comprises 556 markers, only 120 markers have been detected in homologous regions of A. thaliana.5Go

Single nucleotide polymorphisms (SNPs) are the most common types of DNA polymorphism in genomes. In plants, SNP frequencies vary widely: 1 SNP/124 bp in coding region and 1 SNP/31 bp in the non-coding region of 36 inbred lines of maize;16Go 1 SNP/72 bp in expressed genes and 1 SNP/58 bp in non-coding sequences among a panel of 13 lines of sugar beet (Beta vulgaris L.);17Go and 1 SNP/2.1 kb to 1 SNP/1.2 kb between two cultivars of B. napus, Tapidor and Ningyou 7.18Go For genotyping of SNPs, gel electrophoresis, fluorometry, DNA microarrays, MALDI-TOF mass spectrometry, and labeled oligonucleotide hybridization are the most commonly used techniques.19Go,20Go Among them, for a small-scale operation in common laboratories, the dot-blot-SNP technique20Go,21Go based on labeled oligonucleotide hybridization is considered to be a low-cost, labor-saving technique for analysis of a large number of individuals. SNPs have been extensively applied to genetic studies, including association analysis of candidate genes with phenotypic variation,22Go fine mapping of QTLs,14Go EST-mapping,17Go linkage disequilibrium-based association mapping,23Go and genetic diversity assessment.24Go Whereas, in Brassica, very limited SNP markers, ~50, have been developed,14Go,18Go and no SNP linkage map has been published.

As ESTs are highly conserved in Brassica and A. thaliana, the linkage map based on ESTs can eventually provide a direct genomic comparison of macro- and micro-collinearity across these species. In the present study, we mainly used EST-based SNP markers to construct a B. rapa linkage map, which was then compared with A. thaliana linkage map based on homologous loci and revealed a fine synteny relationship. Flowering time and leaf morphological traits are important characteristics as vegetables, and therefore ones of breeding objectives. It has been indicated that multiple loci are involved in the variation of these traits,7Go,8Go,25Go,26Go but limited gene information has been obtained. We analyzed QTLs for these traits and inferred some candidate genes from the corresponding syntenic regions of A. thaliana.


   2. Materials and methods
Top
 Abstract
 1. Introduction
 2. Materials and methods
3. Results and discussion
 Supplementary data
 Funding
 Acknowledgements
 References
 
2.1. Plant materials and growth conditions
Brassica rapa cv. ‘Yellow Sarson’ C634, an early-flowering Indian oilseed rape, has been maintained by self-pollination for more than five generations. A doubled haploid line P11 of B. rapa cv. ‘Osome’, a Japanese commercial variety of late flowering leafy vegetable, was provided by Dr Kuginuki of Asahi Noen Seed Co. Ltd. The two lines were crossed to produce an F2 population. As temperature is an important factor for flowering, the F2 population was divided into two subpopulations of 132 and 134 individuals, which were cultivated using 24-cm-diameter pots in an unheated greenhouse and a heated greenhouse, respectively, in Sendai, Japan (38°16'N, 140°52'E), from October 2007 to April 2008. The average highest and lowest temperature of each month in Sendai is listed in Supplementary Table S1. The room temperature of the heated greenhouse was maintained above 10°C. Genomic DNA was prepared from leaves by a modified CTAB (cetyltrimethylammonium bromide) method.27Go

2.2. Investigation of phenotypes
The leaf characteristics were scored 60 days after sowing using the largest leaf of each plant. Lamina width (LW) and sinus width (SW) were measured as illustrated in Fig. 1, and lobe depth ratio (LDR) was calculated as (LW – SW)/LW. Leaf hairiness (LH) was scored with 0–5 scale: 0, hairless; 1, some hairs along the leaf edge; 2, some hairs on the margin of abaxial surface of leaf; 3, dense hairs on the margin of abaxial surface of leaf; 4, dense hairs covering ca. 1/3 of abaxial surface of leaf and some hairs on the adaxial veins; 5, dense hairs covering more than 1/3 of abaxial surface of leaf and dense hairs on the adaxial veins. Days to bolting (DBO), budding (DBU), and flowering (DFL) were recorded as the number of days from sowing to achievement of a 3-cm-high main flower stalk, emergence of buds, and opening of the first flower, respectively. Days from bolting to bud formation (DBF) was the difference between DBU and DBO.


Figure 1
View larger version (12K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		Figure 1. Diagram of traits of fully expanded leaves. LW and SW indicate lamina width and sinus width, respectively.

 
2.3. Primer design
In total, 533 sequences of genes or hypothetical genes in Brassica from BrassicaDB (http://brassica.bbsrc.ac.uk/BrassicaDB), 80 BAC sequences from the B. rapa Genome Project (http://www.brassica-rapa.org/BRGP/chromosomeSequence.jsp), 6 sequences of A. thaliana genes related to leaf shape from NCBI (http://www.ncbi.nlm.nih.gov/), and 308 EST sequences from the Raphanus sativus database (http://radish.plantbiology.msu.edu/index.php/Sequences:All) were selected to design PCR primers using the Primer3 software (http://frodo.wi.mit.edu/), expected sizes of amplified fragments being between 600 and 1500 bp.

2.4. SCAR and CAPS analyses
PCR was performed in a reaction mixture of 20 µl, consisting of 40 ng genomic DNA as a PCR template, 0.5 µM of primers, 1 U of Taq DNA polymerase (TaKaRa Ex Taq® or TaKaRa TaqTM, Takara Biomedicals, Japan), 1xEx Taq buffer or 1xPCR buffer, and 200 µM of dNTPs. The thermal cycle of PCR was set to be as follows: 1 min denaturation at 94°C, 40 cycles of 30 s denaturation at 94°C, 30 s annealing at 58°C, and 1 min 30 s extension at 72°C, and 1 min 30 s final extension at 72°C. PCR products amplified by the primers of the Brassica sequence were separated on 0.8% agarose gel for SCAR analysis. For CAPS analysis, PCR products were digested by a restriction enzyme MboI or MspI and then were separated by 6% polyacrylamide gel in 1xTBE buffer. The resulting DNA bands of SCAR and CAPS were stained with ethidium bromide.

2.5. SSR analysis
From Brassica microsatellite information exchange (http://www.brassica.info/resource/markers/ssr-exchange.php), 80 SSRs were selected to analyze polymorphism between the two parental genotypes. The PCR products were separated by 8% polyacrylamide gel in 1xTBE buffer.

2.6. Nucleotide sequencing and dot-blot-SNP analysis
When SCAR and CAPS analyses did not differentiate the two parental genotypes, the amplification products were sequenced with a DNA analyzer (CEQ2000, Beckman Coulter, http://www.beckmancoulter.com/Default.asp?bhfv=6) and aligned using SEQUENCHER version 4.7 to identify SNPs. PCR products amplified by the primers of the R. sativus EST sequence were directly sequenced to identify SNPs without SCAR and CAPS analyses, and the SNP frequency between P11 and C634 was estimated.

5'-Biotin-labeled probes and bridge probes21Go were designed for dot-blot-SNP analysis (Supplementary Table S2). DNA fragments containing SNPs were amplified by PCR using the primers listed in Supplementary Table S2. The PCR products were mixed with an equal volume of a denaturation solution containing 0.4 N NaOH and 10 mM EDTA and dot-blotted onto a nylon membrane by Multi-pin Blotter (ATTO, Japan). After UV crosslinking using GS Gene Linker UV Chamber (Bio-Rad Laboratories, USA), the membrane was hybridized for 3 h or overnight with the probes and washed with washing buffer (Supplementary Table S2). Digoxigenin and biotin were reacted with anti-DIG-IgG alkaline phosphatase conjugate (Roche Diagnostics, Switzerland) and streptavidin-alkaline phosphatase conjugate (Promega Corp., Madison, WI, USA), respectively, and alkaline phosphatase activity was detected using CSPD (Roche Diagnostics).

2.7. Linkage analysis and map construction
All 241 DNA markers were scored in the unheated-greenhouse subpopulation. To investigate the environmental effect on QTLs, about half of the markers covering the whole genome of B. rapa were scored in the heated-greenhouse subpopulation. Linkage analysis and map construction were performed using Antmap version 1.2.28Go Linkage groups were identified in the LOD (logarithm of odds) threshold range of 3.0, and the Kosambi mapping function was used to convert recombination frequencies into map distances (cM).

QTL analysis was performed using a composite interval-mapping analysis with Windows QTL Cartographer v2.5.29Go A permutation test was applied to each data set (1000 repetitions) to determine the LOD thresholds (P = 0.05). LOD values of 3.6 for the unheated-greenhouse subpopulation and 3.8 for the heated-greenhouse subpopulation were used as significant thresholds for the presence of a candidate QTL.

2.8. Locus nomenclature
The markers based on the genes or hypothetical genes were named <Br><Name 2-6 letter code>, and those based on ESTs of R. sativus were named <Br><EST name>. The marker assay types were designated with single letters, following the recommendations of De Vicente et al.30Go as follows: m, SSR; p, CAPS; c, SCAR; s, SNP.

2.9. Comparison with the A. thaliana genome sequence
The sequences of sequence-tagged SNP, SCAR, and CAPS loci in the map were aligned with A. thaliana genome sequences using the BLASTN program of the National Centre for Biotechnology Information (http://www.ncbi.nlm.nih.gov/). We regarded the sequences as homologous loci of A. thaliana genome with a threshold value of E < 10–20, and the homologous loci of A. thaliana were located on a physical map of A. thaliana, according to the TAIR database (http://www.arabidopsis.org/). The regions having conserved collinearity with A. thaliana were regarded as homologous syntenic regions. In the A. thaliana syntenic regions corresponding to QTLs in the present study, loci defined as transcribed units (e.g. AT2G03340) in the A. thaliana genome were searched in the TAIR database (http://www.arabidopsis.org/) to find the related candidate genes.

2.10. Sequence and expression analysis of candidate genes
Young leaves were collected from 30-day-old plants in the greenhouse at 20°C for RNA extraction. Total RNA was extracted from 30 mg of leaf using the SV Total RNA Isolation System (Promega Corp.). First-strand cDNA was synthesized from 1 µg of total RNA by reverse-transcription using a GE Healthcare first-strand cDNA Synthesis kit (GE Healthcare UK Limited, Little Chalfont, Buckinghamshire, UK). RT–PCR was performed by specific primers of BrFLCs, listed in Supplementary Table S3, by 30 cycles of 94°C for 30 s, 58°C for 30 s, and 72°C for 1 min. Real-time PCR was performed to semi-quantify mRNAs of BrTTG1 and BrGA20OX3 using a LightCycler (Roche Diagnostics) with SYBR Premix Ex Taq (Takara Biomedicals) according to Fujimoto et al.31Go Specific primers of BrTTG1 and BrGA20OX3 are listed in Supplementary Table S3.

Nucleotide sequences of three BAC clones, KBrH080A08 (accession number in NCBI: AC155344), KBrH004D11 (accession number in NCBI: AC155341), and KBrB004B12 (accession number in NCBI: AC189195 [GenBank] ), which carry BrFLC1, BrFLC2, and BrGA20OX3, respectively, and that of the gene BnTTG1 (accession number in NCBI: EF175930 [GenBank] ) were used to design primer sequences for nucleotide sequencing.


   3. Results and discussion
Top
 Abstract
 1. Introduction
 2. Materials and methods
 3. Results and discussion
Supplementary data
 Funding
 Acknowledgements
 References
 
3.1. Production of DNA markers and mapping
Of the 533 primer pairs designed from Brassica EST and BAC sequences and A. thaliana EST sequences, single DNA fragments were amplified by 383, in which 8 showed polymorphism between the parents of F2 population, i.e. C634 and P11, by SCAR analysis and 70 showed such polymorphism by CAPS analysis (Supplementary Table S4). Of the remaining 305 primers, 160 were used for PCR amplification and the PCR products of C634 and P11 were sequenced and aligned. As a result, 95 SNP markers were produced (Supplementary Table S5). Among the 288 primer pairs from ESTs of R. sativus, 132 amplified single DNA fragments. Of these, 100 PCR amplicons covering ~34 kb were sequenced, and 71 showed nucleotide variation with a total of 427 SNPs (12.6 SNPs/kb) between the parents, 56 of which were used as SNP markers. Surveying 80 SSR primers, 12 revealed polymorphism between the parents. In total, 241 markers, including 170 EST-based markers (120 being SNP markers), 59 BAC sequence-based markers (31 being SNP markers), and 12 SSR markers, were generated (Supplementary Table S6).

A total of 241 markers were assigned to 10 linkage groups (Fig. 2), designated as A01–A10, corresponding to the linkage groups of the JWF3p reference map (http://www.brassica-rapa.org/BRGP/chromosomeSequence.jsp), based on the SSR markers and the BAC sequence-based markers. The linkage group names A01–A10 correspond to R01–R10, a previously determined classification.2Go,5Go,15Go,32Go The linkage map had a total length of 1396 cM, and the average distance between ordered adjacent markers was 5.79 cM. The largest linkage group consisted of 38 markers and had a length of 176.7 cM, and the smallest consisted of 12 markers with a length of 73.2 cM. From the physical length of B. rapa, i.e. 717 Mb,33Go the relationship between the physical and genetic distances was 1 cM = ca. 514 kb.


Figure 2
View larger version (55K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		Figure 2. Linkage map of B. rapa with comparative maps of A. thaliana. Ten linkage groups are represented by vertical bars. Linkage group numbers correspond to the Brassica A genome linkage. Map distances are given in centimorgans (cM) on the left side of the linkage groups, and the locus designations are given on the right side. Loci tested for homology with A. thaliana using BLAST are shown to the right of the B. rapa linkage groups as colored vertical bars, which represent different chromosomes of A. thaliana. The locus names of the homologous regions of A. thaliana are given on the right side of the colored vertical bars and the megabase distances are given to the left. Groups of two or more markers showing homology with A. thaliana in collinearity are regarded as synteny regions. Except for both ends of each synteny region, the names and megabase distances of identified homologous loci are not shown, but are shown in the Supplementary Table S5 and Fig. S1. The A. thaliana chromosomes are colored as the bottom right of the figure.

 
Nucleotide sequences of ESTs for primer design in the present study were mainly derived from B. rapa, B. oleracea, and B. napus. Most of these ESTs have known or hypothetical functions (Supplementary Table S4). The nucleotides sequences of ESTs of R. sativus were also used to design primers and the sequencing results showed that PCR products of B. rapa and R. sativus had high nucleotide similarity, suggesting the feasibility of genome comparison between the two species using markers generated from the PCR products amplified by the same primers. The linkage map also contains 59 BAC sequence-based markers, which are located across all 10 linkage groups, and the order of these markers in each linkage group was almost the same as that of the JWF3p reference map (http://www.brassica-rapa.org/BRGP/chromosomeSequence.jsp).

A conspicuous characteristic of the present map is that more than half of the markers, 151 markers, are SNP markers, which greatly contribute to increasing the number of sequence-based markers. The SNP frequency between the parents was estimated to be 12.6 SNPs/kb. Such high frequencies of SNPs have also been reported in other species such as maize,16Go sugar beet,17Go and forage grass species Lolium perenne.34Go Therefore, almost all of the EST markers mapped in the present study can be applied to segregating populations of other parental lines in B. rapa, even in different species of Brassica or Raphanus, if they are used as SNP markers.

3.2. Synteny between B. rapa and A. thaliana
All the markers except the SSR markers in the present map were sequence-tagged markers. BLASTN was used to test homology between the B. rapa sequence and the A. thaliana genome sequence, the results being shown in Supplementary Table S7. By a significance threshold of E < 10–20, 223 markers had homologous regions in A. thaliana, and these homologous loci covered nearly the whole genome of A. thaliana, except that only one homolog was located on the short arm of chromosome 2 (Supplementary Fig. S1). The five chromosomes of A. thaliana were divided into many segments and distributed to various regions of the B. rapa map (Fig. 2). There were 33 large syntenic regions containing at least three markers and having conserved collinearity with the homologues in A. thaliana. By comparing the syntenic regions from the same chromosome of A. thaliana, doubled and triplicated regions were identified. The region from 3.5 to 6.0 Mb in chromosome 1 was triplicated in A06, A08, and A09 of B. rapa, and that from 3.0 to 5.5 Mb in chromosome 5 was triplicated in A02, A03, and A10. The region between 12.9 and 14.9 Mb in chromosome 2 was doubled in A04 and A05, that between 2.8 and 7.5 Mb in chromosome 3 was doubled in A03 and A05, and that between 8.5 and 13.2 Mb in chromosome 4 was doubled in A01 and A08. Furthermore, these doubled regions also appeared to be partly conserved in other chromosomes of B. rapa. Segment duplications within the same chromosome were shown in A01, A02, A03, A05, A07, and A09. Insertion and rearrangement were the most apparent events and occurred in every chromosome of B. rapa.

Brassica and A. thaliana are considered to have originated from a common ancestor and diverged 14.5–20.4 million years ago.35Go Arabidopsis thaliana has a small genome (~146 Mb) with relatively little repetitive DNA and a high gene density,36Go whereas the diploid Brassica genome is significantly larger and has a lower gene density (B. rapa, 717 Mb).33Go,37Go Comparative genetic mapping between the diploid Brassica species and A. thaliana using RFLP probes has identified genome duplication and triplication in the Brassica genome and collinearity disrupted by multiple rearrangements.3Go,38Go,39Go Suwabe et al.2Go and Choi et al.5Go using 74 and 120 sequence-tagged markers of B. rapa, which have corresponding homologous regions in A. thaliana, have also identified some syntenic regions and illustrated the complex nature of the chromosomal rearrangements that may have occurred before/after differentiation of A. thaliana and B. rapa. Physical mapping and microsynteny analysis by nucleotide sequencing have further corroborated these findings.37Go,40Go The present comparative map of B. rapa and A. thaliana using 223 sequence-based markers revealed 33 large syntenic regions, which comprise the previously reported syntenic regions, and shows consistency.2Go,5Go Calculating physical distances of the syntenic regions in A. thaliana provided genome-wide synteny and micro-synteny, enabling detailed identification of duplicated and triplicated regions in every chromosome of A. thaliana. This increased detail will contribute to further study of genome structures and identification of candidate genes in B. rapa.

3.3. Analysis of QTLs and candidate genes for flowering time
There were significant differences in flowering traits, i.e. DBO, DBU, and DFL, between P11 and C634. Budding time of P11 was almost equal to bolting time; however, bolting time of C634 was significantly earlier than budding time. Flowering traits of F1 were similar to those of P11, indicating that late flowering trait was dominant over early-flowering trait. The traits including DBO, DBU, DFL, and DBF showed continuous phenotypic distribution in the F2 population under both heated and unheated conditions (Supplementary Fig. S2).

Under the two conditions, three QTLs of bolting time (DBOQTL1–3) were identified on A02, A07, and A10 (Table 1, Fig. 3). DBOQTL1 and DBOQTL2 were detected in both heated and unheated conditions, and the explained phenotypic variance was from 14.0% to 27.1%. DBOQTL3 was detected only in the heated condition with an explained phenotypic variance of 19.9%. Two QTLs (DBUQTL1/DFLQTL1 and DBUQTL2/DFLQTL2) for budding time and flowering time were detected on the same regions in A02 and A10 as DBOQTL1 and DBOQTL3 in two different conditions (Table 1, Fig. 3). One QTL for DBF (DBFQTL1) was detected in both conditions (Table 1, Fig. 3), and this QTL was co-located with DBOQTL2 on A07, but with an opposite effect, i.e. a P11 allele of this locus delayed the bolting time and accelerated bud formation so that the bolting time and the budding time were almost equal. Bolting time and budding time of B. rapa have also been analyzed previously, and some QTLs have been identified.25Go,26Go However, these linkage groups have not been assigned to the reference linkage groups and therefore it is not possible to compare these QTLs with the three QTLs detected in the present study.


Figure 3
View larger version (33K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		Figure 3. The genetic locations of QTL for days to bolting (DBO), budding (DBU), flowering (DFL), days from bolting to bud formation (DBF), lobe depth ratio (LDR), and leaf hairiness (LH). Syntenic regions of A. thaliana are shown to the left of the B. rapa linkage groups as patterned vertical bars, which represent different chromosomes of A. thaliana. Candidate genes and locus names of both ends of each syntenic region corresponding to each QTL are shown on the left side of the patterned vertical bars. The lengths of the arrows indicate the 2-LOD support intervals. The traits are indicated above each column. The direction of the arrowheads indicates the allelic effect: upward, P11 increases and C634 decreases for the F2 population; downward: P11 decreases and C634 increases for the F2 population. The numbers in the brackets under the arrows indicate the order numbers of the QTLs for each trait shown in Table 1.

 


View this table:
[in this window]
[in a new window]

  		Table 1. QTL analyses of flowering and leaf morphological traits in B. rapa under two treatments

 
BrFLC2 and BrFLC141Go were located in the regions of DBOQTL1/DBUQTL1/DFLQTL1 and DBOQTL3/DBUQTL2/DFLQTL2, respectively. These two loci have also been detected in previous study on QTL of flowering time.8Go,41Go,42Go In A. thaliana, FLC encodes a MADS box protein that acts as a dosage-dependent repressor of flowering.43Go,44Go Evidence has been presented that the FLC homologues in Brassica species act similar to AtFLC and play a central role in the repression of flowering.8Go,41Go,42Go,45Go In the present study, QTL with BrFLC2 was detected in both heated and unheated conditions, and that of BrFLC1 was sensitive to environment, which is in accord with previous findings.8Go,42Go,46Go

Expressions of BrFLC1, BrFLC3, and BrFLC5 were all detected in the leaves of P11 and C634 (Fig. 4) when the RT–PCR was performed for 30 cycles with the primer pairs of exons 1–4. In contrast, BrFLC2 expression in P11 was strong, whereas very faint in C634. Sequence analysis of BrFLC1 in C634 and P11 revealed three nucleotide polymorphisms (Fig. 5A). One was a nonsynonymous SNP in the first exon, one was in the second intron, and the other was at the 5' splicing site of the sixth intron. RT–PCR using a primer pair of exons 4–7 revealed two different transcripts of BrFLC1 in C634 (Supplementary Fig. S4), suggesting abnormal splicing of the C634 allele. This splicing site polymorphism has recently been reported to contribute greatly to flowering-time variation in B. rapa.47Go For BrFLC2, eight SNP sites were identified between C634 and P11 (Fig. 5B). Four of these SNPs were in the first and sixth introns, one was a synonymous SNP in the third exon, and the remaining three occurred in the region –800 bp upstream of the translation start site of BrFLC2, which might cause the difference of BrFLC2 expression between C634 and P11.


Figure 4
View larger version (32K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		Figure 4. RT–PCR analysis of BrFLCs transcripts in B. rapa lines, P11 and C634, using primer pairs of exons 1–4. Leaf samples were collected from plants grown for 30 days in the green house at 20°C. Actin was used as a control to demonstrate equal RNA loading.

 


Figure 5
View larger version (14K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		Figure 5. Nucleotide polymorphisms of BrFLC1 (A) and BrFLC2 (B) between P11 and C634. Exons (E), introns, and 5' sequences upstream of translational start site (+1) are represented by black, white, and gray boxes, respectively. Numbers above the boxes indicate SNP sites connected by dotted arrows showing nucleotides of P11 and C634 at the left and right of the slashes, respectively.

 
One QTL (DBFQTL1) for DBF with minus effect co-located with DBOQTL2. In DBFQTL1/DBOQTL2, a Brassica ortholog of CLAVATA1 (BrCLV1)48Go was located. In A. thaliana, the CLV1 gene is involved in controlling shoot apical and floral meristem size and contributes to the establishment and maintenance of floral meristem identity.49Go,50Go Opposite effects between DBFQTL1 and DBOQTL2 indicated that at least two related genes located and were responsible for bolting and bud formation, respectively. As the region of DBFQTL1/DBOQTL2 showed a fine syntenic relationship with A. thaliana, we searched the genes in the syntenic region from 27.7 to 30.2 Mb of chromosome 1 in A. thaliana (http://www.arabidopsis.org/) and revealed 705 genes, among which, besides CLV1, a gene involved in the control of flowering time FLOWERING LOCUS M (FLM) (AT1G77080)51Go was found. FLM, also called MAF1 (MADS AFFECTING FLOWERING1), encodes a floral repressor closely related to FLC.51Go Repression of FLM or non-functional FLM contribute to acceleration of flowering.52Go A BAC clone KBrH80C09, containing a tandem array of three MAF genes53Go and mapped on the bottom of A02, has been suggested to correspond to the QTL ‘VFR1’, which has been identified by Osborn et al.42Go The speculated homolog of MAF1 in A07 of B. rapa needs to be cloned and characterized.

3.4. Analysis of QTLs and candidate genes for leaf morphological traits
3.4.1. Lobe depth ratio
Two QTLs for LDR (LDRQTL1 and 2) were detected. A large percentage of phenotypic variance of 35.8% and 43.7% in the unheated-greenhouse subpopulation and the heated-greenhouse subpopulation, respectively, were explained by LDRQTL1 on A10. A minor QTL LDRQTL2, linked with a Brassica ortholog of ASYMETRIC LEAVES 1 (BrAS1) on A03, was only detected in the unheated-greenhouse subpopulation. In A. thaliana, expression of the KNOTTED-like homeobox (KNOX) genes in the shoot apical meristem is required for maintenance of a functional meristem, but suppressed in cells destined for leaf primordia.54Go,55Go AS1 maintains the down-regulation of KNOX expression in leaf primordia, and the as1 mutation leads to lobed leaves.56Go In Brassica, QTLs in A02 and A038Go and an AFLP maker on A0857Go have been indicated to be associated with leaf edge shape (LES). In the present study, we analyzed the trait of LDR calculated by (LW – SW)/LW. A similar trait, LES, scored as 1–4, has been analyzed by Lou et al.,8Go and a major QTL has been detected on the bottom of A03, where is seemed to be near to the LDRQTL2 in the present study.

We searched the genes in the syntenic region of LDRQTL1 from 2.1 to 3.2 Mb of chromosome 5 in A. thaliana (http://www.arabidopsis.org/) and found 348 genes, among which AtGA20OX3 (AT5G07200) encoding a gibberellin 20-oxidase was noted. It has been reported that reduction in GA signaling or biosynthesis increases the number of lobes per leaf of as1 mutation, and conversely, exogenous GA application or constitutive GA signaling suppresses ectopic expression of KNOX and both the number and depth of lobes.58Go Blast analysis showed a B. rapa homologue of AtGA20OX3 (BrGA20OX3) in a B. rapa BAC clone KBrB004B12, which contains the sequence of the BrA9s marker located in LDRQTL1. Identities of nucleotide sequences and deduced amino acid sequences between BrGA20OX and AtGA20OX3 were 92% and 86%, respectively (Supplementary Fig. S3). Semi-quantitative RT–PCR analysis showed that BrGA20OX3 expression in leaves of P11 was more than six times as high as that in leaves of C634 (Fig. 6A). The lobed leaves of C634 may be caused by the low expression of BrGA20OX3. Nucleotide sequence analysis of BrGA20OX3 revealed four SNPs in the promoter region of BrGA20OX3 between P11 and C634 (Fig. 6B), which might be the cause of different expression levels of BrGA20OX3 between P11 and C634.


Figure 6
View larger version (14K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		Figure 6. Gene expression and SNPs of BrGA20OX3. (A) Semi-quantitative RT–PCR analysis of BrGA20OX3 in the leaves of P11 and C634. Data were normalized to C634; mean ± SE (n = 3). (B) Nucleotide polymorphisms of a promoter region of BrGA20OX3 between P11 and C634. The 5' sequence upstream of translational start site (+1) is represented by a white box. Numbers above the box indicate SNP sites connected by dotted arrows showing nucleotides of P11 and C634 at the left and right of the slashes, respectively.

 
3.4.2. Leaf hairiness
A major QTL for LH (LHQTL1) was detected on A06, which explained the phenotypic variances of 66.8% and 50.2% in the unheated-greenhouse subpopulation and the heated-greenhouse subpopulation, respectively. A minor QTL, LHQTL2, co-located with LDRQTL2, was detected only in the unheated-greenhouse subpopulation. The LH trait has been reported to be conditioned by a major QTL and a minor QTL in B. rapa.7Go Association mapping of leaf trichome detected an associated AFLP marker in A05.57Go BrTTG1, encoding a WD40 repeat protein, TRANSPARENT TESTA GLABRA01, which has been recently cloned from A06 and verified to simultaneously control both hairiness and seed coat color,14Go was located in LHQTL1. However, P11 has hairless leaves and brown seeds, whereas C634 has hairy leaves and yellow seeds.

BrTTG1 expression in leaves of P11 was over twice of that in leaves of C634 (Supplementary Fig. S5), which is not in accord with phenotypic difference that P11 was hairless and C634 was hairy. Nucleotide sequences of the coding region of BrTTG1 in P11 and C634 were determined and aligned together with the sequences of hairy DH lines ‘Y177-12’ and ‘P-8’.14Go There were 11 nucleotide polymorphisms between P11 and C634, which were not specific to P11 in multiple alignments (Supplementary Fig. S6). It can be inferred that BrTTG1 of P11 is not responsible for the leaf hairless trait of P11.

Synteny analysis revealed that GLABRA01 (GL1) (AT3G27920), which is a well-known gene controlling A. thaliana trichome development,59Go,60Go was located just in the syntenic region of a central site of LHQTL1. Blast analysis of GL1 showed that a homologue of GL1 with a maximum nucleotide identity of 90% was present in the B. rapa BAC clone KBrS004A14, which was located in LHQTL1. We determined nucleotide sequences of GL1 alleles of C634 and P11 and found a 5-bp deletion in exon 3 of a P11 allele (Fig. 7). This deletion is considered to be a frameshift mutation causing the loss of the function of the P11 allele.


Figure 7
View larger version (20K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		Figure 7. Identification of a 5-bp deletion in exon 3 of GL1 in P11. Exons 1–3 are indicated by different boxes, and nucleotide sequences around the deletion are shown on the bottom.

 
3.5. Conclusion
As the protein-coding regions of the genomes of Brassica species show high sequence conservation with those of A. thaliana, and the SNPs constitute the most common DNA sequence variations found in genomes of most organisms, EST-based SNP markers contribute to fine comparative mapping of Brassica species and A. thaliana. The present study shows that complex rearrangements and overlaps of the A. thaliana genome exist in the genome of B. rapa; however, genomic collinearity is conserved in many segments of B. rapa, which implies that we can use the positional information from the A. thaliana genome. After QTL mapping in B. rapa, we inferred some candidate genes from the corresponding region in A. thaliana. However, these candidate genes need further analyses, such as complementation or association studies.


   Supplementary data
Top
 Abstract
 1. Introduction
 2. Materials and methods
 3. Results and discussion
 Supplementary data
Funding
 Acknowledgements
 References
 
Supplementary Material is available at www.dnaresearch.oxfordjournals.org.


   Funding
Top
 Abstract
 1. Introduction
 2. Materials and methods
 3. Results and discussion
 Supplementary data
 Funding
Acknowledgements
 References
 
This work was supported in part by a Grant-in-Aid for Scientific Research (A) (19208001) from the Japan Society for the Promotion of Science (JSPS) and the Program for Promotion of Basic and Applied Researches for Innovations in Bio-oriented Industry (BRAIN). F.L. is a recipient of the Monbukagakusho fellowship from the Ministry of Education, Culture, Sports, Science and Technology, Japan.


   Acknowledgements
Top
 Abstract
 1. Introduction
 2. Materials and methods
 3. Results and discussion
 Supplementary data
 Funding
 Acknowledgements
References
 
We are grateful to Dr Kuginuki for providing a plant material.


   Footnotes
 
* To whom correspondence should be addressed. Tel. +81 22-717-8650. Fax. +81 22-717-8654. E-mail: nishio{at}bios.tohoku.ac.jp

Edited by Satoshi Tabata

{dagger} Present address: Tea Research Center, Shizuoka Prefectural Research Institute of Agriculture and Forestry, Kurasawa 1706-11, Kikukawa, Shizuoka 439-0002, Japan. Back


   References
Top
 Abstract
 1. Introduction
 2. Materials and methods
 3. Results and discussion
 Supplementary data
 Funding
 Acknowledgements
 References
 

  1. Lagercrantz U., Lydiate D.J. Comparative genome mapping in Brassica. Genetics (1996) 144:1903–10.[Web of Science][Medline]
  2. Suwabe K., Tsukazaki H., Iketani H., et al. Simple sequence repeat-based comparative genomics between Brassica rapa and Arabidopsis thaliana: the genetic origin of clubroot resistance. Genetics (2006) 173:309–19.[CrossRef][Web of Science][Medline]
  3. Parkin I.A., Gulden S.M., Sharpe A.G., et al. Segmental structure of the Brassica napus genome based on comparative analysis with Arabidopsis thaliana. Genetics (2005) 171:765–81.[CrossRef][Web of Science][Medline]
  4. Cheng X.M., Xu J.S., Xia S., et al. Development and genetic mapping of microsatellite markers from genome survey sequences in Brassica napus. Theor. Appl. Genet. (2009) 118:1121–31.[CrossRef][Web of Science][Medline]
  5. Choi S.R., Teakle G.R., Plaha P., et al. The reference genetic linkage map for the multinational Brassica rapa genome sequencing project. Theor. Appl. Genet. (2007) 115:777–92.[CrossRef][Web of Science][Medline]
  6. Lowe A.J., Moule C., Trick M., Edwards K.J. Efficient large-scale development of microsatellites for marker and mapping applications in Brassica crop species. Theor. Appl. Genet. (2004) 108:1103–12.[CrossRef][Web of Science][Medline]
  7. Song K., Slocum M., Osborn T.C. Molecular marker analysis of genes controlling morphological variation in Brassica rapa (syn. campestris). Theor. Appl. Genet. (1995) 90:1–10.[Web of Science]
  8. Lou P., Zhao J., Kim J.S., et al. Quantitative trait loci for flowering time and morphological traits in multiple populations of Brassica rapa. J. Exp. Bot. (2007) 58:4005–16.[Abstract/Free Full Text]
  9. Tanhuanpää P.K., Vilkki J.P., Vilkki H.J. Mapping of a QTL for oleic acid concentration in spring turnip rape (Brassica rapa ssp. oleifera). Theor. Appl. Genet. (1996) 92:952–6.[CrossRef][Web of Science]
  10. Lou P., Zhao J., He H., et al. Quantitative trait loci for glucosinolate accumulation in Brassica rapa leaves. New Phytol. (2008) 179:1017–32.[CrossRef][Web of Science][Medline]
  11. Zhao J., Jamar D.C., Lou P., et al. Quantitative trait loci analysis of phytate and phosphate concentrations in seeds and leaves of Brassica rapa. Plant Cell Environ. (2008) 31:887–900.[Medline]
  12. Saito M., Kubo N., Matsumoto S., Suwabe K., Tsukada M., Hirai M. Fine mapping of the clubroot resistance gene, Crr3, in Brassica rapa. Theor. Appl. Genet. (2006) 114:81–91.[CrossRef][Web of Science][Medline]
  13. Soengas P., Hand P., Vicente J.G., Pole J.M., Pink D.A. Identification of quantitative trait loci for resistance to Xanthomonas campestris pv. campestris in Brassica rapa. Theor. Appl. Genet. (2007) 114:637–45.[CrossRef][Web of Science][Medline]
  14. Zhang J., Lu Y., Yuan Y., et al. Map-based cloning and characterization of a gene controlling hairiness and seed coat color traits in Brassica rapa. Plant Mol. Biol. (2009) 69:553–63.[CrossRef][Web of Science][Medline]
  15. Kim J.S., Chung T.Y., King G.J., et al. A sequence-tagged linkage map of Brassica rapa. Genetics (2006) 174:29–39.[CrossRef][Web of Science][Medline]
  16. Ching A., Caldwell K.S., Jung M., et al. SNP frequency, haplotype structure and linkage disequilibrium in elite maize inbred lines. BMC Genet. (2002) 3:19.[CrossRef][Medline]
  17. Schneider K., Kulosa D., Soerensen T.R., et al. Analysis of DNA polymorphisms in sugar beet (Beta vulgaris L.) and development of an SNP-based map of expressed genes. Theor. Appl. Genet. (2007) 115:601–15.[CrossRef][Web of Science][Medline]
  18. Trick M., Long Y., Meng J., Bancroft I. Single nucleotide polymorphism (SNP) discovery in the polyploid Brassica napus using solexa transcriptome sequencing. Plant Biotechnol. J. (2009) 7:334–46.[CrossRef][Web of Science][Medline]
  19. Khlestkina E.K., Salina E.A. SNP markers: methods of analysis, ways of development, and comparison on an example of common wheat. Russ. J. Genet. (2006) 42:725–36.
  20. Shirasawa K., Shiokai S., Yamaguchi M., Kishitani S., Nishio T. Dot-blot-SNP analysis for practical plant breeding and cultivar identification in rice. Theor. Appl. Genet. (2006) 113:147–55.[CrossRef][Web of Science][Medline]
  21. Shiokai S., Shirasawa K., Sato Y., Nishio T. Improvement of the dot-blot-SNP technique for efficient and cost-effective genotyping. Mol. Breed. (2009) doi:10.1007/s11032-009-9315-6.
  22. Bao J.S., Corke H., Sun M. Nucleotide diversity in starch synthase IIa and validation of single nucleotide polymorphisms in relation to starch gelatinization temperature and other physicochemical properties in rice (Oryza sativa L.). Theor. Appl. Genet. (2006) 113:1171–83.[CrossRef][Web of Science][Medline]
  23. Achenbach U., Paulo J., Ilarionova E., et al. Using SNP markers to dissect linkage disequilibrium at a major quantitative trait locus for resistance to the potato cyst nematode Globodera pallida on potato chromosome V. Theor. Appl. Genet. (2009) 118:619–29.[CrossRef][Web of Science][Medline]
  24. Kawuki R.S., Ferguson M., Labuschagne M., Herselman L., Kim D.J. Identification, characterisation and application of single nucleotide polymorphisms for diversity assessment in cassava (Manihot esculenta Crantz). Mol. Breed. (2009) 23:669–84.[CrossRef]
  25. Ajisaka H., Kuginuki Y., Yui S., Enomoto S., Hirai M. Identification and mapping of a quantitative trait locus controlling extreme late bolting in Chinese cabbage (Brassica rapa L. ssp. pekinensis syn. campestris L.) using bulked segregant analysis. Euphytica (2001) 118:75–81.[CrossRef][Web of Science]
  26. Nishioka M., Tamura K., Hayashi M., et al. Mapping of QTL for bolting time in Brassica rapa (syn. campestris) under different environmental conditions. Breed. Sci. (2005) 55:127–33.[CrossRef]
  27. Escaravage N., Questiau S., Pornon A., Doche B., Taberlet P. Clonal diversity in a Rhododendron ferrugineum L. (Ericaceae) population inferred from AFLP markers. Mol. Ecol. (1998) 7:975–82.[CrossRef]
  28. Iwata H., Ninomiya S. Ant colony optimization for linkage grouping and locus ordering in genome mapping. (2004) San Diego, USA: Plant and Animal Genome XII.
  29. Wang S., Basten C.J., Zeng Z.B. Windows QTL Cartographer 2.5 (2007) Raleigh, NC: Department of Statistics, North Carolina State University.
  30. De Vicente C., Metz T., Adriana A. Descriptors for genetic markers technologies. (2004) Rome, Italy: IPGRI Report.
  31. Fujimoto R., Sasaki T., Nishio T. Characterization of DNA methyltransferase genes in Brassica rapa. Genes Genet. Syst. (2006) 81:235–42.[CrossRef][Web of Science][Medline]
  32. Parkin I.A., Sharpe A.G., Keith D.J., Lydiate D.J. Identification of the A and C genomes of amphidiploid Brassica napus (oilseed rape). Genome (1995) 38:1122–31.[Medline]
  33. Mun J.H., Kwon S.J., Yang T.J., et al. The first generation of a BAC-based physical map of Brassica rapa. BMC Genomics (2008) 9:280.[CrossRef][Medline]
  34. Cogan N.O., Ponting R.C., Vecchies A.C., et al. Gene-associated single nucleotide polymorphism discovery in perennial ryegrass (Lolium perenne L.). Mol. Genet. Genomics (2006) 276:101–12.[CrossRef][Web of Science][Medline]
  35. Yang Y.W., Lai K.N., Tai P.Y., Li W.H. Rates of nucleotide substitution in angiosperm mitochondrial DNA sequences and dates of divergence between Brassica and other angiosperm lineages. J. Mol. Evol. (1999) 48:597–604.[CrossRef][Web of Science][Medline]
  36. Bevan M., Walsh S. The Arabidopsis genome: a foundation for plant research. Genome Res. (2005) 15:1632–42.[Abstract/Free Full Text]
  37. Qiu D., Gao M., Li G., Quiros C. Comparative sequence analysis for Brassica oleracea with similar sequences in B. rapa and Arabidopsis thaliana. Plant Cell Rep. (2009) 28:649–61.[CrossRef][Web of Science][Medline]
  38. Lagercrantz U. Comparative mapping between Arabidopsis thaliana and Brassica nigra indicates that Brassica genomes have evolved through extensive genome replication accompanied by chromosome fusions and frequent rearrangements. Genetics (1998) 150:1217–28.[Web of Science][Medline]
  39. Lan T.H., DelMonte T.A., Reischmann K.P., et al. An EST-enriched comparative map of Brassica oleracea and Arabidopsis thaliana. Genome Res. (2000) 10:776–88.[Abstract/Free Full Text]
  40. Park J.Y., Koo D.H., Hong C.P., et al. Physical mapping and microsynteny of Brassica rapa ssp. pekinensis genome corresponding to a 222 kbp gene-rich region of Arabidopsis chromosome 4 and partially duplicated on chromosome 5. Mol. Genet. Genomics (2005) 274:579–88.[CrossRef][Web of Science][Medline]
  41. Schranz M.E., Quijada P., Sung S.B., Lukens L., Amasino R., Osborn T.C. Characterization and effects of the replicated flowering time gene FLC in Brassica rapa. Genetics (2002) 162:1457–68.[Web of Science][Medline]
  42. Osborn T.C., Kole C., Parkin I.A., et al. Comparison of flowering time genes in Brassica rapa, B. napus and Arabidopsis thaliana. Genetics (1997) 146:1123–9.[Web of Science][Medline]
  43. Michaels S.D., Amasino R.M. FLOWERING LOCUS C encodes a novel MADS domain protein that acts as a repressor of flowering. Plant Cell (1999) 11:949–56.[Abstract/Free Full Text]
  44. Sheldon C.C., Rouse D.T., Finnegan E.J., Peacock W.J., Dennis E.S. The molecular basis of vernalization: the central role of FLOWERING LOCUS C (FLC). Proc. Natl Acad. Sci. USA (2000) 97:3753–8.[Abstract/Free Full Text]
  45. Okazaki K., Sakamoto K., Kikuchi R., et al. Mapping and characterization of FLC homologs and QTL analysis of flowering time in Brassica oleracea. Theor. Appl. Genet. (2007) 114:595–608.[CrossRef][Web of Science][Medline]
  46. Teutonico R.A., Osborn T.C. Mapping loci controlling vernalization requirement in Brassica rapa. Theor. Appl. Genet. (1995) 91:1279–83.[Web of Science]
  47. Yuan Y.X., Wu J., Sun R.F., et al. A naturally occurring splicing site mutation in the Brassica rapa FLC1 gene is associated with variation in flowering time. J. Exp. Bot. (2009) 60:1299–308.[Abstract/Free Full Text]
  48. Martynov V.V., Tsvetkov I.L., Khavkin E.E. Orthologs of Arabidopsis CLAVATA1 gene in cultivated Brassicaceae plants. Russ. J. Dev. Biol. (2004) 35:41–6.
  49. Clark S.E., Williams R.W., Meyerowitz E.M. The CLAVATA1 gene encodes a putative receptor kinase that controls shoot and floral meristem size in Arabidopsis. Cell (1997) 89:575–85.[CrossRef][Web of Science][Medline]
  50. Dievart A., Dalal M., Tax F.E., et al. CLAVATA1 dominant-negative alleles reveal functional overlap between multiple receptor kinases that regulate meristem and organ development. Plant Cell (2003) 15:1198–211.[Abstract/Free Full Text]
  51. Scortecci K.C., Michaels S.D., Amasino R.M. Identification of a MADS-box gene, FLOWERING LOCUS M, that represses flowering. Plant J. (2001) 26:229–36.[CrossRef][Web of Science][Medline]
  52. Werner J.D., Borevitz J.O., Warthmann N., et al. Quantitative trait locus mapping and DNA array hybridization identify an FLM deletion as a cause for natural flowering-time variation. Proc. Natl Acad. Sci. USA (2005) 102:2460–5.[Abstract/Free Full Text]
  53. Yang T.J., Kim J.S., Kwon S.J., et al. Sequence-level analysis of the diploidization process in the triplicated FLOWERING LOCUS C region of Brassica rapa. Plant Cell (2006) 18:1339–47.[Abstract/Free Full Text]
  54. Lincoln C., Long J., Yamaguchi J., Serikawa K., Hake S. A knotted1-like homeobox gene in Arabidopsis is expressed in the vegetative meristem and dramatically alters leaf morphology when overexpressed in transgenic plants. Plant Cell (1994) 6:1859–76.[Abstract/Free Full Text]
  55. Uchida N., Kimura S., Koenig D., Sinha N. Coordination of leaf development via regulation of KNOX1 genes. J. Plant Res. (2009) doi:10.1007/s10265-009-0248-2.
  56. Byrne M.E., Barley R., Curtis M., et al. Asymmetric leaves1 mediates leaf patterning and stem cell function in Arabidopsis. Nature (2000) 408:967–71.[CrossRef][Medline]
  57. Zhao J., Paulo M.J., Jamar D., et al. Association mapping of leaf traits, flowering time, and phytate content in Brassica rapa. Genome (2007) 50:963–73.[Medline]
  58. Hay A., Kaur H., Phillips A., Hedden P., Hake S., Tsiantis M. The gibberellin pathway mediates KNOTTED1-type homeobox function in plants with different body plans. Curr. Biol. (2002) 12:1557–65.[CrossRef][Web of Science][Medline]
  59. Oppenheimer D.G., Herman P.L., Sivakumaran S., Esch J., Marks M.D. A myb gene required for leaf trichome differentiation in Arabidopsis is expressed in stipules. Cell (1991) 67:483–93.[CrossRef][Web of Science][Medline]
  60. Larkin J.C., Oppenheimer D.G., Marks M.D. The GL1 gene and the trichome developmental pathway in Arabidopsis thaliana. Results Probl. Cell Differ. (1994) 20:259–75.[Medline]

Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?



This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary Data
Right arrowOA All Versions of this Article:
16/6/311    most recent
dsp020v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Request Permissions
Google Scholar
Right arrow Articles by Li, F.
Right arrow Articles by Nishio, T.
PubMed
Right arrow PubMed Citation
Right arrow Articles by Li, F.
Right arrow Articles by Nishio, T.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?