Genetic Diversity and Population Structure Analysis of Three Deep-Sea Amphipod Species from Geographically Isolated Hadal Trenches in the Pacific Ocean

Amphipods of the superfamily Lysianassoidea that inhabit the hadal zone ( > 6000 m) have large bathymetric ranges and play a key role in deep ocean ecosystems. The endemism of these amphipod species makes them a good model for investigating potent natural selection and restricted dispersal in deep ocean trenches. Here, we describe genetic diversity and intraspecific population differentiation among three amphipod species from four Pacific trenches based on a mtDNA concatenated dataset (CO Ι and 16S rRNA genes) from 150 amphipod individuals. All amphipod populations had low genetic diversity, as indicated by haplotype and nucleotide diversity values. Population geographic relationship analysis of two Alicella gigantea populations revealed no genetic differentiation between these two localities (pairwise genetic differentiation coefficient = 0.00032, gene flow = 784.58), and the major variation (99.97%) was derived from variation within the populations. Historical demographic events were investigated using Tajima’s D and Fu’s F neutrality tests and analysis of mismatch distribution. Consistent results provided strong evidence to support the premise that demographic expansion occurred only for the Mariana population of Hirondellea gigas, possibly within the last 2.1–3.4 million years. These findings suggest that the formation of amphipod population structure might be the result of multiple factors including high hydrostatic pressure, food distribution, trench topographic forcing and potential ecological interactions.