A central goal of molecular studies on ancient lake faunas is to resolve the origin and phylogeny of their strikingly diverse endemic species flocks. Another equally intriguing goal is to understand the integrity of individual morphologically diagnosed species, which should help to perceive the nature and speed of the speciation process, and the true biological species diversity. In the uniquely diverse Lake Baikal amphipod crustaceans, molecular data from shallow-water species have often disclosed their cryptic subdivision into geographically segregated genetic lineages, but the evidence so far is mainly based on mitochondrial DNA. We now present a lake-wide parallel survey of both mitochondrial and multilocus nuclear genetic structuring in the common shoreline amphipod Eulimnogammarus verrucosus, known to comprise three deep, parapatric mtDNA lineages. Allele frequencies of seven nuclear allozyme loci divide the data into three main groups whose distributions exactly match the distributions of the main mitochondrial lineages S, W, and E and involve a further division of the W cluster into two subgroups. The inter-group differences involve one to four diagnostic loci and additional group-specific alleles. The transition zones are either abrupt (1 km), occur over a long segment of uninhabitable shoreline, or may be gradual with non-coincident clinal change at different loci. Mitochondrial variation is hierarchically structured, each main lineage further subdivided into 2–4 parapatric sublineages or phylogroups, and patterns of further local segregation are seen in some of them. Despite the recurring observations of cryptic diversity in Baikalian amphipods, the geographical subdivisions and clade depths do not match in different taxa, defying a common explanation for the diversification in environmental history.