The high‐quality genome of Brassica napus cultivar ‘ZS11’ reveals the introgression history in semi‐winter morphotype

Allotetraploid oilseed rape (Brassica napus L.) is an agriculturally important crop. Cultivation and breeding of B. napus by humans has resulted in numerous genetically diverse morphotypes with optimized agronomic traits and ecophysiological adaptation. To further understand the genetic basis of diversification and adaptation, we report a draft genome of an Asian semi‐winter oilseed rape cultivar ‘ZS11’ and its comprehensive genomic comparison with the genomes of the winter‐type cultivar ‘Darmor‐bzh’ as well as two progenitors. The integrated BAC‐to‐BAC and whole‐genome shotgun sequencing strategies were effective in the assembly of repetitive regions (especially young long terminal repeats) and resulted in a high‐quality genome assembly of B. napus ‘ZS11’. Within a short evolutionary period (~6700 years ago), semi‐winter‐type ‘ZS11’ and the winter‐type ‘Darmor‐bzh’ maintained highly genomic collinearity. Even so, certain genetic differences were also detected in two morphotypes. Relative to ‘Darmor‐bzh’, both two subgenomes of ‘ZS11’ are closely related to its progenitors, and the ‘ZS11’ genome harbored several specific segmental homoeologous exchanges (HEs). Furthermore, the semi‐winter‐type ‘ZS11’ underwent potential genomic introgressions with B. rapa (A_r). Some of these genetic differences were associated with key agronomic traits. A key gene of A03.FLC3 regulating vernalization‐responsive flowering time in ‘ZS11’ was first experienced HE, and then underwent genomic introgression event with A_r, which potentially has led to genetic differences in controlling vernalization in the semi‐winter types. Our observations improved our understanding of the genetic diversity of different B. napus morphotypes and the cultivation history of semi‐winter oilseed rape in Asia.     

Precise editing of CLAVATA genes in Brassica napus L. regulates multilocular silique development

Multilocular silique is a desirable agricultural trait with great potential for the development of high‐yield varieties of Brassica. To date, no spontaneous or induced multilocular mutants have been reported in Brassica napus, which likely reflects its allotetraploid nature and the extremely low probability of the simultaneous random mutagenesis of multiple gene copies with functional redundancy. Here, we present evidence for the efficient knockout of rapeseed homologues of CLAVATA3 (CLV3) for a secreted peptide and its related receptors CLV1 and CLV2 in the CLV signalling pathway using the CRISPR/Cas9 system and achieved stable transmission of the mutations across three generations. Each BnCLV gene has two copies located in two subgenomes. The multilocular phenotype can be recovered only in knockout mutations of both copies of each BnCLV gene, illustrating that the simultaneous alteration of multiple gene copies by CRISPR/Cas9 mutagenesis has great potential in generating agronomically important mutations in rapeseed. The mutagenesis efficiency varied widely from 0% to 48.65% in T₀ with different single‐guide RNAs (sgRNAs), indicating that the appropriate selection of the sgRNA is important for effectively generating indels in rapeseed. The double mutation of BnCLV3 produced more leaves and multilocular siliques with a significantly higher number of seeds per silique and a higher seed weight than the wild‐type and single mutant plants, potentially contributing to increased seed production. We also assessed the efficiency of the horizontal transfer of Cas9/gRNA cassettes by pollination. Our findings reveal the potential for plant breeding strategies to improve yield traits in currently cultivated rapeseed varieties.     

Rapid genome‐wide evolution in Brassica rapa populations following drought revealed by sequencing of ancestral and descendant gene pools

There is increasing evidence that evolution can occur rapidly in response to selection. Recent advances in sequencing suggest the possibility of documenting genetic changes as they occur in populations, thus uncovering the genetic basis of evolution, particularly if samples are available from both before and after selection. Here, we had a unique opportunity to directly assess genetic changes in natural populations following an evolutionary response to a fluctuation in climate. We analysed genome‐wide differences between ancestors and descendants of natural populations of Brassica rapa plants from two locations that rapidly evolved changes in multiple phenotypic traits, including flowering time, following a multiyear late‐season drought in California. These ancestor‐descendant comparisons revealed evolutionary shifts in allele frequencies in many genes. Some genes showing evolutionary shifts have functions related to drought stress and flowering time, consistent with an adaptive response to selection. Loci differentiated between ancestors and descendants (F_ST outliers) were generally different from those showing signatures of selection based on site frequency spectrum analysis (Tajima's D), indicating that the loci that evolved in response to the recent drought and those under historical selection were generally distinct. Very few genes showed similar evolutionary responses between two geographically distinct populations, suggesting independent genetic trajectories of evolution yielding parallel phenotypic changes. The results show that selection can result in rapid genome‐wide evolutionary shifts in allele frequencies in natural populations, and highlight the usefulness of combining resurrection experiments in natural populations with genomics for studying the genetic basis of adaptive evolution.     

Induction of telomere‐mediated chromosomal truncation and behavior of truncated chromosomes in Brassica napus

Engineered minichromosomes could be stably inherited and serve as a platform for simultaneously transferring and stably expressing multiple genes. Chromosomal truncation mediated by repeats of telomeric sequences is a promising approach for the generation of minichromosomes. In the present work, direct repetitive sequences of Arabidopsis telomere were used to study telomere‐mediated truncation of chromosomes in Brassica napus. Transgenes containing alien Arabidopsis telomere were successfully obtained, and Southern blotting and fluorescence in situ hybridization (FISH) results show that the transgenes resulted in successful chromosomal truncation in B. napus. In addition, truncated chromosomes were inherited at rates lower than that predicted by Mendelian rules. To determine the potential manipulations and applications of the engineered chromosomes, such as the stacking of multiple transgenes and the Cre/lox and FRT/FLP recombination systems, both amenable to genetic manipulations through site‐specific recombination in somatic cells, were tested for their ability to undergo recombination in B. napus. These results demonstrate that alien Arabidopsis telomere is able to mediate chromosomal truncation in B. napus. This technology would be feasible for chromosomal engineering and for studies on chromosome structure and function in B. napus.     

Thiol‐based redox proteins in abscisic acid and methyl jasmonate signaling in Brassica napus guard cells

Reversibly oxidized cysteine sulfhydryl groups serve as redox sensors or targets of redox sensing that are important in various physiological processes. However, little is known about redox‐sensitive proteins in guard cells and how they function in stomatal signaling. In this study, Brassica napus guard‐cell proteins altered by redox in response to abscisic acid (ABA) or methyl jasmonate (MeJA) were identified by complementary proteomics approaches, saturation differential in‐gel electrophoresis and isotope‐coded affinity tagging. In total, 65 and 118 potential redox‐responsive proteins were identified in ABA‐ and MeJA‐treated guard cells, respectively. All the proteins contain at least one cysteine, and over half of them are predicted to form intra‐molecular disulfide bonds. Most of the proteins fall into the functional groups of ‘energy’, ‘stress and defense’ and ‘metabolism’. Based on the peptide sequences identified by mass spectrometry, 30 proteins were common to ABA‐ and MeJA‐treated samples. A total of 44 cysteines were mapped in the identified proteins, and their levels of redox sensitivity were quantified. Two of the proteins, a sucrose non‐fermenting 1‐related protein kinase and an isopropylmalate dehydrogenase, were confirmed to be redox‐regulated and involved in stomatal movement. This study creates an inventory of potential redox switches, and highlights a protein redox regulatory mechanism in ABA and MeJA signal transduction in guard cells.     

Maternal components of RNA‐directed DNA methylation are required for seed development in Brassica rapa

Small RNAs trigger repressive DNA methylation at thousands of transposable elements in a process called RNA‐directed DNA methylation (RdDM). The molecular mechanism of RdDM is well characterized in Arabidopsis, yet the biological function remains unclear, as loss of RdDM in Arabidopsis causes no overt defects, even after generations of inbreeding. It is known that 24 nucleotide Pol IV‐dependent siRNAs, the hallmark of RdDM, are abundant in flowers and developing seeds, indicating that RdDM might be important during reproduction. Here we show that, unlike Arabidopsis, mutations in the Pol IV‐dependent small RNA pathway cause severe and specific reproductive defects in Brassica rapa. High rates of abortion occur when seeds have RdDM mutant mothers, but not when they have mutant fathers. Although abortion occurs after fertilization, RdDM function is required in maternal somatic tissue, not in the female gametophyte or the developing zygote, suggesting that siRNAs from the maternal soma might function in filial tissues. We propose that recently outbreeding species such as B. rapa are key to understanding the role of RdDM during plant reproduction.     

Diploidization and genome size change in allopolyploids is associated with differential dynamics of low‐ and high‐copy sequences

Recent advances have highlighted the ubiquity of whole‐genome duplication (polyploidy) in angiosperms, although subsequent genome size change and diploidization (returning to a diploid‐like condition) are poorly understood. An excellent system to assess these processes is provided by Nicotiana section Repandae, which arose via allopolyploidy (approximately 5 million years ago) involving relatives of Nicotiana sylvestris and Nicotiana obtusifolia. Subsequent speciation in Repandae has resulted in allotetraploids with divergent genome sizes, including Nicotiana repanda and Nicotiana nudicaulis studied here, which have an estimated 23.6% genome expansion and 19.2% genome contraction from the early polyploid, respectively. Graph‐based clustering of next‐generation sequence data enabled assessment of the global genome composition of these allotetraploids and their diploid progenitors. Unexpectedly, in both allotetraploids, over 85% of sequence clusters (repetitive DNA families) had a lower abundance than predicted from their diploid relatives; a trend seen particularly in low‐copy repeats. The loss of high‐copy sequences predominantly accounts for the genome downsizing in N. nudicaulis. In contrast, N. repanda shows expansion of clusters already inherited in high copy number (mostly chromovirus‐like Ty3/Gypsy retroelements and some low‐complexity sequences), leading to much of the genome upsizing predicted. We suggest that the differential dynamics of low‐ and high‐copy sequences reveal two genomic processes that occur subsequent to allopolyploidy. The loss of low‐copy sequences, common to both allopolyploids, may reflect genome diploidization, a process that also involves loss of duplicate copies of genes and upstream regulators. In contrast, genome size divergence between allopolyploids is manifested through differential accumulation and/or deletion of high‐copy‐number sequences.     

Plant‐mediated interactions between two herbivores differentially affect a subsequently arriving third herbivore in populations of wild cabbage

• Plants are part of biodiverse communities and frequently suffer from attack by multiple herbivorous insects. Plant responses to these herbivores are specific for insect feeding guilds: aphids and caterpillars induce different plant phenotypes. Moreover, plants respond differentially to single or dual herbivory, which may cascade into a chain of interactions in terms of resistance to other community members. Whether differential responses to single or dual herbivory have consequences for plant resistance to yet a third herbivore is unknown.
• We assessed the effects of single or dual herbivory by Brevicoryne brassicae aphids and/or Plutella xylostella caterpillars on resistance of plants from three natural populations of wild cabbage to feeding by caterpillars of Mamestra brassicae. We measured plant gene expression and phytohormone concentrations to illustrate mechanisms involved in induced responses.
• Performance of both B. brassicae and P. xylostella was reduced when feeding simultaneously with the other herbivore, compared to feeding alone. Gene expression and phytohormone concentrations in plants exposed to dual herbivory were different from those found in plants exposed to herbivory by either insect alone. Plants previously induced by both P. xylostella and B. brassicae negatively affected growth of the subsequently arriving M. brassicae. Furthermore, induced responses varied between wild cabbage populations.
• Feeding by multiple herbivores differentially activates plant defences, which has plant‐mediated negative consequences for a subsequently arriving herbivore. Plant population‐specific responses suggest that plant populations adapt to the specific communities of insect herbivores. Our study contributes to the understanding of plant defence plasticity in response to multiple insect attacks.

     

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.     

Development of diagnostic microsatellite markers from whole‐genome sequences of Ammodramus sparrows for assessing admixture in a hybrid zone

Studies of hybridization and introgression and, in particular, the identification of admixed individuals in natural populations benefit from the use of diagnostic genetic markers that reliably differentiate pure species from each other and their hybrid forms. Such diagnostic markers are often infrequent in the genomes of closely related species, and genomewide data facilitate their discovery. We used whole‐genome data from Illumina HiSeqS2000 sequencing of two recently diverged (600,000 years) and hybridizing, avian, sister species, the Saltmarsh (Ammodramus caudacutus) and Nelson's (A. nelsoni) Sparrow, to develop a suite of diagnostic markers for high‐resolution identification of pure and admixed individuals. We compared the microsatellite repeat regions identified in the genomes of the two species and selected a subset of 37 loci that differed between the species in repeat number. We screened these loci on 12 pure individuals of each species and report on the 34 that successfully amplified. From these, we developed a panel of the 12 most diagnostic loci, which we evaluated on 96 individuals, including individuals from both allopatric populations and sympatric individuals from the hybrid zone. Using simulations, we evaluated the power of the marker panel for accurate assignments of individuals to their appropriate pure species and hybrid genotypic classes (F1, F2, and backcrosses). The markers proved highly informative for species discrimination and had high accuracy for classifying admixed individuals into their genotypic classes. These markers will aid future investigations of introgressive hybridization in this system and aid conservation efforts aimed at monitoring and preserving pure species. Our approach is transferable to other study systems consisting of closely related and incipient species.     

Identification and comprehensive analysis of the CLV3/ESR‐related (CLE) gene family in Brassica napus L.

• The CLE (CLAVATA3/ESR) gene family, encoding a group of small secretory peptides, plays important roles in cell‐to‐cell communication, thereby controlling a broad spectrum of development processes. The CLE family has been systematically characterized in some plants, but not in Brassica napus.
• In the present study, 116 BnCLE genes were identified in the B. napus genome, including seven unannotated, six incorrectly predicted and five multi‐CLE domain‐encoding genes. These BnCLE members were separated into seven distinct groups based on phylogenetic analysis, which might facilitate the functional characterization of the peptides.
• Further characterization of CLE pre‐propeptides revealed 31 unique CLE peptides from 45 BnCLE genes, which may give rise to distinct roles of BnCLE and expansion of the gene family. The biological activity of these unique CLE dodecamer peptides was tested further through in vitro peptide assays. Variations in several important residues were identified as key contributors to the functional differentiation of BnCLE and expansion of the gene family in B. napus. Expression profile analysis helped to characterize possible functional redundancy and sub‐functionalization among the BnCLE members.
• This study presents a comprehensive overview of the CLE gene family in B. napus and provides a foundation for future evolutionary and functional studies.

     

Flip‐flop organization in the chloroplast genome of Capsosiphon fulvescens (Ulvophyceae,Chlorophyta)

To better understand organelle genome evolution of the ulvophycean green alga Capsosiphon fulvescens, we sequenced and characterized its complete chloroplast genome. The circular chloroplast genome was 111,561 bp in length with 31.3% GC content that contained 108 genes including 77 protein‐coding genes, two copies of rRNA operons, and 27 tRNAs. In this analysis, we found the two types of isoform, called heteroplasmy, were likely caused by a flip‐flop organization. The flip‐flop mechanism may have caused structural variation and gene conversion in the chloroplast genome of C. fulvescens. In a phylogenetic analysis based on all available ulvophycean chloroplast genome data, including a new C. fulvescens genome, we found three major conflicting signals for C. fulvescens and its sister taxon Pseudoneochloris marina within 70 individual genes: (i) monophyly with Ulotrichales, (ii) monophyly with Ulvales, and (iii) monophyly with the clade of Ulotrichales and Ulvales. Although the 70‐gene concatenated phylogeny supported monophyly with Ulvales for both species, these complex phylogenetic signals of individual genes need further investigations using a data‐rich approach (i.e., organelle genome data) from broader taxon sampling.     

Knockout of two BnaMAX1 homologs by CRISPR/Cas9‐targeted mutagenesis improves plant architecture and increases yield in rapeseed (Brassica napus L.)

Plant height and branch number are essential components of rapeseed plant architecture and are directly correlated with its yield. Presently, improvement of plant architecture is a major challenge in rapeseed breeding. In this study, we first verified that the two rapeseed BnaMAX1 genes had redundant functions resembling those of Arabidopsis MAX1, which regulates plant height and axillary bud outgrowth. Therefore, we designed two sgRNAs to edit these BnaMAX1 homologs using the CRISPR/Cas9 system. The T₀ plants were edited very efficiently (56.30%–67.38%) at the BnaMAX1 target sites resulting in homozygous, heterozygous, bi‐allelic and chimeric mutations. Transmission tests revealed that the mutations were passed on to the T₁ and T₂ progeny. We also obtained transgene‐free lines created by the CRISPR/Cas9 editing, and no mutations were detected in potential off‐target sites. Notably, simultaneous knockout of all four BnaMAX1 alleles resulted in semi‐dwarf and increased branching phenotypes with more siliques, contributing to increased yield per plant relative to wild type. Therefore, these semi‐dwarf and increased branching characteristics have the potential to help construct a rapeseed ideotype. Significantly, the editing resources obtained in our study provide desirable germplasm for further breeding of high yield in rapeseed.     

Development of iFOX‐hunting as a functional genomic tool and demonstration of its use to identify early senescence‐related genes in the polyploid Brassica napus

Functional genomic studies of many polyploid crops, including rapeseed (Brassica napus), are constrained by limited tool sets. Here we report development of a gain‐of‐function platform, termed ‘iFOX (inducible Full‐length cDNA OvereXpressor gene)‐Hunting’, for inducible expression of B. napus seed cDNAs in Arabidopsis. A Gateway‐compatible plant gene expression vector containing a methoxyfenozide‐inducible constitutive promoter for transgene expression was developed. This vector was used for cloning of random cDNAs from developing B. napus seeds and subsequent Agrobacterium‐mediated transformation of Arabidopsis. The inducible promoter of this vector enabled identification of genes upon induction that are otherwise lethal when constitutively overexpressed and to control developmental timing of transgene expression. Evaluation of a subset of the resulting ~6000 Arabidopsis transformants revealed a high percentage of lines with full‐length B. napus transgene insertions. Upon induction, numerous iFOX lines with visible phenotypes were identified, including one that displayed early leaf senescence. Phenotypic analysis of this line (rsl‐1327) after methoxyfenozide induction indicated high degree of leaf chlorosis. The integrated B. napuscDNA was identified as a homolog of an Arabidopsis acyl‐CoA binding protein (ACBP) gene designated BnACBP1‐like. The early senescence phenotype conferred by BnACBP1‐like was confirmed by constitutive expression of this gene in Arabidopsis and B. napus. Use of the inducible promoter in the iFOX line coupled with RNA‐Seq analyses allowed mechanistic clues and a working model for the phenotype associated with BnACBP1‐like expression. Our results demonstrate the utility of iFOX‐Hunting as a tool for gene discovery and functional characterization of Brassica napus genome.