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.     

Genome-wide comparative analysis of NBS-encoding genes between Brassica species and Arabidopsis thaliana

Background

Plant disease resistance (R) genes with the nucleotide binding site (NBS) play an important role in offering resistance to pathogens. The availability of complete genome sequences of Brassica oleracea and Brassica rapa provides an important opportunity for researchers to identify and characterize NBS-encoding R genes in Brassica species and to compare with analogues in Arabidopsis thaliana based on a comparative genomics approach. However, little is known about the evolutionary fate of NBS-encoding genes in the Brassica lineage after split from A. thaliana.

Results

Here we present genome-wide analysis of NBS-encoding genes in B. oleracea, B. rapa and A. thaliana. Through the employment of HMM search and manual curation, we identified 157, 206 and 167 NBS-encoding genes in B. oleracea, B. rapa and A. thaliana genomes, respectively. Phylogenetic analysis among 3 species classified NBS-encoding genes into 6 subgroups. Tandem duplication and whole genome triplication (WGT) analyses revealed that after WGT of the Brassica ancestor, NBS-encoding homologous gene pairs on triplicated regions in Brassica ancestor were deleted or lost quickly, but NBS-encoding genes in Brassica species experienced species-specific gene amplification by tandem duplication after divergence of B. rapa and B. oleracea. Expression profiling of NBS-encoding orthologous gene pairs indicated the differential expression pattern of retained orthologous gene copies in B. oleracea and B. rapa. Furthermore, evolutionary analysis of CNL type NBS-encoding orthologous gene pairs among 3 species suggested that orthologous genes in B. rapa species have undergone stronger negative selection than those in B .oleracea species. But for TNL type, there are no significant differences in the orthologous gene pairs between the two species.

Conclusion

This study is first identification and characterization of NBS-encoding genes in B. rapa and B. oleracea based on whole genome sequences. Through tandem duplication and whole genome triplication analysis in B. oleracea, B. rapa and A. thaliana genomes, our study provides insight into the evolutionary history of NBS-encoding genes after divergence of A. thaliana and the Brassica lineage. These results together with expression pattern analysis of NBS-encoding orthologous genes provide useful resource for functional characterization of these genes and genetic improvement of relevant crops.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-3) contains supplementary material, which is available to authorized users.     

Evidence for homology of flowering-time genes VFR2 from Brassica rapa and FLC from Arabidopsis thaliana

Variation in flowering time is important for the adaptation of plant species to different natural and agricultural environments. We previously identified VFR2 as one of two major quantitative trait loci (QTLs) controlling vernalization-responsive flowering time in a segregating population derived from a cross of annual and biennial Brassica rapa. The region containing VFR2 is homologous to a region in Brassica napus that controls the same trait, and also to a region on chromosome 5 of Arabidopsis thaliana that contains several flowering-time loci. In order to determine precisely the allelic effects and map position of VFR2, we backcrossed the late allele into an early flowering line and obtained monogenic segregation for flowering time in a BC₃S₁ population. The two homozygous genotypic classes differed by 43 and 95 days to flowering in the field and growth chamber, respectively; and the effect of the late allele was almost completely additive. DNA probes that were previously shown to detect RFLP loci in the VFR2 region, or in the homologous regions of B. napus or A. thaliana (including two DNA clones of flowering-time genes) were used to construct a high-resolution map around VFR2. An RFLP detected by an A. thaliana cDNA clone of flowering locus C (FLC) co-segregated exactly with VFR2 in 414 gametes analyzed (<0.24 cM). FLC is a repressor of flowering and is required for the winter-annual habit of late-flowering ecotypes of A. thaliana. The regulation of FLC RNA in B. rapa was consistent with that seen in A. thaliana; RNA levels were up-regulated in the biennial parent and down-regulated by cold treatment. Thus, VFR2 appears to be homologous to FLC and may control flowering time though a similar mechanism as in A. thaliana .     

A Survey of the Brassica rapa genome by BAC-end sequence analysis and comparison with Arabidopsis thaliana

Brassica rapa ssp. pekinensis (Chinese cabbage) is an economically important crop and a model plant for studies on polyploidization and phenotypic evolution. To gain an insight into the structure of the B. rapa genome we analyzed 12,017 BAC-end sequences for the presence of transposable elements (TEs), SSRs, centromeric satellite repeats and genes, and similarity to the closely related genome of Arabidopsis thaliana. TEs were estimated to occupy 14% of the genome, with 12.3% of the genome represented by retrotransposons. It was estimated that the B. rapa genome contains 43,000 genes, 1.6 times greater than the genome of A. thaliana. A number of centromeric satellite sequences, representing variations of a 176-bp consensus sequence, were identified. This sequence has undergone rapid evolution within the B. rapa genome and has diverged among the related species of Brassicaceae. A study of SSRs demonstrated a non-random distribution with a greater abundance within predicted intergenic regions. Our results provide an initial characterization of the genome of B. rapa and provide the basis for detailed analysis through whole-genome sequencing.     

Comparative sequence analysis for Brassica oleracea with similar sequences in B. rapa and Arabidopsis thaliana

We sequenced five BAC clones of Brassica oleracea doubled haploid ‘Early Big' broccoli containing major genes in the aliphatic glucosinolate pathway, and comparatively analyzed them with similar sequences in A. thaliana and B. rapa. Additionally, we included in the analysis published sequences from three other B. oleracea BAC clones and a contig of this species corresponding to segments in A. thaliana chromosomes IV and V. A total of 2,946 kb of B. oleracea, 1,069 kb of B. rapa sequence and 2,607 kb of A. thaliana sequence were compared and analyzed. We found conserved collinearity for gene order and content restricted to specific chromosomal segments, but breaks in collinearity were frequent resulting in gene absence likely not due to gene loss but rearrangements. B. oleracea has the lowest gene density of the three species, followed by B. rapa. The genome expansion of the Brassica species, B. oleracea in particular, is due to larger introns and gene spacers resulting from frequent insertion of DNA transposons and retrotransposons. These findings are discussed in relation to the possible origin and evolution of the Brassica genomes. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Linkage mapping of genes controlling resistance to white rust (Albugo candida) in Brassica rapa (syn. campestris) and comparative mapping to Brassica napus and Arabidopsis thaliana.

Genes for resistance to white rust (Albugo candida) in oilseed Brassica rapa were mapped using a recombinant inbred (RI) population and a genetic linkage map consisting of 144 restriction fragment length polymorphism (RFLP) markers and 3 phenotypic markers. Young seedlings were evaluated by inoculating cotyledons with A. candida race 2 (AC2) and race 7 (AC7) and scoring the interaction phenotype (IP) on a 0-9 scale. The IP of each line was nearly identical for the two races and the population showed bimodal distributions, suggesting that a single major gene (or tightly linked genes) controlled resistance to the two races. The IP scores were converted to categorical resistant and susceptible scores, and these data were used to map a single Mendelian gene controlling resistance to both races on linkage group 4 where resistance to race 2 had been mapped previously. A quantitative trait loci (QTL) mapping approach using the IP scores detected the same major resistance locus for both races, plus a second minor QTL effect for AC2 on linkage group 2. These results indicate that either a dominant allele at a single locus (Acal) or two tightly linked loci control seedling resistance to both races of white rust in the biennial turnip rape cultivar Per. The map positions of white rust resistance genes in B. rapa and Brassica napus were compared and the results indicate where additional loci that have not been mapped may be located. Alignment of these maps to the physical map of the Arabidopsis genome identified regions to target for comparative fine mapping using this model organism.     

Molecular mapping of Arabidopsis thaliana lipid-related orthologous genes in Brassica napus

Quantitative Trait Loci (QTL) for oil content has been previously analyzed in a SG-DH population from a cross between a Chinese cultivar and a European cultivar of Brassica napus. Eight QTL with additive and epistatic effects, and with environmental interactions were evaluated. Here we present an integrated linkage map of this population predominantly based on informative markers derived from Brassica sequences, including 249 orthologous A. thaliana genes, where nearly half (112) are acyl lipid metabolism related genes. Comparative genomic analysis between B. napus and A. thaliana revealed 33 colinearity regions. Each of the conserved A. thaliana segments is present two to six?times in the B. napus genome. Approximately half of the mapped lipid-related orthologous gene loci (76/137) were assigned in these conserved colinearity regions. QTL analysis for seed oil content was performed using the new map and phenotypic data from 11 different field trials. Nine significant QTL were identified on linkage groups A1, A5, A7, A9, C2, C3, C6 and C8, together explaining 57.79% of the total phenotypic variation. A total of 14 lipid related candidate gene loci were located in the confidence intervals of six of these QTL, of which ten were assigned in the conserved colinearity regions and felled in the most frequently overlapped QTL intervals. The information obtained from this study demonstrates the potential role of the suggested candidate genes in rapeseed kernel oil accumulation.     

Anthocyanin biosynthetic genes in Brassica rapa

Background

Anthocyanins are a group of flavonoid compounds. As a group of important secondary metabolites, they perform several key biological functions in plants. Anthocyanins also play beneficial health roles as potentially protective factors against cancer and heart disease. To elucidate the anthocyanin biosynthetic pathway in Brassica rapa, we conducted comparative genomic analyses between Arabidopsis thaliana and B. rapa on a genome-wide level.

Results

In total, we identified 73 genes in B. rapa as orthologs of 41 anthocyanin biosynthetic genes in A. thaliana. In B. rapa, the anthocyanin biosynthetic genes (ABGs) have expanded and most genes exist in more than one copy. The anthocyanin biosynthetic structural genes have expanded through whole genome and tandem duplication in B. rapa. More structural genes located upstream of the anthocyanin biosynthetic pathway have been retained than downstream. More negative regulatory genes are retained in the anthocyanin biosynthesis regulatory system of B. rapa.

Conclusions

These results will promote an understanding of the genetic mechanism of anthocyanin biosynthesis, as well as help the improvement of the nutritional quality of B. rapa through the breeding of high anthocyanin content varieties.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-426) contains supplementary material, which is available to authorized users.     

Glucosinolate biosynthetic genes in Brassica rapa

Glucosinolates (GS) are a group of amino acid-derived secondary metabolites found throughout the Cruciferae family. Glucosinolates and their degradation products play important roles in pathogen and insect interactions, as well as in human health. In order to elucidate the glucosinolate biosynthetic pathway in Brassica rapa, we conducted comparative genomic analyses of Arabidopsis thaliana and B. rapa on a genome-wide level. We identified 102 putative genes in B. rapa as the orthologs of 52 GS genes in A. thaliana. All but one gene was successfully mapped on 10 chromosomes. Most GS genes exist in more than one copy in B. rapa. A high co-linearity in the glucosinolate biosynthetic pathway between A. thaliana and B. rapa was also established. The homologous GS genes in B. rapa and A. thaliana share 59-91% nucleotide sequence identity and 93% of the GS genes exhibit synteny between B. rapa and A. thaliana. Moreover, the structure and arrangement of the B. rapa GS (BrGS) genes correspond with the known evolutionary divergence of B. rapa, and may help explain the profiles and accumulation of GS in B. rapa.     

Aspartic proteinase genes in the Brassicaceae Arabidopsis thaliana and Brassica napus

Active aspartic proteinase is isolated from Brassica napus seeds and the peptide sequence is used to generate primers for PCR. We present here cDNA and genomic clones for aspartic proteinases from the closely related Brassicaceae Arabidopsis thaliana and Brassica napus. The Arabidopsis cDNA represents a single gene, while Brassica has at least 4 genes. Like other plant aspartic proteases, the two Brassicaceae enzymes contain an extra protein domain of about 100 amino acids relative to the mammalian forms. The intron/exon arrangement in the Brassica genomic clone is significantly different from that in mammalian genes. As the proteinase is isolated from seeds, the same tissue where 2S albumins are processed, this implies expression of one of the aspartic proteinase genes there.     

Pong-like elements in Arabidopsis and Brassica rapa: its regulation of F-box protein gene in different ecotypes of Arabidopsis thaliana

Pong elements are active class 2 transposons in rice. We found Pong-like elements in Arabidopsis thaliana and named them as AtPong (Arabidopsis thaliana Pong) elements. The AtPong elements carried 13 bp of TIRs and two ORFs in which NAM DNA binding domain and DDE catalytic domain are encoded. Ping and mPing (miniature Ping) elements are deficient Pong elements and we found AtPong derived deficient AtPing and AtmPing elements. We also found a homologous element of the AtPong element in Brassica rapa. This element was a deficient element by internal deletion, but showed high sequence similarity with the AtmPing so that it was named as BrmPing (Brassica rapa miniature Ping). The AtPong, AtPing, and AtmPing elements were present in intergenic regions except one element, AtPing1, which was present in an exon of a F-box protein gene. Among the different A. thaliana ecotypes, the AtPing1 showed polymorphisms of present and absent in the F-box protein gene. The excision of the AtPing1 restored the expression of the F-box protein gene, indicating that the expression of the F-box protein genes is regulated by the AtPing1.     

Segmental structure of the Brassica napus genome based on comparative analysis with Arabidopsis thaliana

Over 1000 genetically linked RFLP loci in Brassica napus were mapped to homologous positions in the Arabidopsis genome on the basis of sequence similarity. Blocks of genetically linked loci in B. napus frequently corresponded to physically linked markers in Arabidopsis. This comparative analysis allowed the identification of a minimum of 21 conserved genomic units within the Arabidopsis genome, which can be duplicated and rearranged to generate the present-day B. napus genome. The conserved regions extended over lengths as great as 50 cM in the B. napus genetic map, equivalent to approximately 9 Mb of contiguous sequence in the Arabidopsis genome. There was also evidence for conservation of chromosome landmarks, particularly centromeric regions, between the two species. The observed segmental structure of the Brassica genome strongly suggests that the extant Brassica diploid species evolved from a hexaploid ancestor. The comparative map assists in exploiting the Arabidopsis genomic sequence for marker and candidate gene identification within the larger, intractable genomes of the Brassica polyploids.     

Characterization of DNA methyltransferase genes in Brassica rapa

DNA methylation is essential for normal development and plays important roles in regulating gene expression in plants. Analysis of the key enzymes catalyzing DNA methylation is important to understand epigenetic phenomena. In this study, three putative methyltransferase genes, BrMET1a, BrMET1b, and BrCMT, were isolated from a genome library of Brassica rapa. Structural conservation of the amino acid sequence between BrMET1a/BrMET1b and AtMET1 and that between BrCMT and AtCMT3 suggests that they may function as DNA methyltransferase. BrMET1a was expressed in vegetative and reproductive organs, while BrMET1b was expressed only in pistils, indicating that these two genes have different functions. BrCMT was expressed especially in stamens at the stage of 2-4 days before anthesis. We isolated three DNA methyltransferase genes in Brassica rapa and indicated differences of expression patterns of these DNA methyltransferase genes and expression levels in different tissues and developmental stages, suggesting that these genes might play important roles in epigenetic gene regulation in B. rapa.     

Metallochaperone-like genes in Arabidopsis thaliana

A complete inventory of metallochaperone-like proteins containing a predicted HMA domain in Arabidopsis revealed a large family of 67 proteins. 45 proteins, the HIPPs, have a predicted isoprenylation site while 22 proteins, the HPPs, do not. Sequence comparisons divided the proteins into seven major clusters (I-VII). Cluster IV is notable for the presence of a conserved Asp residue before the CysXXCys, metal binding motif, analogous to the Zn binding motif in E. coli ZntA. HIPP20, HIPP21, HIPP22, HIPP26 and HIPP27 in Cluster IV were studied in more detail. All but HIPP21 could rescue the Cd-sensitive, ycf1 yeast mutant but failed to rescue the growth of zrt1zrt2, zrc1cot1 and atx1 mutants. In Arabidopsis, single and double mutants did not show a phenotype but the hipp20/21/22 triple mutant was more sensitive to Cd and accumulated less Cd than the wild-type suggesting the HIPPs can have a role in Cd-detoxification, possibly by binding Cd. Promoter-GUS reporter expression studies indicated variable expression of these HIPPs. For example, in roots, HIPP22 and HIPP26 are only expressed in lateral root tips while HIPP20 and HIPP25 show strong expression in the root vasculature.     

Arabidopsis thaliana: a source of candidate disease-resistance genes for Brassica napus.

Common structural and amino acid motifs among cloned plant disease-resistance genes (R genes), have made it possible to identify putative disease-resistance sequences based on DNA sequence identity. Mapping of such R-gene homologues will identify candidate disease-resistance loci to expedite map-based cloning strategies in complex crop genomes. Arabidopsis thaliana expressed sequence tags (ESTs) with homology to cloned plant R genes (R-ESTs), were mapped in both A. thaliana and Brassica napus to identify candidate R-gene loci and investigate intergenomic collinearity. Brassica R-gene homologous sequences were also mapped in B. napus. In total, 103 R-EST loci and 36 Brassica R-gene homologous loci were positioned on the N-fo-61-9 B. napus genetic map, and 48 R-EST loci positioned on the Columbia x Landsberg A. thaliana map. The mapped loci identified collinear regions between Arabidopsis and Brassica which had been observed in previous comparative mapping studies; the detection of syntenic genomic regions indicated that there was no apparent rapid divergence of the identified genomic regions housing the R-EST loci.     

Isolation and characterization of microsatellites in Brassica rapa L

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 F(2) population from a cross of doubled-haploid lines DH27 x 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.     

Gravity control of growth form in Brassica rapa and Arabidopsis thaliana (Brassicaceae): Consequences for secondary metabolism

How gravity influences the growth form and flavor components of plants is of interest to the space program because plants could be used for food and life support during prolonged missions away from the planet, where that constant feature of Earth's environment does not prevail. We used plant growth hardware from prior experiments on the space shuttle to grow Brassica rapa and Arabidopsis thaliana plants during 16-d or 11-d hypergravity treatments on large-diameter centrifuge rotors. Both species showed radical changes in growth form, becoming more prostrate with increasing g-loads (2-g and 4-g). In Brassica, height decreased and stems thickened in a linear relationship with increasing g-load. Glucosinolates, secondary compounds that contribute flavor to Brassica, decreased by 140% over the range of micro to 4-g, while the structural secondary compound, lignin, remained constant at ～15% (w/w) cell wall dry mass. Stem thickening at 4-g was associated with substantial increases in cell size (47%, 226%, and 33% for pith, cortex, and vascular tissue), rather than any change in cell number. The results, which demonstrate the profound effect of gravity on plant growth form and secondary metabolism, are discussed in the context of similar thigmostresses such as touch and wind.     

Mutants of circadian-associated PRR genes display a novel and visible phenotype as to light responses during de-etiolation of Arabidopsis thaliana seedlings

In Arabidopsis thaliana, it is currently accepted that certain mutants with lesions in clock-associated genes commonly display hallmarked phenotypes with regard to three characteristic biological events: (i) altered rhythmic expression of circadian-controlled genes, (ii) changes in flowering time, and (iii) altered sensitivity to red light in elongation of hypocotyls. During the course of examination of the clock-associated mutants of PSEUDO-RESPONSE REGULATORS, PRRs, including TOC1 (PRR1), we found that they commonly show another visible phenotype of anomalous greening responses upon the onset to light exposure of etiolated seedlings. These findings are indicative of a novel link between circadian rhythms and chloroplast development.