Cardiac function in a large animal model of myocardial infarction at 7 T: deep learning based automatic segmentation increases reproducibility

Cardiac magnetic resonance (CMR) imaging allows precise non-invasive quantification of cardiac function. It requires reliable image segmentation for myocardial tissue. Clinically used software usually offers automatic approaches for this step. These are, however, designed for segmentation of human images obtained at clinical field strengths. They reach their limits when applied to preclinical data and ultrahigh field strength (such as CMR of pigs at 7 T). In our study, eleven animals (seven with myocardial infarction) underwent four CMR scans each. Short-axis cine stacks were acquired and used for functional cardiac analysis. End-systolic and end-diastolic images were labelled manually by two observers and inter- and intra-observer variability were assessed. Aiming to make the functional analysis faster and more reproducible, an established deep learning (DL) model for myocardial segmentation in humans was re-trained using our preclinical 7 T data (n = 772 images and labels). We then tested the model on n = 288 images. Excellent agreement in parameters of cardiac function was found between manual and DL segmentation: For ejection fraction (EF) we achieved a Pearson's r of 0.95, an Intraclass correlation coefficient (ICC) of 0.97, and a Coefficient of variability (CoV) of 6.6%. Dice scores were 0.88 for the left ventricle and 0.84 for the myocardium.

Semi-automated sequence curation for reliable reference datasets in ITS2 vascular plant DNA (meta-)barcoding

One of the most critical steps for accurate taxonomic identification in DNA (meta)-barcoding is to have an accurate DNA reference sequence dataset for the marker of choice. Therefore, developing such a dataset has been a long-term ambition, especially in the Viridiplantae kingdom. Typically, reference datasets are constructed with sequences downloaded from general public databases, which can carry taxonomic and other relevant errors. Herein, we constructed a curated (i) global dataset, (ii) European crop dataset, and (iii) 27 datasets for the EU countries for the ITS2 barcoding marker of vascular plants. To that end, we first developed a pipeline script that entails (i) an automated curation stage comprising five filters, (ii) manual taxonomic correction for misclassified taxa, and (iii) manual addition of newly sequenced species. The pipeline allows easy updating of the curated datasets. With this approach, 13% of the sequences, corresponding to 7% of species originally imported from GenBank, were discarded. Further, 259 sequences were manually added to the curated global dataset, which now comprises 307,977 sequences of 111,382 plant species.

Opportunistic Bacteria of Grapevine Crown Galls Are Equipped with the Genomic Repertoire for Opine Utilization

Young grapevines (Vitis vinifera) suffer and eventually can die from the crown gall disease caused by the plant pathogen Allorhizobium vitis (Rhizobiaceae). Virulent members of A. vitis harbor a tumor-inducing plasmid and induce formation of crown galls due to the oncogenes encoded on the transfer DNA. The expression of oncogenes in transformed host cells induces unregulated cell proliferation and metabolic and physiological changes. The crown gall produces opines uncommon to plants, which provide an important nutrient source for A. vitis harboring opine catabolism enzymes. Crown galls host a distinct bacterial community, and the mechanisms establishing a crown gall-specific bacterial community are currently unknown. Thus, we were interested in whether genes homologous to those of the tumor-inducing plasmid coexist in the genomes of the microbial species coexisting in crown galls. We isolated 8 bacterial strains from grapevine crown galls, sequenced their genomes, and tested their virulence and opine utilization ability in bioassays. In addition, the 8 genome sequences were compared with 34 published bacterial genomes, including closely related plant-associated bacteria not from crown galls. Homologous genes for virulence and opine anabolism were only present in the virulent Rhizobiaceae. In contrast, homologs of the opine catabolism genes were present in all strains including the nonvirulent members of the Rhizobiaceae and non-Rhizobiaceae. Gene neighborhood and sequence identity of the opine degradation cluster of virulent and nonvirulent strains together with the results of the opine utilization assay support the important role of opine utilization for cocolonization in crown galls, thereby shaping the crown gall community.

Recognition and reconstruction of cell differentiation patterns with deep learning

Cell lineage decisions occur in three-dimensional spatial patterns that are difficult to identify by eye. There is an ongoing effort to replicate such patterns using mathematical modeling. One approach uses long ranging cell-cell communication to replicate common spatial arrangements like checkerboard and engulfing patterns. In this model, the cell-cell communication has been implemented as a signal that disperses throughout the tissue. On the other hand, machine learning models have been developed for pattern recognition and pattern reconstruction tasks. We combined synthetic data generated by the mathematical model with spatial summary statistics and deep learning algorithms to recognize and reconstruct cell fate patterns in organoids of mouse embryonic stem cells. Application of Moran's index and pair correlation functions for in vitro and synthetic data from the model showed local clustering and radial segregation. To assess the patterns as a whole, a graph neural network was developed and trained on synthetic data from the model. Application to in vitro data predicted a low signal dispersion value. To test this result, we implemented a multilayer perceptron for the prediction of a given cell fate based on the fates of the neighboring cells. The results show a 70% accuracy of cell fate imputation based on the nine nearest neighbors of a cell. Overall, our approach combines deep learning with mathematical modeling to link cell fate patterns with potential underlying mechanisms.

MyoPS: A benchmark of myocardial pathology segmentation combining three-sequence cardiac magnetic resonance images

Assessment of myocardial viability is essential in diagnosis and treatment management of patients suffering from myocardial infarction, and classification of pathology on the myocardium is the key to this assessment. This work defines a new task of medical image analysis, i.e., to perform myocardial pathology segmentation (MyoPS) combining three-sequence cardiac magnetic resonance (CMR) images, which was first proposed in the MyoPS challenge, in conjunction with MICCAI 2020. Note that MyoPS refers to both myocardial pathology segmentation and the challenge in this paper. The challenge provided 45 paired and pre-aligned CMR images, allowing algorithms to combine the complementary information from the three CMR sequences for pathology segmentation. In this article, we provide details of the challenge, survey the works from fifteen participants and interpret their methods according to five aspects, i.e., preprocessing, data augmentation, learning strategy, model architecture and post-processing. In addition, we analyze the results with respect to different factors, in order to examine the key obstacles and explore the potential of solutions, as well as to provide a benchmark for future research. The average Dice scores of submitted algorithms were 0.614±0.231 and 0.644±0.153 for myocardial scars and edema, respectively. We conclude that while promising results have been reported, the research is still in the early stage, and more in-depth exploration is needed before a successful application to the clinics. MyoPS data and evaluation tool continue to be publicly available upon registration via its homepage (www.sdspeople.fudan.edu.cn/zhuangxiahai/0/myops20/).

Corrigendum: Honeybees are buffered against undernourishment during larval stages

[This corrects the article DOI: 10.3389/finsc.2022.951317.].

Novel integrative elements and genomic plasticity in ocean ecosystems

Horizontal gene transfer accelerates microbial evolution. The marine picocyanobacterium Prochlorococcus exhibits high genomic plasticity, yet the underlying mechanisms are elusive. Here, we report a novel family of DNA transposons-"tycheposons"-some of which are viral satellites while others carry cargo, such as nutrient-acquisition genes, which shape the genetic variability in this globally abundant genus. Tycheposons share distinctive mobile-lifecycle-linked hallmark genes, including a deep-branching site-specific tyrosine recombinase. Their excision and integration at tRNA genes appear to drive the remodeling of genomic islands-key reservoirs for flexible genes in bacteria. In a selection experiment, tycheposons harboring a nitrate assimilation cassette were dynamically gained and lost, thereby promoting chromosomal rearrangements and host adaptation. Vesicles and phage particles harvested from seawater are enriched in tycheposons, providing a means for their dispersal in the wild. Similar elements are found in microbes co-occurring with Prochlorococcus, suggesting a common mechanism for microbial diversification in the vast oligotrophic oceans.

Honeybees are buffered against undernourishment during larval stages

The negative impact of juvenile undernourishment on adult behavior has been well reported for vertebrates, but relatively little is known about invertebrates. In honeybees, nutrition has long been known to affect task performance and timing of behavioral transitions. Whether and how a dietary restriction during larval development affects the task performance of adult honeybees is largely unknown. We raised honeybees in-vitro, varying the amount of a standardized diet (150 µl, 160 µl, 180 µl in total). Emerging adults were marked and inserted into established colonies. Behavioral performance of nurse bees and foragers was investigated and physiological factors known to be involved in the regulation of social organization were quantified. Surprisingly, adult honeybees raised under different feeding regimes did not differ in any of the behaviors observed. No differences were observed in physiological parameters apart from weight. Honeybees were lighter when undernourished (150 µl), while they were heavier under the overfed treatment (180 µl) compared to the control group raised under a normal diet (160 µl). These data suggest that dietary restrictions during larval development do not affect task performance or physiology in this social insect despite producing clear effects on adult weight. We speculate that possible effects of larval undernourishment might be compensated during the early period of adult life.

Efficient permutation-based genome-wide association studies for normal and skewed phenotypic distributions

MOTIVATION

Genome-wide association studies (GWAS) are an integral tool for studying the architecture of complex genotype and phenotype relationships. Linear mixed models (LMMs) are commonly used to detect associations between genetic markers and a trait of interest, while at the same time allowing to account for population structure and cryptic relatedness. Assumptions of LMMs include a normal distribution of the residuals and that the genetic markers are independent and identically distributed-both assumptions are often violated in real data. Permutation-based methods can help to overcome some of these limitations and provide more realistic thresholds for the discovery of true associations. Still, in practice, they are rarely implemented due to the high computational complexity.

RESULTS

We propose permGWAS, an efficient LMM reformulation based on 4D tensors that can provide permutation-based significance thresholds. We show that our method outperforms current state-of-the-art LMMs with respect to runtime and that permutation-based thresholds have lower false discovery rates for skewed phenotypes compared to the commonly used Bonferroni threshold. Furthermore, using permGWAS we re-analyzed more than 500 Arabidopsis thaliana phenotypes with 100 permutations each in less than 8 days on a single GPU. Our re-analyses suggest that applying a permutation-based threshold can improve and refine the interpretation of GWAS results.

AVAILABILITY AND IMPLEMENTATION

permGWAS is open-source and publicly available on GitHub for download: https://github.com/grimmlab/permGWAS.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Self-configuring nnU-net pipeline enables fully automatic infarct segmentation in late enhancement MRI after myocardial infarction

PURPOSE

To fully automatically derive quantitative parameters from late gadolinium enhancement (LGE) cardiac MR (CMR) in patients with myocardial infarction and to investigate if phase sensitive or magnitude reconstructions or a combination of both results in best segmentation accuracy.

METHODS

In this retrospective single center study, a convolutional neural network with a U-Net architecture with a self-configuring framework ("nnU-net") was trained for segmentation of left ventricular myocardium and infarct zone in LGE-CMR. A database of 170 examinations from 78 patients with history of myocardial infarction was assembled. Separate fitting of the model was performed, using phase sensitive inversion recovery, the magnitude reconstruction or both contrasts as input channels. Manual labelling served as ground truth. In a subset of 10 patients, the performance of the trained models was evaluated and quantitatively compared by determination of the Sørensen-Dice similarity coefficient (DSC) and volumes of the infarct zone compared with the manual ground truth using Pearson's r correlation and Bland-Altman analysis.

RESULTS

The model achieved high similarity coefficients for myocardium and scar tissue. No significant difference was observed between using PSIR, magnitude reconstruction or both contrasts as input (PSIR and MAG; mean DSC: 0.83 ± 0.03 for myocardium and 0.72 ± 0.08 for scars). A strong correlation for volumes of infarct zone was observed between manual and model-based approach (r = 0.96), with a significant underestimation of the volumes obtained from the neural network.

CONCLUSION

The self-configuring nnU-net achieves predictions with strong agreement compared to manual segmentation, proving the potential as a promising tool to provide fully automatic quantitative evaluation of LGE-CMR.

Sensitivity analysis for interpretation of machine learning based segmentation models in cardiac MRI

BACKGROUND

Image segmentation is a common task in medical imaging e.g., for volumetry analysis in cardiac MRI. Artificial neural networks are used to automate this task with performance similar to manual operators. However, this performance is only achieved in the narrow tasks networks are trained on. Performance drops dramatically when data characteristics differ from the training set properties. Moreover, neural networks are commonly considered black boxes, because it is hard to understand how they make decisions and why they fail. Therefore, it is also hard to predict whether they will generalize and work well with new data. Here we present a generic method for segmentation model interpretation. Sensitivity analysis is an approach where model input is modified in a controlled manner and the effect of these modifications on the model output is evaluated. This method yields insights into the sensitivity of the model to these alterations and therefore to the importance of certain features on segmentation performance.

RESULTS

We present an open-source Python library (misas), that facilitates the use of sensitivity analysis with arbitrary data and models. We show that this method is a suitable approach to answer practical questions regarding use and functionality of segmentation models. We demonstrate this in two case studies on cardiac magnetic resonance imaging. The first case study explores the suitability of a published network for use on a public dataset the network has not been trained on. The second case study demonstrates how sensitivity analysis can be used to evaluate the robustness of a newly trained model.

CONCLUSIONS

Sensitivity analysis is a useful tool for deep learning developers as well as users such as clinicians. It extends their toolbox, enabling and improving interpretability of segmentation models. Enhancing our understanding of neural networks through sensitivity analysis also assists in decision making. Although demonstrated only on cardiac magnetic resonance images this approach and software are much more broadly applicable.

A systematic comparison of chloroplast genome assembly tools

BACKGROUND

Chloroplasts are intracellular organelles that enable plants to conduct photosynthesis. They arose through the symbiotic integration of a prokaryotic cell into an eukaryotic host cell and still contain their own genomes with distinct genomic information. Plastid genomes accommodate essential genes and are regularly utilized in biotechnology or phylogenetics. Different assemblers that are able to assess the plastid genome have been developed. These assemblers often use data of whole genome sequencing experiments, which usually contain reads from the complete chloroplast genome.

RESULTS

The performance of different assembly tools has never been systematically compared. Here, we present a benchmark of seven chloroplast assembly tools, capable of succeeding in more than 60% of known real data sets. Our results show significant differences between the tested assemblers in terms of generating whole chloroplast genome sequences and computational requirements. The examination of 105 data sets from species with unknown plastid genomes leads to the assembly of 20 novel chloroplast genomes.

CONCLUSIONS

We create docker images for each tested tool that are freely available for the scientific community and ensure reproducibility of the analyses. These containers allow the analysis and screening of data sets for chloroplast genomes using standard computational infrastructure. Thus, large scale screening for chloroplasts within genomic sequencing data is feasible.

A systematic comparison of chloroplast genome assembly tools

BACKGROUND

Chloroplasts are intracellular organelles that enable plants to conduct photosynthesis. They arose through the symbiotic integration of a prokaryotic cell into an eukaryotic host cell and still contain their own genomes with distinct genomic information. Plastid genomes accommodate essential genes and are regularly utilized in biotechnology or phylogenetics. Different assemblers that are able to assess the plastid genome have been developed. These assemblers often use data of whole genome sequencing experiments, which usually contain reads from the complete chloroplast genome.

RESULTS

The performance of different assembly tools has never been systematically compared. Here, we present a benchmark of seven chloroplast assembly tools, capable of succeeding in more than 60% of known real data sets. Our results show significant differences between the tested assemblers in terms of generating whole chloroplast genome sequences and computational requirements. The examination of 105 data sets from species with unknown plastid genomes leads to the assembly of 20 novel chloroplast genomes.

CONCLUSIONS

We create docker images for each tested tool that are freely available for the scientific community and ensure reproducibility of the analyses. These containers allow the analysis and screening of data sets for chloroplast genomes using standard computational infrastructure. Thus, large scale screening for chloroplasts within genomic sequencing data is feasible.

Open Science principles for accelerating trait-based science across the Tree of Life

Synthesizing trait observations and knowledge across the Tree of Life remains a grand challenge for biodiversity science. Species traits are widely used in ecological and evolutionary science, and new data and methods have proliferated rapidly. Yet accessing and integrating disparate data sources remains a considerable challenge, slowing progress toward a global synthesis to integrate trait data across organisms. Trait science needs a vision for achieving global integration across all organisms. Here, we outline how the adoption of key Open Science principles-open data, open source and open methods-is transforming trait science, increasing transparency, democratizing access and accelerating global synthesis. To enhance widespread adoption of these principles, we introduce the Open Traits Network (OTN), a global, decentralized community welcoming all researchers and institutions pursuing the collaborative goal of standardizing and integrating trait data across organisms. We demonstrate how adherence to Open Science principles is key to the OTN community and outline five activities that can accelerate the synthesis of trait data across the Tree of Life, thereby facilitating rapid advances to address scientific inquiries and environmental issues. Lessons learned along the path to a global synthesis of trait data will provide a framework for addressing similarly complex data science and informatics challenges.

BCdatabaser: on-the-fly reference database creation for (meta-)barcoding

Summary

DNA barcoding and meta-barcoding have become irreplaceable in research and applications, where identification of taxa alone or within a mixture, respectively, becomes relevant. Pioneering studies were in the microbiological context, yet nowadays also plants and animals become targeted. Given the variety of markers used, formatting requirements for classifiers and constant growth of primary databases, there is a need for dedicated reference database creation. We developed a web and command-line interface to generate such on-the-fly for any applicable marker and taxonomic group with optional filtering, formatting and restriction specific for (meta-)barcoding purposes. Also, databases optionally receive a DOI, making them well-documented with meta-data, publicly sharable and citable.

Availability and implementation

source code: https://www.github.com/molbiodiv/bcdatabaser, webservice: https://bcdatabaser.molecular.eco, documentation: https://molbiodiv.github.io/bcdatabaser.

Linking pollen foraging of megachilid bees to their nest bacterial microbiota

Solitary bees build their nests by modifying the interior of natural cavities, and they provision them with food by importing collected pollen. As a result, the microbiota of the solitary bee nests may be highly dependent on introduced materials. In order to investigate how the collected pollen is associated with the nest microbiota, we used metabarcoding of the ITS2 rDNA and the 16S rDNA to simultaneously characterize the pollen composition and the bacterial communities of 100 solitary bee nest chambers belonging to seven megachilid species. We found a weak correlation between bacterial and pollen alpha diversity and significant associations between the composition of pollen and that of the nest microbiota, contributing to the understanding of the link between foraging and bacteria acquisition for solitary bees. Since solitary bees cannot establish bacterial transmission routes through eusociality, this link could be essential for obtaining bacterial symbionts for this group of valuable pollinators.

OPEN RESEARCH BADGES

This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ebi.ac.uk/ena/data/view/PRJEB27223, https://www.ebi.ac.uk/ena/data/view/PRJEB31610, and https://doi.org/10.5061/dryad.qk36k8q.

Wild bees and their nests host Paenibacillus bacteria with functional potential of avail

BACKGROUND

In previous studies, the gram-positive firmicute genus Paenibacillus was found with significant abundances in nests of wild solitary bees. Paenibacillus larvae is well-known for beekeepers as a severe pathogen causing the fatal honey bee disease American foulbrood, and other members of the genus are either secondary invaders of European foulbrood or considered a threat to honey bees. We thus investigated whether Paenibacillus is a common bacterium associated with various wild bees and hence poses a latent threat to honey bees visiting the same flowers.

RESULTS

We collected 202 samples from 82 individuals or nests of 13 bee species at the same location and screened each for Paenibacillus using high-throughput sequencing-based 16S metabarcoding. We then isolated the identified strain Paenibacillus MBD-MB06 from a solitary bee nest and sequenced its genome. We did find conserved toxin genes and such encoding for chitin-binding proteins, yet none specifically related to foulbrood virulence or chitinases. Phylogenomic analysis revealed a closer relationship to strains of root-associated Paenibacillus rather than strains causing foulbrood or other accompanying diseases. We found anti-microbial evidence within the genome, confirmed by experimental bioassays with strong growth inhibition of selected fungi as well as gram-positive and gram-negative bacteria.

CONCLUSIONS

The isolated wild bee associate Paenibacillus MBD-MB06 is a common, but irregularly occurring part of wild bee microbiomes, present on adult body surfaces and guts and within nests especially in megachilids. It was phylogenetically and functionally distinct from harmful members causing honey bee colony diseases, although it shared few conserved proteins putatively toxic to insects that might indicate ancestral predisposition for the evolution of insect pathogens within the group. By contrast, our strain showed anti-microbial capabilities and the genome further indicates abilities for chitin-binding and biofilm-forming, suggesting it is likely a useful associate to avoid fungal penetration of the bee cuticula and a beneficial inhabitant of nests to repress fungal threats in humid and nutrient-rich environments of wild bee nests.

biojs-io-biom, a BioJS component for handling data in Biological Observation Matrix (BIOM) format

The Biological Observation Matrix (BIOM) format is widely used to store data from high-throughput studies. It aims at increasing interoperability of bioinformatic tools that process this data. However, due to multiple versions and implementation details, working with this format can be tricky. Currently, libraries in Python, R and Perl are available, whilst such for JavaScript are lacking. Here, we present a BioJS component for parsing BIOM data in all format versions. It supports import, modification, and export via a unified interface. This module aims to facilitate the development of web applications that use BIOM data. Finally, we demonstrate its usefulness by two applications that already use this component. Availability: https://github.com/molbiodiv/biojs-io-biom, https://dx.doi.org/10.5281/zenodo.218277.

TBro: visualization and management of de novo transcriptomes

RNA sequencing (RNA-seq) has become a powerful tool to understand molecular mechanisms and/or developmental programs. It provides a fast, reliable and cost-effective method to access sets of expressed elements in a qualitative and quantitative manner. Especially for non-model organisms and in absence of a reference genome, RNA-seq data is used to reconstruct and quantify transcriptomes at the same time. Even SNPs, InDels, and alternative splicing events are predicted directly from the data without having a reference genome at hand. A key challenge, especially for non-computational personnal, is the management of the resulting datasets, consisting of different data types and formats. Here, we present TBro, a flexible de novo transcriptome browser, tackling this challenge. TBro aggregates sequences, their annotation, expression levels as well as differential testing results. It provides an easy-to-use interface to mine the aggregated data and generate publication-ready visualizations. Additionally, it supports users with an intuitive cart system, that helps collecting and analysing biological meaningful sets of transcripts. TBro’s modular architecture allows easy extension of its functionalities in the future. Especially, the integration of new data types such as proteomic quantifications or array-based gene expression data is straightforward. Thus, TBro is a fully featured yet flexible transcriptome browser that supports approaching complex biological questions and enhances collaboration of numerous researchers.

Database URL: tbro.carnivorom.com

bcgTree: automatized phylogenetic tree building from bacterial core genomes

The need for multi-gene analyses in scientific fields such as phylogenetics and DNA barcoding has increased in recent years. In particular, these approaches are increasingly important for differentiating bacterial species, where reliance on the standard 16S rDNA marker can result in poor resolution. Additionally, the assembly of bacterial genomes has become a standard task due to advances in next-generation sequencing technologies. We created a bioinformatic pipeline, bcgTree, which uses assembled bacterial genomes either from databases or own sequencing results from the user to reconstruct their phylogenetic history. The pipeline automatically extracts 107 essential single-copy core genes, found in a majority of bacteria, using hidden Markov models and performs a partitioned maximum-likelihood analysis. Here, we describe the workflow of bcgTree and, as a proof-of-concept, its usefulness in resolving the phylogeny of 293 publically available bacterial strains of the genus Lactobacillus. We also evaluate its performance in both low- and high-level taxonomy test sets. The tool is freely available at github ( https://github.com/iimog/bcgTree ) and our institutional homepage ( http://www.dna-analytics.biozentrum.uni-wuerzburg.de ).

Increased efficiency in identifying mixed pollen samples by meta-barcoding with a dual-indexing approach

BACKGROUND

Meta-barcoding of mixed pollen samples constitutes a suitable alternative to conventional pollen identification via light microscopy. Current approaches however have limitations in practicability due to low sample throughput and/or inefficient processing methods, e.g. separate steps for amplification and sample indexing.

RESULTS

We thus developed a new primer-adapter design for high throughput sequencing with the Illumina technology that remedies these issues. It uses a dual-indexing strategy, where sample-specific combinations of forward and reverse identifiers attached to the barcode marker allow high sample throughput with a single sequencing run. It does not require further adapter ligation steps after amplification. We applied this protocol to 384 pollen samples collected by solitary bees and sequenced all samples together on a single Illumina MiSeq v2 flow cell. According to rarefaction curves, 2,000-3,000 high quality reads per sample were sufficient to assess the complete diversity of 95% of the samples. We were able to detect 650 different plant taxa in total, of which 95% were classified at the species level. Together with the laboratory protocol, we also present an update of the reference database used by the classifier software, which increases the total number of covered global plant species included in the database from 37,403 to 72,325 (93% increase).

CONCLUSIONS

This study thus offers improvements for the laboratory and bioinformatical workflow to existing approaches regarding data quantity and quality as well as processing effort and cost-effectiveness. Although only tested for pollen samples, it is furthermore applicable to other research questions requiring plant identification in mixed and challenging samples.