The reference genetic linkage map for the multinational Brassica rapa genome sequencing project

We describe the construction of a reference genetic linkage map for the Brassica A genome, which will form the backbone for anchoring sequence contigs for the Multinational Brassica rapa Genome Sequencing Project. Seventy-eight doubled haploid lines derived from anther culture of the F₁ of a cross between two diverse Chinese cabbage (B. rapa ssp. pekinensis) inbred lines, ‘Chiifu-401-42’ (C) and ‘Kenshin-402-43’ (K) were used to construct the map. The map comprises a total of 556 markers, including 278 AFLP, 235 SSR, 25 RAPD and 18 ESTP, STS and CAPS markers. Ten linkage groups were identified and designated as R1–R10 through alignment and orientation using SSR markers in common with existing B. napus reference linkage maps. The total length of the linkage map was 1,182 cM with an average interval of 2.83 cM between adjacent loci. The length of linkage groups ranged from 81 to 161 cM for R04 and R06, respectively. The use of 235 SSR markers allowed us to align the A-genome chromosomes of B. napus with those of B. rapa ssp. pekinensis. The development of this map is vital to the integration of genome sequence and genetic information and will enable the international research community to share resources and data for the improvement of B. rapa and other cultivated Brassica species. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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.     

Integration of Brassica A genome genetic linkage map between Brassica napus and B. rapa.

An integrated linkage map between B. napus and B. rapa was constructed based on a total of 44 common markers comprising 41 SSR (33 BRMS, 6 Saskatoon, and 2 BBSRC) and 3 SNP/indel markers. Between 3 and 7 common markers were mapped onto each of the linkage groups A1 to A10. The position and order of most common markers revealed a high level of colinearity between species, although two small regions on A4, A5, and A10 revealed apparent local inversions between them. These results indicate that the A genome of Brassica has retained a high degree of colinearity between species, despite each species having evolved independently after the integration of the A and C genomes in the amphidiploid state. Our results provide a genetic integration of the Brassica A genome between B. napus and B. rapa. As the analysis employed sequence-based molecular markers, the information will accelerate the exploitation of the B. rapa genome sequence for the improvement of oilseed rape.     

A near-complete genome assembly of Brassica rapa provides new insights into the evolution of centromeres

Brassica rapa comprises many important cultivated vegetables and oil crops. However, Chiifu v3.0, the current B. rapa reference genome, still contains hundreds of gaps. Here, we presented a near-complete genome assembly of B. rapa Chiifu v4.0, which was 424.59 Mb with only two gaps, using Oxford Nanopore Technology (ONT) ultralong-read sequencing and Hi-C technologies. The new assembly contains 12 contigs, with a contig N50 of 38.26 Mb. Eight of the ten chromosomes were entirely reconstructed in a single contig from telomere to telomere. We found that the centromeres were mainly invaded by ALE and CRM long terminal repeats (LTRs). Moreover, there is a high divergence of centromere length and sequence among B. rapa genomes. We further found that centromeres are enriched for Copia invaded at 0.14 MYA on average, while pericentromeres are enriched for Gypsy LTRs invaded at 0.51 MYA on average. These results indicated the different invasion mechanisms of LTRs between the two structures. In addition, a novel repetitive sequence PCR630 was identified in the pericentromeres of B. rapa. Overall, the near-complete genome assembly, B. rapa Chiifu v4.0, offers valuable tools for genomic and genetic studies of Brassica species and provides new insights into the evolution of centromeres.     

白菜型油菜自交亲和性变异分析

对不同来源的85份白菜型油菜品种的自交亲和指数研究表明,白菜型油菜的自交亲和性存在较大幅度的变异,自交亲和指数在0．00～9．28之间,这种变异不仅存在于品种间,而且存在于品种内不同个体间。85份参试材料中,自交亲和指数小于1的有50个品种;35个品种自交亲和指数大于1,自交亲和指数在1．00～9．28之间,其中亲和指数大于4的材料有28、44、45、55x、75、119、123、139号品种等。依据自交亲和指数,可将参试材料分为4种类型:高自交亲和类型(自交亲和指数>4．00);自交亲和类型(自交亲和指数1．00～3．99);自交不亲和类型(0．00< 自交亲和指数<1．00);高自交不亲和类型(自交亲和指数=0．00)。自交亲和性因地区而异,一般西部地区的品种自交亲和性高于其它地区的品种。     

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.     

Intersubgenomic heterosis in seed yield potential observed in a new type of Brassica napus introgressed with partial Brassica rapa genome

This paper reports the observation on the intersubgenomic heterosis for seed yield among hybrids between natural Brassica napus (AⁿAⁿCⁿCⁿ) and a new type of B. napus with introgressions of genomic components of Brassica rapa (A^rA^r). This B. napus was selected from the progeny of B. napus × B. rapa and (B. napus × B. rapa) × B. rapa based on extensive phenotypic and cytological observation. Among the 129 studied partial intersubgenomic hybrids, which were obtained by randomly crossing 13 lines of the new type of B. napus in F₃ or BC₁F₃ to 27 cultivars of B. napus from different regions as tester lines, about 90% of combinations exceeded the yield of their respective tester lines, whereas about 75% and 25% of combinations surpassed two elite Chinese cultivars, respectively. This strong heterosis was further confirmed by reevaluating 2 out of the 129 combinations in a successive year and by surveying hybrids between 20 lines of the new type of B. napus in BC₁F₅ and its parental B. napus in two locations. Some DNA segments from B. rapa were identified with significant effects on seed yield and yield components of the new type of B. napus in BC₁F₅ and intersubgenomic hybrids in positive or negative direction. It seems that the genomic components introgressed from B. rapa contributed to improvement of seed yield of rapeseed.     

Genome-wide comparative analysis of the Brassica rapa gene space reveals genome shrinkage and differential loss of duplicated genes after whole genome triplication

Background

Brassica rapa is one of the most economically important vegetable crops worldwide. Owing to its agronomic importance and phylogenetic position, B. rapa provides a crucial reference to understand polyploidy-related crop genome evolution. The high degree of sequence identity and remarkably conserved genome structure between Arabidopsis and Brassica genomes enables comparative tiling sequencing using Arabidopsis sequences as references to select the counterpart regions in B. rapa, which is a strong challenge of structural and comparative crop genomics.

Results

We assembled 65.8 megabase-pairs of non-redundant euchromatic sequence of B. rapa and compared this sequence to the Arabidopsis genome to investigate chromosomal relationships, macrosynteny blocks, and microsynteny within blocks. The triplicated B. rapa genome contains only approximately twice the number of genes as in Arabidopsis because of genome shrinkage. Genome comparisons suggest that B. rapa has a distinct organization of ancestral genome blocks as a result of recent whole genome triplication followed by a unique diploidization process. A lack of the most recent whole genome duplication (3R) event in the B. rapa genome, atypical of other Brassica genomes, may account for the emergence of B. rapa from the Brassica progenitor around 8 million years ago.

Conclusions

This work demonstrates the potential of using comparative tiling sequencing for genome analysis of crop species. Based on a comparative analysis of the B. rapa sequences and the Arabidopsis genome, it appears that polyploidy and chromosomal diploidization are ongoing processes that collectively stabilize the B. rapa genome and facilitate its evolution.     

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.     

Genetic transformation of Brassica campestris var. rapa protoplasts with an engineered cauliflower mosaic virus genome

Summary A hybrid Cauliflower Mosaic Virus (CaMV) genome containing a selectable marker gene was constructed by replacing the gene VI coding region with the aminoglycoside (neomycin) phosphotransferase type II [APH(3)II] gene from Tn5. This modified viral genome was tested for its infectivity both in planta and in a protoplast transformation system of Brassica campestris var. rapa. Stable, genetically transformed cell lines of B. campestris var. rapa were obtained after transformation. DNA of the hybrid CaMV genome was found to be integrated into high molecular weight plant genomic DNA. Transformation was achieved only when the hybrid genome was supplied together with wild type viral DNA. A possible complementation of the modified CaMV genome with the wild type viral DNA as a helper molecule in planta and in the protoplast system is discussed.     

Toward unraveling membrane biogenesis in mammalian autophagy

Autophagy is a unique form of membrane trafficking, which delivers macromolecules and organelles from the cytoplasm to lysosomes for degradation. This fundamental and ubiquitous process in eukaryotic cells is mediated by the double-membrane-bound structures called autophagosomes, which transiently emerge in the cytoplasm. The recent remarkable explosion of knowledge of autophagy has revealed its multiple roles, especially in mammals; in addition to its basic role in turnover and reuse of cellular constituents, the process unexpectedly functions in elimination of invading bacteria and antigen presentation. Analysis of mammalian homologs of the autophagy-related (Atg) proteins identified in yeast has shed light on not only the common molecular mechanisms underlying autophagosome formation, but also specialized mechanisms that are related to the diverse functions and complex regulation of autophagy in higher organisms.     

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.     

Phylogenetic Analysis of Brassica rapa MATH-Domain Proteins

The MATH (meprin and TRAF-C homology) domain is a fold of seven anti-parallel β-helices involved in protein-protein interaction. Here, we report the identification and characterization of 90 MATH-domain proteins from the Brassica rapa genome. By sequence analysis together with MATH-domain proteins from other species, the B. rapa MATH-domain proteins can be grouped into 6 classes. Class-I protein has one or several MATH domains without any other recognizable domain; Class-II protein contains a MATH domain together with a conserved BTB (Broad Complex, Tramtrack, and Bric-a-Brac ) domain; Class-III protein belongs to the MATH/Filament domain family; Class-IV protein contains a MATH domain frequently combined with some other domains; Class-V protein has a relative long sequence but contains only one MATH domain; Class-VI protein is characterized by the presence of Peptidase and UBQ (Ubiquitinylation) domains together with one MATH domain. As part of our study regarding seed development of B. rapa, six genes are screened by SSH (Suppression Subtractive Hybridization) and their expression levels are analyzed in combination with seed developmental stages, and expression patterns suggested that Bra001786, Bra03578 and Bra036572 may be seed development specific genes, while Bra001787, Bra020541 and Bra040904 may be involved in seed and flower organ development. This study provides the first characterization of the MATH domain proteins in B. rapa     

白菜型油菜黄子资源的初步遗传研究

对随机选取的国内外22份黄子白菜型油菜和22份褐子白菜型油菜进行了种皮色泽的显隐性关系鉴定、黄子性状的等位性测验以及遗传多样性分析。结果表明,黄子白菜型油菜与褐子白菜型油菜配组的杂交组合中,部分组合的F1种皮色泽呈现父本的种皮色泽,表现出花粉直感现象;自然界中存在多种白菜型油菜黄子类型,鉴定出的3种黄子白菜型油菜与褐色白菜型油菜的F2种皮色泽均为褐色,表明黄子性状对于褐子性状为隐性;分子标记方差分析结果显示,白菜型油菜的生长习性所解释的遗传变异大于种皮颜色所解释的遗传变异,表明国外不同生长习性的黄子白菜型油菜资源可用于国内黄子白菜型油菜遗传基础的拓宽。     

扫帚菜-白菜轮作对白菜镉吸收的影响

通过盆栽和大田试验研究扫帚菜Kochia scoparia对土壤中镉(Cd)的富集效率,并利用盆栽试验将其与4个白菜品种进行轮作试验,验证扫帚菜对土壤Cd污染的修复效果。结果表明,扫帚菜各部位富集能力表现为叶>根>茎,富集系数分别为15.07、5.44和2.96;种植扫帚菜后土壤总Cd降低6.02%-13.60%;土壤脲酶和酸性磷酸酶活性也有所提高。盆栽轮作结果表明,扫帚菜-白菜轮作系统中白菜地上部Cd含量与未轮作的对照组相比平均降低17.21%,生物量有略微增加,地上部对Cd的转运系数无明显变化。结果显示,通过扫帚菜与白菜轮作不仅可以增加白菜产量,而且可以有效降低白菜可食部Cd含量,实现边生产边治理的绿色农业理念。     

A sequence-tagged linkage map of Brassica rapa

A detailed genetic linkage map of Brassica rapa has been constructed containing 545 sequence-tagged loci covering 1287 cM, with an average mapping interval of 2.4 cM. The loci were identified using a combination of 520 RFLP and 25 PCR-based markers. RFLP probes were derived from 359 B. rapa EST clones and amplification products of 11 B. rapa and 26 Arabidopsis. Including 21 SSR markers provided anchors to previously published linkage maps for B. rapa and B. napus and is followed as the referenced mapping of R1-R10. The sequence-tagged markers allowed interpretation of the pattern of chromosome duplications within the B. rapa genome and comparison with Arabidopsis. A total of 62 EST markers showing a single RFLP band were mapped through 10 linkage groups, indicating that these can be valuable anchoring markers for chromosome-based genome sequencing of B. rapa. Other RFLP probes gave rise to 2-5 loci, inferring that B. rapa genome duplication is a general phenomenon through 10 chromosomes. The map includes five loci of FLC paralogues, which represent the previously reported BrFLC-1, -2, -3, and -5 and additionally identified BrFLC3 paralogues derived from local segmental duplication on R3.     

芜菁的诱导同源四倍体减数分裂分析

细胞学上芜菁 (Brassicarapa)为正常的二倍体 ,具有 1 1的联会和 1 - 1的染色体分离。它的花粉母细胞 (PMC’s )在小孢子发生过程中的减数分裂前期 (M1 )无“B染色体” ,也无任何次级缔合。用 0 0 2 %的秋水仙碱水溶液处理 6h对多倍体的诱导形成明显。诱导出的同源四倍体构型中无联会基因 ,而成 1或 2次交叉的二价体、三价体和四价体。从多倍体逆转的有趣现象中发现 ,4 6 7%的诱导多倍体后来恢复成二倍体。不平衡的染色体分离普遍存在于诱导出的同源四倍体中 ,这就是导致生存能力不同的异常配子体出现的主要原因。异常配子体的结合 ,在芸苔属 (Brassica)现存种的演化中起到了重要作用