I see that setting --quadrature-rule canonical in FFC enforces Gauss--Legendre--Jacobi rule.
Is there a way of using Gauss-Lobatto quadrature rule for the integration of the ufl forms?
Yes, although one may possibly get into troubles when quadrature points on element boundaries are used with floating-point arithmetics.