Nils Arrigo

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference

Restriction site‐associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single‐nucleotide polymorphisms. As an empirical example, we use a double‐digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high‐altitude mountains in Mexico.

Rarely successful polyploids and their legacy in plant genomes

Polyploidy, or whole genome duplication, is recognized as an important feature of eukaryotic genome evolution. Among eukaryotes, polyploidy has probably had the largest evolutionary impact on vascular plants where many contemporary species are of recent polyploid origin. Genomic analyses have uncovered evidence of at least one round of polyploidy in the ancestry of most plants, fueling speculation that genome duplications lead to increases in net diversity. In spite of the frequency of ancient polyploidy, recent analyses have found that recently formed polyploid species have higher extinction rates than their diploid relatives. These results suggest that despite leaving a substantial legacy in plant genomes, only rare polyploids survive over the long term and most are evolutionary dead-ends.

The complex history of the olive tree: from Late Quaternary diversification of Mediterranean lineages to primary domestication in the northern Levant

The location and timing of domestication of the olive tree, a key crop in Early Mediterranean societies, remain hotly debated. Here, we unravel the history of wild olives (oleasters), and then infer the primary origins of the domesticated olive. Phylogeography and Bayesian molecular dating analyses based on plastid genome profiling of 1263 oleasters and 534 cultivated genotypes reveal three main lineages of pre-Quaternary origin. Regional hotspots of plastid diversity, species distribution modelling and macrofossils support the existence of three long-term refugia; namely the Near East (including Cyprus), the Aegean area and the Strait of Gibraltar. These ancestral wild gene pools have provided the essential foundations for cultivated olive breeding. Comparison of the geographical pattern of plastid diversity between wild and cultivated olives indicates the cradle of first domestication in the northern Levant followed by dispersals across the Mediterranean basin in parallel with the expansion of civilizations and human exchanges in this part of the world.

Naturalised Vitis Rootstocks in Europe and Consequences to Native Wild Grapevine

The genus Vitis is represented by several coexisting species in Europe. Our study focuses on naturalised rootstocks that originate in viticulture. The consequences of their presence to the landscape and to native European species (Vitis vinifera ssp. silvestris) are evaluated. This study compares ecological traits (seven qualitative and quantitative descriptors) and the genetic diversity (10 SSR markers) of populations of naturalised rootstocks and native wild grapevines. 18 large naturalised rootstock populations were studied in the Rhône watershed. Wild European grapevines are present in four main habitats (screes, alluvial forests, hedges, and streamside hedges). In contrast, naturalised rootstock populations are mainly located in alluvial forests, but they clearly take advantage of alluvial system dynamics and connectivity at the landscape level. These latter populations appear to reproduce sexually, and show a higher genetic diversity than Vitis vinifera ssp. silvestris. The regrouping of naturalised rootstocks in interconnected populations tends to create active hybrid swarms of rootstocks. The rootstocks show characters of invasive plants. The spread of naturalised rootstocks in the environment, the acceleration of the decline of the European wild grapevine, and the propagation of genes of viticultural interest in natural populations are potential consequences that should be kept in mind when undertaking appropriate management measures.

Hybridization Capture Using RAD Probes (hyRAD), a New Tool for Performing Genomic Analyses on Collection Specimens.

In the recent years, many protocols aimed at reproducibly sequencing reduced-genome subsets in non-model organisms have been published. Among them, RAD-sequencing is one of the most widely used. It relies on digesting DNA with specific restriction enzymes and performing size selection on the resulting fragments. Despite its acknowledged utility, this method is of limited use with degraded DNA samples, such as those isolated from museum specimens, as these samples are less likely to harbor fragments long enough to comprise two restriction sites making possible ligation of the adapter sequences (in the case of double-digest RAD) or performing size selection of the resulting fragments (in the case of single-digest RAD). Here, we address these limitations by presenting a novel method called hybridization RAD (hyRAD). In this approach, biotinylated RAD fragments, covering a random fraction of the genome, are used as baits for capturing homologous fragments from genomic shotgun sequencing libraries. This simple and cost-effective approach allows sequencing of orthologous loci even from highly degraded DNA samples, opening new avenues of research in the field of museum genomics. Not relying on the restriction site presence, it improves among-sample loci coverage. In a trial study, hyRAD allowed us to obtain a large set of orthologous loci from fresh and museum samples from a non-model butterfly species, with a high proportion of single nucleotide polymorphisms present in all eight analyzed specimens, including 58-year-old museum samples. The utility of the method was further validated using 49 museum and fresh samples of a Palearctic grasshopper species for which the spatial genetic structure was previously assessed using mtDNA amplicons. The application of the method is eventually discussed in a wider context. As it does not rely on the restriction site presence, it is therefore not sensitive to among-sample loci polymorphisms in the restriction sites that usually causes loci dropout. This should enable the application of hyRAD to analyses at broader evolutionary scales.

Gene flow among wild and domesticated almond species: insights from chloroplast and nuclear markers

Hybridization has played a central role in the evolutionary history of domesticated plants. Notably, several breeding programs relying on gene introgression from the wild compartment have been performed in fruit tree species within the genus Prunus but few studies investigated spontaneous gene flow among wild and domesticated Prunus species. Consequently, a comprehensive understanding of genetic relationships and levels of gene flow between domesticated and wild Prunus species is needed. Combining nuclear and chloroplastic microsatellites, we investigated the gene flow and hybridization among two key almond tree species, the cultivated Prunus dulcis and one of the most widespread wild relative Prunus orientalis in the Fertile Crescent. We detected high genetic diversity levels in both species along with substantial and symmetric gene flow between the domesticated P. dulcis and the wild P. orientalis. These results were discussed in light of the cultivated species diversity, by outlining the frequent spontaneous genetic contributions of wild species to the domesticated compartment. In addition, crop‐to‐wild gene flow suggests that ad hoc transgene containment strategies would be required if genetically modified cultivars were introduced in the northwestern Mediterranean.

Genetic structure and evolution of Alpine polyploid complexes: Ranunculus kuepferi (Ranunculaceae) as a case study

The alpine white‐flowered buttercup, Ranunculus kuepferi Greuter & Burdet, is a polyploid complex with diploids endemic to the southwestern Alps and polyploids – which have been previously described as apomictic – widespread throughout European mountains. Due to the polymorphic status of both its ploidy level and its reproductive mode, R. kuepferi represents a key species for understanding the evolution of polyploid lineages in alpine habitats. To disentangle the phylogeography of this polyploid taxon, we used cpDNA sequences and AFLP (amplified fragment length polymorphism) markers in 33 populations of R. kuepferi representative of its ploidy level and distribution area. Polyploid individuals were shown to be the result of at least two polyploidization events that may have taken place in the southwestern Alps. From this region, one single main migration of tetraploids colonized the entire Alpine range, the Apennines and Corsica. Genetic recombination among tetraploids was also observed, revealing the facultative nature of the apomictic reproductive mode in R. kuepferi polyploids. Our study shows the contrasting role played by diploid lineages mostly restricted to persistent refugia and by tetraploids, whose dispersal abilities have permitted their range extension all over the previously glaciated Alpine area and throughout neighbouring mountain massifs.

Origin and expansion of the allotetraploid Aegilops geniculata, a wild relative of wheat

*This study reconstructs the phylogeography of Aegilops geniculata, an allotetraploid relative of wheat, to discuss the impact of past climate changes and recent human activities (e.g. the early expansion of agriculture) on the genetic diversity of ruderal plant species. *We combined chloroplast DNA (cpDNA) sequencing, analysed using statistical parsimony network, with nonhierarchical K-means clustering of amplified fragment length polymorphism (AFLP) genotyping, to unravel patterns of genetic structure across the native range of Ae. geniculata. The AFLP dataset was further explored by measurement of the regional genetic diversity and the detection of isolation by distance patterns. *Both cpDNA and AFLP suggest an eastern Mediterranean origin of Ae. geniculata. Two lineages have spread independently over northern and southern Mediterranean areas. Northern populations show low genetic diversity but strong phylogeographical structure among the main peninsulas, indicating a major influence of glacial cycles. By contrast, low genetic structuring and a high genetic diversity are detected in southern Mediterranean populations. Finally, we highlight human-mediated dispersal resulting in substantial introgression between resident and migrant populations. *We have shown that the evolutionary trajectories of ruderal plants can be similar to those of wild species, but are interfered by human activities, promoting range expansions through increased long-distance dispersal and the creation of suitable habitats.

Methods for studying polyploid diversification and the dead end hypothesis: a reply to Soltis et al. (2014)

The fate of polyploid lineages has been of long-standing interest to evolutionary biologists. In our previous work (Mayrose et al., 2011a; reviewed in Arrigo & Barker, 2012), we used likelihood methods to estimate the effects of recent polyploidization events on diversification rates (speciation minus extinction). Our goal was to ask whether, across groups, there is a preponderance of evidence for or against the hypothesis that polyploid species diversify at different rates than diploid species. Our results demonstrated that polyploid lineages, when compared with their diploid relatives, tend to form new species more slowly and go extinct more quickly, for a combined effect of lower diversification rates. Soltis et al. (2014) recently called into question our conclusions and raised concerns about our analyses. Some of these concerns reflect open questions and limitations in the available data, but others are based on misinterpretations of the methods used and their implications. In the spirit of furthering the fascinating debate about the macroevolutionary consequences of polyploidization, we offer this response with the hope that it will clarify what has been shown, demonstrate the utility of the methods applied, and inspire future work. Soltis et al. (2014) frame their concerns as a series of arguments, and we follow their framework in responding to the points raised. Before embarking, it is worth reemphasizing what we defined as a ‘polyploid lineage’. Indeed, all angiosperms have one or more polyploidization events in their history if one looks back far enough (Jiao et al., 2011), and thus polyploidymust be definedwith respect to a reference time point. Because our focus was on the impacts on diversification of recent polyploidization events, we explicitly defined polyploids ‘as those lineages that underwent a polyploidization event since divergence from their generic ancestor’ (Mayrose et al., 2011a). This definition enabled us to perform a large-scale comparative study using chromosome number changes to infer polyploid transitions in 63 clades of plants. Thus, when asking whether polyploidization impacts subsequent diversification, it must be borne in mind that our analyses were conducted on rather short evolutionary timescales, comparing those taxa that had undergone fairly recent polyploidization (‘neopolyploids’) with those that had not. As we concluded in Mayrose et al. (2011a), it remains an interesting open question whether earlier polyploidy events had different impacts on speciation and extinction, as has been suggested, for example, in the context of mass extinction events (Fawcett et al., 2009; Vanneste et al., 2014).

Gene flow between wheat and wild relatives: empirical evidence from Aegilops geniculata, Ae. neglecta and Ae. triuncialis.

Gene flow between domesticated species and their wild relatives is receiving growing attention. This study addressed introgression between wheat and natural populations of its wild relatives (Aegilops species). The sampling included 472 individuals, collected from 32 Mediterranean populations of three widespread Aegilops species (Aegilops geniculata, Ae. neglecta and Ae. triuncialis) and compared wheat field borders to areas isolated from agriculture. Individuals were characterized with amplified fragment length polymorphism fingerprinting, analysed through two computational approaches (i.e. Bayesian estimations of admixture and fuzzy clustering), and sequences marking wheat‐specific insertions of transposable elements. With this combined approach, we detected substantial gene flow between wheat and Aegilops species. Specifically, Ae. neglecta and Ae. triuncialis showed significantly more admixed individuals close to wheat fields than in locations isolated from agriculture. In contrast, little evidence of gene flow was found in Ae. geniculata. Our results indicated that reproductive barriers have been regularly bypassed during the long history of sympatry between wheat and Aegilops.