The transcriptional regulatory structure of plant genomes remains poorly defined relative to animals. It is unclear how many -regulatory elements exist, where these elements lie relative to promoters, and how these features are conserved across plant species. We employed the assay for transposase-accessible chromatin (ATAC-seq) in four plant species (, ,, and ) to delineate open chromatin regions and transcription factor (TF) binding sites across each genome. Despite 10-fold variation in intergenic space among species, the majority of open chromatin regions lie within 3 kb upstream of a transcription start site in all species. We find a common set of four TFs that appear to regulate conserved gene sets in the root tips of all four species, suggesting that TF-gene networks are generally conserved. Comparative ATAC-seq profiling of Arabidopsis root hair and non-hair cell types revealed extensive similarity as well as many cell-type-specific differences. Analyzing TF binding sites in differentially accessible regions identified a MYB-driven regulatory module unique to the hair cell, which appears to control both cell fate regulators and abiotic stress responses. Our analyses revealed common regulatory principles among species and shed light on the mechanisms producing cell-type-specific transcriptomes during development.

ORCID ID: 0000-0002-8568-7125 (J.B-S.).Agrobacterium rhizogenes (or Rhizobium rhizogenes) is able to transform plant genomes and induce the production of hairy roots. We describe the use of A. rhizogenes in tomato (Solanum spp.) to rapidly assess gene expression and function. Gene expression of reporters is indistinguishable in plants transformed by Agrobacterium tumefaciens as compared with A. rhizogenes. A root cell type-and tissue-specific promoter resource has been generated for domesticated and wild tomato (Solanum lycopersicum and Solanum pennellii, respectively) using these approaches. Imaging of tomato roots using A. rhizogenes coupled with laser scanning confocal microscopy is facilitated by the use of a membrane-tagged protein fused to a red fluorescent protein marker present in binary vectors. Tomatooptimized isolation of nuclei tagged in specific cell types and translating ribosome affinity purification binary vectors were generated and used to monitor associated messenger RNA abundance or chromatin modification. Finally, transcriptional reporters, translational reporters, and clustered regularly interspaced short palindromic repeats-associated nuclease9 genome editing demonstrate that SHORT-ROOT and SCARECROW gene function is conserved between Arabidopsis (Arabidopsis thaliana) and tomato.

Parasitic flowering plants are one of the most destructive agricultural pests and have major impact on crop yields throughout the world. Being dependent on finding a host plant for growth, parasitic plants penetrate their host using specialized organs called haustoria. Haustoria establish vascular connections with the host, which enable the parasite to steal nutrients and water. The underlying molecular and developmental basis of parasitism by plants is largely unknown. In order to investigate the process of parasitism, RNAs from different stages (i.e. seed, seedling, vegetative strand, prehaustoria, haustoria, and flower) were used to de novo assemble and annotate the transcriptome of the obligate plant stem parasite dodder (Cuscuta pentagona). The assembled transcriptome was used to dissect transcriptional dynamics during dodder development and parasitism and identified key gene categories involved in the process of plant parasitism. Host plant infection is accompanied by increased expression of parasite genes underlying transport and transporter categories, response to stress and stimuli, as well as genes encoding enzymes involved in cell wall modifications. By contrast, expression of photosynthetic genes is decreased in the dodder infective stages compared with normal stem. In addition, genes relating to biosynthesis, transport, and response of phytohormones, such as auxin, gibberellins, and strigolactone, were differentially expressed in the dodder infective stages compared with stems and seedlings. This analysis sheds light on the transcriptional changes that accompany plant parasitism and will aid in identifying potential gene targets for use in controlling the infestation of crops by parasitic weeds.

Flooding due to extreme weather threatens crops and ecosystems. To understand variation in gene regulatory networks activated by submergence, we conducted a high-resolution analysis of chromatin accessibility and gene expression at three scales of transcript control in four angiosperms, ranging from a dryland-adapted wild species to a wetland crop. The data define a cohort of conserved submergence-activated genes with signatures of overlapping cis regulation by four transcription factor families. Syntenic genes are more highly expressed than nonsyntenic genes, yet both can have the cis motifs and chromatin accessibility associated with submergence up-regulation. Whereas the flexible circuitry spans the eudicot-monocot divide, the frequency of specific cis motifs, extent of chromatin accessibility, and degree of submergence activation are more prevalent in the wetland crop and may have adaptive importance.

Next Generation Sequencing (NGS) is driving rapid advancement in biological understanding and RNA-sequencing (RNA-seq) has become an indispensable tool for biology and medicine. There is a growing need for access to these technologies although preparation of NGS libraries remains a bottleneck to wider adoption. Here we report a novel method for the production of strand specific RNA-seq libraries utilizing the terminal breathing of double-stranded cDNA to capture and incorporate a sequencing adapter. Breath Adapter Directional sequencing (BrAD-seq) reduces sample handling and requires far fewer enzymatic steps than most available methods to produce high quality strand-specific RNA-seq libraries. The method we present is optimized for 3-prime Digital Gene Expression (DGE) libraries and can easily extend to full transcript coverage shotgun (SHO) type strand-specific libraries and is modularized to accommodate a diversity of RNA and DNA input materials. BrAD-seq offers a highly streamlined and inexpensive option for RNA-seq libraries.

Leaf shape is mutable, changing in ways modulated by both development and environment within genotypes. A complete model of leaf phenotype would incorporate the changes in leaf shape during juvenile-to-adult phase transitions and the ontogeny of each leaf. Here, we provide a morphometric description of >33,000 leaflets from a set of tomato (Solanum spp) introgression lines grown under controlled environment conditions. We first compare the shape of these leaves, arising during vegetative development, with >11,000 previously published leaflets from a field setting and >11,000 leaflets from wild tomato relatives. We then quantify the changes in shape, across ontogeny, for successive leaves in the heteroblastic series. Using principal component analysis, we then separate genetic effects modulating (1) the overall shape of all leaves versus (2) the shape of specific leaves in the series, finding the former more heritable than the latter and comparing quantitative trait loci regulating each. Our results demonstrate that phenotype is highly contextual and that unbiased assessments of phenotype, for quantitative genetic or other purposes, would ideally sample the many developmental and environmental factors that modulate it.

Variation in gene expression, in addition to sequence polymorphisms, is known to influence developmental, physiological, and metabolic traits in plants. Genetic mapping populations have facilitated identification of expression quantitative trait loci (eQTL), the genetic determinants of variation in gene expression patterns. We used an introgression population developed from the wild desertadapted Solanum pennellii and domesticated tomato (Solanum lycopersicum) to identify the genetic basis of transcript level variation. We established the effect of each introgression on the transcriptome and identified approximately 7,200 eQTL regulating the steady-state transcript levels of 5,300 genes. Barnes-Hut t-distributed stochastic neighbor embedding clustering identified 42 modules revealing novel associations between transcript level patterns and biological processes. The results showed a complex genetic architecture of global transcript abundance pattern in tomato. Several genetic hot spots regulating a large number of transcript level patterns relating to diverse biological processes such as plant defense and photosynthesis were identified. Important eQTL regulating transcript level patterns were related to leaf number and complexity as well as hypocotyl length. Genes associated with leaf development showed an inverse correlation with photosynthetic gene expression, but eQTL regulating genes associated with leaf development and photosynthesis were dispersed across the genome. This comprehensive eQTL analysis details the influence of these loci on plant phenotypes and will be a valuable community resource for investigations on the genetic effects of eQTL on phenotypic traits in tomato.

Plants sense the foliar shade of competitors and alter their developmental programs through the shade-avoidance response. Internode and petiole elongation, and changes in overall leaf area and leaf mass per area, are the stereotypical architectural responses to foliar shade in the shoot. However, changes in leaf shape and complexity in response to shade remain incompletely, and qualitatively, described. Using a meta-analysis of more than 18,000 previously published leaflet outlines, we demonstrate that shade avoidance alters leaf shape in domesticated tomato (Solanum lycopersicum) and wild relatives. The effects of shade avoidance on leaf shape are subtle with respect to individual traits but are combinatorially strong. We then seek to describe the developmental origins of shade-induced changes in leaf shape by swapping plants between light treatments. Leaf size is light responsive late into development, but patterning events, such as stomatal index, are irrevocably specified earlier. Observing that shade induces increases in shoot apical meristem size, we then describe gene expression changes in early leaf primordia and the meristem using laser microdissection. We find that in leaf primordia, shade avoidance is not mediated through canonical pathways described in mature organs but rather through the expression of KNOTTED1-LIKE HOMEOBOX and other indeterminacy genes, altering known developmental pathways responsible for patterning leaf shape. We also demonstrate that shade-induced changes in leaf primordium gene expression largely do not overlap with those found in successively initiated leaf primordia, providing evidence against classic hypotheses that shaded leaf morphology results from the prolonged production of juvenile leaf types.