Several gene-based association tests for time-to-event traits have been proposed recently, to detect whether a gene region (containing multiple variants), as a set, is associated with the survival outcome. However, for bivariate survival outcomes, to the best of our knowledge, there is no statistical method that can be directly applied for gene-based association analysis. Motivated by a genetic study to discover gene regions associated with the progression of a bilateral eye disease, Age-related Macular Degeneration (AMD), we implement a novel functional regression method under the copula framework. Specifically, the effects of variants within a gene region are modeled through a functional linear model, which then contributes to the marginal survival functions within the copula.Generalized score test and likelihood ratio test statistics are derived to test for the association between bivariate survival traits and the genetic region. Extensive simulation studies are conducted to evaluate the type-I error control and power performance of the proposed approach, with comparisons to several existing methods for a single survival trait, as well as the marginal Cox functional regression model using the robust sandwich estimator for bivariate survival traits. Finally, we apply our method to a large AMD study, the Age-related Eye Disease Study (AREDS), to identify gene regions that are associated with AMD progression. The method has been implemented and added into a newly developed R package {CopulaCenR}.