We numerically investigate the polarization characteristics of photonic crystal fibers selectively filled with metal wires into cladding air holes, through a full-vector modal solver based on the finite-element method (FEM). Firstly, we investigate the fundamental coupling properties between the core guided light and surface plasmon polaritons (SPPs) excited on the surface of metal wire. Secondly, we show that we can obtain highly polarization-dependent transmission characteristics in PCFs by introducing several metal wires closely aligned into the cladding, and reveal the strongly polarization-dependent coupling properties between the core guided modes and the SPP supermodes, which consist of discrete SPP modes. Finally, we show the importance of arranging the metal wires close to each other for high polarization-dependence.