Multiple Pleistocene refugia and Holocene range expansion of an abundant southwestern American desert plant species (Melampodium leucanthum, Asteraceae)

Origin and diversification of the saguaro cactus (Carnegiea gigantea): a within-species phylogenomic analysis

Reconstructing accurate historical relationships within a species poses numerous challenges, not least in many plant groups in which gene flow is high enough to extend well beyond species boundaries. Nonetheless, the extent of tree-like history within a species is an empirical question on which it is now possible to bring large amounts of genome sequence to bear. We assess phylogenetic structure across the geographic range of the saguaro cactus, an emblematic member of Cactaceae, a clade known for extensive hybridization and porous species boundaries. Using 200 Gb of whole genome resequencing data from 20 individuals sampled from 10 localities, we assembled two data sets comprising 150,000 biallelic single nucleotide polymorphisms (SNPs) from protein coding sequences. From these we inferred within-species trees and evaluated their significance and robustness using five qualitatively different inference methods. Despite the low sequence diversity, large census population sizes, and presence of wide-ranging pollen and seed dispersal agents, phylogenetic trees were well resolved and highly consistent across both data sets and all methods. We inferred that the most likely root, based on marginal likelihood comparisons, is to the east and south of the region of highest genetic diversity, which lies along the coast of the Gulf of California in Sonora, Mexico. Together with striking decreases in marginal likelihood found to the north, this supports hypotheses that saguaro's current range reflects post-glacial expansion from the refugia in the south of its range. We conclude with observations about practical and theoretical issues raised by phylogenomic data sets within species, in which SNP-based methods must be used rather than gene tree methods that are widely used when sequence divergence is higher. These include computational scalability, inference of gene flow, and proper assessment of statistical support in the presence of linkage effects.

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

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 C₄ 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.

Fluctuating fortunes: genomes and habitat reconstructions reveal global climate-mediated changes in bats' genetic diversity

Over the last approximately 2.6 Myr, Earth's climate has been dominated by cyclical ice ages that have profoundly affected species' population sizes, but the impact of impending anthropogenic climate change on species' extinction potential remains a worrying problem. We investigated 11 bat species from different taxonomic, ecological and geographical backgrounds using combined information from palaeoclimatic habitat reconstructions and genomes to analyse biotic impacts of historic climate change. We discover tightly correlated fluctuations between species' historic distribution and effective population size, identify frugivores as particularly susceptible to global warming, pinpoint large insectivores as having overall low effective population size and flag the onset of the Holocene (approx. 10-12 000 years ago) as the period with the generally lowest effective population sizes across the last approximately 1 Myr. Our study shows that combining genomic and palaeoclimatological approaches reveals effects of climatic shifts on genetic diversity and may help predict impacts of future climate change.

Population genomics and climate adaptation of a C4 perennial grass, Panicum hallii (Poaceae)

Background

Understanding how and why genetic variation is partitioned across geographic space is of fundamental importance to understanding the nature of biological species. How geographical isolation and local adaptation contribute to the formation of ecotypically differentiated groups of plants is just beginning to be understood through population genomic studies. We used whole genome sequencing combined with association study of climate to discover the drivers of differentiation in the perennial C4 grass Panicum hallii.

Results

Sequencing of 89 natural accessions of P.hallii revealed complex population structure across the species range. Major population genomic separation was found between subspecies P.hallii var. hallii and var. filipes as well as between at least four major unrecognized subgroups within var. hallii. At least 139 genomic SNPs were significantly associated with temperature or precipitation across the range and these SNPs were enriched for non-synonymous substitutions. SNPs associated with temperature and aridity were more often found in or near genes than expected by chance and enriched for putative involvement in dormancy processes, seed maturation, response to hyperosmosis and salinity, abscisic acid metabolism, hormone metabolism, and drought recovery.

Conclusions

Both geography and climate adaptation contribute significantly to patterns of genome-wide variation in P.hallii. Population subgroups within P.hallii may represent early stages in the formation of ecotypes. Climate associated loci identified here represent promising targets for future research in this and other perennial grasses.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5179-7) contains supplementary material, which is available to authorized users.

Geographical features are the predominant driver of molecular diversification in widely distributed North American whipsnakes

Allopatric divergence following the formation of geographical features has been implicated as a major driver of evolutionary diversification. Widespread species complexes provide opportunities to examine allopatric divergence across varying degrees of isolation in both time and space. In North America, several geographical features may play such a role in diversification, including the Mississippi River, Pecos River, Rocky Mountains, Cochise Filter Barrier, Gulf of California and Isthmus of Tehuantepec. We used thousands of nuclear single nucleotide polymorphisms (SNPs) and mitochondrial DNA from several species of whipsnakes (genera Masticophis and Coluber) distributed across North and Central America to investigate the role that these geographical features have played on lineage divergence. We hypothesize that these features restrict gene flow and separate whipsnakes into diagnosable genomic clusters. We performed genomic clustering and phylogenetic reconstructions at the species and population levels using Bayesian and likelihood analyses and quantified migration levels across geographical features to assess the degree of genetic isolation due to allopatry. Our analyses suggest that (i) major genetic divisions are often consistent with isolation by geographical features, (ii) migration rates between clusters are asymmetrical across major geographical features, and (iii) areas that receive proportionally more migrants possess higher levels of genetic diversity. Collectively, our findings suggest that multiple features of the North American landscape contributed to allopatric divergence in this widely distributed snake group.

Shifting Quaternary migration patterns in the Bahamian archipelago: Evidence from the Zamia pumila complex at the northern limits of the Caribbean island biodiversity hotspot

PREMISE OF THE STUDY

The Bahamas archipelago is formed by young, tectonically stable carbonate banks that harbor direct geological evidence of global ice-volume changes. We sought to detect signatures of major changes on gene flow patterns and reconstruct the phylogeographic history of the monophyletic Zamia pumila complex across the Bahamas.

METHODS

Nuclear molecular markers with both high and low mutation rates were used to capture two different time scale signatures and test several gene flow and demographic hypotheses.

KEY RESULTS

Single-copy nuclear genes unveiled apparent ancestral admixture on Andros, suggesting a significant role of this island as main hub of diversity of the archipelago. We detected demographic and spatial expansion of the Zamia pumila complex on both paleo-provinces around the Piacenzian (Pliocene)/Gelasian (Pleistocene). Populations evidenced signatures of different migration models that have occurred at two different times. Populations on Long Island (Z. lucayana) may either represent a secondary colonization of the Bahamas by Zamia or a rapid and early-divergence event of at least one population on the Bahamas.

CONCLUSIONS

Despite changes in migration patterns with global climate, expected heterozygosity with both marker systems remains within the range reported for cycads, but with significant levels of increased inbreeding detected by the microsatellites. This finding is likely associated with reduced gene flow between and within paleo-provinces, accompanied by genetic drift, as rising seas enforced isolation. Our study highlights the importance of the maintenance of the predominant direction of genetic exchange and the role of overseas dispersion among the islands during climate oscillations.

Colonization in North American Arid Lands: The Journey of Agarito (Berberis trifoliolata) Revealed by Multilocus Molecular Data and Packrat Midden Fossil Remains

Here we conduct research to understand the evolutionary history of a shrubby species known as Agarito (Berberis trifoliolata), an endemic species to the Chihuahuan Desert. We identify genetic signatures based on plastid DNA and AFLP markers and perform niche modelling and spatial connectivity analyses as well as niche modelling based on records in packrats to elucidate whether orogenic events such as mountain range uplift in the Miocene or the contraction/expansion dynamics of vegetation in response to climate oscillations in the Pliocene/Pleistocene had an effect on evolutionary processes in Agarito. Our results of current niche modelling and palaeomodelling showed that the area currently occupied by Berberis trifoliolata is substantially larger than it was during the Last Interglacial period and the Last Glacial Maximum. Agarito was probably confined to small areas in the Northeastern and gradually expanded its distribution just after the Last Glacial Maximum when the weather in the Chihuahuan Desert and adjacent regions became progressively warmer and drier. The most contracted range was predicted for the Interglacial period. Populations remained in stable areas during the Last Glacial Maximum and expanded at the beginning of the Holocene. Most genetic variation occured in populations from the Sierra Madre Oriental. Two groups of haplotypes were identified: the Mexican Plateau populations and certain Northeastern populations. Haplogroups were spatially connected during the Last Glacial Maximum and separated during interglacial periods. The most important prediction of packrat middens palaeomodelling lies in the Mexican Plateau, a finding congruent with current and past niche modelling predictions for agarito and genetic results. Our results corroborate that these climate changes in the Pliocene/Pleistocene affected the evolutionary history of agarito. The journey of agarito in the Chihuahuan Desert has been dynamic, expanding and contracting its distribution range and currently occupying the largest area in its history.

New insights on the origin of the woody flora of the Chihuahuan Desert: The case of Lindleya

PREMISE OF THE STUDY

We reconstructed aspects of the evolutionary history of Lindleya, a shrubby element in the rose family, as a model for identifying the ancestral area of the biota of the Chihuahuan Desert and for understanding the effect that pre-Pleistocene or Pleistocene climate fluctuations had on the evolution of the plants of this desert.

METHODS

We conducted phylogenetic, biogeographic, and phylogeographic analyses based on one nuclear DNA and two plastid markers sequenced for 20 populations comprising 153 accessions and representing the species' entire distribution range. We carried out ecological niche modeling to forecast Lindleya's potential distribution under Current, Last Interglacial, and Last Glacial Maximum conditions.

KEY RESULTS

The ancestral area of Lindleya populations was probably in the central and southern parts of the Chihuahuan Desert. Levels of genetic and haplotype variation were the highest in the northernmost and southernmost areas and low levels of genetic variation in the central part. Extended Bayesian skyline plots including all populations identified a slight demographic expansion ∼1 Ma. The predicted potential distribution during the Last Glacial Maximum was very restricted.

CONCLUSIONS

We suggest that pre-Pleistocene changes had an influence on the divergence of Lindleya populations. We also suggest that the ancestral area for Lindleya was the Chihuahuan Desert and that the uplift of the Trans-Mexican Volcanic Belt isolated the Tehuacán Valley populations. Genetic and demographic analyses as well as and ecological niche modeling indicate that populations of Lindleya experienced genetic bottlenecks and that they have expanded since the Last Glacial Maximum.

Ecological distributions, phenological isolation, and genetic structure in sympatric and parapatric populations of the Larrea tridentata polyploid complex

PREMISE OF THE STUDY

Polyploidy is widely recognized as a mechanism of diversification. Contributions of polyploidy to specific pre- and postzygotic barriers-and classifications of polyploid speciation as "ecological" vs. "non-ecological"-are more contentious. Evaluation of these issues requires comprehensive studies that test ecological characteristics of cytotypes as well as the coincidence of genetic structure with cytotype distributions.

METHODS

We investigated a classical example of autopolyploid speciation, Larrea tridentata, at multiple areas of cytotype co-occurrence. Habitat and phenological differences were compared between diploid, tetraploid, and hexaploid populations on the basis of edaphic, community composition, and flowering time surveys. Frequency of hybridization between diploids and tetraploids was investigated using a diploid-specific chloroplast DNA (cpDNA) marker; genetic structure for all cytotypes was assessed using amplified fragment length polymorphisms (AFLPs).

KEY RESULTS

Across contact zones, we found cytotypes in habitats distinguished by soil and vegetation. We observed modest differences in timing and production of flowers, indicating a degree of assortative mating that was asymmetric between cytotypes. Nonetheless, cpDNA analyses in diploid-tetraploid contact zones suggested that ∼5% of tetraploid plants had hybrid origins involving unilateral sexual polyploidization. Genetic structure of AFLPs largely coincided with cytotype distributions in diploid-tetraploid contact zones. In contrast, there was little structure in areas of contact between tetraploids and hexaploids, suggesting intercytotype gene flow or recurrent hexaploid formation.

CONCLUSIONS

Diploid, tetraploid, and hexaploid cytotypes of L. tridentata are segregated by environmental distributions and flowering phenology in contact zones, with diploid and tetraploid populations having corresponding differences in genetic structure.

Historical Biogeography Using Species Geographical Ranges

Spatial variation in biodiversity is the result of complex interactions between evolutionary history and ecological factors. Methods in historical biogeography combine phylogenetic information with current species locations to infer the evolutionary history of a clade through space and time. A major limitation of most methods for historical biogeographic inference is the requirement of single locations for terminal lineages, reducing contemporary species geographical ranges to a point in two-dimensional space. In reality, geographic ranges usually show complex geographic patterns, irregular shapes, or discontinuities. In this article, we describe a method for phylogeographic analysis using polygonal species geographic ranges of arbitrary complexity. By integrating the geographic diversification process across species ranges, we provide a method to infer the geographic location of ancestors in a Bayesian framework. By modeling migration conditioned on a phylogenetic tree, this approach permits reconstructing the geographic location of ancestors through time. We apply this new method to the diversification of two neotropical bird genera, Trumpeters (Psophia) and Cinclodes ovenbirds. We demonstrate the usefulness of our method (called rase) in phylogeographic reconstruction of species ancestral locations and contrast our results with previous methods that compel researchers to reduce the distribution of species to one point in space. We discuss model extensions to enable a more general, spatially explicit framework for historical biogeographic analysis.

Dealing with discordant genetic signal caused by hybridisation, incomplete lineage sorting and paucity of primary nucleotide homologies: a case study of closely related members of the genus Picris subsection Hieracioides (Compositae)

We investigated genetic variation and evolutionary history of closely related taxa of Picris subsect. Hieracioides with major focus on the widely distributed P. hieracioides and its closely related congeners, P. hispidissima, P. japonica, P. olympica, and P. nuristanica. Accessions from 140 sample sites of the investigated Picris taxa were analyzed on the infra- and the inter-specific level using nuclear (ITS1-5.8S-ITS2 region) and chloroplast (rpl32-trnL^(UAG) region) DNA sequences. Genetic patterns of P. hieracioides, P. hispidissima, and P. olympica were shown to be incongruent and, in several cases, both plastid and nuclear alleles transcended borders of the taxa and genetic lineages. The widespread P. hieracioides was genetically highly variable and non-monophyletic across both markers, with allele groups having particular geographic distributions. Generally, all gene trees and networks displayed only a limited and statistically rather unsupported resolution among ingroup taxa causing their phylogenetic relationships to remain rather unresolved. More light on these intricate evolutionary relationships was cast by the Bayesian coalescent-based analysis, although some relationships were still left unresolved. A combination of suite of phylogenetic analyses revealed the ingroup taxa to represent a complex of genetically closely related and morphologically similar entities that have undergone a highly dynamic and recent evolution. This has been especially affected by the extensive and recurrent gene flow among and within the studied taxa and/or by the maintenance of ancestral variation. Paucity of phylogenetically informative signal further hampers the reconstruction of relationships on the infra- as well as on the inter-specific level. In the present study, we have demonstrated that a combination of various phylogenetic analyses of datasets with extremely complex and incongruent phylogenetic signal may shed more light on the interrelationships and evolutionary history of analysed species groups.

Ecological studies of polyploidy in the 100 years following its discovery

Polyploidy is a mutation with profound phenotypic consequences and thus hypothesized to have transformative effects in plant ecology. This is most often considered in the context of geographical and environmental distributions-as achieved from divergence of physiological and life-history traits-but may also include species interactions and biological invasion. This paper presents a historical overview of hypotheses and empirical data regarding the ecology of polyploids. Early researchers of polyploidy (1910 s-1930 s) were geneticists by training but nonetheless savvy to its phenotypic effects, and speculated on the importance of genome duplication to adaptation and crop improvement. Cytogenetic studies in the 1930 s-1950 s indicated that polyploids are larger (sturdier foliage, thicker stems and taller stature) than diploids while cytogeographic surveys suggested that polyploids and diploids have allopatric or parapatric distributions. Although autopolyploidy was initially regarded as common, influential writings by North American botanists in the 1940 s and 1950 s argued for the principle role of allopolyploidy; according to this view, genome duplication was significant for providing a broader canvas for hybridization rather than for its phenotypic effects per se. The emphasis on allopolyploidy had a chilling effect on nascent ecological work, in part due to taxonomic challenges posed by interspecific hybridization. Nonetheless, biosystematic efforts over the next few decades (1950s-1970s) laid the foundation for ecological research by documenting cytotype distributions and identifying phenotypic correlates of polyploidy. Rigorous investigation of polyploid ecology was achieved in the 1980s and 1990 s by population biologists who leveraged flow cytometry for comparative work in autopolyploid complexes. These efforts revealed multi-faceted ecological and phenotypic differences, some of which may be direct consequences of genome duplication. Several classical hypotheses about the ecology of polyploids remain untested, however, and allopolyploidy--regarded by most botanists as the primary mode of genome duplication--is largely unstudied in an ecological context.

Relationship of genetic diversity and niche centrality: a survey and analysis

The distribution of genetic diversity within and among populations in relation to species' geographic ranges is important to understanding processes of evolution, speciation, and biogeography. One hypothesis predicts that natural populations at geographic range margins will have lower genetic diversity relative to those located centrally in species' distributions owing to a link between geographic and environmental marginality; alternatively, genetic variation may be unrelated with geographic marginality via decoupling of geographic and environmental marginality. We investigate the predictivity of geographic patterns of genetic variation based on geographic and environmental marginality using published genetic diversity data for 40 species (insects, plants, birds, mammals, worms). Only about half of species showed positive relationships between geographic and environmental marginality. Three analyses (sign test, multiple linear regression, and meta-analysis of correlation effect sizes) showed a negative relationship between genetic diversity and distance to environmental niche centroid, but no consistent relationship of genetic diversity with distance to geographic range center.

Species delimitation and morphological divergence in the scorpion Centruroides vittatus (Say, 1821): insights from phylogeography

Scorpion systematics and taxonomy have recently shown a need for revision, partially due to insights from molecular techniques. Scorpion taxonomy has been difficult with morphological characters as disagreement exists among researchers with character choice for adequate species delimitation in taxonomic studies. Within the family Buthidae, species identification and delimitation is particularly difficult due to the morphological similarity among species and extensive intraspecific morphological diversity. The genus Centruroides in the western hemisphere is a prime example of the difficulty in untangling the taxonomic complexity within buthid scorpions. In this paper, we present phylogeographic, Ecological Niche Modeling, and morphometric analyses to further understand how population diversification may have produced morphological diversity in Centruroides vittatus (Say, 1821). We show that C. vittatus populations in the Big Bend and Trans-Pecos region of Texas, USA are phylogeographically distinct and may predate the Last Glacial Maximum (LGM). In addition, we suggest the extended isolation of Big Bend region populations may have created the C. vittatus variant once known as C. pantheriensis.

Extensive range persistence in peripheral and interior refugia characterizes Pleistocene range dynamics in a widespread Alpine plant species (Senecio carniolicus, Asteraceae)

Recent evidence suggests that survival of arctic-alpine organisms in peripheral or interior glacial refugia are not mutually exclusive and may both be involved in shaping an organism’s Pleistocene history, yet potentially at different time levels. Here, we test this hypothesis in a high-mountain plant (diploid lineage of Senecio carniolicus, Asteraceae) from the Eastern European Alps, in which patterns of morphological variation and current habitat requirements suggest survival in both types of refugia. To this end, we used AFLPs, nuclear and plastid DNA sequences and analysed them, among others, within a graph theoretic framework and using novel Bayesian methods of phylogeographic inference. On the basis of patterns of genetic diversity, occurrence of rare markers, distribution of distinct genetic lineages and patterns of range connectivity both interior refugia in the formerly strongly glaciated central Alps and peripheral refugia along the southern margin of the Alps were identified. The presence of refugia congruently inferred by markers resolving at different time levels suggests that these refugia acted as such throughout several glacial cycles. The high degree of range persistence together with gradual range expansion, which contrasts with the extent of range shifts implied for other Alpine species, is likely responsible for incipient lineage differentiation evident from the genetic data. Replacing a simplistic peripheral vs. interior refugia dualism by more complex models involving both types of refugia and considering different time levels will help identifying common phylogeographic patterns with respect to, for instance, location of refugia and colonization routes and elucidating their underlying genetic and/or ecological causes.

Species delimitation and phylogeography of Aphonopelma hentzi (Araneae, Mygalomorphae, Theraphosidae): cryptic diversity in North American tarantulas

BACKGROUND

The primary objective of this study is to reconstruct the phylogeny of the hentzi species group and sister species in the North American tarantula genus, Aphonopelma, using a set of mitochondrial DNA markers that include the animal "barcoding gene". An mtDNA genealogy is used to consider questions regarding species boundary delimitation and to evaluate timing of divergence to infer historical biogeographic events that played a role in shaping the present-day diversity and distribution. We aimed to identify potential refugial locations, directionality of range expansion, and test whether A. hentzi post-glacial expansion fit a predicted time frame.

METHODS AND FINDINGS

A Bayesian phylogenetic approach was used to analyze a 2051 base pair (bp) mtDNA data matrix comprising aligned fragments of the gene regions CO1 (1165 bp) and ND1-16S (886 bp). Multiple species delimitation techniques (DNA tree-based methods, a "barcode gap" using percent of pairwise sequence divergence (uncorrected p-distances), and the GMYC method) consistently recognized a number of divergent and genealogically exclusive groups.

CONCLUSIONS

The use of numerous species delimitation methods, in concert, provide an effective approach to dissecting species boundaries in this spider group; as well they seem to provide strong evidence for a number of nominal, previously undiscovered, and cryptic species. Our data also indicate that Pleistocene habitat fragmentation and subsequent range expansion events may have shaped contemporary phylogeographic patterns of Aphonopelma diversity in the southwestern United States, particularly for the A. hentzi species group. These findings indicate that future species delimitation approaches need to be analyzed in context of a number of factors, such as the sampling distribution, loci used, biogeographic history, breadth of morphological variation, ecological factors, and behavioral data, to make truly integrative decisions about what constitutes an evolutionary lineage recognized as a "species".