test loading and parsing of Recon3D

test/data_downloaded.jl

@@ -63,4 +63,9 @@ model_paths = Dict{String,String}(
         df("iJR904.mat"),
         "d17be86293d4caafc32b829da4e2d0d76eb45e1bb837e0138327043a83e20c6e",
     ),
+    "Recond3D.json" => download_data_file(
+        "http://bigg.ucsd.edu/static/models/Recon3D.json",
+        df("Recon3D.json"),
+        "aba925f17547a42f9fdb4c1f685d89364cbf4979bbe7862e9f793af7169b26d5",
+    ),
 )

test/io/io.jl

@@ -16,4 +16,11 @@
           Set("r_" .* lowercase.(reactions(matlabmodel)))
     @test Set(lowercase.(reactions(sbmlmodel))) ==
           Set("r_" .* lowercase.(reactions(jsonmodel)))
+
+    # specifically test parsing of gene-reaction associations in Recon
+    reconmodel = load_model(model_paths["Recon3D.json"])
+    @test n_reactions(reconmodel) == 2583
+    recon_grrs = [r.grr for (i, r) in m.reactions if !isnothing(r.grr)]
+    @test length(recon_grrs) == 5938
+    @test sum(length.(recon_grrs)) == 13903
 end