Umeå University's logo

Biodiversity loss threatens ecosystems and human well-being, making accurate, large-scale monitoring crucial. Environmental DNA (eDNA) has enabled species detection from substrates such as water, without the need for direct observation. Lately, airborne eDNA has been showing promise for tracking organisms from insects to mammals in terrestrial ecosystems. Conventional biodiversity assessments are often labor-intensive and limited in scope, leaving gaps in our understanding of ecosystem response to environmental change. Here, we demonstrate that airborne eDNA can detect organisms across the tree of life, quantify changes in abundance congruent with traditional monitoring, and reveal land-use induced regional decline of diversity in a northern boreal ecosystem over more than three decades. By analyzing 34 years of archived aerosol filters, we reconstruct weekly temporal relative abundance data for more than 2700 genera using non-targeted methods. This study provides unified, ecosystem-scale biodiversity surveillance spanning multiple decades, with data collected at weekly intervals on both the individual species and community level. Previously, large scale analyses of ecosystem changes, targeting all types of organisms, has been prohibitively expensive and difficult to attempt. Here, we present a way of holistically doing this type of analysis in a single framework.

Crop pests and diseases increasingly challenge the global food system. To prepare for and detect outbreaks, surveillance plays an important role. Traditional monitoring methods are often organism-specific, making large-scale monitoring of crop pathogens and pests impractical. We here investigate the potential for using shotgun sequencing of airborne eDNA for large-scale surveillance of crop pathogens and pests. We show that it is possible to detect DNA from all types of organisms in air, and that DNA can be classified down to species level. However, the accuracy of the identification is highly dependent on the quality of reference genomes of both the pathogens or pests, and their close relatives present in the region. Finally, we find that observed degree of crop damages correlate with amount of DNA from crop pathogens and pests in air, showing the promise of this approach for surveillance of all types of crop pathogens and pests.

The Flexible Taxonomy Database framework provides a method for modification and merging official and custom taxonomic databases to create improved databases. Using such databases will increase accuracy and precision of existing methods to classify sequence reads.

The number of national reference populations that are whole-genome sequenced are rapidly increasing. Partly driving this development is the fact that genetic disease studies benefit from knowing the genetic variation typical for the geographical area of interest. A whole-genome sequenced Swedish national reference population (n = 1000) has been recently published but with few samples from northern Sweden. In the present study we have whole-genome sequenced a control population (n = 300) (ACpop) from Västerbotten County, a sparsely populated region in northern Sweden previously shown to be genetically different from southern Sweden. The aggregated variant frequencies within ACpop are publicly available (DOI 10.17044/NBIS/G000005) to function as a basic resource in clinical genetics and for genetic studies. Our analysis of ACpop, representing approximately 0.11% of the population in Västerbotten, indicates the presence of a genetic substructure within the county. Furthermore, a demographic analysis showed that the population from which samples were drawn was to a large extent geographically stationary, a finding that was corroborated in the genetic analysis down to the level of municipalities. Including ACpop in the reference population when imputing unknown variants in a Västerbotten cohort resulted in a strong increase in the number of high-confidence imputed variants (up to 81% for variants with minor allele frequency < 5%). ACpop was initially designed for cancer disease studies, but the genetic structure within the cohort will be of general interest for all genetic disease studies in northern Sweden.

Whole-genome sequencing is a promising approach for human autosomal dominant disease studies. However, the vast number of genetic variants observed by this method constitutes a challenge when trying to identify the causal variants. This is often handled by restricting disease studies to the most damaging variants, e.g. those found in coding regions, and overlooking the remaining genetic variation. Such a biased approach explains in part why the genetic causes of many families with dominantly inherited diseases, in spite of being included in whole-genome sequencing studies, are left unsolved today. Here we explore the use of a geographically matched control population to minimize the number of candidate disease-causing variants without excluding variants based on assumptions on genomic position or functional predictions. To exemplify the benefit of the geographically matched control population we apply a typical disease variant filtering strategy in a family with an autosomal dominant form of colorectal cancer. With the use of the geographically matched control population we end up with 26 candidate variants genome wide. This is in contrast to the tens of thousands of candidates left when only making use of available public variant datasets. The effect of the local control population is dual, it (1) reduces the total number of candidate variants shared between affected individuals, and more importantly (2) increases the rate by which the number of candidate variants are reduced as additional affected family members are included in the filtering strategy. We demonstrate that the application of a geographically matched control population effectively limits the number of candidate disease-causing variants and may provide the means by which variants suitable for functional studies are identified genome wide.

Background: Selecting the proper parameter settings for bioinformatic software tools is challenging. Not only will each parameter have an individual effect on the outcome, but there are also potential interaction effects between parameters. Both of these effects may be difficult to predict. To make the situation even more complex, multiple tools may be run in a sequential pipeline where the final output depends on the parameter configuration for each tool in the pipeline. Because of the complexity and difficulty of predicting outcomes, in practice parameters are often left at default settings or set based on personal or peer experience obtained in a trial and error fashion. To allow for the reliable and efficient selection of parameters for bioinformatic pipelines, a systematic approach is needed.

Results: We present doepipeline, a novel approach to optimizing bioinformatic software parameters, based on core concepts of the Design of Experiments methodology and recent advances in subset designs. Optimal parameter settings are first approximated in a screening phase using a subset design that efficiently spans the entire search space, then optimized in the subsequent phase using response surface designs and OLS modeling. Doepipeline was used to optimize parameters in four use cases; 1) de-novo assembly, 2) scaffolding of a fragmented genome assembly, 3) k-mer taxonomic classification of Oxford Nanopore Technologies MinION reads, and 4) genetic variant calling. In all four cases, doepipeline found parameter settings that produced a better outcome with respect to the characteristic measured when compared to using default values. Our approach is implemented and available in the Python package doepipeline.

Conclusions: Our proposed methodology provides a systematic and robust framework for optimizing software parameter settings, in contrast to labor- and time-intensive manual parameter tweaking. Implementation in doepipeline makes our methodology accessible and user-friendly, and allows for automatic optimization of tools in a wide range of cases. The source code of doepipeline is available at https://github.com/clicumu/doepipeline and it can be installed through conda-forge.

We present the complete genome sequence of Francisella guangzhouensis strain 08HL01032(T), which consists of one chromosome (1,658,482 bp) and one plasmid ( 3,045 bp) with G+C contents of 32.0% and 28.7%, respectively.