Genome-Wide Effects of Postglacial Colonization in Arabidopsis Lyrata
Total Page:16
File Type:pdf, Size:1020Kb
Heredity (2008) 100, 47–58 & 2008 Nature Publishing Group All rights reserved 0018-067X/08 $30.00 www.nature.com/hdy ORIGINAL ARTICLE Genome-wide effects of postglacial colonization in Arabidopsis lyrata M-H Muller1,2, J Leppa¨la¨1 and O Savolainen1 1Department of Biology, University of Oulu, Oulu, Finland and 2UMR DIA-PC 1097, INRA Montpellier, Domaine de Melgueil, Mauguio, France The perennial outcrossing Arabidopsis lyrata is becoming a taken into account in studies aiming at elucidating the genetic plant model species for molecular ecology and evolution. basis of local adaptation. As shown by sequence data, most However, its evolutionary history, and especially the impact of the loci used in this study do not evolve like typical of the climatic oscillations of the Pleistocene on its genetic microsatellite loci and show variable levels of homoplasy: this diversity and population structure, is not well known. We mode of evolution makes these markers less suitable to analyzed the broad-scale population structure of the species investigate the between-continent divergence and more based on microsatellite variation at 22 loci. A wide sample in generally the worldwide evolution of the species. Finally, a Europe revealed that glaciations and postglacial colonization strong negative correlation was detected between levels of have caused high divergence and high variation in variability within-population diversity and indices of differentiation such between populations. Colonization from Central Europe to as FST. We discuss the causes of this correlation as well as Iceland and Scandinavia was associated with a strong the potential bias it induces on the quantification and decrease of genetic diversity from South to North. On the interpretation of population structure. other hand, the Russian population included in our data set Heredity (2008) 100, 47–58; doi:10.1038/sj.hdy.6801057; may originate from a different refugium probably located published online 17 October 2007 more to the East. These genome-wide patterns must be Keywords: Arabidopsis lyrata; microsatellites; postglacial history; population structure; FST Introduction ssp. lyrata (North America), A. lyrata ssp. petraea (Europe and Siberia) and A. lyrata ssp. kamchatica (from Eastern Arabidopsis lyrata is a close relative of Arabidopsis thaliana, Asia to Northwest America). Within Europe, ssp. petraea the model organism for plant molecular biology, and itself has a patchy distribution extending from Central Europe is becoming a model for ecology and evolution (Clauss and to Norway (Jalas and Suominen, 1994; Jonsell et al., 1995). Koch, 2006). It benefits from a wealth of molecular tools It occurs under a variety of climatic and ecological and information accumulated for its congener, and the conditions, but is most often cold-tolerant and grows in genomic resources for A. lyrata are also rapidly improving low competition habitats (Jonsell et al., 1995; Clauss and (Kuittinen et al., 2002, 2004; http://www.jgi.doe.gov/ Mitchell-Olds, 2006). Studying the large-scale population sequencing/why/CSP2006/AlyrataCrubella.html). In con- structure and history of this species is needed as a trast to A. thaliana, it is a perennial, outcrossing species prerequisite for the study of potentially selected varia- and forms stable populations (Clauss and Koch, 2006). tion (Wright and Gaut, 2005). It is also interesting As such, it is a more tractable organism for population because A. lyrata constitutes a good example of the genetics and allows addressing evolutionary questions evolution of a species complex affected by the climatic and testing the validity of theoretical models (Riihima¨ki oscillations of the Pleistocene (Abbott and Brochmann, et al., 2005). For instance, A. lyrata, together with other 2003). species of the genus Arabidopsis, is intensively used for Indeed, the Pleistocene has been a period of intense the study of self-incompatibility evolution (Mable et al., and recurrent climatic changes that have affected the 2005; Schierup et al., 2006) or the genetic basis of distribution and genetic diversity of many species, adaptation at different geographic scales (for example, leading to the evolution of differentiated genetic entities Ka¨rkka¨inen et al., 2004; Riihima¨ki et al., 2005, Kivima¨ki and in some cases to speciation (Hewitt, 2004). Cold et al. 2007). periods displaced the distribution ranges southwards Despite this interest, the evolutionary history of the and European refugia have been described for a variety species is not well known. Follwing three subspecies are of organisms in Spain, Italy and the Balkans (Taberlet distinguished (Al-Shehbaz and O’Kane, 2002): A. lyrata et al., 1998). During interglacial times, species were able to expand northward; in Europe, different routes of colonization have contributed to the present-day dis- Correspondence: Dr M-H Muller, UMR DIA-PC 1097, INRA, Domaine tribution of many different species (Taberlet et al., 1998). de Melgueil, Mauguio F34130, France. E-mail: [email protected] In addition to these North–South movements, the Received 28 March 2007; revised 12 June 2007; accepted 3 August evolution of species having a ‘circumboreal’ distribution 2007; published online 17 October 2007 such as A. lyrata, has been affected by East–West Postglacial history of Arabidopsis lyrata M-H Muller et al 48 movements, from Europe to America, across the Atlantic, For all populations except Bohemia, we analyzed field- as well as through Russia and Asia (Abbott and collected seeds sampled from different families. Bohemia Brochmann, 2003). In plants, the vast majority of studies seeds were a population sample propagated in the of glacial and postglacial history have concerned tree laboratory, kindly provided by Mark McNair. Reykjavik species (for example, Heuertz et al., 2004), and a few seeds were kindly provided by MH Schierup and other circumboreal species (Alsos et al., 2005). Mayodan seeds by CH Langley. Previous studies on neutral variation, based on We used 22 microsatellite loci among which 21 were isozymes, microsatellites or nucleotide polymorphism, previously described (Table 1). The primers used for have revealed a strong differentiation between A.l. ssp. SOC1 were 50 TGTCAAATGTATTCGAGCAAGA and petraea and A.l. ssp. lyrata, and—generally but not 50 TGTAAGAGCAAGCACAAGAGGA. Some of these always—a lower level of diversity in ssp. lyrata, loci are mapped (Kuittinen et al., 2004; J Leppa¨la, potentially due to population bottlenecks (Van Treuren unpublished). Except for ELF3 and SOC1, all loci have et al., 1997; Savolainen et al., 2000; Ramos-Onsins et al., been derived from A. thaliana. 2004; Wright et al., 2006). Within Europe, genetic diversity appears to be strongly influenced by popula- Laboratory methods tion history and demography (Ka¨rkka¨inen et al., 2004; DNA was extracted from about 200 mg of plant leaves Ramos-Onsins et al., 2004). A recent study (Clauss and using the FastDNA kit and the FastPrep instrument Mitchell-Olds, 2006) analyzed the genetic diversity of (Qbiogene Inc, Carlsbad, CA, USA). The amplification ssp. petraea in Central Europe and suggested that the reaction consisted of 10 ng DNA, 0.25 mM each primer, species may have persisted in this region during the last one of which was fluorescently labelled with TET, HEX glacial maximum. The relationships of Northern Euro- or FAM, 1 Â reaction buffer (Promega, Madison, WI, pean, especially Scandinavian populations, with this USA), 2.5 mM MgCl , 100 mM dNTP, 0.15 U Taq DNA potential refugium, the importance of founder effects 2 polymerase (Promega), in a total volume of 15 ml. The during recolonization and the level of differentiation amplification method was as follows: 94 1C for 3 min, 35 within Europe compared to the between continent cycles of 94 1C for 20 s, 50 1C for 30 s, and 72 1C for 10s, differentiation (that is, differentiation with the subspe- followed by a final extension for 45 min at 72 1C. cies ssp. lyrata) remain to be studied. Electrophoresis was performed on an ABI Prism 377 With their high level of diversity, ease of use and sequencer. The GENESCAN 3.1 and GENOTYPER 2.0 codominant nature, microsatellite markers have already (Applied Biosystems, Foster City, CA, USA) were used to proven useful in small-scale analyses of A. lyrata (Clauss analyze the DNA fragments and to score the genotypes. et al., 2002; Clauss and Mitchell-Olds, 2006). In the We reported all observed size classes. The final extension present study, we surveyed microsatellite variation at 22 time was long to ensure that one T was added to all loci in a sample including seven European populations fragments by the Taq polymerase. of ssp. petraea and one American population of ssp. As most of the microsatellite loci had been derived lyrata. We asked the following main questions: how has from A. thaliana, their molecular structure and thus their postglacial colonization impacted the patterns of within mode of evolution in A. lyrata could differ from that of population diversity, and divergence of populations? typical microsatellites. For 18 loci, we thus sequenced Have the populations reached a genetic equilibrium after one allele from each of the populations Karhuma¨ki, the colonization? What are the implications of the Plech, Spiterstulen and Mayodan. For homozygote demographic history for studies of molecular ecology genotypes, we directly sequenced the polymerase chain and functional genomics? reaction product. For heterozygotes, we cloned the We show that populations of A. lyrata exhibit highly polymerase chain reaction products into TOPO TA PCR variable levels of diversity and are strongly differentiated, 2.1 cloning vector (Invitrogen, Carlsbad, CA, USA) and especially within Europe, as a consequence of the different sequenced one single clone per individual. The sequen- routes of postglacial colonization. Further, we illustrate cing reactions were run in the ABI PRISM 377 sequencer how the variable within-population diversity can impact and the contigs were read and aligned using the Staden- the methods of analyzing population differentiation.