SUMMARYNitrate acts as a potent signal to control global gene expression in Arabidopsis. Using an integrative bioinformatics approach we identified TGA1 and TGA4 as putative regulatory factors that mediate nitrate responses in Arabidopsis roots. We showed that both TGA1 and TGA4 mRNAs accumulate strongly after nitrate treatments in roots. Global gene expression analysis revealed 97% of the genes with altered expression in tga1 tga4 double mutant plants respond to nitrate treatments, indicating that these transcription factors have a specific role in nitrate responses in Arabidopsis root organs. We found TGA1 and TGA4 regulate the expression of nitrate transporter genes NRT2.1 and NRT2.2. Specific binding of TGA1 to its cognate DNA sequence on NRT2.1 and NRT2.2 promoters was confirmed by chromatin immunoprecipitation assays. The tga1 tga4 double mutant plants exhibit nitrate-dependent lateral and primary root phenotypes. Lateral root initiation is affected in both tga1 tga4 and nrt1.2 nrt2.2 double mutants, suggesting TGA1 and TGA4 regulate lateral root development at least partly via NRT2.1 and NRT2.2. Additional root phenotypes of tga1 tga4 double mutants indicate that these transcription factors play an important role in root developmental responses to nitrate. These results identify TGA1 and TGA4 as important regulatory factors of the nitrate response in Arabidopsis roots.

Auxin is a key phytohormone regulating central processes in plants. Although the mechanism by which auxin triggers changes in gene expression is well understood, little is known about the specific role of the individual members of the TIR1/AFB auxin receptors, Aux/IAA repressors, and ARF transcription factors and/or molecular pathways acting downstream leading to plant responses to the environment. We previously reported a role for AFB3 in coordinating primary and lateral root growth to nitrate availability. In this work, we used an integrated genomics, bioinformatics, and molecular genetics approach to dissect regulatory networks acting downstream of AFB3 that are activated by nitrate in roots. We found that the NAC4 transcription factor is a key regulatory element controlling a nitrate-responsive network, and that nac4 mutants have altered lateral root growth but normal primary root growth in response to nitrate. This finding suggests that AFB3 is able to activate two independent pathways to control root system architecture. Our systems approach has unraveled key components of the AFB3 regulatory network leading to changes in lateral root growth in response to nitrate.systems biology | nitrogen | Sungear

Root hairs are specialized cells that are important for nutrient uptake. It is well established that nutrients such as phosphate have a great influence on root hair development in many plant species. Here we investigated the role of nitrate on root hair development at a physiological and molecular level. We showed that nitrate increases root hair density in Arabidopsis thaliana. We found that two different root hair defective mutants have significantly less nitrate than wild-type plants, suggesting that in A. thaliana root hairs have an important role in the capacity to acquire nitrate. Nitrate reductase-null mutants exhibited nitrate-dependent root hair phenotypes comparable with wild-type plants, indicating that nitrate is the signal that leads to increased formation of root hairs. We examined the role of two key regulators of root hair cell fate, CPC and WER, in response to nitrate treatments. Phenotypic analyses of these mutants showed that CPC is essential for nitrate-induced responses of root hair development. Moreover, we showed that NRT1.1 and TGA1/TGA4 are required for pathways that induce root hair development by suppression of longitudinal elongation of trichoblast cells in response to nitrate treatments. Our results prompted a model where nitrate signaling via TGA1/TGA4 directly regulates the CPC root hair cell fate specification gene to increase formation of root hairs in A. thaliana.

The Atacama Desert in Chile—hyperarid and with high–ultraviolet irradiance levels—is one of the harshest environments on Earth. Yet, dozens of species grow there, including Atacama-endemic plants. Herein, we establish the Talabre–Lejía transect (TLT) in the Atacama as an unparalleled natural laboratory to study plant adaptation to extreme environmental conditions. We characterized climate, soil, plant, and soil–microbe diversity at 22 sites (every 100 m of altitude) along the TLT over a 10-y period. We quantified drought, nutrient deficiencies, large diurnal temperature oscillations, and pH gradients that define three distinct vegetational belts along the altitudinal cline. We deep-sequenced transcriptomes of 32 dominant plant species spanning the major plant clades, and assessed soil microbes by metabarcoding sequencing. The top-expressed genes in the 32 Atacama species are enriched in stress responses, metabolism, and energy production. Moreover, their root-associated soils are enriched in growth-promoting bacteria, including nitrogen fixers. To identify genes associated with plant adaptation to harsh environments, we compared 32 Atacama species with the 32 closest sequenced species, comprising 70 taxa and 1,686,950 proteins. To perform phylogenomic reconstruction, we concatenated 15,972 ortholog groups into a supermatrix of 8,599,764 amino acids. Using two codon-based methods, we identified 265 candidate positively selected genes (PSGs) in the Atacama plants, 64% of which are located in Pfam domains, supporting their functional relevance. For 59/184 PSGs with an Arabidopsis ortholog, we uncovered functional evidence linking them to plant resilience. As some Atacama plants are closely related to staple crops, these candidate PSGs are a “genetic goldmine” to engineer crop resilience to face climate change.

The rapid increase in the availability of transcriptomics data generated by RNA sequencing represents both a challenge and an opportunity for biologists without bioinformatics training. The challenge is handling, integrating, and interpreting these data sets. The opportunity is to use this information to generate testable hypothesis to understand molecular mechanisms controlling gene expression and biological processes (Fig. 1). A successful strategy to generate tractable hypotheses from transcriptomics data has been to build undirected network graphs based on patterns of gene co-expression. Many examples of new hypothesis derived from network analyses can be found in the literature, spanning different organisms including plants and specific fields such as root developmental biology.In order to make the process of constructing a gene co-expression network more accessible to biologists, here we provide step-by-step instructions using published RNA-seq experimental data obtained from a public database. Similar strategies have been used in previous studies to advance root developmental biology. This guide includes basic instructions for the operation of widely used open source platforms such as Bio-Linux, R, and Cytoscape. Even though the data we used in this example was obtained from Arabidopsis thaliana, the workflow developed in this guide can be easily adapted to work with RNA-seq data from any organism.

Nitrate is a nutrient and a potent signal that impacts global gene expression in plants. However, the regulatory factors controlling temporal and cell type–specific nitrate responses remain largely unknown. We assayed nitrate-responsive transcriptome changes in five major root cell types of the Arabidopsis thaliana root as a function of time. We found that gene-expression response to nitrate is dynamic and highly localized and predicted cell type–specific transcription factor (TF)–target interactions. Among cell types, the endodermis stands out as having the largest and most connected nitrate-regulatory gene network. ABF2 and ABF3 are major hubs for transcriptional responses in the endodermis cell layer. We experimentally validated TF–target interactions for ABF2 and ABF3 by chromatin immunoprecipitation followed by sequencing and a cell-based system to detect TF regulation genome-wide. Validated targets of ABF2 and ABF3 account for more than 50% of the nitrate-responsive transcriptome in the endodermis. Moreover, ABF2 and ABF3 are involved in nitrate-induced lateral root growth. Our approach offers an unprecedented spatiotemporal resolution of the root response to nitrate and identifies important components of cell-specific gene regulatory networks.

Estimating total plant diversity in extreme or hyperarid environments can be challenging, as adaptations to pronounced climate variability include evading prolonged stress periods through seeds or specialized underground organs. Short-term surveys of these ecosystems are thus likely poor estimators of actual diversity. Here we develop a multimethod strategy to obtain a more complete understanding of plant diversity from a community in the Atacama Desert. We explicitly test environmental DNA-based techniques (eDNA) to see if they can reveal the observed and 'hidden' (dormant or locally rare) species. To estimate total plant diversity, we performed long-term traditional surveys during eight consecutive years, including El Niño and La Niña events, we then analyzed eDNA from soil samples using high-throughput sequencing. We further used soil pollen analysis and soil seed bank germination assays to identify 'hidden' species. Each approach offers different subsets of current biodiversity at different taxonomic, spatial and temporal resolution, with a total of 92 taxa identified along the transect. Traditional field surveys identified 77 plant species over eight consecutive years. Observed community composition greatly varies interannually, with only 22 species seen every year. eDNA analysis revealed 37 taxa, eight of which were 'hidden' in our field surveys. Soil samples contain a viable seed bank of 21 taxa. Soil pollen (27 taxa) and eDNA analysis show affinities with vegetation at the landscape scale but a weak relationship to local plot diversity. Multimethod approaches (including eDNA) in deserts are valuable tools that add to a comprehensive assessment of biodiversity in such extreme environments, where using a single method or observations over a few years is insufficient. Our results can also explain the resilience of Atacama plant communities as 'hidden' taxa may have been active in the recent past or could even emerge in the future as accelerated global environmental change continues unabated.

With long-read sequencing we have entered an era where individual genomes are routinely assembled to near-completion and where complex genetic variation can efficiently be resolved. Here we demonstrate that long reads can be applied also to study the genomic architecture of individual human cells. Clonally expanded CD8+ T-cells from a human donor were used as starting material for a droplet-based multiple displacement amplification (dMDA) method designed to ensure long molecule lengths and minimal amplification bias. Sequencing of two single cells was performed on the PacBio Sequel II system, generating over 2.5 million reads and ~20Gb HiFi data (>QV20) per cell, achieving up to 40% genome coverage. This data allowed for single nucleotide variant (SNV) detection, including in genomic regions inaccessible by short reads. Over 1000 high-confidence structural variants (SVs) per cell were discovered in the PacBio data, which is four times more than the number of SVs detected in Illumina dMDA data from clonally related cells. In addition, several putative clone-specific somatic SV events could be identified. Single-cell de novo assembly resulted in 454-598 Mb assembly sizes and 35-42 kb contig N50 values. 1762 (12.8%) of expected gene models were found to be complete in the best single-cell assembly. The de novo constructed mitochondrial genomes were 100% identical for the two single cells subjected to PacBio sequencing, although mitochondrial heteroplasmy was also observed. In summary, the work presented here demonstrates the utility of long-read sequencing towards understanding the extent and distribution of complex genetic variation at the single cell level.