Water ice and gas hydrates can coexist in the permafrost and polar regions on Earth and in the universe. However, the role of ice in the mechanical response of ice-contained methane hydrates is still unclear. Here, we conduct direct million-atom molecular simulations of ice-contained polycrystalline methane hydrates and identify a crossover in the tensile strength and average compressive flow stress due to the presence of ice. The average mechanical shear strengths of hydrate-hydrate bicrystals are about three times as large as those of hydrate-ice bicrystals. The ice content, especially below 70%, shows a significant effect on the mechanical strengths of the polycrystals, which is mainly governed by the proportions of the hydrate-hydrate grain boundaries (HHGBs), the hydrate-ice grain boundaries (HIGBs), and the ice-ice grain boundaries (IIGBs). Quantitative analysis of the microstructure of the water cages in the polycrystals reveals the dissociation and reformation of various water cages due to mechanical deformation. These findings provide molecular insights into the mechanical behavior and microscopic deformation mechanisms of ice-contained methane hydrate systems on Earth and in the universe.@en