BACKGROUND
the conventional process of removing nitrogen from municipal wastewater streams involves nitrification followed by denitrification.
nitrification in the wastewater treatment process is generally performed in two consecutive steps.
the first step is ammonia oxidation by either ammonia-oxidizing bacteria  or ammonia-oxidizing archaea , and the second step is nitrite oxidation by nitrite-oxidizing bacteria .
representative aob genera include nitrosomonas and nitrosospira while nitrososphaera and nitrosopumilus would account as two examples of aoa.
representative nob include nitrospira and nitrobacter.
in recent years, three members within the genus nitrospira, namely candidatus nitrospira inopinata, candidatus nitrospira nitrosa and candidatus nitrospira nitrificans, were shown to perform complete nitrification in a single step and are referred to as comammox  <cit> .
although comammox are not yet isolated from municipal wastewater treatment plants  and their contribution in terms of full nitrification in such systems remain unknown, comammox had been shown to use urea as ammonium source for nitrification and that this trait should in theory enable them to thrive in wwtp environments where urea is often present  <cit> .
the commensal interaction between aob/aoa and nob, and the potential importance of nitrospira-affiliated comammox suggest that the nitrification process can be prone to failure when there is a suboptimal ratio of aob/aoa to nob, particularly that of nitrospira.
this is primarily due to the effect of hydroxylamine, an intermediate of the conversion of ammonia to nitrite, which has a strong inhibitory effect on nitrite oxidizers when released by aob/aoa  <cit> .
failure of the nitrification system in a wwtp can result in discharge of effluent that contains total nitrogen  higher than the permissible limit .
this can result in potentially detrimental environmental impacts in the form of eutrophication of the receiving water bodies.
hence, establishing an optimal ratio and abundance of aob/aoa to nob is important in ensuring the functionality of the nitrification process so that treated effluent will meet the discharge limits for tn.
in such instance, a quantitative method capable of efficiently targeting aob/aoa and nob relative abundances would be particularly useful in tracking the proportional ratio between these two groups.
various methods like quantitative pcr  and fluorescent in-situ hybridization  were used to directly quantify the abundances of aob and nob.
through these approaches, it was determined that nitrosomonas likely accounts for the largest fraction of aob in conventional activated sludge systems.
additionally, in conventional municipal wastewater streams in which ammonium concentrations ranged from  <dig> to 50 mg/l, nitrospira was the predominant nob while nitrobacter was generally in low abundance .
although qpcr and fish have been routinely used in earlier studies, both approaches have substantial limitations.
for example, qpcr requires standard curves to be established for individual targets.
new generation of pcr  does not require standard curves to be generated but similar to conventional qpcr, remains limited by the number of fluorophores available in the ultraviolet spectrum.
hence, both quantitative pcr approaches have limited multiplexing throughput capability.
furthermore, both qpcr approaches can be prone to non-specific amplification when sybr® green fluorescence reporters are utilized, and would require an additional melting curve analysis or verification step to differentiate between true-positive and false-positive amplifications.
although taqman® probes can be developed for use in qpcr and dpcr, such probes are relatively expensive and remain challenging to design for a specific microbial target within the internal region of two flanking qpcr primers.
digital pcr being a considerably new technology, also remains rather costly till this date compared to traditional qpcr.
fish, on the other hand, can be affected by the in-situ accessibility of labeled probes through the cell membranes and the copy numbers of rrna genes inside aob and nob.
for example, species of nitrobacter, nitrosomonas, nitrosospira and nitrospira are reported to have one copy of 16s rrna gene per cell.
this copy number of 16s rrna genes per genome is lower than the mean copy number of  <dig>  among the  <dig> microbial genomes currently available in the same ribosomal rna operon copy number database  <cit> .
hierarchical oligonucleotide primer extension  was developed as a method to complement other quantitative methods like qpcr and fish by providing a high-throughput multiplexing platform capable of addressing some of the previous methods’ limitations.
in this method, oligonucleotide primers are designed to target genes, usually 16s rrna genes, at different hierarchical taxonomical levels .
in the presence of a dna template, the primers anneal at the designed targeting position of the 16s rrna gene and extend with a single fluorescently labeled dideoxynucleoside triphosphate .
the extended primers can then be identified on a dna sequencer based on their anticipated fragment size and fluorescence.
the peak areas of the individual fragments can then be used to determine the relative abundances of bacterial targets of interest  <cit> .
the hope method has been developed to detect specific bacterial populations of interest including host-associated bacteroidales and bifidobacterium, cyanobacteria, and methanogens .
in these studies, hope has been demonstrated to multiplex up to  <dig> primers in a single reaction, allowing for the evaluation of a range of relevant bacterial indicators in a time-efficient manner  <cit> .
the high-throughput capability of hope stems from the use of different fragment sizes of primers and fluorescence colors to distinguish between extended primers.
to exemplify, most primers are synthesized with a length of  <dig> to  <dig> nucleotides  while most dna sequencer platforms are able to detect up to four different fluorescence colors.
hence, in theory, the combination of these two differentiating factors would mean that hope is capable of multiplex up to  <dig> primers in one single reaction tube.
aside from its multiplexing capability, it was demonstrated that hope has a detection limit of  <dig> % of the total pcr-amplified bacterial targets, and that the entire duration required to carry out a single hope reaction and to analyze for the extended primer on a dna sequencer is <1 h, which would be comparable to the time needed to carry out qpcr  <cit> .
given the advantages of hope, the method is a potentially useful tool in tracking the relative abundance of nitrifying populations.
in the present study, we applied the hope approach to target groups of aob and nob  in wastewater environments.
the method was tested against samples collected from a lab-scale biofilm-based trickling bioreactor.
the total microbial community, including the nitrifying groups, was analyzed by 16s rrna gene-based high-throughput sequencing for comparison against the hope data.
to further demonstrate the use of hope in monitoring these relative abundances, the trickling bioreactor was subsequently challenged with high concentrations of incoming ammonium content so as to simulate a toxic shock event that could potentially crash the nitrification system.
the changes in aob and nob proportions were then tracked by the hope approach and correlated to reactor performance.
wastewater samples from a fully operational wastewater treatment plant  were also analyzed in order to further validate the method.
methods
hope primer layout
hope was performed with a total of  <dig> primers targeting ammonia-oxidizing genera nitrosomonas and nitrosospira as well as nitrite-oxidizing genera nitrospira and nitrobacter .
primers were sourced from the probebase database  <cit> .
the mismatch positions of non-targets were moved to the 3′-end of the hope primers as it was previously found that mismatches located at the 3′-end of the hope primers would not facilitate nucleotide base extension of primers, and would hence lead to improved specificity  <cit> .
re-positioning of the mismatches result in minimal change to the original coverage of the targeted bacterial groups.
the coverage of the primers was then verified in-silico against the ribosomal database project  database  <cit> , and shown in additional file  <dig>  primers were arranged into four tubes, with either 27f or 338f primer included in each tube to provide a normalization of the abundance of each nitrifying bacterial target against the total bacteria.
to differentiate primers that extend with the same ddntp and hence color during capillary electrophoresis, some primers were modified with a polya tail at the 5′-end of the primer .
primers in each tube were prepared to form a final concentration of 10 μm.table  <dig> hope primers used to target the different ammonia- and nitrite-oxidizing bacterial groups
nob tube 1: nitrospira
nob tube 2: nitrobacter
aob tube 1: nitrosomonas and nitrosospira
aob tube 2: nitrosomonas and nitrosospira
target coverage for each primer was identified using the rdp probe match function
hope reference template preparation
to prepare the hope reference templates for nitrifying bacterial populations, extracted dna from a separate trickling nitrification biofilter that was set up in year  <dig> was first amplified with 27f  to 1492r  primers and purified for use in 16s rrna gene clone library construction.
approximately  <dig> clones were picked and sequenced to identify nitrifying bacterial species.
two clones associated with candidatus nitrospira defluvii and nitrosomonas eutropha were obtained and used as the reference template for hope primers ntspa <dig>  ntspa <dig>  nso1225b and nse <dig>  the reference standards for ntspa <dig>  nit <dig>  nmocl6a, nmocl6b and nspcl <dig> were obtained by first amplifying the 16s rrna gene portion targeted by that respective hope primer coupled with 27f.
the reference standards for primers nmo <dig> and nspcl <dig> were obtained by amplifying the 16s rrna gene portion with a sense-strand version of that respective hope primer and 1492r.
all amplicons were gel-excised, purified and submitted for sanger sequencing at the kaust genomics core lab.
sequences were checked to ensure perfect matches in regions targeted by the specific primers.
reference templates for nit <dig> and nso <dig> were obtained by checking the above-mentioned reference templates for perfect match regions.
amplicons with regions that perfectly matched the targeted primers were then individually ligated into the pcr <dig>  vector using the ta cloning kit .
the ligated vectors were individually transformed into top <dig> competent cells, and transformed cells were then extracted for the plasmid using pureyield plasmid miniprep system .
the insert gene was resequenced by the sanger approach to verify for perfect match regions.
subsequently, plasmid dna concentrations were measured using an invitrogen qubit  <dig>  fluorometer  and diluted to 10 ng/μl to be used as hope standards in determination of calibration factors.
calibration factors are normalization factors to account for differences in peak areas between two extended hope primers  <cit> , and are calculated based on eq.
below:  <dig> cfa−b=peak area of extended primera,papeak area of extended primerb,pb where primer b is targeting at a higher hierarchical level compared to primer a.
hope reactions and capillary electrophoresis
for each sample that was to be tested with hope, the bacterial 16s rrna genes were first amplified using 27f coupled with 1492r primer in pcr amplifications comprising of 5 ng dna, 25 μl of epicentre biotechnologies failsafe premix f , 200 nm  of forward and reverse primers,  <dig>  u of extaq dna polymerase , and molecular-biology grade water as necessary to reach a total volume of 50 μl.
pcr was performed with 1 cycle of 95 °c for 30 s followed by 25 cycles of thermal cycling  and a final cycle of 72 °c for 10 min.
all amplicons were purified using wizard sv gel and pcr clean-up kit , measured for their concentrations, and diluted to 10 ng/μl to be used as templates in hope reactions.
each hope standard and sample reaction was carried out in a reaction volume that was comprised of  <dig>  μl of applied biosystems snapshot multiplex kit ,  <dig>  μl of 10 μm primer mix,  <dig>  μl of dna template and  <dig>  μl of molecular-biology grade water.
the hope thermal cycling program consisted of 20 cycles of denaturation , annealing  and extension .
after the primer extension reaction, 5 μl of 200 u recombinant shrimp alkaline phosphatase, rsap , was added to each reaction and the mixture was incubated at 37 °c for 60 min prior to denaturation at 80 °c for 10 min.
the rsap-digested samples were prepared for capillary electrophoresis by first adding 1 μl of the sample to 12 μl of applied biosystems hi-di formamide and  <dig>  μl of genescan  <dig> liz size standard, and the mixed samples were denatured at 96 °c for 5 min prior to capillary injection in the applied biosystems  <dig> series genetic analyzer.
extended primers were identified on genemapper v  <dig>  based on the extended nucleotide base  and fragment size.
peak areas of the primers extended in the presence of reference standards were recorded for determination of calibration factors as described previously  <cit>  and in eq.
.
the same was performed for the primers extended in the presence of template obtained from tested samples.
the relative abundance of 16s rrna gene amplicons targeted by primer a with respect to those targeted by primer b can then be calculated based on eq.
below:  <dig> relative abundance of target%=peak area of primerapeak area of primerb×cfa−b×100% where primer b is targeting at a higher hierarchical level compared to primer a.
16s rrna gene-based high-throughput sequencing and data analysis
for all samples quantified by hope, illumina miseq amplicon sequencing was performed to validate the hope data and to provide information on the total microbial community.
to prepare the 16s rrna gene amplicon libraries, 515f  and 907r  primers were modified to encode the overhang adaptor sequences, and used to amplify the 16s rrna genes.
the thermal cycling program included an initial denaturation stage at 95 °c for 3 min, followed by 25 cycles of denaturation at 95 °c for 30 s, annealing at 55 °c for 30 s and extension at 72 °c for 30 s, and then a final extension stage at 72 °c for 5 min.
pcr amplicons were then purified by ampure xp beads  prior to the index pcr.
nextera xt index  was incorporated into each of the individual samples during pcr.
the thermal cycling program included a denaturation stage at 95 °c for 3 min, followed by 8 cycles of denaturation at 95 °c for 30 s, annealing at 55 °c for 30 s and extension at 72 °c for 30 s, and then a final extension stage at 72 °c for 5 min.
the final indexed pcr amplicons were again purified by ampure xp beads and quantified for the concentrations using an invitrogen qubit  <dig>  fluorometer.
the controls for all pcr reactions were negative for amplification.
purified amplicons were submitted to the kaust genomics core lab for unidirectional sequencing on an illumina miseq platform.
raw sequence reads were handled using procedures described previously  <cit> .
to annotate the 16s rrna gene sequences obtained from high-throughput sequencing, rdp classifier was used for taxonomical assignments at a 95% confidence level  <cit> .
primer-e analysis and statistical tests
the taxonomical assignment for each sample obtained by high-throughput amplicon sequencing was calculated for relative abundances of the individual bacterial and archaeal genera.
special emphasis was made to collate the relative abundance of the nitrospira, nitrosomonas, nitrobacter and nitrosospira so as to provide validation against the relative abundances of these genera obtained by hope method.
in addition, the relative abundance of all taxonomical assignments obtained by miseq amplicon sequencing were collated alongside the relative abundances of the bacterial targets identified by hope.
both datasets were individually input into primer e version  <dig>  <cit> , square-root transformed, and then computed for their bray-curtis similarities.
the bray-curtis distance matrix was used for multivariate analysis on a non-metric threshold multidimensional scaling  plot.
the nmds plots utilizes a distance matrix applied to all samples by which each is represented as a point in the two-dimensional space.
the x and y-axes of the nmds plot do not represent any parameters.
instead, the distance between the positioning of two samples would denote the extent of similarity between these two samples.
to illustrate, samples with higher similarity are closer in proximity in the nmds and vice versa.
the stress value measures the scatter in the nmds plot, and ranges from  <dig> to  <dig>  with  <dig> denoting a good representation of the positioning of the samples in the two-dimensional space shown in the nmds plot and vice versa.
bacterial targets that exhibited > <dig>  correlation with the multivariate patterns on the nmds were overlaid as vectors.
vectors are directional lines emanating from a common origin, extending in the directions in which the marked bacterial targets increase in relative abundance, and hence play a dominant role in the subclustering of samples on the nmds plot.
all other significance tests were conducted using two-tailed t-tests on microsoft excel  <dig>
biofilm reactor setup and operation
the samples that were quantified for the respective nitrifying populations by hope and 16s rrna gene-based amplicon sequencing were obtained from a trickling biofilm nitrification reactor.
the trickling biofilm reactor was first packed with sponge cubes, each of dimension  <dig> by  <dig>  cm, to a total volume of  <dig>  l .
500 ml of sludge collected from the aerobic activated sludge tank of the full scale wwtp was inoculated upon commencement of operation and another 500 ml of aerobic sludge was added 2 weeks later.
during the startup phase , basal substrate containing 40 mg/l ammonium  was fed into the reactor.
the components per l of basal substrate were modified from a previous study  <cit>  and include: 153 mg nh4cl, 480 mg nahco <dig>  24 mg mgso <dig> h2o, 20 mg k2hpo <dig>   <dig>  mg caso <dig> h2o and  <dig>  mg fecl <dig> h2o, 100 μl of trace elements that contained  <dig>  g/l na2moo <dig> h2o,  <dig>  g/l znso <dig> h2o,  <dig>  g/l cuso <dig> h2o,  <dig>  g/l cocl <dig> h2o, 7 g/l feso <dig> h2o and  <dig>  g/l mncl <dig> h2o.
upon reaching stabilization and from day  <dig> onwards, the basal substrate was prepared by diluting the anaerobic effluent that was high in ammonium content  <cit> .
the anaerobic effluent contained on average 160 mg/l nh4
+-n and 40 mg/l of chemical oxygen demand , and the ammonium content in the basal feed was diluted accordingly to provide 80 mg/l nh4
+-n from day  <dig> to  <dig> , 110 mg/l nh4
+-n from day  <dig> to  <dig> , and 150 mg/l nh4
+-n from day  <dig> to  <dig> .
ammonium concentration was subsequently increased to 200 mg/l nh4
+-n from day  <dig> to  <dig> , and diluted to 150 mg/l nh4
+-n from day  <dig> to  <dig> .
throughout operation, the ph of trickling biofilm reactor was maintained between  <dig>  and  <dig>  and aeration was provided at a sparging rate of 2 l/min.fig.
<dig> trickling nitrification biofilm reactor column.
the reactor was operated through six phases a through e, with increasing ammonium  concentration in the influent during phase a through e and a decrease to 150 mg/l nh4
+ in the influent at phase d <dig>  samples were collected at the top  and bottom  locations from port u and b, respectively.
influent and effluent were also collected from the marked sampling ports
biofilm reactor sampling and dna extraction
three sponge cubes were sampled individually from both the top and bottom sampling ports , and referred to as u and b samples subsequently in this study.
new replacement sponges were placed back into the reactor after each biomass sampling so as to provide new substrata for the adherent bacteria.
subsequent samplings were carried out such as to avoid sampling the replacement sponges.
sponges with the adherent biomass were placed in 50 ml centrifuge tubes, each containing 20 ml of 1x phosphate buffer saline .
the samples were then vortexed at maximum speed for 10 min and the sponge cubes removed aseptically by sterile forceps.
the biomass suspensions in the 50 ml tubes were then centrifuged at 6800×g for 10 min to collect cell pellets.
<dig>  g of biomass from cell pellet was weighed and extracted for its dna.
dna was extracted using the ultraclean soil dna isolation kit  with slight modifications to the protocol by adding lysozyme and achromopeptidase to the lysis buffer  <cit> .
full-scale wwtp sampling and dna extraction
another set of samples that were quantified for the respective nitrifying populations by hope and 16s rrna gene-based amplicon sequencing were sampled from a full-scale wastewater treatment plant located in kaust.
the wwtp is a full-scale membrane bioreactor  equipped with the following process units:  grid mesh screen,  primary clarifier,  anoxic-oxic activated sludge tanks,  submerged membrane tank, and  holding tank for chlorination.
influent was collected after the primary clarifier, effluent was collected at the mbr discharge point and chlorinated effluent was collected from the holding tank.
activated sludge was collected in the oxic zone of the sludge tank.
sampling was performed on a monthly basis from july to december  <dig>  with the exception of october  <dig> during which two sample sets were collected.
influent samples were prepared for dna extraction by centrifuging  <dig> to 100 ml of influent at 10000×g for 20 min to obtain the biomass pellet.
2 l each of effluent and chlorinated effluent were individually filtered through  <dig>  μm polycarbonate membranes to retain the biomass.
biomass samples were extracted for their dna using the ultraclean soil dna isolation kit  with slight modifications to the protocol by adding lysozyme and achromopeptidase to the lysis buffer  <cit> .
water quality measurement for biofilm reactor samplesand wastewater streams
samples from the fresh basal substrate fed into the trickling biofilm reactor were collected as influent, and samples collected from the discharge port were collected as effluent .
ammonium , nitrite , nitrate  concentrations in influent and effluent samples  collected every 2– <dig> d throughout the reactor operational period were measured using test ‘n tube high range ammonia kit, tntplus  <dig>  and tntplus  <dig>  respectively.
cod was measured using either lck  <dig>  or lck  <dig> cod  cuvette test vials depending on the concentration to be measured.
all measurements were conducted based on protocols specified by the manufacturer .
the water quality data obtained from the trickling biofilm reactor were collated, log-transformed and normalized within primer-e version  <dig>  <cit> .
the normalized data were then used to generate a principal component analysis  plot.
water quality parameters that exhibited > <dig>  correlation with the multivariate pattern on the pca were overlaid as vectors.
for the wastewater samples obtained from the wwtp, nitrate concentrations in all three streams were measured in lab using the tntplus  <dig>  total nitrogen in all streams as well as ammonia in the influents were measured and provided by the plant operator.
RESULTS
aob abundance in biofilm reactor quantified by hope throughout the operational phases
the relative abundances of nitrifying bacterial populations quantified by the hope approach were collated for multivariate analysis on a multidimensional scaling  plot .
the positioning of the samples on the nmds plot showed that operational phase was the dominant factor accounting for the subclustering of samples.
vector-based analysis showed an increase in the abundance of aob in phases e and d <dig> compared to the early operational phases.
to illustrate, the relative abundance of aob groups, specifically nitrosomonas targeted by nse <dig> and nmocl6a, increased from  <dig>  ±  <dig> % and  <dig>  ±  <dig> % of total bacteria, respectively, in phase b to  <dig>  ±  <dig> % and  <dig>  ±  <dig> % of total bacteria, respectively, in phase d <dig> .
the relative abundances of these aob groups were significantly higher in the latter phases  than in earlier phases  .
both nitrosomonas and nitrosospira targeted by nso <dig> also increased from  <dig>  ±  <dig> % in phase b to  <dig>  ±  <dig> % in phase e and  <dig>  ±  <dig> % in phase d <dig> .fig.
<dig> multivariate analysis of the microbial data.
the relative abundance of aob and nob targeted by hope primers revealed differences throughout the operational phases.
a. the relative abundances of the total microbial community obtained from high-throughput sequencing also showed similar clustering based on the operational phases
nitrosomonas
nitrosomonas with nse1472:
ρ =  <dig> , p =  <dig>  × 10−14
nitrosospira
nitrosospira with nsv443:
ρ =  <dig> , p =  <dig>  × 10−5
nitrosomonas and nitrosospira with nso190:
ρ =  <dig> , p =  <dig>  × 10−10
nitrospira
nitrospira with ntspa712:
ρ =  <dig> , p =  <dig>  × 10−4
nitrospira with ntspa572:
ρ =  <dig> , p =  <dig>  × 10−11
nitrobacter
nitrobacter with nit1017:
ρ = − <dig> , p =  <dig>
nitrobacter with nit1000:
ρ =  <dig> , p =  <dig>  × 10−6
the relative abundance of selected targets further correlated with the amplicon sequencing data
however, the increase in the relative abundance of nso190-targets is mainly due to the increase in the relative abundance of nitrosomonas.
this is because nitrosospira detected by primer nsv <dig> remained stable in its relative abundance against the total community throughout the operational phases .
the changes in the relative abundance of nsv443-targeted nitrosospira were not significantly different across phases b to d <dig> .
abundance of nob and nitrospira-affiliated comammox in biofilm reactor quantified by hope throughout the operational phases
vector analysis performed on the same nmds plot also showed that samples collected from the early operational phases a and b were higher in relative abundance of nob targeted by either ntspa <dig> or ntspa <dig>  both primers targeted nitrospirae at the same hierarchical level but with different coverage .
both primers also target at least one of the three recently identified comammox bacteria .
the average relative abundance of nitrospira targeted by ntspa <dig> was  <dig>  ±  <dig> % in phase a and  <dig>  ±  <dig> % in phase b, and was significantly higher  than that detected in phases c , d <dig> and e , and d <dig> .
similarly, the relative abundance of nitrospira targeted by ntspa <dig> was  <dig>  ±  <dig> % in phase a and  <dig>  ±  <dig> % in phase b, and was significantly higher  than that detected in phases c , d <dig> , e  and d <dig> .
in contrast to nitrospira, nitrobacter detected by the nit <dig> primer increased in its relative abundance from  <dig>  ±  <dig> % in phase b to  <dig>  ±  <dig> % in phase d <dig> .
however, the nit <dig> primer did not perform in a similar manner as the nit <dig> primer despite both primers having similar detection coverage of most nitrobacter spp.
.
the nit <dig> primer consistently failed to detect any nitrobacter in the latter operational phases, and in samples that were detected with nitrobacter, the relative abundance was more than 2-fold lower than that detected by the nit <dig> primer.
validation of hope dataset obtained for trickling biofilm reactor with 16s rrna gene-based amplicon sequencing
based on the hope primer coverage shown in additional file  <dig>  the relative abundance of genera nitrosomonas, nitrosospira, nitrobacter and nitrospira obtained by 16s rrna gene-based amplicon sequencing was respectively compared against the relative abundance detected by nse <dig> , nsv <dig> , nit <dig> or nit <dig> , and ntsp <dig> or ntspa <dig>  hope primers.
the relative abundance of nitrosomonas and nitrospira obtained by 16s rrna gene-based amplicon sequencing was also compared against that detected by nso <dig> hope primer.
generally, the relative abundance obtained by hope differed from that detected by amplicon sequencing but both datasets exhibited good spearman’s rank correlation .
to illustrate, the average relative abundance of nitrosomonas reported by 16s rrna gene-based amplicon sequencing was  <dig> %,  <dig> %,  <dig> %,  <dig> %,  <dig> % and  <dig> % in phase a, b, c, d <dig>  e and d <dig>  respectively .
the average relative abundance of nitrosomonas targeted by hope primer nse <dig> was  <dig> %,  <dig> %,  <dig> %,  <dig> %,  <dig> % and  <dig> % in phase a, b, c, d <dig>  e and d <dig>  respectively.
similarly, the summation of relative abundance of nitrosomonas and nitrosospira obtained by amplicon sequencing differed from that detected by nso <dig> hope primer , although the datasets showed a spearman’s rank correlation coefficient of ρ =  <dig> .
for nob, although nit <dig> detected higher relative abundance of nitrobacter than nit <dig> hope primer, the results obtained from nit <dig> did not correlate positively with those detected by amplicon sequencing .
instead, nitrobacter spp.
were detected in a very low abundance of < <dig> % throughout the operational phases by amplicon sequencing and coincides with the range of relative abundance detected by hope using nit <dig> primer .
in addition, the average relative abundance of nitrospira reported by amplicon sequencing decreased from  <dig> % to  <dig> % from phase a to d <dig>  and the same decrease in relative abundance targeted by ntspa <dig> was observed, albeit with different relative abundance that decreased from  <dig> % to  <dig> %.
changes in biofilm reactor microbial community characterized by 16s rrna gene-based amplicon sequencing
the total microbial communities of the trickling biofilm reactor were characterized by 16s rrna gene-based amplicon sequencing .
the nmds plot generated by this data was similar to that which was obtained by hope in that the multivariate clusters were differentiated according to the operational phases of the reactor .
additional information was also obtained from the amplicon sequencing analysis.
to illustrate, nitrososphaera, which represent ammonia-oxidizing archaea, experienced a significant decrease in the relative abundance from  <dig> % during early operational phases  to negligible levels during the late operational phases .
other microbial populations that significantly decreased in their relative abundances across the different operational phases included unclassified actinobacteria, clostridiaceae i, rhodospirillales and deltaproteobacteria .
unclassified xanthomonadaceae had the highest relative abundance of  <dig>  ±  <dig> % during phases c and d <dig>  and averaged a relative abundance of  <dig>  ±  <dig> % in the remaining phases.
unclassified burkholderiales increased significantly from  <dig>  ±  <dig> % in phases a through d <dig> to  <dig>  ±  <dig> % in phases d <dig> and e .
performance of trickling nitrification biofilm reactor
based on the collated water quality data, effluent from the trickling nitrification biofilm reactor clustered into the six operational phases, which corresponded with the influent ammonium content .
as the influent ammonium content was increased stepwise from 40 mg/l to 200 mg/l, the cod in the effluent also increased and was significantly higher in the effluent collected during phase d <dig>  d <dig> and e compared to those collected during phase a and c .
ammonium removal was optimally achieved at an average of 44% reduction per day during days  <dig> through  <dig> of operation , which corresponded with an influent concentration of 80–110 mg/l .
during these phases, ammonium was primarily converted to nitrate  with minimal nitrite accumulation in the effluent .
subsequently, the trickling biofilm reactor experienced deterioration in its nitrification efficacy as influent ammonium content increased to 150 mg/l and 200 mg/l .
ammonium was only partially converted to nitrate and an accumulation of nitrite was observed in the effluent collected during these phases .
nitrite levels in the effluent during the two latter phases were significantly higher than during phases a through c .fig.
<dig> principal component analysis  of the measured water quality parameters in the effluent stream throughout the operational phases.
the pca plot suggests that there was high ammonium nh4
+ oxidation during phases b and c but an accumulation of nitrite and an increase in the chemical oxygen demand  during the late operational phases d <dig> and e
fig.
<dig> correlative trends between effluent quality of the trickling nitrification biofilm reactor and the abundance ratio of predominant aob and nob.
concentrations of ammonium , nitrite , nitrate  and chemical oxygen demand  in the effluent stream of the trickling nitrification biofilm reactor column throughout the different operational phases .
the ratio of nse1472-, nso190- and nsv443-targeted aob against ntspa712-targeted nob, obtained using hope .
the ratio of nitrosospira and nitrosomonas against nitrospira, obtained using amplicon sequencing
correlative trends of aob and nob proportions with reactor performance
the ratios in the relative abundance of these aob  and nob targeted by ntspa <dig> hope primer were further evaluated and compared to the reactor performance based on its water quality parameters .
it was determined that the average ratio of aob targeted by nse <dig> and nso <dig> normalized against ntspa712-targeted nob was  <dig>  and  <dig> , respectively, during the phases that corresponded to stable reactor performance as evidenced by low nitrite content in the effluent .
during phases e and d <dig>  which corresponded with deterioration in the reactor performance, the ratio of aob targeted by nse <dig> and nso <dig> normalized against ntspa712-targeted nob increased to  <dig>  and  <dig> , respectively.
there was, however, a less significant increase in the ratio of aob targeted by nsv <dig> against nit1017-targeted nob .
these ratios obtained using hope were compared against those generated from 16s rrna-based amplicon sequencing data .
the ratio of the relative abundance of nitrosomonas compared against nitrospira  as determined by amplicon sequencing was  <dig>  during the stable operational phases.
this ratio markedly increased to  <dig> when the reactor deteriorated in performance.
likewise, the ratio of nitrosospira normalized against nitrospira was observed to be  <dig>  during phase a to c, and increased to  <dig>  in the latter phases.
detection of aob and nob in wwtp samples using hope and amplicon sequencing
both hope and amplicon sequencing revealed that the activated sludge from the local full-scale wwtp had lower relative abundances of nitrosomonas, nitrosospira and nitrobacter compared to ntspa712-targeted nitrospira .
nitrate concentrations were consequentially higher in the effluent  than in the influent and chlorinated effluent.
hope and amplicon sequencing were further used to evaluate the nitrifying bacterial populations present in these wastewater samples.
similar to the activated sludge, the wastewaters did not have any detectable nitrobacter but instead, had relatively higher abundances of nitrospira compared to nitrosomonas and nitrosospira.
to illustrate, the relative abundance of ntspa712-targeted nob ranged from  <dig> % in the influent to  <dig> % in the chlorinated effluent.
this relative abundance in nitrospira was at least  <dig> times higher than the aob targeted by nse <dig> and nsv <dig> .
amplicon sequencing also showed a relatively higher abundances of nitrospira compared to both aob  in the three types of wastewater streams.table  <dig> water quality for the influent, effluent and chlorinated effluent samples collected from the wastewater treatment plant and the corresponding relative abundances of aob and nob as quantified using the hope and amplicon sequencing approach
nitrosomonas
nitrosospira
nitrobacter
nitrospira
relative abundance  with respect to eub338ia
discussion
this study demonstrates the development of hope as a high-throughput method to quantify for the nitrifying bacterial populations including aob, nob and nitrospira-affiliated comammox.
along with the aoa, which was not targeted by the hope approach in this study, these nitrifying populations play an important role in the activated sludge process to convert ammonia in the untreated wastewater to nitrite and nitrate.
the nitrate is subsequently converted to dinitrogen and nitrogen by denitrifiers.
collectively, the nitrifying and denitrifying populations reduce the total nitrogenous content in wastewater to a level permissible for discharge or reuse.
it has previously been hypothesized that the theoretical aob/nob ratio in a nitrification process would be  <dig> due to the differences in electron generation and biomass yield between aob and nob .
simultaneous nitrification/denitrification can possibly increase this ratio to  <dig>   <cit> .
however, the inverse may have a potentially detrimental impact on overall functionality of the nitrification system, especially since excess aob can lead to accumulation of the hydroxylamine intermediate which is inhibitory towards nob  <cit> .
as such, the hope method is developed with the intention to serve as a practically applicable tool for monitoring the aob and nob ratio and hence used to infer the efficiency of the nitrification process in wastewater treatment systems.
for this purpose, the hope primers were designed to target nitrosomonas, nitrosospira, nitrobacter and nitrospira at the order, subcluster and/or genus level, and the relative abundance at each hierarchical taxonomical level was normalized against that of the total bacteria.
primers were arranged into multiple reaction tubes that involve multiplexing of up to four targets per reaction, and tested against samples collected from a trickling nitrification biofilm reactor.
it was observed from the hope data that the relative abundance of nitrosomonas increased with increasing nh4
+ concentrations while abundances of nitrosospira did not.
in contrast to the aob, the nob genus nitrospira was detrimentally affected by the shock loading of the reactor.
during the stable phases b and c, nitrospira targeted by the ntspa572-primer was more predominant than nitrobacter.
however, with the accumulation of nitrite in the late operational phases, nitrospira was depleted while nitrobacter, targeted by the nit1000-primer, remained relatively stable in its relative abundance against the total bacteria.
this observation coincides with results reported by knapp and graham  <cit> , in which nitrospira was observed to be more perturbed than nitrobacter during a destabilizing experiment  <cit> .
in addition, given that both ntspa <dig> and ntspa <dig> hope primers target at least one of the three recently identified comammox , the higher relative abundance of hope-targeted nitrospira than nitrobacter is likely to suggest that comammox is contributing to an equally important nitrification role in this reactor.
similarly, the depletion in the collective relative abundance of ntspa572-targeted nitrospira suggests the susceptibility of comammox to shock loading events.
the high nh4
+ content exposure appeared to affect specific groups of aob and nob differently.
it is generally thought that nitrospira-like bacteria are k-strategists with a high affinity for no2
− and oxygen, and reach high densities under substrate-limited conditions.
this is in contrast to the nitrobacter species, which are r-strategists with a lower no2
− and oxygen affinity and outcompete nitrospira only at higher substrate concentrations .
an additional postulation to account for the better survival of nitrobacter during shock loading events is its greater metabolic diversity under stressed conditions  <cit> .
this genus exhibits versatile metabolism and is able to grow either mixotrophically or chemoorganotrophically  <cit> .
although nitrobacter was able to withstand shock loading and increase in its relative abundance as quantified by hope, it was observed that the increase in its relative abundance did not correlate with any decrease in nitrite content of the effluent.
furthermore, the increase in nitrobacter did not correlate with a subsequent conversion of nitrite to nitrate in the effluent.
this suggests that nitrobacter may not be playing as important of a role as nitrospira and comammox in nitrite oxidation in the wastewater treatment process.
given that both nitrosomonas and nitrospira seem to correlate with nitrification functionality, the hope-obtained relative abundances of both groups were further expressed as a proportional ratio.
the ratio of nitrosomonas targeted by either nse <dig> or nso <dig> hope primer against ntspa572-targeted nitrospira ranged from  <dig>  to  <dig>  in the nitrification biofilm reactor during the phases that corresponded to stable reactor performance as evidenced by the low nitrite content in the effluent.
however, the exceedingly high ammonium content of the later phases of operation  resulted in a toxic shock response by the nitrifying bacterial populations, which correlated with an observed increase in the aob/nob ratio.
this ratio increased to a range of  <dig>  to  <dig>  when the bioreactor performance deteriorated.
based on energetic calculations, it was proposed that a theoretical aob/nob ratio of  <dig> should be obtained during a stable and functional nitrification process   <cit> .
this theoretical ratio was in agreement with the observation made in this study, and reiterates the need to maintain an optimal proportion of aob to nob in a well-performing nitrification process.
the hope method had been previously validated against other conventional methods such as qpcr and fish  <cit> .
in this study, to validate the data obtained from hope, the relative abundances of nitrifying bacterial populations in the trickling biofilter and in samples collected from a full-scale wwtp were further validated against results from amplicon sequencing .
the observed difference in the abundance values between hope and amplicon sequencing is likely due to differences in the primer coverage for the individual aob and nob groups .
for example, nso <dig> targets  <dig> % of the nitrosomonas and  <dig> % of the nitrosospira sequences collated in the rdp classifier database.
in contrast, the universal primers 515f and 907r used for amplicon sequencing target a higher percentage of nitrosomonas and nitrosospira.
furthermore, the relative abundances reported by hope were expressed by normalizing against total bacteria as detected by universal primers that differed from those used in amplicon sequencing.
although the relative abundances determined by both methods varied in terms of the absolute values, both datasets showed good spearman’s rank correlation.
furthermore, multivariate nmds plots from both hope and amplicon sequencing also showed good correlation at a significant confidence level .
this indicates that the observed trends in the nitrifying bacterial population dynamics revealed by both methods were well aligned.
although similar conclusions were obtained from both hope and amplicon sequencing with regards to the nitrifying population dynamics in the presence of shock loading, a limitation of the hope approach is that it can only detect bacterial targeted by the designed primers.
for instance, through the use of amplicon sequencing and not hope, it was determined that the detrimental effect of the shock loading event on nitrospira can possibly be due to its competition for substrates  with heterotrophic bacteria.
such heterotrophic bacteria included unclassified burkholderiales and unclassified xanthomonadaceae, which increased in relative abundance throughout all operational phase.
furthermore, the hope technique developed in this study does not contain primers that target the aoa.
although aoa had been previously found abundant in wastewater treatment systems  <cit> , both genus nitrososphaera and nitrosopumilus were found to be in low abundance  in the samples when evaluated by high-throughput amplicon sequencing.
further optimization on the hope technique would be required to design a range of aoa-targeting primers, as well as to improve its current detection limits of  <dig> % so as to facilitate the quantification of potentially low abundance species.
regardless, an advantage of the hope approach is that new primers, for example those designed to target the different sublineages of nitrospira and nitrosomonas that were selected for by varying nitrite and ammonium concentrations , can be easily added to the multiplexing tube arrangement when these primers become available.
the additional depth offered by high-throughput amplicon sequencing would be useful as a tool for periodically benchmarking hope-based results obtained on a more frequent basis.
although this study did not attempt to challenge the limits of high-throughput multiplexing capability of hope, the method can in theory perform up to 32-plexes per reaction tube, and modifications to the hope approach can be efficiently made to facilitate the simultaneous detection of the expanding list of nitrifying bacteria.
furthermore, this study demonstrated a good correlation between the results obtained from the hope approach and high-throughput sequencing.
thus, depending on the intended purpose of an experimental or practical study, hope may be particularly useful for the frequent tracking of the occurrence and abundance of selected microbial populations without the need to perform high-throughput amplicon sequencing.
this in turn reduces the need to sieve through a full range of microbial population data for every sampling event.
CONCLUSIONS
in summary, this study demonstrates the applicability and adaptability of hope for assessing abundances of predominant aob and nob groups that have been identified as integral to wastewater treatment systems.
this method allows for the simultaneous monitoring of relative abundances of aob, nob and comammox groups, which can be used to provide indicative data of nitrification performance and efficiency.
coupled with previous application of hope to a variety of target bacteria and sample types , this study demonstrates the versatility and applicability of hope for a wide range of microbial ecology studies.
additional file
additional file 1: target coverage for the different primers used in this study.
abbreviations
aoaammonia-oxidizing archaea
aobammonia-oxidizing bacteria
codchemical oxygen demand
ddntpdideoxynucleoside triphosphate
dpcrdigital polymerase chain reaction
fishfluorescent in-situ hybridization
hopehierarchical oligonucleotide primer extension
mbrmembrane bioreactor
nh4+-nammonium
nmdsnon-metric threshold multidimensional scaling
no2--nnitrite
no3--nnitrate
nobnitrite-oxidizing bacteria
ntnucleotides
pbsphosphate buffer saline
pcaprincipal component analysis
qpcrquantitative polymerase chain reaction
rdpribosomal database project
rsaprecombinant shrimp alkaline phosphatase
tntotal nitrogen
wwtpwastewater treatment plant

 the conventional process of removing nitrogen from municipal wastewater streams involves nitrification followed by denitrification. nitrification in the wastewater treatment process is generally performed in two consecutive steps. the first step is ammonia oxidation by either ammonia-oxidizing bacteria or ammonia-oxidizing archaea, and the second step is nitrite oxidation by nitrite-oxidizing bacteria. representative aob genera include nitrosomonas and nitrosospira while nitrososphaera and nitrosopumilus would account as two examples of aoa. representative nob include nitrospira and nitrobacter. in recent years, three members within the genus nitrospira, namely candidatus nitrospira inopinata, candidatus nitrospira nitrosa and candidatus nitrospira nitrificans, were shown to perform complete nitrification in a single step and are referred to as comammox. although comammox are not yet isolated from municipal wastewater treatment plants and their contribution in terms of full nitrification in such systems remain unknown, comammox had been shown to use urea as ammonium source for nitrification and that this trait should in theory enable them to thrive in wwtp environments where urea is often present. the commensal interaction between aob/aoa and nob, and the potential importance of nitrospira-affiliated comammox suggest that the nitrification process can be prone to failure when there is a suboptimal ratio of aob/aoa to nob, particularly that of nitrospira. this is primarily due to the effect of hydroxylamine, an intermediate of the conversion of ammonia to nitrite, which has a strong inhibitory effect on nitrite oxidizers when released by aob/aoa. failure of the nitrification system in a wwtp can result in discharge of effluent that contains total nitrogen higher than the permissible limit. this can result in potentially detrimental environmental impacts in the form of eutrophication of the receiving water bodies. hence, establishing an optimal ratio and abundance of aob/aoa to nob is important in ensuring the functionality of the nitrification process so that treated effluent will meet the discharge limits for tn. in such instance, a quantitative method capable of efficiently targeting aob/aoa and nob relative abundances would be particularly useful in tracking the proportional ratio between these two groups. various methods like quantitative pcr and fluorescent in-situ hybridization were used to directly quantify the abundances of aob and nob. through these approaches, it was determined that nitrosomonas likely accounts for the largest fraction of aob in conventional activated sludge systems. additionally, in conventional municipal wastewater streams in which ammonium concentrations ranged from to 50 mg/l, nitrospira was the predominant nob while nitrobacter was generally in low abundance. although qpcr and fish have been routinely used in earlier studies, both approaches have substantial limitations. for example, qpcr requires standard curves to be established for individual targets. new generation of pcr does not require standard curves to be generated but similar to conventional qpcr, remains limited by the number of fluorophores available in the ultraviolet spectrum. hence, both quantitative pcr approaches have limited multiplexing throughput capability. furthermore, both qpcr approaches can be prone to non-specific amplification when sybr® green fluorescence reporters are utilized, and would require an additional melting curve analysis or verification step to differentiate between true-positive and false-positive amplifications. although taqman® probes can be developed for use in qpcr and dpcr, such probes are relatively expensive and remain challenging to design for a specific microbial target within the internal region of two flanking qpcr primers. digital pcr being a considerably new technology, also remains rather costly till this date compared to traditional qpcr. fish, on the other hand, can be affected by the in-situ accessibility of labeled probes through the cell membranes and the copy numbers of rrna genes inside aob and nob. for example, species of nitrobacter, nitrosomonas, nitrosospira and nitrospira are reported to have one copy of 16s rrna gene per cell. this copy number of 16s rrna genes per genome is lower than the mean copy number of among the microbial genomes currently available in the same ribosomal rna operon copy number database. hierarchical oligonucleotide primer extension was developed as a method to complement other quantitative methods like qpcr and fish by providing a high-throughput multiplexing platform capable of addressing some of the previous methods’ limitations. in this method, oligonucleotide primers are designed to target genes, usually 16s rrna genes, at different hierarchical taxonomical levels. in the presence of a dna template, the primers anneal at the designed targeting position of the 16s rrna gene and extend with a single fluorescently labeled dideoxynucleoside triphosphate. the extended primers can then be identified on a dna sequencer based on their anticipated fragment size and fluorescence. the peak areas of the individual fragments can then be used to determine the relative abundances of bacterial targets of interest. the hope method has been developed to detect specific bacterial populations of interest including host-associated bacteroidales and bifidobacterium, cyanobacteria, and methanogens. in these studies, hope has been demonstrated to multiplex up to primers in a single reaction, allowing for the evaluation of a range of relevant bacterial indicators in a time-efficient manner. the high-throughput capability of hope stems from the use of different fragment sizes of primers and fluorescence colors to distinguish between extended primers. to exemplify, most primers are synthesized with a length of to nucleotides while most dna sequencer platforms are able to detect up to four different fluorescence colors. hence, in theory, the combination of these two differentiating factors would mean that hope is capable of multiplex up to primers in one single reaction tube. aside from its multiplexing capability, it was demonstrated that hope has a detection limit of % of the total pcr-amplified bacterial targets, and that the entire duration required to carry out a single hope reaction and to analyze for the extended primer on a dna sequencer is <1 h, which would be comparable to the time needed to carry out qpcr. given the advantages of hope, the method is a potentially useful tool in tracking the relative abundance of nitrifying populations. in the present study, we applied the hope approach to target groups of aob and nob in wastewater environments. the method was tested against samples collected from a lab-scale biofilm-based trickling bioreactor. the total microbial community, including the nitrifying groups, was analyzed by 16s rrna gene-based high-throughput sequencing for comparison against the hope data. to further demonstrate the use of hope in monitoring these relative abundances, the trickling bioreactor was subsequently challenged with high concentrations of incoming ammonium content so as to simulate a toxic shock event that could potentially crash the nitrification system. the changes in aob and nob proportions were then tracked by the hope approach and correlated to reactor performance. wastewater samples from a fully operational wastewater treatment plant were also analyzed in order to further validate the method. methods hope primer layout hope was performed with a total of primers targeting ammonia-oxidizing genera nitrosomonas and nitrosospira as well as nitrite-oxidizing genera nitrospira and nitrobacter. primers were sourced from the probebase database. the mismatch positions of non-targets were moved to the 3′-end of the hope primers as it was previously found that mismatches located at the 3′-end of the hope primers would not facilitate nucleotide base extension of primers, and would hence lead to improved specificity. re-positioning of the mismatches result in minimal change to the original coverage of the targeted bacterial groups. the coverage of the primers was then verified in-silico against the ribosomal database project database, and shown in additional file primers were arranged into four tubes, with either 27f or 338f primer included in each tube to provide a normalization of the abundance of each nitrifying bacterial target against the total bacteria. to differentiate primers that extend with the same ddntp and hence color during capillary electrophoresis, some primers were modified with a polya tail at the 5′-end of the primer. primers in each tube were prepared to form a final concentration of 10 μm.table hope primers used to target the different ammonia- and nitrite-oxidizing bacterial groups nob tube 1: nitrospira nob tube 2: nitrobacter aob tube 1: nitrosomonas and nitrosospira aob tube 2: nitrosomonas and nitrosospira target coverage for each primer was identified using the rdp probe match function hope reference template preparation to prepare the hope reference templates for nitrifying bacterial populations, extracted dna from a separate trickling nitrification biofilter that was set up in year was first amplified with 27f to 1492r primers and purified for use in 16s rrna gene clone library construction. approximately clones were picked and sequenced to identify nitrifying bacterial species. two clones associated with candidatus nitrospira defluvii and nitrosomonas eutropha were obtained and used as the reference template for hope primers ntspa ntspa nso1225b and nse the reference standards for ntspa nit nmocl6a, nmocl6b and nspcl were obtained by first amplifying the 16s rrna gene portion targeted by that respective hope primer coupled with 27f. the reference standards for primers nmo and nspcl were obtained by amplifying the 16s rrna gene portion with a sense-strand version of that respective hope primer and 1492r. all amplicons were gel-excised, purified and submitted for sanger sequencing at the kaust genomics core lab. sequences were checked to ensure perfect matches in regions targeted by the specific primers. reference templates for nit and nso were obtained by checking the above-mentioned reference templates for perfect match regions. amplicons with regions that perfectly matched the targeted primers were then individually ligated into the pcr vector using the ta cloning kit. the ligated vectors were individually transformed into top competent cells, and transformed cells were then extracted for the plasmid using pureyield plasmid miniprep system. the insert gene was resequenced by the sanger approach to verify for perfect match regions. subsequently, plasmid dna concentrations were measured using an invitrogen qubit fluorometer and diluted to 10 ng/μl to be used as hope standards in determination of calibration factors. calibration factors are normalization factors to account for differences in peak areas between two extended hope primers, and are calculated based on eq. below: cfa−b=peak area of extended primera,papeak area of extended primerb,pb where primer b is targeting at a higher hierarchical level compared to primer a. hope reactions and capillary electrophoresis for each sample that was to be tested with hope, the bacterial 16s rrna genes were first amplified using 27f coupled with 1492r primer in pcr amplifications comprising of 5 ng dna, 25 μl of epicentre biotechnologies failsafe premix f, 200 nm of forward and reverse primers,  u of extaq dna polymerase, and molecular-biology grade water as necessary to reach a total volume of 50 μl. pcr was performed with 1 cycle of 95 °c for 30 s followed by 25 cycles of thermal cycling and a final cycle of 72 °c for 10 min. all amplicons were purified using wizard sv gel and pcr clean-up kit, measured for their concentrations, and diluted to 10 ng/μl to be used as templates in hope reactions. each hope standard and sample reaction was carried out in a reaction volume that was comprised of  μl of applied biosystems snapshot multiplex kit,  μl of 10 μm primer mix,  μl of dna template and  μl of molecular-biology grade water. the hope thermal cycling program consisted of 20 cycles of denaturation, annealing and extension. after the primer extension reaction, 5 μl of 200 u recombinant shrimp alkaline phosphatase, rsap, was added to each reaction and the mixture was incubated at 37 °c for 60 min prior to denaturation at 80 °c for 10 min. the rsap-digested samples were prepared for capillary electrophoresis by first adding 1 μl of the sample to 12 μl of applied biosystems hi-di formamide and  μl of genescan liz size standard, and the mixed samples were denatured at 96 °c for 5 min prior to capillary injection in the applied biosystems series genetic analyzer. extended primers were identified on genemapper v based on the extended nucleotide base and fragment size. peak areas of the primers extended in the presence of reference standards were recorded for determination of calibration factors as described previously and in eq.. the same was performed for the primers extended in the presence of template obtained from tested samples. the relative abundance of 16s rrna gene amplicons targeted by primer a with respect to those targeted by primer b can then be calculated based on eq. below: relative abundance of target%=peak area of primerapeak area of primerb×cfa−b×100% where primer b is targeting at a higher hierarchical level compared to primer a. 16s rrna gene-based high-throughput sequencing and data analysis for all samples quantified by hope, illumina miseq amplicon sequencing was performed to validate the hope data and to provide information on the total microbial community. to prepare the 16s rrna gene amplicon libraries, 515f and 907r primers were modified to encode the overhang adaptor sequences, and used to amplify the 16s rrna genes. the thermal cycling program included an initial denaturation stage at 95 °c for 3 min, followed by 25 cycles of denaturation at 95 °c for 30 s, annealing at 55 °c for 30 s and extension at 72 °c for 30 s, and then a final extension stage at 72 °c for 5 min. pcr amplicons were then purified by ampure xp beads prior to the index pcr. nextera xt index was incorporated into each of the individual samples during pcr. the thermal cycling program included a denaturation stage at 95 °c for 3 min, followed by 8 cycles of denaturation at 95 °c for 30 s, annealing at 55 °c for 30 s and extension at 72 °c for 30 s, and then a final extension stage at 72 °c for 5 min. the final indexed pcr amplicons were again purified by ampure xp beads and quantified for the concentrations using an invitrogen qubit fluorometer. the controls for all pcr reactions were negative for amplification. purified amplicons were submitted to the kaust genomics core lab for unidirectional sequencing on an illumina miseq platform. raw sequence reads were handled using procedures described previously. to annotate the 16s rrna gene sequences obtained from high-throughput sequencing, rdp classifier was used for taxonomical assignments at a 95% confidence level. primer-e analysis and statistical tests the taxonomical assignment for each sample obtained by high-throughput amplicon sequencing was calculated for relative abundances of the individual bacterial and archaeal genera. special emphasis was made to collate the relative abundance of the nitrospira, nitrosomonas, nitrobacter and nitrosospira so as to provide validation against the relative abundances of these genera obtained by hope method. in addition, the relative abundance of all taxonomical assignments obtained by miseq amplicon sequencing were collated alongside the relative abundances of the bacterial targets identified by hope. both datasets were individually input into primer e version, square-root transformed, and then computed for their bray-curtis similarities. the bray-curtis distance matrix was used for multivariate analysis on a non-metric threshold multidimensional scaling plot. the nmds plots utilizes a distance matrix applied to all samples by which each is represented as a point in the two-dimensional space. the x and y-axes of the nmds plot do not represent any parameters. instead, the distance between the positioning of two samples would denote the extent of similarity between these two samples. to illustrate, samples with higher similarity are closer in proximity in the nmds and vice versa. the stress value measures the scatter in the nmds plot, and ranges from to with denoting a good representation of the positioning of the samples in the two-dimensional space shown in the nmds plot and vice versa. bacterial targets that exhibited > correlation with the multivariate patterns on the nmds were overlaid as vectors. vectors are directional lines emanating from a common origin, extending in the directions in which the marked bacterial targets increase in relative abundance, and hence play a dominant role in the subclustering of samples on the nmds plot. all other significance tests were conducted using two-tailed t-tests on microsoft excel biofilm reactor setup and operation the samples that were quantified for the respective nitrifying populations by hope and 16s rrna gene-based amplicon sequencing were obtained from a trickling biofilm nitrification reactor. the trickling biofilm reactor was first packed with sponge cubes, each of dimension by  cm, to a total volume of  l. 500 ml of sludge collected from the aerobic activated sludge tank of the full scale wwtp was inoculated upon commencement of operation and another 500 ml of aerobic sludge was added 2 weeks later. during the startup phase, basal substrate containing 40 mg/l ammonium was fed into the reactor. the components per l of basal substrate were modified from a previous study and include: 153 mg nh4cl, 480 mg nahco 24 mg mgso h2o, 20 mg k2hpo  mg caso h2o and  mg fecl h2o, 100 μl of trace elements that contained  g/l na2moo h2o,  g/l znso h2o,  g/l cuso h2o,  g/l cocl h2o, 7 g/l feso h2o and  g/l mncl h2o. upon reaching stabilization and from day onwards, the basal substrate was prepared by diluting the anaerobic effluent that was high in ammonium content. the anaerobic effluent contained on average 160 mg/l nh4 +-n and 40 mg/l of chemical oxygen demand, and the ammonium content in the basal feed was diluted accordingly to provide 80 mg/l nh4 +-n from day to, 110 mg/l nh4 +-n from day to, and 150 mg/l nh4 +-n from day to. ammonium concentration was subsequently increased to 200 mg/l nh4 +-n from day to, and diluted to 150 mg/l nh4 +-n from day to. throughout operation, the ph of trickling biofilm reactor was maintained between and and aeration was provided at a sparging rate of 2 l/min.fig. trickling nitrification biofilm reactor column. the reactor was operated through six phases a through e, with increasing ammonium concentration in the influent during phase a through e and a decrease to 150 mg/l nh4 + in the influent at phase d samples were collected at the top and bottom locations from port u and b, respectively. influent and effluent were also collected from the marked sampling ports biofilm reactor sampling and dna extraction three sponge cubes were sampled individually from both the top and bottom sampling ports, and referred to as u and b samples subsequently in this study. new replacement sponges were placed back into the reactor after each biomass sampling so as to provide new substrata for the adherent bacteria. subsequent samplings were carried out such as to avoid sampling the replacement sponges. sponges with the adherent biomass were placed in 50 ml centrifuge tubes, each containing 20 ml of 1x phosphate buffer saline. the samples were then vortexed at maximum speed for 10 min and the sponge cubes removed aseptically by sterile forceps. the biomass suspensions in the 50 ml tubes were then centrifuged at 6800×g for 10 min to collect cell pellets.  g of biomass from cell pellet was weighed and extracted for its dna. dna was extracted using the ultraclean soil dna isolation kit with slight modifications to the protocol by adding lysozyme and achromopeptidase to the lysis buffer. full-scale wwtp sampling and dna extraction another set of samples that were quantified for the respective nitrifying populations by hope and 16s rrna gene-based amplicon sequencing were sampled from a full-scale wastewater treatment plant located in kaust. the wwtp is a full-scale membrane bioreactor equipped with the following process units: grid mesh screen, primary clarifier, anoxic-oxic activated sludge tanks, submerged membrane tank, and holding tank for chlorination. influent was collected after the primary clarifier, effluent was collected at the mbr discharge point and chlorinated effluent was collected from the holding tank. activated sludge was collected in the oxic zone of the sludge tank. sampling was performed on a monthly basis from july to december with the exception of october during which two sample sets were collected. influent samples were prepared for dna extraction by centrifuging to 100 ml of influent at 10000×g for 20 min to obtain the biomass pellet. 2 l each of effluent and chlorinated effluent were individually filtered through  μm polycarbonate membranes to retain the biomass. biomass samples were extracted for their dna using the ultraclean soil dna isolation kit with slight modifications to the protocol by adding lysozyme and achromopeptidase to the lysis buffer. water quality measurement for biofilm reactor samplesand wastewater streams samples from the fresh basal substrate fed into the trickling biofilm reactor were collected as influent, and samples collected from the discharge port were collected as effluent. ammonium, nitrite, nitrate concentrations in influent and effluent samples collected every 2– d throughout the reactor operational period were measured using test ‘n tube high range ammonia kit, tntplus and tntplus respectively. cod was measured using either lck or lck cod cuvette test vials depending on the concentration to be measured. all measurements were conducted based on protocols specified by the manufacturer. the water quality data obtained from the trickling biofilm reactor were collated, log-transformed and normalized within primer-e version. the normalized data were then used to generate a principal component analysis plot. water quality parameters that exhibited > correlation with the multivariate pattern on the pca were overlaid as vectors. for the wastewater samples obtained from the wwtp, nitrate concentrations in all three streams were measured in lab using the tntplus total nitrogen in all streams as well as ammonia in the influents were measured and provided by the plant operator.  aob abundance in biofilm reactor quantified by hope throughout the operational phases the relative abundances of nitrifying bacterial populations quantified by the hope approach were collated for multivariate analysis on a multidimensional scaling plot. the positioning of the samples on the nmds plot showed that operational phase was the dominant factor accounting for the subclustering of samples. vector-based analysis showed an increase in the abundance of aob in phases e and d compared to the early operational phases. to illustrate, the relative abundance of aob groups, specifically nitrosomonas targeted by nse and nmocl6a, increased from  ±  % and  ±  % of total bacteria, respectively, in phase b to  ±  % and  ±  % of total bacteria, respectively, in phase d. the relative abundances of these aob groups were significantly higher in the latter phases than in earlier phases. both nitrosomonas and nitrosospira targeted by nso also increased from  ±  % in phase b to  ±  % in phase e and  ±  % in phase d.fig. multivariate analysis of the microbial data. the relative abundance of aob and nob targeted by hope primers revealed differences throughout the operational phases. a. the relative abundances of the total microbial community obtained from high-throughput sequencing also showed similar clustering based on the operational phases nitrosomonas nitrosomonas with nse1472: ρ = , p =   × 10−14 nitrosospira nitrosospira with nsv443: ρ = , p =   × 10−5 nitrosomonas and nitrosospira with nso190: ρ = , p =   × 10−10 nitrospira nitrospira with ntspa712: ρ = , p =   × 10−4 nitrospira with ntspa572: ρ = , p =   × 10−11 nitrobacter nitrobacter with nit1017: ρ = −, p =  nitrobacter with nit1000: ρ = , p =   × 10−6 the relative abundance of selected targets further correlated with the amplicon sequencing data however, the increase in the relative abundance of nso190-targets is mainly due to the increase in the relative abundance of nitrosomonas. this is because nitrosospira detected by primer nsv remained stable in its relative abundance against the total community throughout the operational phases. the changes in the relative abundance of nsv443-targeted nitrosospira were not significantly different across phases b to d. abundance of nob and nitrospira-affiliated comammox in biofilm reactor quantified by hope throughout the operational phases vector analysis performed on the same nmds plot also showed that samples collected from the early operational phases a and b were higher in relative abundance of nob targeted by either ntspa or ntspa both primers targeted nitrospirae at the same hierarchical level but with different coverage. both primers also target at least one of the three recently identified comammox bacteria. the average relative abundance of nitrospira targeted by ntspa was  ±  % in phase a and  ±  % in phase b, and was significantly higher than that detected in phases c, d and e, and d. similarly, the relative abundance of nitrospira targeted by ntspa was  ±  % in phase a and  ±  % in phase b, and was significantly higher than that detected in phases c, d, e and d. in contrast to nitrospira, nitrobacter detected by the nit primer increased in its relative abundance from  ±  % in phase b to  ±  % in phase d. however, the nit primer did not perform in a similar manner as the nit primer despite both primers having similar detection coverage of most nitrobacter spp.. the nit primer consistently failed to detect any nitrobacter in the latter operational phases, and in samples that were detected with nitrobacter, the relative abundance was more than 2-fold lower than that detected by the nit primer. validation of hope dataset obtained for trickling biofilm reactor with 16s rrna gene-based amplicon sequencing based on the hope primer coverage shown in additional file the relative abundance of genera nitrosomonas, nitrosospira, nitrobacter and nitrospira obtained by 16s rrna gene-based amplicon sequencing was respectively compared against the relative abundance detected by nse, nsv, nit or nit, and ntsp or ntspa hope primers. the relative abundance of nitrosomonas and nitrospira obtained by 16s rrna gene-based amplicon sequencing was also compared against that detected by nso hope primer. generally, the relative abundance obtained by hope differed from that detected by amplicon sequencing but both datasets exhibited good spearman’s rank correlation. to illustrate, the average relative abundance of nitrosomonas reported by 16s rrna gene-based amplicon sequencing was %, %, %, %, % and % in phase a, b, c, d e and d respectively. the average relative abundance of nitrosomonas targeted by hope primer nse was %, %, %, %, % and % in phase a, b, c, d e and d respectively. similarly, the summation of relative abundance of nitrosomonas and nitrosospira obtained by amplicon sequencing differed from that detected by nso hope primer, although the datasets showed a spearman’s rank correlation coefficient of ρ = . for nob, although nit detected higher relative abundance of nitrobacter than nit hope primer, the results obtained from nit did not correlate positively with those detected by amplicon sequencing. instead, nitrobacter spp. were detected in a very low abundance of < % throughout the operational phases by amplicon sequencing and coincides with the range of relative abundance detected by hope using nit primer. in addition, the average relative abundance of nitrospira reported by amplicon sequencing decreased from % to % from phase a to d and the same decrease in relative abundance targeted by ntspa was observed, albeit with different relative abundance that decreased from % to %. changes in biofilm reactor microbial community characterized by 16s rrna gene-based amplicon sequencing the total microbial communities of the trickling biofilm reactor were characterized by 16s rrna gene-based amplicon sequencing. the nmds plot generated by this data was similar to that which was obtained by hope in that the multivariate clusters were differentiated according to the operational phases of the reactor. additional information was also obtained from the amplicon sequencing analysis. to illustrate, nitrososphaera, which represent ammonia-oxidizing archaea, experienced a significant decrease in the relative abundance from % during early operational phases to negligible levels during the late operational phases. other microbial populations that significantly decreased in their relative abundances across the different operational phases included unclassified actinobacteria, clostridiaceae i, rhodospirillales and deltaproteobacteria. unclassified xanthomonadaceae had the highest relative abundance of  ±  % during phases c and d and averaged a relative abundance of  ±  % in the remaining phases. unclassified burkholderiales increased significantly from  ±  % in phases a through d to  ±  % in phases d and e. performance of trickling nitrification biofilm reactor based on the collated water quality data, effluent from the trickling nitrification biofilm reactor clustered into the six operational phases, which corresponded with the influent ammonium content. as the influent ammonium content was increased stepwise from 40 mg/l to 200 mg/l, the cod in the effluent also increased and was significantly higher in the effluent collected during phase d d and e compared to those collected during phase a and c. ammonium removal was optimally achieved at an average of 44% reduction per day during days through of operation, which corresponded with an influent concentration of 80–110 mg/l. during these phases, ammonium was primarily converted to nitrate with minimal nitrite accumulation in the effluent. subsequently, the trickling biofilm reactor experienced deterioration in its nitrification efficacy as influent ammonium content increased to 150 mg/l and 200 mg/l. ammonium was only partially converted to nitrate and an accumulation of nitrite was observed in the effluent collected during these phases. nitrite levels in the effluent during the two latter phases were significantly higher than during phases a through c.fig. principal component analysis of the measured water quality parameters in the effluent stream throughout the operational phases. the pca plot suggests that there was high ammonium nh4 + oxidation during phases b and c but an accumulation of nitrite and an increase in the chemical oxygen demand during the late operational phases d and e fig. correlative trends between effluent quality of the trickling nitrification biofilm reactor and the abundance ratio of predominant aob and nob. concentrations of ammonium, nitrite, nitrate and chemical oxygen demand in the effluent stream of the trickling nitrification biofilm reactor column throughout the different operational phases. the ratio of nse1472-, nso190- and nsv443-targeted aob against ntspa712-targeted nob, obtained using hope. the ratio of nitrosospira and nitrosomonas against nitrospira, obtained using amplicon sequencing correlative trends of aob and nob proportions with reactor performance the ratios in the relative abundance of these aob and nob targeted by ntspa hope primer were further evaluated and compared to the reactor performance based on its water quality parameters. it was determined that the average ratio of aob targeted by nse and nso normalized against ntspa712-targeted nob was and, respectively, during the phases that corresponded to stable reactor performance as evidenced by low nitrite content in the effluent. during phases e and d which corresponded with deterioration in the reactor performance, the ratio of aob targeted by nse and nso normalized against ntspa712-targeted nob increased to and, respectively. there was, however, a less significant increase in the ratio of aob targeted by nsv against nit1017-targeted nob. these ratios obtained using hope were compared against those generated from 16s rrna-based amplicon sequencing data. the ratio of the relative abundance of nitrosomonas compared against nitrospira as determined by amplicon sequencing was during the stable operational phases. this ratio markedly increased to when the reactor deteriorated in performance. likewise, the ratio of nitrosospira normalized against nitrospira was observed to be during phase a to c, and increased to in the latter phases. detection of aob and nob in wwtp samples using hope and amplicon sequencing both hope and amplicon sequencing revealed that the activated sludge from the local full-scale wwtp had lower relative abundances of nitrosomonas, nitrosospira and nitrobacter compared to ntspa712-targeted nitrospira. nitrate concentrations were consequentially higher in the effluent than in the influent and chlorinated effluent. hope and amplicon sequencing were further used to evaluate the nitrifying bacterial populations present in these wastewater samples. similar to the activated sludge, the wastewaters did not have any detectable nitrobacter but instead, had relatively higher abundances of nitrospira compared to nitrosomonas and nitrosospira. to illustrate, the relative abundance of ntspa712-targeted nob ranged from % in the influent to % in the chlorinated effluent. this relative abundance in nitrospira was at least times higher than the aob targeted by nse and nsv. amplicon sequencing also showed a relatively higher abundances of nitrospira compared to both aob in the three types of wastewater streams.table water quality for the influent, effluent and chlorinated effluent samples collected from the wastewater treatment plant and the corresponding relative abundances of aob and nob as quantified using the hope and amplicon sequencing approach nitrosomonas nitrosospira nitrobacter nitrospira relative abundance with respect to eub338ia discussion this study demonstrates the development of hope as a high-throughput method to quantify for the nitrifying bacterial populations including aob, nob and nitrospira-affiliated comammox. along with the aoa, which was not targeted by the hope approach in this study, these nitrifying populations play an important role in the activated sludge process to convert ammonia in the untreated wastewater to nitrite and nitrate. the nitrate is subsequently converted to dinitrogen and nitrogen by denitrifiers. collectively, the nitrifying and denitrifying populations reduce the total nitrogenous content in wastewater to a level permissible for discharge or reuse. it has previously been hypothesized that the theoretical aob/nob ratio in a nitrification process would be due to the differences in electron generation and biomass yield between aob and nob. simultaneous nitrification/denitrification can possibly increase this ratio to. however, the inverse may have a potentially detrimental impact on overall functionality of the nitrification system, especially since excess aob can lead to accumulation of the hydroxylamine intermediate which is inhibitory towards nob. as such, the hope method is developed with the intention to serve as a practically applicable tool for monitoring the aob and nob ratio and hence used to infer the efficiency of the nitrification process in wastewater treatment systems. for this purpose, the hope primers were designed to target nitrosomonas, nitrosospira, nitrobacter and nitrospira at the order, subcluster and/or genus level, and the relative abundance at each hierarchical taxonomical level was normalized against that of the total bacteria. primers were arranged into multiple reaction tubes that involve multiplexing of up to four targets per reaction, and tested against samples collected from a trickling nitrification biofilm reactor. it was observed from the hope data that the relative abundance of nitrosomonas increased with increasing nh4 + concentrations while abundances of nitrosospira did not. in contrast to the aob, the nob genus nitrospira was detrimentally affected by the shock loading of the reactor. during the stable phases b and c, nitrospira targeted by the ntspa572-primer was more predominant than nitrobacter. however, with the accumulation of nitrite in the late operational phases, nitrospira was depleted while nitrobacter, targeted by the nit1000-primer, remained relatively stable in its relative abundance against the total bacteria. this observation coincides with results reported by knapp and graham, in which nitrospira was observed to be more perturbed than nitrobacter during a destabilizing experiment. in addition, given that both ntspa and ntspa hope primers target at least one of the three recently identified comammox, the higher relative abundance of hope-targeted nitrospira than nitrobacter is likely to suggest that comammox is contributing to an equally important nitrification role in this reactor. similarly, the depletion in the collective relative abundance of ntspa572-targeted nitrospira suggests the susceptibility of comammox to shock loading events. the high nh4 + content exposure appeared to affect specific groups of aob and nob differently. it is generally thought that nitrospira-like bacteria are k-strategists with a high affinity for no2 − and oxygen, and reach high densities under substrate-limited conditions. this is in contrast to the nitrobacter species, which are r-strategists with a lower no2 − and oxygen affinity and outcompete nitrospira only at higher substrate concentrations. an additional postulation to account for the better survival of nitrobacter during shock loading events is its greater metabolic diversity under stressed conditions. this genus exhibits versatile metabolism and is able to grow either mixotrophically or chemoorganotrophically. although nitrobacter was able to withstand shock loading and increase in its relative abundance as quantified by hope, it was observed that the increase in its relative abundance did not correlate with any decrease in nitrite content of the effluent. furthermore, the increase in nitrobacter did not correlate with a subsequent conversion of nitrite to nitrate in the effluent. this suggests that nitrobacter may not be playing as important of a role as nitrospira and comammox in nitrite oxidation in the wastewater treatment process. given that both nitrosomonas and nitrospira seem to correlate with nitrification functionality, the hope-obtained relative abundances of both groups were further expressed as a proportional ratio. the ratio of nitrosomonas targeted by either nse or nso hope primer against ntspa572-targeted nitrospira ranged from to in the nitrification biofilm reactor during the phases that corresponded to stable reactor performance as evidenced by the low nitrite content in the effluent. however, the exceedingly high ammonium content of the later phases of operation resulted in a toxic shock response by the nitrifying bacterial populations, which correlated with an observed increase in the aob/nob ratio. this ratio increased to a range of to when the bioreactor performance deteriorated. based on energetic calculations, it was proposed that a theoretical aob/nob ratio of should be obtained during a stable and functional nitrification process. this theoretical ratio was in agreement with the observation made in this study, and reiterates the need to maintain an optimal proportion of aob to nob in a well-performing nitrification process. the hope method had been previously validated against other conventional methods such as qpcr and fish. in this study, to validate the data obtained from hope, the relative abundances of nitrifying bacterial populations in the trickling biofilter and in samples collected from a full-scale wwtp were further validated against results from amplicon sequencing. the observed difference in the abundance values between hope and amplicon sequencing is likely due to differences in the primer coverage for the individual aob and nob groups. for example, nso targets % of the nitrosomonas and % of the nitrosospira sequences collated in the rdp classifier database. in contrast, the universal primers 515f and 907r used for amplicon sequencing target a higher percentage of nitrosomonas and nitrosospira. furthermore, the relative abundances reported by hope were expressed by normalizing against total bacteria as detected by universal primers that differed from those used in amplicon sequencing. although the relative abundances determined by both methods varied in terms of the absolute values, both datasets showed good spearman’s rank correlation. furthermore, multivariate nmds plots from both hope and amplicon sequencing also showed good correlation at a significant confidence level. this indicates that the observed trends in the nitrifying bacterial population dynamics revealed by both methods were well aligned. although similar conclusions were obtained from both hope and amplicon sequencing with regards to the nitrifying population dynamics in the presence of shock loading, a limitation of the hope approach is that it can only detect bacterial targeted by the designed primers. for instance, through the use of amplicon sequencing and not hope, it was determined that the detrimental effect of the shock loading event on nitrospira can possibly be due to its competition for substrates with heterotrophic bacteria. such heterotrophic bacteria included unclassified burkholderiales and unclassified xanthomonadaceae, which increased in relative abundance throughout all operational phase. furthermore, the hope technique developed in this study does not contain primers that target the aoa. although aoa had been previously found abundant in wastewater treatment systems, both genus nitrososphaera and nitrosopumilus were found to be in low abundance in the samples when evaluated by high-throughput amplicon sequencing. further optimization on the hope technique would be required to design a range of aoa-targeting primers, as well as to improve its current detection limits of % so as to facilitate the quantification of potentially low abundance species. regardless, an advantage of the hope approach is that new primers, for example those designed to target the different sublineages of nitrospira and nitrosomonas that were selected for by varying nitrite and ammonium concentrations, can be easily added to the multiplexing tube arrangement when these primers become available. the additional depth offered by high-throughput amplicon sequencing would be useful as a tool for periodically benchmarking hope-based results obtained on a more frequent basis. although this study did not attempt to challenge the limits of high-throughput multiplexing capability of hope, the method can in theory perform up to 32-plexes per reaction tube, and modifications to the hope approach can be efficiently made to facilitate the simultaneous detection of the expanding list of nitrifying bacteria. furthermore, this study demonstrated a good correlation between the results obtained from the hope approach and high-throughput sequencing. thus, depending on the intended purpose of an experimental or practical study, hope may be particularly useful for the frequent tracking of the occurrence and abundance of selected microbial populations without the need to perform high-throughput amplicon sequencing. this in turn reduces the need to sieve through a full range of microbial population data for every sampling event.  in summary, this study demonstrates the applicability and adaptability of hope for assessing abundances of predominant aob and nob groups that have been identified as integral to wastewater treatment systems. this method allows for the simultaneous monitoring of relative abundances of aob, nob and comammox groups, which can be used to provide indicative data of nitrification performance and efficiency. coupled with previous application of hope to a variety of target bacteria and sample types, this study demonstrates the versatility and applicability of hope for a wide range of microbial ecology studies. additional file additional file 1: target coverage for the different primers used in this study. abbreviations aoaammonia-oxidizing archaea aobammonia-oxidizing bacteria codchemical oxygen demand ddntpdideoxynucleoside triphosphate dpcrdigital polymerase chain reaction fishfluorescent in-situ hybridization hopehierarchical oligonucleotide primer extension mbrmembrane bioreactor nh4+-nammonium nmdsnon-metric threshold multidimensional scaling no2--nnitrite no3--nnitrate nobnitrite-oxidizing bacteria ntnucleotides pbsphosphate buffer saline pcaprincipal component analysis qpcrquantitative polymerase chain reaction rdpribosomal database project rsaprecombinant shrimp alkaline phosphatase tntotal nitrogen wwtpwastewater treatment plant 
the method was tested against biomass obtained from a laboratory-scale biofilm-based trickling reactor, and the findings were validated against those obtained by 16s rrna gene-based amplicon sequencing. overall, the hope method was developed and validated against 16s rrna gene-based amplicon sequencing for the purpose of performing simultaneous monitoring of relative abundance of nitrifying populations. single nucleotide primer extensionquantitative monitoring16s rrna gene-based amplicon sequencingaob/nob ratioshock loading eventhttp://dx.doi.org/ /501100004052king abdullah university of science and technologyfcc/1/1971-06-01issue-copyright-statement© the author 2017 quantitative measurements of these nitrifying populations obtained via hope would be indicative of reactor performance and nitrification functionality. hierarchical oligonucleotide primer extension, previously developed to rapidly quantify relative abundances of specific microbial groups of interest, was applied in this study to track the abundances of the important nitrifying bacterial populations. establishing an optimal proportion of nitrifying microbial populations, including ammonia-oxidizing bacteria, nitrite-oxidizing bacteria, complete nitrite oxidizers and ammonia-oxidizing archaea, is important for ensuring the efficiency of nitrification in water treatment systems. based on the relative abundance of aob and nob obtained from hope, it was determined that the optimal ratio of aob against nob ranged from to during stable reactor performance. in contrast, nitrosospira remained stable in its relative abundance against the total community throughout the operational phases. our findings indicated a good correlation between the relative abundance of nitrifying bacterial populations obtained using both hope and amplicon sequencing. hope showed a significant increase in the relative abundance of aob, specifically nitrosomonas, with increasing ammonium content and shock loading. there was a corresponding significant decrease in the relative abundance of nob, specifically nitrospira and those affiliated to comammox, during the shock loading.
 establishing an optimal proportion of nitrifying microbial populations, including ammonia-oxidizing bacteria, nitrite-oxidizing bacteria, complete nitrite oxidizers and ammonia-oxidizing archaea, is important for ensuring the efficiency of nitrification in water treatment systems. hierarchical oligonucleotide primer extension, previously developed to rapidly quantify relative abundances of specific microbial groups of interest, was applied in this study to track the abundances of the important nitrifying bacterial populations.  the method was tested against biomass obtained from a laboratory-scale biofilm-based trickling reactor, and the findings were validated against those obtained by 16s rrna gene-based amplicon sequencing. our findings indicated a good correlation between the relative abundance of nitrifying bacterial populations obtained using both hope and amplicon sequencing. hope showed a significant increase in the relative abundance of aob, specifically nitrosomonas, with increasing ammonium content and shock loading. in contrast, nitrosospira remained stable in its relative abundance against the total community throughout the operational phases. there was a corresponding significant decrease in the relative abundance of nob, specifically nitrospira and those affiliated to comammox, during the shock loading. based on the relative abundance of aob and nob obtained from hope, it was determined that the optimal ratio of aob against nob ranged from to during stable reactor performance.  overall, the hope method was developed and validated against 16s rrna gene-based amplicon sequencing for the purpose of performing simultaneous monitoring of relative abundance of nitrifying populations. quantitative measurements of these nitrifying populations obtained via hope would be indicative of reactor performance and nitrification functionality. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords single nucleotide primer extensionquantitative monitoring16s rrna gene-based amplicon sequencingaob/nob ratioshock loading eventhttp://dx.doi.org/ /501100004052king abdullah university of science and technologyfcc/1/1971-06-01issue-copyright-statement© the author 2017 