Current status of community resources and priorities for weed genomics research.

Weeds are attractive models for basic and applied research due to their impacts on agricultural systems and capacity to swiftly adapt in response to anthropogenic selection pressures. Currently, a lack of genomic information precludes research to elucidate the genetic basis of rapid adaptation for important traits like herbicide resistance and stress tolerance and the effect of evolutionary mechanisms on wild populations. The International Weed Genomics Consortium is a collaborative group of scientists focused on developing genomic resources to impact research into sustainable, effective weed control methods and to provide insights about stress tolerance and adaptation to assist crop breeding.

Metabolic Exploration for Cyhalofop-Butyl Antagonism in Barnyardgrass [Echinochloa crus-galli (L.) P. Beauv.] Following Pretreatment of Malathion.

The present study investigated the effects of broad-spectrum metabolic inhibitors malathion (cytochrome P450 inhibitor) and/or 4-chloro-7-nitrobenzofurazan (NBD-Cl; glutathione S-transferase inhibitor) on the metabolism of cyhalofop-butyl (CyB) in barnyardgrass [Echinochloa crus-galli (L.) P. Beauv.] biotypes confirmed previously with multiple resistance to two herbicides CyB and florpyrauxifen-benzyl. The metabolic inhibitors were not effective at recovering the sensitivity of resistant barnyardgrass biotypes to CyB treated at the labeled rate (313 g ai ha-1). Rather, treatment with malathion followed by CyB caused antagonism, reducing the efficacy of CyB and promoting the growth of resistant biotypes. Pretreatment with malathion did not influence absorption/translocation of the applied form CyB and its conversion to the active herbicide form cyhalofop-acid (CyA), in both susceptible and resistant biotypes. In contrast, metabolism of the applied form (CyB) decreased 1.5 to 10.5 times by the malathion pretreatment. Taken together, the maintained CyA production against the reduced CyB metabolism could be the mechanism to account for the cause of CyB antagonism observed in barnyardgrass following malathion pretreatment. Additionally, the evolution of CyB resistance in barnyardgrass might be associated with reduced production of CyA in resistant biotypes, independent of activities of cytochrome P450 or GST enzymes.

Exploratory Analysis on Herbicide Metabolism and Very-Long-Chain Fatty Acid Production in Metolachlor-Resistant Palmer Amaranth (Amaranthus palmeri S. Wats.).

A Palmer amaranth (Amaranthus palmeri S. Wats.) biotype resistant to S-metolachlor was confirmed from crop fields in Arkansas, USA. This study investigated the metabolic effects of malathion (cytochrome P450 inhibitor) and 4-chloro-7-nitrobenzofurazan [NBD-Cl; glutathione S-transferase inhibitor] on the S-metolachlor-resistant A. palmeri biotype. Root elongation of the resistant biotype was 20% more inhibited by treatment of NBD-Cl (50 nM) and S-metolachlor (2 µM) in mixture than by treatment of S-metolachlor alone. Metabolites of S-metolachlor were 1.4-12.1 times greater produced in the resistant biotype for 7 d than in the susceptible standard. Production of cerotic acid, one of the very-long-chain fatty acids containing 26 carbons, was more reduced in the susceptible standard (3.8-fold) than in the resistant biotype (1.8-fold) by S-metolachlor treatment. Conclusively, evolution of S-metolachlor resistance observed in this study was likely associated with improved activity of glutathione S-transferases. Further studies are needed to genetically evaluate plant endogenous enzymes involving cerotic acid production.

Detection of subgenome bias using an anchored syntenic approach in Eleusine coracana (finger millet).

BACKGROUND: Finger millet (Eleusine coracana 2n = 4x = 36) is a hardy, nutraceutical, climate change tolerant, orphan crop that is consumed throughout eastern Africa and India. Its genome has been sequenced multiple times, but A and B subgenomes could not be separated because no published genome for E. indica existed. The classification of A and B subgenomes is important for understanding the evolution of this crop and provide a means to improve current and future breeding programs. RESULTS: We produced subgenome calls for 704 syntenic blocks and inferred A or B subgenomic identity for 59,377 genes 81% of the annotated genes. Phylogenetic analysis of a super matrix containing 455 genes shows high support for A and B divergence within the Eleusine genus. Synonymous substitution rates between A and B genes support A and B calls. The repetitive content on highly supported B contigs is higher than that on similar A contigs. Analysis of syntenic singletons showed evidence of biased fractionation showed a pattern of A genome dominance, with 61% A, 37% B and 1% unassigned, and was further supported by the pattern of loss observed among cyto-nuclear interacting genes. CONCLUSION: The evidence of individual gene calls within each syntenic block, provides a powerful tool for inference for subgenome classification. Our results show the utility of a draft genome in resolving A and B subgenomes calls, primarily it allows for the proper polarization of A and B syntenic blocks. There have been multiple calls for the use of phylogenetic inference in subgenome classification, our use of synteny is a practical application in a system that has only one parental genome available.

A novel mutation A212T in chloroplast Protoporphyrinogen oxidase (PPO1) confers resistance to PPO inhibitor Oxadiazon in Eleusine indica.

BACKGROUND: Protoporphyrinogen oxidase (PPO) with two isoforms, chloroplast-targeted (PPO1) and mitochondrial-targeted (PPO2), catalyzes a step in the biosynthesis of chlorophyll and heme. PPO1 and PPO2 are herbicide target sites of PPO-inhibiting herbicides. Target-site mutations conferring resistance to PPO inhibitors have all thus far been in PPO2. Oxadiazon is a unique PPO inhibitor utilized for preemergence Eleusine indica control. In this research, we evaluated the response of two previously confirmed oxadiazon-resistant and susceptible E. indica biotypes to other PPO inhibitors and identified the resistance mechanism in two oxadiazon-resistant E. indica biotypes. RESULTS: Two E. indica biotypes were resistant to oxadiazon, but not to other structurally unrelated PPO inhibitors, such as lactofen, flumioxazin and sulfentrazone. A novel mutation A212T was identified in the chloroplast-targeted PPO1, conferring resistance to oxadiazon in a heterologous expression system. Computational structural modeling provided a mechanistic explanation for reduced herbicide binding to the variant protein: the presence of a methyl group of threonine 212 changes the PPO1 active site and produces repulsive electrostatic interactions that repel oxadiazon from the binding pocket. CONCLUSION: The novel A212T mutation in PPO1 conferring resistance specifically to PPO inhibitor oxadiazon was characterized. This is the first evidence of the direct role of PPO1 in the PPO mode of action, and the first evidence of evolved resistance in PPO1. © 2019 Society of Chemical Industry.

Transcriptome Analysis Reveals Unique Relationships Among Eleusine Species and Heritage of Eleusine coracana.

Relationships in the genus Eleusine were obtained through transcriptome analysis. Eleusine coracana (E. coracana ssp. coracana), also known as finger millet, is an allotetraploid minor crop primarily grown in East Africa and India. Domesticated E. coracana evolved from wild E. africana (E. coracana ssp. africana) with the maternal genome donor largely supported to be E. indica; however, the paternal genome donor remains elusive. We developed transcriptomes for six Eleusine species from fully developed seedlings using Illumina technology and three de novo assemblers (Trinity, Velvet, and SOAPdenovo2) with the redundancy-reducing EvidentialGene pipeline. Mapping E. coracana reads to the chloroplast genes of all Eleusine species detected fewer variants between E. coracana and E. indica compared to all other species. Phylogenetic analysis further supports E. indica as the maternal parent of E. coracana and E. africana, in addition to a close relationship between E. indica and E. tristachya, and between E. floccifolia and E. multiflora, and E. intermedia as a separate group. A close relationship between E. floccifolia and E. multiflora was unexpected considering they are reported to have distinct nuclear genomes, BB and CC, respectively. Further, it was expected that E. intermedia and E. floccifolia would have a closer relationship considering they have similar nuclear genomes, AB and BB, respectively. A rethinking of the labeling of ancestral genomes of E. floccifolia, E. multiflora, and E. intermedia is maybe needed based on this data.

Rapid identification of a candidate nicosulfuron sensitivity gene (Nss) in maize (Zea mays L.) via combining bulked segregant analysis and RNA-seq.

KEY MESSAGE: A candidate nicosulfuron sensitivity gene Nss was identified by combining bulked segregant analysis and RNA-seq. Multiple mutations of this gene were discovered in nicosulfuron-sensitive maize compared with the tolerant. It has been demonstrated that variabilities exist in maize response to nicosulfuron. Two nicosulfuron-sensitive inbred lines (HB39, HB41) and two tolerant inbred lines (HB05, HB09) were identified via greenhouse and field trials. Genetic analysis indicated that the sensitivity to nicosulfuron in maize was controlled by a single, recessive gene. To precisely and rapidly map the nicosulfuron sensitivity gene (Nss), two independent F2 segregating populations, Population A (HB41 × HB09) and Population B (HB39 × HB05), were constructed. By applying bulked segregant RNA-Seq (BSR-Seq), the Nss gene was, respectively, mapped on the short arm of chromosome 5 (chr5: 1.1-15.3 Mb) and (chr5: 0.5-18.2 Mb) using two populations, with 14.2 Mb region in common. Further analysis revealed that there were 43 and 119 differentially expressed genes in the mapping intervals, with 18 genes in common. Gene annotation results showed that a cytochrome P450 gene (CYP81A9) appeared to be the candidate gene of Nss associated with nicosulfuron sensitivity in maize. Sequence analysis demonstrated that two common deletion mutations existed in the sensitive maize, which might lead to the nicosulfuron sensitivity in maize. The results will make valuable contributions to the understanding of molecular mechanism of herbicide sensitivity in maize.

The power and potential of genomics in weed biology and management.

There have been previous calls for, and efforts focused on, realizing the power and potential of weed genomics for better understanding of weeds. Sustained advances in genome sequencing and assembly technologies now make it possible for individual research groups to generate reference genomes for multiple weed species at reasonable costs. Here, we present the outcomes from several meetings, discussions, and workshops focused on establishing an International Weed Genomics Consortium (IWGC) for a coordinated international effort in weed genomics. We review the 'state of the art' in genomics and weed genomics, including technologies, applications, and on-going weed genome projects. We also report the outcomes from a workshop and a global survey of the weed science community to identify priority species, key biological questions, and weed management applications that can be addressed through greater availability of, and access to, genomic resources. Major focus areas include the evolution of herbicide resistance and weedy traits, the development of molecular diagnostics, and the identification of novel targets and approaches for weed management. There is increasing interest in, and need for, weed genomics, and the establishment of the IWGC will provide the necessary global platform for communication and coordination of weed genomics research. © 2018 Society of Chemical Industry.

Complete plastid genome sequence of goosegrass (Eleusine indica) and comparison with other Poaceae.

Eleusine indica, also known as goosegrass, is a serious weed in at least 42 countries. In this paper we report the complete plastid genome sequence of goosegrass obtained by de novo assembly of paired-end and mate-paired reads generated by Illumina sequencing of total genomic DNA. The goosegrass plastome is a circular molecule of 135,151bp in length, consisting of two single-copy regions separated by a pair of inverted repeats (IRs) of 20,919 bases. The large (LSC) and the small (SSC) single-copy regions span 80,667 bases and 12,646 bases, respectively. The plastome of goosegrass has 38.19% GC content and includes 108 unique genes, of which 76 are protein-coding, 28 are transfer RNA, and 4 are ribosomal RNA. The goosegrass plastome sequence was compared to eight other species of Poaceae. Although generally conserved with respect to Poaceae, this genomic resource will be useful for evolutionary studies within this weed species and the genus Eleusine.

Transcriptome Assembly and Comparison of an Allotetraploid Weed Species, Annual Bluegrass, with its Two Diploid Progenitor Species, Schrad and Kunth.

Annual bluegrass ( L.) is one of the most widespread weed species in this world. As a young allotetraploid, has occupied diverse environments from Antarctic area to subtropical regions. To unveil the evolutionary mystery behind 's wide distribution, extensive adaptability and phenotypic plasticity needs collaboration from multiple research scopes from ecology and plant physiology to population genetics and molecular biology. However, the lack of omic data and reference has greatly hampered the study. This is the first comprehensive transcriptome study on species. Total RNA was extracted from and its two proposed diploid parents, Schrad and Kunth, and sequenced in Illumina Hiseq2000. Optimized, nonredundant transcriptome references were generated for each species using four de novo assemblers (Trinity, Velvet, SOAPdenovo, and CLC Genomics Workbench) and a redundancy-reducing pipeline (CD-HIT-EST and EvidentialGene tr2aacds). Using the constructed transcriptomes together with sequencing reads, we found high similarity in nucleotide sequences and homeologous polymorphisms between and the two proposed parents. Comparison of chloroplast and mitochondrion genes further confirmed as the maternal parent. Less nucleotide percentage differences were observed between and homeologs than between and homeologs, indicating a higher nucleotide substitution rates in homeologs than in homeologs. Gene ontology (GO) enrichment analysis suggested the more compatible cytoplasmic environment and cellular apparatus for homeologs as the major cause for this phenomenon.

Optimizing Transcriptome Assemblies for Eleusine indica Leaf and Seedling by Combining Multiple Assemblies from Three De Novo Assemblers.

Due to rapid advances in sequencing technology, increasing amounts of genomic and transcriptomic data are available for plant species, presenting enormous challenges for biocomputing analysis. A crucial first step for a successful transcriptomics-based study is the building of a high-quality assembly. Here, we utilized three different de novo assemblers (Trinity, Velvet, and CLC) and the EvidentialGene pipeline tr2aacds to assemble two optimized transcript sets for the notorious weed species, Eleusine indica. Two RNA sequencing (RNA-seq) datasets from leaf and aboveground seedlings were processed using three assemblers, which resulted in 20 assemblies for each dataset. The contig numbers and N50 values of each assembly were compared to study the effect of read number, k-mer size, and in silico normalization on assembly output. The 20 assemblies were then processed through the tr2aacds pipeline to remove redundant transcripts and to select the transcript set with the best coding potential. Each assembly contributed a considerable proportion to the final transcript combination with the exception of the CLC-k14. Thus each assembler and parameter set did assemble better contigs for certain transcripts. The redundancy, total contig number, N50, fully assembled contig number, and transcripts related to target-site herbicide resistance were evaluated for the EvidentialGene and Trinity assemblies. Comparing the EvidentialGene set with the Trinity assembly revealed improved quality and reduced redundancy in both leaf and seedling EvidentialGene sets. The optimized transcriptome references will be useful for studying herbicide resistance in E. indica and the evolutionary process in the three allotetraploid E. indica offspring.

Utilizing next-generation sequencing to study homeologous polymorphisms and herbicide-resistance-endowing mutations in Poa annua acetolactate synthase genes.

BACKGROUND: Detection of single nucleotide polymorphisms (SNPs) related to herbicide resistance in non-model polyploid weed species is fraught with difficulty owing to the gene duplication and lack of reference sequences. Our research seeks to overcome these obstacles by Illumina HiSeq read mapping, SNP calling and allele frequency determinations. Our focus is on the acetolactate synthase (ALS) gene, the target site of ALS-inhibiting herbicides, in Poa annua, an allotetraploid weed species originating from two diploid parents, P. supina and P. infirma. RESULTS: ALS contigs with complete coding regions of P. supina, P. infirma and P. annua were assembled and compared with ALS genes from other plant species. The ALS infirma-homeolog of P. annua showed higher levels of nucleotide sequence variability than the supina-homeolog. Comparisons of read mappings of P. annua and a simulated P. supina × P. infirma hybrid showed high resemblance. Two homeolog-specific primer pairs were designed and used to amplify a 1860 bp region covering all resistance-conferring codons in the ALS gene. Four P. annua populations, GN, RB, GW and LG, showed high resistance to two ALS inhibitors, bispyribac-sodium and foramsulfuron, and two populations, HD and RS, showed lower resistance in the rate-response trial. Mutations conferring Trp-574-Leu substitution were observed in the infirma-homeolog of GN and RB and in the supina-homeolog of GW and LG, but no resistance-conferring mutation was observed in the two populations of lower resistance, HD and RS. CONCLUSION: In this study we have demonstrated the use of NGS data to study homeologous polymorphisms, parentage and herbicide resistance in an allotetraploid weed species, P. annua. Complete coding sequences of the ALS gene were assembled for P. infirma, P. supina, infirma-homeolog and supina-homeolog in P. annua. A pipeline consisting of read mapping, SNP calling and allele frequency calculation was developed to study the parentage of P. annua, which provided a new perspective on this topic besides the views of morphology, karyotype and phylogeny. Our two homeolog-specific primer pairs can be utilized in future research to separate the homeologs of the ALS gene in P. annua and cover all the codons that have been reported to confer herbicide resistance.

Effects of elevated CO2 on biomass and fungi associated with two ecotypes of ragweed (Ambrosia artemisiifolia L.).

Herbicide resistant weed populations have developed due to the repeated application of herbicides. Elevated concentrations of atmospheric CO2 can have positive effects on weed growth, but how rising CO2 might affect herbicide resistant weeds is not known. Ragweed (Ambrosia artemisiifolia L.) ecotypes known to be resistant or susceptible to glyphosate herbicide were exposed to either ambient or elevated (ambient +200 µ mol mol(-1)) concentrations of CO2 in open top chambers. Plants were harvested following 8 weeks of CO2 exposure; at this time, they had begun to exhibit disease symptoms including spots on leaves and stems. Elevated CO2 significantly increased top, root, and total plant biomass. Also, glyphosate resistant plants had significantly greater top, root, and total biomass than plants susceptible to the herbicide. There were no significant CO2 by ecotype interactions. Fungi from 13 genera were associated with ragweed, several of which can be either pathogens (i.e., Alternaria, Fusarium, Rhizoctonia), aiding the decline in health of the ragweed plants, or saprophytes existing on dead plant tissues. The common foliar disease powdery mildew was significantly higher on susceptible compared with resistant ragweed. Susceptible plants also showed an increased frequency of Rhizoctonia on leaves and Alternaria on stems; however, Fusarium occurred more frequently on resistant ragweed leaves. Fungi were not affected by CO2 concentration or its interaction with ecotype. This study reports the first information on the effects of elevated CO2 on growth of herbicide resistant weeds. This is also the first study examining the impact of herbicide resistance and elevated CO2 on fungi associated with weeds. What effects herbicide resistance might have on plant diseases and how rising atmospheric CO2 might impact these effects needs to be addressed, not only with important weeds but also with crops.

Analysis of the Citrullus colocynthis transcriptome during water deficit stress.

Citrullus colocynthis is a very drought tolerant species, closely related to watermelon (C. lanatus var. lanatus), an economically important cucurbit crop. Drought is a threat to plant growth and development, and the discovery of drought inducible genes with various functions is of great importance. We used high throughput mRNA Illumina sequencing technology and bioinformatic strategies to analyze the C. colocynthis leaf transcriptome under drought treatment. Leaf samples at four different time points (0, 24, 36, or 48 hours of withholding water) were used for RNA extraction and Illumina sequencing. qRT-PCR of several drought responsive genes was performed to confirm the accuracy of RNA sequencing. Leaf transcriptome analysis provided the first glimpse of the drought responsive transcriptome of this unique cucurbit species. A total of 5038 full-length cDNAs were detected, with 2545 genes showing significant changes during drought stress. Principle component analysis indicated that drought was the major contributing factor regulating transcriptome changes. Up regulation of many transcription factors, stress signaling factors, detoxification genes, and genes involved in phytohormone signaling and citrulline metabolism occurred under the water deficit conditions. The C. colocynthis transcriptome data highlight the activation of a large set of drought related genes in this species, thus providing a valuable resource for future functional analysis of candidate genes in defense of drought stress.

Increase in nutritionally important sweet corn kernel carotenoids following Mesotrione and atrazine applications.

The herbicide mesotrione inhibits a critical enzyme, phytoene desaturase, in plant carotenoid biosynthesis. Mesotrione is currently labeled for selective weed control in sweet corn ( Zea mays var. rugosa). Mesotrione applied alone, or in mixtures with the photosystem II inhibitor atrazine, acted to increase concentrations of kernel antheraxanthin, lutein, and zeaxanthin carotenoids in several sweet corn genotypes. Kernel lutein and zeaxanthin levels significantly increased 15.6% after mesotrione + atrazine early postemergence applications, as compared to the control treatment. It appears that mesotrione applications resulted in greater pools of kernel carotenoids once the sweet corn genotypes expressing moderate injury overcame the initial herbicidal photo-oxidative stress. This is the first report of herbicides directly up-regulating the carotenoid biosynthetic pathway in corn kernels, which is associated with the nutritional quality of sweet corn. Enhanced accumulation of lutein and zeaxanthin is important because dietary carotenoids function in suppressing aging eye diseases such as macular degeneration, now affecting 1.75 million older Americans.

Mesotrione control and pigment concentration of large crabgrass (Digitaria sanguinalis) under varying environmental conditions.

BACKGROUND: Mesotrione is a carotenoid biosynthesis-inhibiting herbicide currently labeled for crabgrass (Digitaria spp.) control. Mesotrione control of large crabgrass has been reported to vary with temperature and relative humidity; however, the effect of irradiance on mesotrione efficacy has not previously been reported. Likewise, little is known about pigment concentrations of Digitaria spp. The present research investigated the effects of mesotrione on large crabgrass, Digitaria sanguinalis (L.) Scop., control and pigment concentrations under varying irradiance at three temperatures. RESULTS: Mesotrione (0.28 kg ha(-1)) control of large crabgrass did not differ between temperature levels (18, 26 and 32 degrees C). Control was similar at tested irradiance levels (600, 1100 and 1600 micromol m(-2) s(-1)). Mesotrione reduced large crabgrass chlorophyll a, chlorophyll b and total carotenoid concentrations, as well as chlorophyll a to b ratios. Treated plant bleaching was highest 7 days after treatment (DAT) but decreased by 21 DAT. Treated plants were less than 10% necrotic 3 and 7 DAT but nearly 35% necrotic 21 DAT. Treated large crabgrass bleaching was highest and photochemical efficiency was lowest 7 DAT. These results indicate that some plant recovery occurs prior to 21 DAT. CONCLUSION: Although mesotrione efficacy has previously been reported to vary according to environmental factors, mesotrione control of large crabgrass did not vary with measured temperature and irradiance levels in this study. On account of crabgrass convalescence, secondary applications of mesotrione may control large crabgrass more effectively when applied prior to 21 DAT.

Effects of mesotrione on perennial ryegrass (Lolium perenne L.) carotenoid concentrations under varying environmental conditions.

Mesotrione is a carotenoid biosynthesis inhibiting herbicide, which is being evaluated for use in turfgrass. Carotenoids are important light harvesting and photoprotecting pigments that dissipate and quench excess light energy. The effects of mesotrione on carotenoid concentrations in turf and weed species, such as perennial ryegrass (Lolium perenne L.), are poorly understood. Mesotrione injury to perennial ryegrass has been reported, and symptomology may differ due to postapplication environmental factors such as irradiance and temperature. Research was conducted to investigate the effects of mesotrione on perennial ryegrass under varying irradiance (600, 1100, or 1600 micromol/m (2)/s) at three different temperatures (18, 26, and 34 degrees C). Postapplication irradiance and temperature levels did not affect visual injury symptoms in perennial ryegrass. Bleaching of treated plants was highest 7 days after treatment (DAT; 8%) and recovered to nontreated levels by 21 DAT. Mesotrione applications did not decrease perennial ryegrass foliar biomass accumulations. Carotenoid concentrations of nontreated plants were similar to those reported in creeping bentgrass and many green leafy vegetable crops. However, chlorophyll a and b, beta-carotene, lutein, and violaxanthin concentrations decreased due to mesotrione applications, while phytoene and zeaxanthin, a photoprotecting carotenoid, increased. The photochemical efficiency (F v/ F m) of treated plants was lower than nontreated plants at 3 and 7 DAT; however, treated plants recovered to nontreated levels 21 DAT. Results indicate that postapplication irradiance and temperature levels may not affect mesotrione efficacy in perennial ryegrass. Preferential accumulation of zeaxanthin following mesotrione applications may be a stress-related response, which may reduce light harvesting complex size and directly quench excess light energy.

Influence of nitrogen and sulfur on biomass production and carotenoid and glucosinolate concentrations in watercress (Nasturtium officinale R. Br.).

Watercress (Nasturtium officinale R. Br.) is a perennial herb rich in the secondary metabolites of glucosinolates and carotenoids. 2-phenethyl isothiocyanate, the predominate isothiocyanate hydrolysis product in watercress, can reduce carcinogen activation through inhibition of phase I enzymes and induction of phase II enzymes. Sulfur (S) and nitrogen (N) have been shown to influence concentrations of both glucosinolates and carotenoids in a variety of vegetable crops. Our research objectives were to determine how several levels of N and S fertility interact to affect watercress plant tissue biomass production, tissue C/N ratios, concentrations of plant pigments, and glucosinolate concentrations. Watercress was grown using nutrient solution culture under a three by three factorial arrangement, with three S (8, 16, and 32 mg/L) and three N (6, 56, and 106 mg/L) fertility concentrations. Watercress shoot tissue biomass, tissue %N, and tissue C/N ratios were influenced by N but were unaffected by changes in S concentrations or by the interaction of NxS. Tissue pigment concentrations of beta-carotene, lutein, 5,6-epoxylutein, neoxanthin, zeaxanthin, and the chlorophyll pigments responded to changes in N treatment concentrations but were unaffected by S concentrations or through N x S interactions. Watercress tissue concentrations of aromatic, indole, and total glucosinolate concentrations responded to changes in N treatments; whereas aliphatic, aromatic, and total glucosinolates responded to changes in S treatment concentrations. Individual glucosinolates of glucobrassicin, 4-methoxyglucobrassicin, and gluconasturriin responded to N fertility treatments, while gluconapin, glucobrassicin, and gluconasturiin responded to changes in S fertility concentrations. Increases in carotenoid and glucosinolate concentrations through fertility management would be expected to influence the nutritional value of watercress in human diets.