diff --git a/src/reactionsystem.jl b/src/reactionsystem.jl
index a60bc2b818..7991a6053e 100644
--- a/src/reactionsystem.jl
+++ b/src/reactionsystem.jl
@@ -607,7 +607,8 @@ function Base.convert(::Type{<:ODESystem}, rs::ReactionSystem;
                              include_zero_odes=include_zero_odes)                                 
     eqs,sts,ps = addconstraints!(eqs, fullrs)
     ODESystem(eqs, get_iv(fullrs), sts, ps; name=name, defaults=get_defaults(fullrs), 
-                                            checks=checks, kwargs...)
+                                            observed=get_observed(fullrs), checks=checks, 
+                                            kwargs...)
 end
 
 """
@@ -632,7 +633,8 @@ function Base.convert(::Type{<:NonlinearSystem}, rs::ReactionSystem;
     error_if_constraint_odes(NonlinearSystem, fullrs)
     eqs,sts,ps = addconstraints!(eqs, fullrs)
     NonlinearSystem(eqs, sts, ps; name=name, defaults=get_defaults(fullrs), 
-                                  checks = checks, kwargs...)
+                                  observed=get_observed(fullrs), checks = checks, 
+                                  kwargs...)
 end
 
 """
@@ -682,6 +684,7 @@ function Base.convert(::Type{<:SDESystem}, rs::ReactionSystem;
               (noise_scaling===nothing) ? get_ps(flatrs) : union(get_ps(flatrs), toparam(noise_scaling));
               name=name, 
               defaults=get_defaults(flatrs),
+              observed=get_observed(fullrs), 
               checks = checks,
               kwargs...)
 end
@@ -707,7 +710,8 @@ function Base.convert(::Type{<:JumpSystem},rs::ReactionSystem;
 
     eqs = assemble_jumps(flatrs; combinatoric_ratelaws=combinatoric_ratelaws)
     JumpSystem(eqs, get_iv(flatrs), get_states(flatrs), get_ps(flatrs); name=name, 
-               defaults=get_defaults(flatrs), checks = checks, kwargs...)
+               defaults=get_defaults(flatrs), observed=get_observed(fullrs), 
+               checks = checks, kwargs...)
 end