ddradseqtools: a software package for in silico simulation and testing of double-digest RADseq experiments.

Double-digested RADseq (ddRADseq) is a NGS methodology that generates reads from thousands of loci targeted by restriction enzyme cut sites, across multiple individuals. To be statistically sound and economically optimal, a ddRADseq experiment has a preliminary design stage that needs to consider issues related to the selection of enzymes, particular features of the genome of the focal species, possible modifications to the library construction protocol, coverage needed to minimize missing data, and the potential sources of error that may impact upon the coverage. We present ddradseqtools, a software package to help ddRADseq experimental design by (i) the generation of in silico double-digested fragments; (ii) the construction of modified ddRADseq libraries using adapters with either one or two indexes and degenerate base regions (DBRs) to quantify PCR duplicates; and (iii) the initial steps of the bioinformatics preprocessing of reads. ddradseqtools generates single-end (SE) or paired-end (PE) reads that may bear SNPs and/or indels. The effect of allele dropout and PCR duplicates on coverage is also simulated. The resulting output files can be submitted to pipelines of alignment and variant calling, to allow the fine-tuning of parameters. The software was validated with specific tests for the correct operability of the program. The correspondence between in silico settings and parameters from ddRADseq in vitro experiments was assessed to provide guidelines for the reliable performance of the software. ddradseqtools is cost-efficient in terms of execution time, and can be run on computers with standard CPU and RAM configuration.

Signatures of volcanism and aridity in the evolution of an insular pine (Pinus canariensis Chr. Sm. Ex DC in Buch).

Oceanic islands of volcanic origin provide useful templates for the study of evolution because they are subjected to recurrent perturbations that generate steep environmental gradients that may promote adaptation. Here we combine population genetic data from nuclear genes with the analysis of environmental variation and phenotypic data from common gardens to disentangle the confounding effects of demography and selection to identify the factors of importance for the evolution of the insular pine P. canariensis. Eight nuclear genes were partially sequenced in a survey covering the entire species range, and phenotypic traits were measured in four common gardens from contrasting environments. The explanatory power of population substrate age and environmental indices were assessed against molecular and phenotypic diversity estimates. In addition, neutral genetic variability (FST) and the genetic differentiation of phenotypic variation (QST) were compared in order to identify the evolutionary forces acting on these traits. Two key factors in the evolution of the species were identified: (1) recurrent volcanic activity has left an imprint in the genetic diversity of the nuclear genes; (2) aridity in southern slopes promotes local adaptation in the driest localities of P. canariensis, despite high levels of gene flow among populations.

Assessment of spatial discordance of primary and effective seed dispersal of European beech (Fagus sylvatica L.) by ecological and genetic methods.

Spatial discordance between primary and effective dispersal in plant populations indicates that postdispersal processes erase the seed rain signal in recruitment patterns. Five different models were used to test the spatial concordance of the primary and effective dispersal patterns in a European beech (Fagus sylvatica) population from central Spain. An ecological method was based on classical inverse modelling (SSS), using the number of seed/seedlings as input data. Genetic models were based on direct kernel fitting of mother-to-offspring distances estimated by a parentage analysis or were spatially explicit models based on the genotype frequencies of offspring (competing sources model and Moran-Clark's Model). A fully integrated mixed model was based on inverse modelling, but used the number of genotypes as input data (gene shadow model). The potential sources of error and limitations of each seed dispersal estimation method are discussed. The mean dispersal distances for seeds and saplings estimated with these five methods were higher than those obtained by previous estimations for European beech forests. All the methods show strong discordance between primary and effective dispersal kernel parameters, and for dispersal directionality. While seed rain was released mostly under the canopy, saplings were established far from mother trees. This discordant pattern may be the result of the action of secondary dispersal by animals or density-dependent effects; that is, the Janzen-Connell effect.

The Atlantic-Mediterranean watershed, river basins and glacial history shape the genetic structure of Iberian poplars.

Recent phylogeographic studies have elucidated the effects of Pleistocene glaciations and of Pre-Pleistocene events on populations from glacial refuge areas. This study investigates those effects in riparian trees (Populus spp.), whose particular features may convey enhanced resistance to climate fluctuations. We analysed the phylogeographic structure of 44 white (Populus alba), 13 black (Populus nigra) and two grey (Populus x canescens) poplar populations in the Iberian Peninsula using plastid DNA microsatellites and sequences. We also assessed fine-scale spatial genetic structure and the extent of clonality in four white and one grey poplar populations using nuclear microsatellites and we determined quantitative genetic differentiation (Q(ST) ) for growth traits in white poplar. Black poplar displayed higher regional diversity and lower differentiation than white poplar, reflecting its higher cold-tolerance. The dependence of white poplar on phreatic water was evidenced by strong differentiation between the Atlantic and Mediterranean drainage basins and among river basins, and by weaker isolation by distance within than among river basins. Our results suggest confinement to the lower river courses during glacial periods and moderate interglacial gene exchange along coastlines. In northern Iberian river basins, white poplar had lower diversity, fewer private haplotypes and larger clonal assemblies than in southern basins, indicating a stronger effect of glaciations in the north. Despite strong genetic structure and frequent asexual propagation in white poplar, some growth traits displayed adaptive divergence between drainage and river basins (Q(ST) >F(ST)), highlighting the remarkable capacity of riparian tree populations to adapt to regional environmental conditions.