Unexpected differences in the population genetics of phasmavirids (Bunyavirales) from subarctic ponds

Little is known of the evolution of RNA viruses in aquatic systems. Here, we assess the genetic connectivity of two bunyaviruses (Kigluaik phantom orthophasmavirus or KIGV and Nome phantom orthophasmavirus or NOMV) with zooplanktonic hosts from subarctic ponds. We expected weak genetic structure among populations as the hosts (phantom midges) have a terrestrial winged dispersal stage. To test whether their respective viruses mirror this structure, we collected and analyzed population datasets from 21 subarctic freshwater ponds and obtained sequences from all four genes in the viral genomes. Prevalence averaged 66 per cent for 514 host specimens and was not significantly different between recently formed thaw ponds and glacial ponds. Unexpectedly, KIGV from older ponds showed pronounced haplotype divergence with little evidence of genetic connectivity. However, KIGV populations from recent thaw ponds appeared to be represented by a closely related haplotype group, perhaps indicating a genotypic dispersal bias. Unlike KIGV, NOMV had modest structure and diversity in recently formed thaw ponds. For each virus, we found elevated genetic diversity relative to the host, but similar population structures to the host. Our results suggest that non-random processes such as virus-host interactions, genotypic bias, and habitat effects differ among polar aquatic RNA viruses.