Allelic variants confer Arabidopsis adaptation to small regional environmental differences.

Natural populations of Arabidopsis thaliana provide powerful systems to study the adaptation of wild plant species. Previous research has predominantly focused on global populations or accessions collected from regions with diverse climates. However, little is known about the genetics underlying adaptation in regions with mild environmental clines. We have examined a diversity panel consisting of 192 A. thaliana accessions collected from the Netherlands, a region with limited climatic variation. Despite the relatively uniform climate, we identified evidence of local adaptation within this population. Notably, semidwarf accessions, due to mutation of the GIBBERELLIC ACID REQUIRING 5 (GA5) gene, occur at a relatively high frequency near the coast and these displayed enhanced tolerance to high wind velocities. Additionally, we evaluated the performance of the population under iron deficiency conditions and found that allelic variation in the FE SUPEROXIDE DISMUTASE 3 (FSD3) gene affects tolerance to low iron levels. Moreover, we explored patterns of local adaptation to environmental clines in temperature and precipitation, observing that allelic variation at LA RELATED PROTEIN 1C (LARP1c) likely affects drought tolerance. Not only is the genetic variation observed in a diversity panel of A. thaliana collected in a region with mild environmental clines comparable to that in collections sampled over larger geographic ranges but it is also sufficiently rich to elucidate the genetic and environmental factors underlying natural plant adaptation.

Lactuca super-pangenome reduces bias towards reference genes in lettuce research.

BACKGROUND: Breeding of lettuce (Lactuca sativa L.), the most important leafy vegetable worldwide, for enhanced disease resistance and resilience relies on multiple wild relatives to provide the necessary genetic diversity. In this study, we constructed a super-pangenome based on four Lactuca species (representing the primary, secondary and tertiary gene pools) and comprising 474 accessions. We include 68 newly sequenced accessions to improve cultivar coverage and add important foundational breeding lines. RESULTS: With the super-pangenome we find substantial presence/absence variation (PAV) and copy-number variation (CNV). Functional enrichment analyses of core and variable genes show that transcriptional regulators are conserved whereas disease resistance genes are variable. PAV-genome-wide association studies (GWAS) and CNV-GWAS are largely congruent with single-nucleotide polymorphism (SNP)-GWAS. Importantly, they also identify several major novel quantitative trait loci (QTL) for resistance against Bremia lactucae in variable regions not present in the reference lettuce genome. The usability of the super-pangenome is demonstrated by identifying the likely origin of non-reference resistance loci from the wild relatives Lactuca serriola, Lactuca saligna and Lactuca virosa. CONCLUSIONS: The super-pangenome offers a broader view on the gene repertoire of lettuce, revealing relevant loci that are not in the reference genome(s). The provided methodology and data provide a strong basis for research into PAVs, CNVs and other variation underlying important biological traits of lettuce and other crops.

Architecture and dynamics of the abscisic acid gene regulatory network.

The plant hormone abscisic acid (ABA) regulates essential processes in plant development and responsiveness to abiotic and biotic stresses. ABA perception triggers a post-translational signaling cascade that elicits the ABA gene regulatory network (GRN), encompassing hundreds of transcription factors (TFs) and thousands of transcribed genes. To further our knowledge of this GRN, we performed an RNA-seq time series experiment consisting of 14 time points in the 16 h following a one-time ABA treatment of 5-week-old Arabidopsis rosettes. During this time course, ABA rapidly changed transcription levels of 7151 genes, which were partitioned into 44 coexpressed modules that carry out diverse biological functions. We integrated our time-series data with publicly available TF-binding site data, motif data, and RNA-seq data of plants inhibited in translation, and predicted (i) which TFs regulate the different coexpression clusters, (ii) which TFs contribute the most to target gene amplitude, (iii) timing of engagement of different TFs in the ABA GRN, and (iv) hierarchical position of TFs and their targets in the multi-tiered ABA GRN. The ABA GRN was found to be highly interconnected and regulated at different amplitudes and timing by a wide variety of TFs, of which the bZIP family was most prominent, and upregulation of genes encompassed more TFs than downregulation. We validated our network models in silico with additional public TF-binding site data and transcription data of selected TF mutants. Finally, using a drought assay we found that the Trihelix TF GT3a is likely an ABA-induced positive regulator of drought tolerance.

Integrating taxonomic signals from MAGs and contigs improves read annotation and taxonomic profiling of metagenomes.

Metagenomic analysis typically includes read-based taxonomic profiling, assembly, and binning of metagenome-assembled genomes (MAGs). Here we integrate these steps in Read Annotation Tool (RAT), which uses robust taxonomic signals from MAGs and contigs to enhance read annotation. RAT reconstructs taxonomic profiles with high precision and sensitivity, outperforming other state-of-the-art tools. In high-diversity groundwater samples, RAT annotates a large fraction of the metagenomic reads, calling novel taxa at the appropriate, sometimes high taxonomic ranks. Thus, RAT integrative profiling provides an accurate and comprehensive view of the microbiome from shotgun metagenomics data. The package of Contig Annotation Tool (CAT), Bin Annotation Tool (BAT), and RAT is available at https://github.com/MGXlab/CAT_pack (from CAT pack v6.0). The CAT pack now also supports Genome Taxonomy Database (GTDB) annotations.

A low-cost open-source imaging platform reveals spatiotemporal insight into leaf elongation and movement.

Plant organs move throughout the diurnal cycle, changing leaf and petiole positions to balance light capture, leaf temperature, and water loss under dynamic environmental conditions. Upward movement of the petiole, called hyponasty, is one of several traits of the shade avoidance syndrome (SAS). SAS traits are elicited upon perception of vegetation shade signals such as far-red light (FR) and improve light capture in dense vegetation. Monitoring plant movement at a high temporal resolution allows studying functionality and molecular regulation of hyponasty. However, high temporal resolution imaging solutions are often very expensive, making this unavailable to many researchers. Here, we present a modular and low-cost imaging setup, based on small Raspberry Pi computers that can track leaf movements and elongation growth with high temporal resolution. We also developed an open-source, semiautomated image analysis pipeline. Using this setup, we followed responses to FR enrichment, light intensity, and their interactions. Tracking both elongation and the angle of the petiole, lamina, and entire leaf in Arabidopsis (Arabidopsis thaliana) revealed insight into R:FR sensitivities of leaf growth and movement dynamics and the interactions of R:FR with background light intensity. The detailed imaging options of this system allowed us to identify spatially separate bending points for petiole and lamina positioning of the leaf.

The effect of developmental variation on expression QTLs in a multi parental Caenorhabditis elegans population.

Regulation of gene expression plays a crucial role in developmental processes and adaptation to changing environments. expression quantitative trait locus (eQTL) mapping is a technique used to study the genetic regulation of gene expression using the transcriptomes of recombinant inbred lines (RILs). Typically, the age of the inbred lines at the time of RNA sampling is carefully controlled. This is necessary because the developmental process causes changes in gene expression, complicating the interpretation of eQTL mapping experiments. However, due to genetics and variation in ambient micro-environments, organisms can differ in their "developmental age," even if they are of the same chronological age. As a result, eQTL patterns are affected by developmental variation in gene expression. The model organism Caenorhabditis elegans is particularly suited for studying the effect of developmental variation on eQTL mapping patterns. In a span of days, C. elegans transitions from embryo through 4 larval stages to adult while undergoing massive changes to its transcriptome. Here, we use C. elegans to investigate the effect of developmental age variation on eQTL patterns and present a normalization procedure. We used dynamical eQTL mapping, which includes the developmental age as a cofactor, to separate the variation in development from genotypic variation and explain variation in gene expression levels. We compare classical single marker eQTL mapping and dynamical eQTL mapping using RNA-seq data of ∼200 multi-parental RILs of C. elegans. The results show that (1) many eQTLs are caused by developmental variation, (2) most trans-bands are developmental QTLs, and (3) dynamical eQTL mapping detects additional eQTLs not found with classical eQTL mapping. We recommend that correction for variation in developmental age should be strongly considered in eQTL mapping studies given the large impact of processes like development on the transcriptome.

Cryptic genetic variation of expression quantitative trait locus architecture revealed by genetic perturbation in Caenorhabditis elegans.

Genetic perturbation in different genetic backgrounds can cause a range of phenotypes within a species. These phenotypic differences can be the result of the interaction between the genetic background and the perturbation. Previously, we reported that perturbation of gld-1, an important player in the developmental control of Caenorhabditis elegans, released cryptic genetic variation (CGV) affecting fitness in different genetic backgrounds. Here, we investigated the change in transcriptional architecture. We found 414 genes with a cis-expression quantitative trait locus (eQTL) and 991 genes with a trans-eQTL that were specifically found in the gld-1 RNAi treatment. In total, we detected 16 eQTL hotspots, of which 7 were only found in the gld-1 RNAi treatment. Enrichment analysis of those 7 hotspots showed that the regulated genes were associated with neurons and the pharynx. Furthermore, we found evidence of accelerated transcriptional aging in the gld-1 RNAi-treated nematodes. Overall, our results illustrate that studying CGV leads to the discovery of hidden polymorphic regulators.

Plasticity of maternal environment-dependent expression-QTLs of tomato seeds.

Seeds are essential for plant reproduction, survival, and dispersal. Germination ability and successful establishment of young seedlings strongly depend on seed quality and on environmental factors such as nutrient availability. In tomato (Solanum lycopersicum) and many other species, seed quality and seedling establishment characteristics are determined by genetic variation, as well as the maternal environment in which the seeds develop and mature. The genetic contribution to variation in seed and seedling quality traits and environmental responsiveness can be estimated at transcriptome level in the dry seed by mapping genomic loci that affect gene expression (expression QTLs) in contrasting maternal environments. In this study, we applied RNA-sequencing to construct a linkage map and measure gene expression of seeds of a tomato recombinant inbred line (RIL) population derived from a cross between S. lycopersicum (cv. Moneymaker) and S. pimpinellifolium (G1.1554). The seeds matured on plants cultivated under different nutritional environments, i.e., on high phosphorus or low nitrogen. The obtained single-nucleotide polymorphisms (SNPs) were subsequently used to construct a genetic map. We show how the genetic landscape of plasticity in gene regulation in dry seeds is affected by the maternal nutrient environment. The combined information on natural genetic variation mediating (variation in) responsiveness to the environment may contribute to knowledge-based breeding programs aiming to develop crop cultivars that are resilient to stressful environments.

Local light signaling at the leaf tip drives remote differential petiole growth through auxin-gibberellin dynamics.

Although plants are immobile, many of their organs are flexible to move in response to environmental cues. In dense vegetation, plants detect neighbors through far-red light perception with their leaf tip. They respond remotely, with asymmetrical growth between the abaxial and adaxial sides of the leafstalk, the petiole. This results in upward movement that brings the leaf blades into better lit zones of the canopy. The plant hormone auxin is required for this response, but it is not understood how non-differential leaf tip-derived auxin can remotely regulate movement. Here, we show that remote signaling of far-red light promotes auxin accumulation in the abaxial petiole. This local auxin accumulation is facilitated by reinforcing an intrinsic directionality of the auxin transport protein PIN3 on the petiole endodermis, as visualized with a PIN3-GFP line. Using an auxin biosensor, we show that auxin accumulates in all cell layers from endodermis to epidermis in the abaxial petiole, upon far-red light signaling in the remote leaf tip. In the petiole, auxin elicits a response to both auxin itself as well as a second growth promoter; gibberellin. We show that this dual regulation is necessary for hyponastic leaf movement in response to light. Our data indicate that gibberellin is required to permit cell growth, whereas differential auxin accumulation determines which cells can grow. Our results reveal how plants can spatially relay information about neighbor proximity from their sensory leaf tips to the petiole base, thus driving adaptive growth.

Coumarin biosynthesis genes are required after foliar pathogen infection for the creation of a microbial soil-borne legacy that primes plants for SA-dependent defenses.

Plants deposit photosynthetically-fixed carbon in the rhizosphere, the thin soil layer directly around the root, thereby creating a hospitable environment for microbes. To manage the inhabitants of this nutrient-rich environment, plant roots exude and dynamically adjust microbe-attracting and -repelling compounds to stimulate specific members of the microbiome. Previously, we demonstrated that foliar infection of Arabidopsis thaliana by the biotrophic downy mildew pathogen Hyaloperonospora arabidopsidis (Hpa) leads to a disease-induced modification of the rhizosphere microbiome. Soil conditioned with Hpa-infected plants provided enhanced protection against foliar downy mildew infection in a subsequent population of plants, a phenomenon dubbed the soil-borne legacy (SBL). Here, we show that for the creation of the SBL, plant-produced coumarins play a prominent role as coumarin-deficient myb72 and f6'h1 mutants were defective in creating a Hpa-induced SBL. Root exudation profiles changed significantly in Col-0 upon foliar Hpa infection, and this was accompanied by a compositional shift in the root microbiome that was significantly different from microbial shifts occurring on roots of Hpa-infected coumarin-deficient mutants. Our data further show that the Hpa-induced SBL primes Col-0 plants growing in SBL-conditioned soil for salicylic acid (SA)-dependent defenses. The SA-signaling mutants sid2 and npr1 were unresponsive to the Hpa-induced SBL, suggesting that the protective effect of the Hpa-induced shift in the root microbiome results from an induced systemic resistance that requires SA-signaling in the plant.

A metagenomic portrait of the microbial community responsible for two decades of bioremediation of poly-contaminated groundwater.

Biodegradation of pollutants is a sustainable and cost-effective solution to groundwater pollution. Here, we investigate microbial populations involved in biodegradation of poly-contaminants in a pipeline for heavily contaminated groundwater. Groundwater moves from a polluted park to a treatment plant, where an aerated bioreactor effectively removes the contaminants. While the biomass does not settle in the reactor, sediment is collected afterwards and used to seed the new polluted groundwater via a backwash cycle. The pipeline has successfully operated since 1999, but the biological components in the reactor and the contaminated park groundwater have never been described. We sampled seven points along the pipeline, representing the entire remediation process, and characterized the changing microbial communities using genome-resolved metagenomic analysis. We assembled 297 medium- and high-quality metagenome-assembled genome sequences representing on average 46.3% of the total DNA per sample. We found that the communities cluster into two distinct groups, separating the anaerobic communities in the park groundwater from the aerobic communities inside the plant. In the park, the community is dominated by members of the genus Sulfuricurvum, while the plant is dominated by generalists from the order Burkholderiales. Known aromatic compound biodegradation pathways are four times more abundant in the plant-side communities compared to the park-side. Our findings provide a genome-resolved portrait of the microbial community in a highly effective groundwater treatment system that has treated groundwater with a complex contamination profile for two decades.

Disentangling the genetic basis of rhizosphere microbiome assembly in tomato.

Microbiomes play a pivotal role in plant growth and health, but the genetic factors involved in microbiome assembly remain largely elusive. Here, we map the molecular features of the rhizosphere microbiome as quantitative traits of a diverse hybrid population of wild and domesticated tomato. Gene content analysis of prioritized tomato quantitative trait loci suggests a genetic basis for differential recruitment of various rhizobacterial lineages, including a Streptomyces-associated 6.31 Mbp region harboring tomato domestication sweeps and encoding, among others, the iron regulator FIT and the water channel aquaporin SlTIP2.3. Within metagenome-assembled genomes of root-associated Streptomyces and Cellvibrio, we identify bacterial genes involved in metabolism of plant polysaccharides, iron, sulfur, trehalose, and vitamins, whose genetic variation associates with specific tomato QTLs. By integrating 'microbiomics' and quantitative plant genetics, we pinpoint putative plant and reciprocal rhizobacterial traits underlying microbiome assembly, thereby providing a first step towards plant-microbiome breeding programs.

The genetic architecture underlying body-size traits plasticity over different temperatures and developmental stages in Caenorhabditis elegans.

Most ectotherms obey the temperature-size rule, meaning they grow larger in a colder environment. This raises the question of how the interplay between genes and temperature affects the body size of ectotherms. Despite the growing body of literature on the physiological life-history and molecular genetic mechanism underlying the temperature-size rule, the overall genetic architecture orchestrating this complex phenotype is not yet fully understood. One approach to identify genetic regulators of complex phenotypes is quantitative trait locus (QTL) mapping. Here, we explore the genetic architecture of body-size phenotypes, and plasticity of body-size phenotypes at different temperatures using Caenorhabditis elegans as a model ectotherm. We used 40 recombinant inbred lines (RILs) derived from N2 and CB4856, which were reared at four different temperatures (16, 20, 24, and 26 °C) and measured at two developmental stages (L4 and adult). The animals were measured for body length, width at vulva, body volume, length/width ratio, and seven other body-size traits. The genetically diverse RILs varied in their body-size phenotypes with heritabilities ranging from 0.0 to 0.99. We detected 18 QTL underlying the body-size traits across all treatment combinations, with the majority clustering on Chromosome X. We hypothesize that the Chromosome X QTL could result from a known pleiotropic regulator-npr-1-known to affect the body size of C. elegans through behavioral changes. We also found five plasticity QTL of body-size traits where three colocalized with body-size QTL. In conclusion, our findings shed more light on multiple loci affecting body-size plasticity and the possibility of co-regulation of traits and traits plasticity by the same loci under different environments.

Embodiment in distributed information processing: "Solid" plants versus "liquid" ant colonies.

Information processing is an essential part of biology, enabling coordination of intra-organismal processes such as development, environmental adaptation and inter-organismal communication. Whilst in animals with specialised brain tissue a substantial amount of information processing occurs in a centralised manner, most biological computing is distributed across multiple entities, such as cells in a tissue, roots in a root system or ants in a colony. Physical context, called embodiment, also affects the nature of biological computing. While plants and ant colonies both perform distributed computing, in plants the units occupy fixed positions while individual ants move around. This distinction, solid versus liquid brain computing, shapes the nature of computations. Here we compare information processing in plants and ant colonies, highlighting how similarities and differences originate in, as well as make use of, the differences in embodiment. We end with a discussion on how this embodiment perspective may inform the debate on plant cognition.

2D morphometric analysis of Arabidopsis thaliana nuclei reveals characteristic profiles of different cell types and accessions.

Functional changes of cells upon developmental switches and in response to environmental cues are often reflected in nuclear phenotypes, showing distinctive chromatin states corresponding to transcriptional changes. Such characteristic nuclear shapes have been microscopically monitored and can be quantified after differential staining of euchromatin and heterochromatin domains. Here, we examined several nuclear parameters (size, DNA content, DNA density, chromatin compaction, relative heterochromatin fraction (RHF), and number of chromocenters) in relation to spatial distribution of genes and transposon elements (TEs), using standard 2D fluorescence microscopy. We provide nuclear profiles for different cell types and different accessions of Arabidopsis thaliana. A variable, yet significant, fraction of TEs was found outside chromocenters in all cell types, except for guard cells. The latter cell type features nuclei with the highest level of chromatin compaction, while their chromocenters seem to contain gene-rich regions. The highest number of parameter correlations was found in the accession Cvi, whereas Ler showed only few correlations. This may point at differences in phenotype robustness between accessions. The significantly high association of NOR chromocenters in accessions Ws and Cvi corresponds to their low RHF level.

Mechanisms of far-red light-mediated dampening of defense against Botrytis cinerea in tomato leaves.

Plants detect neighboring competitors through a decrease in the ratio between red and far-red light (R:FR). This decreased R:FR is perceived by phytochrome photoreceptors and triggers shade avoidance responses such as shoot elongation and upward leaf movement (hyponasty). In addition to promoting elongation growth, low R:FR perception enhances plant susceptibility to pathogens: the growth-defense tradeoff. Although increased susceptibility in low R:FR has been studied for over a decade, the associated timing of molecular events is still unknown. Here, we studied the chronology of FR-induced susceptibility events in tomato (Solanum lycopersicum) plants pre-exposed to either white light (WL) or WL supplemented with FR light (WL+FR) prior to inoculation with the necrotrophic fungus Botrytis cinerea (B.c.). We monitored the leaf transcriptional changes over a 30-h time course upon infection and followed up with functional studies to identify mechanisms. We found that FR-induced susceptibility in tomato is linked to a general dampening of B.c.-responsive gene expression, and a delay in both pathogen recognition and jasmonic acid-mediated defense gene expression. In addition, we found that the supplemental FR-induced ethylene emissions affected plant immune responses under the WL+FR condition. This study improves our understanding of the growth-immunity tradeoff, while simultaneously providing leads to improve tomato resistance against pathogens in dense cropping systems.

The genetics of gene expression in a Caenorhabditis elegans multiparental recombinant inbred line population.

Studying genetic variation of gene expression provides a powerful way to unravel the molecular components underlying complex traits. Expression quantitative trait locus (eQTL) studies have been performed in several different model species, yet most of these linkage studies have been based on the genetic segregation of two parental alleles. Recently, we developed a multiparental segregating population of 200 recombinant inbred lines (mpRILs) derived from four wild isolates (JU1511, JU1926, JU1931, and JU1941) in the nematode Caenorhabditis elegans. We used RNA-seq to investigate how multiple alleles affect gene expression in these mpRILs. We found 1789 genes differentially expressed between the parental lines. Transgression, expression beyond any of the parental lines in the mpRILs, was found for 7896 genes. For expression QTL mapping almost 9000 SNPs were available. By combining these SNPs and the RNA-seq profiles of the mpRILs, we detected almost 6800 eQTLs. Most trans-eQTLs (63%) co-locate in six newly identified trans-bands. The trans-eQTLs found in previous two-parental allele eQTL experiments and this study showed some overlap (17.5-46.8%), highlighting on the one hand that a large group of genes is affected by polymorphic regulators across populations and conditions, on the other hand, it shows that the mpRIL population allows identification of novel gene expression regulatory loci. Taken together, the analysis of our mpRIL population provides a more refined insight into C. elegans complex trait genetics and eQTLs in general, as well as a starting point to further test and develop advanced statistical models for detection of multiallelic eQTLs and systems genetics studying the genotype-phenotype relationship.

Transcriptome Signatures in Pseudomonas simiae WCS417 Shed Light on Role of Root-Secreted Coumarins in Arabidopsis-Mutualist Communication.

Pseudomonas simiae WCS417 is a root-colonizing bacterium with well-established plant-beneficial effects. Upon colonization of Arabidopsis roots, WCS417 evades local root immune responses while triggering an induced systemic resistance (ISR) in the leaves. The early onset of ISR in roots shows similarities with the iron deficiency response, as both responses are associated with the production and secretion of coumarins. Coumarins can mobilize iron from the soil environment and have a selective antimicrobial activity that impacts microbiome assembly in the rhizosphere. Being highly coumarin-tolerant, WCS417 induces the secretion of these phenolic compounds, likely to improve its own niche establishment, while providing growth and immunity benefits for the host in return. To investigate the possible signaling function of coumarins in the mutualistic Arabidopsis-WCS417 interaction, we analyzed the transcriptome of WCS417 growing in root exudates of coumarin-producing Arabidopsis Col-0 and the coumarin-biosynthesis mutant f6'h1. We found that coumarins in F6'H1-dependent root exudates significantly affected the expression of 439 bacterial genes (8% of the bacterial genome). Of those, genes with functions related to transport and metabolism of carbohydrates, amino acids, and nucleotides were induced, whereas genes with functions related to cell motility, the bacterial mobilome, and energy production and conversion were repressed. Strikingly, most genes related to flagellar biosynthesis were down-regulated by F6'H1-dependent root exudates and we found that application of selected coumarins reduces bacterial motility. These findings suggest that coumarins' function in the rhizosphere as semiochemicals in the communication between the roots and WCS417. Collectively, our results provide important novel leads for future functional analysis of molecular processes in the establishment of plant-mutualist interactions.

The chemical compound 'Heatin' stimulates hypocotyl elongation and interferes with the Arabidopsis NIT1-subfamily of nitrilases.

Temperature passively affects biological processes involved in plant growth. Therefore, it is challenging to study the dedicated temperature signalling pathways that orchestrate thermomorphogenesis, a suite of elongation growth-based adaptations that enhance leaf-cooling capacity. We screened a chemical library for compounds that restored hypocotyl elongation in the pif4-2-deficient mutant background at warm temperature conditions in Arabidopsis thaliana to identify modulators of thermomorphogenesis. The small aromatic compound 'Heatin', containing 1-iminomethyl-2-naphthol as a pharmacophore, was selected as an enhancer of elongation growth. We show that ARABIDOPSIS ALDEHYDE OXIDASES redundantly contribute to Heatin-mediated hypocotyl elongation. Following a chemical proteomics approach, the members of the NITRILASE1-subfamily of auxin biosynthesis enzymes were identified among the molecular targets of Heatin. Our data reveal that nitrilases are involved in promotion of hypocotyl elongation in response to high temperature and Heatin-mediated hypocotyl elongation requires the NITRILASE1-subfamily members, NIT1 and NIT2. Heatin inhibits NIT1-subfamily enzymatic activity in vitro and the application of Heatin accordingly results in the accumulation of NIT1-subfamily substrate indole-3-acetonitrile in vivo. However, levels of the NIT1-subfamily product, bioactive auxin (indole-3-acetic acid), were also significantly increased. It is likely that the stimulation of hypocotyl elongation by Heatin might be independent of its observed interaction with NITRILASE1-subfamily members. However, nitrilases may contribute to the Heatin response by stimulating indole-3-acetic acid biosynthesis in an indirect way. Heatin and its functional analogues present novel chemical entities for studying auxin biology.