Aquarium Nitrification Revisited: Thaumarchaeota Are the Dominant Ammonia Oxidizers in Freshwater Aquarium Biofilters

Ammonia-oxidizing archaea (AOA) outnumber ammonia-oxidizing bacteria (AOB) in many terrestrial and aquatic environments. Although nitrification is the primary function of aquarium biofilters, very few studies have investigated the microorganisms responsible for this process in aquaria. This study used quantitative real-time PCR (qPCR) to quantify the ammonia monooxygenase (amoA) and 16S rRNA genes of Bacteria and Thaumarchaeota in freshwater aquarium biofilters, in addition to assessing the diversity of AOA amoA genes by denaturing gradient gel electrophoresis (DGGE) and clone libraries. AOA were numerically dominant in 23 of 27 freshwater biofilters, and in 12 of these biofilters AOA contributed all detectable amoA genes. Eight saltwater aquaria and two commercial aquarium nitrifier supplements were included for comparison. Both thaumarchaeal and bacterial amoA genes were detected in all saltwater samples, with AOA genes outnumbering AOB genes in five of eight biofilters. Bacterial amoA genes were abundant in both supplements, but thaumarchaeal amoA and 16S rRNA genes could not be detected. For freshwater aquaria, the proportion of amoA genes from AOA relative to AOB was inversely correlated with ammonium concentration. DGGE of AOA amoA genes revealed variable diversity across samples, with nonmetric multidimensional scaling (NMDS) indicating separation of freshwater and saltwater fingerprints. Composite clone libraries of AOA amoA genes revealed distinct freshwater and saltwater clusters, as well as mixed clusters containing both freshwater and saltwater amoA gene sequences. These results reveal insight into commonplace residential biofilters and suggest that aquarium biofilters may represent valuable biofilm microcosms for future studies of AOA ecology

PANDAseq: paired-end assembler for illumina sequences

Abstract Background Illumina paired-end reads are used to analyse microbial communities by targeting amplicons of the 16S rRNA gene. Publicly available tools are needed to assemble overlapping paired-end reads while correcting mismatches and uncalled bases; many errors could be corrected to obtain higher sequence yields using quality information. Results PANDAseq assembles paired-end reads rapidly and with the correction of most errors. Uncertain error corrections come from reads with many low-quality bases identified by upstream processing. Benchmarks were done using real error masks on simulated data, a pure source template, and a pooled template of genomic DNA from known organisms. PANDAseq assembled reads more rapidly and with reduced error incorporation compared to alternative methods. Conclusions PANDAseq rapidly assembles sequences and scales to billions of paired-end reads. Assembly of control libraries showed a 4-50% increase in the number of assembled sequences over naïve assembly with negligible loss of "good" sequence.</p

BAMQL: a query language for extracting reads from BAM files

Abstract Background It is extremely common to need to select a subset of reads from a BAM file based on their specific properties. Typically, a user unpacks the BAM file to a text stream using SAMtools, parses and filters the lines using AWK, then repacks them using SAMtools. This process is tedious and error-prone. In particular, when working with many columns of data, mix-ups are common and the bit field containing the flags is unintuitive. There are several libraries for reading BAM files, such as Bio-SamTools for Perl and pysam for Python. Both allow access to the BAM’s read information and can filter reads, but require substantial boilerplate code; this is high overhead for mostly ad hoc filtering. Results We have created a query language that gathers reads using a collection of predicates and common logical connectives. Queries run faster than equivalents and can be compiled to native code for embedding in larger programs. Conclusions BAMQL provides a user-friendly, powerful and performant way to extract subsets of BAM files for ad hoc analyses or integration into applications. The query language provides a collection of predicates beyond those in SAMtools, and more flexible connectives

BAMQL:A query language for extracting reads from BAM files

BACKGROUND: It is extremely common to need to select a subset of reads from a BAM file based on their specific properties. Typically, a user unpacks the BAM file to a text stream using SAMtools, parses and filters the lines using AWK, then repacks them using SAMtools. This process is tedious and error-prone. In particular, when working with many columns of data, mix-ups are common and the bit field containing the flags is unintuitive. There are several libraries for reading BAM files, such as Bio-SamTools for Perl and pysam for Python. Both allow access to the BAM’s read information and can filter reads, but require substantial boilerplate code; this is high overhead for mostly ad hoc filtering. RESULTS: We have created a query language that gathers reads using a collection of predicates and common logical connectives. Queries run faster than equivalents and can be compiled to native code for embedding in larger programs. CONCLUSIONS: BAMQL provides a user-friendly, powerful and performant way to extract subsets of BAM files for ad hoc analyses or integration into applications. The query language provides a collection of predicates beyond those in SAMtools, and more flexible connectives. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s12859-016-1162-y) contains supplementary material, which is available to authorized users