SNP4OrphanSpecies: A bioinformatics pipeline to isolate molecular markers for studying genetic diversity of orphan species.

Background: For several decades, an increase in disease or pest emergences due to anthropogenic introduction or environmental changes has been recorded. This increase leads to serious threats to the genetic and species diversity of numerous ecosystems. Many of these events involve species with poor or no genomic resources (called here "orphan species"). This lack of resources is a serious limitation to our understanding of the origin of emergent populations, their ability to adapt to new environments and to predict future consequences to biodiversity. Analyses of genetic diversity are an efficient method to obtain this information rapidly, but require available polymorphic genetic markers. New information: We developed a generic bioinformatics pipeline to rapidly isolate such markers with the goal for the pipeline to be applied in studies of invasive taxa from different taxonomic groups, with a special focus on forest fungal pathogens and insect pests. This pipeline is based on: 1) an automated de novo genome assembly obtained from shotgun whole genome sequencing using paired-end Illumina technology; 2) the isolation of single-copy genes conserved in species related to the studied emergent organisms; 3) primer development for multiplexed short sequences obtained from these conserved genes. Previous studies have shown that intronic regions of these conserved genes generally contain several single nucleotide polymorphisms within species. The pipeline's functionality was evaluated with sequenced genomes of five invasive or expanding pathogen and pest species in Europe (Armillariaostoyae (Romagn.) Herink 1973, Bursaphelenchusxylophilus Steiner & Buhrer 1934, Sphaeropsissapinea (fr.) Dicko & B. Sutton 1980, Erysiphealphitoides (Griffon & Maubl.) U. Braun & S. Takam. 2000, Thaumetopoeapityocampa Denis & Schiffermüller, 1775). We successfully isolated several pools of one hundred short gene regions for each assembled genome, which can be amplified in multiplex. The bioinformatics pipeline is user-friendly and requires little computational resources. This easy-to-set-up and run method for genetic marker identification will be useful for numerous laboratories studying biological invasions, but with limited resources and expertise in bioinformatics.

Genetic basis of growth, spring phenology, and susceptibility to biotic stressors in maritime pine.

Forest ecosystems are increasingly challenged by extreme events, for example, drought, storms, pest attacks, and fungal pathogen outbreaks, causing severe ecological and economic losses. Understanding the genetic basis of adaptive traits in tree species is of key importance to preserve forest ecosystems, as genetic variation in a trait (i.e., heritability) determines its potential for human-mediated or evolutionary change. Maritime pine (Pinus pinaster Aiton), a conifer widely distributed in southwestern Europe and northwestern Africa, grows under contrasted environmental conditions promoting local adaptation. Genetic variation at adaptive phenotypes, including height, spring phenology, and susceptibility to two fungal pathogens (Diplodia sapinea and Armillaria ostoyae) and an insect pest (Thaumetopoea pityocampa), was assessed in a range-wide clonal common garden of maritime pine. Broad-sense heritability was significant for height (0.219), spring phenology (0.165-0.310), and pathogen susceptibility (necrosis length caused by D. sapinea, 0.152; and by A. ostoyae, 0.021, measured on inoculated, excised branches under controlled conditions), but not for pine processionary moth incidence in the common garden. The correlations of trait variation among populations revealed contrasting trends for pathogen susceptibility to D. sapinea and A. ostoyae with respect to height. Taller trees showed longer necrosis length caused by D. sapinea while shorter trees were more affected by A. ostoyae. Moreover, maritime pine populations from areas with high summer temperatures and frequent droughts were less susceptible to D. sapinea but more susceptible to A. ostoyae. Finally, an association study using 4227 genome-wide SNPs revealed several loci significantly associated with each trait (range of 3-26), including a possibly disease-induced translation initiation factor, eIF-5, associated with needle discoloration caused by D. sapinea. This study provides important insights to develop genetic conservation and breeding strategies integrating species responses to biotic stressors.

Fast sequence-based microsatellite genotyping development workflow.

Application of high-throughput sequencing technologies to microsatellite genotyping (SSRseq) has been shown to remove many of the limitations of electrophoresis-based methods and to refine inference of population genetic diversity and structure. We present here a streamlined SSRseq development workflow that includes microsatellite development, multiplexed marker amplification and sequencing, and automated bioinformatics data analysis. We illustrate its application to five groups of species across phyla (fungi, plant, insect and fish) with different levels of genomic resource availability. We found that relying on previously developed microsatellite assay is not optimal and leads to a resulting low number of reliable locus being genotyped. In contrast, de novo ad hoc primer designs gives highly multiplexed microsatellite assays that can be sequenced to produce high quality genotypes for 20-40 loci. We highlight critical upfront development factors to consider for effective SSRseq setup in a wide range of situations. Sequence analysis accounting for all linked polymorphisms along the sequence quickly generates a powerful multi-allelic haplotype-based genotypic dataset, calling to new theoretical and analytical frameworks to extract more information from multi-nucleotide polymorphism marker systems.

Whole-genome sequencing reveals recent and frequent genetic recombination between clonal lineages of Cryphonectria parasitica in western Europe.

Changes in the mode of reproduction are frequently observed in invasive fungal populations. The ascomycete Cryphonectria parasitica, which causes Chestnut Blight, was introduced to Europe from North America and Asia in the 20th century. Previous genotyping studies based on ten microsatellite markers have identified several clonal lineages which have spread throughout western Europe, suggesting that asexuality was the main reproductive mode of this species during colonization, although occasional sexual reproduction is not excluded. Based on the whole-genome sequences alignment of 46 C. parasitica isolates from France, North America and Asia, genealogy and population structure analyses mostly confirmed these lineages as clonal. However, one of these clonal lineages showed a signal of strong recombination, suggesting different strategies of reproduction in western Europe. Signatures of several recent recombination events within all the French clonal lineages studied here were also identified, indicating that gene flow is regular between these lineages. In addition, haplotype identification of seven French clonal lineages revealed that emergences of new clonal lineages during colonization were the result of hybridization between the main expanding clonal lineages and minor haplotypes non-sequenced in the present study. This whole-genome sequencing study underlines the importance of recombination events in the invasive success of these clonal populations, and suggests that sexual reproduction may be more frequent within and between the western European clonal lineages of C. parasitica than previously assumed using few genetic markers.

Genetic analysis reveals efficient sexual spore dispersal at a fine spatial scale in Armillaria ostoyae, the causal agent of root-rot disease in conifers.

Armillaria ostoyae (sometimes named Armillaria solidipes) is a fungal species causing root diseases in numerous coniferous forests of the northern hemisphere. The importance of sexual spores for the establishment of new disease centres remains unclear, particularly in the large maritime pine plantations of southwestern France. An analysis of the genetic diversity of a local fungal population distributed over 500 ha in this French forest showed genetic recombination between genotypes to be frequent, consistent with regular sexual reproduction within the population. The estimated spatial genetic structure displayed a significant pattern of isolation by distance, consistent with the dispersal of sexual spores mostly at the spatial scale studied. Using these genetic data, we inferred an effective density of reproductive individuals of 0.1-0.3 individuals/ha, and a second moment of parent-progeny dispersal distance of 130-800 m, compatible with the main models of fungal spore dispersal. These results contrast with those obtained for studies of A. ostoyae over larger spatial scales, suggesting that inferences about mean spore dispersal may be best performed at fine spatial scales (i.e. a few kilometres) for most fungal species.

What was old is new again: thermal adaptation within clonal lineages during range expansion in a fungal pathogen.

Range-expanding species are expected to gain an increasing importance in the context of global change. They provide a great opportunity to study contemporary evolutionary changes and to unravel the mechanisms of evolution. Cryphonectria parasitica, the causal agent of chestnut blight, originating from Asia, has been spread since the beginning of the 20th century into different continents. We took advantage of the C. parasitica recent emergence in northern France to study the changes in population genetic structure and in phenotypic traits along this colonization and climatic gradient. Four hundred twenty-seven C. parasitica isolates were sampled in 47 chestnut sites in northern France. The C. parasitica outbreak in the north was found to be due to the expansion of five dominant clonal groups from southern France and to the emergence of a few rare recombined genotypes. The evolutionary changes during C. parasitica range expansion were studied by analysing phenotypic changes in isolates from the same clonal lineage, with or without a geographic shift. Growth rates were assessed in vitro, at four temperatures. The northern isolates grew faster at 12 and 15 °C and more slowly at 28 and 32 °C than the southern isolates. These results strongly suggest local adaptation to low temperatures in C. parasitica, with a trade-off of slower growth at high temperatures. They also reflect the high evolutionary potential of C. parasitica along a colonization gradient and show that clonal evolution is not a limitation for the rapid thermal adaptation of this invasive fungal species.

Multiple introductions and recombination in Cryphonectria hypovirus 1: perspective for a sustainable biological control of chestnut blight.

Cryphonectria hypovirus 1 (CHV1) is a mycovirus which decreases the virulence of its fungal host Cryphonectria parasitica, the causal agent of chestnut blight recently introduced in Europe. The understanding of the evolutionary processes which have shaped CHV1 populations in Europe is required to develop a sustainable biocontrol strategy targeting chestnut blight and effective in European chestnut forests. To retrace the evolutionary history of CHV1, we analyzed sequences from two genomic regions on a collection of 55 CHV1 strains from France and northern Spain, two countries where multiple introductions of C. parasitica occurred. Several recombination events and variable selection pressures contributed to CHV1 evolution, agreeing with a non-clock-like diversification rate. These two mechanisms may be at the origin of CHV1 population diversity observed in western Europe. Considering the actual prevalence of CHV1 and its association with host genotypes, multiple introductions of CHV1 may have occurred in Europe, some of them directly from Asia and some of them through North America. Although some viral strains remained with low frequency in their introduction area, multiple infections might have allowed homologous recombination within parental sequences. Some of these recombinant lineages are associated with the spread of CHV1 in European regions.

Exploring the role of asexual multiplication in poplar rust epidemics: impact on diversity and genetic structure.

Fungal plant pathogens, especially rust fungi (Pucciniales), are well known for their complex life cycles, which include phases of sexual and asexual reproduction. The effect of asexual multiplication on population genetic diversity has been investigated in the poplar rust fungus Melampsora larici-populina using a nested hierarchical sampling scheme. Four hierarchical levels were considered: leaf, twig, tree and site. Both cultivated and wild poplar stands were sampled at two time points at the start and end of rust epidemics. A total of 641 fungal isolates was analysed using nine microsatellite markers. This study revealed that the genetic signature of asexual multiplication in the wild poplar stand was seen only at lower hierarchical levels (leaf and twig). Moreover, we observed an erosion of clonal structure through time, with an increase in both gene and genotypic diversity. New genotypes contributed to host infection over time, which demonstrates the importance of allo-infection in the epidemic process in this host-pathogen system. Compared with the wild stands, the nearly lack of detection of clonal structure in the cultivated stands reflects the higher infection level on cultivated poplars. More generally, this genetic analysis illustrates the utility of population genetics approach for elucidating the proportion of asexual reproduction in the multiplication of isolates during an epidemic, and for proper quantification of asexual dispersal in plant pathogens.

Absence of isolation by distance patterns at the regional scale in the fungal plant pathogen Leptosphaeria maculans.

Outcomes of host-pathogen coevolution are influenced by migration rates of the interacting species. Reduced gene flow with increasing spatial distance between populations leads to spatial genetic structure, as predicted by the isolation by distance (IBD) model. In wind-dispersed plant-pathogenic fungi, a significant spatial genetic structure is theoretically expected if local spore dispersal is more frequent than long-distance dispersal, but this remains to be documented by empirical data. For 29 populations of the oilseed rape fungus Leptosphaeria maculans sampled from two French regions, genetic structure was determined using eight minisatellite markers. Gene diversity (H = 0.62-0.70) and haplotypic richness (R = 0.96-1) were high in all populations. No linkage disequilibrium was detected between loci, suggesting the prevalence of panmictic sexual reproduction. Analysis of molecular variance showed that > 97% of genetic diversity was observed within populations. Genetic differentiation was low among populations (F(st) < 0.05). Although direct methods previously revealed short-distance dispersal for L. maculans, our findings of no correlation between genetic and geographic distances among populations illustrate that the IBD model does not account for dispersal of the fungus at the spatial scale we examined. These results indicate high gene flow among French populations of L. maculans, suggesting high dispersal rates and/or large effective population sizes, two characteristics giving the pathogen high evolutionary potential against the deployment of resistant oilseed rape cultivars.

Finding single copy genes out of sequenced genomes for multilocus phylogenetics in non-model fungi.

Historically, fungal multigene phylogenies have been reconstructed based on a small number of commonly used genes. The availability of complete fungal genomes has given rise to a new wave of model organisms that provide large number of genes potentially useful for building robust gene genealogies. Unfortunately, cross-utilization of these resources to study phylogenetic relationships in the vast majority of non-model fungi (i.e. "orphan" species) remains an unexamined question. To address this problem, we developed a method coupled with a program named "PHYLORPH" (PHYLogenetic markers for ORPHans). The method screens fungal genomic databases (107 fungal genomes fully sequenced) for single copy genes that might be easily transferable and well suited for studies at low taxonomic levels (for example, in species complexes) in non-model fungal species. To maximize the chance to target genes with informative regions, PHYLORPH displays a graphical evaluation system based on the estimation of nucleotide divergence relative to substitution type. The usefulness of this approach was tested by developing markers in four non-model groups of fungal pathogens. For each pathogen considered, 7 to 40% of the 10-15 best candidate genes proposed by PHYLORPH yielded sequencing success. Levels of polymorphism of these genes were compared with those obtained for some genes traditionally used to build fungal phylogenies (e.g. nuclear rDNA, ß-tubulin, γ-actin, Elongation factor EF-1α). These genes were ranked among the best-performing ones and resolved accurately taxa relationships in each of the four non-model groups of fungi considered. We envision that PHYLORPH will constitute a useful tool for obtaining new and accurate phylogenetic markers to resolve relationships between closely related non-model fungal species.

Genetic structure of the poplar rust fungus Melampsora larici-populina: evidence for isolation by distance in Europe and recent founder effects overseas.

Dispersal has a great impact on the genetic structure of populations, but remains difficult to estimate by direct measures. In particular, gradual and stochastic dispersal are often difficult to assess and to distinguish, although they have different evolutionary consequences. Plant pathogens, especially rust fungi, are suspected to display both dispersal modes, though on different spatial scales. In this study, we inferred dispersal capacities of the poplar rust fungus Melampsora larici-populina by examining the genetic diversity and structure of 13 populations from eight European and two overseas countries in the Northern hemisphere. M. larici-populina was sampled from both cultivated hybrid poplars and on the wild host, Populus nigra. The populations were analyzed with 11 microsatellite and 8 virulence markers. Although isolates displayed different virulence profiles according to the host plant, neutral markers revealed little population differentiation with respect to the type of host. This suggests an absence of reproductive isolation between populations sampled from cultivated and wild poplars. Conversely, studying the relationship between geographic and genetic structure allowed us to distinguish between isolation by distance (IBD) patterns and long distance dispersal (LDD) events. The European populations exhibited a significant IBD pattern, suggesting a regular and gradual dispersal of the pathogen over this spatial scale. Nonetheless, the genetic differentiation between these populations was low, suggesting an important gene flow on a continental scale. The two overseas populations from Iceland and Canada were shown to result from rare LDD events, and exhibited signatures of strong founder effects. Furthermore, the high genetic differentiation between both populations suggested that these two recent introductions were independent. This study illustrated how the proper use of population genetics methods can enable contrasted dispersal modes to be revealed.

Challenges of microsatellite isolation in fungi.

Although they represent powerful genetic markers in many fields of biology, microsatellites have been isolated in few fungal species. The aim of this study was to assess whether obtaining microsatellite markers with an acceptable level of polymorphism is generally harder from fungi than in other organisms. We therefore surveyed the number, nature and polymorphism level of published microsatellite markers in fungi from the literature and from our own data on seventeen fungal microsatellite-enriched libraries, and in five other phylogroups (angiosperms, insects, fishes, birds and mammals). Fungal microsatellites indeed appeared both harder to isolate and to exhibit lower polymorphism than in other organisms. This appeared to be due, at least in part, to genomic specificities, such as scarcity and shortness of fungal microsatellite loci. A correlation was observed between mean repeat number and mean allele number in the published fungal microsatellite loci. The cross-species transferability of fungal microsatellites also appeared lower than in other phylogroups. However, microsatellites have been useful in some fungal species. Thus, the considerable advantages of these markers make their development worthwhile, and this study provides some guidelines for their isolation.

Fine-scale genetic structure and gene dispersal inferences in 10 neotropical tree species.

The extent of gene dispersal is a fundamental factor of the population and evolutionary dynamics of tropical tree species, but directly monitoring seed and pollen movement is a difficult task. However, indirect estimates of historical gene dispersal can be obtained from the fine-scale spatial genetic structure of populations at drift-dispersal equilibrium. Using an approach that is based on the slope of the regression of pairwise kinship coefficients on spatial distance and estimates of the effective population density, we compare indirect gene dispersal estimates of sympatric populations of 10 tropical tree species. We re-analysed 26 data sets consisting of mapped allozyme, SSR (simple sequence repeat), RAPD (random amplified polymorphic DNA) or AFLP (amplified fragment length polymorphism) genotypes from two rainforest sites in French Guiana. Gene dispersal estimates were obtained for at least one marker in each species, although the estimation procedure failed under insufficient marker polymorphism, limited sample size, or inappropriate sampling area. Estimates generally suffered low precision and were affected by assumptions regarding the effective population density. Averaging estimates over data sets, the extent of gene dispersal ranged from 150 m to 1200 m according to species. Smaller gene dispersal estimates were obtained in species with heavy diaspores, which are presumably not well dispersed, and in populations with high local adult density. We suggest that limited seed dispersal could indirectly limit effective pollen dispersal by creating higher local tree densities, thereby increasing the positive correlation between pollen and seed dispersal distances. We discuss the potential and limitations of our indirect estimation procedure and suggest guidelines for future studies.

Genetic diversity of the Chestnut blight fungus Cryphonectria parasitica in four French populations assessed by microsatellite markers.

Microsatellites are powerful markers to infer population genetic parameters. Here, 13 microsatellite loci isolated from a genomic and a cDNA library of Cryphonectria parasitica were used to characterize the genetic diversity and structure of four French populations. Twelve of these loci were polymorphic within populations, and average gene diversity (H(e)) was estimated to be 0.35. There was a lower genetic diversity in a south-eastern population relative to three south-western populations. In these three populations, microsatellite genotypic diversity was higher than vegetative compatibility type diversity. A high genetic differentiation (G(ST) = 0.27) suggested a low gene flow and/or founder effects of French populations which are in agreement with low dispersal of spores and different introductions of this species in southern France. This study demonstrates the significance of these microsatellite loci to assess gene flow and reproductive system in this important pathogen.

Gene flow and fine-scale genetic structure in a wind-pollinated tree species, Quercus lobata (Fagaceaee).

California Valley oak (Quercus lobata), one of the state's most distinctive oak species, has experienced serious demographic attrition since the 19th century, due to human activities. Recent estimates of pollen dispersal suggest a small reproductive neighborhood. Whether small neighborhood size is a recent phenomenon, a consequence of reduced gene flow caused by demographic changes, or whether it has been historically restricted, remains unclear. To examine this question, we have characterized the spatial genetic structure of N = 191 Q. lobata individuals, spread over an area of 230 ha, using eight microsatellite loci. The observed autocorrelogram suggests an historical standard deviation of gene flow distance of about 350 m per generation, higher than contemporary pollen dispersal estimates. To determine whether our estimates were affected by strong prevailing winds from the west-northwest, we developed and utilized a novel anisotropic autocorrelation analysis. We detected no more than a hint of anisotropy, and we concluded that adult spatial structure is indicative of strong historical signature of "isolation by distance." This historical estimate provides a useful reference value against which to gauge the future gene flow consequences of ongoing anthropogenic disturbance.

Using genetic markers to estimate the pollen dispersal curve.

Pollen dispersal is a critical process that shapes genetic diversity in natural populations of plants. Estimating the pollen dispersal curve can provide insight into the evolutionary dynamics of populations and is essential background for making predictions about changes induced by perturbations. Specifically, we would like to know whether the dispersal curve is exponential, thin-tailed (decreasing faster than exponential), or fat-tailed (decreasing slower than the exponential). In the latter case, rare events of long-distance dispersal will be much more likely. Here we generalize the previously developed TWOGENER method, assuming that the pollen dispersal curve belongs to particular one- or two-parameter families of dispersal curves and estimating simultaneously the parameters of the dispersal curve and the effective density of reproducing individuals in the population. We tested this method on simulated data, using an exponential power distribution, under thin-tailed, exponential and fat-tailed conditions. We find that even if our estimates show some bias and large mean squared error (MSE), we are able to estimate correctly the general trend of the curve - thin-tailed or fat-tailed - and the effective density. Moreover, the mean distance of dispersal can be correctly estimated with low bias and MSE, even if another family of dispersal curve is used for the estimation. Finally, we consider three case studies based on forest tree species. We find that dispersal is fat-tailed in all cases, and that the effective density estimated by our model is below the measured density in two of the cases. This latter result may reflect the difficulty of estimating two parameters, or it may be a biological consequence of variance in reproductive success of males in the population. Both the simulated and empirical findings demonstrate the strong potential of TWOGENER for evaluating the shape of the dispersal curve and the effective density of the population (d(e)).

Tracking a genetic signal of extinction-recolonization events in a neotropical tree species: Vouacapoua americana Aublet in French Guiana.

Drier periods from the late Pleistocene and early Holocene have been hypothesized to have caused the disappearance of various rainforest species over large geographical areas in South America and restricted the extant populations to mesic sites. Subsequent improvement in climatic conditions has been associated with recolonization. Changes in population size associated with these extinction-recolonization events should have affected genetic diversity within species. However, these historical hypotheses and their genetic consequences have rarely been tested in South America. Here, we examine the diversity of the chloroplast and nuclear genomes in a Neotropical rainforest tree species, Vouacapoua americana (Leguminosae, Caesalpinioideae) in French Guiana. The chloroplast diversity was analyzed using a polymerase chain reaction-restriction fragment length polymorphism method (six pairs of primers) in 29 populations distributed over most of French Guiana, and a subset of 17 populations was also analyzed at nine polymorphic microsatellite loci. To determine whether this species has experienced extinction-recolonization, we sampled populations in areas supposedly not or only slightly affected by climatic changes, where the populations would not have experienced frequent extinction, and in areas that appear to have been recently recolonized. In the putatively recolonized areas, we found patches of several thousands of hectares homogeneous for chloroplast variation that can be interpreted as the effect of recolonization processes from several geographical origins. In addition, we observed that, for both chloroplast and nuclear genomes, the populations in newly recolonized areas exhibited a significantly smaller allelic richness than others. Controlling for geographic distance, we also detected a significant correlation between chloroplast and nuclear population differentiation. This result indicates a cytonuclear disequilibrium that can be interpreted as a historical signal of a genetic divergence between fragmented populations. In conclusion, the spatial genetic structure of contemporary V. americana populations shows evidence that this species has experienced large extinction-recolonization events, which were possibly caused by past climatic change.