A next-generation sequencing method for overcoming the multiple gene copy problem in polyploid phylogenetics, applied to Poa grasses

Insight into the herbicide resistance patterns in Lolium rigidum populations in Tunisian and Moroccan wheat regions

Rigid ryegrass (Lolium rigidum Gaud.) is one of the most troublesome weeds in Moroccan and Tunisian cereal crop fields. In total, 19 rigid ryegrass field populations were randomly selected in northern wheat crop areas of Morocco and Tunisia to examine the patterns of herbicide resistance to acetolactate synthase (ALS)- and acetyl-CoA carboxylase (ACCase)-inhibiting herbicides. Greenhouse experiments confirmed reduced sensitivity to ALS- and/or ACCase-inhibiting herbicides in all L. rigidum populations. The occurrence of target-site resistance (TSR) was tested using high-throughput genotyping. The advent of next-generation sequencing (NGS) has enabled easy identification of causal mutations and confirmed the presence of ALS and ACCase mutations at specific codons conferring TSR. Thirteen populations showed resistance to ALS-inhibiting herbicides associated with point mutations in positions Pro-197-Thr, Pro-197-Ser, Pro-197-Leu, Pro-197-Gln and Trp-574-Leu, while resistance to ACCase-inhibiting herbicides was detected in 18 populations in positions Asp-2078-Val, Trp-2027-Cys, Ile-1781-Leu, Gly-2096-Ala, and Ile-2041-Asn of the enzymes conferring TSR. Additionally, dose-response experiments with pyroxsulam applied after the inhibition of cytochrome P450 monooxygenase by malathion showed an increase in sensitivity in two out of seven highly resistant (HR) rigid ryegrass populations. This demonstrates the presence of non-target-site resistance (NTSR) in some ryegrass populations. Further evidence of NTSR was investigated in dose-response experiments with pyroxsulam, following pretreatment with the glutathione S-transferase (GST) inhibitor 4-chloro-7-nitrobenzoxadiazole (NBD-Cl), which partially reversed resistance in only a few individuals of two L. rigidum populations. Hence, our study confirms the existence of multiple and cross-resistance to ALS- and ACCase-inhibiting herbicides in L. rigidum from Morocco and Tunisia with both TSR and NTSR mechanisms. These results emphasize local resistance management as an important tool to detect and mitigate gene flow from rigid ryegrass populations where resistance has evolved.

Population Genetic Analysis in Persimmons (Diospyros kaki Thunb.) Based on Genome-Wide Single-Nucleotide Polymorphisms

This study investigated the genetic diversity and population structure of a persimmon (Diospyros kaki Thunb., 2n = 6x = 90) collection in South Korea by evaluating 9751 genome-wide single-nucleotide polymorphisms (SNPs) detected using genotyping-by-sequencing in 93 cultivars. The results of neighbor-joining clustering, principal component analysis, and STRUCTURE analysis based on SNPs indicated clear separation between cultivar groups (pollination-constant nonastringent (PCNA, 40 cultivars), pollination-constant astringent (PCA, 19), pollination-variant nonastringent (PVNA, 23), and the pollination-variant astringent type (PVA, 9)) based on the astringency types, while separation between PVA and PVNA-type cultivars was unclear. Population genetic diversity based on SNPs showed that the proportions of polymorphic SNPs within each group ranged from 99.01% (PVNA) to 94.08% (PVA), and the PVNA group exhibited the highest genetic diversity (He = 3.86 and uHe = 0.397). F (fixation index) values were low ranging from -0.024 (PVA) to 0.176 (PCA) with an average of 0.089, indicating a deficiency of heterozygosity. Analysis of molecular variance (AMOVA) and Fst among cultivar groups indicated that variation within individuals was higher than that among the groups. Pairwise Fst values among the groups ranged from 0.01566 (between PVA and PVNA) to 0.09416 (between PCA and PCNA), indicating a low level of cultivar type differentiation. These findings highlight the potential application of biallelic SNPs in population genetics studies of allopolyploids species and provide valuable insights that may have significant implications for breeding and cultivar identification in persimmon.

PURC Provides Improved Sequence Inference for Polyploid Phylogenetics and Other Manifestations of the Multiple-Copy Problem

Inferring the true biological sequences from amplicon mixtures remains a difficult bioinformatics problem. The traditional approach is to cluster sequencing reads by similarity thresholds and treat the consensus sequence of each cluster as an "operational taxonomic unit" (OTU). Recently, this approach has been improved by model-based methods that correct PCR and sequencing errors in order to infer "amplicon sequence variants" (ASVs). To date, ASV approaches have been used primarily in metagenomics, but they are also useful for determining homeologs in polyploid organisms. To facilitate the usage of ASV methods among polyploidy researchers, we incorporated ASV inference alongside OTU clustering in PURC v2.0, a major update to PURC (Pipeline for Untangling Reticulate Complexes). In addition, PURC v2.0 features faster demultiplexing than the original version and has been updated to be compatible with Python 3. In this chapter we present results indicating that using the ASV approach is more likely to infer the correct biological sequences in comparison to the earlier OTU-based PURC and describe how to prepare sequencing data, run PURC v2.0 under several different modes, and interpret the output.

Reticulate Evolution in the Western Mediterranean Mountain Ranges: The Case of the Leucanthemopsis Polyploid Complex

Polyploidization is one of the most common speciation mechanisms in plants. This is particularly relevant in high mountain environments and/or in areas heavily affected by climatic oscillations. Although the role of polyploidy and the temporal and geographical frameworks of polyploidization have been intensively investigated in the alpine regions of the temperate and arctic biomes, fewer studies are available with a specific focus on the Mediterranean region. Leucanthemopsis (Asteraceae) consists of six to ten species with several infraspecific entities, mainly distributed in the western Mediterranean Basin. It is a polyploid complex including montane, subalpine, and strictly alpine lineages, which are locally distributed in different mountain ranges of Western Europe and North Africa. We used a mixed approach including Sanger sequencing and (Roche-454) high throughput sequencing of amplicons to gather information from single-copy nuclear markers and plastid regions. Nuclear regions were carefully tested for recombinants/PCR artifacts and for paralogy. Coalescent-based methods were used to infer the number of polyploidization events and the age of formation of polyploid lineages, and to reconstruct the reticulate evolution of the genus. Whereas the polyploids within the widespread Leucanthemopsis alpina are autopolyploids, the situation is more complex among the taxa endemic to the western Mediterranean. While the hexaploid, L. longipectinata, confined to the northern Moroccan mountain ranges (north-west Africa), is an autopolyploid, the Iberian polyploids are clearly of allopolyploid origins. At least two different polyploidization events gave rise to L. spathulifolia and to all other tetraploid Iberian taxa, respectively. The formation of the Iberian allopolyploids took place in the early Pleistocene and was probably caused by latitudinal and elevational range shifts that brought into contact previously isolated Leucanthemopsis lineages. Our study thus highlights the importance of the Pleistocene climatic oscillations and connected polyploidization events for the high plant diversity in the Mediterranean Basin.

Sequence coverage required for accurate genotyping by sequencing in polyploid species

Polyploidy plays an important role in the evolution of eukaryotes, especially for flowering plants. Many of ecologically or agronomically important plant or crop species are polyploids, including sycamore maple (tetraploid), the world second and third largest food crops wheat (hexaploid) and potato (tetraploid) as well as economically important aquaculture animals such as Atlantic salmon and trout. The next generation sequencing data enables to allocate genotype at a sequence variant site, known as genotyping by sequencing (GBS). GBS has stimulated enormous interests in population based genomics studies in almost all diploid and many polyploid organisms. DNA sequence polymorphisms are codominant and thus fully informative about the underlying genotype at the polymorphic site, making GBS a straightforward task in diploids. However, sequence data may usually be uninformative in polyploid species, making GBS a far more challenging task in polyploids. This paper presents novel and rigorous statistical methods for predicting the number of sequence reads needed to ensure accurate GBS at a polymorphic site bared by the reads in polyploids and shows that a dozen of reads can ensure a probability of 95% to recover all constituent alleles of any tetraploid genotype but several hundreds of reads are needed to accurately uncover the genotype with probability confidence of 90%, subverting the proposition of GBS using low coverage sequence data in the literature. The theoretical prediction was tested by use of RAD-seq data from tetraploid potato cultivars. The paper provides polyploid experimentalists with theoretical guides and methods for designing and conducting their sequence-based studies.

Unravelling hybridization in Phytophthora using phylogenomics and genome size estimation

The genus Phytophthora comprises many economically and ecologically important plant pathogens. Hybrid species have previously been identified in at least six of the 12 phylogenetic clades. These hybrids can potentially infect a wider host range and display enhanced vigour compared to their progenitors. Phytophthora hybrids therefore pose a serious threat to agriculture as well as to natural ecosystems. Early and correct identification of hybrids is therefore essential for adequate plant protection but this is hampered by the limitations of morphological and traditional molecular methods. Identification of hybrids is also important in evolutionary studies as the positioning of hybrids in a phylogenetic tree can lead to suboptimal topologies. To improve the identification of hybrids we have combined genotyping-by-sequencing (GBS) and genome size estimation on a genus-wide collection of 614 Phytophthora isolates. Analyses based on locus- and allele counts and especially on the combination of species-specific loci and genome size estimations allowed us to confirm and characterize 27 previously described hybrid species and discover 16 new hybrid species. Our method was also valuable for species identification at an unprecedented resolution and further allowed correct naming of misidentified isolates. We used both a concatenation- and a coalescent-based phylogenomic method to construct a reliable phylogeny using the GBS data of 140 non-hybrid Phytophthora isolates. Hybrid species were subsequently connected to their progenitors in this phylogenetic tree. In this study we demonstrate the application of two validated techniques (GBS and flow cytometry) for relatively low cost but high resolution identification of hybrids and their phylogenetic relations.

Polyploid phylogenetics

Polyploidy is a dominant feature of extant plant diversity. However, major research questions, including whether polyploidy is important to long-term evolution or is just 'evolutionary noise', remain unresolved due to difficulties associated with the generation and analysis of data from polyploid lineages. Many of these difficulties have been recently overcome, such that it is now often relatively straightforward to infer the full and often reticulate phylogenetic history of groups with recently formed polyploids. This nascent field of 'polyploid phylogenetics' allows researchers to tackle long-standing questions of polyploid macroevolution, supplies the foundation for mechanistic models of ploidy change, and provides the opportunity to include a more complete and representative sample of plant taxa in our analyses in general.

Sampling Variation of RAD-Seq Data from Diploid and Tetraploid Potato (Solanum tuberosum L.)

The new sequencing technology enables identification of genome-wide sequence-based variants at a population level and a competitively low cost. The sequence variant-based molecular markers have motivated enormous interest in population and quantitative genetic analyses. Generation of the sequence data involves a sophisticated experimental process embedded with rich non-biological variation. Statistically, the sequencing process indeed involves sampling DNA fragments from an individual sequence. Adequate knowledge of sampling variation of the sequence data generation is one of the key statistical properties for any downstream analysis of the data and for implementing statistically appropriate methods. This paper reports a thorough investigation on modeling the sampling variation of the sequence data from the optimized RAD-seq (Restriction sit associated DNA sequencing) experiments with two parents and their offspring of diploid and autotetraploid potato (Solanum tuberosum L.). The analysis shows significant dispersion in sampling variation of the sequence data over that expected under multinomial distribution as widely assumed in the literature and provides statistical methods for modeling the variation and calculating the model parameters, which may be easily implemented in real sequence datasets. The optimized design of RAD-seq experiments enabled effective control of presentation of undesirable chloroplast DNA and RNA genes in the sequence data generated.

Optimizing whole-genomic prediction for autotetraploid blueberry breeding

Blueberry (Vaccinium spp.) is an important autopolyploid crop with significant benefits for human health. Apart from its genetic complexity, the feasibility of genomic prediction has been proven for blueberry, enabling a reduction in the breeding cycle time and increasing genetic gain. However, as for other polyploid crops, sequencing costs still hinder the implementation of genome-based breeding methods for blueberry. This motivated us to evaluate the effect of training population sizes and composition, as well as the impact of marker density and sequencing depth on phenotype prediction for the species. For this, data from a large real breeding population of 1804 individuals were used. Genotypic data from 86,930 markers and three traits with different genetic architecture (fruit firmness, fruit weight, and total yield) were evaluated. Herein, we suggested that marker density, sequencing depth, and training population size can be substantially reduced with no significant impact on model accuracy. Our results can help guide decisions toward resource allocation (e.g., genotyping and phenotyping) in order to maximize prediction accuracy. These findings have the potential to allow for a faster and more accurate release of varieties with a substantial reduction of resources for the application of genomic prediction in blueberry. We anticipate that the benefits and pipeline described in our study can be applied to optimize genomic prediction for other diploid and polyploid species.

Multiple auto- and allopolyploidisations marked the Pleistocene history of the widespread Eurasian steppe plant Astragalus onobrychis (Fabaceae)

The Eurasian steppes occupy a significant portion of the worldwide land surface and their biota have been affected by specific past range dynamics driven by ice ages-related climatic fluctuations. The dynamic alterations in conditions during the Pleistocene often triggered reticulate evolution and whole genome duplication events. Employing genomic, genetic and cytogenetic tools as well as morphometry we investigate the intricate evolution of Astragalus onobrychis, a widespread Eurasian steppe plant with diploid, tetraploid and octoploid cytotypes. To analyse the heteroploid RADseq dataset we employ both genotype-based and genotype-free methods that result in highly consistent results, and complement our inference with information from the plastid ycf1 region. We uncover a complex and reticulate evolutionary history, including at least one auto-tetraploidization event and two allo-octoploidization events; one of them involved also genetic contributions from other species, most likely A. goktschaicus. The present genetic structure points to the existence of four main clades within A. onobrychis, which only partly correspond to different ploidies. Time-calibrated diffusion models suggest that diversification within A. onobrychis was associated with ice age-related climatic fluctuations during the last million years. We finally argue for the usefulness of uniparentally inherited plastid markers, even in the genomic era, especially when investigating heteroploid systems.

Recent polyploidization events in three Saccharum founding species

The complexity of polyploid Saccharum genomes hindered progress of genome research and crop improvement in sugarcane. To understand their genome structure, transcriptomes of 59 F₁ individuals derived from S. officinarumLA Purple and S. robustum Molokai 5829 (2n = 80, x = 10 for both) were sequenced, yielding 11 157 and 8998 SNPs and 83 and 105 linkage groups, respectively. Most markers in each linkage group aligned to single sorghum chromosome. However, 71 interchromosomal rearrangements were detected between sorghum and S. officinarum or S. robustum, and 24 (33.8%) of them were shared between S. officinarum and S. robustum, indicating their occurrence before the speciation event that separated these two species. More than 2000 gene pairs from S. spontaneum, S. officinarum and S. robustum were analysed to estimate their divergence time. Saccharum officinarum and S. robustum diverged about 385 thousand years ago, and the whole-genome duplication events occurred after the speciation event because of shared interchromosomal rearrangements. The ancestor of these two species diverged from S. spontaneum about 769 thousand years ago, and the reduction in basic chromosome number from 10 to 8 in S. spontaneum occurred after the speciation event but before the two rounds of whole-genome duplication. Our results proved that S. officinarum is a legitimate species in its own right and not a selection from S. robustum during the domestication process in the past 10 000 years. Our findings rejected a long-standing hypothesis and clarified the timing of speciation and whole-genome duplication events in Saccharum.

Genetic diversity assessed by genotyping by sequencing (GBS) and for phenological traits in blueberry cultivars

Blueberry is a small fruit crop which includes a complex group of different Vaccinium species of various ploidy levels. Commercial blueberries have been grown in Europe most recently, so there is not much information available about their adaptation into new regions. In this work we investigated adaptation to the environmental conditions of northern Spain, in terms of flowering and ripening seasons, of a set of 70 blueberry cultivars including several of the most important cultivated American species (V. corymbosum, V. virgatum, V. macrocarpon and V. uliginosum) in order to identify which types are best-suited in this geographical area of Europe. Most materials showed high chilling requirements for flowering under local conditions, while materials with low-chilling requirements showed problems in the maturation process of the flowers. Most cultivars were early or mid-season while a relative lack of late-season cultivars was observed. GBS was used for the analysis of genetic diversity in this sample of 70 cultivars. A total of 5255 SNP markers were obtained and a cluster analysis revealed three main groups associated with the ploidy level of the species. A Principal Component Analysis revealed a grouping of the V. corymbosum cultivars according to their chilling requirements. A total of 29 SNPs were identified as being highly informative for diversity analysis and potentially useful for cultivar identification and for breeding purposes. The results obtained from this research should contribute to the expansion of this crop, as well as providing data about genetic diversity useful for the preservation of genetic resources or for future breeding programs.

History vs. legend: Retracing invasion and spread of Oxalis pes-caprae L. in Europe and the Mediterranean area

Oxalis pes-caprae L. is a South African geophyte that behaves as an invasive in the eurimediterranean area. According to a long-established hypothesis, O. pes-caprae may have invaded Europe and the Mediterranean area starting from a single plant introduced in the Botanical Garden of Malta at the beginning of the 19^th century. The aim of this work was to test this hypothesis, to track the arrival of O. pes-caprae in different countries of the Euro-Mediterranean area and to understand the pathways of spreading and particularly its starting point(s). Historical data attesting the presence of the plant in the whole Euro-Mediterranean region were collected from different sources: herbarium specimens, Floras and other botanical papers, plant lists of gardens, catalogs of plant nurseries and plant dealers. First records of the plant (both cultivated and wild) for each Territorial Unit (3rd level of NUTS) were selected and used to draw up a diachronic map and an animated graphic. Both documents clearly show that oldest records are scattered throughout the whole area, proving that the plant arrived in Europe and in the Mediterranean region more times independently and that its spreading started in different times from several different centers of invasion. Botanical gardens and other public or private gardens, nurseries and plant dealers, and above all seaside towns and harbors seemingly played a strategic role as a source of either intentional and unintentional introduction or spread. A geographic profiling analysis was performed to analyse the data. We used also techniques (Silhouette, Kmeans and Voronoi tessellation) capable of verifying the presence of more than one independent clusters of data on the basis of their geographical distribution. Microsatellites were employed for a preliminary analysis of genetic variation in the Mediterranean. Even if the sampling was insufficient, particularly among the populations of the original area, our data supported three main groups of populations, one of them corresponding to the central group of populations identified by GP analysis, and the other two corresponding, respectively, to the western and the eastern cluster of data. The most probable areas of origin of the invasion in the three clusters of observations are characterized by the presence of localities where the invasive plant was cultivated, with the exception of the Iberian cluster of observation where the observations in the field predate the data about known cultivation localities. Alternative possible reasons are also suggested, to explain the current prevalence of pentaploid short-styled plants in the Euro-Mediterranean area.

New specific primers for amplification of the Internal Transcribed Spacer region in Clitellata (Annelida)

Nuclear molecular evidence, for example, the rapidly evolving Internal Transcribed Spacer region (ITS), integrated with maternally inherited (mitochondrial) COI barcodes, has provided new insights into the diversity of clitellate annelids. PCR amplification and sequencing of ITS, however, are often hampered by poor specificity of primers used. Therefore, new clitellate‐specific primers for amplifying the whole ITS region (ITS: 29F/1084R) and a part of it (ITS2: 606F/1082R) were developed on the basis of a collection of previously published ITS sequences with flanking rDNA coding regions. The specificity of these and other ITS primers used for clitellates were then tested in silico by evaluating their mismatches with all assembled and annotated sequences (STD, version r127) from EMBL, and the new primers were also tested in vitro for a taxonomically broad sample of clitellate species (71 specimens representing 11 families). The in silico analyses showed that the newly designed primers have a better performance than the universal ones when amplifying clitellate ITS sequences. In vitro PCR and sequencing using the new primers were successful, in particular, for the 606F/1082R pair, which worked well for 65 of the 71 specimens. Thus, using this pair for amplifying the ITS2 will facilitate further molecular systematic investigation of various clitellates. The other pair (29F/1084R), will be a useful complement to existing ITS primers, when amplifying ITS as a whole.

Evaluating allopolyploid origins in strawberries (Fragaria) using haplotypes generated from target capture sequencing

BACKGROUND

Hybridization is observed in many eukaryotic lineages and can lead to the formation of polyploid species. The study of hybridization and polyploidization faces challenges both in data generation and in accounting for population-level phenomena such as coalescence processes in phylogenetic analysis. Genus Fragaria is one example of a set of plant taxa in which a range of ploidy levels is observed across species, but phylogenetic origins are unknown.

RESULTS

Here, using 20 diploid and polyploid Fragaria species, we combine approaches from NGS data analysis and phylogenetics to infer evolutionary origins of polyploid strawberries, taking into account coalescence processes. We generate haplotype sequences for 257 low-copy nuclear markers assembled from Illumina target capture sequence data. We then identify putative hybridization events by analyzing gene tree topologies, and further test predicted hybridizations in a coalescence framework. This approach confirms the allopolyploid ancestry of F. chiloensis and F. virginiana, and provides new allopolyploid ancestry hypotheses for F. iturupensis, F. moschata, and F. orientalis. Evidence of gene flow between diploids F. bucharica and F. vesca is also detected, suggesting that it might be appropriate to consider these groups as conspecifics.

CONCLUSIONS

This study is one of the first in which target capture sequencing followed by computational deconvolution of individual haplotypes is used for tracing origins of polyploid taxa. The study also provides new perspectives on the evolutionary history of Fragaria.

Deep reticulation and incomplete lineage sorting obscure the diploid phylogeny of rain-lilies and allies (Amaryllidaceae tribe Hippeastreae)

Hybridization is a frequent and important force in plant evolution. Next-generation sequencing (NGS) methods offer new possibilities for clade resolution and ambitious sampling of gene genealogies, yet difficulty remains in detecting deep reticulation events using currently available methods. We reconstructed the phylogeny of diploid representatives of Amaryllidaceae tribe Hippeastreae to test the hypothesis of ancient hybridizations preceding the radiation of its major subclade, Hippeastrinae. Through hybrid enrichment of DNA libraries and NGS, we obtained data for 18 nuclear loci through a curated assembly approach and nearly complete plastid genomes for 35 ingroup taxa plus 5 outgroups. Additionally, we obtained alignments for 39 loci through an automated assembly algorithm. These data were analyzed with diverse phylogenetic methods, including concatenation, coalescence-based species tree estimation, Bayesian concordance analysis, and network reconstructions, to provide insights into the evolutionary relationships of Hippeastreae. Causes for gene tree heterogeneity and cytonuclear discordance were examined through a Bayesian posterior predictive approach (JML) and coalescent simulations. Two major clades were found, Hippeastrinae and Traubiinae, as previously reported. Our results suggest the presence of two major nuclear lineages in Hippeastrinae characterized by different chromosome numbers: (1) Tocantinia and Hippeastrum with 2n=22, and (2) Eithea, Habranthus, Rhodophiala, and Zephyranthes mostly with 2n=12, 14, and 18. Strong cytonuclear discordance was confirmed in Hippeastrinae, and a network scenario with at least six hybridization events is proposed to reconcile nuclear and plastid signals, along a backbone that may also have been affected by incomplete lineage sorting at the base of each major subclade.

Next-generation polyploid phylogenetics: rapid resolution of hybrid polyploid complexes using PacBio single-molecule sequencing

Difficulties in generating nuclear data for polyploids have impeded phylogenetic study of these groups. We describe a high-throughput protocol and an associated bioinformatics pipeline (Pipeline for Untangling Reticulate Complexes (Purc)) that is able to generate these data quickly and conveniently, and demonstrate its efficacy on accessions from the fern family Cystopteridaceae. We conclude with a demonstration of the downstream utility of these data by inferring a multi-labeled species tree for a subset of our accessions. We amplified four c. 1-kb-long nuclear loci and sequenced them in a parallel-tagged amplicon sequencing approach using the PacBio platform. Purc infers the final sequences from the raw reads via an iterative approach that corrects PCR and sequencing errors and removes PCR-mediated recombinant sequences (chimeras). We generated data for all gene copies (homeologs, paralogs, and segregating alleles) present in each of three sets of 50 mostly polyploid accessions, for four loci, in three PacBio runs (one run per set). From the raw sequencing reads, Purc was able to accurately infer the underlying sequences. This approach makes it easy and economical to study the phylogenetics of polyploids, and, in conjunction with recent analytical advances, facilitates investigation of broad patterns of polyploid evolution.

Bioinformatics Approach in Plant Genomic Research

The advance in genomics technology leads to the dramatic change in plant biology research. Plant biologists now easily access to enormous genomic data to deeply study plant high-density genetic variation at molecular level. Therefore, fully understanding and well manipulating bioinformatics tools to manage and analyze these data are essential in current plant genome research. Many plant genome databases have been established and continued expanding recently. Meanwhile, analytical methods based on bioinformatics are also well developed in many aspects of plant genomic research including comparative genomic analysis, phylogenomics and evolutionary analysis, and genome-wide association study. However, constantly upgrading in computational infrastructures, such as high capacity data storage and high performing analysis software, is the real challenge for plant genome research. This review paper focuses on challenges and opportunities which knowledge and skills in bioinformatics can bring to plant scientists in present plant genomics era as well as future aspects in critical need for effective tools to facilitate the translation of knowledge from new sequencing data to enhancement of plant productivity.

Genomics-based strategies for the use of natural variation in the improvement of crop metabolism

Next-generation genomics holds great potential in the study of plant phenotypic variation. With several crop reference genomes now available, the affordable costs of de novo genome assembly or target resequencing offer the opportunity to mine the enormous amount of genetic diversity hidden in crop wild relatives. Wide introgressions from these wild ancestors species or land races represent a possible strategy to improve cultivated varieties. In this review, we discuss the mechanisms underlying metabolic diversity within plant species and the possible strategies (and barriers) to introgress novel metabolic traits into cultivated varieties. We show how deep genomic surveys uncover various types of structural variants from extended gene pools of major crops and highlight how this variation may be used for the improvement of crop metabolism.

Diversity and population structure of northern switchgrass as revealed through exome capture sequencing

Panicum virgatum L. (switchgrass) is a polyploid, perennial grass species that is native to North America, and is being developed as a future biofuel feedstock crop. Switchgrass is present primarily in two ecotypes: a northern upland ecotype, composed of tetraploid and octoploid accessions, and a southern lowland ecotype, composed of primarily tetraploid accessions. We employed high-coverage exome capture sequencing (~2.4 Tb) to genotype 537 individuals from 45 upland and 21 lowland populations. From these data, we identified ~27 million single-nucleotide polymorphisms (SNPs), of which 1 590 653 high-confidence SNPs were used in downstream analyses of diversity within and between the populations. From the 66 populations, we identified five primary population groups within the upland and lowland ecotypes, a result that was further supported through genetic distance analysis. We identified conserved, ecotype-restricted, non-synonymous SNPs that are predicted to affect the protein function of CONSTANS (CO) and EARLY HEADING DATE 1 (EHD1), key genes involved in flowering, which may contribute to the phenotypic differences between the two ecotypes. We also identified, relative to the near-reference Kanlow population, 17 228 genes present in more copies than in the reference genome (up-CNVs), 112 630 genes present in fewer copies than in the reference genome (down-CNVs) and 14 430 presence/absence variants (PAVs), affecting a total of 9979 genes, including two upland-specific CNV clusters. In total, 45 719 genes were affected by an SNP, CNV, or PAV across the panel, providing a firm foundation to identify functional variation associated with phenotypic traits of interest for biofuel feedstock production.

Species-Level Phylogeny and Polyploid Relationships in Hordeum (Poaceae) Inferred by Next-Generation Sequencing and In Silico Cloning of Multiple Nuclear Loci

Polyploidization is an important speciation mechanism in the barley genus Hordeum. To analyze evolutionary changes after allopolyploidization, knowledge of parental relationships is essential. One chloroplast and 12 nuclear single-copy loci were amplified by polymerase chain reaction (PCR) in all Hordeum plus six out-group species. Amplicons from each of 96 individuals were pooled, sheared, labeled with individual-specific barcodes and sequenced in a single run on a 454 platform. Reference sequences were obtained by cloning and Sanger sequencing of all loci for nine supplementary individuals. The 454 reads were assembled into contigs representing the 13 loci and, for polyploids, also homoeologues. Phylogenetic analyses were conducted for all loci separately and for a concatenated data matrix of all loci. For diploid taxa, a Bayesian concordance analysis and a coalescent-based dated species tree was inferred from all gene trees. Chloroplast matK was used to determine the maternal parent in allopolyploid taxa. The relative performance of different multilocus analyses in the presence of incomplete lineage sorting and hybridization was also assessed. The resulting multilocus phylogeny reveals for the first time species phylogeny and progenitor-derivative relationships of all di- and polyploid Hordeum taxa within a single analysis. Our study proves that it is possible to obtain a multilocus species-level phylogeny for di- and polyploid taxa by combining PCR with next-generation sequencing, without cloning and without creating a heavy load of sequence data.

Molecular phylogeny of anoplocephalid tapeworms (Cestoda: Anoplocephalidae) infecting humans and non-human primates

Anoplocephalid tapeworms of the genus Bertiella Stiles and Hassall, 1902 and Anoplocephala Blanchard, 1848, found in the Asian, African and American non-human primates are presumed to sporadic ape-to-man transmissions. Variable nuclear (5.8S-ITS2; 28S rRNA) and mitochondrial genes (cox1; nad1) of isolates of anoplocephalids originating from different primates (Callicebus oenanthe, Gorilla beringei, Gorilla gorilla, Pan troglodytes and Pongo abelii) and humans from various regions (South America, Africa, South-East Asia) were sequenced. In most analyses, Bertiella formed a monophyletic group within the subfamily Anoplocephalinae, however, the 28S rRNA sequence-based analysis indicated paraphyletic relationship between Bertiella from primates and Australian marsupials and rodents, which should thus be regarded as different taxa. Moreover, isolate determined as Anoplocephala cf. gorillae from mountain gorilla clustered within the Bertiella clade from primates. This either indicates that A. gorillae deserves to be included into the genus Bertiella, or, that an unknown Bertiella species infects also mountain gorillas. The analyses allowed the genetic differentiation of the isolates, albeit with no obvious geographical or host-related patterns. The unexpected genetic diversity of the isolates studied suggests the existence of several Bertiella species in primates and human and calls for revision of the whole group, based both on molecular and morphological data.

Single Nucleotide Polymorphism Identification in Polyploids: A Review, Example, and Recommendations

Understanding the relationship between genotype and phenotype is a major biological question and being able to predict phenotypes based on molecular genotypes is integral to molecular breeding. Whole-genome duplications have shaped the history of all flowering plants and present challenges to elucidating the relationship between genotype and phenotype, especially in neopolyploid species. Although single nucleotide polymorphisms (SNPs) have become popular tools for genetic mapping, discovery and application of SNPs in polyploids has been difficult. Here, we summarize common experimental approaches to SNP calling, highlighting recent polyploid successes. To examine the impact of software choice on these analyses, we called SNPs among five peanut genotypes using different alignment programs (BWA-mem and Bowtie 2) and variant callers (SAMtools, GATK, and Freebayes). Alignments produced by Bowtie 2 and BWA-mem and analyzed in SAMtools shared 24.5% concordant SNPs, and SAMtools, GATK, and Freebayes shared 1.4% concordant SNPs. A subsequent analysis of simulated Brassica napus chromosome 1A and 1C genotypes demonstrated that, of the three software programs, SAMtools performed with the highest sensitivity and specificity on Bowtie 2 alignments. These results, however, are likely to vary among species, and we therefore propose a series of best practices for SNP calling in polyploids.

Parallel tagged amplicon sequencing of relatively long PCR products using the Illumina HiSeq platform and transcriptome assembly

In phylogenetics and population genetics, a large number of loci are often needed to accurately resolve species relationships. Normally, loci are enriched by PCR and sequenced by Sanger sequencing, which is expensive when the number of amplicons is large. Next-generation sequencing (NGS) techniques are increasingly used for parallel amplicon sequencing, which reduces sequencing costs tremendously, but has not reduced preparation costs very much. Moreover, for most current NGS methods, amplicons need to be purified and quantified before sequencing and their lengths are also restricted (normally <700 bp). Here, we describe an approach to sequence pooled amplicons of any length using the Illumina platform. Using this method, amplicons are pooled at equal volume rather than at equal concentration, thus eliminating the laborious purification and quantification steps. We then shear the pooled amplicons, repair the ends, add sample identifying linkers and pool multiple samples prior to Illumina library preparation. Data are then assembled using the transcriptome assembly program trinity, which is optimized to deal with templates of highly varying quantities. We demonstrated the utility of our approach by recovering 93.5% of the target amplicons (size up to 1650 bp) in full length for a 16 taxa × 101 loci project, using ~2.0 GB of Illumina HiSeq paired-end 90-bp data. Overall, we validate a rapid, cost-effective and scalable approach to sequence a large number of targeted loci from a large number of samples that is particularly suitable for both phylogenetics and population genetics studies that require a modest scale of data.

Detection of induced mutations in CaFAD2 genes by next-generation sequencing leading to the production of improved oil composition in Crambe abyssinica

Crambe abyssinica is a hexaploid oil crop for industrial applications. An increase of erucic acid (C22:1) and reduction of polyunsaturated fatty acid (PUFA) contents in crambe oil is a valuable improvement. An increase in oleic acid (C18:1), a reduction in PUFA and possibly an increase in C22:1 can be obtained by down-regulating the expression of fatty acid desaturase2 genes (CaFAD2), which code for the enzyme that converts C18:1 into C18:2. We conducted EMS-mutagenesis in crambe, followed by Illumina sequencing, to screen mutations in three expressed CaFAD2 genes. Two novel analysis strategies were used to detect mutation sites. In the first strategy, mutation detection targeted specific sequence motifs. In the second strategy, every nucleotide position in a CaFAD2 fragment was tested for the presence of mutations. Seventeen novel mutations were detected in 1100 one-dimensional pools (11 000 individuals) in three expressed CaFAD2 genes, including non-sense mutations and mis-sense mutations in CaFAD2-C1, -C2 and -C3. The homozygous non-sense mutants for CaFAD2-C3 resulted in a 25% higher content of C18:1 and 25% lower content of PUFA compared to the wild type. The mis-sense mutations only led to small changes in oil composition. Concluding, targeted mutation detection using NGS in a polyploid was successfully applied and it was found that a non-sense mutation in even a single CaFAD2 gene can lead to changes in crambe oil composition. Stacking the mutations in different CaFAD2 may gain additional changes in C18:1 and PUFA contents.

Biological invasions, climate change and genomics

The rate of biological invasions is expected to increase as the effects of climate change on biological communities become widespread. Climate change enhances habitat disturbance which facilitates the establishment of invasive species, which in turn provides opportunities for hybridization and introgression. These effects influence local biodiversity that can be tracked through genetic and genomic approaches. Metabarcoding and metagenomic approaches provide a way of monitoring some types of communities under climate change for the appearance of invasives. Introgression and hybridization can be followed by the analysis of entire genomes so that rapidly changing areas of the genome are identified and instances of genetic pollution monitored. Genomic markers enable accurate tracking of invasive species' geographic origin well beyond what was previously possible. New genomic tools are promoting fresh insights into classic questions about invading organisms under climate change, such as the role of genetic variation, local adaptation and climate pre-adaptation in successful invasions. These tools are providing managers with often more effective means to identify potential threats, improve surveillance and assess impacts on communities. We provide a framework for the application of genomic techniques within a management context and also indicate some important limitations in what can be achieved.

Genomic profile of the plants with pharmaceutical value

There is an ample genetic diversity of plants with medicinal importance around the globe and this pool of genetic variation serves as the base for selection as well as for plant improvement. Thus, identification, characterization and documentation of the gene pool of medicinal plants are essential for this purpose. Genomic information of many a medicinal plant species has increased rapidly since the past decade and genetic resources available for domestication and improvement programs include genome sequencing, expressed sequence tags sequencing, transcript profiling, gene transmit, molecular markers in favor of mapping and breeding. In recent years, multiple endeavors have been undertaken for genomic characterization of medicinal plant species with the aid of molecular markers for sustainable utilization of gene pool, its conservation and future studies. Recent advancement in genomics is so fast that only some researches have been published till date and to a large extent documentation is restricted to electronic resources. Whole genome profiling of the identified medicinal plant species, carried out by several researchers, based on the DNA fingerprinting, is well documented in the present review. This review will facilitate preparing a database of the widely used, economically important medicinal plant species, based on their genomic organization.

Limited genetic divergence among Australian alpine Poa tussock grasses coupled with regional structuring points to ongoing gene flow and taxonomic challenges

BACKGROUND AND AIMS

While molecular approaches can often accurately reconstruct species relationships, taxa that are incompletely differentiated pose a challenge even with extensive data. Such taxa are functionally differentiated, but may be genetically differentiated only at small and/or patchy regions of the genome. This issue is considered here in Poa tussock grass species that dominate grassland and herbfields in the Australian alpine zone.

METHODS

Previously reported tetraploidy was confirmed in all species by sequencing seven nuclear regions and five microsatellite markers. A Bayesian approach was used to co-estimate nuclear and chloroplast gene trees with an overall dated species tree. The resulting species tree was used to examine species structure and recent hybridization, and intertaxon fertility was tested by experimental crosses.

KEY RESULTS

Species tree estimation revealed Poa gunnii, a Tasmanian endemic species, as sister to the rest of the Australian alpine Poa. The taxa have radiated in the last 0·5-1·2 million years and the non-gunnii taxa are not supported as genetically distinct. Recent hybridization following past species divergence was also not supported. Ongoing gene flow is suggested, with some broad-scale geographic structure within the group.

CONCLUSIONS

The Australian alpine Poa species are not genetically distinct despite being distinguishable phenotypically, suggesting recent adaptive divergence with ongoing intertaxon gene flow. This highlights challenges in using conventional molecular taxonomy to infer species relationships in recent, rapid radiations.

Identification of putative orthologous genes for the phylogenetic reconstruction of temperate woody bamboos (Poaceae: Bambusoideae)

The temperate woody bamboos (Arundinarieae) are highly diverse in morphology but lack a substantial amount of genetic variation. The taxonomy of this lineage is intractable, and the relationships within the tribe have not been well resolved. Recent studies indicated that this tribe could have a complex evolutionary history. Although phylogenetic studies of the tribe have been carried out, most of these phylogenetic reconstructions were based on plastid data, which provide lower phylogenetic resolution compared with nuclear data. In this study, we intended to identify a set of desirable nuclear genes for resolving the phylogeny of the temperate woody bamboos. Using two different methodologies, we identified 209 and 916 genes, respectively, as putative single copy orthologous genes. A total of 112 genes was successfully amplified and sequenced by next-generation sequencing technologies in five species sampled from the tribe. As most of the genes exhibited intra-individual allele heterozygotes, we investigated phylogenetic utility by reconstructing the phylogeny based on individual genes. Discordance among gene trees was observed and, to resolve the conflict, we performed a range of analyses using BUCKy and HybTree. While caution should be taken when inferring a phylogeny from multiple conflicting genes, our analysis indicated that 74 of the 112 investigated genes are potential markers for resolving the phylogeny of the temperate woody bamboos.

DNA fingerprinting in botany: past, present, future

Almost three decades ago Alec Jeffreys published his seminal Nature papers on the use of minisatellite probes for DNA fingerprinting of humans (Jeffreys and colleagues Nature 1985, 314:67-73 and Nature 1985, 316:76-79). The new technology was soon adopted for many other organisms including plants, and when Hilde Nybom, Kurt Weising and Alec Jeffreys first met at the very First International Conference on DNA Fingerprinting in Berne, Switzerland, in 1990, everybody was enthusiastic about the novel method that allowed us for the first time to discriminate between humans, animals, plants and fungi on the individual level using DNA markers. A newsletter coined "Fingerprint News" was launched, T-shirts were sold, and the proceedings of the Berne conference filled a first book on "DNA fingerprinting: approaches and applications". Four more conferences were about to follow, one on each continent, and Alec Jeffreys of course was invited to all of them. Since these early days, methodologies have undergone a rapid evolution and diversification. A multitude of techniques have been developed, optimized, and eventually abandoned when novel and more efficient and/or more reliable methods appeared. Despite some overlap between the lifetimes of the different technologies, three phases can be defined that coincide with major technological advances. Whereas the first phase of DNA fingerprinting ("the past") was dominated by restriction fragment analysis in conjunction with Southern blot hybridization, the advent of the PCR in the late 1980s gave way to the development of PCR-based single- or multi-locus profiling techniques in the second phase. Given that many routine applications of plant DNA fingerprinting still rely on PCR-based markers, we here refer to these methods as "DNA fingerprinting in the present", and include numerous examples in the present review. The beginning of the third phase actually dates back to 2005, when several novel, highly parallel DNA sequencing strategies were developed that increased the throughput over current Sanger sequencing technology 1000-fold and more. High-speed DNA sequencing was soon also exploited for DNA fingerprinting in plants, either in terms of facilitated marker development, or directly in the sense of "genotyping-by-sequencing". Whereas these novel approaches are applied at an ever increasing rate also in non-model species, they are still far from routine, and we therefore treat them here as "DNA fingerprinting in the future".

Using multi-locus allelic sequence data to estimate genetic divergence among four Lilium (Liliaceae) cultivars

Next Generation Sequencing (NGS) may enable estimating relationships among genotypes using allelic variation of multiple nuclear genes simultaneously. We explored the potential and caveats of this strategy in four genetically distant Lilium cultivars to estimate their genetic divergence from transcriptome sequences using three approaches: POFAD (Phylogeny of Organisms from Allelic Data, uses allelic information of sequence data), RAxML (Randomized Accelerated Maximum Likelihood, tree building based on concatenated consensus sequences) and Consensus Network (constructing a network summarizing among gene tree conflicts). Twenty six gene contigs were chosen based on the presence of orthologous sequences in all cultivars, seven of which also had an orthologous sequence in Tulipa, used as out-group. The three approaches generated the same topology. Although the resolution offered by these approaches is high, in this case there was no extra benefit in using allelic information. We conclude that these 26 genes can be widely applied to construct a species tree for the genus Lilium.

Genetic and genomic interactions of animals with different ploidy levels

Polyploid animals have independently evolved from diploids in diverse taxa across the tree of life. We review a few polyploid animal species or biotypes where recently developed molecular and cytogenetic methods have significantly improved our understanding of their genetics, reproduction and evolution. Mitochondrial sequences that target the maternal ancestor of a polyploid show that polyploids may have single (e.g. unisexual salamanders in the genus Ambystoma) or multiple (e.g. parthenogenetic polyploid lizards in the genus Aspidoscelis) origins. Microsatellites are nuclear markers that can be used to analyze genetic recombinations, reproductive modes (e.g. Ambystoma) and recombination events (e.g. polyploid frogs such as Pelophylax esculentus). Hom(e)ologous chromosomes and rare intergenomic exchanges in allopolyploids have been distinguished by applying genome-specific fluorescent probes to chromosome spreads. Polyploids arise, and are maintained, through perturbations of the 'normal' meiotic program that would include pre-meiotic chromosome replication and genomic integrity of homologs. When possible, asexual, unisexual and bisexual polyploid species or biotypes interact with diploid relatives, and genes are passed from diploid to polyploid gene pools, which increase genetic diversity and ultimately evolutionary flexibility in the polyploid. When diploid relatives do not exist, polyploids can interact with another polyploid (e.g. species of African Clawed Frogs in the genus Xenopus). Some polyploid fish (e.g. salmonids) and frogs (Xenopus) represent independent lineages whose ancestors experienced whole genome duplication events. Some tetraploid frogs (P. esculentus) and fish (Squaliusalburnoides) may be in the process of becoming independent species, but diploid and triploid forms of these 'species' continue to genetically interact with the comparatively few tetraploid populations. Genetic and genomic interaction between polyploids and diploids is a complex and dynamic process that likely plays a crucial role for the evolution and persistence of polyploid animals. See also other articles in this themed issue.

A next-generation sequencing method for genotyping-by-sequencing of highly heterozygous autotetraploid potato

Assessment of genomic DNA sequence variation and genotype calling in autotetraploids implies the ability to distinguish among five possible alternative allele copy number states. This study demonstrates the accuracy of genotyping-by-sequencing (GBS) of a large collection of autotetraploid potato cultivars using next-generation sequencing. It is still costly to reach sufficient read depths on a genome wide scale, across the cultivated gene pool. Therefore, we enriched cultivar-specific DNA sequencing libraries using an in-solution hybridisation method (SureSelect). This complexity reduction allowed to confine our study to 807 target genes distributed across the genomes of 83 tetraploid cultivars and one reference (DM 1–3 511). Indexed sequencing libraries were paired-end sequenced in 7 pools of 12 samples using Illumina HiSeq2000. After filtering and processing the raw sequence data, 12.4 Gigabases of high-quality sequence data was obtained, which mapped to 2.1 Mb of the potato reference genome, with a median average read depth of 63× per cultivar. We detected 129,156 sequence variants and genotyped the allele copy number of each variant for every cultivar. In this cultivar panel a variant density of 1 SNP/24 bp in exons and 1 SNP/15 bp in introns was obtained. The average minor allele frequency (MAF) of a variant was 0.14. Potato germplasm displayed a large number of relatively rare variants and/or haplotypes, with 61% of the variants having a MAF below 0.05. A very high average nucleotide diversity (π = 0.0107) was observed. Nucleotide diversity varied among potato chromosomes. Several genes under selection were identified. Genotyping-by-sequencing results, with allele copy number estimates, were validated with a KASP genotyping assay. This validation showed that read depths of ∼60–80× can be used as a lower boundary for reliable assessment of allele copy number of sequence variants in autotetraploids. Genotypic data were associated with traits, and alleles strongly influencing maturity and flesh colour were identified.

Reprint of: using nuclear gene data for plant phylogenetics: progress and prospects

The paper reviews the current state of low and single copy nuclear markers that have been applied successfully in plant phylogenetics to date, and discusses case studies highlighting the potential of massively parallel high throughput or next-generation sequencing (NGS) approaches for molecular phylogenetic and evolutionary investigations. The current state, prospects and challenges of specific single- or low-copy plant nuclear markers as well as phylogenomic case studies are presented and evaluated.

Spatiotemporal reconstruction of the Aquilegia rapid radiation through next-generation sequencing of rapidly evolving cpDNA regions

Aquilegia is a well-known model system in the field of evolutionary biology, but obtaining a resolved and well-supported phylogenetic reconstruction for the genus has been hindered by its recent and rapid diversification. Here, we applied 454 next-generation sequencing to PCR amplicons of 21 of the most rapidly evolving regions of the plastome to generate c. 24 kb of sequences from each of 84 individuals from throughout the genus. The resulting phylogeny has well-supported resolution of the main lineages of the genus, although recent diversification such as in the European taxa remains unresolved. By producing a chronogram of the whole Ranunculaceae family based on published data, we inferred calibration points for dating the Aquilegia radiation. The genus originated in the upper Miocene c. 6.9 million yr ago (Ma) in Eastern Asia, and diversification occurred c. 4.8 Ma with the split of two main clades, one colonizing North America, and the other Western Eurasia through the mountains of Central Asia. This was followed by a back-to-Asia migration, originating from the European stock using a North Asian route. These results provide the first backbone phylogeny and spatiotemporal reconstruction of the Aquilegia radiation, and constitute a robust framework to address the adaptative nature of speciation within the group.

A DNA barcoding method to discriminate between the model plant Brachypodium distachyon and its close relatives B. stacei and B. hybridum (Poaceae)

BACKGROUND

Brachypodium distachyon s. l. has been widely investigated across the world as a model plant for temperate cereals and biofuel grasses. However, this annual plant shows three cytotypes that have been recently recognized as three independent species, the diploids B. distachyon (2n = 10) and B. stacei (2n = 20) and their derived allotetraploid B. hybridum (2n = 30).

METHODOLOGY/PRINCIPAL FINDINGS

We propose a DNA barcoding approach that consists of a rapid, accurate and automatable species identification method using the standard DNA sequences of complementary plastid (trnLF) and nuclear (ITS, GI) loci. The highly homogenous but largely divergent B. distachyon and B. stacei diploids could be easily distinguished (100% identification success) using direct trnLF (2.4%), ITS (5.5%) or GI (3.8%) sequence divergence. By contrast, B. hybridum could only be unambiguously identified through the use of combined trnLF+ITS sequences (90% of identification success) or by cloned GI sequences (96.7%) that showed 5.4% (ITS) and 4% (GI) rate divergence between the two parental sequences found in the allopolyploid.

CONCLUSION/SIGNIFICANCE

Our data provide an unbiased and effective barcode to differentiate these three closely-related species from one another. This procedure overcomes the taxonomic uncertainty generated from methods based on morphology or flow cytometry identifications that have resulted in some misclassifications of the model plant and its allies. Our study also demonstrates that the allotetraploid B. hybridum has resulted from bi-directional crosses of B. distachyon and B. stacei plants acting either as maternal or paternal parents.

Parallel tagged amplicon sequencing reveals major lineages and phylogenetic structure in the North American tiger salamander (Ambystoma tigrinum) species complex

Modern analytical methods for population genetics and phylogenetics are expected to provide more accurate results when data from multiple genome-wide loci are analysed. We present the results of an initial application of parallel tagged sequencing (PTS) on a next-generation platform to sequence thousands of barcoded PCR amplicons generated from 95 nuclear loci and 93 individuals sampled across the range of the tiger salamander (Ambystoma tigrinum) species complex. To manage the bioinformatic processing of this large data set (344 330 reads), we developed a pipeline that sorts PTS data by barcode and locus, identifies high-quality variable nucleotides and yields phased haplotype sequences for each individual at each locus. Our sequencing and bioinformatic strategy resulted in a genome-wide data set with relatively low levels of missing data and a wide range of nucleotide variation. structure analyses of these data in a genotypic format resulted in strongly supported assignments for the majority of individuals into nine geographically defined genetic clusters. Species tree analyses of the most variable loci using a multi-species coalescent model resulted in strong support for most branches in the species tree; however, analyses including more than 50 loci produced parameter sampling trends that indicated a lack of convergence on the posterior distribution. Overall, these results demonstrate the potential for amplicon-based PTS to rapidly generate large-scale data for population genetic and phylogenetic-based research.

Patterns of linkage disequilibrium and association mapping in diploid alfalfa (M. sativa L.)

Association mapping enables the detection of marker-trait associations in unstructured populations by taking advantage of historical linkage disequilibrium (LD) that exists between a marker and the true causative polymorphism of the trait phenotype. Our first objective was to understand the pattern of LD decay in the diploid alfalfa genome. We used 89 highly polymorphic SSR loci in 374 unimproved diploid alfalfa (Medicago sativa L.) genotypes from 120 accessions to infer chromosome-wide patterns of LD. We also sequenced four lignin biosynthesis candidate genes (caffeoyl-CoA 3-O-methyltransferase (CCoAoMT), ferulate-5-hydroxylase (F5H), caffeic acid-O-methyltransferase (COMT), and phenylalanine amonialyase (PAL 1)) to identify single nucleotide polymorphisms (SNPs) and infer within gene estimates of LD. As the second objective of this study, we conducted association mapping for cell wall components and agronomic traits using the SSR markers and SNPs from the four candidate genes. We found very little LD among SSR markers implying limited value for genomewide association studies. In contrast, within gene LD decayed within 300 bp below an r (2) of 0.2 in three of four candidate genes. We identified one SSR and two highly significant SNPs associated with biomass yield. Based on our results, focusing association mapping on candidate gene sequences will be necessary until a dense set of genome-wide markers is available for alfalfa.

A tailed PCR procedure for cost-effective, two-order multiplex sequencing of candidate genes in polyploid plants

Complex polyploid crop genomes can be recalcitrant towards conventional DNA sequencing approaches for allele mining in candidate genes for valuable traits. In the past, this has greatly complicated the transfer of knowledge on promising candidate genes from model plants to even closely related polyploid crops. Next-generation sequencing offers diverse solutions to overcome such difficulties. Here, we present a method for multiplexed 454 sequencing in gene-specific PCR amplicons that can simultaneously address multiple homologues of given target genes. We devised a simple two-step PCR procedure employing a set of barcoded M13/T7 universal fusion primers that enable a cost-effective and efficient amplification of large numbers of target gene amplicons. Sequencing-ready amplicons are generated that can be simultaneously sequenced in pools comprising multiple amplicons from multiple genotypes. High-depth sequencing allows resolution of the resulting sequence reads into contigs representing multiple homologous loci, with only insignificant off-target capture of paralogues or PCR artefacts. In a case study, the procedure was tested in the complex polyploid genome of Brassica napus for a set of nine genes identified in Arabidopsis as candidates for regulation of seed development and oil content. Up to six copies of these genes were expected in B. napus. SNP discovery was performed by pooled multiplex sequencing of 30 amplicons in 20 diverse B. napus accessions with interesting trait variation for oil content, providing a basis for comparative mapping to relevant quantitative trait loci and for subsequent marker-assisted breeding.

Targeted capture of homoeologous coding and noncoding sequence in polyploid cotton

Targeted sequence capture is a promising technology in many areas in biology. These methods enable efficient and relatively inexpensive sequencing of hundreds to thousands of genes or genomic regions from many more individuals than is practical using whole-genome sequencing approaches. Here, we demonstrate the feasibility of target enrichment using sequence capture in polyploid cotton. To capture and sequence both members of each gene pair (homeologs) of wild and domesticated Gossypium hirsutum, we created custom hybridization probes to target 1000 genes (500 pairs of homeologs) using information from the cotton transcriptome. Two widely divergent samples of G. hirsutum were hybridized to four custom NimbleGen capture arrays containing probes for targeted genes. We show that the two coresident homeologs in the allopolyploid nucleus were efficiently captured with high coverage. The capture efficiency was similar between the two accessions and independent of whether the samples were multiplexed. A significant amount of flanking, nontargeted sequence (untranslated regions and introns) was also captured and sequenced along with the targeted exons. Intraindividual heterozygosity is low in both wild and cultivated Upland cotton, as expected from the high level of inbreeding in natural G. hirsutum and bottlenecks accompanying domestication. In addition, levels of heterozygosity appeared asymmetrical with respect to genome (A_T or D_T) in cultivated cotton. The approach used here is general, scalable, and may be adapted for many different research inquiries involving polyploid plant genomes.

Next-generation hybridization and introgression

Hybridization has a major role in evolution-from the introgression of important phenotypic traits between species, to the creation of new species through hybrid speciation. Molecular studies of hybridization aim to understand the class of hybrids and the frequency of introgression, detect the signature of ancient hybridization, and understand the behaviour of introgressed loci in their new genomic background. This often involves a large investment in the design and application of molecular markers, leading to a compromise between the depth and breadth of genomic data. New techniques designed to assay a large sub-section of the genome, in association with next-generation sequencing (NGS) technologies, will allow genome-wide hybridization and introgression studies in organisms with no prior sequence data. These detailed genotypic data will unite the breadth of sampling of loci characteristic of population genetics with the depth of sequence information associated with molecular phylogenetics. In this review, we assess the theoretical and methodological constraints that limit our understanding of natural hybridization, and promote the use of NGS for detecting hybridization and introgression between non-model organisms. We also make recommendations for the ways in which emerging techniques, such as pooled barcoded amplicon sequencing and restriction site-associated DNA tags, should be used to overcome current limitations, and enhance our understanding of this evolutionary significant process.