Determining population structure and hybridization for two iris species

Complete plastid genome of Iris orchioides and comparative analysis with 19 Iris plastomes

Iris is a cosmopolitan genus comprising approximately 280 species distributed throughout the Northern Hemisphere. Although Iris is the most diverse group in the Iridaceae, the number of taxa is debatable owing to various taxonomic issues. Plastid genomes have been widely used for phylogenetic research in plants; however, only limited number of plastid DNA markers are available for phylogenetic study of the Iris. To understand the genomic features of plastids within the genus, including its structural and genetic variation, we newly sequenced and analyzed the complete plastid genome of I. orchioides and compared it with those of 19 other Iris taxa. Potential plastid markers for phylogenetic research were identified by computing the sequence divergence and phylogenetic informativeness. We then tested the utility of the markers with the phylogenies inferred from the markers and whole-plastome data. The average size of the plastid genome was 152,926 bp, and the overall genomic content and organization were nearly identical among the 20 Iris taxa, except for minor variations in the inverted repeats. We identified 10 highly informative regions (matK, ndhF, rpoC2, ycf1, ycf2, rps15-ycf, rpoB-trnC, petA-psbJ, ndhG-ndhI and psbK-trnQ) and inferred a phylogeny from each region individually, as well as from their concatenated data. Remarkably, the phylogeny reconstructed from the concatenated data comprising three selected regions (rpoC2, ycf1 and ycf2) exhibited the highest congruence with the phylogeny derived from the entire plastome dataset. The result suggests that this subset of data could serve as a viable alternative to the complete plastome data, especially for molecular diagnoses among closely related Iris taxa, and at a lower cost.

Revisiting the comparative phylogeography of unglaciated eastern North America: 15 years of patterns and progress

In a landmark comparative phylogeographic study, “Comparative phylogeography of unglaciated eastern North America,” Soltis et al. (Molecular Ecology, 2006, 15, 4261) identified geographic discontinuities in genetic variation shared across taxa occupying unglaciated eastern North America and proposed several common biogeographical discontinuities related to past climate fluctuations and geographic barriers. Since 2006, researchers have published many phylogeographical studies and achieved many advances in genotyping and analytical techniques; however, it is unknown how this work has changed our understanding of the factors shaping the phylogeography of eastern North American taxa. We analyzed 184 phylogeographical studies of eastern North American taxa published between 2007 and 2019 to evaluate: (1) the taxonomic focus of studies and whether a previously detected taxonomic bias towards studies focused on vertebrates has changed over time, (2) the extent to which studies have adopted genotyping technologies that improve the resolution of genetic groups (i.e., NGS DNA sequencing) and analytical approaches that facilitate hypothesis‐testing (i.e., divergence time estimation and niche modeling), and (3) whether new studies support the hypothesized biogeographic discontinuities proposed by Soltis et al. (Molecular Ecology, 2006, 15, 4261) or instead support new, previously undetected discontinuities. We observed little change in taxonomic focus over time, with studies still biased toward vertebrates. Although many technological and analytical advances became available during the period, uptake was slow and they were employed in only a small proportion of studies. We found variable support for previously identified discontinuities and identified one new recurrent discontinuity. However, the limited resolution and taxonomic breadth of many studies hindered our ability to clarify the most important climatological or geographical factors affecting taxa in the region. Broadening the taxonomic focus to include more non‐vertebrate taxa, employing technologies that improve genetic resolution, and using analytical approaches that improve hypothesis testing are necessary to strengthen our inference of the forces shaping the phylogeography of eastern North America.

Genetic structure in Louisiana Iris species reveals patterns of recent and historical admixture

PREMISE

When divergent lineages come into secondary contact, reproductive isolation may be incomplete, thus providing an opportunity to investigate how speciation is manifested in the genome. The Louisiana Irises (Iris, series Hexagonae) comprise a group of three or more ecologically and reproductively divergent lineages that can produce hybrids where they come into contact. We estimated standing genetic variation to understand the current distribution of population structure in the Louisiana Irises.

METHODS

We used genotyping-by-sequencing techniques to sample the genomes of Louisiana Iris species across their ranges. We sampled 20 populations (n = 632 individuals) across 11,249 loci and used Entropy and PCA models to assess population genetic data.

RESULTS

We discovered evidence for interspecific gene flow in parts of the range. Our analysis revealed patterns of population structure at odds with widely accepted nominal taxonomy. We discovered undescribed hybrid populations, designated as belonging to the I. brevicaulis lineage. Iris nelsonii shared significant ancestry with only one of the purported parent species, I. fulva, evidence inconsistent with a hybrid origin.

CONCLUSIONS

This study provides several key findings important to the investigation of standing genetic variation in the Louisiana Iris species complex. Compared to the other nominal species, I. brevicaulis contains a large amount of genetic diversity. In addition, we discovered a previously unknown hybrid zone between I. brevicaulis and I. hexagona along the Texas coast. Finally, our results do not support the long-standing hypothesis that I. nelsonii has mixed ancestry from three parental taxa.

Molecular signatures of long-distance oceanic dispersal and the colonization of Pacific islands in Lycium carolinianum

PREMISE

Long-distance dispersal has been important in explaining the present distributions of many plant species. Despite being infrequent, such dispersal events have considerable evolutionary consequences, because bottlenecks during colonization can result in reduced genetic diversity. We examined the phylogeographic history of Lycium carolinianum, a widespread taxon that ranges from southeastern North America to several Pacific islands, with intraspecific diversity in sexual and mating systems.

METHODS

We used Bayesian, likelihood, and coalescent approaches with nuclear and plastid sequence data and genome-wide single nucleotide polymorphisms to reconstruct the dispersal history of this species. We also compared patterns of genetic variation in mainland and island populations using single nucleotide polymorphisms and allelic diversity at the S-RNase mating system gene.

RESULTS

Lycium carolinianum is monophyletic and dispersed once from the North American mainland, colonizing the Pacific islands ca. 40,100 years ago. This dispersal was accompanied by a loss of genetic diversity in SNPs and the S-RNase locus due to a colonization bottleneck and the loss of self-incompatibility. Additionally, we documented at least two independent transitions to gynodioecy: once following the colonization of the Hawaiian Islands and loss of self-incompatibility, and a second time associated with polyploidy in the Yucatán region of Mexico.

CONCLUSIONS

Long-distance dispersal via fleshy, bird dispersed fruits best explains the unusually widespread distribution of L. carolinianum. The collapse of diversity at the S-RNase locus in island populations suggests that self-fertilization may have facilitated the subsequent colonization of Pacific islands following a single dispersal from mainland North America.

From pollen dispersal to plant diversification: genetic consequences of pollination mode

Pollinators influence patterns of plant speciation, and one intuitive hypothesis is that pollinators affect rates of plant diversification through their effects on pollen dispersal. By specifying mating events and pollen flow across the landscape, distinct types of pollinators may cause different opportunities for allopatric speciation. This pollen dispersal-dependent speciation hypothesis predicts that pollination mode has effects on the spatial context of mating events that scale up to impact population structure and rates of species formation. Here I consider recent comparative studies, including genetic analyses of plant mating events, population structure and comparative phylogenetic analyses, to examine evidence for this model. These studies suggest that highly mobile pollinators conduct greater gene flow within and among populations, compared to less mobile pollinators. These differences influence patterns of population structure across the landscape. However, the effects of pollination mode on speciation rates is less predictable. In some contexts, the predicted effects of pollen dispersal are outweighed by other factors that govern speciation rates. A multiscale approach to examine effects of pollination mode on plant mating system, population structure and rates of diversification is key to determining the role of pollen dispersal on plant speciation for model clades.

Population Genetics of Odontarrhena (Brassicaceae) from Albania: The Effects of Anthropic Habitat Disturbance, Soil, and Altitude on a Ni-Hyperaccumulator Plant Group from a Major Serpentine Hotspot

Albanian taxa and populations of the genus Odontarrhena are most promising candidates for research on metal tolerance and Ni-agromining, but their genetic structure remains unknown. We investigated phylogenetic relationships and genetic differentiation in relation to distribution and ploidy of the taxa, anthropic site disturbance, elevation, soil type, and trace metals at each population site. After performing DNA sequencing of selected accessions, we applied DNA-fingerprinting to analyze the genetic structure of 32 populations from ultramafic and non-ultramafic outcrops across Albania. Low sequence divergence resulted in poorly resolved phylograms, but supported affinity between the two diploid serpentine endemics O. moravensis and O. rigida. Analysis of molecular variance (AMOVA) revealed significant population differentiation, but no isolation by distance. Among-population variation was higher in polyploids than in diploids, in which genetic distances were lower. Genetic admixing at population and individual level occurred especially in the polyploids O. chalcidica, O. decipiens, and O. smolikana. Admixing increased with site disturbance. Outlier loci were higher in serpentine populations but decreased along altitude with lower drought and heat stress. Genetic variability gained by gene flow and hybridization at contact zones with "resident" species of primary ultramafic habitats promoted expansion of the tetraploid O. chalcidica across anthropogenic sites.

Morphological and genome-wide evidence for natural hybridisation within the genus Stipa (Poaceae)

Hybridisation in the wild between closely related species is a common mechanism of speciation in the plant kingdom and, in particular, in the grass family. Here we explore the potential for natural hybridisation in Stipa (one of the largest genera in Poaceae) between genetically distant species at their distribution edges in Mountains of Central Asia using integrative taxonomy. Our research highlights the applicability of classical morphological and genome reduction approaches in studies on wild plant species. The obtained results revealed a new nothospecies, Stipa × lazkovii, which exhibits intermediate characters to S. krylovii and S. bungeana. A high-density DArTseq assay disclosed that S. × lazkovii is an F1 hybrid, and established that the plastid and mitochondrial DNA was inherited from S. bungeana. In addition, molecular markers detected a hybridisation event between morphologically and genetically distant species S. bungeana and probably S. glareosa. Moreover, our findings demonstrated an uncertainty on the taxonomic status of S. bungeana that currently belongs to the section Leiostipa, but it is genetically closer to S. breviflora from the section Barbatae. Finally, we noticed a discrepancy between the current molecular data with the previous findings on S. capillata and S. sareptana.

Biological action in Read-Write genome evolution

Many of the most important evolutionary variations that generated phenotypic adaptations and originated novel taxa resulted from complex cellular activities affecting genome content and expression. These activities included (i) the symbiogenetic cell merger that produced the mitochondrion-bearing ancestor of all extant eukaryotes, (ii) symbiogenetic cell mergers that produced chloroplast-bearing ancestors of photosynthetic eukaryotes, and (iii) interspecific hybridizations and genome doublings that generated new species and adaptive radiations of higher plants and animals. Adaptive variations also involved horizontal DNA transfers and natural genetic engineering by mobile DNA elements to rewire regulatory networks, such as those essential to viviparous reproduction in mammals. In the most highly evolved multicellular organisms, biological complexity scales with 'non-coding' DNA content rather than with protein-coding capacity in the genome. Coincidentally, 'non-coding' RNAs rich in repetitive mobile DNA sequences function as key regulators of complex adaptive phenotypes, such as stem cell pluripotency. The intersections of cell fusion activities, horizontal DNA transfers and natural genetic engineering of Read-Write genomes provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

Targeted capture to assess neutral genomic variation in the narrow-leaf hopbush across a continental biodiversity refugium

The Adelaide geosyncline, a mountainous region in central southern Australia, is purported to be an important continental refugium for Mediterranean and semi-arid Australian biota, yet few population genetic studies have been conducted to test this theory. Here, we focus on a plant species distributed widely throughout the region, the narrow-leaf hopbush, Dodonaea viscosa ssp. angustissima, and examine its genetic diversity and population structure. We used a hybrid-capture target enrichment technique to selectively sequence over 700 genes from 89 individuals across 17 sampling locations. We compared 815 single nucleotide polymorphisms among individuals and populations to investigate population genetic structure. Three distinct genetic clusters were identified; a Flinders/Gammon ranges cluster, an Eastern cluster, and a Kangaroo Island cluster. Higher genetic diversity was identified in the Flinders/Gammon Ranges cluster, indicating that this area is likely to have acted as a refugium during past climate oscillations. We discuss these findings and consider the historical range dynamics of these populations. We also provide methodological considerations for population genomics studies that aim to use novel genomic approaches (such as target capture methods) on non-model systems. The application of our findings to restoration of this species across the region are also considered.

Genotyping-by-Sequencing in a Species Complex of Australian Hummock Grasses (Triodia): Methodological Insights and Phylogenetic Resolution

Next-generation sequencing is becoming increasingly accessible to researchers asking biosystematic questions, but current best practice in both choosing a specific approach and effectively analysing the resulting data set is still being explored. We present a case study for the use of genotyping-by-sequencing (GBS) to resolve relationships in a species complex of Australian arid and semi-arid grasses (Triodia R.Br.), highlighting our solutions to methodological challenges in the use of GBS data. We merged overlapping paired-end reads then optimised locus assembly in the program PyRAD to generate GBS data sets for phylogenetic and distance-based analyses. In addition to traditional concatenation analyses in RAxML, we also demonstrate the novel use of summary species tree analyses (taking gene trees as input) with GBS loci. We found that while species tree analyses were relatively robust to variation in PyRAD assembly parameters, our RAxML analyses resulted in well-supported but conflicting topologies under different assembly settings. Despite this conflict, multiple clades in the complex were consistently supported as distinct across analyses. Our GBS data assembly and analyses improve the resolution of taxa and phylogenetic relationships in the Triodia basedowii complex compared to our previous study based on Sanger sequencing of nuclear (ITS/ETS) and chloroplast (rps16-trnK spacer) markers. The genomic results also partly support previous evidence for hybridization between species in the complex. Our methodological insights for analysing GBS data will assist researchers using similar data to resolve phylogenetic relationships within species complexes.

Niche Conservatism for Ecological Preference in the Louisiana iris species complex

Spatial and temporal environmental variation influences evolutionary processes such as divergence among populations and species. Here, we investigate the patterns of niche evolution for the Louisiana irises as well as understanding the phylogenetic relationship between these species. Using BEAST, a species phylogeny was generated for the Louisiana irises in order to test the hypothesis of whether niche conservatism has played an important role for this species complex. Species Distribution Models were constructed for present day distributions to determine the environmental factors, which contribute to species ranges. Tests of niche similarity were performed in order to evaluate if niche conservatism is apparent within this species complex. We demonstrate that the Louisiana iris species complex is a monophyletic clade with I. brevicaulis and I. fulva as being sister to each other. The differences observed among the iris species, in regard to associated environmental factors suggest an effect from these components on the distributions and habitats occupied. Furthermore, tests of niche similarity indicate niche conservatism for all species comparisons. Working at the species level and assessing various factors that can influence differentiation, this study assessed a more complete picture of the ecological and evolutionary history of this species complex across their geographic and ecological range.

Single nucleotide polymorphism discovery via genotyping by sequencing to assess population genetic structure and recurrent polyploidization in Andropogon gerardii

PREMISE OF THE STUDY

Autopolyploidy, genome duplication within a single lineage, can result in multiple cytotypes within a species. Geographic distributions of cytotypes may reflect the evolutionary history of autopolyploid formation and subsequent population dynamics including stochastic (drift) and deterministic (differential selection among cytotypes) processes. Here, we used a population genomic approach to investigate whether autopolyploidy occurred once or multiple times in Andropogon gerardii, a widespread, North American grass with two predominant cytotypes.

METHODS

Genotyping by sequencing was used to identify single nucleotide polymorphisms (SNPs) in individuals collected from across the geographic range of A. gerardii. Two independent approaches to SNP calling were used: the reference-free UNEAK pipeline and a reference-guided approach based on the sequenced Sorghum bicolor genome. SNPs generated using these pipelines were analyzed independently with genetic distance and clustering.

KEY RESULTS

Analyses of the two SNP data sets showed very similar patterns of population-level clustering of A. gerardii individuals: a cluster of A. gerardii individuals from the southern Plains, a northern Plains cluster, and a western cluster. Groupings of individuals corresponded to geographic localities regardless of cytotype: 6x and 9x individuals from the same geographic area clustered together.

CONCLUSIONS

SNPs generated using reference-guided and reference-free pipelines in A. gerardii yielded unique subsets of genomic data. Both data sets suggest that the 9x cytotype in A. gerardii likely evolved multiple times from 6x progenitors across the range of the species. Genomic approaches like GBS and diverse bioinformatics pipelines used here facilitate evolutionary analyses of complex systems with multiple ploidy levels.

GIbPSs: a toolkit for fast and accurate analyses of genotyping-by-sequencing data without a reference genome

Genotyping-by-sequencing (GBS) and related methods are increasingly used for studies of non-model organisms from population genetic to phylogenetic scales. We present GIbPSs, a new genotyping toolkit for the analysis of data from various protocols such as RAD, double-digest RAD, GBS, and two-enzyme GBS without a reference genome. GIbPSs can handle paired-end GBS data and is able to assign reads from both strands of a restriction fragment to the same locus. GIbPSs is most suitable for population genetic and phylogeographic analyses. It avoids genotyping errors due to indel variation by identifying and discarding affected loci. GIbPSs creates a genotype database that offers rich functionality for data filtering and export in numerous formats. We performed comparative analyses of simulated and real GBS data with GIbPSs and another program, pyRAD. This program accounts for indel variation by aligning homologous sequences. GIbPSs performed better than pyRAD in several aspects. It required much less computation time and displayed higher genotyping accuracy. GIbPSs retained smaller numbers of loci overall in analyses of real GBS data. It nevertheless delivered more complete genotype matrices with greater locus overlap between individuals and greater numbers of loci sampled in all individuals.

Genetic differentiation and reduced genetic diversity at the northern range edge of two species with different dispersal modes

Theory predicts that genetic variation should be reduced at range margins, but empirical support is equivocal. Here, we used genotyping-by-sequencing technology to investigate genetic variation in central and marginal populations of two species in the marine gastropod genus Crepidula. These two species have different development and dispersal types and might therefore show different spatial patterns of genetic variation. Both allelic richness and the proportion of private alleles were highest in the most central populations of both species, and lower at the margin. The species with low dispersal, Crepidula convexa, showed high degrees of structure throughout the range that conform to the pattern found in previous studies using other molecular markers. The northernmost populations of the high-dispersing species, Crepidula fornicata, are distinct from more central populations, although this species has been previously observed to have little genetic structure over much of its range. Although genetic diversity was significantly lower at the range margin, the absolute reduction in diversity observed with these genomewide markers was slight, and it is not yet known whether there are functional consequences for the marginal populations.

Multiple hydroxyphenethyl glucosinolate isomers and their tandem mass spectrometric distinction in a geographically structured polymorphism in the crucifer Barbarea vulgaris

Two distinct glucosinolate (GSL) chemotypes (P and G-types) of Barbarea vulgaris (Brassicaceae) were known from southern Scandinavia, but whether the types were consistent in a wider geographic area was not known. Populations (26) from Eastern and Central Europe were analyzed for GSLs in order to investigate whether the two types were consistent in this area. Most (21) could be attributed to one of the previously described GSL profile types, the P-type (13 populations) and the G-type (8 populations), based on differences in the stereochemistry of 2-hydroxylation, presence or absence of phenolic glucobarbarin derivatives, and qualitative differences in indole GSL decoration (tested for a subset of 8+6 populations only). The distinction agreed with previous molecular genetic analysis of the same individuals. Geographically, the P-type typically occurred in Eastern Europe while the G-type mainly occurred in Central Europe. Of the remaining five populations, minor deviations were observed in some individuals from two populations genetically assigned to the G-type, and a hybrid population from Finland contained an additional dihydroxyphenethyl GSL isomer attributed to a combinatorial effect of P-type and G-type genes. Major exceptions to the typical GSL profiles were observed in two populations: (1) A G-type population from Slovenia deviated by a high frequency of a known variant in glucobarbarin biosynthesis ('NAS form') co-occurring with usual G-type individuals. (2) A population from Caucasus exhibited a highly deviating GSL profile dominated by p-hydroxyphenethyl GSL that was insignificant in other accessions, as well as two GSLs investigated by NMR, m-hydroxyphenethylGSL and a partially identified m,p disubstituted hydroxy-methoxy derivative of phenethylGSL. Tandem HPLC-MS of seven NMR-identified desulfoGSLs was carried out and interpreted for increased certainty in peak identification and as a tool for partial structure elucidation. The distinct, geographically separated chemotypes and rare variants are discussed in relation to future taxonomic revision and the genetics and ecology of GSLs in B. vulgaris.