Brassica taxonomy based on nuclear restriction fragment length polymorphisms (RFLPs)

Summary Preliminary analysis using nuclear RFLPs provided evidence that subspecies within Brassica rapa originated from two different centers. One center is in Europe, represented by turnip and turnip rape from which the oilseed sarson was derived. A second center is in South China containing a variety of Chinese vegetables of which pak choi and narinosa seem to be the most ancient forms. Based on RFLP data, the accessions of B. oleracea examined could be divided into three distinct groups, represented by thousand head kale, broccoli and cabbage. Thousand head kale and Chinese kale appear to be the primitive types. Observations of parallel variation among subspecies of both species are discussed.     

Brassica taxonomy based on nuclear restriction fragment length polymorphisms (RFLPs)

Summary RFLPs were used to study genome evolution and phylogeny in Brassica and related genera. Thirtyeight accessions, including 10 accessions of B. rapa (syn. campestris), 9 cultivated types of B. oleracea, 13 nine-chromosome wild brassicas related to B. oleracea, and 6 other species in Brassica and allied genera, were examined with more then 30 random genomic DNA probes, which identified RFLPs mapping to nine different linkage groups of the B. rapa genome. Based on the RFLP data, phylogenetic trees were constructed using the PAUP microcomputer program. Within B. rapa, accessions of pak choi, narinosa, and Chinese cabbage from East Asia constituted a group distinct from turnip and wild European populations, consistent with the hypothesis that B. rapa had two centers of domestication. A wild B. rapa accession from India was positioned in the tree between European types and East Asian types, suggesting an evolutionary pathway from Europe to India, then to South China. Cultivated B. oleracea morphotypes showed monophyletic origin with wild B. oleracea or B. alboglabra as possible ancestors. Various kales constitute a highly diverse group, and represent the primitive morphotypes of cultivated B. oleracea from which cabbage, broccoli, cauliflower, etc. probably have evolved. Cauliflower was found to be closely related to broccoli, whereas cabbage was closely related to leafy kales. A great diversity existed among the 13 collections of nine-chromosome wild brassicas related to B. oleracea, representing various taxonomic states from subspecies to species. Results from these studies suggested that two basic evolutionary pathways exist for the diploid species examined. One pathway gave rise to B. fruticulosa, B. nigra, and Sinapis arvensis, with B. adpressa or a close relative as the initial ancestor. Another pathway gave rise to B. oleracea and B. rapa, with Diplotaxis erucoides or a close relative as the initial ancestor. Raphanus sativus and Eruca sativus represented intermediate types between the two lineages, and might have been derived from introgression or hybridization between species belonging to different lineages. Molecular evidence for an ascending order of chromosome numbers in the evolution of Brassica and allied genera was obtained on the basis of RFLP data and phylogenetic analysis.     

Isolation and Characterization of an Anther-Specific Polygalacturonase Gene, BcMF16, in Brassica campestris ssp. chinensis

Many genes in the genic male sterile A/B line (Bajh97-01A/B) of Chinese cabbage pak choi (Brassica campestris L. subsp. chinensis Makino) are expressed differentially, and some play critical roles in the formation of pollen walls. In this study, one of these genes, Brassica campestris Male Fertility 16 (BcMF16), has been isolated and characterized. The BcMF16 gene shares approximately 85% nucleotide sequence homology with two exopolygalacturonase (EC3.2.1.67) genes of Arabidopsis thaliana. Cluster analysis of polygalacturonase peptides indicate that BcMF16 belongs to the pollen polygalacturonase clade. Quantitative real-time PCR analysis has revealed that BcMF16 is specifically expressed in reproductive tissues of the fertile line of genic male sterile A/B line of Chinese cabbage pak choi, and that expression levels dramatically increased during later stages of pollen development. In situ hybridization has demonstrated that BcMF16 is specifically and transiently expressed in both tapetum and pollen following microspore separation at the tetrad stage.     

BrpSPL9 (Brassica rapa ssp. pekinensis SPL9) controls the earliness of heading time in Chinese cabbage

The leafy heads of cabbage (Brassica oleracea), Chinese cabbage (Brassica rapa ssp. pekinensis), Brussels sprouts (B. oleracea ssp. gemmifera) and lettuce (Lactuca sativa) comprise extremely incurved leaves that are edible vegetable products. The heading time is important for high quality and yield of these crops. Here, we report that BrpSPL9‐2 (B. rapa ssp. pekinensis SQUAMOSA PROMOTER BINDING‐LIKE 9‐2), a target gene of microRNA brp‐miR156, controls the heading time of Chinese cabbage. Quantitative measurements of leaf shapes, sizes, colour and curvature indicated that heading is a late adult phase of vegetative growth. During the vegetative period, miR156 levels gradually decreased from the seedling stage to the heading one, whereas BrpSPL9‐2 and BrpSPL15‐1 mRNAs increased progressively and reached the highest levels at the heading stage. Overexpression of a mutated miR156‐resistant form of BrpSPL9‐2 caused the significant earliness of heading, concurrent with shortening of the seedling and rosette stages. By contrast, overexpression of miR156 delayed the folding time, concomitant with prolongation of the seedling and rosette stages. Morphological analysis reveals that the significant earliness of heading in the transgenic plants overexpressing BrpSPL9‐2 gene was produced because the juvenile phase was absent and the early adult phase shortened, whereas the significant delay of folding in the transgenic plants overexpressing Brp‐MIR156a was due to prolongation of the juvenile and early adult phases. Thus, miR156 and BrpSPL9 genes are potentially important for genetic improvement of earliness of Chinese cabbage and other crops.     

Beyond genomic variation - comparison and functional annotation of three Brassica rapa genomes: a turnip,a rapid cycling and a Chinese cabbage

Background

Brassica rapa is an economically important crop species. During its long breeding history, a large number of morphotypes have been generated, including leafy vegetables such as Chinese cabbage and pakchoi, turnip tuber crops and oil crops.

Results

To investigate the genetic variation underlying this morphological variation, we re-sequenced, assembled and annotated the genomes of two B. rapa subspecies, turnip crops (turnip) and a rapid cycling. We then analysed the two resulting genomes together with the Chinese cabbage Chiifu reference genome to obtain an impression of the B. rapa pan-genome. The number of genes with protein-coding changes between the three genotypes was lower than that among different accessions of Arabidopsis thaliana, which can be explained by the smaller effective population size of B. rapa due to its domestication. Based on orthology to a number of non-brassica species, we estimated the date of divergence among the three B. rapa morphotypes at approximately 250,000 YA, far predating Brassica domestication (5,000-10,000 YA).

Conclusions

By analysing genes unique to turnip we found evidence for copy number differences in peroxidases, pointing to a role for the phenylpropanoid biosynthesis pathway in the generation of morphological variation. The estimated date of divergence among three B. rapa morphotypes implies that prior to domestication there was already considerably divergence among B. rapa genotypes. Our study thus provides two new B. rapa reference genomes, delivers a set of computer tools to analyse the resulting pan-genome and uses these to shed light on genetic drivers behind the rich morphological variation found in B. rapa.

Electronic supplementary material

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

Variation in abundance of predicted resistance genes in the Brassica oleracea pangenome

Brassica oleracea is an important agricultural species encompassing many vegetable crops including cabbage, cauliflower, broccoli and kale; however, it can be susceptible to a variety of fungal diseases such as clubroot, blackleg, leaf spot and downy mildew. Resistance to these diseases is meditated by specific disease resistance genes analogs (RGAs) which are differently distributed across B. oleracea lines. The sequenced reference cultivar does not contain all B. oleracea genes due to gene presence/absence variation between individuals, which makes it necessary to search for RGA candidates in the B. oleracea pangenome. Here we present a comparative analysis of RGA candidates in the pangenome of B. oleracea. We show that the presence of RGA candidates differs between lines and suggests that in B. oleracea, SNPs and presence/absence variation drive RGA diversity using separate mechanisms. We identified 59 RGA candidates linked to Sclerotinia, clubroot, and Fusarium wilt resistance QTL, and these findings have implications for crop breeding in B. oleracea, which may also be applicable in other crops species.     

Multiple copies of eukaryotic translation initiation factors in Brassica rapa facilitate redundancy,enabling diversification through variation in splicing and broad‐spectrum virus resistance

Recessive strain‐specific resistance to a number of plant viruses in the Potyvirus genus has been found to be based on mutations in the eukaryotic translation initiation factor 4E (eIF4E) and its isoform, eIF(iso)4E. We identified three copies of eIF(iso)4E in a number of Brassica rapa lines. Here we report broad‐spectrum resistance to the potyvirus Turnip mosaic virus (TuMV) due to a natural mechanism based on the mis‐splicing of the eIF(iso)4E allele in some TuMV‐resistant B. rapa var. pekinensis lines. Of the splice variants, the most common results in a stop codon in intron 1 and a much truncated, non‐functional protein. The existence of multiple copies has enabled redundancy in the host plant's translational machinery, resulting in diversification and emergence of the resistance. Deployment of the resistance is complicated by the presence of multiple copies of the gene. Our data suggest that in the B. rapa subspecies trilocularis, TuMV appears to be able to use copies of eIF(iso)4E at two loci. Transformation of different copies of eIF(iso)4E from a resistant B. rapa line into an eIF(iso)4E knockout line of Arabidopsis thaliana proved misleading because it showed that, when expressed ectopically, TuMV could use multiple copies which was not the case in the resistant B. rapa line. The inability of TuMV to access multiple copies of eIF(iso)4E in B. rapa and the broad spectrum of the resistance suggest it may be durable.     

Genome ancestry mosaics reveal multiple and cryptic contributors to cultivated banana

Hybridizations between closely related species commonly occur in the domestication process of many crops. Banana cultivars are derived from such hybridizations between species and subspecies of the Musa genus that have diverged in various tropical Southeast Asian regions and archipelagos. Among the diploid and triploid hybrids generated, those with seedless parthenocarpic fruits were selected by humans and thereafter dispersed through vegetative propagation. Musa acuminata subspecies contribute to most of these cultivars. We analyzed sequence data from 14 M. acuminata wild accessions and 10 M. acuminata‐based cultivars, including diploids and one triploid, to characterize the ancestral origins along their chromosomes. We used multivariate analysis and single nucleotide polymorphism clustering and identified five ancestral groups as contributors to these cultivars. Four of these corresponded to known M. acuminata subspecies. A fifth group, found only in cultivars, was defined based on the ‘Pisang Madu’ cultivar and represented two uncharacterized genetic pools. Diverse ancestral contributions along cultivar chromosomes were found, resulting in mosaics with at least three and up to five ancestries. The commercially important triploid Cavendish banana cultivar had contributions from at least one of the uncharacterized genetic pools and three known M. acuminata subspecies. Our results highlighted that cultivated banana origins are more complex than expected – involving multiple hybridization steps – and also that major wild banana ancestors have yet to be identified. This study revealed the extent to which admixture has framed the evolution and domestication of a crop plant.     

A rich TILLING resource for studying gene function in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Background  

The Brassicaceae family includes the model plant Arabidopsis thaliana as well as a number of agronomically important species such as oilseed crops (in particular Brassica napus, B. juncea and B. rapa) and vegetables (eg. B. rapa and B. oleracea).     

A new avirulence gene of Leptosphaeria maculans,AvrLm14, identifies a resistance source in American broccoli (Brassica oleracea) genotypes

In many cultivated crops, sources of resistance to diseases are sparse and rely on introgression from wild relatives. Agricultural crops often are allopolyploids resulting from interspecific crosses between related species, which are sources of diversity for resistance genes. This is the case for Brassica napus (oilseed rape, canola), an interspecific hybrid between Brassica rapa (turnip) and Brassica oleracea (cabbage). B. napus has a narrow genetic basis and few effective resistance genes against stem canker (blackleg) disease, caused by the fungus Leptosphaeria maculans, are currently available. B. rapa diversity has proven to be a valuable source of resistance (Rlm, LepR) genes, while B. oleracea genotypes were mostly considered susceptible. Here we identified a new resistance source in B. oleracea genotypes from America, potentially effective against French L. maculans isolates under both controlled and field conditions. Genetic analysis of fungal avirulence and subsequent cloning and validation identified a new avirulence gene termed AvrLm14 and suggested a typical gene-for-gene interaction between AvrLm14 and the postulated Rlm14 gene. AvrLm14 shares all the usual characteristics of L. maculans avirulence genes: it is hosted in a genomic region enriched in transposable elements and heterochromatin marks H3K9me3, its expression is repressed during vegetative growth but shows a strong overexpression 5–9 days following cotyledon infection, and it encodes a small secreted protein enriched in cysteine residues with few matches in databases. Similar to the previously cloned AvrLm10-A, AvrLm14 contributes to reduce lesion size on susceptible cotyledons, pointing to a complex interplay between effectors promoting or reducing lesion development.     

High‐throughput multiplex cpDNA resequencing clarifies the genetic diversity and genetic relationships among Brassica napus,Brassica rapa and Brassica oleracea

Brassica napus (rapeseed) is a recent allotetraploid plant and the second most important oilseed crop worldwide. The origin of B. napus and the genetic relationships with its diploid ancestor species remain largely unresolved. Here, chloroplast DNA (cpDNA) from 488 B. napus accessions of global origin, 139 B. rapa accessions and 49 B. oleracea accessions were populationally resequenced using Illumina Solexa sequencing technologies. The intraspecific cpDNA variants and their allelic frequencies were called genomewide and further validated via EcoTILLING analyses of the rpo region. The cpDNA of the current global B. napus population comprises more than 400 variants (SNPs and short InDels) and maintains one predominant haplotype (Bncp1). Whole‐genome resequencing of the cpDNA of Bncp1 haplotype eliminated its direct inheritance from any accession of the B. rapa or B. oleracea species. The distribution of the polymorphism information content (PIC) values for each variant demonstrated that B. napus has much lower cpDNA diversity than B. rapa; however, a vast majority of the wild and cultivated B. oleracea specimens appeared to share one same distinct cpDNA haplotype, in contrast to its wild C‐genome relatives. This finding suggests that the cpDNA of the three Brassica species is well differentiated. The predominant B. napus cpDNA haplotype may have originated from uninvestigated relatives or from interactions between cpDNA mutations and natural/artificial selection during speciation and evolution. These exhaustive data on variation in cpDNA would provide fundamental data for research on cpDNA and chloroplasts.     

Diamondback moth egg susceptibility to rainfall: effects of host plant and oviposition behavior

Oviposition patterns of the diamondback moth (DBM), Plutella xylostella L. (Lepidoptera: Plutellidae), differ between common cabbage (Brassica oleracea L. var. capitata) and Chinese cabbage (Brassica rapa L. var. pekinensis) (Brassicaceae) host plants. This study shows that the moth prefers to oviposit on adaxial rather than abaxial leaf surfaces and petioles of both host plants. More eggs were laid in leaf veins than on leaf laminas of both host plants, especially in Chinese cabbage, where 94.6% of eggs were laid in veins. On Chinese cabbage, very few eggs were laid in clusters (≥2 eggs), whereas on common cabbage approximately 30% of eggs were laid in groups of 2 or more eggs. Removal of wax from common cabbage leaves dramatically increased the number of eggs laid singly on the leaf lamina of treated plants, suggesting that leaf waxes affect how eggs are distributed by ovipositing DBM. Eggs were most susceptible to removal by rainfall from the plant surface immediately (<1 h) after oviposition and when close to hatching (>72h old) whereas they were least susceptible 24 h after oviposition. Eggs laid on common cabbage plants were more susceptible to simulated rainfall than eggs laid on Chinese cabbage plants. On common cabbage plants, egg susceptibility to rainfall on different plant parts ranked adaxial leaf surfaces>petioles = abaxial leaf surfaces>stem, but there was no difference in egg susceptibility to rainfall on the various plant parts of Chinese cabbage. Furthermore, on common cabbage plants, eggs laid on both adaxial and abaxial leaf surfaces were afforded significant protection from the effects of rainfall by leaves higher in the plant canopy. On common cabbage plants, oviposition patterns reduce the potential impact of rainfall on eggs, possibly reducing the effect of this important abiotic mortality factor in the field.     

Gene transferability from transgenic Brassica napus L. to various subspecies and varieties of Brassica rapa

Gene transferability from transgenic rapeseed to various subspecies and varieties of Brassica rapa was assessed in this study. Artificial crossability was studied in 118 cultivars of 7 B. rapa subspecies and varieties with the transgenic rapeseed GT73 (Brassica napus) as the pollen donor. On average 5.7 seeds were obtained per pollination, with a range from 0.05 to 19.4. The heading type of B. rapa L. showed significantly higher crossability than non-heading types of B. rapa. The spontaneous outcrossing rate between B. rapa (female) and the transgenic rapeseed Ms8 × Rf3 (B. napus) (male) ranged from 0.039 to 0.406%, with an average of 0.19%. The fertilization process and the development of the hybrid seeds as shown by fluorescent staining techniques indicated that the number of adhered pollens on the stigma was reduced by 80%, the number of pollen tubes in the style was reduced by 2/3 and the fertilization time was delayed by over 20 h when pollinated with the transgenic rapeseed Ms8 × Rf3 in comparison with the bud self-pollination of B. rapa as control. About 10–70% of the interspecific hybrid embryos were aborted in the course of development. Some seeds looked cracked in mature pods, which showed germination abilities lower than 10%. The spontaneous outcrossing rates were much lower than the artificial crossability, and their survival fitness of the interspecific hybrid was very low, indicating that it should be possible to keep the adventitious presence of the off-plants under the allowed threshold, if proper measures are taken.     

Diamondback moth oviposition: effects of host plant and herbivory

The oviposition behaviour of Plutella xylostella L. (Lepidoptera: Plutellidae) on Chinese cabbage (Brassica rapa L. Pekinensis, cv. Wombok), canola (Brassica napus L. cv. Thunder TT), and cabbage (Brassica oleracea L. Capitata, cv. sugarloaf) (Brassicaceae) was studied in the laboratory. In no‐choice experiments moths laid most eggs on the stems and lower three leaves of cabbage plants, the lower three leaves of canola plants, but on the upper three leaves of Chinese cabbage plants. The effects of conspecific herbivore damage to foliage could be replicated by mechanical damage. When foliage was damaged, injured cabbage and canola plants were preferred for oviposition over intact conspecifics, whereas injured Chinese cabbage plants were less preferred than intact conspecifics. However, when root tissue was damaged, intact cabbage and canola plants were preferred over injured conspecifics, whereas moths did not discriminate between root‐damaged and intact Chinese cabbage plants. Injury to upper leaves significantly affected the intra‐plant distribution of eggs. In cabbage and canola plants, injury to leaf 6 significantly increased the number of eggs laid on this leaf, resulting in a significant decrease in the number of eggs laid on the lower foliage/stem of plants, whereas in Chinese cabbage plants it significantly decreased the number of eggs laid on leaf 6. Following oviposition on intact plants, neonate larvae established the vast majority of feeding sites on leaves 5–8 in all three host plants, indicating that larvae moved a considerable distance from preferred oviposition sites in cabbage and canola plants. The growth rate of neonates fed on leaf‐6 tissue was significantly greater than that of those fed on leaf‐1 tissue; >90% of larvae completed development when fed exclusively on leaf‐6 tissue but no larvae completed development when fed exclusively on leaf‐1 tissue. The study demonstrates the complex and unpredictable interactions between P. xylostella and its host plants and provides a basis from which we can begin to understand observed distributions of the pest in Brassica crops.     

Development of public immortal mapping populations,molecular markers and linkage maps for rapid cycling <Emphasis Type="Italic">Brassica rapa</Emphasis> and <Emphasis Type="Italic">B. oleracea</Emphasis>

Publicly available genomic tools help researchers integrate information and make new discoveries. In this paper, we describe the development of immortal mapping populations of rapid cycling, self-compatible lines, molecular markers, and linkage maps for Brassica rapa and B. oleracea and make the data and germplasm available to the Brassica research community. The B. rapa population consists of 160 recombinant inbred (RI) lines derived from the cross of highly inbred lines of rapid cycling and yellow sarson B. rapa. The B. oleracea population consists of 155 double haploid (DH) lines derived from an F1 cross between two DH lines, rapid cycling and broccoli. A total of 120 RFLP probes, 146 SSR markers, and one phenotypic trait (flower color) were used to construct genetic linkage maps for both species. The B. rapa map consists of 224 molecular markers distributed along 10 linkage groups (A1–A10) with a total distance of 1125.3 cM and a marker density of 5.7 cM/marker. The B. oleracea genetic map consists of 279 molecular markers and one phenotypic marker distributed along nine linkage groups (C1–C9) with a total distance of 891.4 cM and a marker density of 3.2 cM/marker. A syntenic analysis with Arabidopsis thaliana identified collinear genomic blocks that are in agreement with previous studies, reinforcing the idea of conserved chromosomal regions across the Brassicaceae.     

Genetic changes in a novel breeding population of Brassica napus synthesized from hundreds of crosses between B. rapa and B. carinata

Introgression of genomic variation between and within related crop species is a significant evolutionary approach for population differentiation, genome reorganization and trait improvement. Using the Illumina Infinium Brassica 60K SNP array, we investigated genomic changes in a panel of advanced generation new‐type Brassica napus breeding lines developed from hundreds of interspecific crosses between 122 Brassica rapa and 74 Brassica carinata accessions, and compared them with representative accessions of their three parental species. The new‐type B. napus population presented rich genetic diversity and abundant novel genomic alterations, consisting of introgressions from B. rapa and B. carinata, novel allelic combinations, reconstructed linkage disequilibrium patterns and haplotype blocks, and frequent deletions and duplications (nonrandomly distributed), particularly in the C subgenome. After a much shorter, but very intensive, selection history compared to traditional B. napus, a total of 15 genomic regions with strong selective sweeps and 112 genomic regions with putative signals of selective sweeps were identified. Some of these regions were associated with important agronomic traits that were selected for during the breeding process, while others were potentially associated with restoration of genome stability and fertility after interspecific hybridization. Our results demonstrate how a novel method for population‐based crop genetic improvement can lead to rapid adaptation, restoration of genome stability and positive responses to artificial selection.     

Brassica rapa Domestication: Untangling Wild and Feral Forms and Convergence of Crop Morphotypes

The study of domestication contributes to our knowledge of evolution and crop genetic resources. Human selection has shaped wild Brassica rapa into diverse turnip, leafy, and oilseed crops. Despite its worldwide economic importance and potential as a model for understanding diversification under domestication, insights into the number of domestication events and initial crop(s) domesticated in B. rapa have been limited due to a lack of clarity about the wild or feral status of conspecific noncrop relatives. To address this gap and reconstruct the domestication history of B. rapa, we analyzed 68,468 genotyping-by-sequencing-derived single nucleotide polymorphisms for 416 samples in the largest diversity panel of domesticated and weedy B. rapa to date. To further understand the center of origin, we modeled the potential range of wild B. rapa during the mid-Holocene. Our analyses of genetic diversity across B. rapa morphotypes suggest that noncrop samples from the Caucasus, Siberia, and Italy may be truly wild, whereas those occurring in the Americas and much of Europe are feral. Clustering, tree-based analyses, and parameterized demographic inference further indicate that turnips were likely the first crop type domesticated, from which leafy types in East Asia and Europe were selected from distinct lineages. These findings clarify the domestication history and nature of wild crop genetic resources for B. rapa, which provides the first step toward investigating cases of possible parallel selection, the domestication and feralization syndrome, and novel germplasm for Brassica crop improvement.     

Risk assessment of Bacillus thuringiensis in wild Brassica rapa: A field simulation of introgressed transgenes

Abstract

Brassica crops are able to hybridize with closely related wild and weedy species such as Brassica rapa and therefore these species could become recipients of transgenes from GM brassica crops, such as oilseed rape. Transgenes which protect against herbivory, such as a gene conferring the production of a Bt toxin, could increase the recipient's fitness and therefore enhance its competitiveness. We used microbial Bt and several other pesticides to exclude several guilds of herbivores and thus simulate the effect of transgene introgression on the performance of wild B. rapa. There were only minor negative effects of excluding insect herbivores on the performance of B. rapa and it appears that vertebrate herbivory has a more significant effect. The advantages and disadvantages of simulating the risks associated with GM plants are discussed. If we apply this simulation to GM risk assessment and make inferences from these findings, the implications of this for gene flow from insect-resistant transgenic plants are that it will be of little significance should it occur.