SNiPloid: A Utility to Exploit High-Throughput SNP Data Derived from RNA-Seq in Allopolyploid Species

Allopolyploid subgenome identification and implications for evolutionary analysis

Whole-genome duplications (WGDs) are widespread genomic events in eukaryotes that are hypothesized to contribute to the evolutionary success of many lineages, including flowering plants, Saccharomyces yeast, and vertebrates. WGDs generally can be classified into autopolyploids (ploidy increase descended from one species) or allopolyploids (ploidy increase descended from multiple species). Assignment of allopolyploid progenitor species (called subgenomes in the polyploid) is important to understanding the biology and evolution of polyploids, including the asymmetric subgenome evolution following hybridization (biased fractionation). Here, I review the different methodologies used to identify the ancestors of allopolyploid subgenomes, discuss the advantages and disadvantages of these methods, and outline the implications of how these methods affect the subsequent evolutionary analysis of these genomes.

Species delimitation 4.0: integrative taxonomy meets artificial intelligence

Although species are central units for biological research, recent findings in genomics are raising awareness that what we call species can be ill-founded entities due to solely morphology-based, regional species descriptions. This particularly applies to groups characterized by intricate evolutionary processes such as hybridization, polyploidy, or asexuality. Here, challenges of current integrative taxonomy (genetics/genomics + morphology + ecology, etc.) become apparent: different favored species concepts, lack of universal characters/markers, missing appropriate analytical tools for intricate evolutionary processes, and highly subjective ranking and fusion of datasets. Now, integrative taxonomy combined with artificial intelligence under a unified species concept can enable automated feature learning and data integration, and thus reduce subjectivity in species delimitation. This approach will likely accelerate revising and unraveling eukaryotic biodiversity.

Meeting in Liguria: hybridisation between Apennine endemic Euphorbia barrelieri and western Mediterranean E. nicaeensis led to the allopolyploid origin of E. ligustica

The Mediterranean Basin is renowned for its extremely rich biota and is considered as one of the 25 Global Biodiversity Hotspots, but its diversity is not homogeneously distributed. Outstanding in the number of (endemic) species are the Ligurian Alps (Italy). At the foot of the Ligurian Alps, little above the Mediterranean Sea, a disjunct occurrence of Italian endemic Euphorbia barrelieri was reported. Using an array of integrative methods ranging from cytogenetic (chromosome number and relative genome size estimation), over phylogenetic approaches (plastid, ITS and RAD sequencing) to multivariate morphometrics we disentangled the origin of these populations that were shown to be tetraploid. We performed phylogenetic analyses of the nuclear ITS and plastid regions of a broad taxonomic sampling of Euphorbia sect. Pithyusa to identify possible species involved in the origin of the tetraploid populations and then applied various analyses of RADseq data to identify the putative parental species. Our results have shown that the Ligurian populations of E. barrelieri are of allotetraploid origin that involved E. barrelieri and western Mediterranean E. nicaeensis as parental species. We thus describe a new species, E. ligustica, and hypothesise that its adaptation to similar environments in which E. barrelieri occurs, triggered development of similar morphology, whereas its genetic composition appears to be closer to that of E. nicaeensis. Our study emphasises the importance of polyploidisation for plant diversification, highlights the value of the Ligurian Alps as a hotspot of biodiversity and endemism and underlines the importance of integrative taxonomic approaches in uncovering cryptic diversity.

Untying Gordian knots: unraveling reticulate polyploid plant evolution by genomic data using the large Ranunculus auricomus species complex

Speciation via hybridization and polyploidization is a major evolutionary force in plant evolution but is still poorly understood for neopolyploid groups. Challenges are attributed to high heterozygosity, low genetic divergence, and missing information on progenitors, ploidy, and reproduction. We study the large Eurasian Ranunculus auricomus species complex and use a comprehensive workflow integrating reduced-representation sequencing (RRS) genomic data to unravel reticulate evolution, genome diversity and composition of polyploids. We rely on 97 312 restriction site-associated DNA sequencing (RAD-Seq) loci, 576 targeted nuclear genes (48 phased), and 71 plastid regions derived from 78 polyploid apomictic taxa and four diploid and one tetraploid putative sexual progenitor species. We applied (phylo)genomic structure, network, and single nucleotide polymorphism (SNP)-origin analyses. Results consistently showed only 3-5 supported and geographically structured polyploid genetic groups, each containing extant sexual and one unknown progenitor species. Combined analyses demonstrated predominantly allopolyploid origins, each involving 2-3 different diploid sexual progenitor species. Young allotetraploids were characterized by subgenome dominance and nonhybrid SNPs, suggesting substantial post-origin but little lineage-specific evolution. The biodiversity of neopolyploid complexes can result from multiple hybrid origins involving different progenitors and substantial post-origin evolution (e.g. homoeologous exchanges, hybrid segregation, gene flow). Reduced-representation sequencing genomic data including multi-approach information is efficient to delimit shallow reticulate relationships.

Novel Approaches for Species Concepts and Delimitation in Polyploids and Hybrids

Hybridization and polyploidization are important processes for plant evolution. However, classification of hybrid or polyploid species has been notoriously difficult because of the complexity of processes and different evolutionary scenarios that do not fit with classical species concepts. Polyploid complexes are formed via combinations of allopolyploidy, autopolyploidy and homoploid hybridization with persisting sexual reproduction, resulting in many discrete lineages that have been classified as species. Polyploid complexes with facultative apomixis result in complicated net-work like clusters, or rarely in agamospecies. Various case studies illustrate the problems that apply to traditional species concepts to hybrids and polyploids. Conceptual progress can be made if lineage formation is accepted as an inevitable consequence of meiotic sex, which is established already in the first eukaryotes as a DNA restoration tool. The turnaround of the viewpoint that sex forms species as lineages helps to overcome traditional thinking of species as "units". Lineage formation and self-sustainability is the prerequisite for speciation and can also be applied to hybrids and polyploids. Species delimitation is aided by the improved recognition of lineages via various novel -omics methods, by understanding meiosis functions, and by recognizing functional phenotypes by considering morphological-physiological-ecological adaptations.

Subgenome evolution in allotetraploid plants

Polyploidization is a well-known speciation and adaptation mechanism. Traces of former polyploidization events were discovered within many genomes, and especially in plants. Allopolyploidization by interspecific hybridization between two species is common. Among hybrid plants, many are domesticated species of agricultural interest and many of their genomes and of their presumptive parents have been sequenced. Hybrid genomes remain challenging to analyse because of the presence of multiple subgenomes. The genomes of hybrids often undergo rearrangement and degradation over time. Based on 10 hybrid plant genomes from six different genera, with hybridization dating from 10,000 to 5 million years ago, we assessed subgenome degradation, subgenomic intermixing and biased subgenome fractionation. The restructuring of hybrid genomes does not proceed proportionally with the age of the hybrid. The oldest hybrids in our data set display completely different fates: whereas the subgenomes of the tobacco plant Nicotiana benthamiana are in an advanced stage of degradation, the subgenomes of quinoa (Chenopodium quinoa) are exceptionally well conserved by structure and sequence. We observed statistically significant biased subgenome fractionation in seven out of 10 hybrids, which had different ages and subgenomic intermixing levels. Hence, we conclude that no correlation exists between biased fractionation and subgenome intermixing. Lastly, domestication may encourage or hinder subgenome intermixing, depending on the evolutionary context. In summary, comparative analysis of hybrid genomes and their presumptive parents allowed us to determine commonalities and differences between their evolutionary fates. In order to facilitate the future analysis of further hybrid genomes, we automated the analysis steps within manticore, which is publicly available at https://github.com/MatteoSchiavinato/manticore.git.

Phylogenomic Relationships and Evolution of Polyploid Salix Species Revealed by RAD Sequencing Data

Polyploidy is common in the genus Salix. However, little is known about the origin, parentage and genomic composition of polyploid species because of a lack of suitable molecular markers and analysis tools. We established a phylogenomic framework including species of all described sections of Eurasian shrub willows. We analyzed the genomic composition of seven polyploid willow species in comparison to putative diploid parental species to draw conclusions on their origin and the effects of backcrossing and post-origin evolution. We applied recently developed programs like SNAPP, HyDe, and SNiPloid to establish a bioinformatic pipeline for unravelling the complexity of polyploid genomes. RAD sequencing revealed 23,393 loci and 320,010 high quality SNPs for the analysis of relationships of 35 species of Eurasian shrub willows (Salix subg. Chamaetia/Vetrix). Polyploid willow species appear to be predominantly of allopolyploid origin. More ancient allopolyploidization events were observed for two hexaploid and one octoploid species, while our data suggested a more recent allopolyploid origin for the included tetraploids and identified putative parental taxa. SNiPloid analyses disentangled the different genomic signatures resulting from hybrid origin, backcrossing, and secondary post-origin evolution in the polyploid species. Our RAD sequencing data demonstrate that willow genomes are shaped by ancient and recent reticulate evolution, polyploidization, and post-origin divergence of species.

Homoeologous gene expression and co-expression network analyses and evolutionary inference in allopolyploids

Polyploidy is a widespread phenomenon throughout eukaryotes. Due to the coexistence of duplicated genomes, polyploids offer unique challenges for estimating gene expression levels, which is essential for understanding the massive and various forms of transcriptomic responses accompanying polyploidy. Although previous studies have explored the bioinformatics of polyploid transcriptomic profiling, the causes and consequences of inaccurate quantification of transcripts from duplicated gene copies have not been addressed. Using transcriptomic data from the cotton genus (Gossypium) as an example, we present an analytical workflow to evaluate a variety of bioinformatic method choices at different stages of RNA-seq analysis, from homoeolog expression quantification to downstream analysis used to infer key phenomena of polyploid expression evolution. In general, EAGLE-RC and GSNAP-PolyCat outperform other quantification pipelines tested, and their derived expression dataset best represents the expected homoeolog expression and co-expression divergence. The performance of co-expression network analysis was less affected by homoeolog quantification than by network construction methods, where weighted networks outperformed binary networks. By examining the extent and consequences of homoeolog read ambiguity, we illuminate the potential artifacts that may affect our understanding of duplicate gene expression, including an overestimation of homoeolog co-regulation and the incorrect inference of subgenome asymmetry in network topology. Taken together, our work points to a set of reasonable practices that we hope are broadly applicable to the evolutionary exploration of polyploids.

Merging genotypes: graft union formation and scion-rootstock interactions

Grafting has been utilised for at least the past 7000 years. Historically, grafting has been developed by growers without particular interest beyond the agronomical and ornamental effects, and thus knowledge about grafting has remained largely empirical. Much of the commercial production of fruit, and increasingly vegetables, relies upon grafting with rootstocks to provide resistance to soil-borne pathogens and abiotic stresses as well as to influence scion growth and performance. Although there is considerable agronomic knowledge about the use and selection of rootstocks for many species, we know little of the molecular mechanisms underlying rootstock adaptation to different soil environments and rootstock-conferred modifications of scion phenotypes. Furthermore, the processes involved in the formation of the graft union and graft compatibility are poorly understood despite over a hundred years of scientific study. In this paper, we provide an overview of what is known about grafting and the mechanisms underlying rootstock-scion interactions. We highlight recent studies that have advanced our understanding of graft union formation and outline subjects that require further development.

Reference Transcriptomes and Detection of Duplicated Copies in Hexaploid and Allododecaploid Spartina Species (Poaceae)

In this study, we report the assembly and annotation of five reference transcriptomes for the European hexaploid Spartina species (S. maritima, S. alterniflora and their homoploid hybrids S. x townsendii and S. x neyrautii) and the allododecaploid invasive species S. anglica These transcriptomes were constructed from various leaf and root cDNA libraries that were sequenced using both Roche-454 and Illumina technologies. Considering the high ploidy levels of the Spartina genomes under study, and considering the absence of diploid reference genome and the need of an appropriate analytical strategy, we developed generic bioinformatics tools to (1) detect different haplotypes of each gene within each species and (2) assign a parental origin to haplotypes detected in the hexaploid hybrids and the neo-allopolyploid. The approach described here allows the detection of putative homeologs from sets of short reads. Synonymous substitution rate (KS) comparisons between haplotypes from the hexaploid species revealed the presence of one KS peak (likely resulting from the tetraploid duplication event). The procedure developed in this study can be applied for future differential gene expression or genomics experiments to study the fate of duplicated genes in the invasive allododecaploid S. anglica.

De-novo transcriptome assembly for gene identification, analysis, annotation, and molecular marker discovery in Onobrychis viciifolia

BACKGROUND

Sainfoin (Onobrychis viciifolia) is a highly nutritious tannin-containing forage legume. In the diet of ruminants sainfoin can have anti-parasitic effects and reduce methane emissions under in vitro conditions. Many of these benefits have been attributed to condensed tannins or proanthocyanidins in sainfoin. A combination of increased use of industrially produced nitrogen fertilizer, issues with establishment and productivity in the first year and more reliable alternatives, such as red clover led to a decline in the use of sainfoin since the middle of the last century. In recent years there has been a resurgence of interest in sainfoin due to its potential beneficial nutraceutical and environmental attributes. However, genomic resources are scarce, thus hampering progress in genetic analysis and improvement. To address this we have used next generation RNA sequencing technology to obtain the first transcriptome of sainfoin. We used the library to identify gene-based simple sequence repeats (SSRs) and potential single nucleotide polymorphisms (SNPs).

RESULTS

One genotype from each of five sainfoin accessions was sequenced. Paired-end (PE) sequences were generated from cDNA libraries of RNA extracted from 7 day old seedlings. A combined assembly of 92,772 transcripts was produced de novo using the Trinity programme. About 18,000 transcripts were annotated with at least one GO (gene ontology) term. A total of 63 transcripts were annotated as involved in the tannin biosynthesis pathway. We identified 3786 potential SSRs. SNPs were identified by mapping the reads of the individual assemblies against the combined assembly. After stringent filtering a total of 77,000 putative SNPs were identified. A phylogenetic analysis of single copy number genes showed that sainfoin was most closely related to red clover and Medicago truncatula, while Lotus japonicus, bean and soybean are more distant relatives.

CONCLUSIONS

This work describes the first transcriptome assembly in sainfoin. The 92 K transcripts provide a rich source of SNP and SSR polymorphisms for future use in genetic studies of this crop. Annotation of genes involved in the condensed tannin biosynthesis pathway has provided the basis for further studies of the genetic control of this important trait in sainfoin.

Gene coexpression network analysis of oil biosynthesis in an interspecific backcross of oil palm

Global demand for vegetable oils is increasing at a dramatic rate, while our understanding of the regulation of oil biosynthesis in plants remains limited. To gain insights into the mechanisms that govern oil synthesis and fatty acid (FA) composition in the oil palm fruit, we used a multilevel approach combining gene coexpression analysis, quantification of allele-specific expression and joint multivariate analysis of transcriptomic and lipid data, in an interspecific backcross population between the African oil palm, Elaeis guineensis, and the American oil palm, Elaeis oleifera, which display contrasting oil contents and FA compositions. The gene coexpression network produced revealed tight transcriptional coordination of fatty acid synthesis (FAS) in the plastid with sugar sensing, plastidial glycolysis, transient starch storage and carbon recapture pathways. It also revealed a concerted regulation, along with FAS, of both the transfer of nascent FA to the endoplasmic reticulum, where triacylglycerol assembly occurs, and of the production of glycerol-3-phosphate, which provides the backbone of triacylglycerols. Plastid biogenesis and auxin transport were the two other biological processes most tightly connected to FAS in the network. In addition to WRINKLED1, a transcription factor (TF) known to activate FAS genes, two novel TFs, termed NF-YB-1 and ZFP-1, were found at the core of the FAS module. The saturated FA content of palm oil appeared to vary above all in relation to the level of transcripts of the gene coding for β-ketoacyl-acyl carrier protein synthase II. Our findings should facilitate the development of breeding and engineering strategies in this and other oil crops.

Unidirectional diploid-tetraploid introgression among British birch trees with shifting ranges shown by restriction site-associated markers

Hybridization may lead to introgression of genes among species. Introgression may be bidirectional or unidirectional, depending on factors such as the demography of the hybridizing species, or the nature of reproductive barriers between them. Previous microsatellite studies suggested bidirectional introgression between diploid Betula nana (dwarf birch) and tetraploid B. pubescens (downy birch) and also between B. pubescens and diploid B. pendula (silver birch) in Britain. Here, we analyse introgression among these species using 51 237 variants in restriction site‐associated (RAD) markers in 194 individuals, called with allele dosages in the tetraploids. In contrast to the microsatellite study, we found unidirectional introgression into B. pubescens from both of the diploid species. This pattern fits better with the expected nature of the reproductive barrier between diploids and tetraploids. As in the microsatellite study, introgression into B. pubescens showed clear clines with increasing introgression from B. nana in the north and from B. pendula in the south. Unlike B. pendula alleles, introgression of B. nana alleles was found far from the current area of sympatry or allopatry between B. nana and B. pubescens. This pattern fits a shifting zone of hybridization due to Holocene reduction in the range of B. nana and expansion in the range of B. pubescens.

Haplotype Detection from Next-Generation Sequencing in High-Ploidy-Level Species: 45S rDNA Gene Copies in the Hexaploid Spartina maritima

Gene and whole-genome duplications are widespread in plant nuclear genomes, resulting in sequence heterogeneity. Identification of duplicated genes may be particularly challenging in highly redundant genomes, especially when there are no diploid parents as a reference. Here, we developed a pipeline to detect the different copies in the ribosomal RNA gene family in the hexaploid grass Spartina maritima from next-generation sequencing (Roche-454) reads. The heterogeneity of the different domains of the highly repeated 45S unit was explored by identifying single nucleotide polymorphisms (SNPs) and assembling reads based on shared polymorphisms. SNPs were validated using comparisons with Illumina sequence data sets and by cloning and Sanger (re)sequencing. Using this approach, 29 validated polymorphisms and 11 validated haplotypes were reported (out of 34 and 20, respectively, that were initially predicted by our program). The rDNA domains of S. maritima have similar lengths as those found in other Poaceae, apart from the 5′-ETS, which is approximately two-times longer in S. maritima. Sequence homogeneity was encountered in coding regions and both internal transcribed spacers (ITS), whereas high intragenomic variability was detected in the intergenic spacer (IGS) and the external transcribed spacer (ETS). Molecular cytogenetic analysis by fluorescent in situ hybridization (FISH) revealed the presence of one pair of 45S rDNA signals on the chromosomes of S. maritima instead of three expected pairs for a hexaploid genome, indicating loss of duplicated homeologous loci through the diploidization process. The procedure developed here may be used at any ploidy level and using different sequencing technologies.

Evolutionary dynamics of retrotransposons following autopolyploidy in the Buckler Mustard species complex

Long terminal repeat retrotransposons (LTR-RTs) represent a major fraction of plant genomes, but processes leading to transposition bursts remain elusive. Polyploidy expectedly leads to LTR-RT proliferation, as the merging of divergent diploids provokes a genome shock activating LTR-RTs and/or genetic redundancy supports the accumulation of active LTR-RTs through relaxation of selective constraints. Available evidence supports interspecific hybridization as the main trigger of genome dynamics, but few studies have addressed the consequences of intraspecific polyploidy (i.e. autopolyploidy), where the genome shock is expectedly minimized. The dynamics of LTR-RTs was thus here evaluated through low coverage 454 sequencing of three closely related diploid progenitors and three independent autotetraploids from the young Biscutella laevigata species complex. Genomes from this early diverging Brassicaceae lineage presented a minimum of 40% repeats and a large diversity of transposable elements. Differential abundances and patterns of sequence divergence among genomes for 37 LTR-RT families revealed contrasted dynamics during species diversification. Quiescent LTR-RT families with limited genetic variation among genomes were distinguished from active families (37.8%) having proliferated in specific taxa. Specific families proliferated in autopolyploids only, but most transpositionally active families in polyploids were also differentiated among diploids. Low expression levels of transpositionally active LTR-RT families in autopolyploids further supported that genome shock and redundancy are non-mutually exclusive triggers of LTR-RT proliferation. Although reputed stable, autopolyploid genomes show LTR-RT fractions presenting analogies with polyploids between widely divergent genomes.

Methods for mapping and categorization of DNA sequence reads from allopolyploid organisms

Genome read categorization determines the genome of origin for sequence reads from an allopolyploid organism. Different techniques have been used to perform read categorization, mostly based on homoeo-SNPs identified between extant diploid relatives of allopolyploids. We present a novel technique for read categorization implemented by the software PolyDog. We demonstrate its accuracy and improved categorization relative to other methods. We discuss the situations in which one method or another might be most appropriate.

Genome rearrangements derived from homoeologous recombination following allopolyploidy speciation in coffee

Allopolyploidization is widespread and has played a major role in flowering plant diversification. Genomic changes are common consequences of allopolyploidization, but their mechanisms of occurrence and dynamics over time are still poorly understood. Coffea arabica, a recently formed allotetraploid, was chosen as a model to investigate genetic changes in allopolyploid using an approach that exploits next-generation sequencing technologies. Genes affected by putative homoeolog loss were inferred by comparing the numbers of single-nucleotide polymorphisms detected using RNA-seq in individual accessions of C. arabica, and between accessions of its two diploid progenitor species for common sequence positions. Their physical locations were investigated and clusters of genes exhibiting homoeolog loss were identified. To validate these results, genome sequencing data were generated from one accession of C. arabica and further analyzed. Genomic rearrangements involving homoeologous exchanges appear to occur in C. arabica and to be a major source of genetic diversity. At least 5% of the C. arabica genes were inferred to have undergone homoeolog loss. The detection of a large number of homoeologous exchange events (HEEs) shared by all accessions of C. arabica strongly reinforces the assumption of a single allopolyploidization event. Furthermore, HEEs were specific to one or a few accessions, suggesting that HEE accumulates gradually. Our results provide evidence for the important role of HEE in allopolyploid genome evolution.