AbstractA clear insight into large-scale community structure of planktonic copepods is critical to understanding mechanisms controlling diversity and biogeography of marine taxa, owing to their high abundance, ubiquity, and sensitivity to environmental changes. Here, we applied a 28S metabarcoding approach to large-scale communities of epipelagic and mesopelagic copepods at 70 stations across the Pacific Ocean and three stations in the Arctic Ocean. Major patterns of community structure and diversity, influenced by water mass structures, agreed with results from previous morphology-based studies. However, large-scale metabarcoding approach could detected community changes even under stable environmental conditions, including changes in the north/south subtropical gyres and east/west areas within each subtropical gyre. There were strong effects of epipelagic environment on mesopelagic communities, and community subdivisions were observed in the environmentally-stable mesopelagic layer. In each sampling station, higher operational taxonomic unit (OTU) numbers and lower phylogenetic diversity were observed in the mesopelagic layer than in the epipelagic layer, indicating a recent rapid increase of species numbers in the mesopelagic layer. The phylogenetic analysis utilizing representative sequences of OTUs revealed trends of recent emergence of cold-water OTUs mainly distributed at high latitudes with low water temperatures. Conversely, high diversity of copepods at low latitudes was suggested to have been formed through long evolutionary history under high water temperature. The metabarcoding results suggest that evolutionary processes have strong impacts on current patterns of copepod diversity, and support the “out of the tropics” theory explaining latitudinal diversity gradients of copepods. Both diversity patterns in epipelagic and mesopelagic showed high correlations to sea surface temperature; thus, predicted global warming may have a significant impact on copepod diversity in both layers.Author SummaryMarine planktonic copepods are highly dominant and diverse, and revealing their community structure and diversity is important to understanding marine ecosystems. We used molecular-based metabarcoding to reveal a total of 205 copepod communities in the ‘sunlight’ or epipelagic layer (0– 200 m) and the ‘twilight’ or mesopelagic layer (200–500 m and 500–1,000 m), mainly in the Pacific Ocean (data for 70 stations), but also in the Arctic Ocean (data for three stations). Different copepod communities were found in each geographical region with different environmental conditions, including tropical, subtropical, transition, Kuroshio Current, California Current, subarctic and arctic areas. The metabarcoding method sensitively detected small changes of copepod community even in environmentally-stable subtropical ocean systems and the mesopelagic layer. A high diversity of copepods was detected at low latitudes, and copepod diversity was higher in the mesopelagic layer than in the epipelagic layer in each area. These diversity patterns were influenced by both evolutionary history and present environmental conditions. The copepod community in the mesopelagic layer was strongly influenced by environmental conditions in the epipelagic layer. Thus, predicted climate changes may affect marine ecosystems not only in the epipelagic layer but also in the mesopelagic layer.