Transcriptional Interactions of Single B-Subgenome Chromosome with C-Subgenome in B. oleracea-nigra Additional Lines.

Serial monosomic alien addition lines (MAALs) provide an ideal system to elucidate the transcriptomic interactions between the alien chromosomes and recipient genome under aneuploidy. Herein, five available Brassica oleracea-nigra MAALs (CCB1, CCB4, CCB5, CCB6, CCB8), their derived B. oleracea plants (non-MAALs), and two parents were analyzed for their gene expressions by using high-throughput technology. Compared to parental B. oleracea, all MAALs showed various numbers of DEGs, but CCB8 gave much higher DEGs; the number of downregulated DEGs was slightly higher than the number of upregulated ones, except for in relation to CCB8. All derived B. oleracea plants also gave certain numbers of DEGs, despite these being much lower than in the respective MAALs. Compared to B. nigra, in all five MAALs more DEGs were downregulated than upregulated. Trans-effects were likely more prevailing than cis-effects, and these DEGs were predominantly associated with material transport by dysregulating the cellular component. Meanwhile, the orthologous genes on alien chromosomes could only play a feeble compensatory role for those gene pairs in C-subgenome, and different levels of the expressed genes had a greater tendency towards downregulation. These results revealed transcriptional aneuploidy response patterns between two genomes and suggested that cis- and trans-mechanisms synergistically regulated alien gene transcriptions after distant hybridization.

Fine mapping of genes controlling pigment accumulation in oilseed rape (Brassica napus L.).

Purple/red appearance is one of the common phenotypic variations in leaves, stems, and siliques of oilseed rape (Brassica napus L.) but very rare in flowers. In this study, the causal genes for the purple/red traits in stems and flowers in two accessions of oilseed rape (DH_PR and DH_GC001, respectively) derived from the wide hybridization were fine mapped, and candidate genes were determined by methods combined with bulked segregant analysis (BSA) and RNA-seq analysis. Both traits of purple stem and red flowers were mapped to the locus as AtPAP2 homologous genes (BnaPAP2.C6a and BnaPAP2.A7b, respectively) belonging to the R2R3-MYB family. Sequence comparisons of full-length allelic genes revealed several InDels and SNPs in intron 1 as well as exons, and completely different promoter region of BnaPAP2.C6a and a 211 bp insertion was identified in the promoter region of BnaPAP2.A7b of DH_GC001. Our results not only contribute to a better understanding of anthocyanin inheritance in B. napus, but also provide a useful toolbox for future breeding of cultivars with purple/red traits through the combination of different functional alleles and homologs. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-023-01365-5.

Unraveling Large and Polyploidy Genome of the Crucifer Orychophragmus violaceus in China, a Potential Oil Crop.

The genus Orychophragmus in the Brassicaceae family includes the types with 2n = 20, 22, 24, and 48. The species O. violaceus (L.) O. E. Schulz has 2n = 24 and is widely cultivated as an ornamental plant in China. This review summarizes the research progress of its genome structure and evolution in the context of cytogenetics and genome sequencing. This species has a large genome size of ~1 Gb and longer chromosomes than those of Brassica species, which is attributable to the burst of TE insertions. Even more, one tetraploidization event from about 600-800 million years ago is elucidated to occur during its genome evolution, which is consistent with the polyploidy nature of its genome revealed by the meiotic pairing patterns. Its chromosomes are still characterized by a larger size and deeper staining than those from Brassica species in their intergeneric hybrids, which is likely related to their inherent differences between genome structures and cytology. Its genome is dissected by the development of additional alien lines, and some traits are located on individual chromosomes. Due to the abundant dihydroxy fatty acids in its seed oil with superior lubricant properties and wide environmental adaptations, this plant promises to be utilized as one new oil crop in the future.

Sinapis genomes provide insights into whole-genome triplication and divergence patterns within tribe Brassiceae.

Sinapis alba and Sinapis arvensis are mustard crops within the Brassiceae tribe of the Brassicaceae family, and represent an important genetic resource for crop improvement. We performed the de novo assembly of Brassica nigra, S. alba, and S. arvensis, and conducted comparative genomics to investigate the pattern of genomic evolution since an ancient whole-genome triplication event. Both Sinapis species retained evidence of the Brassiceae whole-genome triplication approximately 20.5 million years ago (Mya), with subgenome dominance observed in gene density, gene expression, and selective constraint. While S. alba diverged from the ancestor of Brassica and Raphanus at approximately 12.5 Mya, the divergence time of S. arvensis and B. nigra was approximately 6.5 Mya. S. arvensis and B. nigra had greater collinearity compared with their relationship to either Brassica rapa or Brassica oleracea. Two chromosomes of S. alba (Sal03 and Sal08) were completely collinear with two ancestral chromosomes proposed in the Ancestral Crucifer Karyotype (ACK) genomic block model, the first time this has been observed in the Brassiceae. These results are consistent with S. alba representing a relatively ancient lineage of the species evolved from the common ancestor of tribe Brassiceae, and suggest that the phylogeny of the Brassica and Sinapis genera requires some revision. Our study provides new insights into the genome evolution and phylogenetic relationships of Brassiceae and provides genomic information for genetic improvement of these plants.

Genome assembly of the Brassicaceae diploid Orychophragmus violaceus reveals complex whole-genome duplication and evolution of dihydroxy fatty acid metabolism.

Orychophragmus violaceus is a Brassicaceae species widely cultivated in China, particularly as a winter cover crop in northern China because of its low-temperature tolerance and low water demand. Recently, O. violaceus has also been cultivated as a potential industrial oilseed crop because of its abundant 24-carbon dihydroxy fatty acids (diOH-FAs), which contribute to superior high-temperature lubricant properties. In this study, we performed de novo assembly of the O. violaceus genome. Whole-genome synteny analysis of the genomes of its relatives demonstrated that O. violaceus is a diploid that has undergone an extra whole-genome duplication (WGD) after the Brassicaceae-specific α-WGD event, with a basic chromosome number of x = 12. Formation of diOH-FAs is hypothesized to have occurred after the WGD event. Based on the genome and the transcriptome data from multiple stages of seed development, we predicted that OvDGAT1-1 and OvDGAT1-2 are candidate genes for the regulation of diOH-FA storage in O. violaceus seeds. These results may greatly facilitate the development of heat-tolerant and eco-friendly plant-based lubricants using O. violaceus seed oil and improve our understanding of the genomic evolution of Brassicaceae.

Transcriptomic and iTRAQ-Based Quantitative Proteomic Analyses of inap CMS in Brassica napus L.

Brassica napus inap cytoplasmic male sterility (CMS) is a novel sterile line with potential application in rapeseed hybrid breeding. Sterile cytoplasm was obtained from Isatis indigotica through somatic fusion and then recurrent backcrossing with B. napus. Previous studies have shown that inap CMS abortion occurred before the stamen primordia (stage 4-5), but the genetic mechanism of sterility needs to be studied. RNA-seq analyses were performed on the floral buds at two stages (0-5 and 6-8), before and after the formation of stamen primordium. As a result, a total of 1769 and 594 differentially expressed genes (DEGs) were detected in the CMS line compared to its maintainer line at the two stages, respectively. In accordance with the CMS phenotype, the up- and downstream regulators of the stamen identity genes AP3 and PI were up- and downregulated in the CMS line, respectively. Furthermore, isobaric tags for relative and absolute quantitation (iTRAQ) analysis showed that a total of 760 differentially abundant proteins (DAPs) were identified in flower buds at stages 0-8, and most of the proteins related to the anther development, oxidative phosphorylation, and programmed cell death (PCD) were downregulated in inap CMS. In combined transcriptomic and proteomic analysis, a total of 32 DEGs/DAPs were identified, of which 7 common DEGs/DAPs had the same expression trend at stage 0-8 of flower development. The downregulation of genes related to the energy deficiency, hormone signal transduction, and the maintenance of mitochondrial metabolic homeostasis at stage 0-5 might disturb the normal differentiation of stamen primordium, resulting in carpelloid stamen of inap CMS. The study will help provide insights into the molecular mechanism of this new male sterility.

Production and cytology of Brassica autoallohexaploids with two and four copies of two subgenomes.

KEY MESSAGE: Different digenomic Brassica autoallohexaploids were produced from the crosses of three allotetraploids and ancestral diploids and characterized for the cytological behavior of two subgenomes with two and four copies. Interspecific hybridization and allopolyploidization present an important pathway for plant evolution and breeding. In this study, different types of digenomic autoallohexaploids with two or four copies of two subgenomes (AAAACC, AACCCC, AAAABB, BBBBCC, BBCCCC) were synthesized by the crosses between three Brassica allotetraploids and their diploid progenitors and the chromosome doubling, and their meiotic behaviors were analyzed by fluorescence in situ hybridization (FISH). These autoallohexaploids showed some variations in pollen fertility and seed-sets and produced both euploid and aneuploid progenies with some chromosomes lost. Two subgenomes in these autoallohexaploids showed some aberrant pairings and segregations, and the degrees of meiotic regularity were negatively associated with the genome affinities. The chromosomes of the subgenome with four copies formed few quadrivalents with the average number < 2, and mainly paired as bivalents, and majority of the chromosomes from the subgenome with two copies gave the expected bivalents. The different extents of the equal and unequal segregations corresponded to the chromosome pairings. The development and cytological investigation of these autoallohexaploids provide not only the new germplasm for genetic research and breeding but also the new clues for the genome behavior and interplay between these subgenomes with different copies.

R gene triplication confers European fodder turnip with improved clubroot resistance.

Clubroot is one of the most important diseases for many important cruciferous vegetables and oilseed crops worldwide. Different clubroot resistance (CR) loci have been identified from only limited species in Brassica, making it difficult to compare and utilize these loci. European fodder turnip ECD04 is considered one of the most valuable resources for CR breeding. To explore the genetic and evolutionary basis of CR in ECD04, we sequenced the genome of ECD04 using de novo assembly and identified 978 candidate R genes. Subsequently, the 28 published CR loci were physically mapped to 15 loci in the ECD04 genome, including 62 candidate CR genes. Among them, two CR genes, CRA3.7.1 and CRA8.2.4, were functionally validated. Phylogenetic analysis revealed that CRA3.7.1 and CRA8.2.4 originated from a common ancestor before the whole-genome triplication (WGT) event. In clubroot susceptible Brassica species, CR-gene homologues were affected by transposable element (TE) insertion, resulting in the loss of CR function. It can be concluded that the current functional CR genes in Brassica rapa and non-functional CR genes in other Brassica species were derived from a common ancestral gene before WGT. Finally, a hypothesis for CR gene evolution is proposed for further discussion.

QTL Identification for Stem Fiber, Strength and Rot Resistance in a DH Population from an Alien Introgression of Brassica napus.

Stem fiber, stem strength and stem-rot resistance are important agronomic traits in Brassica napus. To understand the molecular mechanism that controls the stem-related traits, we investigated the stem lignin (ADL), cellulose (Cel), hemicellulose (Hem) content, S/G monolignol ratio (SG), stem breaking force (BF), breaking strength (F) and Sclerotinia sclerotiorum resistance (SSR). Each trait was significantly positively or negatively correlated with more than three of the other six traits. QTL mapping for ADL, Cel, Hem, SG, BF, F and SSR were performed using a doubled haploid population derived from an intertribal B. napus introgression line 'Y689' crossed with B. napus cv. 'Westar'. A total of 67 additive QTL were identified and integrated into 55 consensus QTL by meta-analysis. Among the 55 consensus QTL, 23 (41.8%) QTL were co-located and were integrated into 11 unique QTL. The QTL by environment (Q × E) interactions were analyzed and 22 combined QTL were identified. In addition, candidate genes within the QTL intervals were proposed based on the known function of Arabidopsis orthologs. These results provided valuable information for improving lodging resistance, S. sclerotiorum resistance and mechanized harvesting of B. napus.

Integrating GWAS, linkage mapping and gene expression analyses reveal the genetic control of first branch height in Brassica napus L.

Rapeseed (Brassica napus L.) is a crucial oil crop cultivated worldwide. First branch height, an essential component of rapeseed plant architecture, has an important effect on yield and mechanized harvesting; however, the underlying genetic mechanism remains unclear. In this study, based on the 60K single nucleotide polymorphism array and a recombinant inbred lines population derived from M083 and 888-5, a total of 19 QTLs were detected in five environments, distributed on linkage groups A02, A09, A10, C06, and C07, which explained phenotypic variation ranging from 4.87 to 29.87%. Furthermore, 26 significant SNPs were discovered on Chr.A02 by genome-wide association study in a diversity panel of 324 re-sequencing accessions. The major QTL of the first branch height trait was co-located on Chr.A02 by integrating linkage mapping and association mapping. Eleven candidate genes were screened via allelic variation analysis, inter-subgenomic synteny analysis, and differential expression of genes in parental shoot apical meristem tissues. Among these genes, BnaA02g13010D, which encodes a TCP transcription factor, was confirmed as the target gene according to gene function annotation, haplotype analysis, and full-length gene sequencing, which revealed that TATA insertion/deletion in the promoter region was closely linked to significantly phenotypic differences BnaA02.TCP1 M083 overexpression resulted in decreased branch height and increased branch number in Arabidopsis. These results provide a genetic basis for first branch height and the ideal architecture of B. napus.

Dosage-Dependent Gynoecium Development and Gene Expression in Brassica napus-Orychophragmus violaceus Addition Lines.

Distant hybridization usually leads to female sterility of the hybrid but the mechanism behind this is poorly understood. Complete pistil abortion but normal male fertility was shown by one Brassica napus-Orychophragmus violaceus monosomic alien addition line (MA, AACC + 1 IO, 2n = 39) produced previously. To study the effect of a single O. violaceus chromosome addition on pistil development in different genetic backgrounds, hybrids between the MA and B. carinata (BBCC), B. juncea (AABB), and two synthetic hexaploids (AABBCC) were firstly produced in this study which show complete female sterility. A microspore culture was further performed to produce the haploid monosomic alien addition line (HMA, AC + 1 IO, 2n = 20) and disomic addition line (DA, AACC + 2 IO, 2n = 40) together with haploid (H, AC, 2n = 19) and double haploid (DH, AACC, 2n = 38) plants of B. napus from MA to investigate the dosage effect of the alien O. violaceus chromosome on pistil development and gene expression. Compared to MA, the development of the pistils of DA and HMA was completely or partially recovered, in which the pistils could swell and elongate to a normal shape after open pollination, although no seeds were produced. Comparative RNA-seq analyses revealed that the numbers of the differentially expressed genes (DEGs) were significantly different, dosage-dependent, and consistent with the phenotypic difference in pairwise comparisons of HMA vs. H, DA vs. DH, MA vs. DH, MA vs. DA, and MA vs. HMA. The gene ontology (GO) enrichment analysis of DEGs showed that a number of genes involved in the development of the gynoecium, embryo sac, ovule, and integuments. Particularly, several common DEGs for pistil development shared in HMA vs. H and DA vs. DH showed functions in genotoxic stress response, auxin transport, and signaling and adaxial/abaxial axis specification. The results provided updated information for the molecular mechanisms behind the gynoecium development of B. napus responding to the dosage of alien O. violaceus chromosomes.

3-Indoleacetonitrile Is Highly Effective in Treating Influenza A Virus Infection In Vitro and In Vivo.

Influenza A viruses are serious zoonotic pathogens that continuously cause pandemics in several animal hosts, including birds, pigs, and humans. Indole derivatives containing an indole core framework have been extensively studied and developed to prevent and/or treat viral infection. This study evaluated the anti-influenza activity of several indole derivatives, including 3-indoleacetonitrile, indole-3-carboxaldehyde, 3-carboxyindole, and gramine, in A549 and MDCK cells. Among these compounds, 3-indoleacetonitrile exerts profound antiviral activity against a broad spectrum of influenza A viruses, as tested in A549 cells. Importantly, in a mouse model, 3-indoleacetonitrile with a non-toxic concentration of 20 mg/kg effectively reduced the mortality and weight loss, diminished lung virus titers, and alleviated lung lesions of mice lethally challenged with A/duck/Hubei/WH18/2015 H5N6 and A/Puerto Rico/8/1934 H1N1 influenza A viruses. The antiviral properties enable the potential use of 3-indoleacetonitrile for the treatment of IAV infection.

Functional homoeologous alleles of CONSTANS contribute to seasonal crop type in rapeseed.

KEY MESSAGE: Two CO paralogs in Brassica napus were confirmed and shown distinct expression pattern and function in promoting flowering and allelic variation s within BnaCO.A10 were found closely associated with ecotype divergence. CONSTANS (CO) is a key gene that responds to photoperiod and in Arabidopsis can promote flowering under long-day (LD) conditions. Brassica napus L. is a major oil crop and close relative of Arabidopsis, and arose via allopolyploidization from the diploids B. rapa (A genome) and B. oleracea (C genome). In this study, we confirmed that B. napus has two CO genes located on the A10 (BnaCO.A10) and C9 (BnaCO.C9) chromosomes. Significant differences in level and temporal pattern of transcription, as well as in protein function, of these homoeologous may have resulted from sequence variation in the promoter as well as in the coding region. Apart from two insertions of 527 bp and 2002 bp in the promoter of BnaCO.C9 that function as transcriptional enhancers, this gene is otherwise highly conserved in both promoter and coding region. However, BnaCO.A10 was classified into two haplotypes and transgene analysis in Arabidopsis and backcross analysis in rapeseed indicated that the winter-type haplotype had a greater effect in promoting flowering than the spring type. We discuss the contribution of CO alleles to species evolution, and for eco-geographic radiation following crop domestication, alongside scope for managing this locus in future breeding.

Genomic signatures of vegetable and oilseed allopolyploid Brassica juncea and genetic loci controlling the accumulation of glucosinolates.

Allopolyploid Brassica juncea crops in Brassicaceae are becoming increasingly revitalized as vegetables and oilseeds owing to wide adaptability and significant economic values. However, the genomic differentiation of diversified vegetables and oilseed B. juncea and the genetic basis underlying glucosinolates accumulation have yet to be elucidated. To address this knowledge gap, we report the sequencing of pairwise genomes of vegetable and oilseed B. juncea at chromosome scale. Comparative genomics analysis unveils panoramic structural variation footprints, particularly the genetic loci of HSP20 and TGA1 associated with abiotic and biotic stresses responses between oilseed and vegetable subgroups. We anchored two major loci of MYB28 (HAG1) orthologues caused by copy number variations on A02 and A09 chromosomes using scored genomic SNPs-based GWAS that are responsible for seed oil quality-determining glucosinolates biosynthesis. These findings will provide valuable repertories of polyploidy genomic information enabling polyploidy genome evolution studies and precise genomic selections for crucial traits like functional components of glucosinolates in B. juncea crops and beyond.

Global gene expression perturbations in rapeseed due to the introduction of alien radish chromosomes.

The aneuploidy chromosome addition lines with individual chromosomes of one species into the whole genome of another species developed from interspecific or intergeneric crosses are ideal for studying the gene expression regulation under asymmetrical genome interactions. Here, we chose two rapeseed-radish addition lines (2n=40, AACC+2R) with all chromosomes of rapeseed and one pair of radish chromosomes carrying 25S rDNA loci and both 5S and 25S rDNA loci, respectively. Transcriptome sequencing of these two rapeseed-radish addition lines, together with two parents (Brassica napus, 2n=38, AACC; Raphanus sativus, 2n=18, RR), was performed to assess gene expression changes due to the aneuploidy effect. Our results showed that the global gene expression perturbations in two addition lines showed asymmetric distributions between A and C subgenomes, for more downregulated genes located in the C subgenome. Moreover, several dysregulated domains occurred in the addition lines and majority of them were clustered in C subgenome, further revealing that C subgenome was more inclined to be repressed by the aneuploidy. Homoeolog expression was slightly biased toward C subgenome in the parent B. napus, but had no preference to either of A or C subgenome in the addition lines. The total number of biased genes increased sharply and most of them were shared only by two addition lines, indicating that aneuploidy could change the extent and direction of homoeolog expression. The triplicated subgenomes which were orthologs to the chromosomes of Arabidopsis thaliana also exhibited dramatic alternations in two additions. Together, our results revealed that the addition of individual radish chromosomes could induce dramatically transcriptomic disturbances and those expression changes gave asymmetric distributions between two subgenomes in B. napus. The possible mechanisms including the uniparental expression of rRNA genes in the allopolyploids and additions for the findings are discussed.

Alternatively Spliced BnaPAP2.A7 Isoforms Play Opposing Roles in Anthocyanin Biosynthesis of Brassica napus L.

Brassica napus L. (rapeseed, oilseed rape, and canola) and varieties of its two diploid parents, B. oleracea and B. rapa, display a large amount of variation in anthocyanin pigmentation of the leaf, stem, and fruit. Here, we demonstrate that BnaPAP2.A7, an ortholog of the B. oleracea anthocyanin activator BoMYB2 that confers purple traits, positively regulates anthocyanin biosynthesis in leaves of B. napus. Sequencing of BnaPAP2.A7 and transgenic analysis suggests that activation of this gene in purple rapeseed may result from a single nucleotide and/or 2bp insertion in its promoter region. BnaPAP2.A7 gives rise to three splice variants, designated BnaPAP2.A7-744, BnaPAP2.A7-910, and BnaPAP2.A7-395 according to the length of the transcripts. While BnaPAP2.A7-744 encodes a full-length R2R3-MYB, both BnaPAP2.A7-910 and BnaPAP2.A7-395 encode truncated proteins that lack both a partial R3 repeat and the complete C terminal domain, and so in vitro are unable to interact with the Arabidopsis bHLH protein AtTT8. Although expression of either BnaPAP2.A7-910 or BnaPAP2.A7-395 in green rapeseed does not result in purple leaves, both genes do modify genome-wide gene expression, with a strong repression of anthocyanin-related genes. We have demonstrated that BnaPAP.A7 regulates anthocyanin accumulation in leaves of B. napus and propose a potential mechanism for modulation of anthocyanin biosynthesis by alternative splicing.

Author Correction: Integration of metabolome and transcriptome reveals flavonoid accumulation in the intergeneric hybrid between Brassica rapa and Raphanus sativus.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Association of Clubroot Resistance Locus PbBa8.1 With a Linkage Drag of High Erucic Acid Content in the Seed of the European Turnip.

Clubroot caused by Plasmodiophora brassicae is a severe threat to the production of Brassica napus, worldwide. The cultivation of resistant varieties is the most efficient and environmentally friendly way to limit disease spread. We developed a highly resistant B. napus line, ZHE226, containing the resistance locus PbBa8.1. However, ZHE226 seeds contain high erucic acid content, which limits its cultivation owing to its low edible oil quality. A segregation population of BC3F2 was developed by crossing ECD04, a resistant European turnip donor, with Huangshuang5, an elite variety with no erucic acid in its seeds, as a recurrent plant. Fine mapping using the bulk segregation analysis sequencing (BSA-Seq) approach detected PbBa8.1 within a 2.9 MB region on chromosome A08. Interestingly, the previously reported resistance gene Crr1a was found in the same region. Genetic analysis revealed that the CAP-134 marker for Crr1a was closely linked with clubroot resistance (CR). Thus, PbBa8.1 and Crr1a might be allelic for CR. Moreover, comparative and genetic analysis showed that high erucic acid in the seeds of ZHE226 was due to linkage drag of fatty acid elongase 1 (FAE1) in the ECD04 line, which was located in the interval of PbBa8.1 with a physical and genetic distance of 729 Kb and 1.86 cm, respectively. Finally, a clubroot-resistant line with a low erucic acid content was successfully developed through gene-specific molecular marker assistant selection from BC4F4. These results will accelerate CR breeding programs in B. napus.

Transposon insertions within alleles of BnaFLC.A10 and BnaFLC.A2 are associated with seasonal crop type in rapeseed.

In Brassicaceae, the requirement for vernalization is conferred by high expression of FLOWERING LOCUS C (FLC). The expression of FLC is known to be repressed by prolonged exposure to cold. Rapeseed (Brassica napus L.) cultivars can be classified into spring, winter, and semi-winter crop types, depending on their respective vernalization requirements. In addition to two known distinct transposon insertion events, here we identified a 4.422 kb hAT and a 5.625 kb long interspersed nuclear element transposon insertion within BnaFLC.A10, and a 810 bp miniature inverted-repeat transposable element (MITE) in BnaFLC.A2. Quantitative PCR demonstrated that these insertions lead to distinct gene expression patterns and contribute differentially to the vernalization response. Transgenic and haplotype analysis indicated that the known 621 bp MITE in the promoter region of BnaFLC.A10 is a transcriptional enhancer that appears to be the main determinant of rapeseed vernalization, and has contributed to the adaptation of rapeseed in winter cultivation environments. In the absence of this transposon insertion, the functional allele of BnaFLC.A2 is a major determinant of vernalization demand. Thus, the combination of BnaFLC.A10 carrying the 621 bp MITE insertion and a functional BnaFLC.A2 appears necessary to establish the winter rapeseed crop phenotype.

A chromosome-scale genome assembly of Isatis indigotica, an important medicinal plant used in traditional Chinese medicine: An Isatis genome.

Isatis indigotica (2n = 14) is an important medicinal plant in China. Its dried leaves and roots (called Isatidis Folium and Isatidis Radix, respectively) are broadly used in traditional Chinese medicine for curing diseases caused by bacteria and viruses such as influenza and viral pneumonia. Various classes of compounds isolated from this species have been identified as effective ingredients. Previous studies based on transcriptomes revealed only a few candidate genes for the biosynthesis of these active compounds in this medicinal plant. Here, we report a high-quality chromosome-scale genome assembly of I. indigotica with a total size of 293.88 Mb and scaffold N50 = 36.16 Mb using single-molecule real-time long reads and high-throughput chromosome conformation capture techniques. We annotated 30,323 high-confidence protein-coding genes. Based on homolog searching and functional annotations, we identified many candidate genes involved in the biosynthesis of main active components such as indoles, terpenoids, and phenylpropanoids. In addition, we found that some key enzyme-coding gene families related to the biosynthesis of these components were expanded due to tandem duplications, which likely drove the production of these major active compounds and explained why I. indigotica has excellent antibacterial and antiviral activities. Our results highlighted the importance of genome sequencing in identifying candidate genes for metabolite synthesis in medicinal plants.