It is an inherent uncertainty problem that the application of laminar flow technology to the wing of large passenger aircraft is affected by flight conditions. In order to seek a more robust natural laminar flow control effect, it is necessary to develop an effective optimization design method. Meanwhile, attention must be given to the impact of crossflow (CF) instability brought on by the sweep angle. This paper constructs a robust optimization design framework based on discrete adjoint methods and non-intrusive polynomial chaos. Transition prediction is implemented by coupled Reynolds-Averaged Navier-Stokes (RANS) and simplified eN method, which can consider both Tollmien-Schlichting (TS) wave and crossflow vortex instability. We have performed gradient enhancement processing on the general Polynomial Chaos Expansion (PCE), which is advantageous to reduce the computational cost of single uncertainty propagation. This processing takes advantage of the gradient information obtained by solving the coupled adjoint equations considering transition. The statistical moment gradient solution used for the robust optimization design also uses the derivatives of coupled adjoint equations. The framework is applied to the robust design of a 25° swept wing with infinite span in transonic flow. The uncertainty quantification and sensitivity analysis on the baseline wing shows that the uncertainty quantification method in this paper has high accuracy, and qualitatively reveals the factors that dominate in different flow field regions. By the robust optimization design, the mean and standard deviation of the drag coefficient can be reduced by 29% and 45%, respectively, and compared with the deterministic optimization design results, there is less possibility of forming shock waves under flight condition uncertainties. Robust optimization results illustrate the trade-off between the transition delay and the wave drag reduction.