Weak evidence for heritable changes in response to selection by aphids in Arabidopsis accessions

  1. Marc W Schmid
  2. Klara Kropivšek
  3. Samuel E Wuest
  4. Bernhard Schmid
  5. Ueli Grossniklaus  Is a corresponding author
  1. MWSchmid GmbH, Switzerland
  2. Department of Plant and Microbial Biology, University of Zurich, Switzerland
  3. Zurich-Basel Plant Science Center, University of Zurich, ETH Zurich, University of Basel, Switzerland
  4. Department of Evolutionary Biology and Environmental Studies, University of Zurich, Switzerland
  5. Department of Geography, University of Zurich, Switzerland
5 figures, 3 tables and 5 additional files

Figures

Schematic of the original experimental design (top) and the current study (bottom).
Phenotypes measured in the G7 generation averaged by the ID of the maternal lines from the G6 generation.

p-values for contrasts with p<0.05 are indicated (based on the analyses done per accession). p-values with asterisk were not significant when outliers were removed (see Figure 3). Boxplots were drawn with the function graphics::boxplot() in R. The box extends from the 25th to the 75th percentile. The wiskers extend to the most extreme data point which is no more than 1.5 times the interquartile range from the box.

Phenotypes measured in the G7 generation.

Brown and green colors represent plants from blocks 1–3 and 4–5, respectively (slower growth in blocks 4–5). Dots with black outlines were identified as outliers. The number of plants available for the test differed slightly between the measurements. The number of plants given at the bottom of the figure corresponds to the minimal number of plants (i.e., for each phenotype there were at least Nplants plants). It varies by 2–3 plants per accession and treatment as some measures were missing. p-values for contrasts with p<0.05 are indicated (based on the analyses done per accession). p-values with asterisks were not significant when outliers were removed. Boxplots were drawn with the function graphics::boxplot() in R. The box extends from the 25th to the 75th percentile. The wiskers extend to the most extreme data point which is no more than 1.5 times the interquartile range from the box.

Epigenetic differences between ancestral and selected populations using a subset of mothers from the accession Wei-0.

(a) Location of the differentially methylated cytosines along the genome. (b) Epigenetic variation within different populations (selection replicates), i.e., distances to each centroid. Only populations with at least four individuals are shown. Number of individuals per population: A0: 7, BB19: 5, BB27: 4, LE4: 4, LE13: 4, LE28: 4 (not shown: BB14 and BB16 with only two individuals each). No overall and no pairwise p-value was smaller than 0.05. Boxplots were drawn with the function graphics::boxplot() in R. The box extends from the 25th to the 75th percentile. The wiskers extend to the most extreme data point which is no more than 1.5 times the interquartile range from the box.

Genetic differentiation between some treatments in the study of Lin et al., 2024.

(a) Genetic distances between the 190 RNA-seq samples from Lin et al., 2024. The two genetic clusters diverging from the rest of their accession group are highlighted. (b) Within accession-label genetic distances, drawn to scale. Maximal distance given below the label, colored according to treatment and generation. The two treatments that have been assigned the wrong accessions are highlighted. (c) Pairwise genetic distances within lines (plants of the same lineage, i.e., single-seed descendants), within accessions, and between accessions using either the original accession labels from Lin et al., 2024 (left; 26) or the accessions assigned based on genetic similarity to the data available at 1001genomes.org (right). The red line demarcates the distance equal to the average of within-line distances plus three standard deviations. Boxplots were drawn with the function graphics::boxplot() in R. The box extends from the 25th to the 75th percentile. The wiskers extend to the most extreme data point which is no more than 1.5 times the interquartile range from the box.

Tables

Table 1
Percent sum of squares explained by model terms and p-values using all three accessions.

Contrasts of a term are marked with indents. Complete output and further results for individual accessions are given in Supplementary files 1 and 2.

Blade areaTime to boltingRosette diameterNumber of rosette leaves
Model term% SSp-value% SSp-value% SSp-value% SSp-value
Blocks46.97<0.0016.55<0.00143.23<0.00127.1<0.001
>Blocks 1–3 vs. blocks 4–5 (EX)41.85<0.0010.51<0.00139.62<0.00126.42<0.001
>Remaining blocks5.12<0.0016.04<0.0013.61<0.0010.68<0.001
Selection treatments (ST)0.10.8870.240.9250.760.0330.50.782
>Ancestral vs. rest (AN)0.020.7420.060.7350.110.2220.240.484
>Control vs. treated (CO)0.070.4950.050.7570.330.0430.210.506
>B. brassicae vs. L. erysimi (SE)0.010.8520.130.6240.320.0470.050.746
Original replicate population (LI)2.38<0.0018.82<0.0011.18<0.0017.88<0.001
Accession (ET)8.16<0.00123.23<0.0012.2<0.00114.7<0.001
ET:ST0.420.0260.340.410.620.0020.160.655
>ET:AN0.050.3150.090.4270.350.0010.120.238
>ET:CO0.060.2880.220.1530.010.840.010.843
>ET:SE0.310.0050.030.7760.260.0050.030.707
ET:LI0.370.4430.880.0760.30.8640.660.175
Mothers6.070.0078.110.2167.380.0246.810.342
Table 2
Percent sum of squares explained by model terms and p-values using only accessions Ma-0 and Wei-0.

Contrasts of a term are marked with indents. Complete output and further results for individual accessions are given in Supplementary files 1 and 2.

Blade areaTime to boltingRosette diameterNumber of rosette leavesTrichome density
Model term% SSp-value% SSp-value% SSp-value% SSp-value% SSp-value
Blocks48.54<0.0015.91<0.00142.45<0.00129.24<0.00128.97<0.001
>Blocks 1–3 vs. blocks 4–5 (EX)43.45<0.0010.280.02638.73<0.00127.94<0.00125.28<0.001
>Remaining blocks5.09<0.0015.63<0.0013.72<0.0011.3<0.0013.69<0.001
Selection treatments (ST)0.060.9260.530.8270.190.6130.190.9410.180.901
>Ancestral vs. rest (AN)00.9550.170.59700.8880.110.64100.936
>Control vs. treated (CO)00.9080.210.5610.10.3380.070.7150.020.801
>B. brassicae vs. L. erysimi (SE)0.060.5150.150.6280.090.3670.010.8740.160.49
Original replicate population (LI)2.16<0.0019.5<0.0011.670.0027.96<0.0014.99<0.001
Accession (ET)5.78<0.00118.5<0.0012.46<0.00111.78<0.00113.98<0.001
ET:ST0.120.4240.420.3490.380.0360.040.7990.070.572
>ET:AN0.020.4710.190.2250.130.060.020.510.040.304
>ET:CO0.030.3940.180.2410.010.5880.020.5220.020.422
>ET:SE0.070.2220.050.5150.240.01700.7750.010.66
ET:LI0.340.2410.890.0480.210.7750.30.6090.280.597
Mothers7.690.02311.270.3199.160.210.20.1068.250.46
Table 3
Gene Ontology (GO) enrichment among genes mapped by differentially methylated cytosines (DMCs).

‘Found’ indicates the number of different genes with the term that were mapped by DMCs. ‘Expected’ refers to the number of genes that were expected to be mapped by DMCs if six DMCs were randomly distributed.

GO.IDTermFoundExpectedp-value
GO:0070646Protein modification by small protein removal71.530.00083
GO:0016571Histone methylation103.070.00101
GO:0008284Positive regulation of cell proliferation50.80.00119
GO:0042732D-xylose metabolic process30.30.00321
GO:0006487Protein N-linked glycosylation51.030.0036
GO:0071702Organic substance transport2511.890.00423
GO:0015691Cadmium ion transport20.110.00455
GO:0009834Plant-type secondary cell wall biogenesis30.350.00482
GO:0006306DNA methylation72.090.00492
GO:0006342Chromatin silencing82.860.00794
GO:0071705Nitrogen compound transport2311.740.00931
GO:0050665Hydrogen peroxide biosynthetic process40.830.00945
GO:1905039Carboxylic acid transmembrane transport20.150.00945

Additional files

Supplementary file 1

Collection of workbooks with results from the different analyses.

The files are named as <phenotype > _<data set>.xlsx. For each phenotype, there are analyses with all three accessions (‘all’), only the two accessions with trichomes (‘NoWil3’), or with only one accession (‘OnlyWei0’, ‘OnlyMa0’, and ‘OnlyWil3’). For each of them, there is one version with all data included and one in which the extreme values were removed (‘.outliersRemoved’). For example: bladeArea_all.outliersRemoved.xlsx refers to the analysis of the blade area using all three accessions and removing extreme values. Within each workbook, there are three sheets for untransformed data (‘y’), log-transformed data (‘logy’), and square-root-transformed data (‘sqrty’). Within each sheet, the output of the two models tested either with the regular linear model (lm) or the random effects model (lmer) is given. In case of the regular lm output, alternative F- and p-values (F_altern and P_altern) refer to the tests against the correct error stratum (i.e., nested model). The columns F and P contain the F- and p-values from tests against the residuals. Regarding the model terms: ‘__Blocks’ is the five-level factor of experimental blocks that was split into two contrasts (‘EX’ and ‘BL’); ‘__Selection’ is the four-level factor with the selection treatments that was split into three contrasts (‘AN’, ‘CO’, and ‘SE’). ‘__Selection_X_accession/experiment’ is accordingly the interaction of the four-level factor with the selection treatments and the accession (‘ET’) or the ‘EX’ contrast. Abbreviations: Df: degrees-of-freedom, SS: sum of squares, MS: mean squares, F: F-value, P: p-value, F_altern: F-value of nested model, P_altern: p-value of nested model, perc_SS: % sum of squares explained by the term.

https://cdn.elifesciences.org/articles/106930/elife-106930-supp1-v1.zip
Supplementary file 2

A workbook containing tables from the power estimations.

For each phenotype, power was estimated using all three accessions or only the two accessions with trichomes. In case of trichome density, only the latter is available. Each table contains the percentage of significant outcomes for a given simulation and factor/contrast of interest (based on 1,000 simulations). The simulations are named as <type > _<effect size>. There were five different types of simulations with effect sizes ranging from 1% to 20%: (‘step’) differences between all four groups with step sizes of 1%, 5%, 10%, 15%, and 20% between each group (with the level following AN<CO<BB<LE, note that the largest difference here is 60%); (‘onlyAn’) only different in ancestors; (‘onlyCo’) only different in controls; (‘onlySe’) different in both selection treatments; (‘onlyLe’) different in only one selection treatment.

https://cdn.elifesciences.org/articles/106930/elife-106930-supp2-v1.xlsx
Supplementary file 3

Table with the samples from Lin et al., 2024.

Given are the sample IDs, line ID, treatment, generation, accession label (used by Lin et al., 2024), and the accession assignments based on the genotypes inferred from the RNA-seq data and the publicly available data from 1001 Arabidopsis accessions (this study).

https://cdn.elifesciences.org/articles/106930/elife-106930-supp3-v1.csv
Supplementary file 4

Table with the phenotype data.

Rosette diameter in mm, blade area in mm2. Factors are labeled as described in ‘Materials and methods’.

https://cdn.elifesciences.org/articles/106930/elife-106930-supp4-v1.xlsx
MDAR checklist
https://cdn.elifesciences.org/articles/106930/elife-106930-mdarchecklist1-v1.docx

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Marc W Schmid
  2. Klara Kropivšek
  3. Samuel E Wuest
  4. Bernhard Schmid
  5. Ueli Grossniklaus
(2025)
Weak evidence for heritable changes in response to selection by aphids in Arabidopsis accessions
eLife 14:RP106930.
https://doi.org/10.7554/eLife.106930.3