Responses of Transposable Elements to Polyploidy

Studying Whole-Genome Duplication Using Experimental Evolution of Chlamydomonas

In this chapter, we present the use of Chlamydomonas reinhardtii in experiments designed to study the evolutionary impacts of whole genome duplication. We shortly introduce the algal species and depict why it is an excellent model for experimental evolution. Subsequently, we discuss the most relevant steps and methods in the design of a ploidy-related Chlamydomonas experiment. These steps include strain selection, ploidy determination, different methods of making diplo- and polyploid Chlamydomonas cells, replication, culturing conditions, preservation, and the ways to quantify phenotypic and genotypic change.

TEMM: A Curated Data Resource for Transposon Element-Based Molecular Markers in Plants

Transposon elements (TEs) are mobile genetic elements that can insert themselves into new locations and modify the plant genome. In recent years, they have been used as molecular markers in plant breeding programs. TE-based molecular markers (TE-markers) are divided into two categories depending on the transcription mechanism of the TEs. The first category is retrotransposon-based molecular markers, which include RBIP, IRAP, REMAP, and iPBS. The second group is DNA-based-TE-markers, which include MITE, TE-junction, and CACTA TE-markers. These markers are a good tool for studying genetic diversity and can provide information on plants' phylogenetic and evolutionary history. They can help improve breeding programs to increase agronomic traits and develop new varieties. Overall, TE-markers play an important role in plant genetics and plant breeding and contribute to a better understanding of plant biology. Here, we present TEMM, a curated data resource for TE-markers in plants. Relevant research articles were screened to collect primer sequences and related information. Only articles containing primer sequences are added to the present data resource. TEMM contains 784 primers with their associated PCR reaction programs and their applications in various crops. These include 203 IPBS, 191 RBIP, 140 IRAP, 78 TE-junction, 76 IRAPS, 47 RBIP-IRAP, 16 IRAP-REMAP, 12 REMAP, 12 REMA-IRAP, 6 REMA, and 3 ISBP primers. The data resource is freely available at https://bioinformatics.um6p.ma/TEMM .

Power and Weakness of Repetition - Evaluating the Phylogenetic Signal From Repeatomes in the Family Rosaceae With Two Case Studies From Genera Prone to Polyploidy and Hybridization (Rosa and Fragaria)

Plant genomes consist, to a considerable extent, of non-coding repetitive DNA. Several studies showed that phylogenetic signals can be extracted from such repeatome data by using among-species dissimilarities from the RepeatExplorer2 pipeline as distance measures. Here, we advanced this approach by adjusting the read input for comparative clustering indirectly proportional to genome size and by summarizing all clusters into a main distance matrix subjected to Neighbor Joining algorithms and Principal Coordinate Analyses. Thus, our multivariate statistical method works as a "repeatomic fingerprint," and we proved its power and limitations by exemplarily applying it to the family Rosaceae at intrafamilial and, in the genera Fragaria and Rosa, at the intrageneric level. Since both taxa are prone to hybridization events, we wanted to show whether repeatome data are suitable to unravel the origin of natural and synthetic hybrids. In addition, we compared the results based on complete repeatomes with those from ribosomal DNA clusters only, because they represent one of the most widely used barcoding markers. Our results demonstrated that repeatome data contained a clear phylogenetic signal supporting the current subfamilial classification within Rosaceae. Accordingly, the well-accepted major evolutionary lineages within Fragaria were distinguished, and hybrids showed intermediate positions between parental species in data sets retrieved from both complete repeatomes and rDNA clusters. Within the taxonomically more complicated and particularly frequently hybridizing genus Rosa, we detected rather weak phylogenetic signals but surprisingly found a geographic pattern at a population scale. In sum, our method revealed promising results at larger taxonomic scales as well as within taxa with manageable levels of reticulation, but success remained rather taxon specific. Since repeatomes can be technically easy and comparably inexpensively retrieved even from samples of rather poor DNA quality, our phylogenomic method serves as a valuable alternative when high-quality genomes are unavailable, for example, in the case of old museum specimens.

Ancient gene duplications, rather than polyploidization, facilitate diversification of petal pigmentation patterns in Clarkia gracilis (Onagraceae)

It has been suggested that gene duplication and polyploidization create opportunities for the evolution of novel characters. However, the connections between the effects of polyploidization and morphological novelties have rarely been examined. In this study, we investigated whether petal pigmentation patterning in an allotetraploid Clarkia gracilis has evolved as a result of polyploidization. C. gracilis is thought to be derived through a recent polyploidization event with two diploid species, C. amoena huntiana and an extinct species that is closely related to C. lassenensis. We reconstructed phylogenetic relationships of the R2R3-MYBs (the regulators of petal pigmentation) from two subspecies of C. gracilis and the two purported progenitors, C. a. huntiana and C. lassenensis. The gene tree reveals that these R2R3-MYB genes have arisen through duplications that occurred before the divergence of the two progenitor species, i.e., before polyploidization. After polyploidization and subsequent gene loss, only one of the two orthologous copies inherited from the progenitors was retained in the polyploid, turning it to diploid inheritance. We examined evolutionary changes in these R2R3-MYBs and in their expression, which reveals that the changes affecting patterning (including expression domain contraction, loss-of-function mutation, cis-regulatory mutation) occurred after polyploidization within the C. gracilis lineages. Our results thus suggest that polyploidization itself is not necessary in producing novel petal color patterns. By contrast, duplications of R2R3-MYB genes in the common ancestor of the two progenitors have apparently facilitated diversification of petal pigmentation patterns.

Characterization and Dynamics of Repeatomes in Closely Related Species of Hieracium (Asteraceae) and Their Synthetic and Apomictic Hybrids

The repetitive content of the plant genome (repeatome) often represents its largest fraction and is frequently correlated with its size. Transposable elements (TEs), the main component of the repeatome, are an important driver in the genome diversification due to their fast-evolving nature. Hybridization and polyploidization events are hypothesized to induce massive bursts of TEs resulting, among other effects, in an increase of copy number and genome size. Little is known about the repeatome dynamics following hybridization and polyploidization in plants that reproduce by apomixis (asexual reproduction via seeds). To address this, we analyzed the repeatomes of two diploid parental species, Hieracium intybaceum and H. prenanthoides (sexual), their diploid F1 synthetic and their natural triploid hybrids (H. pallidiflorum and H. picroides, apomictic). Using low-coverage next-generation sequencing (NGS) and a graph-based clustering approach, we detected high overall similarity across all major repeatome categories between the parental species, despite their large phylogenetic distance. Medium and highly abundant repetitive elements comprise ∼70% of Hieracium genomes; most prevalent were Ty3/Gypsy chromovirus Tekay and Ty1/Copia Maximus-SIRE elements. No TE bursts were detected, neither in synthetic nor in natural hybrids, as TE abundance generally followed theoretical expectations based on parental genome dosage. Slight over- and under-representation of TE cluster abundances reflected individual differences in genome size. However, in comparative analyses, apomicts displayed an overabundance of pararetrovirus clusters not observed in synthetic hybrids. Substantial deviations were detected in rDNAs and satellite repeats, but these patterns were sample specific. rDNA and satellite repeats (three of them were newly developed as cytogenetic markers) were localized on chromosomes by fluorescence in situ hybridization (FISH). In a few cases, low-abundant repeats (5S rDNA and certain satellites) showed some discrepancy between NGS data and FISH results, which is due partly to the bias of low-coverage sequencing and partly to low amounts of the satellite repeats or their sequence divergence. Overall, satellite DNA (including rDNA) was markedly affected by hybridization, but independent of the ploidy or reproductive mode of the progeny, whereas bursts of TEs did not play an important role in the evolutionary history of Hieracium.

Copy Number Variation of Transposable Elements in Thinopyrum intermedium and Its Diploid Relative Species

Diploid and polyploid wild species of Triticeae have complex relationships, and the understanding of their evolution and speciation could help to increase the usability of them in wheat breeding as a source of genetic diversity. The diploid species Pseudoroegneria spicata (St), Thinopyrum bessarabicum (Jb), Dasypyrum villosum (V) derived from a hypothetical common ancestor are considered to be possible subgenome donors in hexaploid species Th. intermedium (JrJvsSt, where indices r, v, and s stand for the partial relation to the genomes of Secale, Dasypyrum, and Pseudoroegneria, respectively). We quantified 10 families of transposable elements (TEs) in P. spicata, Th. bessarabicum, D. villosum (per one genome), and Th. intermedium (per one average subgenome) using the quantitative real time PCR assay and compared their abundance within the studied genomes as well as between them. Sabrina was the most abundant among all studied elements in P. spicata, D. villosum, and Th. intermedium, and among Ty3/Gypsy elements in all studied species. Among Ty1/Copia elements, Angela-A and WIS-A showed the highest and close abundance with the exception of D. villosum, and comprised the majority of all studied elements in Th. bessarabicum. Sabrina, BAGY2, and Angela-A showed similar abundance among diploids and in Th. intermedium hexaploid; Latidu and Barbara demonstrated sharp differences between diploid genomes. The relationships between genomes of Triticeae species based on the studied TE abundance and the role of TEs in speciation and polyploidization in the light of the current phylogenetic models is discussed.

Relaxed purifying selection in autopolyploids drives transposable element over-accumulation which provides variants for local adaptation

Polyploidization is frequently associated with increased transposable element (TE) content. However, what drives TE dynamics following whole genome duplication (WGD) and the evolutionary implications remain unclear. Here, we leverage whole-genome resequencing data available for ~300 individuals of Arabidopsis arenosa, a well characterized natural diploid-autotetraploid plant species, to address these questions. Based on 43,176 TE insertions we detect in these genomes, we demonstrate that relaxed purifying selection rather than transposition bursts is the main driver of TE over-accumulation after WGD. Furthermore, the increased pool of TE insertions in tetraploids is especially enriched within or near environmentally responsive genes. Notably, we show that the major flowering-time repressor gene FLC is disrupted by a TE insertion specifically in the rapid-cycling tetraploid lineage that colonized mainland railways. Together, our findings indicate that tetrasomy leads to an enhanced accumulation of genic TE insertions, some of which likely contribute to local adaptation.

Why transposable elements (TEs) accumulate in polyploids and the evolutionary implications remain unclear. Here, the authors show that following whole genome duplication, relaxed purifying selection is the main driver of TE over-accumulation, which provides variants for rapid local adaptation.

Comparative effect of allopolyploidy on transposable element composition and gene expression between Gossypium hirsutum and its two diploid progenitors

An allopolyploidization event formed allotetraploid Gossypium species from an A-genome diploid species and a D-genome diploid species. To explore the responses of transposable elements (TEs) to allopolyploidy, we assembled parallel TE datasets from G. hirsutum, G. arboreum and G. raimondii and analyzed the TE types and the effects of TEs on orthologous gene expression in the three Gossypium genomes. Gypsy was the most abundant TE type and most TEs were located ∼500 bp from genes in all three genomes. In G. hirsutum, 35.6% of genes harbored TE insertions, whereas insertions were more frequent in G. arboreum and G. raimondii. G. hirsutum had the highest proportion of uniquely matching 24-nt small interfering RNAs (siRNAs) that targeted TEs. TEs, particularly those targeted by 24-nt siRNAs, were associated with reduced gene expression, but the effect of TEs on orthologous gene expression varied substantially among species. Orthologous gene expression levels in G. hirsutum were intermediate between those of G. arboreum and G. raimondii, which did not experience TE expansion or reduction resulting from allopolyploidization. This study underscores the diversity of TEs co-opted by host genes and provides insights into the roles of TEs in regulating gene expression in Gossypium.

Dating the Species Network: Allopolyploidy and Repetitive DNA Evolution in American Daisies (Melampodium sect. Melampodium, Asteraceae)

Allopolyploidy has played an important role in the evolution of the flowering plants. Genome mergers are often accompanied by significant and rapid alterations of genome size and structure via chromosomal rearrangements and altered dynamics of tandem and dispersed repetitive DNA families. Recent developments in sequencing technologies and bioinformatic methods allow for a comprehensive investigation of the repetitive component of plant genomes. Interpretation of evolutionary dynamics following allopolyploidization requires both the knowledge of parentage and the age of origin of an allopolyploid. Whereas parentage is typically inferred from cytogenetic and phylogenetic data, age inference is hampered by the reticulate nature of the phylogenetic relationships. Treating subgenomes of allopolyploids as if they belonged to different species (i.e., no recombination among subgenomes) and applying cross-bracing (i.e., putting a constraint on the age difference of nodes pertaining to the same event), we can infer the age of allopolyploids within the framework of the multispecies coalescent within BEAST2. Together with a comprehensive characterization of the repetitive DNA fraction using the RepeatExplorer pipeline, we apply the dating approach in a group of closely related allopolyploids and their progenitor species in the plant genus Melampodium (Asteraceae). We dated the origin of both the allotetraploid, Melampodium strigosum, and its two allohexaploid derivatives, Melampodium pringlei and Melampodium sericeum, which share both parentage and the direction of the cross, to the Pleistocene ($<$1.4 Ma). Thus, Pleistocene climatic fluctuations may have triggered formation of allopolyploids possibly in short intervals, contributing to difficulties in inferring the precise temporal order of allopolyploid species divergence of M. sericeum and M. pringlei. The relatively recent origin of the allopolyploids likely played a role in the near-absence of major changes in the repetitive fraction of the polyploids' genomes. The repetitive elements most affected by the postpolyploidization changes represented retrotransposons of the Ty1-copia lineage Maximus and, to a lesser extent, also Athila elements of Ty3-gypsy family.

Patterns of Evolutionary Trajectories and Domestication History within the Genus Hordeum Assessed by REMAP Markers

The patterns of genetic diversity related to the taxonomy and domestication history of 85 accessions representing the main four species of the genus Hordeum were examined by retrotransposon-microsatellite amplified polymorphism (REMAP) markers based on the retrotransposon BARE-1. A substantial level of genetic polymorphisms at among- and within-species level was observed showing that this retrotransposon family and its adjacent genomic regions has been a target for genome dynamics during the evolution and domestication of barley. The obtained data are consistent with the current taxonomic status within the genus Hordeum. Similar level of genetic diversity was observed between the wild and the domesticated barley accessions suggesting that transposable elements` activity and accumulation may counteract the decrease of genome-wide diversity following domestication. In addition, eco-geographical sub-genome pools of the cultivated barley were identified in support to the theory of multiple origins of domestication within the genus Hordeum. We also provide conclusions about the relationship between accessions of different species and the putative routes of barley domestication. In conclusion, the retrotransposon BARE-1 stands as a reliable and perspective DNA marker for the assessment of the phylogenetic and domestication history in the genus Hordeum and other crop species.

Heritable gene expression differences between apomictic clone members in Taraxacum officinale: Insights into early stages of evolutionary divergence in asexual plants

BACKGROUND

Asexual reproduction has the potential to enhance deleterious mutation accumulation and to constrain adaptive evolution. One source of mutations that can be especially relevant in recent asexuals is activity of transposable elements (TEs), which may have experienced selection for high transposition rates in sexual ancestor populations. Predictions of genomic divergence under asexual reproduction therefore likely include a large contribution of transposable elements but limited adaptive divergence. For plants empirical insight into genome divergence under asexual reproduction remains limited. Here, we characterize expression divergence between clone members of a single apomictic lineage of the common dandelion (Taraxacum officinale) to contribute to our knowledge of genome evolution under asexuality.

RESULTS

Using RNA-Seq, we show that about one third of heritable divergence within the apomictic lineage is driven by TEs and TE-related gene activity. In addition, we identify non-random transcriptional differences in pathways related to acyl-lipid and abscisic acid metabolisms which might reflect functional divergence within the apomictic lineage. We analyze SNPs in the transcriptome to assess genetic divergence between the apomictic clone members and reveal that heritable expression differences between the accessions are not explained simply by genome-wide genetic divergence.

CONCLUSION

The present study depicts a first effort towards a more complete understanding of apomictic plant genome evolution. We identify abundant TE activity and ecologically relevant functional genes and pathways affecting heritable within-lineage expression divergence. These findings offer valuable resources for future work looking at epigenetic silencing and Cis-regulation of gene expression with particular emphasis on the effects of TE activity on asexual species' genome.

Genome reorganization in F1 hybrids uncovers the role of retrotransposons in reproductive isolation

Interspecific hybridization leads to new interactions among divergent genomes, revealing the nature of genetic incompatibilities having accumulated during and after the origin of species. Conflicts associated with misregulation of transposable elements (TEs) in hybrids expectedly result in their activation and genome-wide changes that may be key to species boundaries. Repetitive genomes of wild wheats have diverged under differential dynamics of specific long terminal repeat retrotransposons (LTR-RTs), offering unparalleled opportunities to address the underpinnings of plant genome reorganization by selfish sequences. Using reciprocal F1 hybrids between three Aegilops species, restructuring and epigenetic repatterning was assessed at random and LTR-RT sequences with amplified fragment length polymorphism and sequence-specific amplified polymorphisms as well as their methylation-sensitive counterparts, respectively. Asymmetrical reorganization of LTR-RT families predicted to cause conflicting interactions matched differential survival of F1 hybrids. Consistent with the genome shock model, increasing divergence of merged LTR-RTs yielded higher levels of changes in corresponding genome fractions and lead to repeated reorganization of LTR-RT sequences in F1 hybrids. Such non-random reorganization of hybrid genomes is coherent with the necessary repression of incompatible TE loci in support of hybrid viability and indicates that TE-driven genomic conflicts may represent an overlooked factor supporting reproductive isolation.

The striking and unexpected cytogenetic diversity of genus Tanacetum L. (Asteraceae): a cytometric and fluorescent in situ hybridisation study of Iranian taxa

BACKGROUND

Although karyologically well studied, the genus Tanacetum (Asteraceae) is poorly known from the perspective of molecular cytogenetics. The prevalence of polyploidy, including odd ploidy warranted an extensive cytogenetic study. We studied several species native to Iran, one of the most important centres of diversity of the genus. We aimed to characterise Tanacetum genomes through fluorochrome banding, fluorescent in situ hybridisation (FISH) of rRNA genes and the assessment of genome size by flow cytometry. We appraise the effect of polyploidy and evaluate the existence of intraspecific variation based on the number and distribution of GC-rich bands and rDNA loci. Finally, we infer ancestral genome size and other cytogenetic traits considering phylogenetic relationships within the genus.

RESULTS

We report first genome size (2C) estimates ranging from 3.84 to 24.87 pg representing about 11 % of those recognised for the genus. We found striking cytogenetic diversity both in the number of GC-rich bands and rDNA loci. There is variation even at the population level and some species have undergone massive heterochromatic or rDNA amplification. Certain morphometric data, such as pollen size or inflorescence architecture, bear some relationship with genome size. Reconstruction of ancestral genome size, number of CMA+ bands and number of rDNA loci show that ups and downs have occurred during the evolution of these traits, although genome size has mostly increased and the number of CMA+ bands and rDNA loci have decreased in present-day taxa compared with ancestral values.

CONCLUSIONS

Tanacetum genomes are highly unstable in the number of GC-rich bands and rDNA loci, although some patterns can be established at the diploid and tetraploid levels. In particular, aneuploid taxa and some odd ploidy species show greater cytogenetic instability than the rest of the genus. We have also confirmed a linked rDNA arrangement for all the studied Tanacetum species. The labile scenario found in Tanacetum proves that some cytogenetic features previously regarded as relatively constant, or even diagnostic, can display high variability, which is better interpreted within a phylogenetic context.

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.

Polyploidy and novelty: Gottlieb's legacy

Nearly four decades ago, Roose & Gottlieb (Roose & Gottlieb 1976 Evolution 30, 818-830. (doi:10.2307/2407821)) showed that the recently derived allotetraploids Tragopogon mirus and T. miscellus combined the allozyme profiles of their diploid parents (T. dubius and T. porrifolius, and T. dubius and T. pratensis, respectively). This classic paper addressed the link between genotype and biochemical phenotype and documented enzyme additivity in allopolyploids. Perhaps more important than their model of additivity, however, was their demonstration of novelty at the biochemical level. Enzyme multiplicity-the production of novel enzyme forms in the allopolyploids-can provide an extensive array of polymorphism for a polyploid individual and may explain, for example, the expanded ranges of polyploids relative to their diploid progenitors. In this paper, we extend the concept of evolutionary novelty in allopolyploids to a range of genetic and ecological features. We observe that the dynamic nature of polyploid genomes-with alterations in gene content, gene number, gene arrangement, gene expression and transposon activity-may generate sufficient novelty that every individual in a polyploid population or species may be unique. Whereas certain combinations of these features will undoubtedly be maladaptive, some unique combinations of newly generated variation may provide tremendous evolutionary potential and adaptive capabilities.

The causes and molecular consequences of polyploidy in flowering plants

Polyploidy is an important force shaping plant genomes. All flowering plants are descendants of an ancestral polyploid species, and up to 70% of extant vascular plant species are believed to be recent polyploids. Over the past century, a significant body of knowledge has accumulated regarding the prevalence and ecology of polyploid plants. In this review, we summarize our current understanding of the causes and molecular consequences of polyploidization in angiosperms. We also provide a discussion on the relationships between polyploidy and adaptation and suggest areas where further research may provide a better understanding of polyploidy.

Contrasting evolutionary trajectories of multiple retrotransposons following independent allopolyploidy in wild wheats

Transposable elements (TEs) are expectedly central to genome evolution. To assess the impact of TEs in driving genome turnover, we used allopolyploid genomes, showing considerable deviation from the predicted additivity of their diploid progenitors and thus having undergone major restructuring. Genome survey sequencing was used to select 17 putatively active families of long terminal repeat retrotransposons. Genome-wide TE insertions were genotyped with sequence-specific amplified polymorphism (SSAP) in diploid progenitors and their derived polyploids, and compared with changes in random sequences to assess restructuring of four independent Aegilops allotetraploid genomes. Generally, TEs with different evolutionary trajectories from those of random sequences were identified. Thus, TEs presented family-specific and species-specific dynamics following polyploidy, as illustrated by Sabine showing proliferation in particular polyploids, but massive elimination in others. Contrasting with that, only a few families (BARE1 and Romani) showed proliferation in all polyploids. Overall, TE divergence between progenitors was strongly correlated with the degree of restructuring in polyploid TE fractions. TE families present evolutionary trajectories that are decoupled from genome-wide changes after allopolyploidy and have a pervasive impact on their restructuring.

Evolutionary dynamics of retrotransposons assessed by high-throughput sequencing in wild relatives of wheat

Transposable elements (TEs) represent a major fraction of plant genomes and drive their evolution. An improved understanding of genome evolution requires the dynamics of a large number of TE families to be considered. We put forward an approach bypassing the required step of a complete reference genome to assess the evolutionary trajectories of high copy number TE families from genome snapshot with high-throughput sequencing. Low coverage sequencing of the complex genomes of Aegilops cylindrica and Ae. geniculata using 454 identified more than 70% of the sequences as known TEs, mainly long terminal repeat (LTR) retrotransposons. Comparing the abundance of reads as well as patterns of sequence diversity and divergence within and among genomes assessed the dynamics of 44 major LTR retrotransposon families of the 165 identified. In particular, molecular population genetics on individual TE copies distinguished recently active from quiescent families and highlighted different evolutionary trajectories of retrotransposons among related species. This work presents a suite of tools suitable for current sequencing data, allowing to address the genome-wide evolutionary dynamics of TEs at the family level and advancing our understanding of the evolution of nonmodel genomes.