Symbolic regression problems can be solved using grammatical evolution (GE), an evolutionary computation (EC) method, to find a function that coincides satisfactorily with the given datasets. The evolutional approach of GE is based on the grammar learning paradigm, which can translate the genotype (binary digit) into the phenotype (terminals and non-terminals). Unlike traditional codons in a genotype, the fittest codons in phenotype represented by the Backus-Naur form (BNF) are difficult for next generation genes to inherit the traits of parents, accounting for crossover and mutation. For this issue, this article presents a proposal of an advanced improvement to GE using a two-dimensional gene (GE2DG). In contrast to multi-chromosomal GE (GEMC), our proposal not only encloses the two-dimensional gene-expression for symbolic regression, but also introduces one independent gene defined as a conditional statement to express a new BNF grammar of an if-then (-else) branch. In the experiments described herein, continuous/discontinuous non-branch functions and continuous/discontinuous branch functions, four testing patterns, are considered as numerical examples. Results show that GE2DG has better performance than the original GE or GEMC. Especially for the case of branch functions, GE with hybrid chromosome (GEHC), where GE2DG is incorporated with GEMC, has faster convergence in symbolic regression than other methods.