Modeling genetic regulatory networks is an important research topic in genomic research and computational systems biology. This paper considers the problem of constructing a genetic regulatory network (GRN) using the discrete dynamic system (DDS) model approach. Although considerable research has been devoted to building GRNs, many of the works did not consider the time-delay effect. Here, the authors propose a time-delay DDS model composed of linear difference equations to represent temporal interactions among significantly expressed genes. The authors also introduce interpolation scheme and re-sampling method for equalizing the non-uniformity of sampling time points. Statistical significance plays an active role in obtaining the optimal interaction matrix of GRNs. The constructed genetic network using linear multiple regression matches with the original data very well. Simulation results are given to demonstrate the effectiveness of the proposed method and model.