The schema describes our use of parallel (A) and convergent (B) evolution for within-species chemotypic variation. The letters in the blue box represent the state of the source/ancestral haplotypes. …
Short names and structures of the GSLs are in black. Genes encoding the causal enzyme for each reaction (arrow) are in gray. GS-OX is a gene family of five or more genes. OH-But: 2-OH-3-Butenyl.
(A) The accessions are plotted on the map based on their collection site and colored based on their principal component (PC)1 score. (B) Manhattan plot of genome-wideassociation analyses using PC1. …
(A) Percentage of variance explained by each PC. (B, C) Contribution of the individual GSLs to PC1 (B) and PC2 (C). Red bars: contribution of four carbon GSLs; blue bars: contribution of three …
(A) The accessions were plotted on the map based on their collection site and colored based on their principal component (PC)2 score. (B) Manhattan plot of genome-wideassociation analyses using PC2. …
Horizontal lines represent 5% significance thresholds using Bonferroni (red) and Benjamini–Hochberg (blue).
(A) Using the GSL accumulation, each accession was classified to one of seven aliphatic short-chained GSL chemotypes based on the enzyme functions as follows: MAM2, AOP null: classified as 3MSO …
Environmental conditions differentially associate with MAM status in the north versus the south.
Linear model for MAM status (carbon side chain length) was conducted with the indicated environmental parameters, for the northern and southern collection, separately (for more details, see Methods). The tables show p values for each term from the linear model. For the interaction with geography, the linear model was run using the total dataset, and the geography parameter (north or south) was added to the model.
Accessions were classified based on the side chain length of the aliphatic short-chained glucosinolates (GSLs). Accessions with a majority of GSLs containing three carbons in their side chains are …
Relative amounts of alkenyl glucosinolates (GSLs), alkyl GSLs and methylsulfinyl (MSO) GSLs were calculated in respect to the total short-chained aliphatic GSLs as described in the Methods section. …
The ratio between 2-OH-3-Butenyl to 3-Butenyl glucosinolate (GSL) was calculated only for MAM1-dominant accessions (accessions with GSLs containing four carbons in their side chain). Accessions with …
(A) The accessions were divided to two collections using the following chain of mountains: the Pyrenees between Spain and France, the Alps between Italy and Germany, and the Carpathians in the …
(A) MAM3 phylogeny of Arabidopsis thaliana accessions, rooted by Arabidopsis lyrata MAMb, which is not shown because of distance. Tree tips are colored based on the accession chemotype. (B) The …
(A) MAM3 phylogeny of 637 Arabidopsis thaliana accessions, rooted by Arabidopsis lyrata MAMb, excluding accessions with low-quality sequences. (B–E) MAM3 phylogeny of different combinations of A. …
The GS-Elong locus from different accessions was sequenced, and the MAM1 and MAM2 structures were analyzed. The table indicates the dominant chemotype of each accession, the MAM status of each …
Sequences of MAM locus.
Sequencing for the accessions in Figure 5—figure supplement 2 were generated. The extracted regions are 30000bps upstream of AT5G23000 (MYB37) and 60000 bps downstream of AT5G23020 (MAM3).
Domain (A) and full sequence (B) amino acid phylogenies of the MAM/IPMS gene family. Sequences were taken from Abrahams et al., 2020, which uses Arabidopsis thalina Col-0 genome and the MAM2 amino …
(A) All accessions from Iberia were plotted, colored and shaped based on the side chain length of the aliphatic short chained GSLs. Accessions with a majority of GSLs containing 3 carbons in their …
The MAM phylogeny is split by the major clades/haplotypes and each sub-clade’s phylogeny is reflected on the map. Tree tips are colored based on the accessions chemotype.
The genomic structure and causality of the major AOP2/AOP3 haplotypes are illustrated. Pink arrows show the AOP2 gene while yellow arrows represent AOP3. The black arrows represent the direction of …
Linear model for the two major chemotypes, Allyl and 2-OH-3-Butenyl, was conducted with the indicated environmental parameters, for the northern and southern collection, separately (for more …
Environmental parameter | Effect on chemotype – north | Effect on chemotype – south | Interaction with geography |
---|---|---|---|
Genomic group | <0.0001 | <0.0001 | <0.0001 |
Max temperature of warmest month | 0.0382 | <0.0001 | 0.3574 |
Min temperature of coldest month | 0.0007 | <0.0001 | 0.0049 |
Precipitation of wettest month | 0.1645 | 0.0003 | 0.0094 |
Precipitation of driest month | 0.0665 | 0.2026 | 0.47425 |
Distance to the coast | 0.2781 | 0.02680 | 0.1279 |
The structures of GS-OH in the 3-Butenyl accessions are illustrated. Gray boxes represent exons, and blue lines represent introns. Black line represents a mutation, and gray lines represent unknown …
Accession | Type of mutation | Allele structure | Fraction (out of C4 Alkenyl accessions) | Observed frequency (out of non-C4 Alkenyl accessions) |
---|---|---|---|---|
Sorbo, Pien | Polymorphism at SNP10831302 | ![]() | 0.009 (2/226) | 0.067 (38/564) |
Cvi-0 | Active site mutation | ![]() | 0.004 (1/226) | 0.025 (14/564) |
IP-Mot-0, IP-Tri-0 | Gene deletion | ![]() | 0.009 (2/226) | 0.055 (31/564) |
Multiple accessions (T670, FlyA-3, Ting-1, T880, T710, T850) | Unidentified mutations | ![]() | 0.026 (6/226) | Unknown |
GSLs data.
(A) List of GSLs and structures. (B) Accessions and glucosinolate (GSL) data – raw data. (C) Heritability values. (D) Accessions and GSL data – emmeans.
SNPs in glucosinolates (GSLs)-related loci under different genome-wide association (GWA) studies: GSL values were used as traits to conduct GWA studies.
The number of significant SNPs in the GSLs related loci (columns c to o) was counted for each GWA study separately. Rows 2–3: common name and AT number of gene/s in the loci. Rows 4–5: upstream and downstream positions of the relevant loci (10 kb were added upstream and downstream of the genes). Rows 6–33: GSLs traits used for GWA studies. In black: number of SNPs with p value between 0.00001 and 0.0000001. In red: number of SNPs with p value equal or smaller than 0.0000001.