diff --git a/DESCRIPTION b/DESCRIPTION
index ff7678c4a..03954b96c 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -2,7 +2,7 @@ Package: rstanarm
 Type: Package
 Title: Bayesian Applied Regression Modeling via Stan
 Version: 2.21.1
-Date: 2020-07-07
+Date: 2020-07-20
 Encoding: UTF-8
 Authors@R: c(person("Jonah", "Gabry", email = "jsg2201@columbia.edu", role = "aut"),
              person("Imad", "Ali", role = "ctb"),
diff --git a/NEWS.md b/NEWS.md
index 103e3b671..46b548637 100644
--- a/NEWS.md
+++ b/NEWS.md
@@ -1,3 +1,27 @@
+# rstanarm 2.21.1
+
+* Consistent with [Regression and Other Stories](https://statmodeling.stat.columbia.edu/2020/07/08/regression-and-other-stories-is-available/)
+
+### Backwards incompatible changes
+
+* `autoscale` argument to various prior functions not defaults to `FALSE`, although 
+  by default `normal` is now called with `autoscale = TRUE` in `stan_glm`, `stan_glmer`, etc.
+* The default prior on the intercept is different than it was in rstanarm <= 2.19.3
+  for `stan_glm`, `stan_glmer`, etc.
+* `stan_jm` is not available for 32bit Windows
+
+### New functions
+
+* `posterior_epred` returns the posterior distribution of the conditional expectation,
+  which is previously accomplished via `posterior_linpred` with `transform = TRUE`
+* `predict` produces predictions in more cases where it previously threw errors
+
+### Bug fixes
+
+* `singular.ok` now rules out singular design matrices in `stan_lm`
+* `newdata` now works when the family was `mgcv::betar`
+* now works better with `data.table`s
+
 # rstanarm 2.19.3
 
 ### Bug fixes