The primary objective of this study is to compare the impact of CT80, CTVT and IT on exercise tolerance measured after the 12-week intervention phase. Preliminary data analysis of the first 18 subjects (6 in each arm) showed an average improvement in endurance time to constant-load cycling (ES=0.32). Due to this interim analysis, α level for sample size calculation should be adjusted from an α spending function. Using Pocock, α spending function at t=0.1765 (18 from 102 subjects projected), the loss of α calculated is 0.01324 resulting in a new α of 0.03676. For a power of 0.80 and an α of 0.03676, we would need 36 patients/group for a total sample size of 108 to see a statistical difference across the groups. Due to expected 10% loss of subjects during follow-up, 40 subjects/groups (N=120) will be recruited. To examine between-group differences in exercise tolerance from week 0 to week 12, 2-way repeated-measures ANOVA will be performed, with treatment group as the between-subjects factor with 3 levels and time as the within-subjects factor with 2 levels. If a significant interaction between treatment group and time is detected, simple main effect tests will be conducted to identify the time point for which differences between treatment groups are significant. If a significant simple main effect is obtained, pairwise comparisons with TukeyB corrections will be conducted to identify between which treatments groups the differences occurred. Lost cases will be managed using the "bringing the last value forward" method. Analyses will be conducted with SPSS 19.0 (Chicago, IL). Between groups differences in exercise tolerance from week 12 to week 52 will be compared with 2-way repeated-measures ANOVA. To compare the 3 arms on functional, psychological, cognition and HRQL variables from week 0 to week 12 to week 52, mixed model repeated measure analysis will be applied. If a significant interaction between treatment group and time is detected, simple main effect tests will be conducted to identify the time point(s) for which differences between treatment groups are significant. Mean compliance to the 12-week program and mean adherence at week 52 will be compared between treatment arms with 1-way ANOVA, followed up, if needed, by pairwise comparisons to locate the differences. Patient-related covariates and potential co-interventions (pulmonary function, medications and smoking status, exacerbations, hospitalisations) will be compared between treatment groups. If significant group differences are found on these variables, ANOVAs (or Mixed models) will be replaced by analyses of covariance (ANCOVAs or Mixed models) with the variable(s) in question as the covariate(s).