Fabienne Justy

Comparative allozyme and microsatellite population structure in a narrow endemic plant species, Centaurea corymbosa Pourret (Asteraceae)

Centaurea corymbosa Pourret (Asteraceae) is a narrow endemic species known only from six populations located in a 3‐km2 area in the south of France. Earlier field experiments have suggested that pollen and seed dispersal were highly restricted within and among populations. Consistent with the field results, populations were highly differentiated for five allozyme loci and among‐population variation fitted an isolation‐by‐distance model. In the present study, we investigated the genetic structure of C. corymbosa using six microsatellite loci. As with allozymes, microsatellites revealed no within‐population structure and a large differentiation among populations. However, allozyme loci were less powerful than microsatellites in detecting the extent of gene flow assessed by assignment tests. The patterns of structuration greatly varied among loci for both types of marker; we suggest that differences in single‐locus pattern could mainly be an effect of stochastic variation for allozymes and an effect of variation in mutation rate for microsatellites. In contrast to the multilocus results, the two most polymorphic microsatellite loci did not show any isolation‐by‐distance pattern. Our results suggest that highly variable loci might not always be the best suited markers to quantify levels of gene flow among populations.

THE INFLUENCE OF SELF-FERTILIZATION AND POPULATION DYNAMICS ON THE GENETIC STRUCTURE OF SUBDIVIDED POPULATIONS: A CASE STUDY USING MICROSATELLITE MARKERS IN THE FRESHWATER SNAIL BULINUS TRUNCATUS

The distribution of neutral genetic variability within and among sets of populations results from the combined actions of genetic drift, migration, extinction and recolonization processes, mutation, and the mating system. We here analyzed these factors in 38 populations of the hermaphroditic snail Bulinus truncatus. The sampling area covered a large part of the species range. The variability was analyzed using four polymorphic microsatellite loci. A very large number of alleles (up to 55) was found at the level of the whole study. Observed heterozygote deficiencies within populations are consistent with very high selfing rates, generally above 0.80, in all populations. These should depress the variability within populations, because of low effective size, genetic hitchhiking, and background selection, whatever the model of mutation assumed. However, that some populations exhibit much more variability than others suggests that historical demographic processes (e.g., population size variation, bottlenecks, or founding events) may play a significant role. A hierarchical analysis of the distribution of the variability across populations indicates a strong pattern of isolation by distance, whatever the geographical scale considered. Our analysis also illustrates how the mutation rate may affect population differentiation, as different mutation rates result in different levels of homoplasy at microsatellite loci. The effects of both genetic drift and gene flow vary with the temporal and spatial scales considered in B. truncatus populations.

Conservation genetics and population history of the threatened European mink Mustela lutreola , with an emphasis on the west European population

In species of great conservation concern, special attention must be paid to their phylogeography, in particular the origin of animals for captive breeding and reintroduction. The endangered European mink lives now in at least three well‐separated populations in northeast, southeast and west Europe. Our aim is to assess the genetic structure of these populations to identify ‘distinct population segments’ (DPS) and advise captive breeding programmes. First, the mtDNA control region was completely sequenced in 176 minks and 10 polecats. The analysis revealed that the western population is characterized by a single mtDNA haplotype that is closely related to those in eastern regions but nevertheless, not found there to date. The northeast European animals are much more variable (π = 0.012, h = 0.939), with the southeast samples intermediate (π = 0.0012, h = 0.469). Second, 155 European mink were genotyped using six microsatellites. The latter display the same trends of genetic diversity among regions as mtDNA [gene diversity and allelic richness highest in northeast Europe (HE = 0.539, RS = 3.76), lowest in west Europe (HE = 0.379, RS = 2.12)], and provide evidences that the southeast and possibly the west populations have undergone a recent bottleneck. Our results indicate that the western population derives from a few animals which recently colonized this region, possibly after a human introduction. Microsatellite data also reveal that isolation by distance occurs in the western population, causing some inbreeding because related individuals mate. As genetic data indicate that the three populations have not undergone independent evolutionary histories for long (no phylogeographical structure), they should not be considered as distinct DPS. In conclusion, the captive breeding programme should use animals from different parts of the species’ present distribution area.

Permanent genetic resources added to molecular ecology resources database 1 January 2009-30 April 2009.

This article documents the addition of 283 microsatellite marker loci to the Molecular Ecology Resources Database. Loci were developed for the following species: Agalinis acuta; Ambrosia artemisiifolia; Berula erecta; Casuarius casuarius; Cercospora zeae‐maydis; Chorthippus parallelus; Conyza canadensis; Cotesia sesamiae; Epinephelus acanthistius; Ficedula hypoleuca; Grindelia hirsutula; Guadua angustifolia; Leucadendron rubrum; Maritrema novaezealandensis; Meretrix meretrix; Nilaparvata lugens; Oxyeleotris marmoratus; Phoxinus neogaeus; Pristomyrmex punctatus; Pseudobagrus brevicorpus; Seiridium cardinale; Stenopsyche marmorata; Tetranychus evansi and Xerus inauris. These loci were cross‐tested on the following species: Agalinis decemloba; Agalinis tenella; Agalinis obtusifolia; Agalinis setacea; Agalinis skinneriana; Cercospora zeina; Cercospora kikuchii; Cercospora sorghi; Mycosphaerella graminicola; Setosphaeria turcica; Magnaporthe oryzae; Cotesia flavipes; Cotesia marginiventris; Grindelia Xpaludosa; Grindelia chiloensis; Grindelia fastigiata; Grindelia lanceolata; Grindelia squarrosa; Leucadendron coniferum; Leucadendron salicifolium; Leucadendron tinctum; Leucadendron meridianum; Laodelphax striatellus; Sogatella furcifera; Phoxinus eos; Phoxinus rigidus; Phoxinus brevispinosus; Phoxinus bicolor; Tetranychus urticae; Tetranychus turkestani; Tetranychus ludeni; Tetranychus neocaledonicus; Tetranychus amicus; Amphitetranychus viennensis; Eotetranychus rubiphilus; Eotetranychus tiliarium; Oligonychus perseae; Panonychus citri; Bryobia rubrioculus; Schizonobia bundi; Petrobia harti; Xerus princeps; Spermophilus tridecemlineatus and Sciurus carolinensis.

Next-generation sequencing and phylogenetic signal of complete mitochondrial genomes for resolving the evolutionary history of leaf-nosed bats (Phyllostomidae)

Leaf-nosed bats (Phyllostomidae) are one of the most studied groups within the order Chiroptera mainly because of their outstanding species richness and diversity in morphological and ecological traits. Rapid diversification and multiple homoplasies have made the phylogeny of the family difficult to solve using morphological characters. Molecular data have contributed to shed light on the evolutionary history of phyllostomid bats, yet several relationships remain unresolved at the intra-familial level. Complete mitochondrial genomes have proven useful to deal with this kind of situation in other groups of mammals by providing access to a large number of molecular characters. At present, there are only two mitogenomes available for phyllostomid bats hinting at the need for further exploration of the mitogenomic approach in this group. We used both standard Sanger sequencing of PCR products and next-generation sequencing (NGS) of shotgun genomic DNA to obtain new complete mitochondrial genomes from 10 species of phyllostomid bats, including representatives of major subfamilies, plus one outgroup belonging to the closely-related mormoopids. We then evaluated the contribution of mitogenomics to the resolution of the phylogeny of leaf-nosed bats and compared the results to those based on mitochondrial genes and the RAG2 and VWF nuclear makers. Our results demonstrate the advantages of the Illumina NGS approach to efficiently obtain mitogenomes of phyllostomid bats. The phylogenetic signal provided by entire mitogenomes is highly comparable to the one of a concatenation of individual mitochondrial and nuclear markers, and allows increasing both resolution and statistical support for several clades. This enhanced phylogenetic signal is the result of combining markers with heterogeneous evolutionary rates representing a large number of nucleotide sites. Our results illustrate the potential of the NGS mitogenomic approach for resolving the evolutionary history of phyllostomid bats based on a denser species sampling.

Ascidian Mitogenomics: Comparison of Evolutionary Rates in Closely Related Taxa Provides Evidence of Ongoing Speciation Events

Ascidians are a fascinating group of filter-feeding marine chordates characterized by rapid evolution of both sequences and structure of their nuclear and mitochondrial genomes. Moreover, they include several model organisms used to investigate complex biological processes in chordates. To study the evolutionary dynamics of ascidians at short phylogenetic distances, we sequenced 13 new mitogenomes and analyzed them, together with 15 other available mitogenomes, using a novel approach involving detailed whole-mitogenome comparisons of conspecific and congeneric pairs. The evolutionary rate was quite homogeneous at both intraspecific and congeneric level, and the lowest congeneric rates were found in cryptic (morphologically undistinguishable) and in morphologically very similar species pairs. Moreover, congeneric nonsynonymous rates (dN) were up to two orders of magnitude higher than in intraspecies pairs. Overall, a clear-cut gap sets apart conspecific from congeneric pairs. These evolutionary peculiarities allowed easily identifying an extraordinary intraspecific variability in the model ascidian Botryllus schlosseri, where most pairs show a dN value between that observed at intraspecies and congeneric level, yet consistently lower than that of the Ciona intestinalis cryptic species pair. These data suggest ongoing speciation events producing genetically distinct B. schlosseri entities. Remarkably, these ongoing speciation events were undetectable by the cox1 barcode fragment, demonstrating that, at low phylogenetic distances, the whole mitogenome has a higher resolving power than cox1. Our study shows that whole-mitogenome comparative analyses, performed on a suitable sample of congeneric and intraspecies pairs, may allow detecting not only cryptic species but also ongoing speciation events.

Genetic, ecological, behavioral and geographic differentiation of populations in a thistle weevil: implications for speciation and biocontrol.

Because weevils are used as biocontrol agents against thistles, it is important to document and understand host shifts and the evolution of host‐specificity in these insects. Furthermore, such host shifts are of fundamental interest to mechanisms of speciation. The mediterranean weevil Larinus cynarae normally parasitizes either one of two thistle genera, Onopordum and Cynara, being locally monophagous. In Sardinia, however, both host genera are used. We used three types of data to help understand this complex host use: (i) weevil attack rates on the two host genera among 53 different populations in Sardinia and nearby Corsica, (ii) host preference in a lab setting, and (iii) genetic (allozyme) differentiation among weevil populations exploiting the same or different hosts. Using a subset of populations from northern Sardinia, we attempted to relate interpopulation differences in host preference to gene flow among populations by comparing pairwise differences in oviposition preference (Qst) and in allozyme frequencies (Fst). Overall, Qst and Fst were positively correlated. Fst was positively correlated with geographic distance among pairs of populations using the same host, but not among different‐host population pairs. As mating occurs on the hosts, this result suggests reinforcement. Genetic evidence indicates Cynara as the ancestral host of the weevils from both islands and our current studies suggest repeated attempts to colonize Onopordum, with a successful shift in Corsica and a partial shift in Sardinia. This scenario would explain why in Sardinia the level of attack was higher on Cynara than on Onopordum and why, when given a choice in the laboratory, Sardinian weevils preferred Cynara even when sampled from Onopordum. The lability of host shifts in L. cynarae supports caution in using these or related weevils as biocontrol agents of exotic thistles.

Mitogenomic phylogeny, diversification, and biogeography of South American spiny rats

Echimyidae is one of the most speciose and ecologically diverse rodent families in the world, occupying a wide range of habitats in the Neotropics. However, a resolved phylogeny at the genus-level is still lacking for these 22 genera of South American spiny rats, including the coypu (Myocastorinae), and 5 genera of West Indian hutias (Capromyidae) relatives. Here, we used Illumina shotgun sequencing to assemble 38 new complete mitogenomes, establishing Echimyidae, and Capromyidae as the first major rodent families to be completely sequenced at the genus-level for their mitochondrial DNA. Combining mitogenomes and nuclear exons, we inferred a robust phylogenetic framework that reveals several newly supported nodes as well as the tempo of the higher level diversification of these rodents. Incorporating the full generic diversity of extant echimyids leads us to propose a new higher level classification of two subfamilies: Euryzygomatomyinae and Echimyinae. Of note, the enigmatic Carterodon displays fast-evolving mitochondrial and nuclear sequences, with a long branch that destabilizes the deepest divergences of the echimyid tree, thereby challenging the sister-group relationship between Capromyidae and Euryzygomatomyinae. Biogeographical analyses involving higher level taxa show that several vicariant and dispersal events impacted the evolutionary history of echimyids. The diversification history of Echimyidae seems to have been influenced by two major historical factors, namely (1) recurrent connections between Atlantic and Amazonian Forests and (2) the Northern uplift of the Andes.

A cost-effective straightforward protocol for shotgun Illumina libraries designed to assemble complete mitogenomes from non-model species

The mitogenome is an inescapable tool in conservation biology studies. Yet, its routine sequencing may remain tricky despite next-generation sequencing technologies. An enrichment step is often necessary but not always straightforward depending on the initial DNA quality or quantity. Furthermore, the availability of close mitochondrial DNA reference sequences for non-model species limits the primer design for long-range PCR or bait synthesis. Here we propose an easy and cost-effective protocol without enrichment step for building and sequencing multiplexed Illumina libraries from small quantities of either high-quality or degraded genomic DNA. We validated the approach through the successful assembly of the complete mitogenome of 60 bats and 7 tunicates. Our protocol allows the sequencing and assembly of mitochondrial genomes from non-model species with sufficient coverage for applications in conservation genetics.

Interaction of climate, demography and genetics: a ten-year study of Brassica insularis, a narrow endemic Mediterranean species

Long-term demographic surveys, needed to obtain accurate information on population dynamics and efficiently manage rare species, are still very scarce. Matrix population models are useful tools to identify key demographic transitions and thus help setting up conservation actions. Furthermore, the combination of ecological, demographic and genetic data is likely to improve the identification of the threats acting upon populations and help conservation decisions. In this paper we illustrate the power of this approach on Brassica insularis, a Mediterranean endemic plant species, rare and endangered in Corsica (France). In four populations of this species, a long-term demographic survey (2000–2009), genetic analyses (in 2000 and 2009) and survey of ecological variables (climatic variables, competition and herbivory) were performed. By using both deterministic and stochastic matrix model analyses, we assessed the viability of each population and tested for both spatial and temporal variations in demographic vital rates. Populations exhibited differing demographic behaviours and environmental stochasticity occurred in populations. Significant correlations between climatic variables and vital rates were detected. Stochastic simulations suggested that three out of the four populations studied might present a high risk of extinction on the short-term and should actively be managed, or at least surveyed. It could be, however, that two of these populations are experiencing density-dependent regulation, rather than being declining. Microsatellite diversity was slightly reduced in a single population and similar in the three others, consistently with expectations based on population census size and geographic area, as well as with diversity at the S-locus observed in 2000. The combination of all data led to specific recommendations for managing each population. We discuss the implications for conservation of such a general approach.