LeafGo: Leaf to Genome, a quick workflow to produce high-quality de novo plant genomes using long-read sequencing technology

Chromosome-level haplotype-resolved genome assembly of bread wheat's wild relative Aegilops mutica

Bread wheat (Triticum aestivum) is a vital staple crop, with an urgent need for increased production to help feed the world's growing population. Aegilops mutica (2n = 2x = 14; T genome) is a diploid wild relative of wheat carrying valuable agronomic traits resulting in its extensive exploitation for wheat improvement. This paper reports a chromosome-scale, haplotype-resolved genome assembly of Ae. mutica using HiFi reads and Omni-C data. The final lengths for the curated genomes were ~4.65 Gb (haplotype 1) and 4.56 Gb (haplotype 2), featuring a contig N50 of ~4.35 Mb and ~4.60 Mb, respectively. Genome annotation predicted 96,723 gene models and repeats. In summary, the genome assembly of Ae. mutica provides a valuable resource for the wheat breeding community, facilitating faster and more efficient pre-breeding of wheat to enhance food security.

Third generation sequencing transforming plant genome research: Current trends and challenges

In recent years, third-generation sequencing (TGS) technologies have transformed genomics and transcriptomics research, providing novel opportunities for significant discoveries. The long-read sequencing platforms, with their unique advantages over next-generation sequencing (NGS), including a definitive protocol, reduced operational time, and real-time sequencing, possess the potential to transform plant genomics. TGS optimizes and enhances the efficiency of data analysis by removing the necessity for time-consuming assembly tools. The current review examines the development and application of bioinformatics tools for data analysis and annotation, driven by the rapid advancement of TGS platforms like Oxford Nanopore Technologies and Pacific Biosciences. Transcriptome analysis utilizing TGS has been extensively employed to elucidate complex plant transcriptomes and genomes, particularly those characterized by high frequencies of duplicated genomes and repetitive sequences. As a result, current methodologies that allow for generating transcriptomes and comprehensive whole-genome sequences of complex plant genomes employing tailored hybrid sequencing techniques that integrate NGS and TGS technologies have been emphasized herein. This paper, thus, articulates a vision for a future in which TGS effectively addresses the challenges faced in plant research, offering a comprehensive understanding of its advantages, applications, limitations, and promising prospects.

Chromosome-level genome assembly and annotation of largemouth bronze gudgeon (Coreius guichenoti)

Coreius guichenoti, mainly distributed in upstream regions of the Yangtze River China, is currently on the brink of extinction and listed as national secondary protected animal. In this study, we aimed to obtain the chromosome-level genome of C. guichenoti using PacBio and Hi-C techniques. According to the PacBio sequencing, C. guichenoti genome was successfully assembled to 1100.1 Mb size, with a Contig N50 size of 25.0 Mb, and containing 731.0 Mb of repeats. Hi-C sequencing data was utilized for chromosome assembly and 25 chromosome sequences were ultimately yielded, with a total length of 1076.8 Mb. Moreover, a total of 22,506 protein-coding genes were predicted with average intron length of 2293 bp. Evolutionary analysis and divergence time prediction revealed that C. guichenoti was closely related to C. heterodon and they phylogenetically diverged from common ancestor ~20.7 million years ago (Mya), following the separation of Cyprinidae at 28.3 Mya. In the future, the utilization of comparative genomics research is important in elucidating the molecular mechanisms of Ichthyophthirius disease and ensuring the conservation of biological resources.

Genomic Variability Survey in Ilex aquifolium L., with Reference to Four Insular Populations from Eastern Europe

Cosmopolitan in the western areas of Europe as well as on other continents, the Ilex genus is interesting for its genetic, phenotypic, and biogeographic variabilities. Its insular/local distribution, according to existing data on the periphery of the central and southern European areas, represents a suitable case study with reference to the adaptive plasticity or acclimatization of the Ilex aquifolium L. species to new climatic conditions. The aim of the present study was to analyze the genetic variability at the genome level in four insular populations of Ilex aquifolium L., i.e., in three spontaneous populations from Romania (RO), Serbia (SR), and Bulgaria (BG) and a cultivated population from Hungary (HU). According to the obtained results, the most genetically similar populations among the four considered in this study were those from SR and RO. Genetic variation overlapped genes that were generally associated with metabolic regulation/transport factors, water, and abiotic stress factors. The analysis of single-nucleotide polymorphisms (SNPs) at the levels of the chloroplast and mitochondrion, from the point of view of their distributions at the gene level, identified two clusters: one that includes the native populations (BG, SR, and RO) and a second one including the cultured population from HU.

Genome assembly of a diversity panel of Chenopodium quinoa

Quinoa (Chenopodium quinoa) is an important crop for the future challenges of food and nutrient security. Deep characterization of quinoa diversity is needed to support the agronomic improvement and adaptation of quinoa as its worldwide cultivation expands. In this study, we report the construction of chromosome-scale genome assemblies of eight quinoa accessions covering the range of phenotypic and genetic diversity of both lowland and highland quinoas. The assemblies were produced from a combination of PacBio HiFi reads and Bionano Saphyr optical maps, with total assembly sizes averaging 1.28 Gb with a mean N50 of 71.1 Mb. Between 43,733 and 48,564 gene models were predicted for the eight new quinoa genomes, and on average, 66% of each quinoa genome was classified as repetitive sequences. Alignment between the eight genome assemblies allowed the identification of structural rearrangements including inversions, translocations, and duplications. These eight novel quinoa genome assemblies provide a resource for association genetics, comparative genomics, and pan-genome analyses for the discovery of genetic components and variations underlying agriculturally important traits.

Relics of interspecific hybridization retained in the genome of a drought-adapted peanut cultivar

Peanut (Arachis hypogaea L.) is a globally important oil and food crop frequently grown in arid, semi-arid, or dryland environments. Improving drought tolerance is a key goal for peanut crop improvement efforts. Here, we present the genome assembly and gene model annotation for "Line8," a peanut genotype bred from drought-tolerant cultivars. Our assembly and annotation are the most contiguous and complete peanut genome resources currently available. The high contiguity of the Line8 assembly allowed us to explore structural variation both between peanut genotypes and subgenomes. We detect several large inversions between Line8 and other peanut genome assemblies, and there is a trend for the inversions between more genetically diverged genotypes to have higher gene content. We also relate patterns of subgenome exchange to structural variation between Line8 homeologous chromosomes. Unexpectedly, we discover that Line8 harbors an introgression from A.cardenasii, a diploid peanut relative and important donor of disease resistance alleles to peanut breeding populations. The fully resolved sequences of both haplotypes in this introgression provide the first in situ characterization of A.cardenasii candidate alleles that can be leveraged for future targeted improvement efforts. The completeness of our genome will support peanut biotechnology and broader research into the evolution of hybridization and polyploidy.

Expanding the Triangle of U: Comparative analysis of the Hirschfeldia incana genome provides insights into chromosomal evolution, phylogenomics and high photosynthesis-related traits

BACKGROUND AND AIMS

The Brassiceae tribe encompasses many economically important crops and exhibits high intraspecific and interspecific phenotypic variation. After a shared whole-genome triplication (WGT) event (Br-α, ~15.9 million years ago), differential lineage diversification and genomic changes contributed to an array of divergence in morphology, biochemistry, and physiology underlying photosynthesis-related traits. Here, the C3 species Hirschfeldia incana is studied as it displays high photosynthetic rates under high-light conditions. Our aim was to elucidate the evolution that gave rise to the genome of H. incana and its high-photosynthesis traits.

METHODS

We reconstructed a chromosome-level genome assembly for H. incana (Nijmegen, v2.0) using nanopore and chromosome conformation capture (Hi-C) technologies, with 409Mb in size and an N50 of 52Mb (a 10× improvement over the previously published scaffold-level v1.0 assembly). The updated assembly and annotation was subsequently employed to investigate the WGT history of H. incana in a comparative phylogenomic framework from the Brassiceae ancestral genomic blocks and related diploidized crops.

KEY RESULTS

Hirschfeldia incana (x=7) shares extensive genome collinearity with Raphanus sativus (x=9). These two species share some commonalities with Brassica rapa and B. oleracea (A genome, x=10 and C genome, x=9, respectively) and other similarities with B. nigra (B genome, x=8). Phylogenetic analysis revealed that H. incana and R. sativus form a monophyletic clade in between the Brassica A/C and B genomes. We postulate that H. incana and R. sativus genomes are results of hybridization or introgression of the Brassica A/C and B genome types. Our results might explain the discrepancy observed in published studies regarding phylogenetic placement of H. incana and R. sativus in relation to the "Triangle of U" species. Expression analysis of WGT retained gene copies revealed sub-genome expression divergence, likely due to neo- or sub-functionalization. Finally, we highlighted genes associated with physio-biochemical-anatomical adaptive changes observed in H. incana which likely facilitate its high-photosynthesis traits under high light.

CONCLUSIONS

The improved H. incana genome assembly, annotation and results presented in this work will be a valuable resource for future research to unravel the genetic basis of its ability to maintain a high photosynthetic efficiency in high-light conditions and thereby improve photosynthesis for enhanced agricultural production.

Unlocking plant genetics with telomere-to-telomere genome assemblies

Contiguous genome sequence assemblies will help us to realize the full potential of crop translational genomics. Recent advances in sequencing technologies, especially long-read sequencing strategies, have made it possible to construct gapless telomere-to-telomere (T2T) assemblies, thus offering novel insights into genome organization and function. Plant genomes pose unique challenges, such as a continuum of ancient to recent polyploidy and abundant highly similar and long repetitive elements. Owing to progress in sequencing approaches, for most crop plants, chromosome-scale reference genome assemblies are available, but T2T assembly construction remains challenging. Here we describe methods for haplotype-resolved, gapless T2T assembly construction in plants, including various crop species. We outline the impact of T2T assemblies in elucidating the roles of repetitive elements in gene regulation, as well as in pangenomics, functional genomics, genome-assisted breeding and targeted genome manipulation. In conjunction with sequence-enriched germplasm repositories, T2T assemblies thus hold great promise for basic and applied plant sciences.

Origin and evolution of the bread wheat D genome

Bread wheat (Triticum aestivum) is a globally dominant crop and major source of calories and proteins for the human diet. Compared with its wild ancestors, modern bread wheat shows lower genetic diversity, caused by polyploidisation, domestication and breeding bottlenecks1,2. Wild wheat relatives represent genetic reservoirs, and harbour diversity and beneficial alleles that have not been incorporated into bread wheat. Here we establish and analyse extensive genome resources for Tausch's goatgrass (Aegilops tauschii), the donor of the bread wheat D genome. Our analysis of 46 Ae. tauschii genomes enabled us to clone a disease resistance gene and perform haplotype analysis across a complex disease resistance locus, allowing us to discern alleles from paralogous gene copies. We also reveal the complex genetic composition and history of the bread wheat D genome, which involves contributions from genetically and geographically discrete Ae. tauschii subpopulations. Together, our results reveal the complex history of the bread wheat D genome and demonstrate the potential of wild relatives in crop improvement.

Chromosome-level genome assembly of Cnidium monnieri, a highly demanded traditional Chinese medicine

Cnidium monnieri, a medicinal herb of the Cnidium genus and the Apiaceae family, is among the most important traditional Chinese medicines and is widely distributed in China. However, to date, no C. monnieri-related genomic information has been described. In this study, we assembled the C. monnieri genome of approximately 1210.23 Mb with a contig N50 of 83.14 Mb. Using PacBio HiFi and Hi-C sequencing data, we successfully anchored 93.86% of the assembled sequences to 10 pseudochromosomes (2n = 20). We predicted a total of 37,460 protein-coding genes, with 97.02% of them being functionally annotated in Non-Redundant, Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, and other databases. In addition, we identified 2,778 tRNAs, 4,180 rRNAs, 258 miRNAs, and 1,700 snRNAs in the genome. This is the first reported C. monnieri genome. Hopefully, the availability of this chromosome-level reference genome provides a significant basis for upcoming natural product-related biosynthetic pathway assessment in C. monnieri.

Genomic, transcriptomic, and metabolomic analyses provide insights into the evolution and development of a medicinal plant Saposhnikovia divaricata (Apiaceae)

Saposhnikovia divaricata, 2n = 2x = 16, as a perennial species, is widely distributed in China, Mongolia, Russia, etc. It is a traditional Chinese herb used to treat tetanus, rubella pruritus, rheumatic arthralgia, and other diseases. Here, we assembled a 2.07 Gb and N50 scaffold length of 227.67 Mb high-quality chromosome-level genome of S. divaricata based on the PacBio Sequel II sequencing platform. The total number of genes identified was 42 948, and 42 456 of them were functionally annotated. A total of 85.07% of the genome was composed of repeat sequences, comprised mainly of long terminal repeats (LTRs) which represented 73.7% of the genome sequence. The genome size may have been affected by a recent whole-genome duplication event. Transcriptional and metabolic analyses revealed bolting and non-bolting S. divaricata differed in flavonoids, plant hormones, and some pharmacologically active components. The analysis of its genome, transcriptome, and metabolome helped to provide insights into the evolution of bolting and non-bolting phenotypes in wild and cultivated S. divaricata and lays the basis for genetic improvement of the species.

Chromosome-scale genome assembly of bread wheat's wild relative Triticum timopheevii

Wheat (Triticum aestivum) is one of the most important food crops with an urgent need for increase in its production to feed the growing world. Triticum timopheevii (2n = 4x = 28) is an allotetraploid wheat wild relative species containing the At and G genomes that has been exploited in many pre-breeding programmes for wheat improvement. In this study, we report the generation of a chromosome-scale reference genome assembly of T. timopheevii accession PI 94760 based on PacBio HiFi reads and chromosome conformation capture (Hi-C). The assembly comprised a total size of 9.35 Gb, featuring a contig N50 of 42.4 Mb and included the mitochondrial and plastid genome sequences. Genome annotation predicted 166,325 gene models including 70,365 genes with high confidence. DNA methylation analysis showed that the G genome had on average more methylated bases than the At genome. In summary, the T. timopheevii genome assembly provides a valuable resource for genome-informed discovery of agronomically important genes for food security.

Plasticity of repetitive sequences demonstrated by the complete mitochondrial genome of Eucalyptus camaldulensis

The tree Eucalyptus camaldulensis is a ubiquitous member of the Eucalyptus genus, which includes several hundred species. Despite the extensive sequencing and assembly of nuclear genomes from various eucalypts, the genus has only one fully annotated and complete mitochondrial genome (mitogenome). Plant mitochondria are characterized by dynamic genomic rearrangements, facilitated by repeat content, a feature that has hindered the assembly of plant mitogenomes. This complexity is evident in the paucity of available mitogenomes. This study, to the best of our knowledge, presents the first E. camaldulensis mitogenome. Our findings suggest the presence of multiple isomeric forms of the E. camaldulensis mitogenome and provide novel insights into minor rearrangements triggered by nested repeat sequences. A comparative sequence analysis of the E. camaldulensis and E. grandis mitogenomes unveils evolutionary changes between the two genomes. A significant divergence is the evolution of a large repeat sequence, which may have contributed to the differences observed between the two genomes. The largest repeat sequences in the E. camaldulensis mitogenome align well with significant yet unexplained structural variations in the E. grandis mitogenome, highlighting the adaptability of repeat sequences in plant mitogenomes.

A haplotype resolved chromosome-scale assembly of North American wild apple Malus fusca and comparative genomics of the fire blight Mfu10 locus

The Pacific crabapple (Malus fusca) is a wild relative of the commercial apple (Malus × domestica). With a range extending from Alaska to Northern California, M. fusca is extremely hardy and disease resistant. The species represents an untapped genetic resource for the development of new apple cultivars with enhanced stress resistance. However, gene discovery and utilization of M. fusca have been hampered by the lack of genomic resources. Here, we present a high‐quality, haplotype‐resolved, chromosome‐scale genome assembly and annotation for M. fusca. The genome was assembled using high‐fidelity long‐reads and scaffolded using genetic maps and high‐throughput chromatin conformation capture sequencing, resulting in one of the most contiguous apple genomes to date. We annotated the genome using public transcriptomic data from the same species taken from diverse plant structures and developmental stages. Using this assembly, we explored haplotypic structural variation within the genome of M. fusca, identifying thousands of large variants. We further showed high sequence co‐linearity with other domesticated and wild Malus species. Finally, we resolve a known quantitative trait locus associated with resistance to fire blight (Erwinia amylovora). Insights gained from the assembly of a reference‐quality genome of this hardy wild apple relative will be invaluable as a tool to facilitate DNA‐informed introgression breeding.

Chromosome-scale assembly of the wild wheat relative Aegilops umbellulata

Wild wheat relatives have been explored in plant breeding to increase the genetic diversity of bread wheat, one of the most important food crops. Aegilops umbellulata is a diploid U genome-containing grass species that serves as a genetic reservoir for wheat improvement. In this study, we report the construction of a chromosome-scale reference assembly of Ae. umbellulata accession TA1851 based on corrected PacBio HiFi reads and chromosome conformation capture. The total assembly size was 4.25 Gb with a contig N50 of 17.7 Mb. In total, 36,268 gene models were predicted. We benchmarked the performance of hifiasm and LJA, two of the most widely used assemblers using standard and corrected HiFi reads, revealing a positive effect of corrected input reads. Comparative genome analysis confirmed substantial chromosome rearrangements in Ae. umbellulata compared to bread wheat. In summary, the Ae. umbellulata assembly provides a resource for comparative genomics in Triticeae and for the discovery of agriculturally important genes.

An unusual tandem kinase fusion protein confers leaf rust resistance in wheat

The introgression of chromosome segments from wild relatives is an established strategy to enrich crop germplasm with disease-resistance genes¹. Here we use mutagenesis and transcriptome sequencing to clone the leaf rust resistance gene Lr9, which was introduced into bread wheat from the wild grass species Aegilops umbellulata². We established that Lr9 encodes an unusual tandem kinase fusion protein. Long-read sequencing of a wheat Lr9 introgression line and the putative Ae. umbellulata Lr9 donor enabled us to assemble the ~28.4-Mb Lr9 translocation and to identify the translocation breakpoint. We likewise cloned Lr58, which was reportedly introgressed from Aegilops triuncialis³, but has an identical coding sequence compared to Lr9. Cytogenetic and haplotype analyses corroborate that the two genes originate from the same translocation event. Our work sheds light on the emerging role of kinase fusion proteins in wheat disease resistance, expanding the repertoire of disease-resistance genes for breeding.

Isolating Linum usitatissimum L. Nuclear DNA Enabled Assembling High-Quality Genome

High-quality genome sequences help to elucidate the genetic basis of numerous biological processes and track species evolution. For flax (Linum usitatissimum L.)-a multifunctional crop, high-quality assemblies from Oxford Nanopore Technologies (ONT) data were unavailable, largely due to the difficulty of isolating pure high-molecular-weight DNA. This article proposes a scheme for gaining a contiguous L. usitatissimum assembly using Nanopore data. We developed a protocol for flax nuclei isolation with subsequent DNA extraction, which allows obtaining about 5 μg of pure high-molecular-weight DNA from 0.5 g of leaves. Such an amount of material can be collected even from a single plant and yields more than 30 Gb of ONT data in two MinION runs. We performed a comparative analysis of different genome assemblers and polishers on the gained data and obtained the final 447.1-Mb assembly of L. usitatissimum line 3896 genome using the Canu-Racon (two iterations)-Medaka combination. The genome comprised 1695 contigs and had an N50 of 6.2 Mb and a completeness of 93.8% of BUSCOs from eudicots_odb10. Our study highlights the impact of the chosen genome construction strategy on the resulting assembly parameters and its eligibility for future genomic studies.

Establishing MinION Sequencing and Genome Assembly Procedures for the Analysis of the Rooibos (Aspalathus linearis) Genome

While plant genome analysis is gaining speed worldwide, few plant genomes have been sequenced and analyzed on the African continent. Yet, this information holds the potential to transform diverse industries as it unlocks medicinally and industrially relevant biosynthesis pathways for bioprospecting. Considering that South Africa is home to the highly diverse Cape Floristic Region, local establishment of methods for plant genome analysis is essential. Long-read sequencing is becoming standard procedure for plant genome research, as these reads can span repetitive regions of the DNA, substantially facilitating reassembly of a contiguous genome. With the MinION, Oxford Nanopore offers a cost-efficient sequencing method to generate long reads; however, DNA purification protocols must be adapted for each plant species to generate ultra-pure DNA, essential for these analyses. Here, we describe a cost-effective procedure for the extraction and purification of plant DNA and evaluate diverse genome assembly approaches for the reconstruction of the genome of rooibos (Aspalathus linearis), an endemic South African medicinal plant widely used for tea production. We discuss the pros and cons of nine tested assembly programs, specifically Redbean and NextDenovo, which generated the most contiguous assemblies, and Flye, which produced an assembly closest to the predicted genome size.

A High-Quality Genome Assembly of Striped Catfish (Pangasianodon hypophthalmus) Based on Highly Accurate Long-Read HiFi Sequencing Data

The HiFi sequencing technology yields highly accurate long-read data with accuracies greater than 99.9% that can be used to improve results for complex applications such as genome assembly. Our study presents a high-quality chromosome-scale genome assembly of striped catfish (Pangasianodon hypophthalmus), a commercially important species cultured mainly in Vietnam, integrating HiFi reads and Hi-C data. A 788.4 Mb genome containing 381 scaffolds with an N50 length of 21.8 Mb has been obtained from HiFi reads. These scaffolds have been further ordered and clustered into 30 chromosome groups, ranging from 1.4 to 57.6 Mb, based on Hi-C data. The present updated assembly has a contig N50 of 14.7 Mb, representing a 245-fold and 4.2-fold improvement over the previous Illumina and Illumina-Nanopore-Hi-C based version, respectively. In addition, the proportion of repeat elements and BUSCO genes identified in our genome is remarkably higher than in the two previously released striped catfish genomes. These results highlight the power of using HiFi reads to assemble the highly repetitive regions and to improve the quality of genome assembly. The updated, high-quality genome assembled in this work will provide a valuable genomic resource for future population genetics, conservation biology and selective breeding studies of striped catfish.

De novo chromosome level assembly of a plant genome from long read sequence data

Recent advances in the sequencing and assembly of plant genomes have allowed the generation of genomes with increasing contiguity and sequence accuracy. Chromosome level genome assemblies using sequence contigs generated from long read sequencing have involved the use of proximity analysis (Hi-C) or traditional genetic maps to guide the placement of sequence contigs within chromosomes. The development of highly accurate long reads by repeated sequencing of circularized DNA (HiFi; PacBio) has greatly increased the size of contigs. We now report the use of HiFiasm to assemble the genome of Macadamia jansenii, a genome that has been used as a model to test sequencing and assembly. This achieved almost complete chromosome level assembly from the sequence data alone without the need for higher level chromosome map information. Eight of the 14 chromosomes were represented by a single large contig (six with telomere repeats at both ends) and the other six assembled from two to four main contigs. The small number of chromosome breaks appears to be the result of highly repetitive regions including ribosomal genes that cannot be assembled by these approaches. De novo assembly of near complete chromosome level plant genomes now appears possible using these sequencing and assembly tools. Further targeted strategies might allow these remaining gaps to be closed.

Green plant genomes: What we know in an era of rapidly expanding opportunities

Green plants play a fundamental role in ecosystems, human health, and agriculture. As de novo genomes are being generated for all known eukaryotic species as advocated by the Earth BioGenome Project, increasing genomic information on green land plants is essential. However, setting standards for the generation and storage of the complex set of genomes that characterize the green lineage of life is a major challenge for plant scientists. Such standards will need to accommodate the immense variation in green plant genome size, transposable element content, and structural complexity while enabling research into the molecular and evolutionary processes that have resulted in this enormous genomic variation. Here we provide an overview and assessment of the current state of knowledge of green plant genomes. To date fewer than 300 complete chromosome-scale genome assemblies representing fewer than 900 species have been generated across the estimated 450,000 to 500,000 species in the green plant clade. These genomes range in size from 12 Mb to 27.6 Gb and are biased toward agricultural crops with large branches of the green tree of life untouched by genomic-scale sequencing. Locating suitable tissue samples of most species of plants, especially those taxa from extreme environments, remains one of the biggest hurdles to increasing our genomic inventory. Furthermore, the annotation of plant genomes is at present undergoing intensive improvement. It is our hope that this fresh overview will help in the development of genomic quality standards for a cohesive and meaningful synthesis of green plant genomes as we scale up for the future.

Comparing whole-genome shotgun sequencing and DNA metabarcoding approaches for species identification and quantification of pollen species mixtures

Molecular identification of mixed-species pollen samples has a range of applications in various fields of research. To date, such molecular identification has primarily been carried out via amplicon sequencing, but whole-genome shotgun (WGS) sequencing of pollen DNA has potential advantages, including (1) more genetic information per sample and (2) the potential for better quantitative matching. In this study, we tested the performance of WGS sequencing methodology and publicly available reference sequences in identifying species and quantifying their relative abundance in pollen mock communities. Using mock communities previously analyzed with DNA metabarcoding, we sequenced approximately 200Mbp for each sample using Illumina HiSeq and MiSeq. Taxonomic identifications were based on the Kraken k-mer identification method with reference libraries constructed from full-genome and short read archive data from the NCBI database. We found WGS to be a reliable method for taxonomic identification of pollen with near 100% identification of species in mixtures but generating higher rates of false positives (reads not identified to the correct taxon at the required taxonomic level) relative to rbcL and ITS2 amplicon sequencing. For quantification of relative species abundance, WGS data provided a stronger correlation between pollen grain proportion and sequence read proportion, but diverged more from a 1:1 relationship, likely due to the higher rate of false positives. Currently, a limitation of WGS-based pollen identification is the lack of representation of plant diversity in publicly available genome databases. As databases improve and costs drop, we expect that eventually genomics methods will become the methods of choice for species identification and quantification of mixed-species pollen samples.