甘蓝型油菜BnFLA基因的克隆及表达分析

利用同源克隆法从甘蓝型油菜中获得了1个类成束阿拉伯半乳聚糖蛋白基因(FLA),命名为BnFLA。BnFLA基因开放阅读框长为1 200bp,编码399个氨基酸,分子量为42 885.9Da,等电点为6.37。预测的BnFLA蛋白包含N-端信号肽、2个AGP-like结构域、2个fasciclin-like结构域和C-端GPI-anchor序列。系统进化分析表明BnFLA氨基酸序列与BrFLA17和AtFLA2进化关系较近,一致性分别为98%和87%。qRT-PCR分析表明,BnFLA基因在油菜各组织均有表达,并以下胚轴中表达量最高,其次为子叶,茎秆中表达最少;BnFLA基因的表达受到GA3、BR、IAA、ABA和NaCl的诱导,但受6-BA、蔗糖、低温和PEG抑制。研究认为,油菜中BnFLA基因可能参与激素信号转导途径和非生物胁迫应答。     

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.     

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.     

<Emphasis Type="Italic">BcMF11</Emphasis> and its homologous sequences may form a lncRNA family in <Emphasis Type="Italic">Brassica</Emphasis> diploids

BcMF11 is a long non-coding RNA that has been identified in Brassica rapa and shown to be involved in pollen development. Here, when re-cloned the gene sequence, multiple paralogous copies of BcMF11 were identified in B. rapa (A genome). Multiple paralogous copies of BcMF11 were also found in B. nigra (B genome) and Brassica oleracea (C genome), the other two primary diploids of Brassica U triangle. While in the early diverging Brassicaceae lineage including Arabidopsis thaliana, no BcMF11 homolog was found. Phylogenetic analysis showed that the BcMF11 homologous sequences cloned from A genome or C genome could be clustered into a separate branch, respectively. However, there was no distinct cluster defined for BcMF11 homologous sequences cloned from B genome. The expression of BcMF11 in B. rapa was investigated and revealed a different result in the previous study. In addition, 12 expressed sequence tags from B. napus and B. rapa showing high similarities with BcMF11 were identified in the NCBI database, which further verified that rather than the useless repeat fragments in the genome, the BcMF11 homologous genes could transcribe. It is possible that BcMF11 and its homologous sequences may form a large gene family which might be originated in the recent ancestral lineage of Brassica.     

The putative phytocyanin genes in Chinese cabbage (Brassica rapa L.): genome-wide identification, classification and expression analysis

Phytocyanins (PCs) are a plant-specific family of small copper-containing electron transfer proteins. PCs may bind with a single copper atom to function as electron transporters in various biological systems, such as copper trafficking and plant photosynthesis. Evidence indicates that PCs may also be involved in plant developmental processes and stress responses. Many PCs possess arabinogalactan protein-like regions and are therefore termed chimeric arabinogalactan proteins (CAGPs). Previously, 38 and 62 PC genes have been identified in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), respectively. The recent release of the Chinese cabbage genome (B. rapa ssp. Pekinensis line Chiifu-401-42) enabled us to perform a genome-wide identification and analysis. In this study we identified 84 putative PC genes in the B. rapa genome. All of the Brassica rapa phytocyanins (BrPCs) described here could be divided, based on motif constitution, into the following three main subclasses: 52 early nodulin-like proteins (ENODLs), 16 uclacyanin-like proteins (UCLs), and 11 stellacyanin-like proteins (SCLs). A structural analysis predicted that 71 BrPCs contained N-terminal secretion signals and 45 BrPCs may be glycosylphosphatidylinositol-anchored to the plasma membrane. Glycosylation prediction revealed that 48 BrPCs were CAGPs with putative arabinogalactan glycomodules, and 57 BrPCs had N-glycosylation sites. Additionally, gene duplication analysis demonstrated that almost all of the duplicated BrPC genes shared the same conserved collinear blocks and that segmental duplications play an important role in the diversification of this gene family. Surprisingly, all BrUCL genes were duplicated except for BrUCL16. Expression analyses indicated that BrENODL22/27 and BrSCL8/9 were highly expressed in reproductive organs; BrUCL6/16 was strongly expressed in roots and even more strongly expressed in stems. The genome-wide identification, classification and expression analysis of BrPCs will provide a fundamental basis for the evolution and modification of the gene family after a polyploidy event and enable the functional study of PC genes in a polyploid crop species.     

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

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.     

Carotenoid biosynthetic genes in Brassica rapa: comparative genomic analysis,phylogenetic analysis,and expression profiling

Background

Carotenoids are isoprenoid compounds synthesized by all photosynthetic organisms. Despite much research on carotenoid biosynthesis in the model plant Arabidopsis thaliana, there is a lack of information on the carotenoid pathway in Brassica rapa. To better understand its carotenoid biosynthetic pathway, we performed a systematic analysis of carotenoid biosynthetic genes at the genome level in B. rapa.

Results

We identified 67 carotenoid biosynthetic genes in B. rapa, which were orthologs of the 47 carotenoid genes in A. thaliana. A high level of synteny was observed for carotenoid biosynthetic genes between A. thaliana and B. rapa. Out of 47 carotenoid biosynthetic genes in A. thaliana, 46 were successfully mapped to the 10 B. rapa chromosomes, and most of the genes retained more than one copy in B. rapa. The gene expansion was caused by the whole-genome triplication (WGT) event experienced by Brassica species. An expression analysis of the carotenoid biosynthetic genes suggested that their expression levels differed in root, stem, leaf, flower, callus, and silique tissues. Additionally, the paralogs of each carotenoid biosynthetic gene, which were generated from the WGT in B. rapa, showed significantly different expression levels among tissues, suggesting differentiated functions for these multi-copy genes in the carotenoid pathway.

Conclusions

This first systematic study of carotenoid biosynthetic genes in B. rapa provides insights into the carotenoid metabolic mechanisms of Brassica crops. In addition, a better understanding of carotenoid biosynthetic genes in B. rapa will contribute to the development of conventional and transgenic B. rapa cultivars with enriched carotenoid levels in the future.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1655-5) contains supplementary material, which is available to authorized users.     

Identification of Novel QTLs for Isolate-Specific Partial Resistance to Plasmodiophora brassicae in Brassica rapa

Plasmodiophora brassicae, the causal agent of clubroot disease of the Brassica crops, is widespread in the world. Quantitative trait loci (QTLs) for partial resistance to 4 different isolates of P. brassicae (Pb2, Pb4, Pb7, and Pb10) were investigated using a BC₁F₁ population from a cross between two subspecies of Brassica rapa, i.e. Chinese cabbage inbred line C59-1 as a susceptible recurrent parent and turnip inbred line ECD04 as a resistant donor parent. The BC₁F₂ families were assessed for resistance under controlled conditions. A linkage map constructed with simple sequence repeats (SSR), unigene-derived microsatellite (UGMS) markers, and specific markers linked to published clubroot resistance (CR) genes of B. rapa was used to perform QTL mapping. A total of 6 QTLs residing in 5 CR QTL regions of the B. rapa chromosomes A01, A03, and A08 were identified to account for 12.2 to 35.2% of the phenotypic variance. Two QTL regions were found to be novel except for 3 QTLs in the respective regions of previously identified Crr1, Crr2, and Crr3. QTL mapping results indicated that 1 QTL region was common for partial resistance to the 2 isolates of Pb2 and Pb7, whereas the others were specific for each isolate. Additionally, synteny analysis between B. rapa and Arabidopsis thaliana revealed that all CR QTL regions were aligned to a single conserved crucifer blocks (U, F, and R) on 3 Arabidopsis chromosomes where 2 CR QTLs were detected in A. thaliana. These results suggest that some common ancestral genomic regions were involved in the evolution of CR genes in B. rapa.     

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 .     

Sequence and structure of Brassica rapa chromosome A3

Background

The species Brassica rapa includes important vegetable and oil crops. It also serves as an excellent model system to study polyploidy-related genome evolution because of its paleohexaploid ancestry and its close evolutionary relationships with Arabidopsis thaliana and other Brassica species with larger genomes. Therefore, its genome sequence will be used to accelerate both basic research on genome evolution and applied research across the cultivated Brassica species.

Results

We have determined and analyzed the sequence of B. rapa chromosome A3. We obtained 31.9 Mb of sequences, organized into nine contigs, which incorporated 348 overlapping BAC clones. Annotation revealed 7,058 protein-coding genes, with an average gene density of 4.6 kb per gene. Analysis of chromosome collinearity with the A. thaliana genome identified conserved synteny blocks encompassing the whole of the B. rapa chromosome A3 and sections of four A. thaliana chromosomes. The frequency of tandem duplication of genes differed between the conserved genome segments in B. rapa and A. thaliana, indicating differential rates of occurrence/retention of such duplicate copies of genes. Analysis of 'ancestral karyotype' genome building blocks enabled the development of a hypothetical model for the derivation of the B. rapa chromosome A3.

Conclusions

We report the near-complete chromosome sequence from a dicotyledonous crop species. This provides an example of the complexity of genome evolution following polyploidy. The high degree of contiguity afforded by the clone-by-clone approach provides a benchmark for the performance of whole genome shotgun approaches presently being applied in B. rapa and other species with complex genomes.     

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.     

Chloroplast and mitochondrial SSR help to distinguish allo-cytoplasmic male sterile types in cabbage (Brassica oleracea L. var. capitata)

The profiles of single sequence repeat (SSR) in six distinct allo-cytoplasmic male sterile (CMS) types of cabbage (Brassica oleracea L. var. capitata) were generated using 32 SSR primer pairs derived from the Arabidopsis thaliana chloroplast (cp) genome and another 21 SSR primers from the B. napus mitochondrial (mt) genome sequences. In total, 11 cpSSR and 4 mtSSR primers revealed polymorphism among the six cabbage CMS types, namely NigCMS, OguCMSR₁, OguCMSR₂, OguCMSR₃, OguCMSHY and PolCMS. Through cluster analysis, six cabbage CMS types could be unambiguously differentiated with just three sets of primers (ACP43, ACP47, mtSSR2). Analysis of the selected amplicon sequences showed high identity to that of the corresponding sequences in A. thaliana, B. rapa and B. napus. The aligned cluster analysis revealed that the polymorphism mainly included SSR number variation, single nucleotide polymorphism (SNP), and sequence insertion or deletion (InDel). Our results demonstrated that specific mitochondrial or chloroplast SSR analysis could be a feasible alternative means for cabbage CMS type identification.     

Comparative Analysis between Homoeologous Genome Segments of Brassica napus and Its Progenitor Species Reveals Extensive Sequence-Level Divergence

Homoeologous regions of Brassica genomes were analyzed at the sequence level. These represent segments of the Brassica A genome as found in Brassica rapa and Brassica napus and the corresponding segments of the Brassica C genome as found in Brassica oleracea and B. napus. Analysis of synonymous base substitution rates within modeled genes revealed a relatively broad range of times (0.12 to 1.37 million years ago) since the divergence of orthologous genome segments as represented in B. napus and the diploid species. Similar, and consistent, ranges were also identified for single nucleotide polymorphism and insertion-deletion variation. Genes conserved across the Brassica genomes and the homoeologous segments of the genome of Arabidopsis thaliana showed almost perfect collinearity. Numerous examples of apparent transduplication of gene fragments, as previously reported in B. oleracea, were observed in B. rapa and B. napus, indicating that this phenomenon is widespread in Brassica species. In the majority of the regions studied, the C genome segments were expanded in size relative to their A genome counterparts. The considerable variation that we observed, even between the different versions of the same Brassica genome, for gene fragments and annotated putative genes suggest that the concept of the pan-genome might be particularly appropriate when considering Brassica genomes.     

Sequence‐level comparative analysis of the Brassica napus genome around two stearoyl‐ACP desaturase loci

We conducted a sequence‐level comparative analyses, at the scale of complete bacterial artificial chromosome (BAC) clones, between the genome of the most economically important Brassica species, Brassica napus (oilseed rape), and those of Brassica rapa, the genome of which is currently being sequenced, and Arabidopsis thaliana. We constructed a new B. napus BAC library and identified and sequenced clones that contain homoeologous regions of the genome including stearoyl‐ACP desaturase‐encoding genes. We sequenced the orthologous region of the genome of B. rapa and conducted comparative analyses between the Brassica sequences and those of the orthologous region of the genome of A. thaliana. The proportion of genes conserved (～56%) is lower than has been reported previously between A. thaliana and Brassica (～66%). The gene models for sets of conserved genes were used to determine the extent of nucleotide conservation of coding regions. This was found to be 84.2 ± 3.9% and 85.8 ± 3.7% between the B. napus A and C genomes, respectively, and that of A. thaliana, which is consistent with previous results for other Brassica species, and 97.5 ± 3.1% between the B. napus A genome and B. rapa, and 93.1 ± 4.9% between the B. napus C genome and B. rapa. The divergence of the B. napus genes from the A genome and the B. rapa genes was greater than anticipated and indicates that the A genome ancestor of the B. napus cultivar studied was relatively distantly related to the cultivar of B. rapa selected for genome sequencing.     

Extensive Chromosome Homoeology among Brassiceae Species Were Revealed by Comparative Genetic Mapping with High-Density EST-Based SNP Markers in Radish (Raphanus sativus L.)

A linkage map of expressed sequence tag (EST)-based markers in radish (Raphanus sativus L.) was constructed using a low-cost and high-efficiency single-nucleotide polymorphism (SNP) genotyping method named multiplex polymerase chain reaction–mixed probe dot-blot analysis developed in this study. Seven hundred and forty-six SNP markers derived from EST sequences of R. sativus were assigned to nine linkage groups with a total length of 806.7 cM. By BLASTN, 726 markers were found to have homologous genes in Arabidopsis thaliana, and 72 syntenic regions, which have great potential for utilizing genomic information of the model species A. thaliana in basic and applied genetics of R. sativus, were identified. By construction and analysis of the genome structures of R. sativus based on the 24 genomic blocks within the Brassicaceae ancestral karyotype, 23 of the 24 genomic blocks were detected in the genome of R. sativus, and half of them were found to be triplicated. Comparison of the genome structure of R. sativus with those of the A, B, and C genomes of Brassica species and that of Sinapis alba L. revealed extensive chromosome homoeology among Brassiceae species, which would facilitate transfer of the genomic information from one Brassiceae species to another.