Brassicaceae Species Research Articles

The Brassicaceae genome resource (TBGR): A comprehensive genome platform for Brassicaceae plants.

The Brassicaceae is an important plant family. We built a user-friendly, web-based, comparative, and functional genomic database, The Brassicaceae Genome Resource (TBGR, http://www.tbgr.org.cn), based on 82 released genomes from 27 Brassicaceae species. The TBGR database contains a large number of important functional genes, including 4,096 glucosinolate genes, 6,625 auxin genes, 13,805 flowering genes, 36,632 resistance genes, 1,939 anthocyanin genes, and 1,231 m6A genes. A total of 1,174,049 specific guide sequences for clustered regularly interspaced short palindromic repeats and 5,856,479 transposable elements were detected in Brassicaceae. TBGR also provides information on synteny, duplication, and orthologs for 27 Brassicaceae species. The TBGR database contains 1,183,851 gene annotations obtained using the TrEMBL, Swiss-Prot, Nr, GO, and Pfam databases. The BLAST, Synteny, Primer Design, Seq_fetch, and JBrowse tools are provided to help users perform comparative genomic analyses. All the genome assemblies, gene models, annotations, and bioinformatics results can be easily downloaded from the TBGR database. We plan to improve and continuously update the database with newly assembled genomes and comparative genomic studies. We expect the TBGR database to become a key resource for the study of the Brassicaceae.

Effects of Phospholipase Dε Overexpression on Soybean Response to Nitrogen and Nodulation.

Nitrogen is a key macronutrient to plant growth. We found previously that increased expression of phospholipase Dε (PLDε), which hydrolyzes phospholipids into phosphatidic acid (PA), enhanced plant growth under nitrogen deficiency in Brassicaceae species Arabidopsis and canola. The present study investigated the effect of AtPLDε-overexpression (OE) on soybean (Glycine max), a species capable of symbiotic nitrogen fixation. AtPLDε-OE soybean plants displayed increased root length and leaf size, and the effect of AtPLDε-ΟΕ on leaf size was greater under nitrogen-deficient than -sufficient condition. Under nitrogen deficiency, AtPLDε-OE soybean plants had a higher chlorophyll content and activity of nitrogen assimilation-related enzymes than wild-type soybean plants. AtPLDε-OE led to a higher level of specific PA species in roots after rhizobium inoculation than wild type. AtPLDε-OE soybean plants also increased seed production under nitrogen deprivation with and without nodulation and decreased seed germination in response to high humidity storage and artificial aging. These results suggest that PLDε promotes nitrogen response and affects adversely seed viability during storage.

Divergence in the ABA gene regulatory network underlies differential growth control.

The phytohormone abscisic acid (ABA) is a central regulator of acclimation to environmental stress; however, its contribution to differences in stress tolerance between species is unclear. To establish a comparative framework for understanding how stress hormone signalling pathways diverge across species, we studied the growth response of four Brassicaceae species to ABA treatment and generated transcriptomic and DNA affinity purification and sequencing datasets to construct a cross-species gene regulatory network (GRN) for ABA. Comparison of genes bound directly by ABA-responsive element binding factors suggests that cis-factors are most important for determining the target loci represented in the ABA GRN of a particular species. Using this GRN, we reveal how rewiring of growth hormone subnetworks contributes to stark differences in the response to ABA in the extremophyte Schrenkiella parvula. Our study provides a model for understanding how divergence in gene regulation can lead to species-specific physiological outcomes in response to hormonal cues.

Admixture of divergent genomes facilitates hybridization across species in the family Brassicaceae.

Summary Hybridization and polyploidization are pivotal to plant evolution. Genetic crosses between distantly related species are rare in nature due to reproductive barriers but how such hurdles can be overcome is largely unknown. Here we report the hybrid genome structure of xBrassicoraphanus, a synthetic allotetraploid of Brassica rapa and Raphanus sativus.We performed cytogenetic analysis and de novo genome assembly to examine chromosome behaviors and genome integrity in the hybrid. Transcriptome analysis was conducted to investigate expression of duplicated genes in conjunction with epigenome analysis to address whether genome admixture entails epigenetic reconfiguration.Allotetraploid xBrassicoraphanus retains both parental chromosomes without genome rearrangement. Meiotic synapsis formation and chromosome exchange are avoided between nonhomologous progenitor chromosomes. Reconfiguration of transcription network occurs, and less divergent cis‐elements of duplicated genes are associated with convergent expression. Genome‐wide DNA methylation asymmetry between progenitors is largely maintained but, notably, B. rapa‐originated transposable elements are transcriptionally silenced in xBrassicoraphanus through gain of DNA methylation.Our results demonstrate that hybrid genome stabilization and transcription compatibility necessitate epigenome landscape adjustment and rewiring of cis–trans interactions. Overall, this study suggests that a certain extent of genome divergence facilitates hybridization across species, which may explain the great diversification and expansion of angiosperms during evolution.

Plant products with acetylcholinesterase inhibitory activity for insect control

Acetylcholinesterase (AChE) inhibitors are widely used in Alzheimer’s treatment, but they are also crucial for their action on organophosphorus insecticides. The latter exert their toxicity by inhibiting the AChE enzyme in insects, leading to their death. Amaryllidaceae alkaloids have been proven to be potent AChE inhibitors. In the present study methanolic extracts and essential oils being obtained from species of Asteraceae, Lamiaceae, Brassicaceae and Amaryllidaceae were evaluated in vitro for AChE inhibitory activity. Ellman’s colourimetric method, with modifications, was used for AChE activity evaluation. According to the activity level, the tested plant products were divided into three categories. First: plant products with strong activity comparable to that of galanthamine; second: plant products with medium activity, with IC50 value about 1 mg/ml and the last group with low activity, with IC50 value greater than 1 mg/ml. Essential oils of Origanum vulgare subsp. hirtum Ietswaart., Satureja pilosa Vel., Monarda fistulosa L., Thymus longedentatus (Degen & Urum.) Ronniger and the methanolic extract of Leucojum aestivum L. showed the most potent activity and were referred to as the first group. Carvacrol was identified as the main component of the most active essential oils. In L. aestivum extract, galanthamine was found as the main alkaloid. The obtained results indicate that essential oils and alkaloid-rich plant extracts possess the strongest AChE inhibitory activity. This gives us a reason to recommend these plant products to be tested for insecticidal activity in the future.

Hybridization and introgression are prevalent in Southern European Erysimum (Brassicaceae) species.

Hybridization is a common and important force in plant evolution. One of its outcomes is introgression - the transfer of small genomic regions from one taxon to another by hybridization and repeated backcrossing. This process is believed to be common in glacial refugia, where range expansions and contractions can lead to cycles of sympatry and isolation, creating conditions for extensive hybridization and introgression. Polyploidization is another genome-wide process with a major influence on plant evolution. Both hybridization and polyploidization can have complex effects on plant evolution. However, these effects are often difficult to understand in recently evolved species complexes. We combined flow cytometry, analyses of transcriptomic sequences and pollen tube growth assays to investigate the consequences of polyploidization, hybridization and introgression on the recent evolution of several Erysimum (Brassicaceae) species from the South of the Iberian Peninsula, a well-known glacial refugium. This species complex differentiated in the last 2 million years, and its evolution has been hypothesized to be determined mainly by polyploidization, interspecific hybridization and introgression. Our results support a scenario of widespread hybridization involving both extant and 'ghost' taxa. Several taxa studied here, most notably those with purple corollas, are polyploids, probably of allopolyploid origin. Moreover, hybridization in this group might be an ongoing phenomenon, as pre-zygotic barriers appeared weak in many cases. The evolution of Erysimum spp. has been determined by hybridization to a large extent. Species with purple (polyploids) and yellow flowers (mostly diploid) exhibit a strong signature of introgression in their genomes, indicating that hybridization occurred regardless of colour and across ploidy levels. Although the adaptive value of such genomic exchanges remains unclear, our results demonstrate the significance of hybridization for plant diversification, which should be taken into account when studying plant evolution.

Complete chloroplast genome features of the model heavy metal hyperaccumulator Arabis paniculata Franch and its phylogenetic relationships with other Brassicaceae species.

The online version contains supplementary material available at 10.1007/s12298-022-01151-1.

Novel Brassica hybrids with different resistance to Leptosphaeria maculans reveal unbalanced rDNA signal patterns.

Hybridization of Brassica napus with various Brassicaceae species can result in obtaining new forms with increased resistance to blackleg, a dangerous disease caused mainly by Leptosphaeria maculans. In this study, we aimed to correlate the field resistance of selected Brassica hybrids to blackleg with chromosomal structure revealed by Fluorescence in situ hybridization. Tested genotypes varied in the number of chromosomes and rDNA signals. The greatest variation was observed for A1-type chromosomes. Field evaluation also revealed significant differences in L. maculans resistance. Performed analyses allowed to distinguish three B. napus × Brassica fruticulosa genotypes in which variable patterns of chromosomal structure might be connected to field resistance. However, a more thorough study, including the detection of all A-genome chromosomes, is required.

New insight into the roles of lipid transfer protein and seed storage albumin gene families involved in oil and protein accumulation in rapeseed (Brassica napus L.)

In rapeseed ( Brassica napus L.), lipid transfer protein (LTP) and seed storage albumin (SSA) gene families are involved in lipid and protein metabolism. Understanding the mechanism of oil biosynthesis has great value to increase oil production. The LTPs and SSA gene families were identified of Brassicaceae. The LTPs were divided into six groups according to their sequences, and LTPs of all five species in Brassicaceae also had the same pattern. The genes at a higher expression level were evenly distributed across five groups except group ii. The gene groups in these two gene families showed different expression patterns, with most of the LTP genes expressed at a higher level at 25 days after flowering (DAF), but SSA genes were highly expressed at 40 DAF stage. Gene Ontology (GO) enrichment analysis of the regulatory genes for LTPs also was performed, including response to biotic stimulus (GO:0009607), generation of precursor metabolites and energy (GO:0006091), electron transport chain (GO:0022900), oxidative phosphorylation (GO:0006119), adenosine triphosphate metabolic process (GO:0046034), and phosphorylation (GO:0016310). These results may facilitate further research to understand the expression patterns and regulatory mechanism of LTPs and SSAs.

New records and deletions of Cruciferae for Russia and some neighbouring countries

Based on the study of collections of LE, RV and RWBG, information on occurrence (or absence) of certain Cruciferae (Brassicaceae) species in Russia and some neighbouring states is updated. Camelina hispida is first recorded from Russia, Camelina sativa from Turkmenistan and Goldbachia torulosa from Ukraine as alien species apparently not naturalized in relevant countries. In contrast, Isatis armena, a novelty for the flora of both Russia and Europe, along with Plagioloba persica for the first time recorded from Kyrgyzstan, appears to be aboriginal elements in the newly revealed localities. At the same time, occurrence of Isatis cardiocarpa in Europe, as well as an introduction ofGoldbachia laevigata in Ukraine does not get support. In addition, the combination Plagioloba perfoliata is validated.

Impact of simulated acid rain on the growth of three species Brassica integrifolia, Brassica rapa, Brassica juncea in Hanoi, Vietnam.

Over the past decade, the extent and magnitude of acid rain in Vietnam and other Asian countries have become more apparent. In this study, the effect of simulated acid rain (pH 5.0, 4.0, and 3.0) and control treatment (pH 6.0) are observed for three species Brassica integrifolia, Brassica rapa, and Brassica juncea in Hanoi. The pot experiment was conducted for 42days and arranged according to a randomized complete block design (RCBD), replicated 3 times with acid rain exposure being supplied every 4days. The results show that acid rain causes direct damage to leaves. Observations reveal white spots on leaves; leaves getting discolored and gradually turning yellow, curling leaf marginals, and turning dark blue, with the most severe symptoms being necrotic leaves. Parameters of the shoot and root length, leaf area, biomass, and chlorophyll content all decrease as pH drops. However, the accumulation of proline content in leaves tends to increase with greater acidity. In conclusion, Brassica rara has the highest resistance capability to acid rain compared with Brassica integrifolia and Brassica juncea, especially its proline content is the highest at pH 3.0 in three Brassicaceae species.

Unearthing the plant-microbe quid pro quo in root associations with beneficial fungi.

Mutualistic symbiotic associations between multicellular eukaryotes and their microbiota are driven by the exchange of nutrients in a quid pro quo manner. In the widespread arbuscular mycorrhizal (AM) symbiosis involving plant roots and Glomeromycotina fungi, the mycobiont is supplied with carbon through photosynthesis, which in return supplies the host plant with essential minerals such as phosphorus (P). Most terrestrial plants are largely dependent on AM fungi for nutrients, which raises the question of how plants that are unable to form a functional AM sustain their P nutrition. AM nonhost plants can form alternative, evolutionarily younger, mycorrhizal associations such as the ectomycorrhiza, ericoid and orchid mycorrhiza. However, it is unclear how plants such as the Brassicaceae species Arabidopsis thaliana, which do not form known mycorrhizal symbioses, have adapted to the loss of these essential mycorrhizal traits. Isotope tracing experiments with root-colonizing fungi have revealed the existence of new 'mycorrhizal-like' fungi capable of transferring nutrients such as nitrogen (N) and P to plants, including Brassicaceae. Here, we provide an overview of the biology of trophic relationships between roots and fungi and how these associations might support plant adaptation to climate change.

Genotyping-by-sequencing-based identification of Arabidopsis pattern recognition receptor RLP32 recognizing proteobacterial translation initiation factor IF1

Activation of plant pattern-triggered immunity (PTI) relies on the recognition of microbe-derived structures, termed patterns, through plant-encoded surface-resident pattern recognition receptors (PRRs). We show that proteobacterial translation initiation factor 1 (IF1) triggers PTI in Arabidopsis thaliana and related Brassicaceae species. Unlike for most other immunogenic patterns, IF1 elicitor activity cannot be assigned to a small peptide epitope, suggesting that tertiary fold features are required for IF1 receptor activation. We have deployed natural variation in IF1 sensitivity to identify Arabidopsis leucine-rich repeat (LRR) receptor-like protein 32 (RLP32) as IF1 receptor using a restriction site-associated DNA sequencing approach. RLP32 confers IF1 sensitivity to rlp32 mutants, IF1-insensitive Arabidopsis accessions and IF1-insensitive Nicotiana benthamiana, binds IF1 specifically and forms complexes with LRR receptor kinases SOBIR1 and BAK1 to mediate signaling. Similar to other PRRs, RLP32 confers resistance to Pseudomonas syringae, highlighting an unexpectedly complex array of bacterial pattern sensors within a single plant species.

Comprehensive analysis of 1R- and 2R-MYBs reveals novel genic and protein features, complex organisation, selective expansion and insights into evolutionary tendencies.

Myeloblastosis (MYB) family, the largest plant transcription factor family, has been subcategorised based on the number and type of repeats in the MYB domain. In spite of several reports, evolution of MYB genes and repeats remains enigmatic. Brassicaceae members are endowed with complex genomes, including dysploidy because of its unique history with multiple rounds of polyploidisation, genomic fractionations and rearrangements. The present study is an attempt to gain insights into the complexities of MYB family diversity, understand impacts of genome evolution on gene families and develop an evolutionary framework to understand the origin of various subcategories of MYB gene family. We identified and analysed 1129 MYBs that included 1R-, 2R-, 3R- and atypical-MYBs across sixteen species representing protists, fungi, animals and plants and exclude MYB identified from Brassicaceae except Arabidopsis thaliana; in addition, a total of 1137 2R-MYB genes from six Brassicaceae species were also analysed. Comparative analysis revealed predominance of 1R-MYBs in protists, fungi, animals and lower plants. Phylogenetic reconstruction and analysis of selection pressure suggested ancestral nature of R1-type repeat containing 1R-MYBs that might have undergone intragenic duplication to form multi-repeat MYBs. Distinct differences in gene structure between 1R-MYB and 2R-MYBs were observed regarding intron number, the ratio of gene length to coding DNA sequence (CDS) length and the length of exons encoding the MYB domain. Conserved as well as novel and lineage-specific intron phases were identified. Analyses of physicochemical properties revealed drastic differences indicating functional diversification in MYBs. Phylogenetic reconstruction of 1R- and 2R-MYB genes revealed a shared structure-function relationship in clades which was supported when transcriptome data was analysed in silico. Comparative genomics to study distribution pattern and mapping of 2R-MYBs revealed congruency and greater degree of synteny and collinearity among closely related species. Micro-synteny analysis of genomic segments revealed high conservation of genes that are immediately flanking the surrounding tandemly organised 2R-MYBs along with instances of local duplication, reorganisations and genome fractionation. In summary, polyploidy, dysploidy, reshuffling and genome fractionation were found to cause loss or gain of 2R-MYB genes. The findings need to be supported with functional validation to understand gene structure-function relationship along the evolutionary lineage and adaptive strategies based on comparative functional genomics in plants.

Identification of Potential Genes Encoding Protein Transporters in Arabidopsis thaliana Glucosinolate (GSL) Metabolism.

Several species in Brassicaceae produce glucosinolates (GSLs) to protect themselves against pests. As demonstrated in A. thaliana, the reallocation of defence compounds, of which GSLs are a major part, is highly dependent on transport processes and serves to protect high-value tissues such as reproductive tissues. This study aimed to identify potential GSL-transporter proteins (TPs) using a network-biology approach. The known A. thaliana GSL genes were retrieved from the literature and pathway databases and searched against several co-expression databases to generate a gene network consisting of 1267 nodes and 14,308 edges. In addition, 1151 co-expressed genes were annotated, integrated, and visualised using relevant bioinformatic tools. Based on three criteria, 21 potential GSL genes encoding TPs were selected. The AST68 and ABCG40 potential GSL TPs were chosen for further investigation because their subcellular localisation is similar to that of known GSL TPs (SULTR1;1 and SULTR1;2) and ABCG36, respectively. However, AST68 was selected for a molecular-docking analysis using AutoDOCK Vina and AutoDOCK 4.2 with the generated 3D model, showing that both domains were well superimposed on the homologs. Both molecular-docking tools calculated good binding-energy values between the sulphate ion and Ser419 and Val172, with the formation of hydrogen bonds and van der Waals interactions, respectively, suggesting that AST68 was one of the sulphate transporters involved in GSL biosynthesis. This finding illustrates the ability to use computational analysis on gene co-expression data to screen and characterise plant TPs on a large scale to comprehensively elucidate GSL metabolism in A. thaliana. Most importantly, newly identified potential GSL transporters can serve as molecular tools in improving the nutritional value of crops.

Hybridization and complex evolution of Barbarea vulgaris and related species (Brassicaceae)

Barbarea, winter-cress, is a genus of 29 species in Brassicaceae, the mustard family, which has emerged as a model for evolution of plant defence and specialised metabolites. Notably, some Barbarea species have evolved the ability to produce triterpenoid saponins as the only ones in Brassicaceae, some of which make plants resistant to important herbivores. Resistance has, however, been lost in a distinct group of plants within B. vulgaris ssp. arcuata, which is genetically strongly diverged from other B. vulgaris plants. This divergence is not reflected present in taxonomy. Thus, a phylogeny is needed to understand evolution and defence in Barbarea. Here, we analysed the nuclear ITS and the plastid matK, ndhF, rps16, and psbA-trnH DNA regions from seven out of 29 Barbarea species, 57 accessions of B. vulgaris, 10 accessions of other Barbarea species, and eight outgroup species, in addition to sequences available from GenBank. All Barbarea species formed a highly supported monophyletic group, separated from sister genera. Several clades seem to have radiated within the genus with no simple branching pattern, and discordant nuclear and plastid DNA phylogenies indicate reticulate evolution and chloroplast capture. One of the complex patterns may have resulted from chloroplast capture of a non-Nordic Barbarea species not included in the study. Two pairs of species were almost identical, B. australis and B. grayi, and B. orthoceras and B. stricta. Despite hybridization, chloroplast capture, and incongruence among the plastid and nuclear DNA data, the high level of intraspecific diversity, coupled with lineage specificity, lead us to recognize three groups of Barbarea vulgaris: G-type (glabrous) and P-type (pubescent) individuals of the current B. vulgaris ssp. arcuata as two distinct groups and the current B. vulgaris ssp. vulgaris as the third. Despite the high molecular diversity below species level, the evolutionary history of the saponin-based resistance remains unsettled due to unresolved basal branching.

Niche breadth and elevational range size: a comparative study on Middle-European Brassicaceae species.

High-elevation species are predicted to have larger elevational ranges compared with species of lower elevations. The reasoning is that temperature variability is greater at higher elevation, selecting for wider niche breadth and more plastic genotypes. We used macroevolutionary comparisons involving 90 Brassicaceae species of the central Alps to test for associations among median elevation of occurrence, elevational range size and thermal variability over space and time on the one hand, and their associations with performance breadth or trait plasticity on the other hand. Performance breadth and trait plasticity were estimated by raising replicate plants per species under three temperature treatments (mild, recurrent frost, recurrent heat). Against prediction, we found that mid-elevation species had the largest elevational ranges, and their ranges were associated with increased spatial thermal variability. Nevertheless, variability in the thermal regime was positively associated neither with niche breadth nor with plasticity. Evidence for adaptive constraints was limited to a trade-off between acclimation-based increases in frost and heat resistance, and phylogenetic niche conservatism for median elevation of occurrence and temporal thermal variability. Results suggest that large elevational range size is associated with divergent adaptation within species, but not with more niche breadth or trait plasticity. This article is part of the theme issue 'Species' ranges in the face of changing environments (part I)'.

An LC-MS/MS assay for enzymatic characterization of methylthioalkylmalate synthase (MAMS) involved in glucosinolate biosynthesis.

Brassicaceae are blessed with specialized metabolites called glucosinolates (GSLs), which along with their degradation products, are beneficial in agriculture and human health. To date, more than 130 GSL structures have been identified, mostly derived from the amino acid methionine. The biosynthesis of methionine-derived aliphatic GSLs starts with a side-chain elongation step involving a recursive three-step cyclic process that incorporates a new methylene group into the 2-oxo acid to form a series of elongated 2-oxo acids. Methylthioalkylmalate synthase (MAMS) catalyzes the first committed step in the side-chain elongation of methionine-derived GSLs. The substrate specificity of MAMS with different 2-oxo acids determines whether reaction products of a given cycle enter for an additional round of chain elongation or enter into core GSLs structure formation. Multiple MAMS encoding genes are present in the Brassicaceae species and are known to play a central role in shaping the diverse profile of aliphatic GSLs. We recently established a highly sensitive LC-MS/MS-based methodology that quantifies the MAMS activity by estimating the amount of the next intermediate of the pathway, the 2-malate derivatives. Overall, this chapter describes the protocol for the expression, purification, and steady-state kinetic analysis of the recombinant MAMS protein.

Agro-Morphological Characterization and Some Seed Characteristics of Wild Crambe (Brassicaceae) Species in Turkey

The seeds of wild Crambe species have potential to be used as a source of industrial oil and animal feed. In this study, 48 genotypes of three Crambe species collected from the flora of Turkey were grown under field conditions in Ankara/Turkey in 2014–2016. The seed protein ratio, plant height, number of branches per plant, number of seeds per plant, seed yield per plant, thousand seed weight and hulless/hulled seed ratio (H/H) were determined. The highest protein ratio was determined as 26.02% in the t18 accession of Crambe tataria species. Variations in the characteristics were analyzed using principal component analysis. In the factor analysis of Crambe maritima, Crambe orientalis, Crambe tataria and the mean of these three species, the first two principal components accounted for 100%, 58.06%, 59.93% and 100% of the total variations, respectively. There were positive correlations between the plant height and number of seeds per plant, seed yield per plant for C. orientalis, and number of branches per plant for C. tataria. Although seed yield per plant was high in C. tataria and C. orientalis, they are not suitable for conventional agriculture due to shell thickness, inhomogeneous plant emergence and shooting. Conventional cultivation of wild Crambe species can be made possible by eliminating these negative features with breeding and agronomic studies.

Analyses of Chloroplast Genome of Eutrema japonicum Provide New Insights into the Evolution of Eutrema Species

Wasabi (Eutrema japonicum) is a vegetable of Brassicaceae family, currently cultivated in Southwest Asia. It is rich in nutritional and has a spicy flavour. It is regarded as a rare condiment worldwide. Its genetic profile for yield improvement and the development of E. japonicum germplasm resources remains unknown. Cognizant of this, this study sequenced and assembled the chloroplast (cp) genome of E. japonicum to enrich our genomic information of wasabi and further understand genetic relationships within the Eutrema species. The structural characteristics, phylogeny, and evolutionary relationship of cp genomes among other Brassicaceae plants were analyzed and compared to those of Eutrema species. The cp genome of E. japonicum has 153,851 bp with a typical quadripartite structure, including 37 tRNA genes, 8 rRNA genes, and 87 protein-coding genes. It contains 290 simple sequence repeats and prefers to end their codons with an A or T, which is the same as other Brassicaceae species. Moreover, the cp genomes of the Eutrema species had a high degree of collinearity and conservation during the evolution process. Nucleotide diversity analysis revealed that genes in the IR regions had higher Pi values than those in LSC (Large single copy) and SSC (Small single copy) regions, making them potential molecular markers for wasabi diversity studies. The analysis of genetic distance between Eutrema plants and other Brassicacea plants showed that intraspecies variation was found to be low, while large differences were found between genera and species. Phylogenetic analysis based on 29 cp genomes revealed the existence of a close relationship amongst the Eutrema species. Overall, this study provides baseline information for cp genome-based molecular breeding and genetic transformation studies of Eutrema plants.