diff --git a/DESCRIPTION b/DESCRIPTION
index 4b0894ba2deb99d357853925b4dfe62dc2d8c3a5..e6f968eb632c1ffc740cab137fa642355ed4eb97 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,8 +1,8 @@
 Package: stressaddition
 Type: Package
 Title: Modelling Tri-Phasic Concentration-Response Relationships
-Version: 3.0.2
-Date: 2020-09-09
+Version: 3.0.3
+Date: 2020-09-21
 Authors@R: c(person("Sebastian", 
                     "Henz", 
                     role = c("aut", "cre"), 
@@ -22,7 +22,8 @@ Description: The stress addition approach is an alternative to the traditional
     concentration addition or effect addition models. It allows the modelling 
     of tri-phasic concentration-response relationships either as single toxicant 
     experiments, in combination with an environmental stressor or as mixtures of 
-    two toxicants. See Liess et al. (2019) <doi:10.1038/s41598-019-51645-4>.
+    two toxicants. See Liess et al. (2019) <doi:10.1038/s41598-019-51645-4> 
+    and Liess et al. (2020) <doi:10.1186/s12302-020-00394-7>.
 License: GPL-3
 Copyright: file inst/COPYRIGHTS
 URL: https://git.ufz.de/oekotox/stressaddition
diff --git a/NEWS.md b/NEWS.md
index 3e75fb8de986ea801c250a63afdd411207a78571..beb7e46ddad9cb1b8a5041e60d198e8341a4b35c 100644
--- a/NEWS.md
+++ b/NEWS.md
@@ -1,3 +1,7 @@
+# stressaddition 3.0.3
+
+* Added references to the journal article about the Multi-Tox model which was recently published.
+
 # stressaddition 3.0.2
 
 * Fixed a bug where the plotting functions printed `NULL` to the console.
diff --git a/R/data.R b/R/data.R
index eb8c1d8bd4dd1f2dbc27b45a6fa0bc299362c0f9..0a4c6cbb2a85a3ead544bfdf40b1cbc32293fa61 100644
--- a/R/data.R
+++ b/R/data.R
@@ -41,6 +41,7 @@
 #'     of Pesticide Mixtures on Daphnia magna. Environ. Sci. Technol. 53,
 #'     12586–12593.}
 #'
-#'   Liess, M., Henz, S., Shahid, N., 2020. Modelling the synergistic effects
-#'     of toxicant mixtures. Manuscript submitted for publication.
+#'   \href{https://doi.org/10.1186/s12302-020-00394-7}{Liess, M., Henz, S. &
+#'   Shahid, N. Modeling the synergistic effects of
+#'   toxicant mixtures. Environ Sci Eur 32, 119 (2020).}
 "multiple_stress"
diff --git a/R/multi_tox.R b/R/multi_tox.R
index 97ae971924e26d36c6e8ca07f34c87d46583c292..a7d0a94fa6c01e63114865797c14ef7b5b1d8502 100644
--- a/R/multi_tox.R
+++ b/R/multi_tox.R
@@ -86,8 +86,9 @@
 #' mt[1:3]  # The concentrations and survival of the 4 mixtures.
 #'
 #'
-#' @references Liess, M., Henz, S., Shahid, N., 2020. Modelling the synergistic
-#'   effects of toxicant mixtures. Manuscript submitted for publication.
+#' @references \href{https://doi.org/10.1186/s12302-020-00394-7}{Liess, M.,
+#'   Henz, S. & Shahid, N. Modeling the synergistic effects of
+#'   toxicant mixtures. Environ Sci Eur 32, 119 (2020).}
 #'
 #' @export
 multi_tox <- function(model_a,
diff --git a/R/stressaddition-package.R b/R/stressaddition-package.R
index c761fcb32f8da660b4885c9c34d1778a9a0114a7..2ac20d4ba949000425572e82b31c7c840beb1f74 100644
--- a/R/stressaddition-package.R
+++ b/R/stressaddition-package.R
@@ -39,8 +39,9 @@
 #'   Henz, S. & Knillmann, S. Predicting low-concentration effects of
 #'   pesticides. Sci Rep 9, 15248 (2019).}
 #'
-#'   Liess, M., Henz, S., Shahid, N., 2020. Modelling the synergistic effects of
-#'   toxicant mixtures. Manuscript submitted for publication.
+#'   \href{https://doi.org/10.1186/s12302-020-00394-7}{Liess, M., Henz, S. &
+#'   Shahid, N. Modeling the synergistic effects of
+#'   toxicant mixtures. Environ Sci Eur 32, 119 (2020).}
 #'
 #' @import stats
 #' @import graphics
diff --git a/README.md b/README.md
index 05c4136f52068caaa2e87a682cb7955a60de490c..0f368387c6349463f7760562dee3b276e2c8f360 100644
--- a/README.md
+++ b/README.md
@@ -4,7 +4,7 @@ This R package makes it possible to model tri-phasic concentration-response rela
 `ecxsys()` implements **EC<sub>x-SyS</sub>**, the tri-phasic concentration-response model introduced in
 [Liess, M., Henz, S. & Knillmann, S. Predicting low-concentration effects of pesticides. Sci Rep 9, 15248 (2019)](https://doi.org/10.1038/s41598-019-51645-4). It is applicable to modelling ecotoxicological experiments with and without environmental stress where the response contains a hormesis effect.
 
-`multi_tox()` implements **Multi-TOX**, a model for binary mixtures of toxicants where each toxicant exhibits a tri-phasic concentration-response relationship. See *Liess, M., Henz, S., Shahid, N. (2020), Modelling the synergistic effects of toxicant mixtures. Manuscript submitted for publication*.
+`multi_tox()` implements **Multi-TOX**, a model for binary mixtures of toxicants where each toxicant exhibits a tri-phasic concentration-response relationship. See [Liess, M., Henz, S. & Shahid, N. Modeling the synergistic effects of toxicant mixtures. Environ Sci Eur 32, 119 (2020)](https://doi.org/10.1186/s12302-020-00394-7).
 
 The EC<sub>x-SyS</sub> and Multi-TOX models are also available as part of the [Indicate app](http://www.systemecology.eu/indicate) which offers an easy to use graphical user interface.
 
diff --git a/man/multi_tox.Rd b/man/multi_tox.Rd
index 7abbb9fd1f38b4ab7559e596f676fd4bc1cae8aa..47994cc34a89744402eaa5ef9193e1d7803ee65c 100644
--- a/man/multi_tox.Rd
+++ b/man/multi_tox.Rd
@@ -89,6 +89,7 @@ mt[1:3]  # The concentrations and survival of the 4 mixtures.
 
 }
 \references{
-Liess, M., Henz, S., Shahid, N., 2020. Modelling the synergistic
-  effects of toxicant mixtures. Manuscript submitted for publication.
+\href{https://doi.org/10.1186/s12302-020-00394-7}{Liess, M.,
+  Henz, S. & Shahid, N. Modeling the synergistic effects of
+  toxicant mixtures. Environ Sci Eur 32, 119 (2020).}
 }
diff --git a/man/multiple_stress.Rd b/man/multiple_stress.Rd
index 8667622d5138803f6d0c1c33163ac7dad14ea3ea..1042d71fc23668247b21d1c18c3825cc9a587fd7 100644
--- a/man/multiple_stress.Rd
+++ b/man/multiple_stress.Rd
@@ -19,8 +19,9 @@ A data frame with 58 rows and 4 variables:
     of Pesticide Mixtures on Daphnia magna. Environ. Sci. Technol. 53,
     12586–12593.}
 
-  Liess, M., Henz, S., Shahid, N., 2020. Modelling the synergistic effects
-    of toxicant mixtures. Manuscript submitted for publication.
+  \href{https://doi.org/10.1186/s12302-020-00394-7}{Liess, M., Henz, S. &
+  Shahid, N. Modeling the synergistic effects of
+  toxicant mixtures. Environ Sci Eur 32, 119 (2020).}
 }
 \usage{
 multiple_stress
diff --git a/man/stressaddition-package.Rd b/man/stressaddition-package.Rd
index 223841f7f118de2b6065b124306e7a21602de48c..224c8fcac0f15f93ecee4bb53715aabbe692885b 100644
--- a/man/stressaddition-package.Rd
+++ b/man/stressaddition-package.Rd
@@ -10,7 +10,8 @@ The stress addition approach is an alternative to the traditional
     concentration addition or effect addition models. It allows the modelling 
     of tri-phasic concentration-response relationships either as single toxicant 
     experiments, in combination with an environmental stressor or as mixtures of 
-    two toxicants. See Liess et al. (2019) <doi:10.1038/s41598-019-51645-4>.
+    two toxicants. See Liess et al. (2019) <doi:10.1038/s41598-019-51645-4> 
+    and Liess et al. (2020) <doi:10.1186/s12302-020-00394-7>.
 }
 \details{
 This R package contains the definitions of the
@@ -26,8 +27,9 @@ This R package contains the definitions of the
   Henz, S. & Knillmann, S. Predicting low-concentration effects of
   pesticides. Sci Rep 9, 15248 (2019).}
 
-  Liess, M., Henz, S., Shahid, N., 2020. Modelling the synergistic effects of
-  toxicant mixtures. Manuscript submitted for publication.
+  \href{https://doi.org/10.1186/s12302-020-00394-7}{Liess, M., Henz, S. &
+  Shahid, N. Modeling the synergistic effects of
+  toxicant mixtures. Environ Sci Eur 32, 119 (2020).}
 }
 \seealso{
 Useful links: