The population demography of Betula maximowicziana,a cool‐temperate tree species in Japan,in relation to the last glacial period: its admixture‐like genetic structure is the result of simple population splitting not admixing

Conservation of the local genetic variation and evolutionary integrity of economically and ecologically important trees is a key aspect of studies involving forest genetics, and a population demographic history of the target species provides valuable information for this purpose. Here, the genetic structure of 48 populations of Betula maximowicziana was assessed using 12 expressed sequence tag–simple sequence repeat (EST‐SSR) markers. Genetic diversity was lower in northern populations than southern ones and structure analysis revealed three groups: northern and southern clusters and an admixed group. Eleven more genomic‐SSR loci were added and the demographic history of these three groups was inferred by approximate Bayesian computation (ABC). The ABC revealed that a simple split scenario was much more likely than isolation with admixture, suggesting that the admixture‐like structure detected in this species was due to ancestral polymorphisms. The ABC analysis suggested that the population growth and divergence of the three groups occurred 96 800 (95% CI, 20 500–599 000) and 28 300 (95% CI, 8700–98 400) years ago, respectively. We need to be aware of several sources of uncertainty in the inference such as assumptions about the generation time, overlapping of generations, confidence intervals of the estimated parameters and the assumed model in the ABC. However, the results of the ABC together with the model‐based maps of reconstructed past species distribution and palaeoecological data suggested that the modern genetic structure of B. maximowicziana originated prior to the last glacial maximum (LGM) and that some populations survived in the northern range even during the LGM.     

Linkage disequilibrium network analysis (LDna) gives a global view of chromosomal inversions,local adaptation and geographic structure

Recent advances in sequencing allow population‐genomic data to be generated for virtually any species. However, approaches to analyse such data lag behind the ability to generate it, particularly in nonmodel species. Linkage disequilibrium (LD, the nonrandom association of alleles from different loci) is a highly sensitive indicator of many evolutionary phenomena including chromosomal inversions, local adaptation and geographical structure. Here, we present linkage disequilibrium network analysis (LDna), which accesses information on LD shared between multiple loci genomewide. In LD networks, vertices represent loci, and connections between vertices represent the LD between them. We analysed such networks in two test cases: a new restriction‐site‐associated DNA sequence (RAD‐seq) data set for Anopheles baimaii, a Southeast Asian malaria vector; and a well‐characterized single nucleotide polymorphism (SNP) data set from 21 three‐spined stickleback individuals. In each case, we readily identified five distinct LD network clusters (single‐outlier clusters, SOCs), each comprising many loci connected by high LD. In A. baimaii, further population‐genetic analyses supported the inference that each SOC corresponds to a large inversion, consistent with previous cytological studies. For sticklebacks, we inferred that each SOC was associated with a distinct evolutionary phenomenon: two chromosomal inversions, local adaptation, population‐demographic history and geographic structure. LDna is thus a useful exploratory tool, able to give a global overview of LD associated with diverse evolutionary phenomena and identify loci potentially involved. LDna does not require a linkage map or reference genome, so it is applicable to any population‐genomic data set, making it especially valuable for nonmodel species.     

Genomic signature of successful colonization of Eurasia by the allopolyploid shepherd's purse (Capsella bursa‐pastoris)

Polyploidization is a dominant feature of flowering plant evolution. However, detailed genomic analyses of the interpopulation diversification of polyploids following genome duplication are still in their infancy, mainly because of methodological limits, both in terms of sequencing and computational analyses. The shepherd's purse (Capsella bursa‐pastoris) is one of the most common weed species in the world. It is highly self‐fertilizing, and recent genomic data indicate that it is an allopolyploid, resulting from hybridization between the ancestors of the diploid species Capsella grandiflora and Capsella orientalis. Here, we investigated the genomic diversity of C. bursa‐pastoris, its population structure and demographic history, following allopolyploidization in Eurasia. To that end, we genotyped 261 C. bursa‐pastoris accessions spread across Europe, the Middle East and Asia, using genotyping‐by‐sequencing, leading to a total of 4274 SNPs after quality control. Bayesian clustering analyses revealed three distinct genetic clusters in Eurasia: one cluster grouping samples from Western Europe and Southeastern Siberia, the second one centred on Eastern Asia and the third one in the Middle East. Approximate Bayesian computation (ABC) supported the hypothesis that C. bursa‐pastoris underwent a typical colonization history involving low gene flow among colonizing populations, likely starting from the Middle East towards Europe and followed by successive human‐mediated expansions into Eastern Asia. Altogether, these findings bring new insights into the recent multistage colonization history of the allotetraploid C. bursa‐pastoris and highlight ABC and genotyping‐by‐sequencing data as promising but still challenging tools to infer demographic histories of selfing allopolyploids.     

Phylogeographic inference using Bayesian model comparison across a fragmented chorus frog species complex

Fragmented species complexes provide an interesting system for investigating biogeographic history and the present distribution of genetic variation. Recent advances in sequencing technology and statistical phylogeography enable the collection and rigorous analysis of large multilocus data sets, but designing studies that produce meaningful phylogeographic inferences remains challenging. We implemented a Bayesian model comparison approach to investigate previous biogeographic hypotheses while simultaneously inferring the presence of genetic structure in a chorus frog species complex. The Illinois chorus frog (Pseudacris illinoensis), originally described as a subspecies of the broadly distributed Strecker's chorus frog (Pseudacris streckeri), occurs in small, disjunct regions associated with scarce sand prairie habitats that have been impacted by human development. We used high‐throughput sequencing to develop and collect a multitiered genetic data set comprised of three different marker types (23 anonymous nuclear sequence loci, four mitochondrial genes and 14 microsatellite loci) designed to address questions across different evolutionary timescales. Phylogenetic analyses uncovered a deep divergence between populations in the Edwards Plateau of central Texas and all other P. streckeri/P. illinoensis populations, but suggest the disjunct distribution of P. illinoensis occurred more recently. Our best‐supported migration model is consistent with the hypothesis that central Texas represented a refugium from which populations expanded via multiple routes. This model also indicates that disjunct northern and southern regions of P. illinoensis should be considered genetically distinct management units. Our study provides an evolutionary context for future studies and conservation efforts in P. illinoensis and demonstrates the utility of model‐based approaches for phylogeographic inference.     

Sequence capture using AFLP‐generated baits: A cost‐effective method for high‐throughput phylogenetic and phylogeographic analysis

Target sequence capture is an efficient technique to enrich specific genomic regions for high‐throughput sequencing in ecological and evolutionary studies. In recent years, many sequence capture approaches have been proposed, but most of them rely on commercial synthetic baits which make the experiment expensive. Here, we present a novel sequence capture approach called AFLP‐based genome sequence capture (AFLP Capture). This method uses the AFLP (amplified fragment length polymorphism) technique to generate homemade capture baits without the need for prior genome information, thus is applicable to any organisms. In this approach, biotinylated AFLP fragments representing a random fraction of the genome are used as baits to capture the homologous fragments from genomic shotgun sequencing libraries. In a trial study, by using AFLP Capture, we successfully obtained 511 orthologous loci (>700,000 bp in total length) from 11 Odorrana species and more than 100,000 single nucleotide polymorphisms (SNPs) in four analyzed individuals of an Odorrana species. This result shows that our method can be used to address questions of various evolutionary depths (from interspecies level to intraspecies level). We also discuss the flexibility in bait preparation and how the sequencing data are analyzed. In summary, AFLP Capture is a rapid and flexible tool and can significantly reduce the experimental cost for phylogenetic studies that require analyzing genome‐scale data (hundreds or thousands of loci).     

Identification of SNP markers for inferring phylogeny in temperate bamboos (Poaceae: Bambusoideae) using RAD sequencing

Phylogenetic relationships among temperate species of bamboo are difficult to resolve, owing to both the challenge of detecting sufficiently variable markers and their polyploid history. Here, we use restriction site–associated DNA sequencing to identify candidate loci with fixed allelic differences segregating between and within two temperate species of bamboos: Arundinaria faberi and Yushania brevipaniculata. Approximately 27 million paired‐end sequencing reads were generated across four samples. From pooled data, we assembled 67 685 and 70 668 de novo contigs from partial overlap among paired‐end reads, with an average length of 240 and 241 bp for the two species, respectively, which were used to investigate functional classification of RAD tags in a blastx search. Analysed separately by population, we recovered 29 443 putatively orthologous RAD tags shared across the four sampled populations, containing 28 023 sequence variants, of which c. 13 000 are segregating between species, and c. 3000 segregating between populations within each species. Analyses based on these RAD tags yielded robust phylogenetic inferences, even with data set constructed from surprisingly few loci. This study illustrates the potential for reduced‐representation genome data to resolve difficult phylogenetic relationships in temperate bamboos.     

Genomic data detect corresponding signatures of population size change on an ecological time scale in two salamander species

Understanding the demography of species over recent history (e.g. <100 years) is critical in studies of ecology and evolution, but records of population history are rarely available. Surveying genetic variation is a potential alternative to census‐based estimates of population size, and can yield insight into the demography of a population. However, to assess the performance of genetic methods, it is important to compare their estimates of population history to known demography. Here, we leveraged the exceptional resources from a wetland with 37 years of amphibian mark–recapture data to study the utility of genetically based demographic inference on salamander species with documented population declines (Ambystoma talpoideum) and expansions (A. opacum), patterns that have been shown to be correlated with changes in wetland hydroperiod. We generated ddRAD data from two temporally sampled populations of A. opacum (1993, 2013) and A. talpoideum (1984, 2011) and used coalescent‐based demographic inference to compare alternate evolutionary models. For both species, demographic model inference supported population size changes that corroborated mark–recapture data. Parameter estimation in A. talpoideum was robust to our variations in analytical approach, while estimates for A. opacum were highly inconsistent, tempering our confidence in detecting a demographic trend in this species. Overall, our robust results in A. talpoideum suggest that genome‐based demographic inference has utility on an ecological scale, but researchers should also be cognizant that these methods may not work in all systems and evolutionary scenarios. Demographic inference may be an important tool for population monitoring and conservation management planning.     

A high frequency of allopolyploid speciation in the gymnospermous genus Ephedra and its possible association with some biological and ecological features

The origin and evolution of polyploids have been studied extensively in angiosperms and ferns but very rarely in gymnosperms. With the exception of three species of conifers, all natural polyploid species of gymnosperms belong to Ephedra, in which more than half of the species show polyploid cytotypes. Here, we investigated the origin and evolution of polyploids of Ephedra distributed in the Qinghai–Tibetan Plateau (QTP) and neighbouring areas. Flow cytometry (FCM) was used to measure the ploidy levels of the sampled species that are represented by multiple individuals from different populations, and then, two single‐copy nuclear genes (LFY and DDB2) and two chloroplast DNA fragments were used to unravel the possible origins and maternal donors of the polyploids. The results indicate that the studied polyploid species are allopolyploids, and suggest that allotetraploidy is a dominant mode of speciation in Ephedra. The high percentage of polyploids in the genus could be related to some of its biological attributes such as vegetative propagation, a relatively high rate of unreduced gamete formation, and a small genome size relative to most other gymnosperms. Significant ecological divergences between allotetraploids and their putative progenitors were detected by PCAs and anova and Tukey's tests, with the exception of E. saxatilis. The overlap of geographical distributions and ecological niches of some diploid species could have provided opportunities for interspecific hybridization and allopolyploid speciation.     

What remains from a 454 run: estimation of success rates of microsatellite loci development in selected newt species (Calotriton asper,Lissotriton helveticus,and Triturus cristatus) and comparison with Illumina‐based approaches

The development of microsatellite loci has become more efficient using next‐generation sequencing (NGS) approaches, and many studies imply that the amount of applicable loci is large. However, few studies have sought to quantify the number of loci that are retained for use out of the thousands of sequence reads initially obtained. We analyzed the success rate of microsatellite loci development for three amphibian species using a 454 NGS approach on tetra‐nucleotide motif‐enriched species‐specific libraries. The number of sequence reads obtained differed strongly between species and ranged from 19,562 for Triturus cristatus to 55,626 for Lissotriton helveticus, with 52,075 reads obtained for Calotriton asper. PHOBOS was used to identify sequences with tetra‐nucleotide repeat motifs with a minimum repeat number of ten and high quality primer binding sites. Of 107 sequences for T. cristatus, 316 for C. asper and 319 for L. helveticus, we tested the amplification success, polymorphism, and degree of heterozygosity for 41 primer combinations each for C. asper and T. cristatus, and 22 for L. helveticus. We found 11 polymorphic loci for T. cristatus, 20 loci for C. asper, and 15 loci for L. helveticus. Extrapolated, the number of potentially amplifiable loci (PALs) resulted in estimated species‐specific success rates of 0.15% (T. cristatus), 0.30% (C. asper), and 0.39% (L. helveticus). Compared with representative Illumina NGS approaches, our applied 454‐sequencing approach on specifically enriched sublibraries proved to be quite competitive in terms of success rates and number of finally applicable loci.     

ABC inference of multi‐population divergence with admixture from unphased population genomic data

Rapidly developing sequencing technologies and declining costs have made it possible to collect genome‐scale data from population‐level samples in nonmodel systems. Inferential tools for historical demography given these data sets are, at present, underdeveloped. In particular, approximate Bayesian computation (ABC) has yet to be widely embraced by researchers generating these data. Here, we demonstrate the promise of ABC for analysis of the large data sets that are now attainable from nonmodel taxa through current genomic sequencing technologies. We develop and test an ABC framework for model selection and parameter estimation, given histories of three‐population divergence with admixture. We then explore different sampling regimes to illustrate how sampling more loci, longer loci or more individuals affects the quality of model selection and parameter estimation in this ABC framework. Our results show that inferences improved substantially with increases in the number and/or length of sequenced loci, while less benefit was gained by sampling large numbers of individuals. Optimal sampling strategies given our inferential models included at least 2000 loci, each approximately 2 kb in length, sampled from five diploid individuals per population, although specific strategies are model and question dependent. We tested our ABC approach through simulation‐based cross‐validations and illustrate its application using previously analysed data from the oak gall wasp, Biorhiza pallida.     

Haploid,diploid, and pooled exome capture recapitulate features of biology and paralogy in two non-model tree species

Despite their suitability for studying evolution, many conifer species have large and repetitive giga-genomes (16–31 Gbp) that create hurdles to producing high coverage SNP data sets that capture diversity from across the entirety of the genome. Due in part to multiple ancient whole genome duplication events, gene family expansion and subsequent evolution within Pinaceae, false diversity from the misalignment of paralog copies creates further challenges in accurately and reproducibly inferring evolutionary history from sequence data. Here, we leverage the cost-saving benefits of pool-seq and exome-capture to discover SNPs in two conifer species, Douglas-fir (Pseudotsuga menziesii var. menziesii (Mirb.) Franco, Pinaceae) and jack pine (Pinus banksiana Lamb., Pinaceae). We show, using minimal baseline filtering, that allele frequencies estimated from pooled individuals show a strong, positive correlation with those estimated by sequencing the same population as individuals (r > .948), on par with such comparisons made in model organisms. Further, we highlight the utility of haploid megagametophyte tissue for identifying sites that are probably due to misaligned paralogs. Together with additional minor filtering, we show that it is possible to remove many of the loci with large frequency estimate discrepancies between individual and pooled sequencing approaches, improving the correlation further (r > .973). Our work addresses bioinformatic challenges in non-model organisms with large and complex genomes, highlights the use of megagametophyte tissue for the identification of paralogous artefacts, and suggests the combination of pool-seq and exome capture to be robust for further evolutionary hypothesis testing in these systems.     

Spatio‐temporal patterns of genome evolution in allotetraploid species of the genus Oryza

Despite knowledge that polyploidy is widespread and a major evolutionary force in flowering plant diversification, detailed comparative molecular studies on polyploidy have been confined to only a few species and families. The genus Oryza is composed of 23 species that are classified into ten distinct ‘genome types’ (six diploid and four polyploid), and is emerging as a powerful new model system to study polyploidy. Here we report the identification, sequence and comprehensive comparative annotation of eight homoeologous genomes from a single orthologous region (Adh1–Adh2) from four allopolyploid species representing each of the known Oryza genome types (BC, CD, HJ and KL). Detailed comparative phylogenomic analyses of these regions within and across species and ploidy levels provided several insights into the spatio‐temporal dynamics of genome organization and evolution of this region in ‘natural’ polyploids of Oryza. The major findings of this study are that: (i) homoeologous genomic regions within the same nucleus experience both independent and parallel evolution, (ii) differential lineage‐specific selection pressures do not occur between polyploids and their diploid progenitors, (iii) there have been no dramatic structural changes relative to the diploid ancestors, (iv) a variation in the molecular evolutionary rate exists between the two genomes in the BC complex species even though the BC and CD polyploid species appear to have arisen <2 million years ago, and (v) there are no clear distinctions in the patterns of genome evolution in the diploid versus polyploid species.     

Cyto‐nuclear discordance suggests complex evolutionary history in the cave‐dwelling salamander,Eurycea lucifuga

Our understanding of the evolutionary history and ecology of cave‐associated species has been driven historically by studies of morphologically adapted cave‐restricted species. Our understanding of the evolutionary history and ecology of nonrestricted cave species, troglophiles, is limited to a few studies, which present differing accounts of troglophiles’ relationship with the cave habitat, and its impact on population dynamics. Here, we used phylogenetics, demographic statistics, and population genetic methods to study lineage divergence, dates of divergence, and population structure in the Cave Salamander, Eurycea lucifuga, across its range. In order to perform these analyses, we sampled 233 individuals from 49 populations, using sequence data from three gene loci as well as genotyping data from 19 newly designed microsatellite markers. We find, as in many other species studied in a phylogeographic context, discordance between patterns inferred from mitochondrial relationships and those inferred by nuclear markers indicating a complicated evolutionary history in this species. Our results suggest Pleistocene‐based divergence among three main lineages within E. lucifuga corresponding to the western, central, and eastern regions of the range, similar to patterns seen in species separated in multiple refugia during climatic shifts. The conflict between mitochondrial and nuclear patterns is consistent with what we would expect from secondary contact between regional populations following expansion from multiple refugia.     

Congruent population structure across paralogous and nonparalogous loci in Salish Sea chum salmon (Oncorhynchus keta)

Whole‐genome duplications are major evolutionary events with a lasting impact on genome structure. Duplication events complicate genetic analyses as paralogous sequences are difficult to distinguish; consequently, paralogs are often excluded from studies. The effects of an ancient whole‐genome duplication (approximately 88 MYA) are still evident in salmonids through the persistence of numerous paralogous gene sequences and partial tetrasomic inheritance. We use restriction site‐associated DNA sequencing on 10 collections of chum salmon from the Salish Sea in the USA and Canada to investigate genetic diversity and population structure in both tetrasomic and rediploidized regions of the genome. We use a pedigree and high‐density linkage map to identify paralogous loci and to investigate genetic variation across the genome. By applying multivariate statistical methods, we show that it is possible to characterize paralogous loci and that they display similar patterns of population structure as the diploidized portion of the genome. We find genetic associations with the adaptively important trait of run‐timing in both sets of loci. By including paralogous loci in genome scans, we can observe evolutionary signals in genomic regions that have routinely been excluded from population genetic studies in other polyploid‐derived species.     

A novel method to infer the origin of polyploids from Amplified Fragment Length Polymorphism data reveals that the alpine polyploid complex of Senecio carniolicus (Asteraceae) evolved mainly via autopolyploidy

Despite its evolutionary and ecological relevance, the mode of polyploid origin has been notoriously difficult to be reconstructed from molecular data. Here, we present a method to identify the putative parents of polyploids and thus to infer the mode of their origin (auto‐ vs. allopolyploidy) from Amplified Fragment Length Polymorphism (AFLP) data. To this end, we use Cohen's d of distances between in silico polyploids, generated within a priori defined scenarios of origin from a priori delimited putative parental entities (e.g. taxa, genetic lineages), and natural polyploids. Simulations show that the discriminatory power of the proposed method increases mainly with increasing divergence between the lower‐ploid putative ancestors and less so with increasing delay of polyploidization relative to the time of divergence. We apply the new method to the Senecio carniolicus aggregate, distributed in the European Alps and comprising two diploid, one tetraploid and one hexaploid species. In the eastern part of its distribution, the S. carniolicus aggregate was inferred to comprise an autopolyploid series, whereas for western populations of the tetraploid species, an allopolyploid origin involving the two diploid species was the most likely scenario. Although this suggests that the tetraploid species has two independent origins, other evidence (ribotype distribution, morphology) is consistent with the hypothesis of an autopolyploid origin with subsequent introgression by the second diploid species. Altogether, identifying the best among alternative scenarios using Cohen's d can be straightforward, but particular scenarios, such as allopolyploid origin vs. autopolyploid origin with subsequent introgression, remain difficult to be distinguished.     

Genome‐wide sequence information reveals recurrent hybridization among diploid wheat wild relatives

Many conflicting hypotheses regarding the relationships among crops and wild species closely related to wheat (the genera Aegilops, Amblyopyrum, and Triticum) have been postulated. The contribution of hybridization to the evolution of these taxa is intensely discussed. To determine possible causes for this, and provide a phylogeny of the diploid taxa based on genome‐wide sequence information, independent data were obtained from genotyping‐by‐sequencing and a target‐enrichment experiment that returned 244 low‐copy nuclear loci. The data were analyzed using Bayesian, likelihood and coalescent‐based methods. D statistics were used to test if incomplete lineage sorting alone or together with hybridization is the source for incongruent gene trees. Here we present the phylogeny of all diploid species of the wheat wild relatives. We hypothesize that most of the wheat‐group species were shaped by a primordial homoploid hybrid speciation event involving the ancestral Triticum and Am. muticum lineages to form all other species except Ae. speltoides. This hybridization event was followed by multiple introgressions affecting all taxa except Triticum. Mostly progenitors of the extant species were involved in these processes, while recent interspecific gene flow seems insignificant. The composite nature of many genomes of wheat‐group taxa results in complicated patterns of diploid contributions when these lineages are involved in polyploid formation, which is, for example, the case for tetraploid and hexaploid wheats. Our analysis provides phylogenetic relationships and a testable hypothesis for the genome compositions in the basic evolutionary units within the wheat group of Triticeae.     

Is RAD‐seq suitable for phylogenetic inference? An in silico assessment and optimization

Inferring phylogenetic relationships between closely related taxa can be hindered by three factors: (1) the lack of informative molecular variation at short evolutionary timescale; (2) the lack of established markers in poorly studied taxa; and (3) the potential phylogenetic conflicts among different genomic regions due to incomplete lineage sorting or introgression. In this context, Restriction site Associated DNA sequencing (RAD‐seq) seems promising as this technique can generate sequence data from numerous DNA fragments scattered throughout the genome, from a large number of samples, and without preliminary knowledge on the taxa under study. However, divergence beyond the within‐species level will necessarily reduce the number of conserved and non‐duplicated restriction sites, and therefore the number of loci usable for phylogenetic inference. Here, we assess the suitability of RAD‐seq for phylogeny using a simulated experiment on the 12 Drosophila genomes, with divergence times ranging from 5 to 63 million years. These simulations show that RAD‐seq allows the recovery of the known Drosophila phylogeny with strong statistical support, even for relatively ancient nodes. Notably, this conclusion is robust to the potentially confounding effects of sequencing errors, heterozygosity, and low coverage. We further show that clustering RAD‐seq data using the BLASTN and SiLiX programs significantly improves the recovery of orthologous RAD loci compared with previously proposed approaches, especially for distantly related species. This study therefore validates the view that RAD sequencing is a powerful tool for phylogenetic inference.     

Molecular phylogeny and reticulate origins of the polyploid Bromus species from section Genea (Poaceae)

The origin of polyploid Bromus species of section Genea was investigated using molecular data. This group of annual species native from the Old-World is composed of three diploids, two tetraploids, one hexaploid, and one octoploid. Molecular cloning, sequencing, and phylogenetic analyses were performed on several accessions per species. We used the low copy nuclear gene Waxy, repeated rDNA spacers ITS1 and ITS2 and chloroplast spacers trnT-trnL and trnL-trnF. Our analyses revealed four different lineages involved in the parentage of the polyploids and confirmed their reticulate origin. Three of these lineages are closely related to the diploid species B. sterilis, B. tectorum, and B. fasciculatus. The fourth lineage could not be related to any diploid according to the available data. Our data gave insights on the origin of all the polyploids of section Genea, and chloroplast data allowed us to identify the maternal lineages. The Waxy gene was the most informative regarding origin of the polyploids. The Waxy copies duplicated by polyploidy appear selectively maintained in the polyploid species. No sequence heterogeneity was encountered in the ITS region, where concerted evolution seems to have occurred toward either maternal or paternal repeats. These results provide new information about the origin and molecular evolution of these polyploids and will allow a more accurate taxonomic treatment of the concerned species, based on their evolutionary history.