Strategies to improve usability and preserve accuracy in biological sequence databases

Next-generation sequencing and PCR technologies in monitoring the hospital microbiome and its drug resistance

The hospital environment significantly contributes to the onset of healthcare-associated infections (HAIs), which represent one of the most frequent complications occurring in healthcare facilities worldwide. Moreover, the increased antimicrobial resistance (AMR) characterizing HAI-associated microbes is one of the human health’s main concerns, requiring the characterization of the contaminating microbial population in the hospital environment. The monitoring of surface microbiota in hospitals is generally addressed by microbial cultural isolation. However, this has some important limitations mainly relating to the inability to define the whole drug-resistance profile of the contaminating microbiota and to the long time period required to obtain the results. Hence, there is an urgent need to implement environmental surveillance systems using more effective methods. Molecular approaches, including next-generation sequencing and PCR assays, may be useful and effective tools to monitor microbial contamination, especially the growing AMR of HAI-associated pathogens. Herein, we summarize the results of our recent studies using culture-based and molecular analyses in 12 hospitals for adults and children over a 5-year period, highlighting the advantages and disadvantages of the techniques used.

Metabolic networks of the Nicotiana genus in the spotlight: content, progress and outlook

Manually curated metabolic databases residing at the Sol Genomics Network comprise two taxon-specific databases for the Solanaceae family, i.e. SolanaCyc and the genus Nicotiana, i.e. NicotianaCyc as well as six species-specific databases for Nicotiana tabacum TN90, N. tabacum K326, Nicotiana benthamiana, N. sylvestris, N. tomentosiformis and N. attenuata. New pathways were created through the extraction, examination and verification of related data from the literature and the aid of external database guided by an expert-led curation process. Here we describe the curation progress that has been achieved in these databases since the first release version 1.0 in 2016, the curation flow and the curation process using the example metabolic pathway for cholesterol in plants. The current content of our databases comprises 266 pathways and 36 superpathways in SolanaCyc and 143 pathways plus 21 superpathways in NicotianaCyc, manually curated and validated specifically for the Solanaceae family and Nicotiana genus, respectively. The curated data have been propagated to the respective Nicotiana-specific databases, which resulted in the enrichment and more accurate presentation of their metabolic networks. The quality and coverage in those databases have been compared with related external databases and discussed in terms of literature support and metabolic content.

Developing a non-destructive metabarcoding protocol for detection of pest insects in bulk trap catches

Metabarcoding has the potential to revolutionise insect surveillance by providing high-throughput and cost-effective species identification of all specimens within mixed trap catches. Nevertheless, incorporation of metabarcoding into insect diagnostic laboratories will first require the development and evaluation of protocols that adhere to the specialised regulatory requirements of invasive species surveillance. In this study, we develop a multi-locus non-destructive metabarcoding protocol that allows sensitive detection of agricultural pests, and subsequent confirmation using traditional diagnostic techniques. We validate this protocol for the detection of tomato potato psyllid (Bactericera cockerelli) and Russian wheat aphid (Diuraphis noxia) within mock communities and field survey traps. We find that metabarcoding can reliably detect target insects within mixed community samples, including specimens that morphological identification did not initially detect, but sensitivity appears inversely related to community size and is impacted by primer biases, target loci, and sample indexing strategy. While our multi-locus approach allowed independent validation of target detection, lack of reference sequences for 18S and 12S restricted its usefulness for estimating diversity in field samples. The non-destructive DNA extraction proved invaluable for resolving inconsistencies between morphological and metabarcoding identification results, and post-extraction specimens were suitable for both morphological re-examination and DNA re-extraction for confirmatory barcoding.

On the complementarity of DNA barcoding and morphology to distinguish benign endemic insects from possible pests: the case of Dirioxa pornia and the tribe Acanthonevrini (Diptera: Tephritidae: Phytalmiinae) in Australia

Fruit flies are considered economically important insects due to some species being agricultural pests. However, morphological identification of fruit fly adults and larvae can be difficult requiring a high level of taxonomic expertise, with misidentifications causing problematic false-positive/negative results. While destructive molecular techniques can assist with the identification process, these often cannot be applied where it is mandatory to retain a voucher reference specimen. In this work, we non-destructively (and partial-destructively) processed larvae and adults mostly belonging to the species Dirioxa pornia (Walker, 1849), of the poorly studied nonpest fruit fly tribe Acanthonevrini (Tephritidae) from Australia, to enable molecular identifications whilst retaining morphological vouchers. By retaining the morphological features of specimens, we confirmed useful characters for genus/species-level identification, contributing to improved accuracy for future diagnostics using both molecular and morphological approaches. We provide DNA barcode information for three species of Acanthonevrini known from Australia, which prior to our study was only available for a single species, D. pornia. Our specimen examinations provide new distribution records for three nonpest species: Acanthonevroides variegatus Permkam and Hancock, 1995 in South Australia, Acanthonevroides basalis (Walker, 1853) and D. pornia in Victoria, Australia; as well as new host plant records for D. pornia, from kangaroo apple, apricot and loquat.

ARGDIT: a validation and integration toolkit for Antimicrobial Resistance Gene Databases

MOTIVATION

Antimicrobial resistance is currently one of the main challenges in public health due to the excessive use of antimicrobials in medical treatments and agriculture. The advancements in high-throughput next-generation sequencing and development of bioinformatics tools allow simultaneous detection and identification of antimicrobial resistance genes (ARGs) from clinical, food and environment samples, to monitor the prevalence and track the dissemination of these ARGs. Such analyses are however reliant on a comprehensive database of ARGs with accurate sequence content and annotation. Most of the current ARG databases are therefore manually curated, but this is a time-consuming process and the resulting curation errors could be hard to detect. Several secondary ARG databases consolidate contents from different source ARG databases, and hence modifications in the primary databases might not be propagated and updated promptly in the secondary ARG databases.

RESULTS

To address these problems, a validation and integration toolkit called ARGDIT was developed to validate ARG database fidelity, and merge multiple primary ARG databases into a single consolidated secondary ARG database with optional automated sequence re-annotation. Experimental results demonstrated the effectiveness of this toolkit in identifying errors such as sequence annotation typos in current ARG databases and generating an integrated non-redundant ARG database with structured annotation. A toolkit-oriented workflow is also proposed to minimize the efforts in validating, curating and merging multiple ARG protein or coding sequence databases. Database developers therefore benefit from faster update cycles and lower costs for database maintenance, while ARG pipeline users can easily evaluate the reference ARG database quality.

AVAILABILITY AND IMPLEMENTATION

ARGDIT is available at https://github.com/phglab/ARGDIT.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Pre- and post-sequencing recommendations for functional annotation of human fecal metagenomes

Background

Shotgun metagenomes are often assembled prior to annotation of genes which biases the functional capacity of a community towards its most abundant members. For an unbiased assessment of community function, short reads need to be mapped directly to a gene or protein database. The ability to detect genes in short read sequences is dependent on pre- and post-sequencing decisions. The objective of the current study was to determine how library size selection, read length and format, protein database, e-value threshold, and sequencing depth impact gene-centric analysis of human fecal microbiomes when using DIAMOND, an alignment tool that is up to 20,000 times faster than BLASTX.

Results

Using metagenomes simulated from a database of experimentally verified protein sequences, we find that read length, e-value threshold, and the choice of protein database dramatically impact detection of a known target, with best performance achieved with longer reads, stricter e-value thresholds, and a custom database. Using publicly available metagenomes, we evaluated library size selection, paired end read strategy, and sequencing depth. Longer read lengths were acheivable by merging paired ends when the sequencing library was size-selected to enable overlaps. When paired ends could not be merged, a congruent strategy in which both ends are independently mapped was acceptable. Sequencing depths of 5 million merged reads minimized the error of abundance estimates of specific target genes, including an antimicrobial resistance gene.

Conclusions

Shotgun metagenomes of DNA extracted from human fecal samples sequenced using the Illumina platform should be size-selected to enable merging of paired end reads and should be sequenced in the PE150 format with a minimum sequencing depth of 5 million merge-able reads to enable detection of specific target genes. Expecting the merged reads to be 180-250 bp in length, the appropriate e-value threshold for DIAMOND would then need to be more strict than the default. Accurate and interpretable results for specific hypotheses will be best obtained using small databases customized for the research question.

Pulse Crop Genetics for a Sustainable Future: Where We Are Now and Where We Should Be Heading

The last decade has witnessed dramatic changes in global food consumption patterns mainly because of population growth and economic development. Food substitutions for healthier eating, such as swapping regular servings of meat for protein-rich crops, is an emerging diet trend that may shape the future of food systems and the environment worldwide. To meet the erratic consumer demand in a rapidly changing world where resources become increasingly scarce due largely to anthropogenic activity, the need to develop crops that benefit both human health and the environment has become urgent. Legumes are often considered to be affordable plant-based sources of dietary proteins. Growing legumes provides significant benefits to cropping systems and the environment because of their natural ability to perform symbiotic nitrogen fixation, which enhances both soil fertility and water-use efficiency. In recent years, the focus in legume research has seen a transition from merely improving economically important species such as soybeans to increasingly turning attention to some promising underutilized species whose genetic resources hold the potential to address global challenges such as food security and climate change. Pulse crops have gained in popularity as an affordable source of food or feed; in fact, the United Nations designated 2016 as the International Year of Pulses, proclaiming their critical role in enhancing global food security. Given that many studies have been conducted on numerous underutilized pulse crops across the world, we provide a systematic review of the related literature to identify gaps and opportunities in pulse crop genetics research. We then discuss plausible strategies for developing and using pulse crops to strengthen food and nutrition security in the face of climate and anthropogenic changes.

Metaxa2 Database Builder: enabling taxonomic identification from metagenomic or metabarcoding data using any genetic marker

Motivation

Correct taxonomic identification of DNA sequences is central to studies of biodiversity using both shotgun metagenomic and metabarcoding approaches. However, no genetic marker gives sufficient performance across all the biological kingdoms, hampering studies of taxonomic diversity in many groups of organisms. This has led to the adoption of a range of genetic markers for DNA metabarcoding. While many taxonomic classification software tools can be re-trained on these genetic markers, they are often designed with assumptions that impair their utility on genes other than the SSU and LSU rRNA. Here, we present an update to Metaxa2 that enables the use of any genetic marker for taxonomic classification of metagenome and amplicon sequence data.

Results

We evaluated the Metaxa2 Database Builder on 11 commonly used barcoding regions and found that while there are wide differences in performance between different genetic markers, our software performs satisfactorily provided that the input taxonomy and sequence data are of high quality.

Availability and implementation

Freely available on the web as part of the Metaxa2 package at http://microbiology.se/software/metaxa2/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Prospects and challenges of implementing DNA metabarcoding for high-throughput insect surveillance

Trap-based surveillance strategies are widely used for monitoring of invasive insect species, aiming to detect newly arrived exotic taxa as well as track the population levels of established or endemic pests. Where these surveillance traps have low specificity and capture non-target endemic species in excess of the target pests, the need for extensive specimen sorting and identification creates a major diagnostic bottleneck. While the recent development of standardized molecular diagnostics has partly alleviated this requirement, the single specimen per reaction nature of these methods does not readily scale to the sheer number of insects trapped in surveillance programmes. Consequently, target lists are often restricted to a few high-priority pests, allowing unanticipated species to avoid detection and potentially establish populations. DNA metabarcoding has recently emerged as a method for conducting simultaneous, multi-species identification of complex mixed communities and may lend itself ideally to rapid diagnostics of bulk insect trap samples. Moreover, the high-throughput nature of recent sequencing platforms could enable the multiplexing of hundreds of diverse trap samples on a single flow cell, thereby providing the means to dramatically scale up insect surveillance in terms of both the quantity of traps that can be processed concurrently and number of pest species that can be targeted. In this review of the metabarcoding literature, we explore how DNA metabarcoding could be tailored to the detection of invasive insects in a surveillance context and highlight the unique technical and regulatory challenges that must be considered when implementing high-throughput sequencing technologies into sensitive diagnostic applications.

Best practice data life cycle approaches for the life sciences

Throughout history, the life sciences have been revolutionised by technological advances; in our era this is manifested by advances in instrumentation for data generation, and consequently researchers now routinely handle large amounts of heterogeneous data in digital formats. The simultaneous transitions towards biology as a data science and towards a 'life cycle' view of research data pose new challenges. Researchers face a bewildering landscape of data management requirements, recommendations and regulations, without necessarily being able to access data management training or possessing a clear understanding of practical approaches that can assist in data management in their particular research domain. Here we provide an overview of best practice data life cycle approaches for researchers in the life sciences/bioinformatics space with a particular focus on 'omics' datasets and computer-based data processing and analysis. We discuss the different stages of the data life cycle and provide practical suggestions for useful tools and resources to improve data management practices.

On expert curation and scalability: UniProtKB/Swiss-Prot as a case study

Motivation

Biological knowledgebases, such as UniProtKB/Swiss-Prot, constitute an essential component of daily scientific research by offering distilled, summarized and computable knowledge extracted from the literature by expert curators. While knowledgebases play an increasingly important role in the scientific community, their ability to keep up with the growth of biomedical literature is under scrutiny. Using UniProtKB/Swiss-Prot as a case study, we address this concern via multiple literature triage approaches.

Results

With the assistance of the PubTator text-mining tool, we tagged more than 10 000 articles to assess the ratio of papers relevant for curation. We first show that curators read and evaluate many more papers than they curate, and that measuring the number of curated publications is insufficient to provide a complete picture as demonstrated by the fact that 8000–10 000 papers are curated in UniProt each year while curators evaluate 50 000–70 000 papers per year. We show that 90% of the papers in PubMed are out of the scope of UniProt, that a maximum of 2–3% of the papers indexed in PubMed each year are relevant for UniProt curation, and that, despite appearances, expert curation in UniProt is scalable.

Availability and implementation

UniProt is freely available at http://www.uniprot.org/.

Supplementary information

Supplementary data are available at Bioinformatics online.

The diversity of uncharacterized antibiotic resistance genes can be predicted from known gene variants-but not always

BACKGROUND

Antibiotic resistance is considered one of the most urgent threats to modern healthcare, and the role of the environment in resistance development is increasingly recognized. It is often assumed that the abundance and diversity of known resistance genes are representative also for the non-characterized fraction of the resistome in a given environment, but this assumption has not been verified. In this study, this hypothesis is tested, using the resistance gene profiles of 1109 metagenomes from various environments.

RESULTS

This study shows that the diversity and abundance of known antibiotic resistance genes can generally predict the diversity and abundance of undescribed resistance genes. However, the extent of this predictability is dependent on the type of environment investigated. Furthermore, it is shown that carefully selected small sets of resistance genes can describe total resistance gene diversity remarkably well.

CONCLUSIONS

The results of this study suggest that knowledge gained from large-scale quantifications of known resistance genes can be utilized as a proxy for unknown resistance factors. This is important for current and proposed monitoring efforts for environmental antibiotic resistance and has implications for the design of risk ranking strategies and the choices of measures and methods for describing resistance gene abundance and diversity in the environment.

Using metagenomics to investigate human and environmental resistomes

Antibiotic resistance is a global health concern declared by the WHO as one of the largest threats to modern healthcare. In recent years, metagenomic DNA sequencing has started to be applied as a tool to study antibiotic resistance in different environments, including the human microbiota. However, a multitude of methods exist for metagenomic data analysis, and not all methods are suitable for the investigation of resistance genes, particularly if the desired outcome is an assessment of risks to human health. In this review, we outline the current state of methods for sequence handling, mapping to databases of resistance genes, statistical analysis and metagenomic assembly. In addition, we provide an overview of important considerations related to the analysis of resistance genes, and recommend some of the currently used tools and methods that are best equipped to inform research and clinical practice related to antibiotic resistance.

Best practice data life cycle approaches for the life sciences

Throughout history, the life sciences have been revolutionised by technological advances; in our era this is manifested by advances in instrumentation for data generation, and consequently researchers now routinely handle large amounts of heterogeneous data in digital formats. The simultaneous transitions towards biology as a data science and towards a 'life cycle' view of research data pose new challenges. Researchers face a bewildering landscape of data management requirements, recommendations and regulations, without necessarily being able to access data management training or possessing a clear understanding of practical approaches that can assist in data management in their particular research domain. Here we provide an overview of best practice data life cycle approaches for researchers in the life sciences/bioinformatics space with a particular focus on 'omics' datasets and computer-based data processing and analysis. We discuss the different stages of the data life cycle and provide practical suggestions for useful tools and resources to improve data management practices.