A method for stacking sequence optimization and aeroelastic tailoring of forward-swept composite wings is presented. It exploits bend-twist coupling to mitigate aeroelastic divergence. The method proposed here is intended for estimating potential weight savings during the preliminary aircraft design stages. A structural beam model of the composite wingbox is derived from anisotropic shell theory and the governing aeroelastic equations are presented for a spanwise discretized forward swept wing. Optimization of the system to reduce wing mass is undertaken for sweep angles of -35 • to 0 • and Mach numbers from 0.7 to 0.9. A subset of lamination parameters (LPs) and the number of laminate plies in each pre-defined direction (restricted to {0 • ,±45 • , 90 • }) serve as design variables. A bi-level hybrid optimization approach is employed, making use of a genetic algorithm (GA) and a subsequent gradientbased optimizer. Constraints are implemented to match lift requirements and prevent aeroelastic divergence, excessive deformations, airfoil stalling and structural failure. A permutation GA is then used to match specific composite ply stacking sequences to the optimum design variables with a limited number of manufacturing constraints considered for demonstration purposes. The optimization results in positive bend-twist coupling and