Global genomic instability caused by reduced expression of DNA polymerase ε in yeast.

SignificanceAlthough most studies of the genetic regulation of genome stability involve an analysis of mutations within the coding sequences of genes required for DNA replication or DNA repair, recent studies in yeast show that reduced levels of wild-type enzymes can also produce a mutator phenotype. By whole-genome sequencing and other methods, we find that reduced levels of the wild-type DNA polymerase ε in yeast greatly increase the rates of mitotic recombination, aneuploidy, and single-base mutations. The observed pattern of genome instability is different from those observed in yeast strains with reduced levels of the other replicative DNA polymerases, Pol α and Pol δ. These observations are relevant to our understanding of cancer and other diseases associated with genetic instability.

Genome-wide mapping of spontaneous genetic alterations in diploid yeast cells.

Genomic alterations including single-base mutations, deletions and duplications, translocations, mitotic recombination events, and chromosome aneuploidy generate genetic diversity. We examined the rates of all of these genetic changes in a diploid strain of Saccharomyces cerevisiae by whole-genome sequencing of many independent isolates (n = 93) subcloned about 100 times in unstressed growth conditions. The most common alterations were point mutations and small (<100 bp) insertion/deletions (n = 1,337) and mitotic recombination events (n = 1,215). The diploid cells of most eukaryotes are heterozygous for many single-nucleotide polymorphisms (SNPs). During mitotic cell divisions, recombination can produce derivatives of these cells that have become homozygous for the polymorphisms, termed loss-of-heterozygosity (LOH) events. LOH events can change the phenotype of the cells and contribute to tumor formation in humans. We observed two types of LOH events: interstitial events (conversions) resulting in a short LOH tract (usually less than 15 kb) and terminal events (mostly cross-overs) in which the LOH tract extends to the end of the chromosome. These two types of LOH events had different distributions, suggesting that they may have initiated by different mechanisms. Based on our results, we present a method of calculating the probability of an LOH event for individual SNPs located throughout the genome. We also identified several hotspots for chromosomal rearrangements (large deletions and duplications). Our results provide insights into the relative importance of different types of genetic alterations produced during vegetative growth.

Hydrogen peroxide, a potent inducer of global genomic instability.

Oxidative stress has been implicated in a variety of human diseases. One plausible mechanism is that reactive active species can induce DNA damages and jeopardize genome integrity. To explore how oxidative stress results in global genomic instability in cells, our current study examined the genomic alterations caused by H2O2 exposure at the whole genome level in yeast. Using SNP microarrays and genome sequencing, we mapped H2O2-induced genomic alterations in the yeast genome ranging from point mutations and mitotic recombination to chromosomal aneuploidy. Our results suggested most H2O2-induced mitotic recombination events were the result of DNA double-stand breaks generated by hydroxyl radicals. Moreover, the mutagenic effect of H2O2 was shown to be largely dependent on DNA polymerase ζ. Lastly, we showed that H2O2 exposure allows rapid phenotypic evolution in yeast strains. Our findings indicate DNA lesions resulting from H2O2 may be general factors that drive genome instability and phenotypic evolution in organisms.

Global Analysis of Furfural-Induced Genomic Instability Using a Yeast Model.

Furfural is an important renewable precursor for multiple commercial chemicals and fuels; a main inhibitor existing in cellulosic hydrolysate, which is used for bioethanol fermentation; and a potential carcinogen, as well. Using a genetic system in Saccharomyces cerevisiae that allows detection of crossover events, we observed that the frequency of mitotic recombination was elevated by 1.5- to 40-fold when cells were treated with 0.1 g/liter to 20 g/liter furfural. Analysis of the gene conversion tracts associated with crossover events suggested that most furfural-induced recombination resulted from repair of DNA double-strand breaks (DSBs) that occurred in the G1 phase. Furfural was incapable of breaking DNA directly in vitro but could trigger DSBs in vivo related to reactive oxygen species accumulation. By whole-genome single nucleotide polymorphism (SNP) microarray and sequencing, furfural-induced genomic alterations that range from single base substitutions, loss of heterozygosity, and chromosomal rearrangements to aneuploidy were explored. At the whole-genome level, furfural-induced events were evenly distributed across 16 chromosomes but were enriched in high-GC-content regions. Point mutations, particularly the C-to-T/G-to-A transitions, were significantly elevated in furfural-treated cells compared to wild-type cells. This study provided multiple novel insights into the global effects of furfural on genomic stability.IMPORTANCE Whether and how furfural affects genome integrity have not been clarified. Using a Saccharomyces cerevisiae model, we found that furfural exposure leads to in vivo DSBs and elevation in mitotic recombination by orders of magnitude. Gross chromosomal rearrangements and aneuploidy events also occurred at a higher frequency in furfural-treated cells. In a genome-wide analysis, we show that the patterns of mitotic recombination and point mutations differed dramatically in furfural-treated cells and wild-type cells.

Global analysis of genomic instability caused by DNA replication stress in Saccharomyces cerevisiae.

DNA replication stress (DRS)-induced genomic instability is an important factor driving cancer development. To understand the mechanisms of DRS-associated genomic instability, we measured the rates of genomic alterations throughout the genome in a yeast strain with lowered expression of the replicative DNA polymerase δ. By a genetic test, we showed that most recombinogenic DNA lesions were introduced during S or G2 phase, presumably as a consequence of broken replication forks. We observed a high rate of chromosome loss, likely reflecting a reduced capacity of the low-polymerase strains to repair double-stranded DNA breaks (DSBs). We also observed a high frequency of deletion events within tandemly repeated genes such as the ribosomal RNA genes. By whole-genome sequencing, we found that low levels of DNA polymerase δ elevated mutation rates, both single-base mutations and small insertions/deletions. Finally, we showed that cells with low levels of DNA polymerase δ tended to accumulate small promoter mutations that increased the expression of this polymerase. These deletions conferred a selective growth advantage to cells, demonstrating that DRS can be one factor driving phenotypic evolution.

Genetic characterization and modification of a bioethanol-producing yeast strain.

Yeast Saccharomyces cerevisiae strains isolated from different sources generally show extensive genetic and phenotypic diversity. Understanding how genomic variations influence phenotypes is important for developing strategies with improved economic traits. The diploid S. cerevisiae strain NY1308 is used for cellulosic bioethanol production. Whole genome sequencing identified an extensive amount of single nucleotide variations and small insertions/deletions in the genome of NY1308 compared with the S288c genome. Gene annotation of the assembled NY1308 genome showed that 43 unique genes are absent in the S288c genome. Phylogenetic analysis suggested most of the unique genes were obtained through horizontal gene transfer from other species. RNA-Seq revealed that some unique genes were not functional in NY1308 due to unidentified intron sequences. During bioethanol fermentation, NY1308 tends to flocculate when certain inhibitors (derived from the pretreatment of cellulosic feedstock) are present in the fermentation medium. qRT-PCR and genetic manipulation confirmed that the novel gene, NYn43, contributed to the flocculation ability of NY1308. Deletion of NYn43 resulted in a faster fermentation rate for NY1308. This work disclosed the genetic characterization of a bioethanol-producing S. cerevisiae strain and provided a useful paradigm showing how the genetic diversity of the yeast population would facilitate the personalized development of desirable traits.

Effects of genome duplication on phenotypes and industrial applications of Saccharomyces cerevisiae strains.

Polyploidy is common in Saccharomyces cerevisiae strains, but the physiological and phenotypic effects of ploidy changes have not been fully clarified. Here, isogenic diploid, triploid, and tetraploid S. cerevisiae strains were constructed from a haploid strain, CEN.PK2-1C. Stress tolerance and ethanol fermentation performance of the four euploid strains were compared. Each euploid strain had strengths and weaknesses in tolerance to certain stressors, and no single strain was tolerant of all stressors. The diploid had higher ethanol production than the other strains in normal fermentation medium, while the triploid strain showed the fastest fermentation rate in the presence of inhibitors found in lignocellulosic hydrolysate. Physiological determination revealed diverse physiological attributes, such as trehalose, ergosterol, glutathione, and anti-oxidative enzymes among the strains. Our analyses suggest that both ploidy parity and number of chromosome sets contribute to changes in physiological status. Using qRT-PCR, different expression patterns of genes involved in the regulation of cell morphology and the biosynthesis of key physiological attributes among strains were determined. Our data provide novel insights into the multiple effects of genome duplication on yeast cells and are a useful reference for breeding excellent strains used in specific industrial applications.

Genomic structural variation contributes to phenotypic change of industrial bioethanol yeast Saccharomyces cerevisiae.

Genomic structural variation (GSV) is a ubiquitous phenomenon observed in the genomes of Saccharomyces cerevisiae strains with different genetic backgrounds; however, the physiological and phenotypic effects of GSV are not well understood. Here, we first revealed the genetic characteristics of a widely used industrial S. cerevisiae strain, ZTW1, by whole genome sequencing. ZTW1 was identified as an aneuploidy strain and a large-scale GSV was observed in the ZTW1 genome compared with the genome of a diploid strain YJS329. These GSV events led to copy number variations (CNVs) in many chromosomal segments as well as one whole chromosome in the ZTW1 genome. Changes in the DNA dosage of certain functional genes directly affected their expression levels and the resultant ZTW1 phenotypes. Moreover, CNVs of large chromosomal regions triggered an aneuploidy stress in ZTW1. This stress decreased the proliferation ability and tolerance of ZTW1 to various stresses, while aneuploidy response stress may also provide some benefits to the fermentation performance of the yeast, including increased fermentation rates and decreased byproduct generation. This work reveals genomic characters of the bioethanol S. cerevisiae strain ZTW1 and suggests that GSV is an important kind of mutation that changes the traits of industrial S. cerevisiae strains.

Two UDP-glucuronic acid decarboxylases involved in the biosynthesis of a bacterial exopolysaccharide in Paenibacillus elgii.

Xylose is described as a component of bacterial exopolysaccharides in only a limited number of bacterial strains. A bacterial strain, Paenibacillus elgii, B69 was shown to be efficient in producing a xylose-containing exopolysaccharide. Sequence analysis was performed to identify the genes encoding the uridine diphosphate (UDP)-glucuronic acid decarboxylase required for the synthesis of UDP-xylose, the precursor of the exopolysaccharide. Two sequences, designated as Peuxs1 and Peuxs2, were found as the candidate genes for such enzymes. The activities of the UDP-glucuronic acid decarboxylases were proven by heterologous expression and real-time nuclear magnetic resonance analysis. The intracellular activity and effect of these genes on the synthesis of exopolysaccharide were further investigated by developing a thymidylate synthase based knockout system. This system was used to substitute the conventional antibiotic resistance gene system in P. elgii, a natural multi-antibiotic resistant strain. Results of intracellular nucleotide sugar analysis showed that the intracellular UDP-xylose and UDP-glucuronic acid levels were affected in Peuxs1 or Peuxs2 knockout strains. The knockout of either Peuxs1 or Peuxs2 reduced the polysaccharide production and changed the monosaccharide ratio. No polysaccharide was found in the Peuxs1/Peuxs2 double knockout strain. Our results show that P. elgii can be efficient in forming UDP-xylose, which is then used for the synthesis of xylose-containing exopolysaccharide.

Genomic structural variations contribute to trait improvement during whole-genome shuffling of yeast.

Whole-genome shuffling (WGS) is a powerful technology of improving the complex traits of many microorganisms. However, the molecular mechanisms underlying the altered phenotypes in isolates were less clarified. Isolates with significantly enhanced stress tolerance and ethanol titer under very-high-gravity conditions were obtained after WGS of the bioethanol Saccharomyces cerevisiae strain ZTW1. Karyotype analysis and RT-qPCR showed that chromosomal rearrangement occurred frequently in genome shuffling. Thus, the phenotypic effects of genomic structural variations were determined in this study. RNA-Seq and physiological analyses revealed the diverse transcription pattern and physiological status of the isolate S3-110 and ZTW1. Our observations suggest that the improved stress tolerance of S3-110 can be largely attributed to the copy number variations in large DNA regions, which would adjust the ploidy of yeast cells and expression levels of certain genes involved in stress response. Overall, this work not only constructed shuffled S. cerevisiae strains that have potential industrial applications but also provided novel insights into the molecular mechanisms of WGS and enhanced our knowledge on this useful breeding strategy.

Comparative functional genomics to reveal the molecular basis of phenotypic diversities and guide the genetic breeding of industrial yeast strains.

An understanding of the genetic basis underlying the phenotypic variations of yeast strains would guide the breeding of this useful microorganism. Here, comparative functional genomics (CFG) of two bioethanol Saccharomyces cerevisiae strains (YJS329 and ZK2) with different stress tolerances and ethanol fermentation performances were performed. Our analysis indicated that different patterns of gene expression in the central carbon metabolism, antioxidative factors, and membrane compositions of these two strains are the main contributors to their various traits. Some of the differently expressed genes were directly caused by the genomic structural variations between YJS329 and ZK2. Moreover, CFG of these two strains also led to novel insights into the mechanism of stress tolerance in yeast. For example, it was found that more oleic acid in the plasma membrane contributes to the acetic acid tolerance of yeast. Based on the genetic information particular to each strain, strategies to improve their adaptability and ethanol fermentation performances were designed and confirmed. Thus, CFG could not only help reveal basis of phenotypic diversities but also guide the genetic breeding of industrial microorganisms.

Genome sequencing and genetic breeding of a bioethanol Saccharomyces cerevisiae strain YJS329.

BACKGROUND: Environmental stresses and inhibitors encountered by Saccharomyces cerevisiae strains are the main limiting factors in bioethanol fermentation. Strains with different genetic backgrounds usually show diverse stress tolerance responses. An understanding of the mechanisms underlying these phenotypic diversities within S. cerevisiae populations could guide the construction of strains with desired traits. RESULTS: We explored the genetic characteristics of the bioethanol S. cerevisiae strain YJS329 and elucidated how genetic variations in its genome were correlated with specified traits compared to similar traits in the S288c-derived strain, BYZ1. Karyotypic electrophoresis combined with array-comparative genomic hybridization indicated that YJS329 was a diploid strain with a relatively constant genome as a result of the fewer Ty elements and lack of structural polymorphisms between homologous chromosomes that it contained. By comparing the sequence with the S288c genome, a total of 64,998 SNPs, 7,093 indels and 11 unique genes were identified in the genome of YJS329-derived haploid strain YJSH1 through whole-genome sequencing. Transcription comparison using RNA-Seq identified which of the differentially expressed genes were the main contributors to the phenotypic differences between YJS329 and BYZ1. By combining the results obtained from the genome sequences and the transcriptions, we predicted how the SNPs, indels and chromosomal copy number variations may affect the mRNA expression profiles and phenotypes of the yeast strains. Furthermore, some genetic breeding strategies to improve the adaptabilities of YJS329 were designed and experimentally verified. CONCLUSIONS: Through comparative functional genomic analysis, we have provided some insights into the mechanisms underlying the specific traits of the bioenthanol strain YJS329. The work reported here has not only enriched the available genetic resources of yeast but has also indicated how functional genomic studies can be used to improve genetic breeding in yeast.

Battacin (Octapeptin B5), a new cyclic lipopeptide antibiotic from Paenibacillus tianmuensis active against multidrug-resistant Gram-negative bacteria.

Hospital-acquired infections caused by drug-resistant bacteria are a significant challenge to patient safety. Numerous clinical isolates resistant to almost all commercially available antibiotics have emerged. Thus, novel antimicrobial agents, specifically those for multidrug-resistant Gram-negative bacteria, are urgently needed. In the current study, we report the isolation, structure elucidation, and preliminary biological characterization of a new cationic lipopeptide antibiotic, battacin or octapeptin B5, produced from a Paenibacillus tianmuensis soil isolate. Battacin kills bacteria in vitro and has potent activity against Gram-negative bacteria, including multidrug-resistant and extremely drug-resistant clinical isolates. Hospital strains of Escherichia coli and Pseudomonas aeruginosa are the pathogens most sensitive to battacin, with MICs of 2 to 4 µg/ml. The ability of battacin to disrupt the outer membrane of Gram-negative bacteria is comparable to that of polymyxin B, the last-line therapy for infections caused by antibiotic-resistant Gram-negative bacteria. However, the capacity of battacin to permeate bacterial plasma membranes is less extensive than that of polymyxin B. The bactericidal kinetics of battacin correlate with the depolarization of the cell membrane, suggesting that battacin kills bacteria by disrupting the cytoplasmic membrane. Other studies indicate that battacin is less acutely toxic than polymyxin B and has potent in vivo biological activity against E. coli. Based on the findings of the current study, battacin may be considered a potential therapeutic agent for the treatment of infections caused by antibiotic-resistant Gram-negative bacteria.

Identification and functional analysis of gene cluster involvement in biosynthesis of the cyclic lipopeptide antibiotic pelgipeptin produced by Paenibacillus elgii.

BACKGROUND: Pelgipeptin, a potent antibacterial and antifungal agent, is a non-ribosomally synthesised lipopeptide antibiotic. This compound consists of a ß-hydroxy fatty acid and nine amino acids. To date, there is no information about its biosynthetic pathway. RESULTS: A potential pelgipeptin synthetase gene cluster (plp) was identified from Paenibacillus elgii B69 through genome analysis. The gene cluster spans 40.8 kb with eight open reading frames. Among the genes in this cluster, three large genes, plpD, plpE, and plpF, were shown to encode non-ribosomal peptide synthetases (NRPSs), with one, seven, and one module(s), respectively. Bioinformatic analysis of the substrate specificity of all nine adenylation domains indicated that the sequence of the NRPS modules is well collinear with the order of amino acids in pelgipeptin. Additional biochemical analysis of four recombinant adenylation domains (PlpD A1, PlpE A1, PlpE A3, and PlpF A1) provided further evidence that the plp gene cluster involved in pelgipeptin biosynthesis. CONCLUSIONS: In this study, a gene cluster (plp) responsible for the biosynthesis of pelgipeptin was identified from the genome sequence of Paenibacillus elgii B69. The identification of the plp gene cluster provides an opportunity to develop novel lipopeptide antibiotics by genetic engineering.

Drug resistance marker-aided genome shuffling to improve acetic acid tolerance in Saccharomyces cerevisiae.

Acetic acid existing in a culture medium is one of the most limiting constraints in yeast growth and viability during ethanol fermentation. To improve acetic acid tolerance in Saccharomyces cerevisiae strains, a drug resistance marker-aided genome shuffling approach with higher screen efficiency of shuffled mutants was developed in this work. Through two rounds of genome shuffling of ultraviolet mutants derived from the original strain 308, we obtained a shuffled strain YZ2, which shows significantly faster growth and higher cell viability under acetic acid stress. Ethanol production of YZ2 (within 60 h) was 21.6% higher than that of 308 when 0.5% (v/v) acetic acid was added to fermentation medium. Membrane integrity, higher in vivo activity of the H+-ATPase, and lower oxidative damage after acetic acid treatment are the possible reasons for the acetic acid-tolerance phenotype of YZ2. These results indicated that this novel genome shuffling approach is powerful to rapidly improve the complex traits of industrial yeast strains.

The coiled-coil neck domain of human pulmonary surfactant protein D drives trimerization and stabilization of thioredoxin, a heterologous non-collagenous protein.

The coiled-coil neck domain of pulmonary surfactant protein D (SP-D) is required for trimeric association and the subsequent assembly of functional dodecamers of SP-D. It is also necessary and sufficient for trimerization of a heterologous collagen sequence. To investigate whether it is capable of driving trimerization of heterologous non-collagenous proteins, we expressed and purified a fusion of a heterologous non-collagenous sequence (thioredoxin) to the coiled-coil neck domain of human SP-D here. While western blot analysis detected a small population of stable trimers of the fusion protein, chemical cross-linking and SEC-HPLC indicated that the fusion protein was predominantly a trimer. In contrast, purified thioredoxin without the fusion was found only as monomers and dimers. We also measured the thermal stabilities (with circular dichroism) and degradation rates of these two proteins. Our data showed that the fusion protein had a melting temperature that was 13 K higher than that of thioredoxin and a longer degradation half life than thioredoxin. Our findings indicate that the coiled-coil neck domain of SP-D enables the trimerization and stabilization of the heterologous non-collagenous thioredoxin. It may provide new clues for further study on the application of this human original coiled-coil domain in protein engineering to construct trimeric functional fusion proteins.

Isolation and identification of newly isolated antagonistic Streptomyces sp. strain AP19-2 producing chromomycins.

A new antagonistic strain of actinomycete, designated AP19-2, was isolated from the feces of giant pandas inhabiting the Foping National Nature Reserve in China. Cultural characteristic studies strongly suggested that this strain is a member of the genus Streptomyces. The nucleotide sequence of the 16S rRNA gene of strain AP19-2 evidenced profound similarity (97-99%) with other Streptomyces strains. Two pure active molecules were isolated from a fermentation broth of Streptomyces sp. strain AP19-2 via extraction, concentration, silica gel G column chromatography, and HPLC. The chemical structures of the two related compounds (referred to as chromomycin A2 and chromomycin A3) were established on the basis of their Infrared spectra (IR), High Resolution Electrospray Ionization Mass Spectrometry (HR-ESI-MS), and Nuclear Magnetic Resonance (NMR) data, and by comparison with published data.

Effects of Temperature on the Meiotic Recombination Landscape of the Yeast Saccharomyces cerevisiae.

Although meiosis in warm-blooded organisms takes place in a narrow temperature range, meiosis in many organisms occurs over a wide variety of temperatures. We analyzed the properties of meiosis in the yeast Saccharomyces cerevisiae in cells sporulated at 14°C, 30°C, or 37°C. Using comparative-genomic-hybridization microarrays, we examined the distribution of Spo11-generated meiosis-specific double-stranded DNA breaks throughout the genome. Although there were between 300 and 400 regions of the genome with high levels of recombination (hot spots) observed at each temperature, only about 20% of these hot spots were found to have occurred independently of the temperature. In S. cerevisiae, regions near the telomeres and centromeres tend to have low levels of meiotic recombination. This tendency was observed in cells sporulated at 14°C and 30°C, but not at 37°C. Thus, the temperature of sporulation in yeast affects some global property of chromosome structure relevant to meiotic recombination. Using single-nucleotide polymorphism (SNP)-specific whole-genome microarrays, we also examined crossovers and their associated gene conversion events as well as gene conversion events that were unassociated with crossovers in all four spores of tetrads obtained by sporulation of diploids at 14°C, 30°C, or 37°C. Although tetrads from cells sporulated at 30°C had slightly (20%) more crossovers than those derived from cells sporulated at the other two temperatures, spore viability was good at all three temperatures. Thus, despite temperature-induced variation in the genetic maps, yeast cells produce viable haploid products at a wide variety of sporulation temperatures.IMPORTANCE In the yeast Saccharomyces cerevisiae, recombination is usually studied in cells that undergo meiosis at 25°C or 30°C. In a genome-wide analysis, we showed that the locations of genomic regions with high and low levels of meiotic recombination (hot spots and cold spots, respectively) differed dramatically in cells sporulated at 14°C, 30°C, and 37°C. Thus, in yeast, and likely in other non-warm-blooded organisms, genetic maps are strongly affected by the environment.

Novel strategy to improve vanillin tolerance and ethanol fermentation performances of Saccharomycere cerevisiae strains.

The aim of this work was to develop a novel strategy for improving the vanillin tolerance and ethanol fermentation performances of Saccharomyces cerevisiae strains. Isogeneic diploid, triploid, and tetraploid S. cerevisiae strains were generated by genome duplication of haploid strain CEN.PK2-1C. Ploidy increments improved vanillin tolerance and diminished proliferation capability. Antimitotic drug methyl benzimidazol-2-ylcarbamate (MBC) was used to introduce chromosomal aberrations into the tetraploid S. cerevisiae strain. Interestingly, aneuploid mutants with DNA contents between triploid and tetraploid were more resistant to vanillin and showed faster ethanol fermentation rates than all euploid strains. The physiological characteristics of these mutants suggest that higher bioconversion capacities of vanillin and ergosterol contents might contribute to improved vanillin tolerance. This study demonstrates that genome duplication and MBC treatment is a powerful strategy to improve the vanillin tolerance of yeast strains.

Evaluation of nutrient limitation in aquatic ecosystems with nitrogen fixing bacteria.

There has always been a great need for simple and accurate bioassays for evaluating nutrient limitation in aquatic ecosystems. Whereas organic carbon is usually considered to be the limiting nutrient for microbial growth in many aquatic ecosystems, there are, however, many water sources that are limited by phosphorus or nitrogen. A method named "nitrogen fixing bacterial growth potential" (NFBGP) test, which is based on pre-culturing of autochthonous (target) microorganisms was described. The method was applied to evaluate phosphorus or nitrogen nutrient limitation in lake and sewage water samples using an isolate of the nitrogen fixing bacterium, Azorhizobium sp. WS6. The results corresponded well to those from the traditional algal growth potential (AGP) test and the bacterial regrowth potential (BRP) test, suggesting that the NFBGP test is a useful supplementary method for evaluating the limiting nutrient, especially phosphorus, in an aquatic environment.