Characterization of the complete chloroplast genome of Brassica oleracea var. italica and phylogenetic relationships in Brassicaceae.

Broccoli (Brassica oleracea var. italica) is an important B. oleracea cultivar, with high economic and agronomic value. However, comparative genome analyses are still needed to clarify variation among cultivars and phylogenetic relationships within the family Brassicaceae. Herein, the complete chloroplast (cp) genome of broccoli was generated by Illumina sequencing platform to provide basic information for genetic studies and to establish phylogenetic relationships within Brassicaceae. The whole genome was 153,364 bp, including two inverted repeat (IR) regions of 26,197 bp each, separated by a small single copy (SSC) region of 17,834 bp and a large single copy (LSC) region of 83,136 bp. The total GC content of the entire chloroplast genome accounts for 36%, while the GC content in each region of SSC,LSC, and IR accounts for 29.1%, 34.15% and 42.35%, respectively. The genome harbored 133 genes, including 88 protein-coding genes, 37 tRNAs, and 8 rRNAs, with 17 duplicates in IRs. The most abundant amino acid was leucine and the least abundant was cysteine. Codon usage analyses revealed a bias for A/T-ending codons. A total of 35 repeat sequences and 92 simple sequence repeats were detected, and the SC-IR boundary regions were variable between the seven cp genomes. A phylogenetic analysis suggested that broccoli is closely related to Brassica oleracea var. italica MH388764.1, Brassica oleracea var. italica MH388765.1, and Brassica oleracea NC_0441167.1. Our results are expected to be useful for further species identification, population genetics analyses, and biological research on broccoli.

BRAD V3.0: an upgraded Brassicaceae database.

The Brassicaceae Database (BRAD version 3.0, BRAD V3.0; http://brassicadb.cn) has evolved from the former Brassica Database (BRAD V2.0), and represents an important community portal hosting genome information for multiple Brassica and related Brassicaceae plant species. Since the last update in 2015, the complex genomes of numerous Brassicaceae species have been decoded, accompanied by many omics datasets. To provide an up-to-date service, we report here a major upgrade of the portal. The Model-View-ViewModel (MVVM) framework of BRAD has been re-engineered to enable easy and sustainable maintenance of the database. The collection of genomes has been increased to 26 species, along with optimization of the user interface. Features of the previous version have been retained, with additional new tools for exploring syntenic genes, gene expression and variation data. In the 'Syntenic Gene @ Subgenome' module, we added features to view the sequence alignment and phylogenetic relationships of syntenic genes. New modules include 'MicroSynteny' for viewing synteny of selected fragment pairs, and 'Polymorph' for retrieval of variation data. The updated BRAD provides a substantial expansion of genomic data and a comprehensive improvement of the service available to the Brassicaceae research community.

Genome evolution of the psammophyte Pugionium for desert adaptation and further speciation.

Deserts exert strong selection pressures on plants, but the underlying genomic drivers of ecological adaptation and subsequent speciation remain largely unknown. Here, we generated de novo genome assemblies and conducted population genomic analyses of the psammophytic genus Pugionium (Brassicaceae). Our results indicated that this bispecific genus had undergone an allopolyploid event, and the two parental genomes were derived from two ancestral lineages with different chromosome numbers and structures. The postpolyploid expansion of gene families related to abiotic stress responses and lignin biosynthesis facilitated environmental adaptations of the genus to desert habitats. Population genomic analyses of both species further revealed their recent divergence with continuous gene flow, and the most divergent regions were found to be centered on three highly structurally reshuffled chromosomes. Genes under selection in these regions, which were mainly located in one of the two subgenomes, contributed greatly to the interspecific divergence in microhabitat adaptation.

Assessment of antioxidant and antimicrobial activities of Crambe spp. during vegetation.

BACKGROUND: The search for new plant raw material as a potential source of antioxidants is still ongoing. This study aimed to evaluate the antioxidant and antimicrobial capacity of the plant raw material of Crambe spp. during vegetation. METHODS: The free radical scavenging activity and molybdenum reducing power of the extracts were used to determine antioxidant activity. The quantification of polyphenol compounds was conducted with Folin-Ciocalteu reagent. Flavonoids and phenolic acids were also determined. The disc diffusion method was used to determine antimicrobial activity. RESULTS: It was determined that the free radical scavenging activity, assessed using the DPPH-method, was 4.38-8.20 mg TE/g DW, the molybdenum reducing power of the extracts was 40.07-129.12 mg TE/g DW, total polyphenol content was 20.24-70.88 mg GAE/g DW, total flavonoid content was 5.73-29.92 mg QE/g DW, and phenolic acid content was 3.00-10.63 mg CAE/g DW. Antimicrobial activity depended on the stage of growth and the part of the plant used. CONCLUSIONS: Crambe spp. possess the antioxidant and antimicrobial potential to mean that they could be used in pharmaceutical studies and the food industry.

Complete chloroplast genome features and phylogenetic analysis of Eruca sativa (Brassicaceae).

Eruca sativa Mill. (Brassicaceae) is an important edible vegetable and a potential medicinal plant due to the antibacterial activity of its seed oil. Here, the complete chloroplast (cp) genome of E. sativa was de novo assembled with a combination of long PacBio reads and short Illumina reads. The E. sativa cp genome had a quadripartite structure that was 153,522 bp in size, consisting of one large single-copy region of 83,320 bp and one small single-copy region of 17,786 bp which were separated by two inverted repeat (IRa and IRb) regions of 26,208 bp. This complete cp genome harbored 113 unique genes: 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Forty-nine long repetitive sequences and 69 simple sequence repeats were identified in the E. sativa cp genome. A codon usage analysis of the E. sativa cp genome showed a bias toward codons ending in A/T. The E. sativa cp genome was similar in size, gene composition, and linearity of the structural region when compared with other Brassicaceae cp genomes. Moreover, the analysis of the synonymous (Ks) and non-synonymous (Ka) substitution rates demonstrated that protein-coding genes generally underwent purifying selection pressure, expect ycf1, ycf2, and rps12. A phylogenetic analysis determined that E. sativa is evolutionarily close to important Brassica species, indicating that it may be possible to transfer favorable E. sativa alleles into other Brassica species. Our results will be helpful to advance genetic improvement and breeding of E. sativa, and will provide valuable information for utilizing E. sativa as an important resource to improve other Brassica species.

Accumulation of potentially toxic elements by plants of North Caucasian Alyssum species and their molecular phylogenetic analysis.

Accumulations of potentially toxic metals were investigated in soils and five North Caucasian Alyssum species from metalliferous areas and non-metalliferous areas in Karachay-Cherkessia, Kabardino-Balkaria, Dagestan and the Krasnodar region. Analyses of field samples showed that chemical features of the soils significantly affected the concentrations of Ni, Co, Zn, but had less effect on Cu and Pb concentrations in the shoots of Alyssum. Variations in the degree of accumulating ability were found in the studied species, including hyperaccumulation of Ni in Alyssum murale (up to 12,100 mg kg-1), and significant accumulation of Zn in A. gehamense (up to 1700 mg kg-1). A comparative molecular genetic analysis of two A. murale populations, both Ni-hyperaccumulating population from Karachay-Cherkessia and non-hyperaccumulating population from Dagestan, indicated considerable genetic difference between them. This result supports the hypothesis that the selection of metal hyperaccumulator species with enhanced phytoremediation efficiency should be considered at the population level.

Brassicaceae transcriptomes reveal convergent evolution of super-accumulation of sinigrin.

Wasabi, horseradish and mustard are popular pungent crops in which the characteristic bioactive hydrolysis of specialized glucosinolates (GSLs) occurs. Although the metabolic pathways of GSLs are well elucidated, how plants have evolved convergent mechanisms to accumulate identical GSL components remains largely unknown. In this study, we discovered that sinigrin is predominantly synthesized in wasabi, horseradish and mustard in Brassicaceae. We de novo assembled the transcriptomes of the three species, revealing the expression patterns of gene clusters associated with chain elongation, side chain modification and transport. Our analysis further revealed that several gene clusters were convergently selected during evolution, exhibiting convergent shifts in amino acid preferences in mustard, wasabi and horseradish. Collectively, our findings provide insights into how unrelated crop species evolve the capacity for sinigrin super-accumulation and thus promise a potent strategy for engineering metabolic pathways at multiple checkpoints to fortify bioactive compounds for condiment or pharmaceutical purposes.

Genome Size Evolution Mediated by Gypsy Retrotransposons in Brassicaceae.

The dynamic activity of transposable elements (TEs) contributes to the vast diversity of genome size and architecture among plants. Here, we examined the genomic distribution and transposition activity of long terminal repeat retrotransposons (LTR-RTs) in Arabidopsis thaliana (Ath) and three of its relatives, Arabidopsis lyrata (Aly), Eutrema salsugineum (Esa), and Schrenkiella parvula (Spa), in Brassicaceae. Our analyses revealed the distinct evolutionary dynamics of Gypsyretrotransposons, which reflects the different patterns of genome size changes of the four species over the past million years. The rate of Gypsy transposition in Aly is approximately five times more rapid than that of Ath and Esa, suggesting an expanding Aly genome. Gypsy insertions in Esa are strictly confined to pericentromeric heterochromatin and associated with dramatic centromere expansion. In contrast, Gypsy insertions in Spa have been largely suppressed over the last million years, likely as a result of a combination of an inherent molecular mechanism of preferential DNA removal and purifying selection at Gypsy elements. Additionally, species-specific clades of Gypsy elements shaped the distinct genome architectures of Aly and Esa.

An updated explanation of ancestral karyotype changes and reconstruction of evolutionary trajectories to form Camelina sativa chromosomes.

BACKGROUND: Belonging to lineage I of Brassicaceae, Camelina sativa is formed by two hybridizations of three species (three sub-genomes). The three sub-genomes were diverged from a common ancestor, likely derived from lineage I (Ancestral Crucifer karyotype, ACK). The karyotype evolutionary trajectories of the C. sativa chromosomes are currently unknown. Here, we managed to adopt a telomere-centric theory proposed previously to explain the karyotype evolution in C. sativa. RESULTS: By characterizing the homology between A. lyrata and C. sativa chromosomes, we inferred ancestral diploid karyotype of C. sativa (ADK), including 7 ancestral chromosomes, and reconstructed the evolutionary trajectories leading to the formation of extant C. sativa genome. The process involved 2 chromosome fusions. We found that sub-genomes Cs-G1 and Cs-G2 may share a closer common ancestor than Cs-G3. Together with other lines of evidence from Arabidopsis, we propose that the Brassicaceae plants, even the eudicots, follow a chromosome fusion mechanism favoring end-end joining of different chromosomes, rather than a mechanism favoring the formation circular chromosomes and nested chromosome fusion preferred by the monocots. CONCLUSIONS: The present work will contribute to understanding the formation of C. sativa chromosomes, providing insight into Brassicaceae karyotype evolution.

Phylogeny of Euclidieae (Brassicaceae) based on plastome and nuclear ribosomal DNA data.

Euclidieae, a morphologically diverse tribe in the family Brassicaceae (Cruciferae), consists of 29 genera and more than 150 species distributed mainly in Asia. Prior phylogenetic analyses on Euclidieae are inadequate. In this study, sequence data from the plastid genome and nuclear ribosomal DNA of 72 species in 27 genera of Euclidieae were used to infer the inter- and intra-generic relationships within. The well-resolved and strongly supported plastome phylogenies revealed that Euclidieae could be divided into five clades. Both Cymatocarpus and Neotorularia are polyphyletic in nuclear and plastome phylogenies. Besides, the conflicts of systematic positions of three species of Braya and two species of Solms-laubachia s.l. indicated that hybridization and or introgression might have happened during the evolutionary history of the tribe. Results from divergence-time analyses suggested an early Miocene origin of Euclidieae, and it probably originated from the Central Asia, Pamir Plateau and West Himalaya. In addition, multiple ndh genes loss and pseudogenization were detected in eight species based on comparative genomic study.

Nested whole-genome duplications coincide with diversification and high morphological disparity in Brassicaceae.

Angiosperms have become the dominant terrestrial plant group by diversifying for ~145 million years into a broad range of environments. During the course of evolution, numerous morphological innovations arose, often preceded by whole genome duplications (WGD). The mustard family (Brassicaceae), a successful angiosperm clade with ~4000 species, has been diversifying into many evolutionary lineages for more than 30 million years. Here we develop a species inventory, analyze morphological variation, and present a maternal, plastome-based genus-level phylogeny. We show that increased morphological disparity, despite an apparent absence of clade-specific morphological innovations, is found in tribes with WGDs or diversification rate shifts. Both are important processes in Brassicaceae, resulting in an overall high net diversification rate. Character states show frequent and independent gain and loss, and form varying combinations. Therefore, Brassicaceae pave the way to concepts of phylogenetic genome-wide association studies to analyze the evolution of morphological form and function.

Life-History Plasticity and Water-Use Trade-Offs Associated with Drought Resistance in a Clade of California Jewelflowers.

Water limitation is a primary driver of plant geographic distributions and individual plant fitness. Drought resistance is the ability to survive and reproduce despite limited water, and numerous studies have explored its physiological basis in plants. However, it is unclear how drought resistance and trade-offs associated with drought resistance evolve within plant clades. We quantified the relationship between water availability and fitness for 13 short-lived plant taxa in the Streptanthus clade that vary in their phenology and the availability of water in the environments where they occur. We derived two parameters from these relationships: plant fitness when water is not limiting and the water inflection point (WIF), the watering level at which additional water is most efficiently turned into fitness. We used phylogenetic comparative methods to explore trade-offs related to drought resistance and trait plasticity and the degree to which water relationship parameters are conserved. Taxa from drier climates produced fruits at the lowest water levels, had a lower WIF, flowered earlier, had shorter life spans, had greater plastic water-use efficiency (WUE), and had lower fitness at nonlimiting water. In contrast, later-flowering Streptanthus taxa from less xeric climates experienced high fitness at nonlimiting water but had no fitness at the lowest water levels. Across the clade, we found a trade-off between drought resistance and fitness at high water, though a single ruderal species was an outlier in this relationship. Our results suggest that drought escape trades off with maximal fitness under nonlimiting water, and both are tied to phenology. We also found that variation in trait plasticity determines how different plant species produce fitness over a water gradient.

Exploring glucosinolates diversity in Brassicaceae: a genomic and chemical assessment for deciphering abiotic stress tolerance.

Brassica is one of the most economically important genus of the Brassicaceae family, encompassing several key crops like Brassica napus (cabbage) and broccoli (Brassica oleraceae var. italica). This family is well known for their high content of characteristic secondary metabolites such as glucosinolates (GLS) compounds, recognize for their beneficial health properties and role in plants defense. In this work, we have looked through gene clusters involved in the biosynthesis of GLS, by combining genomic analysis with biochemical pathways and chemical diversity assessment. A total of 101 Brassicaceae genes involved in GLS biosynthesis were identified, using a multi-database approach. Through a UPGMA and PCA analysis on the 101 GLS genes recorded, revealed a separation between the genes mainly involved in GLS core structure synthesis and genes belonging to the CYP450s and MYBs gene families. After, a detailed phylogenetic analysis was conducted to better understand the disjunction of the aliphatic and indolic genes, by focusing on CYP79F1-F2 and CYP81F1-F4, respectively. Our results point to a recent diversification of the aliphatic CYP79F1 and F2 genes in Brassica crops, while for indolic genes an earliest diversification is observed for CYP81F1-F4 genes. Chemical diversity revealed that Brassica crops have distinct GLS chemo-profiles from other Brassicaceae genera; being highlighted the high contents of GLS found among the Diplotaxis species. Also, we have explored GLS-rich species as a new source of taxa with great agronomic potential, particularly in abiotic stress tolerance, namely Diplotaxis, the closest wild relatives of Brassica crops.

Microscopic investigation of pollen morphology of Brassicaceae from Central Punjab-Pakistan.

The intention of the current study is to provide an account on the palynological features of Brassicaceae from Central Punjab-Pakistan as a basis for future studies. Different morpho-palynological characteristics both qualitative and quantitative were analyzed during this research which includes shape of pollen, diameter of pollen, P/E ratio, exine sculpturing, thickness of exine, type of pollen, shape and size of lumens, and thickness of murus. Taxonomic keys were also constructed based on pollen morphological characters for correct identification of species. This study aims to provide detailed information of pollen diversity and their exine structure based on both qualitative and quantitative characters by using Light microscopy and Scanning electron microscopy. Shape of pollen is mostly prolate, but some species also have sub-prolate to spheroidal prolate types. Exine ornamentation in most species was reticulate, whereas micro reticulate (one species) and coarsely reticulate (one species) exine also observed in some pollen. All the pollen mentioned in this study have tricolpate apertures. Variation found in thickness of exine and other characters proved to be helpful at generic and specific level. The results reinforced the significance of pollen morphological features of family Brassicaceae and aid for valuable taxonomic tool in plant systematics.

The complete chloroplast genome sequence of yellow mustard (Sinapis alba L.) and its phylogenetic relationship to other Brassicaceae species.

As a member of the large Brassicaceae family, yellow mustard (Sinapis alba L.) has been used as an important gene pool for the genetic improvement of cash crops in Brassicaceae. Understanding the phylogenetic relationship between Sinapis alba (S. alba) and other Brassicaceae crops can provide guidance on the introgression of its favorable alleles into related species. The chloroplast (cp) genome is an ideal model for assessing genome evolution and the phylogenetic relationships of complex angiosperm families. Herein, we de novo assembled the complete cp genome of S. alba by integrating the PacBio and Illumina sequencing platforms. A 153,760 bp quadripartite cycle without any gap was obtained, including a pair of inverted repeats (IRa and IRb) of 26,221 bp, separated by a large single copy (LSC) region of 83,506 bp and a small single copy (SSC) region of 17,821 bp. A total of 78 protein-coding genes, 30 tRNA genes, and four rRNA genes were identified in this cp genome, as were 89 simple sequence repeat (SSR) loci of 18 types. The codon usage analysis revealed a preferential use of the Leu codon with the A/U ending. The phylogenetic analysis using 82 Brassicaceae species demonstrated that S. alba had a close relationship with important Brassica and Raphanus species; moreover, it likely originated from a separate evolutionary pathway compared with the congeneric Sinapis arvensis. The synonymous (Ks) and non-synonymous (Ks) substitution rate analysis showed that genes encoding "Subunits of cytochrome b/f complex" were under the lowest purifying selection pressure, whereas those associated with "Maturase", "Subunit of acetyl-CoA", and "Subunits of NADH-dehydrogenase" underwent relatively higher purifying selection pressures. Our results provide valuable information for fully utilizing the S. alba cp genome as a potential genetic resource for the genetic improvement of Brassica and Raphanus species.

Diversification of Biscutella ser. Biscutella (Brassicaceae) followed post-Miocene geologic and climatic changes in the Mediterranean basin.

Biscutella ser. Biscutella (= ser. Lyratae Malin.) is a group of mostly annual or short-lived perennial plants, with petals gradually tapering at the base and lateral intrastaminal nectaries, endemic to the Mediterranean basin and the Middle East. Recent taxonomic work has revealed that a relative morphological homogeneity occurs in Europe and Asia, but a high plasticity is found in N Africa for most of the characters traditionally used for taxonomic arrangements. This fact had generally led to overestimation of the number of taxa, which currently is reduced to ten (namely 7 species and 3 additional varieties), some of them being narrow endemics. In the present contribution, on the basis of a previous detailed morphological study carried out by the authors, the first comprehensive phylogeny based on 47 DNA sequence data including concatenation of two plastid (rpl32-trnL and trnV) and one nuclear (ITS) regions, together with the first time-calibrated phylogenetic tree, allows reappraisal of evolutionary and biogeographic relationships among the accepted taxa in the series. According to all evidence gathered in the present study, the current distribution of B. ser. Biscutella, mostly centred in the southern parts of the Mediterranean basin and the Middle East, suggests that it evolved in relation with the major geological and climatic events occurred in the Mediterranean basin and Eurasia within the last 20 million years. The origin of Biscutella is dated ca 18.75 Mya, and the radiation of the series triggered ca 5.87 Mya with the Messinian Salinity Crisis. Rapid diversification occurred coetaneously to the Intensification of Northern Hemisphere Glaciation (ca 2.86 Mya) onwards, with parallel large-amplitude aridity cycles in Africa and southwestern Asia. In recent times, the divergence of lineages became faster in the W Mediterranean (ca 1.54 to 0.43 Mya), mostly related to geographical and ecological patterns of specialisation. In many cases, the distribution of the current species is apparently linked to ancient glacial refuges in S Mediterranean basin.

Spatially-resolved localization and chemical speciation of nickel and zinc in Noccaea tymphaea and Bornmuellera emarginata.

Hyperaccumulator plants present the ideal model system for studying the physiological regulation of the essential (and potentially toxic) transition elements nickel and zinc. This study used synchrotron X-ray Fluorescence Microscopy (XFM) elemental imaging and spatially resolved X-ray Absorption Spectroscopy (XAS) to elucidate elemental localization and chemical speciation of nickel and zinc in the hyperaccumulators Noccaea tymphaea and Bornmuellera emarginata (synonym Leptoplax emarginata). The results show that in the leaves of N. tymphaea nickel and zinc have contrasting localization, and whereas nickel is present in vacuoles of epidermal cells, zinc occurs mainly in the mesophyll cells. In the seeds Ni and Zn are similarly localized and strongly enriched in the cotyledons in N. tymphaea. Nickel is strongly enriched in the tip of the radicle of B. emarginata. Noccaea tymphaea has an Fe-rich provascular strand network in the cotyledons of the seed. The chemical speciation of Ni in the seeds of N. tymphaea is unequivocally associated with carboxylic acids, whereas Zn is present as the phytate complex. The spatially resolved spectroscopy did not reveal any spatial variation in chemical speciation of Ni and Zn within the N. tymphaea seed. The dissimilar ecophysiological behaviour of Ni and Zn in N. tymphaea and B. emarginata raises questions about the evolution of hyperaccumulation in these species.

Origin and Evolution of Diploid and Allopolyploid Camelina Genomes Were Accompanied by Chromosome Shattering.

Complexes of diploid and polyploid species have formed frequently during the evolution of land plants. In false flax (Camelina sativa), an important hexaploid oilseed crop closely related to Arabidopsis (Arabidopsis thaliana), the putative parental species as well as the origin of other Camelina species remained unknown. By using bacterial artificial chromosome-based chromosome painting, genomic in situ hybridization, and multi-gene phylogenetics, we aimed to elucidate the origin and evolution of the polyploid complex. Genomes of diploid camelinas (Camelina hispida, n = 7; Camelina laxa, n = 6; and Camelina neglecta, n = 6) originated from an ancestral n = 7 genome. The allotetraploid genome of Camelina rumelica (n = 13, N6H) arose from hybridization between diploids related to C. neglecta (n = 6, N6) and C. hispida (n = 7, H), and the N subgenome has undergone a substantial post-polyploid fractionation. The allohexaploid genomes of C. sativa and Camelina microcarpa (n = 20, N6N7H) originated through hybridization between an auto-allotetraploid C. neglecta-like genome (n = 13, N6N7) and C. hispida (n = 7, H), and the three subgenomes have remained stable overall since the genome merger. Remarkably, the ancestral and diploid Camelina genomes were shaped by complex chromosomal rearrangements, resembling those associated with human disorders and resulting in the origin of genome-specific shattered chromosomes.plantcell;31/11/2596/FX1F1fx1.

Viruses infecting cool season crops in the northern Turkey.

Virus diseases of cool season vegetable crops (mainly cabbage, white and red head cabbage, broccoli, kale, radish, rocket salad, garden cress, and turnip) were surveyed in Bafra Plain, Turkey during winter 2017, and 2018. Leaf samples were collected from different species of the Brassicaceae family showing mosaic, mottling, necrotic spots, malformation, and chlorosis symptoms. These samples were tested for the presence of Cauliflower mosaic virus (CaMV), Cucumber mosaic virus (CMV), Beet western yellows virus (BWYV), Radish mosaic virus (RaMV), Turnip mosaic virus (TuMV), Turnip yellow mosaic virus (TYMV), and Turnip yellows virus (TuYV) by biological and serological methods. A total of 455 samples were collected from cole crop fields and tested for the seven viruses by double-antibody sandwich ELISA using specific polyclonal antibodies. According to the results, out of these, 7 % of the samples were infected by at least one of these viruses. TuMV was the most prevalent virus detected in cole crops. TuMV, CaMV, and CMV were detected in 3 %, 2 %, and 2 % of infected samples, respectively, and the infection rate of these three viruses changed significantly among Brassica species.

Interspecies association mapping links reduced CG to TG substitution rates to the loss of gene-body methylation.

Comparative genomics can unravel the genetic basis of species differences; however, successful reports on quantitative traits are still scarce. Here we present genome assemblies of 31 so-far unassembled Brassicaceae plant species and combine them with 16 previously published assemblies to establish the Brassicaceae Diversity Panel. Using a new interspecies association strategy for quantitative traits, we found a so-far unknown association between the unexpectedly high variation in CG to TG substitution rates in genes and the absence of CHROMOMETHYLASE3 (CMT3) orthologues. Low substitution rates were associated with the loss of CMT3, while species with conserved CMT3 orthologues showed high substitution rates. Species without CMT3 also lacked gene-body methylation (gbM), suggesting an evolutionary trade-off between the unknown function of gbM and low substitution rates in Brassicaceae, possibly due to low mutability of non-methylated cytosines.