diff --git a/.gitattributes b/.gitattributes
index 8cece8998908ac7bb479143e780cd9eafb2b1edf..6d4654fb03dc0b97ba836ce39c9579ca219fa7d6 100644
--- a/.gitattributes
+++ b/.gitattributes
@@ -1,2 +1,3 @@
 **/dataset/** filter=lfs diff=lfs merge=lfs -text
 /_publication/De[[:space:]]Jager[[:space:]]et[[:space:]]al.[[:space:]]-[[:space:]]2024[[:space:]]-[[:space:]]Traits[[:space:]]linked[[:space:]]to[[:space:]]natural[[:space:]]variation[[:space:]]of[[:space:]]sulfur[[:space:]]conte.pdf filter=lfs diff=lfs merge=lfs -text
+_publication/erad401_suppl_supplementary_figures_s1-s12.pdf filter=lfs diff=lfs merge=lfs -text
diff --git a/_publication/README.md b/_publication/README.md
new file mode 100644
index 0000000000000000000000000000000000000000..ac22d98cf3af10d20ff781b98d7d0ce3bb5019da
--- /dev/null
+++ b/_publication/README.md
@@ -0,0 +1,42 @@
+## Supplementary Data
+
+The following supplementary data are available at JXB online. 
+ https://academic.oup.com/jxb/article/75/3/1036/7313479#438654548
+
+**erad401_suppl_Supplementary_Figures_S1-S12 - pdf file**
+
+Fig. S1. Natural variation of S content in leaves of Arabidopsis thaliana accessions.
+
+Fig. S2. Quantification of anion content in leaves of three groups of Arabidopsis thaliana accessions, based on the S content in leaf.
+
+Fig. S3. Sulfate uptake traits in Arabidopsis thaliana accessions.
+
+Fig. S4. Linear regression analysis on the relationship of S uptake in opposite groups of Arabidopsis thaliana accessions, based on the S content in leaves.
+
+Fig. S5. Clustering analysis was performed on 383 Arabidopsis thaliana nutrient-related genes in seven contrasting accessions.
+
+Fig. S6. RT-qPCR analysis in 14 A. thaliana contrasting accessions, based on the S content in leaves.
+
+Fig. S7. Quantification of cysteine (Cys) and gluthatione (GSH) in leaves in A. thaliana contrasting accessions, based on the S content in leaves.
+
+Fig. S8. Quantification of glucosinolates (GSLs) in leaves of A. thaliana contrasting accessions, based on the S content in leaves.
+
+Fig. S9. Linear regression analysis on relationship of S uptake and S metabolites.
+
+**erad401_suppl_Supplementary_Tables_S1-S8 - xlsx file**
+
+Table S1. Ionomics data for 174 overlapping accessions between Baxter et al. (2007) and Campos et al. (2021) were used to perform PCA, and select three groups of accessions with low, mid, and high S content.
+
+Table S2. RT-qPCR primers used.
+
+Table S3. Loading matrix and eigenvalues of the first 10 principal components (Prin1–10). Ionomics data for 174 overlapping accessions between Baxter et al. (2007) and Campos et al. (2021) were used to perform PCA.
+
+Table S4. Z-scores of 383 nutrient-related genes in seven accessions with contrasting S content in leaves (expression values transformed from Kawakatsu et al., 2016).
+
+Table S5. P-values between low and high S-containing accession groups, for significantly different GSLs. t-test (P-value <0.05).
+
+Table S6. Correlation matrix of quantified mean traits.
+
+Table S7. P-value matrix of quantified mean traits.
+
+Table S8. Pairwise correlations of all quantified mean traits.
\ No newline at end of file
diff --git a/_publication/erad401_suppl_supplementary_figures_s1-s12.pdf b/_publication/erad401_suppl_supplementary_figures_s1-s12.pdf
new file mode 100644
index 0000000000000000000000000000000000000000..4e22805e0bdf77251113f6274e6c91adeb773a21
--- /dev/null
+++ b/_publication/erad401_suppl_supplementary_figures_s1-s12.pdf
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:5a31e698a04f4502c3eddb0d707b2b4975171dc22fddb7b2527bb88863d0ef4e
+size 2563571
diff --git a/_publication/erad401_suppl_supplementary_tables_s1-s8.xlsx b/_publication/erad401_suppl_supplementary_tables_s1-s8.xlsx
new file mode 100644
index 0000000000000000000000000000000000000000..a22fa2aaf8db76474c63260612f451ddb39af5e8
Binary files /dev/null and b/_publication/erad401_suppl_supplementary_tables_s1-s8.xlsx differ