Assessing the reliability of point mutation as data augmentation for deep learning with genomic data.

BACKGROUND: Deep neural networks (DNNs) have the potential to revolutionize our understanding and treatment of genetic diseases. An inherent limitation of deep neural networks, however, is their high demand for data during training. To overcome this challenge, other fields, such as computer vision, use various data augmentation techniques to artificially increase the available training data for DNNs. Unfortunately, most data augmentation techniques used in other domains do not transfer well to genomic data. RESULTS: Most genomic data possesses peculiar properties and data augmentations may significantly alter the intrinsic properties of the data. In this work, we propose a novel data augmentation technique for genomic data inspired by biology: point mutations. By employing point mutations as substitutes for codons, we demonstrate that our newly proposed data augmentation technique enhances the performance of DNNs across various genomic tasks that involve coding regions, such as translation initiation and splice site detection. CONCLUSION: Silent and missense mutations are found to positively influence effectiveness, while nonsense mutations and random mutations in non-coding regions generally lead to degradation. Overall, point mutation-based augmentations in genomic datasets present valuable opportunities for improving the accuracy and reliability of predictive models for DNA sequences.

Mutate and observe: utilizing deep neural networks to investigate the impact of mutations on translation initiation.

MOTIVATION: The primary regulatory step for protein synthesis is translation initiation, which makes it one of the fundamental steps in the central dogma of molecular biology. In recent years, a number of approaches relying on deep neural networks (DNNs) have demonstrated superb results for predicting translation initiation sites. These state-of-the art results indicate that DNNs are indeed capable of learning complex features that are relevant to the process of translation. Unfortunately, most of those research efforts that employ DNNs only provide shallow insights into the decision-making processes of the trained models and lack highly sought-after novel biologically relevant observations. RESULTS: By improving upon the state-of-the-art DNNs and large-scale human genomic datasets in the area of translation initiation, we propose an innovative computational methodology to get neural networks to explain what was learned from data. Our methodology, which relies on in silico point mutations, reveals that DNNs trained for translation initiation site detection correctly identify well-established biological signals relevant to translation, including (i) the importance of the Kozak sequence, (ii) the damaging consequences of ATG mutations in the 5'-untranslated region, (iii) the detrimental effect of premature stop codons in the coding region, and (iv) the relative insignificance of cytosine mutations for translation. Furthermore, we delve deeper into the Beta-globin gene and investigate various mutations that lead to the Beta thalassemia disorder. Finally, we conclude our work by laying out a number of novel observations regarding mutations and translation initiation. AVAILABILITY AND IMPLEMENTATION: For data, models, and code, visit github.com/utkuozbulak/mutate-and-observe.

Unraveling the MAX2 Protein Network in Arabidopsis thaliana: Identification of the Protein Phosphatase PAPP5 as a Novel MAX2 Interactor.

The F-box protein MORE AXILLARY GROWTH 2 (MAX2) is a central component in the signaling cascade of strigolactones (SLs) as well as of the smoke-derived karrikins (KARs) and the so far unknown endogenous KAI2 ligand (KL). The two groups of molecules are involved in overlapping and unique developmental processes, and signal-specific outcomes are attributed to perception by the paralogous α/ß-hydrolases DWARF14 (D14) for SL and KARRIKIN INSENSITIVE 2/HYPOSENSITIVE TO LIGHT (KAI2/HTL) for KAR/KL. In addition, depending on which receptor is activated, specific members of the SUPPRESSOR OF MAX2 1 (SMAX1)-LIKE (SMXL) family control KAR/KL and SL responses. As proteins that function in the same signal transduction pathway often occur in large protein complexes, we aimed at discovering new players of the MAX2, D14, and KAI2 protein network by tandem affinity purification in Arabidopsis cell cultures. When using MAX2 as a bait, various proteins were copurified, among which were general components of the Skp1-Cullin-F-box complex and members of the CONSTITUTIVE PHOTOMORPHOGENIC 9 signalosome. Here, we report the identification of a novel interactor of MAX2, a type 5 serine/threonine protein phosphatase, designated PHYTOCHROME-ASSOCIATED PROTEIN PHOSPHATASE 5 (PAPP5). Quantitative affinity purification pointed at PAPP5 as being more present in KAI2 rather than in D14 protein complexes. In agreement, mutant analysis suggests that PAPP5 modulates KAR/KL-dependent seed germination under suboptimal conditions and seedling development. In addition, a phosphopeptide enrichment experiment revealed that PAPP5 might dephosphorylate MAX2 in vivo independently of the synthetic SL analog, rac-GR24. Together, by analyzing the protein complexes to which MAX2, D14, and KAI2 belong, we revealed a new MAX2 interactor, PAPP5, that might act through dephosphorylation of MAX2 to control mainly KAR/KL-related phenotypes and, hence, provide another link with the light pathway.

MAX2-dependent competence for callus formation and shoot regeneration from Arabidopsis thaliana root explants.

Although the division of the pericycle cells initiates both lateral root development and root-derived callus formation, these developmental processes are affected differently in the strigolactone and karrikin/KARRIKIN INSENSITIVE 2 (KAI2) ligand signalling mutant more axillary growth 2 (max2). Whereas max2 produces more lateral roots than the wild type, it is defective in the regeneration of shoots from root explants. We suggest that the decreased shoot regeneration of max2 originates from delayed formation of callus primordium, yielding less callus material to regenerate shoots. Indeed, when incubated on callus-inducing medium, the pericycle cell division was reduced in max2 and the early gene expression varied when compared with the wild type, as determined by a transcriptomics analysis. Furthermore, the expression of the LATERAL ORGAN BOUNDARIES DOMAIN genes and of callus-induction genes was modified in correlation with the max2 phenotype, suggesting a role for MAX2 in the regulation of the interplay between cytokinin, auxin, and light signalling in callus initiation. Additionally, we found that the in vitro shoot regeneration phenotype of max2 might be caused by a defect in KAI2, rather than in DWARF14, signalling. Nevertheless, the shoot regeneration assays revealed that the strigolactone biosynthesis mutants max3 and max4 also play a minor role.

Rapid toxicity assessment of six antifouling booster biocides using a microplate-based chlorophyll fluorescence in Undaria pinnatifida gametophytes.

Biocides of antifouling agents can cause problems in marine ecosystems by damaging to non-target algal species. Aquatic bioassays are important means of assessing the quality of water containing mixtures of contaminants and of providing a safety standard for water management in an ecological context. In this study, a rapid, sensitive and inexpensive test method was developed using free-living male and female gametophytes of the brown macroalga Undaria pinnatifida. A conventional fluorometer was employed to evaluate the acute (48 h) toxic effects of six antifouling biocides: 4,5-Dichloro-2-octyl-isothiazolone (DCOIT), diuron, irgarol, medetomidine, tolylfluanid, zinc pyrithione (ZnPT). The decreasing toxicity in male and female gametophytes as estimated by EC50 (effective concentration at which 50% inhibition occurs) values was: diuron (0.037 and 0.128 mg l-1, respectively) > irgarol (0.096 and 0.172 mg l-1, respectively) > tolylfluanid (0.238 and 1.028 mg l-1, respectively) > DCOIT (1.015 and 0.890 mg l-1, respectively) > medetomidine (12.032 and 12.763 mg l-1, respectively). For ZnPT, 50% fluorescence inhibition of U. pinnatifida gametophytes occurred at concentrations above 0.4 mg l-1. The Undaria method is rapid, simple, practical, and cost-effective for the detection of photosynthesis-inhibiting biocides, thus making a useful tool for testing the toxicity of antifouling agents in marine environments.

Towards a multi-bioassay-based index for toxicity assessment of fluvial waters.

Despite their proven reliability for revealing 'acceptable' degrees of toxicity in waste- and reclaimed waters, bioassays are rarely used to assess the toxicity of hazardous contaminants present in natural waters. In this study, we used organisms from different trophic levels to assess the toxicity of water samples collected from four different South Korean rivers. The main objective was to develop a multi-descriptor index of toxicity for undiluted river water. The responses of six test organisms (Aliivibrio fischeri, Pseudokirchneriella subcapitata, Heterocypris incongruens, Moina macrocopa, Danio rerio and Lemna minor) after laboratory exposure to water samples were considered for this index, as well as the frequency of teratologies in diatom assemblages. Each individual test was attributed a toxicity class and score (three levels; no toxicity = 0, low toxicity = 1, confirmed toxicity = 2) based on the organism's response after exposure and a total score was calculated. The proposed index also considers the number of test organisms that received the highest toxicity score (value = 2). An overall toxicity category was then attributed to the water sample based on those two metrics: A = no toxicity, B = slight toxicity, C = moderate toxicity; D = toxicity and E = high toxicity. The susceptibility of the test organisms varied greatly and the sensitivity of their response also differed among bioassays. The combined responses of organisms from different trophic levels and with different life strategies provided multi-level diagnostic information about the intensity and the nature of contamination.

Testing the toxicity of metals, phenol, effluents, and receiving waters by root elongation in Lactuca sativa L.

Phytotoxicity tests using higher plants are among the most simple, sensitive, and cost-effective of the methods available for ecotoxicity testing. In the present study, a hydroponic-based phytotoxicity test using seeds of Lactuca sativa was used to evaluate the water quality of receiving waters and effluents near two industrial sites (Soyo and Daejon) in Korea with respect to the toxicity of 10 metals (As, Cd, Cr, Cu, Fe, Pb, Mn, Hg, Ni, Zn) and phenol, and of the receiving waters and effluents themselves. First, the L. sativa hydroponic bioassay was used to determine whether the receiving water or effluents were toxic; then, the responsible toxicant was identified. The results obtained with the L. sativa bioassay ranked the EC50 toxicities of the investigated metal ions and phenol as: Cd > Ni > Cu > Zn > Hg > phenol > As > Mn > Cr > Pb > Fe. We found that Zn was the toxicant principally responsible for toxicity in Daejeon effluents. The Daejeon field effluent had a higher Zn concentration than permitted by the effluent discharge criteria of the Ministry of Environment of Korea. Our conclusion on the importance of Zn toxicity was supported by the results of the L. sativa hydroponic assay, which showed that the concentration of Zn required to inhibit root elongation in L. sativa by 50% (EC50) was higher in the Daejeon field effluent than that of pure Zn. More importantly, we proved that the L. sativa hydroponic test method can be applied not only as an alternative tool for determining whether a given waste is acceptable for discharge into public water bodies, but also as an alternative method for measuring the safety of aquatic environments using EC20 values, with respect to the water pollutants investigated (i.e., Cd, Cr, Cu, Pb, Mn, Hg, Ni, Zn, and phenol).

Combining linkage and association mapping identifies RECEPTOR-LIKE PROTEIN KINASE1 as an essential Arabidopsis shoot regeneration gene.

De novo shoot organogenesis (i.e., the regeneration of shoots on nonmeristematic tissue) is widely applied in plant biotechnology. However, the capacity to regenerate shoots varies highly among plant species and cultivars, and the factors underlying it are still poorly understood. Here, we evaluated the shoot regeneration capacity of 88 Arabidopsis thaliana accessions and found that the process is blocked at different stages in different accessions. We show that the variation in regeneration capacity between the Arabidopsis accessions Nok-3 and Ga-0 is determined by five quantitative trait loci (QTL): REG-1 to REG-5. Fine mapping by local association analysis identified RECEPTOR-LIKE PROTEIN KINASE1 (RPK1), an abscisic acid-related receptor, as the most likely gene underlying REG-1, which was confirmed by quantitative failure of an RPK1 mutation to complement the high and low REG-1 QTL alleles. The importance of RPK1 in regeneration was further corroborated by mutant and expression analysis. Altogether, our results show that association mapping combined with linkage mapping is a powerful method to discover important genes implicated in a biological process as complex as shoot regeneration.

Molecular genetic framework for protophloem formation.

The phloem performs essential systemic functions in tracheophytes, yet little is known about its molecular genetic specification. Here we show that application of the peptide ligand CLAVATA3/embryo surrounding region 45 (CLE45) specifically inhibits specification of protophloem in Arabidopsis roots by locking the sieve element precursor cell in its preceding developmental state. CLE45 treatment, as well as viable transgenic expression of a weak CLE45(G6T) variant, interferes not only with commitment to sieve element fate but also with the formative sieve element precursor cell division that creates protophloem and metaphloem cell files. However, the absence of this division appears to be a secondary effect of discontinuous sieve element files and subsequent systemically reduced auxin signaling in the root meristem. In the absence of the formative sieve element precursor cell division, metaphloem identity is seemingly adopted by the normally procambial cell file instead, pointing to possibly independent positional cues for metaphloem formation. The protophloem formation and differentiation defects in brevis radix (brx) and octopus (ops) mutants are similar to those observed in transgenic seedlings with increased CLE45 activity and can be rescued by loss of function of a putative CLE45 receptor, barely any meristem 3 (BAM3). Conversely, a dominant gain-of-function ops allele or mild OPS dosage increase suppresses brx defects and confers CLE45 resistance. Thus, our data suggest that delicate quantitative interplay between the opposing activities of BAM3-mediated CLE45 signals and OPS-dependent signals determines cellular commitment to protophloem sieve element fate, with OPS acting as a positive, quantitative master regulator of phloem fate.

Strigolactones spatially influence lateral root development through the cytokinin signaling network.

Strigolactones are important rhizosphere signals that act as phytohormones and have multiple functions, including modulation of lateral root (LR) development. Here, we show that treatment with the strigolactone analog GR24 did not affect LR initiation, but negatively influenced LR priming and emergence, the latter especially near the root-shoot junction. The cytokinin module ARABIDOPSIS HISTIDINE KINASE3 (AHK3)/ARABIDOPSIS RESPONSE REGULATOR1 (ARR1)/ARR12 was found to interact with the GR24-dependent reduction in LR development, because mutants in this pathway rendered LR development insensitive to GR24. Additionally, pharmacological analyses, mutant analyses, and gene expression analyses indicated that the affected polar auxin transport stream in mutants of the AHK3/ARR1/ARR12 module could be the underlying cause. Altogether, the data reveal that the GR24 effect on LR development depends on the hormonal landscape that results from the intimate connection with auxins and cytokinins, two main players in LR development.

Suppression of Arabidopsis protophloem differentiation and root meristem growth by CLE45 requires the receptor-like kinase BAM3.

Peptide signaling presumably occupies a central role in plant development, yet only few concrete examples of receptor-ligand pairs that act in the context of specific differentiation processes have been described. Here we report that second-site null mutations in the Arabidopsis leucine-rich repeat receptor-like kinase gene barely any meristem 3 (BAM3) perfectly suppress the postembryonic root meristem growth defect and the associated perturbed protophloem development of the brevis radix (brx) mutant. The roots of bam3 mutants specifically resist growth inhibition by the CLAVATA3/ENDOSPERM SURROUNDING REGION 45 (CLE45) peptide ligand. WT plants transformed with a construct for ectopic overexpression of CLE45 could not be recovered, with the exception of a single severely dwarfed and sterile plant that eventually died. By contrast, we obtained numerous transgenic bam3 mutants transformed with the same construct. These transgenic plants displayed a WT phenotype, however, supporting the notion that CLE45 is the likely BAM3 ligand. The results correlate with the observation that external CLE45 application represses protophloem differentiation in WT, but not in bam3 mutants. BAM3, BRX, and CLE45 are expressed in a similar spatiotemporal trend along the developing protophloem, up to the end of the transition zone. Induction of BAM3 expression upon CLE45 application, ectopic overexpression of BAM3 in brx root meristems, and laser ablation experiments suggest that intertwined regulatory activity of BRX, BAM3, and CLE45 could be involved in the proper transition of protophloem cells from proliferation to differentiation, thereby impinging on postembryonic growth capacity of the root meristem.

The sensitivity of an hydroponic lettuce root elongation bioassay to metals, phenol and wastewaters.

The root elongation bioassay is one of the most straightforward test methods used for environmental monitoring in terms of simplicity, rapidity and economy since it merely requires filter paper, distilled water and Petri dishes. However, filter paper as a support material is known to be problematic as it can reduce the sensitivity of the test. The newly developed hydroponic method reported here differs from the conventional root elongation method (US EPA filter paper method) in that no support material is used and the exposure time is shorter (48 h in this test versus 120 h in the US EPA test). For metals, the hydroponic test method was 3.3 (for Hg) to 57 (for Cu) times more sensitive than the US EPA method with the rank orders of sensitivity, estimated from EC50 values, being Cu≥Cd>Ni≥Zn≥Hg for the former and Hg≥Cu≥Ni≥Cd≥Zn for the latter methods. For phenol, the results did not differ significantly; EC50 values were 124 mg L(-1) and 108-180 mg L(-1) for the hydroponic and filter paper methods, respectively. Lettuce was less sensitive than daphnids to wastewaters, but the root elongation response appears to be wastewater-specific and is especially sensitive for detecting the presence of fluorine. The new hydroponic test thus provides many practical advantages, especially in terms of cost and time-effectiveness requiring only a well plate, a small volume of distilled water and short exposure period; furthermore, no specialist expertise is required. The method is simpler than the conventional EPA technique in not using filter paper which can influence the sensitivity of the test. Additionally, plant seeds have a long shelf-life and require little or no maintenance.

The IBO germination quantitative trait locus encodes a phosphatase 2C-related variant with a nonsynonymous amino acid change that interferes with abscisic acid signaling.

Natural genetic variation is crucial for adaptability of plants to different environments. Seed dormancy prevents precocious germination in unsuitable conditions and is an adaptation to a major macro-environmental parameter, the seasonal variation in temperature and day length. Here we report the isolation of IBO, a quantitative trait locus (QTL) that governs c. 30% of germination rate variance in an Arabidopsis recombinant inbred line (RIL) population derived from the parental accessions Eilenburg-0 (Eil-0) and Loch Ness-0 (Lc-0). IBO encodes an uncharacterized phosphatase 2C-related protein, but neither the Eil-0 nor the Lc-0 variant, which differ in a single amino acid, have any appreciable phosphatase activity in in vitro assays. However, we found that the amino acid change in the Lc-0 variant of the IBO protein confers reduced germination rate. Moreover, unlike the Eil-0 variant of the protein, the Lc-0 variant can interfere with the activity of the phosphatase 2C ABSCISIC ACID INSENSITIVE 1 in vitro. This suggests that the Lc-0 variant possibly interferes with abscisic acid signaling, a notion that is supported by physiological assays. Thus, we isolated an example of a QTL allele with a nonsynonymous amino acid change that might mediate local adaptation of seed germination timing.

Strigolactones as an auxiliary hormonal defence mechanism against leafy gall syndrome in Arabidopsis thaliana.

Leafy gall syndrome is the consequence of modified plant development in response to a mixture of cytokinins secreted by the biotrophic actinomycete Rhodococcus fascians. The similarity of the induced symptoms with the phenotype of plant mutants defective in strigolactone biosynthesis and signalling prompted an evaluation of the involvement of strigolactones in this pathology. All tested strigolactone-related Arabidopsis thaliana mutants were hypersensitive to R. fascians. Moreover, treatment with the synthetic strigolactone mixture GR24 and with the carotenoid cleavage dioxygenase inhibitor D2 illustrated that strigolactones acted as antagonistic compounds that restricted the morphogenic activity of R. fascians. Transcript profiling of the MORE AXILLARY GROWTH1 (MAX1), MAX2, MAX3, MAX4, and BRANCHED1 (BRC1) genes in the wild-type Columbia-0 accession and in different mutant backgrounds revealed that upregulation of strigolactone biosynthesis genes was triggered indirectly by the bacterial cytokinins via host-derived auxin and led to the activation of BRC1 expression, inhibiting the outgrowth of the newly developing shoots, a typical hallmark of leafy gall syndrome. Taken together, these data support the emerging insight that balances are critical for optimal leafy gall development: the long-lasting biotrophic interaction is possible only because the host activates a set of countermeasures-including the strigolactone response-in reaction to bacterial cytokinins to constrain the activity of R. fascians.

Interactive Effects of Blue Light and Water Turbulence on the Growth of the Green Macroalga Ulva australis (Chlorophyta).

Macroalgal growth and yield are key to sustainable aquaculture. Although light and water turbulence are two important factors that affect algal productivity, research on their interaction is limited. Therefore, in this study, we investigated the effects of different wavelengths of light and the presence or absence of water turbulence on the growth of the green macroalga Ulva australis. Water turbulence was found to enhance the growth of U. australis irrespective of photosynthetic performance, but only in blue light cultures. The quantum dose of blue light required to induce 50% growth promotion was 1.02 mol m-2, which is comparable to the reported values for cryptochrome-mediated effects in other macroalgae. The combined effect of blue light and water turbulence led to the accumulation of photosynthesis-related proteins that support plastid differentiation and facilitate efficient photosynthesis and growth. Our findings thus highlight the potential of harnessing blue light and water turbulence to maximise macroalgal cultivation for sustainable and profitable algal aquaculture.

Strigolactones fine-tune the root system.

Strigolactones were originally discovered to be involved in parasitic weed germination, in mycorrhizal association and in the control of shoot architecture. Despite their clear role in rhizosphere signaling, comparatively less attention has been given to the belowground function of strigolactones on plant development. However, research has revealed that strigolactones play a key role in the regulation of the root system including adventitious roots, primary root length, lateral roots, root hairs and nodulation. Here, we review the recent progress regarding strigolactone regulation of the root system and the antagonism and interplay with other hormones.

Metal Toxicity across Different Thallus Sections of the Green Macroalga, Ulva australis.

We aimed to identify functional differences between different sections of the thallus of Ulva australis and develop tissue-endpoint combinations to assess the toxicity of six metals (i.e., Ag, As, Cd, Cr, Cu, and Ni). EC50 values for these metals in three sections of the thallus of Ulva were obtained for multiple endpoints: relative growth rate (RGR), chlorophyll a fluorescence, pigment contents, and the expression of the photosynthesis-related gene, rbcL. The responses of the endpoints varied across the respective thallus sections; overall, the most toxic metals were Ag and Cu. These endpoints were the best for evaluating metal toxicity: ETRmax of the middle thallus sections for Ag toxicity; RGR of the middle thallus section for As and Cd; ETRmax of the marginal thallus section for Cr; Chl b contents of the marginal thallus section for Cu; RGR of the basal thallus section for Ni. The EC50 values for the inhibition of ETRmax in middle (0.06 mg∙L-1) and Chl b in the marginal thallus sections (0.06 mg∙L-1) were all lower than those of the quality standard for wastewater discharge values of Ag and Cu in Republic of Korea and the US, pointing to the suitability of U. australis-based endpoints for risk assessment.

How to make sense of 3D representations for plant phenotyping: a compendium of processing and analysis techniques.

Computer vision technology is moving more and more towards a three-dimensional approach, and plant phenotyping is following this trend. However, despite its potential, the complexity of the analysis of 3D representations has been the main bottleneck hindering the wider deployment of 3D plant phenotyping. In this review we provide an overview of typical steps for the processing and analysis of 3D representations of plants, to offer potential users of 3D phenotyping a first gateway into its application, and to stimulate its further development. We focus on plant phenotyping applications where the goal is to measure characteristics of single plants or crop canopies on a small scale in research settings, as opposed to large scale crop monitoring in the field.

Giant Duckweed (Spirodela polyrhiza) Root Growth as a Simple and Sensitive Indicator of Copper and Chromium Contamination.

Aquatic environment are often contaminated with heavy metals from various industrial sources. However, physicochemical techniques for pollutant detection are limited, thus prompting the need for additional bioassays. We investigated the use of greater duckweed (Spirodela polyrhiza) as a bioindicator of metal pollution. We exposed S. polyrhiza to four pollutants (namely, silver, cadmium, copper, and chromium) and assessed metal toxicity by measuring its frond area and the length of its regrown roots. The plant displayed significant differences in both frond size and root growth in response to the four metals. Silver was the most toxic (EC50 = 23 µg L-1) while copper the least (EC50 = 365-607 µg L-1). Direct comparisons of metal sensitivity and the reliability of the two endpoint assays showed that root growth was more sensitive (lower in terms of 50% effective concentration) to chromium, cadmium, and copper, and was more reliable (lower in terms of coefficient of variation) than those for frond area. Compared to conventional Lemna-based tests, the S. polyrhiza test is easier to perform (requiring only one 24-well plate, 3 mL of medium and a 72-h exposure). Moreover, it does not require livestock cultivation/maintenance, making it more suitable for repeated measurements. Measurements of S. polyrhiza root length may be suitable for assessment when copper and chromium in municipal and industrial wastewater exceed the environmentally permissible levels.

Evaluating ecotoxicological assays for comprehensive risk assessment of toxic metals present in industrial wastewaters in the Republic of Korea.

Toxicity tests represent a rapid, user-friendly and cost-effective means to assess the impact of wastewater quality on aquatic ecosystems. There are not many cases where wastewater management standards are set based on various bio-based ecotoxicity values. Here, we tested a novel multitaxon approach to compare standard water quality indices to toxicity metrics obtained from ecotoxicity tests, conducted using aquatic organisms representing several trophic levels (Aliivibrio, Ulva, Daphnia, and Lemna), for 99 industrial wastewater samples from South Korea. For five wastewater samples, the concentrations of Se, Zn, or Ni exceeded the permissible limits (1, 5, and 3 mg L-1, respectively). All the four physiochemical water quality indices tested were positively correlated with Se and Pb concentrations. The toxicity unit (TU) scores indicated a declining sensitivity to pollutants, in the order Lemna (2.87) >Daphnia (2.24) >Aliivibrio (1.78) >Ulva (1.42). Significant correlations were observed between (1) Cd and Ni, and Aliivibrio, (2) Cu and Daphnia, (3) Cd, Cu, Zn, and Cr and Lemna, and (4) Cu, Zn, and Ni and Ulva. Daphnia-Lemna and Lemna-Ulva were found to be good indicators of ecologically harmful Se and Ni contents in wastewater, respectively. We suggest that regulatory thresholds based on these bioassays should be set at TU = 1 for all the species or at TU = 1 for Aliivibrio and Ulva and TU = 2 for Daphnia and Lemna, if the number of companies whose wastewater discharge exceeds the allowable TU levels is <1 % or 5 % of the total number of industries, respectively. Taken together, these findings could help in establishing a rapid, ecologically relevant wastewater quality assessment system that would be useful for developing strategies to protect aquatic ecosystems.