the advent of next - generation sequencing and metagenomics has resulted in increasing numbers of ever - larger datasets describing the community structure and function of a variety of different environments , from the human gut ( arumugam et  al . 2011 ; david et  al . 2014 ) to arctic peat soils ( lipson et  al . 2013 ) and deep - sea vents ( xie et  al .
next - generation sequencing technologies have greatly reduced sequencing costs and speed , and researchers can now affordably study whole microbial communities and functions . prior to this , the focus was on community species composition , studied using 16s rrna targeted amplicon sequencing .
amplicon sequencing does not require the dna coverage that metagenomic studies require and can accurately identify which species are present in a sample ( woese and fox 1977 ; lane et  al .
1985 ; hugenholtz 2002 ) , but it does not provide the depth of information , such as gene function , that full metagenome sequencing and annotation provides .
cost is no longer the primary limiting factor for undertaking metagenomic studies , rather it is now bioinformatics and processing power required to process the data produced .
illumina 's hiseq platform , for example , can affordably sequence the most complex of microbial communities , and the challenge now is to interpret the data produced .
( 2014 ) provide an extensive directory of tools available for different tasks involved in a metagenomic project pipeline , related to a range of  omics studies .
megan ( huson et  al . 2007 ) is a popular graphical user interface program for analyzing and visualizing blast results to study the taxonomy of microbial communities .
while megan typically analyses blast results in a few minutes , running blast searches against reference sequences in a database is computationally intensive and slow for metagenomes ( desai et  al .
2012 ; hunter et  al . 2012 ; thomas , gilbert and meyer 2012 ) .
web - based servers are increasingly popular for processing large amounts of data . with an intuitive web interface and a variety of analytical tools to choose from , mg - rast ( meyer et  al .
mg - rast allows users to upload raw sequence files that are processed through quality filters and annotated using a selection of user - defined parameters , such as reference databases , minimum identity cut - off values , maximum e - values or expect - values , and minimum alignment lengths .
details of the processing procedure can be found in the mg - rast technical report ( wilke et  al .
rapsearch2 ( zhao , tang and ye 2012 ) translates nucleotide sequences and aligns them with annotated protein sequences , reporting to be c. 100-fold faster than blastx with only a 1.3%3.2% reduction in sensitivity ( the proportion of sequences annotated ) . with 
pauda ( huson and xie 2014 ) uses a similar approach and claims to be 10  000-fold faster than blastx , although with a significant reduction in sensitivity .
diamond ( buchfink , xie and huson 2015 ) purports to be both fast and accurate , with a 20  000-fold increase in processing speed compared to blastx . in sensitive mode ,
99% of sequences are aligned , with a speed increase of 2000-fold compared to blastx . like blast with megablast , rapsearch2 and
diamond offer fast and sensitive modes , each coming at the cost of the other .
one codex is a web - based program that uses a different technique to blast and mg - rast to classify sequences ( https://onecodex.com/ ) .
the program designers report that it runs 900 times faster than blast while maintaining similar genus - level sensitivity and precision ( the proportion of annotated sequences that are correctly identified ) , taking hours rather than days to classify most metagenomes .
one codex works by comparing k - mers ( sequences of a set length ) from a sequence to a reference database of k - mers ; the greatest number of 100% k - mer matches determines the classification .
blast and mg - rast classify sequences by matching them with the most similar sequences in a database . unlike mg - rast , one codex does not annotate genes for function .
the choice of database , minimum identity cut - off value ( i.e. sequence match stringency ) , minimum alignment length cut - off value and minimum e - value limit ( the probability a match has occurred by chance ) all influence sequence annotation accuracy , which , in turn , affect the reproducibility and interpretation of the data .
an inherent issue with metagenomic studies is that establishing the accuracy of sequence annotation for environmental samples is practically impossible , given that the quantities of organisms and genes are unknown .
therefore , determining the most effective annotation method is fundamental to investigating environmental communities with confidence .
there are a variety of different reference nucleotide and amino acid databases available for annotating gene or protein sequences ( table s1 , supporting information ) .
the m5nr database ( wilke et  al . 2012 ) incorporates information from a selection of different databases ( see table s1 , supporting information ) , increasing the amount of reference data available for annotation . using a single reference database may be the best option in some cases , for example , where 16s rrna amplicons are used as a method to identify taxa , rather than other genes . whereas taxonomic nomenclature is universal , governed by international conventions , there are multiple approaches for functional classification .
two popular methods include clusters of orthologous groups ( cogs ) ( tatusov , koonin and lipman 1997 ) and the kyoto encyclopedia of genes and genomes ( kegg ) ( kanehisa and goto 2000 ) .
cogs comprise orthologous functions that allow for functional description of poorly characterized genomes based on protein orthologs .
cog descriptions are characterized under cellular processes , information storage and processing , metabolism and poorly characterized .
kegg descriptions are characterized under cellular processes , environmental information processing , genetic information processing , human diseases and metabolism . due to the differences in characterization approaches ,
kegg operates on a subscription basis , and mg - rast uses the latest freely available version ( updated in 2008 ) . selecting a minimum identity cut - off value for metagenome analysis is challenging because interspecific sequence identity varies among genes . too high
a value will accurately identify genes with highly conserved regions , such as 16s rrna or highly conserved coding genes with little synonymous substitution , but may fail to identify genes or non - coding regions that are highly variable .
conversely , a value too low will allow for highly variable genes to be identified , but may also incorrectly identify an organism / function , thus providing false community / function profiles . the optimum identity cut - off point for species identification using the 16s rrna gene
is widely accepted as 97% ( stackebrandt and goebel 1994 ; rossell - mora and amann 2001 ; chun et  al . 2007 ; richter and rossell - mra 2009 ; mende et  al .
2013 ; vtrovsk and baldrian 2013 ) , although this value has its limitations . some species , such as certain rickettsia spp .
, have a 16s rrna gene similarity greater than 97% , thus a cut - off value at this level would not differentiate between the species ( fournier et  al . 2003 ) .
stackebrandt and goebel ( 1994 ) suggest that a higher value may be more appropriate , but fewer sequences would be annotated due to sequencing errors and sequence mutations .
typically , lower cut - off values are suitable for metagenomic studies as the multitude of genes that contain varying degrees of conservation are sequenced .
the default value used by mg - rast , and used in many metagenomic studies ( e.g. tatusov , koonin and lipman 1997 ; lipson et  al . 2013 ) , is 60% , as this allows for identification using less conserved genes and non - coding regions .
a lower value allows shorter sequences to be annotated , although the chance of incorrectly annotating a shorter sequence is higher . a higher value will reduce this chance , but may also reduce the number of annotations overall . combining a low minimum alignment length with a strict minimum identity
cut - off value allows shorter sequences to be annotated but with a high match criteria .
setting maximum e - values and minimum alignment lengths allows stringency of annotations to be controlled .
e - values denote the maximum probability that a sequence annotation has occurred by chance .
lower maximum e - values will reduce the number of possible incorrect annotations , although this also reduces number of annotations retained for analysis .
the aim of this study is to evaluate the accuracy of megan , mg - rast and one codex annotation methods while investigating how using different databases and parameters impact the annotation of metagenomes .
to do this , a novel simulated metagenome was generated using the ncbi whole bacterial genome database and annotated using each pipeline and , for mg - rast , with different reference databases , minimum identity cut - off values , minimum alignment lengths and maximum e - values .
using a simulated metagenome comprising known genome abundances allows the accuracy of annotation to be quantified .
the simulated metagenome was also annotated using megablast , a faster variation of blast , to provide a control and so that megan , mg - rast and one codex could be compared to a standard in sequence annotation . comparing the megan
, mg - rast and one codex annotations to the megablast annotations will quantify the accuracy of these programs for annotating sequences from organisms whose genomes are stored in the ncbi databases .
a simulated metagenome , hereafter simmet , was created using nessm ( jia et  al .
2013 ) , comprising the complete ncbi bacterial genome database ( ftp://ftp.ncbi.nlm.nih.gov/genomes/bacteria/all.fna.tar.gz , may 2013 collection , accessed on 29/04/14 ) .
nessm creates synthetic metagenomes from input genomes based on user - defined parameters ( e.g. sequence count , length and abundance distributions ) that aim to simulate real sequencing data , including expected sequencing errors ( i.e. substitutions , insertions and deletions ) based on the chosen sequencing technology simulated ( see step ii : error models and sequencing coverage bias estimation in jia et  al .
a total of 2400  000 sequences with a read length of 450 base pairs were designated for simulation , based on 454 pyrosequencing .
one strain for each of the 1505 species in the ncbi bacterial genome database was randomly selected to be included in the simulation because certain species , e.g. model organisms and human pathogens such as escherichia coli , salmonella enterica , mycobacterium tuberculosis , bacillus cereus and staphylococcus aureus , have been extensively studied and are overrepresented in the databases .
the species abundance distribution used for simulation was derived from the abundance distribution of a pasture soil metagenome ( sequence count : 2378  586 , mg - rast i d : 4554767.3 ) ( see equation  1 ) . 

( 1)\documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{upgreek } 
 \usepackage{mathrsfs } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document } 
 } { } \begin{equation * } y = - 2490\ln \left ( x \right ) + \ 19748 \end{equation*}\end{document}where x is the randomly selected species rank . the sequences were processed with sickle ( joshi and fass 2011 ) to trim low - quality ends , with the average threshold phred score set at 20 ( a base call error rate of 1% ) .
the simmet metagenome file was annotated with megablast ( available from : http://blast.ncbi.nlm.nih.gov/blast.cgi?page_type = blastdocs&doc_type = download ) as a control , using a reference database of the genomes used to create simmet .
the ncbi nucleotide database ( updated 17/11/2014 ) ( ftp://ftp.ncbi.nlm.nih.gov/blast/db/ ) was also used to assess the annotation performance of megablast .
the maximum e - value selected was 1-e and the minimum alignment length , 15 bases .
megablast annotations using the simmet database will be referred to as  control and those using the ncbi nucleotide database will be referred to as 
megablast. the blast results were uploaded to megan ( version 5.2.3 ) and analyzed using the same parameters used in the blast .
the databases investigated within mg - rast were genbank , greengenes , rdp , refseq , seed , swissprot and trembl .
the m5nr and m5rna databases were excluded from individual sequence analysis , as individual sequence annotations were not available for download from mg - rast for these databases .
for both the megablast and mg - rast annotations , which use a minimum sequence alignment match to annotated sequences , the minimum identity cut - off values tested were as follows : 40% , 50% , 60% , 70% , 80% , 90% , 95% and 97% .
the minimum alignment lengths tested were as follows : 10 , 15 , 20 , 25 , 30 , 25 , 40 , 45 , 50 , 55 and 60 bp .
the maximum e - values tested were as follows : 1-e , 1-e , 1-e and 1-e . aside from testing ,
default parameters were used : 60% , 15 bp and 1-e , respectively , for minimum identify cut - off , minimum alignment length and maximum e - value .
the sequence ids and annotations were extracted from the megablast results ( https://github.com/sandyjmacdonald/blast_parser ) , and full taxonomic lineages were generated for each sequence using the ncbi taxonomy database ( available from : ftp://ftp.ncbi.nlm.nih.gov/pub/taxonomy/ ; ncbi database version generated 13/05/2013 ) .
species level was excluded from analysis due to the high variation in annotated species nomenclature and the accepted caveats associated with microbial species classification ( gevers et  al .
2005 ; achtman and wagner 2008 ) , e.g. horizontal gene transfer ( gogarten and townsend 2005 ; bapteste and boucher 2009 ) . discrepancies identified between databases for organism names were corrected for , such as ncbi using the old name chloroflexia ( as of 17/11/14 ) and mg - rast using the new name chloroflexi for the same class .
unidentified and those that were annotated but were either ambiguously annotated or not annotated at all taxonomic levels had the corresponding levels in the lineage replaced with 
unclassified. for megan and one codex , ncbi taxa ids were used to generate the lineages . the taxonomic lineage for each annotated sequence
was compared to the lineage for the corresponding source sequence in simmet to determine the annotation sensitivity and precision at each taxonomic level .
the effect that minimum identity cut - off values , minimum alignment lengths and maximum e - values had on annotation sensitivity and precision were established using megablast and mg - rast .
the correlations between the relative abundances for each taxon in simmet and in the annotations were calculated using pearson 's product moment correlation coefficient .
the natural logarithms of the relative abundance values were calculated for plotting , as the original distributions would not visually convey the variations in low abundance taxa .
the taxa richness values for each taxonomic level were calculated . unlike investigating taxa , correct functional annotations can not be ascertained with 100% confidence . to investigate functional annotation performance , protein sequences associated with the sequences in simmet
were extracted from genbank records and annotated using the kegg automatic annotation server ( kaas ) ( moriya et  al .
both are web - based functional annotation tools independent of those investigated in this study .
they did not contain sequencing errors and thus provided the best possible indication of the functional annotation accuracy , although the caveats associated with sequence annotation ( e.g. possibly incorrectly assigning a function ) are present .
kegg orthology and cog ids were extracted from the kass and webmga results , respectively , for each sequence annotated and compared with the ids assigned by mg - rast .
the parameters set for the taxonomic investigation were used , and the minimum identity cut - off values investigated were as follows : 40% , 50% , 60% , 70% , 80% , 90% and 95% .
the minimum alignment lengths tested were as follows : 10 , 15 , 20 , 25 , 30 , 25 , 40 , 45 , 50 , 55 and 60 base pairs .
the maximum e - values tested were as follows : 1-e , 1-e , 1-e and 1-e .
nessm produced 2399  077 sequences ( length range : 195459 bp , median length 377 bp ) .
s1 , supporting information ) and 98.7% of sequences are between 300 and 400 base pairs long ( fig . s2 , supporting information ) . of the 2399  077 sequences ,
kass annotated 1  341  362 ( 55.9% ) sequences and webmga annotated 1945  674 sequences ( 81.1% ) .
more stringent parameter values resulted in fewer sequence annotations but had a greater precision ; lower values resulted in more annotations being made , but these comprised increases in both correct and incorrect annotations .
for example , with cut - off values of 95% and 40% , mg - rast refseq annotated 40.2% and 90.3% sequences , respectively , with incorrect annotation rates of 2.9% and 34.5% at the genus level .
this was observed for all parameters tested and for both taxonomic and functional annotations ( figs  16 ) .
as the taxonomic level moved up the taxonomic hierarchy , more sequences were correctly annotated ( e.g. 0.5% and 10.2% for mg - rast refseq with cut - off values of 95% and 40% , respectively , at the class level ) .
the effect of changing minimum identity cut - off value on the number of sequences correctly and incorrectly annotated across the taxonomic levels .
the effect of changing minimum identity cut - off value on the number of sequences correctly and incorrectly annotated for functions .
the effect of changing minimum alignment length on the number of sequences correctly and incorrectly annotated across the taxonomic levels .
the effect of changing minimum alignment length on the number of sequences correctly and incorrectly annotated for functions .
the effect of changing maximum e - value on the number of sequences correctly and incorrectly annotated across the taxonomic levels .
the effect of changing maximum e - value value on the number of sequences correctly and incorrectly annotated for functions .
the correlation coefficients between the taxa relative abundances in simmet and in the annotations decreased as parameter stringency increased ( fig .  7 , associated scatter plots in figs s3s5 , supporting information ) .
most databases achieved maximum correlations with a minimum identity cut - off value of 50% , a minimum alignment length of 30 bp and a maximum e - value of 1-e .
greater decreases in correlation coefficients occurred with a minimum identity cut - off value above 70% and a minimum alignment length greater than 40 bp .
the pearson 's product - moment correlation coefficients for the correlations between the genus relative abundances from simmet and those from various annotation methods using different ( a ) minimum identity cut - off values , ( b ) minimum alignment lengths and ( c ) maximum e - values .
this produced the greatest number of correct annotations ( 99.4% ) ( table  1 , fig .  8 , table s2 , supporting information , for all taxonomic levels ) .
one codex annotated all of the sequences , but incorrectly annotated more sequences ( 5.8% ) than megan ( 2.9% ) , megablast ( 2.5% ) and the control ( 0.5% ) .
megablast , megan and one codex correctly annotated 97.3% , 95.7% and 94.2% sequences respectively , significantly more than the next most successful methods : mg - rast refseq ( 55.9% ) , mg - rast trembl ( 54.9% ) and mg - rast genbank ( 52.7% ) .
mg - rast rdp and mg - rast greengenes , both rrna databases , annotated less than 1% of the sequences .
this is consistent with the expected frequency of rrna genes within bacterial genomes ( vtrovsk and baldrian 2013 ) . as the taxonomic level increases , precision increases and becomes more similar across the different databases .
the annotation sensitivity and number of sequences correctly annotated from a variety of methods and databases across the taxonomic levels investigated .
the simmet taxonomic annotation statistics for each method and database at the genus level using default parameters .
mg - rast kegg annotated 63.3% of the sequences and had a precision of 71.7% , with 45.4% of sequences correctly assigned a function and 17.9% incorrectly assigned a function .
mg - rast cog annotated 50.5% of the sequences and had a precision of 91.1% , resulting in 46.0% of sequence being correctly assigned a function and 4.5% being incorrectly assigned a function ( tables s3s5 , supporting information ) .
the portions of sequences correctly annotated by both methods were 81.5% for mg - rast kegg and 55.4% for mg - rast cog .
megablast had the greatest genus - level correlation with simmet after the control ( r = 0.95 ) , while mg - rast seed had the weakest ( r = 0.49 ) .
megan and one codex had genus - level correlations of r = 0.90 and r = 0.93 , respectively .
the greatest correlation achieved , aside from the control , was by megablast at the phylum level ( r > 0.99 ) ( fig . 
the pearson 's product - moment correlation coefficients for the correlations between the relative abundances from simmet and those from the annotation methods .
mg - rast m5nr and mg - rast refseq generated 87 and 56 false positive class identifications , respectively ( table  2 ) .
megan had only two false positive class identifications (  unidentified and insecta ) and one false negative identification ( solibacteres ) .
one codex also had a low abundance of false positive class identifications ( eight ) and no false negative class identifications
. classes with many false positive identifications include eukaryotes , particularly fungi and bacteria such as spartobacteria .
the greatest fold differences for classes can be found in table s6 ( supporting information ) .
false positive and negative class abundances . the false positive and negative classes from the mg - rast m5nr , refseq , one codex and megan annotations .
six of the annotation methods underestimated the genus richness and six overestimated it ( table  3 , fig . 
the next closest estimate was achieved by megan ( 97.7% ) , followed by mg - rast swissprot ( 95.5% ) , mg - rast m5rna ( 95.2% ) , megablast ( 110.2% ) and mg - rast refseq ( 118.2% ) .
mg - rast m5nr produced the most incorrect richness value at 1244 genera ( 180.8% ) .
the methods were inconsistent in response to the taxonomic level . with increasing taxonomic level , some estimates increased in accuracy while others decreased ( fig . 
10 , table s7 , supporting information ) . excluding the control and the domain level ,
where the number of taxa is low , megan achieved the most accurate richness value ( 101.2% ) at the family level .
mg - rast m5nr achieved the most inaccurate richness value ( 253.3% ) at the order level .
megablast and one codex achieved accurate results relative to other methods , but they still overstated taxa richness at every taxonomic level .
the differences between annotated richness values and the actual richness value ( dashed line ) for each taxonomic level .
the genus richness estimates and the differences from simmet for each annotation method . due to the low numbers
richness values at all taxonomic levels can be found in table s7 ( supporting information ) .
in the study , we evaluated the performances of megan , mg - rast , one codex and megablast by determining their sequence annotation accuracies .
all common taxonomic levels above species are studied , building on the work by lindgreen , adair and gardner ( 2016 ) who study several tools at the genus and phylum levels . by studying a range of taxonomic levels ,
we provide a guideline for researchers to establish the annotation accuracy costs of investigating lower taxonomic levels , allowing them to optimize their investigations depending on their requirements for taxonomic resolution .
mg - rast and megablast use a selection of parameters to determine the stringency of matching a sequence with a reference sequence in a database .
less stringent parameters ( i.e. lower minimum identity cut - off values , lower minimum alignment lengths and higher maximum e - values ) annotate more sequences , but more incorrect annotations are made , thus producing an incorrect community profile .
more stringent parameters reduce the number of incorrect annotations , but many fewer annotations are made , resulting in much of the data being rejected .
decreases in sensitivity generally occur from minimum identify cut - off values above 60% , a minimum alignment lengths greater than 30 bp or a maximum e - value below 1-e ; therefore , the default values used by mg - rast maximize sensitivity . according to carr and borenstein ( 2014 ) ,
the impact of parameters such as e - value will vary depending on read length , something that should be considered in future evaluations as newer sequencing technologies produce longer reads ( e.g. nanopore sequencing ; branton et  al .
the sensitivities and the number of sequences correctly annotated are relatively low for mg - rast at the genus and family levels . at the order level , the values are higher , suggesting that this would be the optimum taxonomic level to study , which maximizes the amount of data used without producing too many incorrect annotations . ultimately , there is a trade - off between taxonomic resolution and annotation accuracy , and this must be considered when determining methods for metagenomic studies .
a marginal number of sequences were not annotated by the control and an even smaller number were incorrectly annotated .
we can therefore conclude that 0.5% of intersample difference at the genus level may be attributed to sequencing error , an important consideration when interpreting data obtained from environmental samples using these methods .
this is supported by hoff ( 2009 ) and carr and borenstein ( 2014 ) , who found that increasing error rates decrease gene prediction accuracy .
as the error rates of next - generation sequencing technologies improve , this effect will reduce .
this is likely to be due to a combination of the kmer - based annotation method that it uses and that the simulated metagenome was created using the ncbi genome database , the primary reference source for one codex .
other than the control , megablast correctly annotated the most sequences at the genus level ( 97.3% ) , although the sensitivity of this method was 0.2% less than one codex .
megan had the second highest precision , annotating 98.6% of sequences , with 95.7% correct annotations .
this suggests that megablast is the most reliable method for annotating sequences , and indicates that it is more conservative than one codex when assigning a sequence hit but also less likely to misidentify a sequence .
megan 's performance was similar to megablast , which is expected as megan processed the megablast output .
mg - rast refseq had the fifth greatest annotation sensitivity and the greatest of the mg - rast annotations ( excluding mg - rast m5nr , for which sequence - specific annotation data were unavailable ) , although it also achieved the greatest number of misidentifications . at the genus level , 33.7% of sequences were misidentified and 55.9% were correctly identified , leaving the remainder unassigned despite the fact that all taxa in simmet are fully sequenced .
this would suggest that investigating metagenomes at the genus level would be unreliable , generating many false positives and implying an incorrect community structure and composition .
this supports garcia - etxebarria , garcia - garcer and calafell ( 2014 ) , who found that more annotations are made at higher taxonomic levels and that discrepancies between known frequencies and annotations increase at lower taxonomic levels , and lindgreen , adair and gardner ( 2016 ) , who report decreases in community annotation accuracy at the genus level compared to phylum . at the class level , the proportion of incorrect annotations is reduced to fewer than 10% for mg - rast refseq , with 80% being annotated correctly .
while taxonomic resolution is reduced , it ensures that the confidence in the annotations remains high .
mg - rast kegg correctly annotated a similar number of sequences to mg - rast cog , but incorrectly annotated many more .
kegg offers a more descriptive annotation as it comprises specific gene and pathway annotations , whereas cog provides descriptions based on orthologous sequences
. however , the specificity of kegg classifications may be the cause of the incorrect annotations as there are more annotations to be selected from and there may be more closely related functions , increasing the chance of misidentification . because kegg is now subscription based , and mg - rast uses the last free version ( 2008 ) , it will not contain information added after that date .
our results are in line with those produced by lindgreen , adair and gardner ( 2016 ) , who also conclude that mg - rast 's functional annotation was accurate . the control , one codex , megablast and megan achieved the greatest correlation coefficients between simmet and annotation abundances at the genus level , all above 0.9 . for all mg - rast annotations ,
the greatest correlation of all abundances was achieved at the order level by the m5nr database , closely followed by trembl and refseq .
these correlations inform us about community - wide analyses , but they are not as sequence sensitivity and precision as correlating abundances values may occur from coincidental incorrect annotations .
mg - rast overannotated many more classes than megan and one codex , for which the most abundant feature was the unidentified group .
this supports the sensitivity and precision data in suggesting that one codex is more likely to categorize unknown sequences as unidentified , rather than incorrectly identifying them .
the genus richness estimated by mg - rast m5nr was 81.0% greater that simmet 's actual richness , the highest overstatement , while megan achieved the most accurate genus richness value ( 2.3% lower ) after the control ( 100.2% ) .
this overstatement could be due to the greater number of sequences present in mg - rast m5nr .
mg - rast m5rna produced a relatively accurate estimate of genus richness ( 95.2% ) ; as m5rna is a 16s rrna database , it is unlikely to annotate non-16s rrna sequences , reducing the number of incorrect identifications .
however , the taxa abundance correlations show that mg - rast m5rna achieved the second lowest correlation with simmet at the genus level , and the lowest at all other taxonomic levels
. mg - rast refseq generated the fifth most accurate richness value , greater than one codex , although not as accurate as megablast and megan .
combined with its high abundance correlation with simmet , this suggests that mg - rast refseq provides a relatively accurate representation of both the richness of a community and the abundance of organisms present .
megan and one codex achieve more accurate taxa richness values and taxa abundance correlations than mg - rast refseq at the family level and above , suggesting they would be a viable alternative to mg - rast refseq .
one limitation with evaluating annotations using organism nomenclature , rather than taxon ids ( which were unavailable for mg - rast sequence - specific annotation data ) , is the lack of taxonomic metadata curation in some databases .
some genomes in the ncbi database are stored with the abbreviated species name rather than complete name , thus amycolatopsis mediterranei would not automatically be identified as an amycolatopsis species .
for example , the class chloroflexia has been renamed to chloroflexi , and is called this by mg - rast . however , ncbi is using the old name chloroflexia ( as of 17/11/14 ) , thus sequences identified as chloroflexi would not be correctly matched in simmet .
these issues were corrected for during data processing ; however , there may be other cases of disparities in the plethora of organisms present in the analysis .
a solution to this would be to use the taxon ids instead ; however , these were not available for sequence - specific annotations downloaded from mg - rast . in conclusion
, we found that one codex , megablast and megan are suitable methods for annotating dna sequences that are located in the reference databases that they use for annotation , with one codex offering fast , web - based analyses and megan providing a user - friendly graphical user interface to analyze blast results .
results appear to vary significantly depending on the program and parameters used , a conclusion also drawn by lindgreen , adair and gardner ( 2016 ) .
while mg - rast appears to have a greater rate of incorrect assignments , this is reduced when investigating higher taxonomic levels ( e.g. with refseq : over 33% at the genus level compared to less than 15% and 10% at the order and class levels ) .
the correlations between the annotated taxa abundances are greatest for mg - rast at the order level , using m5nr , trembl or refseq .
in many of the tests , mg - rast m5nr proved to be a reliable database , but the diversity indices suggest that it is less reliable than mg - rast refseq ; at the class , order and family levels , mg - rast m5nr estimates more the double the actual richness values .
therefore , we hypothesize that mg - rast m5nr would generate more false positive sequence annotations than mg - rast refseq .
( 2007 ) , who evaluated different metagenomic processing methods using simulated metagenome developed from 113 isolated genomes , and by pignatelli and moya ( 2011 ) , who used simulated data to study the performances of de novo short - read assembly programs .
it should be noted that the performances of the methods discussed in this study are likely to differ from the reported results when annotating environmental sequence data ; a greater number of sequences are likely to be unidentified due to the multitude of uncultured microorganisms ( streit and schmitz 2004 ) and non - sequenced microbial genomes ( tringe et  al . 2005 ) that are currently absent from the ncbi whole bacterial genome database .
while this research focused on a selection of annotation methods , the overall conclusions drawn should be considered for any pipeline . in this study , we highlight and quantify the annotation errors for a selection of parameters and databases . we show that analysis pipelines are not equivalent and certain parameters can significantly reduce the confidence in results .
these findings should be used as a guideline when determining methods for annotating metagenomic sequences and considered when interpreting metagenomic results .
ultimately , the most appropriate balance between taxonomic resolution , annotation sensitivity and annotation precision needs to be identified for each study conducted .
 the advent of next - generation sequencing has allowed huge amounts of dna sequence data to be produced , advancing the capabilities of microbial ecosystem studies . 
 the current challenge is to identify from which microorganisms and genes the dna originated . 
 several tools and databases are available for annotating dna sequences . 
 the tools , databases and parameters used can have a significant impact on the results : nave choice of these factors can result in a false representation of community composition and function . 
 we use a simulated metagenome to show how different parameters affect annotation accuracy by evaluating the sequence annotation performances of megan , mg - rast , one codex and megablast . 
 this simulated metagenome allowed the recovery of known organism and function abundances to be quantitatively evaluated , which is not possible for environmental metagenomes . 
 the performance of each program and database varied , e.g. one codex correctly annotated many sequences at the genus level , whereas mg - rast refseq produced many false positive annotations . 
 this effect decreased as the taxonomic level investigated increased . 
 selecting more stringent parameters decreases the annotation sensitivity , but increases precision . ultimately , there is a trade - off between taxonomic resolution and annotation accuracy . 
 these results should be considered when annotating metagenomes and interpreting results from previous studies . 
INTRODUCTION
METHODOLOGY
RESULTS
DISCUSSION
SUPPLEMENTARY DATA
FUNDING