Geographic patterns of genomic diversity and structure in the C4 grass Panicum hallii across its natural distribution

Abstract Geographic patterns of within-species genomic diversity are shaped by evolutionary processes, life history and historical and contemporary factors. New genomic approaches can be used to infer the influence of such factors on the current distribution of infraspecific lineages. In this study, we evaluated the genomic and morphological diversity as well as the genetic structure of the C4 grass Panicum hallii across its complex natural distribution in North America. We sampled extensively across the natural range of P. hallii in Mexico and the USA to generate double-digestion restriction-associated DNA (ddRAD) sequence data for 423 individuals from 118 localities. We used these individuals to study the divergence between the two varieties of P. hallii, P. hallii var. filipes and P. hallii var. hallii as well as the genetic diversity and structure within these groups. We also examined the possibility of admixture in the geographically sympatric zone shared by both varieties, and assessed distribution shifts related with past climatic fluctuations. There is strong genetic and morphological divergence between the varieties and consistent genetic structure defining seven genetic clusters that follow major ecoregions across the range. South Texas constitutes a hotspot of genetic diversity with the co-occurrence of all genetic clusters and admixture between the two varieties. It is likely a recolonization and convergence point of populations that previously diverged in isolation during fragmentation events following glaciation periods.

[1]  J. Schmutz,et al.  Population genomics and climate adaptation of a C4 perennial grass, Panicum hallii (Poaceae) , 2018, BMC Genomics.

[2]  Nicolas C Rochette,et al.  Deriving genotypes from RAD-seq short-read data using Stacks , 2017, Nature Protocols.

[3]  Aaron B. A. Shafer,et al.  Bioinformatic processing of RAD‐seq data dramatically impacts downstream population genetic inference , 2017 .

[4]  L. Seeb,et al.  Paralogs are revealed by proportion of heterozygotes and deviations in read ratios in genotyping‐by‐sequencing data from natural populations , 2017, Molecular ecology resources.

[5]  B. Emerson,et al.  ddradseqtools: a software package for in silico simulation and testing of double‐digest RADseq experiments , 2017, Molecular ecology resources.

[6]  Niklaus J Grünwald,et al.  vcfr: a package to manipulate and visualize variant call format data in R , 2017, Molecular ecology resources.

[7]  Anders Albrechtsen,et al.  Next-generation biology: Sequencing and data analysis approaches for non-model organisms. , 2016, Marine genomics.

[8]  Thomas E. Juenger,et al.  The Genetic Basis of Upland/Lowland Ecotype Divergence in Switchgrass (Panicum virgatum) , 2016, G3: Genes, Genomes, Genetics.

[9]  S. Hazen,et al.  Population Structure in the Model Grass Brachypodium distachyon Is Highly Correlated with Flowering Differences across Broad Geographic Areas , 2016, The plant genome.

[10]  H. Ellegren,et al.  Determinants of genetic diversity , 2016, Nature Reviews Genetics.

[11]  A. Márcia Barbosa,et al.  fuzzySim: applying fuzzy logic to binary similarity indices in ecology , 2015 .

[12]  Matthew E. Aiello-Lammens,et al.  spThin: an R package for spatial thinning of species occurrence records for use in ecological niche models , 2015 .

[13]  U. Mueller,et al.  Gone to Texas: phylogeography of two Trachymyrmex (Hymenoptera: Formicidae) species along the southeastern coastal plain of North America , 2015 .

[14]  E. Ortiz,et al.  The grass family (Poaceae) in Coahuila, Mexico: Diversity and distribution , 2015 .

[15]  Jennifer K. Costanza,et al.  How global biodiversity hotspots may go unrecognized: lessons from the North American Coastal Plain , 2015 .

[16]  B. Emerson,et al.  Restriction site‐associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference , 2015, Molecular ecology resources.

[17]  A. von Haeseler,et al.  IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies , 2014, Molecular biology and evolution.

[18]  O. Lepais,et al.  SimRAD: an R package for simulation‐based prediction of the number of loci expected in RADseq and similar genotyping by sequencing approaches , 2014, Molecular ecology resources.

[19]  J. Chapman,et al.  The genetics of divergence and reproductive isolation between ecotypes of Panicum hallii , 2014, The New phytologist.

[20]  James M Sobel,et al.  UNIFICATION OF METHODS FOR ESTIMATING THE STRENGTH OF REPRODUCTIVE ISOLATION , 2014, Evolution; international journal of organic evolution.

[21]  J. Barnett,et al.  Forgotten Grasslands of the South: Natural History and Conservation , 2012 .

[22]  Brian Bushnell,et al.  BBMap: A Fast, Accurate, Splice-Aware Aligner , 2014 .

[23]  Robert A. Boria,et al.  Spatial filtering to reduce sampling bias can improve the performance of ecological niche models , 2014 .

[24]  P. Pearman,et al.  Caves as microrefugia: Pleistocene phylogeography of the troglophilic North American scorpion Pseudouroctonus reddelli , 2014, BMC Evolutionary Biology.

[25]  Boris Schröder,et al.  The importance of correcting for sampling bias in MaxEnt species distribution models , 2013 .

[26]  J. Stinchcombe,et al.  The effect of leaf shape on the thermoregulation and frost tolerance of an annual vine, Ipomoea hederacea (Convolvulaceae). , 2013, American journal of botany.

[27]  A. Anjos,et al.  P-TRAP: a Panicle Trait Phenotyping tool , 2013, BMC Plant Biology.

[28]  M. Beaumont,et al.  The shaping of genetic variation in edge-of-range populations under past and future climate change , 2013, Ecology letters.

[29]  B. Stevens,et al.  Climate and carbon cycle changes from 1850 to 2100 in MPI‐ESM simulations for the Coupled Model Intercomparison Project phase 5 , 2013 .

[30]  Angel Amores,et al.  Stacks: an analysis tool set for population genomics , 2013, Molecular ecology.

[31]  P. Reich,et al.  New handbook for standardised measurement of plant functional traits worldwide , 2013 .

[32]  Christopher Phillips,et al.  An overview of STRUCTURE: applications, parameter settings, and supporting software , 2013, Front. Genet..

[33]  T. Juenger,et al.  A population genetic transect of Panicum hallii (Poaceae). , 2013, American journal of botany.

[34]  Minh Anh Nguyen,et al.  Ultrafast Approximation for Phylogenetic Bootstrap , 2013, Molecular biology and evolution.

[35]  H. Hoekstra,et al.  Double Digest RADseq: An Inexpensive Method for De Novo SNP Discovery and Genotyping in Model and Non-Model Species , 2012, PloS one.

[36]  Shane S. Sturrock,et al.  Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data , 2012, Bioinform..

[37]  T. Juenger,et al.  Microsatellite markers for the native Texas perennial grass, Panicum hallii (Poaceae). , 2012, American journal of botany.

[38]  D. Soltis,et al.  Cytogeography of the Humifusa clade of Opuntia s.s. Mill. 1754 (Cactaceae, Opuntioideae, Opuntieae): correlations with pleistocene refugia and morphological traits in a polyploid complex , 2012, Comparative cytogenetics.

[39]  A. Amores,et al.  Stacks: Building and Genotyping Loci De Novo From Short-Read Sequences , 2011, G3: Genes | Genomes | Genetics.

[40]  Gonçalo R. Abecasis,et al.  The variant call format and VCFtools , 2011, Bioinform..

[41]  P. Schönswetter,et al.  Multiple Pleistocene refugia and Holocene range expansion of an abundant southwestern American desert plant species (Melampodium leucanthum, Asteraceae) , 2010, Molecular ecology.

[42]  Eduardo Estrada-Castillón,et al.  Flora del centro del estado de Chihuahua, México , 2010 .

[43]  M. Araújo,et al.  BIOMOD – a platform for ensemble forecasting of species distributions , 2009 .

[44]  Joe Hereford,et al.  A Quantitative Survey of Local Adaptation and Fitness Trade‐Offs , 2009, The American Naturalist.

[45]  K. Wright,et al.  The strength and genetic basis of reproductive isolating barriers in flowering plants , 2008, Philosophical Transactions of the Royal Society B: Biological Sciences.

[46]  Thibaut Jombart,et al.  adegenet: a R package for the multivariate analysis of genetic markers , 2008, Bioinform..

[47]  Mark Vellend,et al.  Ecological consequences of genetic diversity. , 2008, Ecology letters.

[48]  J. Betancourt,et al.  Inferences about winter temperatures and summer rains from the late Quaternary record of C4 perennial grasses and C3 desert shrubs in the northern Chihuahuan Desert , 2007 .

[49]  Omri Allouche,et al.  Assessing the accuracy of species distribution models: prevalence, kappa and the true skill statistic (TSS) , 2006 .

[50]  Robert P. Anderson,et al.  Maximum entropy modeling of species geographic distributions , 2006 .

[51]  J. L. Parra,et al.  Very high resolution interpolated climate surfaces for global land areas , 2005 .

[52]  T. Hastie,et al.  Using multivariate adaptive regression splines to predict the distributions of New Zealand ’ s freshwater diadromous fish , 2005 .

[53]  W. Thuiller,et al.  Predicting species distribution: offering more than simple habitat models. , 2005, Ecology letters.

[54]  N. Swenson,et al.  Clustering of Contact Zones, Hybrid Zones, and Phylogeographic Breaks in North America , 2005, The American Naturalist.

[55]  M. Vellend,et al.  Connections between species diversity and genetic diversity , 2005 .

[56]  N. Swenson,et al.  DO SUTURE ZONES EXIST? , 2004, Evolution; international journal of organic evolution.

[57]  G. Danabasoglu,et al.  The Community Climate System Model Version 4 , 2011 .

[58]  T. Stuessy,et al.  Diploid and polyploid cytotype distribution in Melampodium cinereum and M. leucanthum (Asteraceae, Heliantheae). , 2004, American journal of botany.

[59]  G. Hewitt Genetic consequences of climatic oscillations in the Quaternary. , 2004, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[60]  Jeffrey K. Conner,et al.  A Primer of Ecological Genetics , 2004 .

[61]  Eric R. Ziegel,et al.  Generalized Linear Models , 2002, Technometrics.

[62]  G. Powell,et al.  Terrestrial Ecoregions of the World: A New Map of Life on Earth , 2001 .

[63]  Leo Breiman,et al.  Random Forests , 2001, Machine Learning.

[64]  G. Hewitt Speciation, hybrid zones and phylogeography — or seeing genes in space and time , 2001, Molecular ecology.

[65]  G. Hewitt The genetic legacy of the Quaternary ice ages , 2000, Nature.

[66]  P. Donnelly,et al.  Inference of population structure using multilocus genotype data. , 2000, Genetics.

[67]  J. Hamrick,et al.  Effects of life history traits on genetic diversity in plant species , 1996 .

[68]  S. Elias Late Quaternary zoogeography of the chihuahuan desert insect fauna, based on fossil records from packrat middens , 1992 .

[69]  F. Smeins,et al.  Composition, Classification and Species Response Patterns of Remnant Tallgrass Prairies in Texas , 1985 .

[70]  S. Renvoize,et al.  The Grasses of Texas , 1976 .

[71]  Floyd R. Waller A New Combination in Panicum (Gramineae) , 1974 .

[72]  H. G. Baker Seed Weight in Relation to Environmental Conditions in California , 1972 .

[73]  Edzer Pebesma,et al.  Simple Features for R: Standardized Support for Spatial Vector Data , 2018, R J..

[74]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[75]  Wei-Yin Loh,et al.  Classification and regression trees , 2011, WIREs Data Mining Knowl. Discov..

[76]  Y. H. Arrieta,et al.  LISTADO FLORISTICO Y ASPECTOS ECOLÓGICOS DE LA FAMILIA POACEAE PARA CHIHUAHUA, DURANGO Y ZACATECAS, MEXICO , 2010 .

[77]  J. Peñuelas,et al.  Environmental change and the option value of genetic diversity. , 2009, Trends in plant science.

[78]  G. De’ath Boosted trees for ecological modeling and prediction. , 2007, Ecology.

[79]  Mark New,et al.  Ensemble forecasting of species distributions. , 2007, Trends in ecology & evolution.

[80]  M. Flores-Vergara,et al.  A modified protocol for rapid DNA isolation from plant tissues using cetyltrimethylammonium bromide , 2006, Nature Protocols.

[81]  Hiroyasu Hasumi,et al.  K-1 Coupled GCM (MIROC) Description , 2004 .

[82]  Michael H. MacRoberts,et al.  ENDEMISM IN THE WEST GULF COASTAL PLAIN: IMPORTANCE OF XERIC HABITATS , 2002 .

[83]  Stefan Sperlich,et al.  Generalized Additive Models , 2014 .

[84]  C. Remington,et al.  Suture-Zones of Hybrid Interaction Between Recently Joined Biotas , 1968 .