1
|
Warson J, Baguette M, Stevens VM, Honnay O, De Kort H. The impact of habitat loss on molecular signatures of coevolution between an iconic butterfly (Alcon blue) and its host plant (Marsh gentian). J Hered 2023; 114:22-34. [PMID: 36749638 DOI: 10.1093/jhered/esac059] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2022] [Accepted: 10/19/2022] [Indexed: 11/30/2022] Open
Abstract
Habitat loss is threatening natural communities worldwide. Small and isolated populations suffer from inbreeding and genetic drift, which jeopardize their long-term survival and adaptive capacities. However, the consequences of habitat loss for reciprocal coevolutionary interactions remain poorly studied. In this study, we investigated the effects of decreasing habitat patch size and connectivity associated with habitat loss on molecular signatures of coevolution in the Alcon blue butterfly (Phengaris alcon) and its most limited host, the marsh gentian (Gentiana pneumonanthe). Because reciprocal coevolution is characterized by negative frequency-dependent selection as a particular type of balancing selection, we investigated how signatures of balancing selection vary along a gradient of patch size and connectivity, using single nucleotide polymorphisms (SNPs). We found that signatures of coevolution were unaffected by patch characteristics in the host plants. On the other hand, more pronounced signatures of coevolution were observed in both spatially isolated and in large Alcon populations, together with pronounced spatial variation in SNPs that are putatively involved in coevolution. These findings suggest that habitat loss can facilitate coevolution in large butterfly populations through limiting swamping of locally beneficial alleles by maladaptive ones. We also found that allelic richness (Ar) of the coevolutionary SNPs is decoupled from neutral Ar in the butterfly, indicating that habitat loss has different effects on coevolutionary as compared with neutral processes. We conclude that this specialized coevolutionary system requires particular conservation interventions aiming at generating a spatial mosaic of both connected and of isolated habitat to maintain coevolutionary dynamics.
Collapse
Affiliation(s)
- Jonas Warson
- Plant Conservation and Population Biology, Department of Biology, University of Leuven, Heverlee, Belgium
- Leuven Plant Institute, Heverlee, Belgium
| | - Michel Baguette
- Centre National de la Recherche Scientifique, SETE Station d'Ecologie Théorique et Expérimentale, UMR 5321, Moulis, France
- Institut Systématique, Evolution, Biodiversité (ISYEB), UMR 7205 Museum National d'HistoireNaturelle, CNRS, Sorbonne Université, EPHE, Université des Antilles, Paris, France
| | - Virginie M Stevens
- Centre National de la Recherche Scientifique, SETE Station d'Ecologie Théorique et Expérimentale, UMR 5321, Moulis, France
| | - Olivier Honnay
- Plant Conservation and Population Biology, Department of Biology, University of Leuven, Heverlee, Belgium
- Leuven Plant Institute, Heverlee, Belgium
| | - Hanne De Kort
- Plant Conservation and Population Biology, Department of Biology, University of Leuven, Heverlee, Belgium
- Leuven Plant Institute, Heverlee, Belgium
| |
Collapse
|
2
|
MacDonald ZG, Snape KL, Roe AD, Sperling F. Host association, environment, and geography underlie genomic differentiation in a major forest pest. Evol Appl 2022; 15:1749-1765. [PMID: 36426133 PMCID: PMC9679251 DOI: 10.1111/eva.13466] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2022] [Accepted: 07/29/2022] [Indexed: 11/30/2022] Open
Abstract
Diverse geographic, environmental, and ecological factors affect gene flow and adaptive genomic variation within species. With recent advances in landscape ecological modelling and high-throughput DNA sequencing, it is now possible to effectively quantify and partition their relative contributions. Here, we use landscape genomics to identify determinants of genomic differentiation in the forest tent caterpillar, Malacosoma disstria, a widespread and irruptive pest of numerous deciduous tree species in North America. We collected larvae from multiple populations across Eastern Canada, where the species experiences a diversity of environmental gradients and feeds on a number of different host tree species, including trembling aspen (Populus tremuloides), sugar maple (Acer saccharum), red oak (Quercus rubra), and white birch (Betula papyrifera). Using a combination of reciprocal causal modelling (RCM) and distance-based redundancy analyses (dbRDA), we show that differentiation of thousands of genome-wide single nucleotide polymorphisms (SNPs) among individuals is best explained by a combination of isolation by distance, isolation by environment (spatial variation in summer temperatures and length of the growing season), and differences in host association. Configuration of suitable habitat inferred from ecological niche models was not significantly related to genomic differentiation, suggesting that M. disstria dispersal is agnostic with respect to habitat quality. Although population structure was not discretely related to host association, our modelling framework provides the first molecular evidence of host-associated differentiation in M. disstria, congruent with previous documentation of reduced growth and survival of larvae moved between natal host species. We conclude that ecologically mediated selection is contributing to variation within M. disstria, and that divergent adaptation related to both environmental conditions and host association should be considered in ongoing research and management of this important forest pest.
Collapse
Affiliation(s)
- Zachary G. MacDonald
- Department of Biological SciencesUniversity of AlbertaEdmontonAlbertaCanada
- UCLA La Kretz Center for California Conservation ScienceUniversity of California Los AngelesLos AngelesCaliforniaUSA
- Institute of the Environmental and SustainabilityUniversity of California Los AngelesLos AngelesCaliforniaUSA
| | - Kyle L. Snape
- Department of Biological SciencesUniversity of AlbertaEdmontonAlbertaCanada
| | - Amanda D. Roe
- Great Lakes Forestry Centre, Canadian Forest ServiceNatural Resources CanadaSault Ste. MarieOntarioCanada
| | | |
Collapse
|
3
|
Serranito B, Cavalazzi M, Vidal P, Taurisson-Mouret D, Ciani E, Bal M, Rouvellac E, Servin B, Moreno-Romieux C, Tosser-Klopp G, Hall SJG, Lenstra JA, Pompanon F, Benjelloun B, Da Silva A. Local adaptations of Mediterranean sheep and goats through an integrative approach. Sci Rep 2021; 11:21363. [PMID: 34725398 PMCID: PMC8560853 DOI: 10.1038/s41598-021-00682-z] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2021] [Accepted: 09/03/2021] [Indexed: 11/20/2022] Open
Abstract
Small ruminants are suited to a wide variety of habitats and thus represent promising study models for identifying genes underlying adaptations. Here, we considered local Mediterranean breeds of goats (n = 17) and sheep (n = 25) from Italy, France and Spain. Based on historical archives, we selected the breeds potentially most linked to a territory and defined their original cradle (i.e., the geographical area in which the breed has emerged), including transhumant pastoral areas. We then used the programs PCAdapt and LFMM to identify signatures of artificial and environmental selection. Considering cradles instead of current GPS coordinates resulted in a greater number of signatures identified by the LFMM analysis. The results, combined with a systematic literature review, revealed a set of genes with potentially key adaptive roles in relation to the gradient of aridity and altitude. Some of these genes have been previously implicated in lipid metabolism (SUCLG2, BMP2), hypoxia stress/lung function (BMPR2), seasonal patterns (SOX2, DPH6) or neuronal function (TRPC4, TRPC6). Selection signatures involving the PCDH9 and KLH1 genes, as well as NBEA/NBEAL1, were identified in both species and thus could play an important adaptive role.
Collapse
Affiliation(s)
- Bruno Serranito
- INRA, EA7500, USC1061 GAMAA, Univ. Limoges, 87000, Limoges, France
- CRESCO, Museum National d'Histoire Naturelle (MNHN), 35800, Dinard, France
| | | | - Pablo Vidal
- Universidad Catolica de Valencia, Valencia, Spain
| | - Dominique Taurisson-Mouret
- GEOLAB, UMR 6042, Univ. Limoges, Limoges, France
- CNRS, UMR 5815, Dynamiques du droit, Université de Montpellier, Montpellier, France
| | - Elena Ciani
- Department of Biosciences, Biotechnologies and Biopharmaceutics, University of Bari, Bari, Italy
| | - Marie Bal
- GEOLAB, UMR 6042, Univ. Limoges, Limoges, France
| | | | - Bertrand Servin
- GenPhySE, INRAE, ENVT, Université de Toulouse, 31326, Castanet-Tolosan, France
| | | | | | - Stephen J G Hall
- Estonian University of Life Sciences, Kreutzwaldi 5, 51014, Tartu, Estonia
| | - Johannes A Lenstra
- Faculty of Veterinary Medicine, Utrecht University, Yalelaan 104, 3584 CM, Utrecht, The Netherlands
| | - François Pompanon
- Univ. Grenoble Alpes, Univ. Savoie Mont Blanc, CNRS, LECA, F-38000, Grenoble, France
| | - Badr Benjelloun
- Univ. Grenoble Alpes, Univ. Savoie Mont Blanc, CNRS, LECA, F-38000, Grenoble, France
- National Institute of Agronomic Research (INRA), Regional Centre of Agronomic Research, Beni-Mellal, Morocco
| | - Anne Da Silva
- INRA, EA7500, USC1061 GAMAA, Univ. Limoges, 87000, Limoges, France.
| |
Collapse
|
4
|
Abstract
Understanding the genetic architecture of complex traits is a major objective in biology. The standard approach for doing so is genome-wide association studies (GWAS), which aim to identify genetic polymorphisms responsible for variation in traits of interest. In human genetics, consistency across studies is commonly used as an indicator of reliability. However, if traits are involved in adaptation to the local environment, we do not necessarily expect reproducibility. On the contrary, results may depend on where you sample, and sampling across a wide range of environments may decrease the power of GWAS because of increased genetic heterogeneity. In this study, we examine how sampling affects GWAS in the model plant species Arabidopsis thaliana. We show that traits like flowering time are indeed influenced by distinct genetic effects in local populations. Furthermore, using gene expression as a molecular phenotype, we show that some genes are globally affected by shared variants, whereas others are affected by variants specific to subpopulations. Remarkably, the former are essentially all cis-regulated, whereas the latter are predominately affected by trans-acting variants. Our result illustrate that conclusions about genetic architecture can be extremely sensitive to sampling and population structure.
Collapse
Affiliation(s)
| | - Stephan Reinert
- Center for Computational and Theoretical Biology, University of Würzburg, Würzburg, Germany
| | - Magnus Nordborg
- Gregor Mendel Institute, Austrian Academy of Sciences, Vienna BioCenter, Vienna, Austria
| | - Arthur Korte
- Center for Computational and Theoretical Biology, University of Würzburg, Würzburg, Germany
| |
Collapse
|
5
|
Helsen K, Acharya KP, Graae BJ, De Kort H, Brunet J, Chabrerie O, Cousins SAO, De Frenne P, Hermy M, Verheyen K, Pélabon C. Earlier onset of flowering and increased reproductive allocation of an annual invasive plant in the north of its novel range. ANNALS OF BOTANY 2020; 126:1005-1016. [PMID: 32582950 PMCID: PMC7596373 DOI: 10.1093/aob/mcaa110] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/25/2019] [Accepted: 06/17/2020] [Indexed: 06/11/2023]
Abstract
BACKGROUND AND AIMS It remains unclear whether invasive species can maintain both high biomass and reproductive output across their invaded range. Along latitudinal gradients, allocation theory predicts that faster flowering onset at high latitudes results in maturation at smaller size and thus reduced reproductive output. For annual invasive plants, more favourable environmental conditions at low latitudes probably result in stronger competition of co-occurring species, potentially driving selection for higher investment in vegetative biomass, while harsher climatic conditions and associated reproductive uncertainty at higher latitudes could reduce selection for vegetative biomass and increased selection for high reproductive investment (stress-gradient hypothesis). Combined, these drivers could result in increased or constant reproductive allocation with increasing latitude. METHODS We quantified life-history traits in the invasive annual plant Impatiens glandulifera along a latitudinal gradient in Europe. By growing two successive glasshouse generations, we assessed genetic differentiation in vegetative growth and reproductive output across six populations, and tested whether onset of flowering drives this divergence. KEY RESULTS Trait variation was mainly caused by genetic differentiation. As expected, flowering onset was progressively earlier in populations from higher latitudes. Plant height and vegetative biomass also decreased in populations from higher latitudes, as predicted by allocation theory, but their variation was independent of the variation in flowering onset. Reproductive output remained constant across latitudes, resulting in increased reproductive allocation towards higher latitudes, supporting the stress-gradient hypothesis. We also observed trait genetic differentiation among populations that was independent of latitude. CONCLUSIONS We show that an annual invasive plant evolved several life-history traits across its invaded range in ~150 years. The evolution of vegetative and reproductive traits seems unconstrained by evolution of flowering onset. This genetic decoupling between vegetative and reproductive traits possibly contributes to the invasion success of this species.
Collapse
Affiliation(s)
- Kenny Helsen
- Plant Conservation and Population Biology, University of Leuven, Arenbergpark 31, Leuven, Belgium
- Department of Biology, Norwegian University of Science and Technology, NTNU, Trondheim, Norway
| | - Kamal Prasad Acharya
- Department of Biology, Norwegian University of Science and Technology, NTNU, Trondheim, Norway
- Department of Sports, Food and Natural Sciences, Western Norway University of Applied Sciences, Bergen, Norway
| | - Bente Jessen Graae
- Department of Biology, Norwegian University of Science and Technology, NTNU, Trondheim, Norway
| | - Hanne De Kort
- Plant Conservation and Population Biology, University of Leuven, Arenbergpark 31, Leuven, Belgium
| | - Jörg Brunet
- Swedish University of Agricultural Sciences, Southern Swedish Forest Research Centre, Alnarp, Sweden
| | - Olivier Chabrerie
- Research Unit ‘Ecologie et Dynamique des Systèmes Anthropisés’, EDYSAN, UMR 7058 CNRS, Université de Picardie Jules Verne, 1 rue des Louvels, Amiens cedex, France
| | - Sara A O Cousins
- Biogeography and Geomatics, Department of Physical Geography, Stockholm University, Stockholm, Sweden
| | - Pieter De Frenne
- Forest & Nature Lab, Ghent University, Geraardsbergsesteenweg 267, Melle-Gontrode, Belgium
| | - Martin Hermy
- Division Forest, Nature and Landscape, University of Leuven, Celestijnenlaan 200E, Leuven, Belgium
| | - Kris Verheyen
- Forest & Nature Lab, Ghent University, Geraardsbergsesteenweg 267, Melle-Gontrode, Belgium
| | - Christophe Pélabon
- Centre for Biodiversity Dynamics, Department of Biology, Norwegian University of Science and Technology, NTNU, 7491 Trondheim, Norway
| |
Collapse
|
6
|
Brousseau L, Fine PVA, Dreyer E, Vendramin GG, Scotti I. Genomic and phenotypic divergence unveil microgeographic adaptation in the Amazonian hyperdominant tree Eperua falcata Aubl. (Fabaceae). Mol Ecol 2020; 30:1136-1154. [PMID: 32786115 DOI: 10.1111/mec.15595] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2018] [Revised: 06/19/2020] [Accepted: 07/31/2020] [Indexed: 01/04/2023]
Abstract
Plant populations can undergo very localized adaptation, allowing widely distributed populations to adapt to divergent habitats in spite of recurrent gene flow. Neotropical trees-whose large and undisturbed populations often span a variety of environmental conditions and local habitats-are particularly good models to study this process. Here, we explore patterns of adaptive divergence from large (i.e., regional) to small (i.e., microgeographic) spatial scales in the hyperdominant Amazonian tree Eperua falcata Aubl. (Fabaceae) under a replicated design involving two microhabitats (~300 m apart) in two study sites (~300 km apart). A three-year reciprocal transplant illustrates that, beyond strong maternal effects and phenotypic plasticity, genetically driven divergence in seedling growth and leaf traits was detected both between seedlings originating from different regions, and between seedlings from different microhabitats. In parallel, a complementary genome scan for selection was carried out through whole-genome sequencing of tree population pools. A set of 290 divergence outlier SNPs was detected at the regional scale (between study sites), while 185 SNPs located in the vicinity of 106 protein-coding genes were detected as replicated outliers between microhabitats within regions. Outlier-surrounding genomic regions are involved in a variety of physiological processes, including plant responses to stress (e.g., oxidative stress, hypoxia and metal toxicity) and biotic interactions. Together with evidence of microgeographic divergence in functional traits, the discovery of genomic candidates for microgeographic adaptive divergence represents a promising advance in our understanding of local adaptation, which probably operates across multiple spatial scales and underpins divergence and diversification in Neotropical trees.
Collapse
Affiliation(s)
- Louise Brousseau
- UMR EcoFoG, AgroParisTech, CIRAD, CNRS, INRAE, Université de Guyane, Université des Antilles, Kourou Cedex, France.,AMAP, Univ. Montpellier, CIRAD, CNRS, INRAE, IRD, Montpellier, France
| | - Paul V A Fine
- Department of Integrative Biology, University of California, Berkeley, Berkeley, CA, USA
| | - Erwin Dreyer
- Université de Lorraine, AgroParisTech, INRAE, Silva, Nancy, France
| | - Giovanni G Vendramin
- Institute of Biosciences and BioResources (IBBR-CNR), National Research Council, Division of Florence, Sesto Fiorentino, Italy
| | - Ivan Scotti
- UR629 Ecologie des Forêts Méditerranéennes (URFM), INRAE, Avignon, France
| |
Collapse
|
7
|
MacDonald ZG, Dupuis JR, Davis CS, Acorn JH, Nielsen SE, Sperling FAH. Gene flow and climate-associated genetic variation in a vagile habitat specialist. Mol Ecol 2020; 29:3889-3906. [PMID: 32810893 DOI: 10.1111/mec.15604] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2020] [Revised: 07/08/2020] [Accepted: 07/16/2020] [Indexed: 12/31/2022]
Abstract
Previous work in landscape genetics suggests that geographic isolation is of greater importance to genetic divergence than variation in environmental conditions. This is intuitive when configurations of suitable habitat are a dominant factor limiting dispersal and gene flow, but has not been thoroughly examined for habitat specialists with strong dispersal capability. Here, we evaluate the effects of geographic and environmental isolation on genetic divergence for a vagile invertebrate with high habitat specificity and a discrete dispersal life stage: Dod's Old World swallowtail butterfly, Papilio machaon dodi. In Canada, P. m. dodi are generally restricted to eroding habitat along major river valleys where their larval host plant occurs. A series of causal and linear mixed effects models indicate that divergence of genome-wide single nucleotide polymorphisms is best explained by a combination of environmental isolation (variation in summer temperatures) and geographic isolation (Euclidean distance). Interestingly, least-cost path and circuit distances through a resistance surface parameterized as the inverse of habitat suitability were not supported. This suggests that, although habitat associations of many butterflies are specific due to reproductive requirements, habitat suitability and landscape permeability are not equivalent concepts due to considerable adult vagility. We infer that divergent selection related to variation in summer temperatures has produced two genetic clusters within P. m. dodi, differing in voltinism and diapause propensity. Within the next century, temperatures are predicted to rise by amounts greater than the present-day difference between regions of the genetic clusters, potentially affecting the persistence of the northern cluster under continued climate change.
Collapse
Affiliation(s)
- Zachary G MacDonald
- Department of Renewable Resources, University of Alberta, Edmonton, Alberta, Canada
| | - Julian R Dupuis
- Department of Entomology, University of Kentucky, Lexington, Kentucky, USA
| | - Corey S Davis
- Department of Biological Sciences, University of Alberta, Edmonton, Alberta, Canada
| | - John H Acorn
- Department of Renewable Resources, University of Alberta, Edmonton, Alberta, Canada
| | - Scott E Nielsen
- Department of Renewable Resources, University of Alberta, Edmonton, Alberta, Canada
| | - Felix A H Sperling
- Department of Biological Sciences, University of Alberta, Edmonton, Alberta, Canada
| |
Collapse
|
8
|
Knowledge status and sampling strategies to maximize cost-benefit ratio of studies in landscape genomics of wild plants. Sci Rep 2020; 10:3706. [PMID: 32111897 PMCID: PMC7048820 DOI: 10.1038/s41598-020-60788-8] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2019] [Accepted: 02/11/2020] [Indexed: 11/27/2022] Open
Abstract
To avoid local extinction due to the changes in their natural ecosystems, introduced by anthropogenic activities, species undergo local adaptation. Landscape genomics approach, through genome–environment association studies, has helped evaluate the local adaptation in natural populations. Landscape genomics, is still a developing discipline, requiring refinement of guidelines in sampling design, especially for studies conducted in the backdrop of stark socioeconomic realities of the rainforest ecologies, which are global biodiversity hotspots. In this study we aimed to devise strategies to improve the cost-benefit ratio of landscape genomics studies by surveying sampling designs and genome sequencing strategies used in existing studies. We conducted meta-analyses to evaluate the importance of sampling designs, in terms of (i) number of populations sampled, (ii) number of individuals sampled per population, (iii) total number of individuals sampled, and (iv) number of SNPs used in different studies, in discerning the molecular mechanisms underlying local adaptation of wild plant species. Using the linear mixed effects model, we demonstrated that the total number of individuals sampled and the number of SNPs used, significantly influenced the detection of loci underlying the local adaptation. Thus, based on our findings, in order to optimize the cost-benefit ratio of landscape genomics studies, we suggest focusing on increasing the total number of individuals sampled and using a targeted (e.g. sequencing capture) Pool-Seq approach and/or a random (e.g. RAD-Seq) Pool-Seq approach to detect SNPs and identify SNPs under selection for a given environmental cline. We also found that the existing molecular evidences are inadequate in predicting the local adaptations to climate change in tropical forest ecosystems.
Collapse
|
9
|
Gamboa M, Watanabe K. Genome-wide signatures of local adaptation among seven stoneflies species along a nationwide latitudinal gradient in Japan. BMC Genomics 2019; 20:84. [PMID: 30678640 PMCID: PMC6346529 DOI: 10.1186/s12864-019-5453-3] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2018] [Accepted: 01/14/2019] [Indexed: 11/16/2022] Open
Abstract
Background Environmental heterogeneity continuously produces a selective pressure that results in genomic variation among organisms; understanding this relationship remains a challenge in evolutionary biology. Here, we evaluated the degree of genome-environmental association of seven stonefly species across a wide geographic area in Japan and additionally identified putative environmental drivers and their effect on co-existing multiple stonefly species. Double-digest restriction-associated DNA (ddRAD) libraries were independently sequenced for 219 individuals from 23 sites across four geographical regions along a nationwide latitudinal gradient in Japan. Results A total of 4251 candidate single nucleotide polymorphisms (SNPs) strongly associated with local adaptation were discovered using Latent mixed models; of these, 294 SNPs showed strong correlation with environmental variables, specifically precipitation and altitude, using distance-based redundancy analysis. Genome–genome comparison among the seven species revealed a high sequence similarity of candidate SNPs within a geographical region, suggesting the occurrence of a parallel evolution process. Conclusions Our results revealed genomic signatures of local adaptation and their influence on multiple, co-occurring species. These results can be potentially applied for future studies on river management and climatic stressor impacts. Electronic supplementary material The online version of this article (10.1186/s12864-019-5453-3) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Maribet Gamboa
- Department of Civil and Environmental Engineering, Ehime University, Matsuyama, 790-0871, Japan.
| | - Kozo Watanabe
- Department of Civil and Environmental Engineering, Ehime University, Matsuyama, 790-0871, Japan
| |
Collapse
|
10
|
Martins K, Gugger PF, Llanderal‐Mendoza J, González‐Rodríguez A, Fitz‐Gibbon ST, Zhao J, Rodríguez‐Correa H, Oyama K, Sork VL. Landscape genomics provides evidence of climate-associated genetic variation in Mexican populations of Quercus rugosa. Evol Appl 2018; 11:1842-1858. [PMID: 30459833 PMCID: PMC6231481 DOI: 10.1111/eva.12684] [Citation(s) in RCA: 40] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2017] [Revised: 05/31/2018] [Accepted: 07/11/2018] [Indexed: 12/30/2022] Open
Abstract
Local adaptation is a critical evolutionary process that allows plants to grow better in their local compared to non-native habitat and results in species-wide geographic patterns of adaptive genetic variation. For forest tree species with a long generation time, this spatial genetic heterogeneity can shape the ability of trees to respond to rapid climate change. Here, we identify genomic variation that may confer local environmental adaptations and then predict the extent of adaptive mismatch under future climate as a tool for forest restoration or management of the widely distributed high-elevation oak species Quercus rugosa in Mexico. Using genotyping by sequencing, we identified 5,354 single nucleotide polymorphisms (SNPs) genotyped from 103 individuals across 17 sites in the Trans-Mexican Volcanic Belt, and, after controlling for neutral genetic structure, we detected 74 F ST outlier SNPs and 97 SNPs associated with climate variation. Then, we deployed a nonlinear multivariate model, Gradient Forests, to map turnover in allele frequencies along environmental gradients and predict areas most sensitive to climate change. We found that spatial patterns of genetic variation were most strongly associated with precipitation seasonality and geographic distance. We identified regions of contemporary genetic and climatic similarities and predicted regions where future populations of Q. rugosa might be at risk due to high expected rate of climate change. Our findings provide preliminary details for future management strategies of Q. rugosa in Mexico and also illustrate how a landscape genomic approach can provide a useful tool for conservation and resource management strategies.
Collapse
Affiliation(s)
- Karina Martins
- Department of Ecology and Evolutionary BiologyUniversity of California, Los AngelesLos AngelesCalifornia
- Departamento de BiologiaUniversidade Federal de São CarlosSorocabaSPBrazil
| | - Paul F. Gugger
- Department of Ecology and Evolutionary BiologyUniversity of California, Los AngelesLos AngelesCalifornia
- Appalachian LaboratoryUniversity of Maryland Center for Environmental ScienceFrostburgMaryland
| | - Jesus Llanderal‐Mendoza
- Instituto de Investigaciones en Ecosistemas y SustentabilidadUniversidad Nacional Autónoma de México (UNAM)MoreliaMichoacánMéxico
- Escuela Nacional de Estudios Superiores Unidad MoreliaUniversidad Nacional Autónoma de México (UNAM)MoreliaMichoacánMéxico
| | - Antonio González‐Rodríguez
- Instituto de Investigaciones en Ecosistemas y SustentabilidadUniversidad Nacional Autónoma de México (UNAM)MoreliaMichoacánMéxico
| | - Sorel T. Fitz‐Gibbon
- Department of Ecology and Evolutionary BiologyUniversity of California, Los AngelesLos AngelesCalifornia
| | - Jian‐Li Zhao
- Key Laboratory of Tropical Forest EcologyXishuangbanna Tropical Botanical GardenChinese Academy of SciencesMenglaYunnanChina
| | - Hernando Rodríguez‐Correa
- Escuela Nacional de Estudios Superiores Unidad MoreliaUniversidad Nacional Autónoma de México (UNAM)MoreliaMichoacánMéxico
| | - Ken Oyama
- Escuela Nacional de Estudios Superiores Unidad MoreliaUniversidad Nacional Autónoma de México (UNAM)MoreliaMichoacánMéxico
| | - Victoria L. Sork
- Department of Ecology and Evolutionary BiologyUniversity of California, Los AngelesLos AngelesCalifornia
- Institute of the Environment and SustainabilityUniversity of California, Los AngelesLos AngelesCalifornia
| |
Collapse
|
11
|
Konar A, Choudhury O, Bullis R, Fiedler L, Kruser JM, Stephens MT, Gailing O, Schlarbaum S, Coggeshall MV, Staton ME, Carlson JE, Emrich S, Romero-Severson J. High-quality genetic mapping with ddRADseq in the non-model tree Quercus rubra. BMC Genomics 2017; 18:417. [PMID: 28558688 PMCID: PMC5450186 DOI: 10.1186/s12864-017-3765-8] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2016] [Accepted: 05/04/2017] [Indexed: 11/10/2022] Open
Abstract
Background Restriction site associated DNA sequencing (RADseq) has the potential to be a broadly applicable, low-cost approach for high-quality genetic linkage mapping in forest trees lacking a reference genome. The statistical inference of linear order must be as accurate as possible for the correct ordering of sequence scaffolds and contigs to chromosomal locations. Accurate maps also facilitate the discovery of chromosome segments containing allelic variants conferring resistance to the biotic and abiotic stresses that threaten forest trees worldwide. We used ddRADseq for genetic mapping in the tree Quercus rubra, with an approach optimized to produce a high-quality map. Our study design also enabled us to model the results we would have obtained with less depth of coverage. Results Our sequencing design produced a high sequencing depth in the parents (248×) and a moderate sequencing depth (15×) in the progeny. The digital normalization method of generating a de novo reference and the SAMtools SNP variant caller yielded the most SNP calls (78,725). The major drivers of map inflation were multiple SNPs located within the same sequence (77% of SNPs called). The highest quality map was generated with a low level of missing data (5%) and a genome-wide threshold of 0.025 for deviation from Mendelian expectation. The final map included 849 SNP markers (1.8% of the 78,725 SNPs called). Downsampling the individual FASTQ files to model lower depth of coverage revealed that sequencing the progeny using 96 samples per lane would have yielded too few SNP markers to generate a map, even if we had sequenced the parents at depth 248×. Conclusions The ddRADseq technology produced enough high-quality SNP markers to make a moderately dense, high-quality map. The success of this project was due to high depth of coverage of the parents, moderate depth of coverage of the progeny, a good framework map, an optimized bioinformatics pipeline, and rigorous premapping filters. The ddRADseq approach is useful for the construction of high-quality genetic maps in organisms lacking a reference genome if the parents and progeny are sequenced at sufficient depth. Technical improvements in reduced representation sequencing (RRS) approaches are needed to reduce the amount of missing data. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3765-8) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Arpita Konar
- Department of Biological Sciences, University of Notre Dame, Notre Dame, IN, 46556, USA
| | - Olivia Choudhury
- Department of Computer Science and Engineering, University of Notre Dame, Notre Dame, IN, 46556, USA
| | - Rebecca Bullis
- Department of Biological Sciences, University of Notre Dame, Notre Dame, IN, 46556, USA
| | - Lauren Fiedler
- Department of Biological Sciences, University of Notre Dame, Notre Dame, IN, 46556, USA
| | | | - Melissa T Stephens
- Department of Biological Sciences, University of Notre Dame, Notre Dame, IN, 46556, USA
| | - Oliver Gailing
- School of Forest Resources and Environmental Science, Michigan Technological University, Houghton, MI, 49931, USA
| | - Scott Schlarbaum
- Department of Forestry, Wildlife and Fisheries, University of Tennessee, Knoxville, TN, 37996, USA
| | - Mark V Coggeshall
- School of Natural Resources, University of Missouri-Columbia, Columbia, MO, 65211, USA.,Hardwood Tree Improvement and Regeneration Center, USDA Forest Service Northern Research Station, West Lafayette, IN, 47907, USA
| | - Margaret E Staton
- Department of Entomology and Plant Pathology, University of Tennessee, Knoxville, TN, 37996, USA
| | - John E Carlson
- Department of Ecosystem Science and Management, Penn State, University Park, State College, PA, 16802, USA
| | - Scott Emrich
- Department of Computer Science and Engineering, University of Notre Dame, Notre Dame, IN, 46556, USA
| | - Jeanne Romero-Severson
- Department of Biological Sciences, University of Notre Dame, Notre Dame, IN, 46556, USA.
| |
Collapse
|
12
|
De Kort H, Mergeay J, Jacquemyn H, Honnay O. Transatlantic invasion routes and adaptive potential in North American populations of the invasive glossy buckthorn, Frangula alnus. ANNALS OF BOTANY 2016; 118:1089-1099. [PMID: 27539599 PMCID: PMC5091722 DOI: 10.1093/aob/mcw157] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/24/2015] [Revised: 05/02/2016] [Accepted: 06/17/2016] [Indexed: 05/31/2023]
Abstract
BACKGROUND AND AIMS Many invasive species severely threaten native biodiversity and ecosystem functioning. One of the most prominent questions in invasion genetics is how invasive populations can overcome genetic founder effects to establish stable populations after colonization of new habitats. High native genetic diversity and multiple introductions are expected to increase genetic diversity and adaptive potential in the invasive range. Our aim was to identify the European source populations of Frangula alnus (glossy buckthorn), an ornamental and highly invasive woody species that was deliberately introduced into North America at the end of the 18th century. A second aim of this study was to assess the adaptive potential as an explanation for the invasion success of this species. METHODS Using a set of annotated single-nucleotide polymorphisms (SNPs) that were assigned a putative function based on sequence comparison with model species, a total of 38 native European and 21 invasive North American populations were subjected to distance-based structure and assignment analyses combined with population genomic tools. Genetic diversity at SNPs with ecologically relevant functions was considered as a proxy for adaptive potential. KEY RESULTS Patterns of invasion coincided with early modern transatlantic trading routes. Multiple introductions through transatlantic trade from a limited number of European port regions to American urban areas led to the establishment of bridgehead populations with high allelic richness and expected heterozygosity, allowing continuous secondary migration to natural areas. CONCLUSIONS Targeted eradication of the urban populations, where the highest genetic diversity and adaptive potential were observed, offers a promising strategy to arrest further invasion of native American prairies and forests.
Collapse
Affiliation(s)
- Hanne De Kort
- Plant Conservation and Population Biology, Biology Department, University of Leuven, Kasteelpark Arenberg 31, B-3001 Heverlee, Belgium Station d'Ecologie Théorique et Expérimentale du CNRS, Centre National de la Recherche Scientifique, 2 Route du CNRS, FR-09200 Moulis, France
| | - Joachim Mergeay
- Research Institute for Nature and Forest, Gaverstraat 4, B-9500 Geraardsbergen, Belgium
| | - Hans Jacquemyn
- Plant Conservation and Population Biology, Biology Department, University of Leuven, Kasteelpark Arenberg 31, B-3001 Heverlee, Belgium
| | - Olivier Honnay
- Plant Conservation and Population Biology, Biology Department, University of Leuven, Kasteelpark Arenberg 31, B-3001 Heverlee, Belgium
| |
Collapse
|