BackgroundThe variation of venom composition with geography is an important aspect of intraspecific variability in the Vipera genus, although causes of this variability remain unclear. The diversity of snake venom is important both for our understanding of venomous snake evolution and for the preparation of relevant antivenoms to treat envenomations. A geographic intraspecific variation in snake venom composition was recently reported for Vipera aspis aspis venom in France. Since 1992, cases of human envenomation after Vipera aspis aspis bites in south-east France involving unexpected neurological signs were regularly reported. The presence of genes encoding PLA2 neurotoxins in the Vaa snake genome led us to investigate any neurological symptom associated with snake bites in other regions of France and in neighboring countries. In parallel, we used several approaches to characterize the venom PLA2 composition of the snakes captured in the same areas.Methodology/Principal FindingsWe conducted an epidemiological survey of snake bites in various regions of France. In parallel, we carried out the analysis of the genes and the transcripts encoding venom PLA2s. We used SELDI technology to study the diversity of PLA2 in various venom samples. Neurological signs (mainly cranial nerve disturbances) were reported after snake bites in three regions of France: Languedoc-Roussillon, Midi-Pyrénées and Provence-Alpes-Côte d'Azur. Genomes of Vipera aspis snakes from south-east France were shown to contain ammodytoxin isoforms never described in the genome of Vipera aspis from other French regions. Surprisingly, transcripts encoding venom neurotoxic PLA2s were found in snakes of Massif Central region. Accordingly, SELDI analysis of PLA2 venom composition confirmed the existence of population of neurotoxic Vipera aspis snakes in the west part of the Massif Central mountains.Conclusions/SignificanceThe association of epidemiological studies to genetic, biochemical and immunochemical analyses of snake venoms allowed a good evaluation of the potential neurotoxicity of snake bites. A correlation was found between the expression of neurological symptoms in humans and the intensity of the cross-reaction of venoms with anti-ammodytoxin antibodies, which is correlated with the level of neurotoxin (vaspin and/or ammodytoxin) expression in the venom. The origin of the two recently identified neurotoxic snake populations is discussed according to venom PLA2 genome and transcriptome data.