Find research datasets worth reusing

opencc-zeroDec 2017View details →

Data from: Genetic dissection of hybrid male sterility across stages of spermatogenesis

Hybrid sterility is a common form of reproductive isolation between nascent species. Although hybrid sterility is routinely documented and genetically dissected in speciation studies, its developmental basis is rarely examined, especially in generations beyond the F1. To identify phenotypic and genetic determinants of hybrid male sterility from a developmental perspective, we characterized testis histology in 312 F2 hybrids generated by intercrossing inbred strains of Mus musculus domesticus and M. m. musculus, two subspecies of house mice. Hybrids display a range of histologic abnormalities that indicate defective spermatogenesis. Among these abnormalities, we quantified decreased testis size, reductions in spermatocyte and spermatid number, increased apoptosis of meiosis I spermatocytes, and more multinucleated syncytia. Collectively, our phenotypic data point to defects in meiosis I as a primary barrier to reproduction. We identified seven quantitative trait loci (QTL) controlling five histologic traits. A region of chromosome 17 that contains Prdm9, a gene known to confer F1 hybrid male sterility, affects multinucleated syncytia and round spermatids, potentially extending the phenotypic outcomes of this incompatibility. The X chromosome also plays a key role, with loci affecting multinucleated syncytia, apoptosis of round spermatids, and round spermatid numbers. We detected an epistatic interaction between QTL on chromosomes 17 and X for multinucleated syncytia. Our results refine the developmental basis of a key reproductive barrier in a classic model system for speciation genetics.

Data from: Genome of the pitcher plant Cephalotus reveals genetic changes associated with carnivory

Carnivorous plants exploit animals as a nutritional source and have inspired long-standing questions about the origin and evolution of carnivory-related traits. To investigate the molecular bases of carnivory, we sequenced the genome of the heterophyllous pitcher plant Cephalotus follicularis, in which we succeeded in regulating the developmental switch between carnivorous and non-carnivorous leaves. Transcriptome comparison of the two leaf types and gene repertoire analysis identified genetic changes associated with prey attraction, capture, digestion and nutrient absorption. Analysis of digestive fluid proteins from C. follicularis and three other carnivorous plants with independent carnivorous origins revealed repeated co-options of stress-responsive protein lineages coupled with convergent amino acid substitutions to acquire digestive physiology. These results imply constraints on the available routes to evolve plant carnivory.

Data from: Metapopulation patterns of additive and nonadditive genetic variance in the sea bass (Dicentrarchus labrax)

Describing and explaining the geographic within-species variation in phenotypes ("phenogeography") in the sea over a species distribution range is central to our understanding of a variety of eco-evolutionary topics. However, phenogeographic studies that have a large potential to investigate adaptive variation are overcome by phylogeographic studies, still mainly focusing on neutral markers. How genotypic and phenotypic data could covary over large geographic scales remains poorly understood in marine species. We crossed 75 noninbred sires (five origins) and 26 dams (two origins; each side of a hybrid zone) in a factorial diallel cross in order to investigate geographic variation for early survival and sex ratio in the metapopulation of the European sea bass (Dicentrarchus labrax), a highly prized marine fish species. Full-sib families (N = 1,950) were produced and reared in a common environment. Parentage assignment of 7,200 individuals was performed with seven microsatellite markers. Generalized linear models showed significant additive effects for both traits and pleiotropy between traits. A significant nonadditive genetic effect was detected. Different expression of traits and distinct relative performances were found for reciprocal crosses involving populations located on each side of the main hybrid zone located at the Almeria-Oran front, illustrating asymmetric reproductive isolation. The poor fitness performance observed for the Western Mediterranean population of sea bass is discussed as it represents the main source of seed hatchery production, but also because it potentially illustrates nonadaptive introgression and maladaptation.

Data from: Evidence of genetic change in the flowering phenology of sea beets along a latitudinal cline within two decades

Sea beets grown from seeds collected in 1989 and 2009 along the coasts of France and adjacent regions were compared for flowering date under controlled conditions. Seeds from both collection years were sown simultaneously and cultivated under the same glasshouse conditions. Date of flowering onset and year of first flowering were recorded. There was an overall northward shift in flowering time of about 0.35° latitude (i.e. 39 km) over the 20 year period. The southern portion of the latitudinal gradient — i.e. from 44.7°N to 47.28°N — flowered significantly later by a mean of 1.78 days, equivalent to a 43.2 km northward shift of phenotypes. In the northern latitudes between 48.6°N and 52°N, flowering date was significantly earlier by a mean of 4.04 days, corresponding to a mean northward shift of 104.9 km, and this shift was apparently due to a diminished requirement of exposure to cold temperatures (i.e. vernalization), for which we found direct and indirect evidence. As all plants were grown from seed under identical conditions, we conclude that genetic changes occurred in the sensitivity to environmental cues that mediate the onset of flowering in both the northern and the southern latitudes of the gradient. Microevolution and gene flow may have contributed to this change. There was no significant change in the frequency of plants that flowered without vernalization. The lack of vernalization requirement may be associated with environmental instability rather than with climate conditions.

opencc-zeroDec 2012View details →

Data from: Parallel tagged next-generation sequencing on pooled samples – a new approach for population genetics in ecology and conservation

Next-generation sequencing (NGS) on pooled samples has already been broadly applied in human medical diagnostics and plant and animal breeding. However, thus far it has been only sparingly employed in ecology and conservation, where it may serve as a useful diagnostic tool for rapid assessment of species genetic diversity and structure at the population level. Here we undertake a comprehensive evaluation of the accuracy, practicality and limitations of parallel tagged amplicon NGS on pooled population samples for estimating species population diversity and structure. We obtained 16S and Cyt b data from 20 populations of Leiopelma hochstetteri, a frog species of conservation concern in New Zealand, using two approaches – parallel tagged NGS on pooled population samples and individual Sanger sequenced samples. Data from each approach were then used to estimate two standard population genetic parameters, nucleotide diversity (π) and population differentiation (FST), that enable population genetic inference in a species conservation context. We found a positive correlation between our two approaches for population genetic estimates, showing that the pooled population NGS approach is a reliable, rapid and appropriate method for population genetic inference in an ecological and conservation context. Our experimental design also allowed us to identify both the strengths and weaknesses of the pooled population NGS approach and outline some guidelines and suggestions that might be considered when planning future projects.

Data from: Creating new evolutionary pathways through bio-invasion: the population genetics of brushtail possums in New Zealand

Rapid increases in global trade and human movement have created novel mixtures of organisms bringing with them the potential to rapidly accelerate the evolution of new forms. The common brushtail possum (Trichosurus vulpecula), introduced into New Zealand from Australia in the 19th Century, is one such species having been sourced from multiple populations in its native range. Here, we combine microsatellite DNA and GIS-based spatial data to show that T.vulpecula originating from at least two different Australian locations exhibit a population structure that is commensurate with their introduction history and which cannot be explained by landscape features alone. Most importantly, we identify a hybrid zone between the two subspecies which appears to function as a barrier to dispersal. When combined with previous genetic, morphological and captive studies, our data suggest that assortative mating between the two subspecies may operate at a behavioural or species recognition level rather than through fertilization, genetic incompatibility, or developmental inhibition. Nevertheless, hybridization between the two subspecies of possum clearly occurs, creating the opportunity for novel genetic combinations that would not occur in their natural ranges and which is especially likely given that multiple contact zones occur in New Zealand. This discovery has implications for wildlife management in New Zealand because multiple contact zones are likely to influence the dispersal patterns of possums and because differential susceptibility to baiting with sodium fluoroacetate between possums of different origins, may promote novel genetic forms.

opencc-zeroDec 2017View details →

Data from: The hidden side of a major marine biogeographic boundary: a wide mosaic hybrid zone at the Atlantic–Mediterranean divide reveals the complex interaction between natural and genetic barriers in mussels

The Almeria-Oran Front (AOF) is a recognised hotspot of genetic differentiation in the sea. It is a barrier to dispersal and an ecological boundary, which explain the position of genetic breaks. However, the maintenance of genetic differentiation is likely reinforced by genetic barriers. A general drawback of previous studies is an insufficient density of sampling sites at the transition zone with a conspicuous lack of samples from the southern coastline. We analysed the genetic structure in the mussel Mytilus galloprovincialis with ancestry-informative loci. We discovered a 600 km wide mosaic hybrid zone eastward of the AOF along the Algerian coasts. This mosaic zone provides a new twist to our understanding of the Atlantic-Mediterranean transition because it demonstrates the two lineages can live in sympatry but hardly interbreed. This implies some form of reproductive isolation must exist to maintain the two genetic backgrounds locally cohesive. The zone ends with an abrupt genetic shift at a barrier to dispersal in the Gulf of Bejaia. Simulations in models that account for hydrodynamic features of the region support the hypothesis that sister hybrid zones could have been differentially trapped at two alternative barriers to dispersal or environmental boundaries. A preponderantly unidirectional north-south gene flow next to the AOF can also maintain a patch of an intrinsically maintained genetic background in the south and the mosaic structure. Our results concur with the coupling hypothesis that suggests natural barriers mostly explain the position of genetic breaks while their maintenance must additionally require genetic barriers.

opencc-zeroDec 2015View details →

Data from: Genetic variation in social environment construction influences the development of aggressive behavior in Drosophila melanogaster

Individuals are not merely subject to their social environments; they choose and create them, through a process called social environment (or social niche) construction. When genotypes differ in social environment-constructing behaviors, different genotypes are expected to experience different social environments. As social experience often affects behavioral development, quantitative genetics and psychology theories predict that genetic variation in social environment construction should have an important role in determining phenotypic variation; however, this hypothesis has not been tested directly. I identify multiple mechanisms of social environment construction that differ among natural genotypes of Drosophila melanogaster and investigate their consequences for the development of aggressive behavior. Male genotypes differed in the group sizes that they preferred and in their aggressive behavior; both of these behaviors influenced social experience, demonstrating that these behaviors function as social environment-constructing traits. Further, the effects of social experience—as determined in part by social environment construction—carried over to affect focal male aggression at a later time and with a new opponent. These results provide manipulative experimental support for longstanding hypotheses in psychology, that genetic variation in social environment construction has a causal role in behavioral development. More broadly, these results imply that studies of the genetic basis of complex traits should be expanded to include mechanisms by which genetic variation shapes the environments that individuals experience.

opencc-zeroDec 2010View details →

Data from: Genetic variation and evolution of secondary compounds in native and introduced populations of the invasive plant Melaleuca quinquenervia

We examined multivariate evolution of 20 leaf terpenoids in the invasive plant Melaleuca quinquenervia in a common garden experiment. While most compounds, including 1,8-Cineole and Viridiflorol, were reduced in home compared with invaded range genotypes, consistent with an evolutionary decrease in defense, one compound (E-Nerolidol) was greater in invaded than home range genotypes. Nerolidol was negatively genetically correlated with Cineole and Viridiflorol, and the increase in this compound in the new range may have been driven by this negative correlation. There was positive selection on all three focal compounds, and a loss of genetic variation in introduced range genotypes. Selection skewers analysis predicted an increase in Cineole and Viridiflorol and a decrease or no change in Nerolidol, in direct contrast to the observed changes in the new range. This discrepancy could be due to differences in patterns of selection, genetic correlations or the herbivore communities in the home versus introduced ranges. While evolutionary changes in most compounds were consistent with the Evolution of Increased Competitive Ability Hypothesis, changes in other compounds as well as selection patterns were not, indicating that it is important to understand selection and the nature of genetic correlations to predict evolutionary change in invasive species.

Data from: Do genetic drift and accumulation of deleterious mutations preclude adaptation? Empirical investigation using RADseq in a northern lacustrine fish

Understanding genomic signatures of divergent selection underlying long-term adaptation in populations located in heterogeneous environments is a key goal in evolutionary biology. In this study, we investigated neutral, adaptive and deleterious genetic variation using 7,192 SNPs in 31 Lake Trout (Salvelinus namaycush) populations (n = 673) from Québec, Canada. Average genetic diversity was low, weakly shared among lakes, and positively correlated to lake size, indicating a major role for genetic drift subsequent to lake isolation. Putatively deleterious mutations were on average at lower frequencies than the other SNPs, and their abundance relative to the entire polymorphism in each population was positively correlated to inbreeding, suggesting that the effectiveness of purifying selection was negatively correlated to inbreeding, as predicted from theory. Despite evidence for pronounced genetic drift and inbreeding, several outlier loci were associated with temperature and found in or close to genes with biologically relevant functions notably related to heat-stress and immune responses. Outcomes of gene-temperature associations were influenced by the inclusion of the most inbred populations, in which allele frequencies deviated the most from model predictions. This result illustrates challenge in identifying gene-environment associations in cases of high genetic drift and restricted gene flow and suggests limited adaptation in populations experiencing higher inbreeding. We discuss the relevance of these findings for the conservation and management, notably regarding stocking and genetic rescue, of Lake Trout populations and other species inhabiting highly fragmented habitats.

opencc-zeroDec 2011View details →

Data from: Genome divergence and the genetic architecture of barriers to gene flow between Lycaeides idas and L. melissa

Genome divergence during speciation is a dynamic process that is affected by various factors, including the genetic architecture of barriers to gene flow. Herein we quantitatively describe aspects of the genetic architecture of two sets of traits, male genitalic morphology and oviposition preference, that putatively function as barriers to gene flow between the butterfly species Lycaeides idas and L. melissa. Our analyses are based on unmapped DNA sequence data and a recently developed Bayesian regression approach that includes variable selection and explicit parameters for the genetic architecture of traits. A modest number of nucleotide polymorphisms explained a small to large proportion of the variation in each trait, and average genetic variant effects were non-negligible. Several genetic regions were associated with variation in multiple traits or with trait variation within- and among-populations. In some instances genetic regions associated with trait variation also exhibited exceptional genetic differentiation between speices or exceptional introgression in hybrids. These results are consistent with the hypothesis that divergent selection on male genitalia has contributed to heterogeneous genetic differentiation, and that both sets of traits affect fitness in hybrids. Although these results are encouraging, we highlight several difficulties related to understanding the genetics of speciation.

opencc-zeroDec 2012View details →

Data from: Adaptive genetic markers discriminate migratory runs of Chinook salmon (Oncorhynchus tshawytscha) amid continued gene flow

Neutral genetic markers are routinely used to define distinct units within species that warrant discrete management. Human-induced changes to gene flow however may reduce the power of such an approach. We tested the efficiency of adaptive vs. neutral genetic markers in differentiating temporally divergent migratory runs of Chinook salmon (Oncorhynchus tshawytscha) amid high gene flow owing to artificial propagation and habitat alteration. We compared seven putative migration timing genes to ten microsatellite loci in delineating three migratory groups of Chinook in the Feather River, CA: offspring of fall run hatchery broodstock that returned as adults to freshwater in fall (fall run), spring run offspring that returned in spring (spring run), and fall run offspring that returned in spring (FRS). We found evidence for significant differentiation between the fall and federally-listed threatened spring groups based on divergence at three circadian clock genes (OtsClock1b, OmyFbxw11 and Omy1009UW), but not neutral markers. We thus demonstrate the importance of genetic marker choice in resolving complex life history types. These findings directly impact conservation management strategies and add to previous evidence from Pacific and Atlantic salmon indicating that circadian clock genes influence migration timing.

opencc-zeroDec 2014View details →

Data from: Do island plant populations really have lower genetic variation than mainland populations? Effects of selection and distribution range on genetic diversity estimates

Ecological and evolutionary studies largely assume that island populations display low levels of neutral genetic variation. However, this notion has only been formally tested in a few cases involving plant taxa, and the confounding effect of selection on genetic diversity (GD) estimates based on putatively neutral markers has typically been overlooked. Here, we generated nuclear microsatellite and plastid DNA sequence data in Periploca laevigata, a plant taxon with an island-mainland distribution area, to (i) investigate whether selection affects GD estimates of populations across contrasting habitats and (ii) test the long-standing idea that island populations have lower GD than their mainland counterparts. Plastid data showed that colonization of the Canary Islands promoted strong lineage divergence within P. laevigata, which was accompanied by selective sweeps at several nuclear microsatellite loci. Inclusion of loci affected by strong divergent selection produced a significant downward bias in the GD estimates of the mainland lineage, but such underestimates were substantial (>14%) only when more than one loci under selection were included in the computations. When loci affected by selection were removed, we did not find evidence that insular Periploca populations have less GD than their mainland counterparts. The analysis of data obtained from a comprehensive literature survey reinforced this result, as overall comparisons of GD estimates between island and mainland populations were not significant across plant taxa (N=66), with the only exception of island endemics with narrow distributions. This study suggests that identification and removal of markers potentially affected by selection should be routinely implemented in estimates of GD, particularly if different lineages are compared. Furthermore, it provides compelling evidence that the expectation of low GD cannot be generalized to island plant populations.

Data from: Temporally isolated lineages of pink salmon reveal unique signatures of selection on distinct pools of standing genetic variation

A species' genetic diversity bears the marks of evolutionary processes that have occurred throughout its history. However, robust detection of selection in wild populations is difficult and often impeded by lack of replicate tests. Here, we investigate selection in pink salmon (Oncorhynchus gorbuscha) using genome scans coupled with inference from a haploid-assisted linkage map. Pink salmon have a strict 2-year semelparous life history which has resulted in temporally isolated (allochronic) lineages that remain sympatric through sharing of spawning habitats in alternate years. The lineages differ in a range of adaptive traits, suggesting different genetic backgrounds. We used genotyping by sequencing of haploids to generate a high-density linkage map with 7035 loci and screened an existing panel of 8036 loci for signatures of selection. The linkage map enabled identification of novel genomic regions displaying signatures of parallel selection shared between lineages. Furthermore, 24 loci demonstrated divergent selection and differences in genetic diversity between lineages, suggesting that adaptation in the 2 lineages has arisen from different pools of standing genetic variation. Findings have implications for understanding asynchronous population abundances as well as predicting future ecosystem impacts from lineage-specific responses to climate change.

Data from: Extinctions, genetic erosion and conservation options for the black rhinoceros (Diceros bicornis)

The black rhinoceros is again on the verge of extinction due to unsustainable poaching in its native range. Despite a wide historic distribution, the black rhinoceros was traditionally thought of as depauperate in genetic variation, and with very little known about its evolutionary history. This knowledge gap has hampered conservation efforts because hunting has dramatically reduced the species' once continuous distribution, leaving five surviving gene pools of unknown genetic affinity. Here we examined the range-wide genetic structure of historic and modern populations using the largest and most geographically representative sample of black rhinoceroses ever assembled. Using both mitochondrial and nuclear datasets, we described a staggering loss of 69% of the species' mitochondrial genetic variation, including the most ancestral lineages that are now absent from modern populations. Genetically unique populations in countries such as Nigeria, Cameroon, Chad, Eritrea, Ethiopia, Somalia, Mozambique, Malawi and Angola no longer exist. We found that the historic range of the West African subspecies (D. b. longipes), declared extinct in 2011, extends into southern Kenya, where a handful of individuals survive in the Masai Mara. We also identify conservation units that will help maintain evolutionary potential. Our results suggest a complete re-evaluation of current conservation management paradigms for the black rhinoceros.

opencc-zeroDec 2015View details →

Data from: Climate-related adaptive genetic variation and population structure in natural stands of Norway spruce in the South-Eastern Alps

Forest trees dominate many Alpine landscapes that are currently exposed to changing climate. Norway spruce is one of the most important conifer species of the Italian Alps, and natural populations are found across steep environmental gradients with large differences in temperature and moisture availability. This study seeks to determine and quantify patterns of genetic diversity in natural populations toward understanding adaptive responses to changing climate. Across the Italian species range, 24 natural stands were sampled with a major focus on the Eastern Italian Alps. Sampled trees were genotyped for 384 selected single nucleotide polymorphisms (SNPs) from 285 genes. A wide array of potential candidate genes was tested for correlation with climatic parameters. To minimize false-positive association between genotype and climate, population structure was investigated. Pairwise F ST estimates between sampled populations ranged between 0.000 and 0.075, with the highest values involving the two disjoint populations, Valdieri, on the western Italian Alps, and Campolino, the most southern population on the Apennines. Despite considerable genetic admixture among populations, both Bayesian and multivariate approach identified four genetic clusters. Selection scans revealed five F ST outliers, and the environmental association analysis detected ten SNPs associated to one or more climatic variables. Overall, 13 potentially adaptive loci were identified, three of which have been reported in a previous study on the same species conducted on a broader geographical scale. In our study, precipitation, more than temperature, was often associated with genotype; therefore, it appears as the most important environmental variable associated with the high sensitivity of Norway spruce to soil water supply. These findings provide relevant information for understanding and quantifying climate change effects on this species and its ability to genetically adapt.

opencc-zeroSep 2019View details →

Data from: Genetic architecture constrains exploitation of siderophore cooperation in the bacterium Burkholderia cenocepacia.

Abstract. Explaining how cooperation can persist in the presence of cheaters, exploiting the cooperative acts, is a challenge for evolutionary biology. Microbial systems have proved extremely useful to test evolutionary theory and identify mechanisms maintaining cooperation. One of the most widely studied system is the secretion and sharing of iron-scavenging siderophores by Pseudomonas bacteria, with many insights gained from this system now being considered as hallmarks of bacterial cooperation. Here, we introduce siderophore secretion by the bacterium Burkholderia cenocepacia H111 as a novel parallel study system, and show that this system behaves differently. For ornibactin, the main siderophore of this species, we discovered a novel mechanism of how cheating can be prevented. Particularly, we found that secreted ornibactin cannot be exploited by ornibactin-defective mutants because ornibactin receptor and synthesis genes are co-expressed from the same operon, such that disruptive mutations in synthesis genes compromise receptor availability required for siderophore uptake and cheating. For pyochelin, the secondary siderophore of this species, we found that cheating was possible, but the relative success of cheaters was positive frequency-dependent, thus diametrically opposite to the Pseudomonas and other microbial systems. Altogether, our results highlight that expanding our repertoire of microbial study systems leads to new discoveries and suggest that there is an enormous diversity of social interactions out there in nature, and we might have only looked at the tip of the iceberg so far.

Data from: Phylogeography in continuous space: coupling species distribution models and circuit theory to assess the effect of contiguous migration at different climatic periods on genetic differentiation in Busseola fusca (Lepidoptera: Noctuidae)

Current population genetic models fail to cope with genetic differentiation for species with large, contiguous and heterogeneous distribution. We show that in such a case, genetic differentiation can be predicted at equilibrium by circuit theory, where conductance corresponds to abundance in species distribution models (SDM). Circuit-SDM approach was used for the phylogeographic study of the lepidopteran cereal stemborer Busseola fusca Füller (Noctuidae) across sub-Saharan Africa. Species abundance was surveyed across its distribution range. SDM models were optimized and selected by cross validation. Relationship between observed matrices of genetic differentiation between individuals, and matrices of resistance distance was assessed through Mantel tests and redundancy discriminant analyses (RDA). A total of 628 individuals from 130 localities in 17 countries were genotyped at 7 microsatellite loci. Six population clusters were found based on a Bayesian analysis. The eastern margin of Dahomey Gap between East and West Africa was the main factor of genetic differentiation. The SDM projections at present, last interglacial and last glacial maximum periods were used for estimation of circuit resistance between locations of genotyped individuals. For all periods of time, when using either all individuals or only East-African individuals, partial Mantel r and RDA analyses conditioning on geographic distance were found significant. Under future projections (year 2080), partial r and RDA significance were different. From this study, it is concluded that analytical solutions provided by circuit theory are useful for the evolutionary management of populations and for phylogeographic analysis when coalescence times are not accessible by approximate Bayesian simulations.

opencc-zeroDec 2015View details →

Data from: Sea-level driven glacial-age refugia and post-glacial mixing on subtropical coasts, a palaeohabitat and genetic study

Using a novel combination of palaeohabitat modelling and genetic mixture analyses, we identify and assess a sea-level-driven recolonization process following the Last Glacial Maximum (LGM). Our palaeohabitat modelling reveals dramatic changes in estuarine habitat distribution along the coast of California (USA) and Baja California (Mexico). At the LGM (approx. 20 kya), when sea level was approximately 130 m lower, the palaeo-shoreline was too steep for tidal estuarine habitat formation, eliminating this habitat type from regions where it is currently most abundant, and limiting such estuaries to a northern and a southern refugium separated by 1000 km. We assess the recolonization of estuaries formed during post-LGM sea-level rise through examination of refugium-associated alleles and approximate Bayesian computation in three species of estuarine fishes. Results reveal sourcing of modern populations from both refugia, which admix in the newly formed habitat between the refuges. We infer a dramatic peak in habitat area between 15 and 10 kya with subsequent decline. Overall, this approach revealed a previously undocumented dynamic and integrated relationship between sea-level change, coastal processes and population genetics. These results extend glacial refugial dynamics to unglaciated subtropical coasts and have significant implications for biotic response to predicted sea-level rise.

Data from: Comparative assessment of SSR and SNP markers for inferring the population genetic structure of the common fungus Armillaria cepistipes

During the last years, simple sequence repeats (SSRs, also known as microsatellites) and single-nucleotide polymorphisms (SNPs) have become the most popular molecular markers for describing neutral genetic variation in populations of a wide range of organisms. However, only a limited number of studies has focused on comparing the performance of these two types of markers for describing the underlying genetic structure of wild populations. Moreover, none of these studies targeted fungi, the group of organisms with one of the most complex reproductive strategies. We evaluated the utility of SSRs and SNPs for inferring the neutral genetic structure of Armillaria cepistipes (basidiomycetes) at different spatial scales. For that, 407 samples were collected across a small (150 km2) area in the Ukrainian Carpathians and a large (41 000 km2) area in the Swiss Alps. All isolates were analyzed at 17 SSR loci distributed throughout the whole genome and at 24 SNP loci located in different single-copy conserved genes. The two markers showed different patterns of structure within the two spatial scales studied. The multi-allelic SSR markers seemed to be best suited for detecting genetic structure in indigenous fungal populations at a rather small spatial scale (radius of ~50-100 km). The pattern observed at SNP markers rather reflected ancient divergence of distant (~1000 km) populations that in addition are separated by mountain ranges. Despite these differences, both marker types were suitable for detecting the weak genetic structure of the two A. cepistipes populations investigated.