The genome of the zebra mussel, Dreissena polymorpha: a resource for comparative genomics, invasion genetics, and biocontrol.

The zebra mussel, Dreissena polymorpha, continues to spread from its native range in Eurasia to Europe and North America, causing billions of dollars in damage and dramatically altering invaded aquatic ecosystems. Despite these impacts, there are few genomic resources for Dreissena or related bivalves. Although the D. polymorpha genome is highly repetitive, we have used a combination of long-read sequencing and Hi-C-based scaffolding to generate a high-quality chromosome-scale genome assembly. Through comparative analysis and transcriptomics experiments, we have gained insights into processes that likely control the invasive success of zebra mussels, including shell formation, synthesis of byssal threads, and thermal tolerance. We identified multiple intact steamer-like elements, a retrotransposon that has been linked to transmissible cancer in marine clams. We also found that D. polymorpha have an unusual 67 kb mitochondrial genome containing numerous tandem repeats, making it the largest observed in Eumetazoa. Together these findings create a rich resource for invasive species research and control efforts.

Spatial heterogeneity in landscape structure influences dispersal and genetic structure: empirical evidence from a grasshopper in an agricultural landscape.

Dispersal may be strongly influenced by landscape and habitat characteristics that could either enhance or restrict movements of organisms. Therefore, spatial heterogeneity in landscape structure could influence gene flow and the spatial structure of populations. In the past decades, agricultural intensification has led to the reduction in grassland surfaces, their fragmentation and intensification. As these changes are not homogeneously distributed in landscapes, they have resulted in spatial heterogeneity with generally less intensified hedged farmland areas remaining alongside streams and rivers. In this study, we assessed spatial pattern of abundance and population genetic structure of a flightless grasshopper species, Pezotettix giornae, based on the surveys of 363 grasslands in a 430-km² agricultural landscape of western France. Data were analysed using geostatistics and landscape genetics based on microsatellites markers and computer simulations. Results suggested that small-scale intense dispersal allows this species to survive in intensive agricultural landscapes. A complex spatial genetic structure related to landscape and habitat characteristics was also detected. Two P. giornae genetic clusters bisected by a linear hedged farmland were inferred from clustering analyses. This linear hedged farmland was characterized by high hedgerow and grassland density as well as higher grassland temporal stability that were suspected to slow down dispersal. Computer simulations demonstrated that a linear-shaped landscape feature limiting dispersal could be detected as a barrier to gene flow and generate the observed genetic pattern. This study illustrates the relevance of using computer simulations to test hypotheses in landscape genetics studies.

No clear effect of admixture between two European invading outbreaks of Diabrotica virgifera virgifera in natura.

In this study, we challenged the hypothesis that admixture may have had a positive impact in the context of the European invasion of the western corn rootworm (WCR), Diabrotica virgifera virgifera, LeConte. This beetle was introduced in Europe from the USA several times since the 1980's. The multiple introductions of this major pest of cultivated corn led to the formation of two major outbreaks in North Western (NW) Italy and in Central and South Eastern (CSE) Europe that eventually merged into a secondary contact zone where insects from both outbreaks interbreed. We collected about 600 insects from this contact zone and genotyped them using 13 microsatellite markers. Three types of information were obtained from the collected individuals: (i) their survival under starvation; (ii) their admixed status, determined through a Bayesian method of genetic clustering and (iii) their mating probability, studied via the detection, isolation and genotyping of sperm in female spermathecae. Twenty six % and 12% of the individuals were assigned to the NW Italy or the CSE Europe parental types, respectively, and 23% and 39% to the F1 and backcross hybrid types, respectively. Globally, our results do not reveal any significant impact of the admixed status on the mating probability and on the choice of mating partners. However the admixed status had a sex- and sampling site-dependent effect on survival in adults under starvation. In addition sex had an effect on survival, with mortality hazard about 3 times larger in males than in females. The consequences of these findings for the evolution of the admixture zone of northern Italy are discussed.

Sequence variability of the MspI satellite DNA family of the pinewood nematode Bursaphelenchus xylophilus at different geographic scales.

Tandemly repeated sequences known as satellite DNA (satDNA) generally exhibit complex evolutionary patterns of concerted evolution in which mutations are homogenized and fixed in a stochastic process of molecular drive. Here, the nucleotidic variability of the MspI satDNA family of the pinewood nematode Bursaphelenchus xylophilus is analyzed in order to understand the evolutionary dynamics of satDNA at the intraspecific level. A total of 425 MspI monomer units, either PCR-amplified from isolates of local (Peninsula of Setúbal, Portugal) or worldwide origin, or retrieved from the B. xylophilus genome sequence, were characterized and compared. Whatever their origin, sliding window analysis of sequence variability patterns among monomers revealed low, moderate and highly variant domains, indicating that variable levels of evolutionary constraint may act upon the entire monomers. The phylogenetic inference based on the different sets of MspI satDNA family for this species shows a broad polymorphism of the individual monomers, which were distributed into four main clusters. However, such clustering appeared independent from the geographic origin of the nematodes, and could not discriminate isolates or groups of geographically close isolates. Rather, the formation of different phylogenetic groups within this satDNA family suggests an a priori embodying of a set of diverging repeats from a common ancestor satDNA library, which have been differently amplified along the evolutionary pathway of this species. The present work improves knowledge on the evolutionary dynamics of satDNA at the intraspecific level, and provides new information on satDNA sequence variability among natural populations sampled at a local geographic scale.

Estimation of the dispersal of a major pest of maize by cline analysis of a temporary contact zone between two invasive outbreaks.

Dispersal is a key factor in invasion and in the persistence and evolution of species. Despite the importance of estimates of dispersal distance, dispersal measurement remains a real methodological challenge. In this study, we characterized dispersal by exploiting a specific case of biological invasion, in which multiple introductions in disconnected areas lead to secondary contact between two differentiated expanding outbreaks. By applying cline theory to this ecological setting, we estimated σ, the standard deviation of the parent-offspring distance distribution, of the western corn rootworm, Diabrotica virgifera virgifera, one of the most destructive pests of maize. This species is currently invading Europe, and the two largest invasive outbreaks, in northern Italy and Central Europe, have recently formed a secondary contact zone in northern Italy. We identified vanishing clines at 12 microsatellite loci throughout the contact zone. By analysing both the rate of change of cline slope and the spatial variation of linkage disequilibrium at these markers, we obtained two σ estimates of about 20 km/generation(1/2). Simulations indicated that these estimates were robust to changes in dispersal kernels and differences in population density between the two outbreaks, despite a systematic weak bias. These estimates are consistent with the results of direct methods for measuring dispersal applied to the same species. We conclude that secondary contact resulting from multiple introductions is very useful for the inference of dispersal parameters and should be more widely used in other species.

First insights into the genetic diversity of the pinewood nematode in its native area using new polymorphic microsatellite loci.

The pinewood nematode, Bursaphelenchus xylophilus, native to North America, is the causative agent of pine wilt disease and among the most important invasive forest pests in the East-Asian countries, such as Japan and China. Since 1999, it has been found in Europe in the Iberian Peninsula, where it also causes significant damage. In a previous study, 94 pairs of microsatellite primers have been identified in silico in the pinewood nematode genome. In the present study, specific PCR amplifications and polymorphism tests to validate these loci were performed and 17 microsatellite loci that were suitable for routine analysis of B. xylophilus genetic diversity were selected. The polymorphism of these markers was evaluated on nematodes from four field origins and one laboratory collection strain, all originate from the native area. The number of alleles and the expected heterozygosity varied between 2 and 11 and between 0.039 and 0.777, respectively. First insights into the population genetic structure of B. xylophilus were obtained using clustering and multivariate methods on the genotypes obtained from the field samples. The results showed that the pinewood nematode genetic diversity is spatially structured at the scale of the pine tree and probably at larger scales. The role of dispersal by the insect vector versus human activities in shaping this structure is discussed.