diff --git a/R/rmReact.R b/R/rmReact.R
index 3e9ff3c5d53bbd48cc31e0b9cd673ecc884bfa9f..9b30a44388c8af1c7788484aa72fc8ea773e845e 100644
--- a/R/rmReact.R
+++ b/R/rmReact.R
@@ -169,7 +169,7 @@ rmReact <- function(model, react, rm_met = TRUE) {
       keepGenes <- ifelse(is.na(keepGenes), FALSE, TRUE)
       #print(keepGenes)
 
-      rxnGeneMat(mod_out)   <- rxnGeneMat(mod_out)[, keepGenes]
+      rxnGeneMat(mod_out)   <- rxnGeneMat(mod_out)[, keepGenes, drop = FALSE]
       #print(dim(rxnGeneMat))
   }
   
diff --git a/inst/NEWS.Rd b/inst/NEWS.Rd
index 0b5d4ac1358f939ad9c0edfa0989ad38f6d0782b..9eb40dc87d2b1237e4718606d8acc07c0d976755 100644
--- a/inst/NEWS.Rd
+++ b/inst/NEWS.Rd
@@ -8,6 +8,7 @@
 
 \section{Changes in version 1.3.1 2015-08-06}{
   \itemize{
+  	\item \code{rmReact} error, if resulting model had only one gene left, corrected.
     \item \code{deadEndMetabolite()} could miss deadEndMetabolites if reactions differ in stoichiometry e.g. 1 and 2.
   }
 }