Regulatory Divergence in Wound-Responsive Gene Expression between Domesticated and Wild Tomato

Nanofabrication of silk microneedles for high-throughput micronutrient delivery and continuous sap monitoring in plants

Biomaterials bridging the biotic-abiotic interface in plants offer the opportunity to precisely deliver agrochemicals and continuously monitor plant health, with the goals of increasing resilience to climate change, enhancing crop production and mitigating environmental impact. In this study we report the manipulation of silk fibroin assembly with inorganics nucleation at their phase front to nanomanufacture porous and hollow microneedles that can be interfaced with plants. Plant growth analysis and quantification of wounding gene expression show a non-significant systemic wounding response to the injection of silk microneedles in tomato plants. Microneedles with a hollow structure enable the systemic delivery of plant micronutrients to treat chlorosis in tomato plants and crop biofortification through transport of human micronutrients injected in the petiole and loaded into tomato fruits. Hollow microneedles also provide access to plant vasculature for sap sampling, enabling continuous monitoring and early detection of phytoaccumulation of environmental contaminants such as cadmium.

Genetic and functional traits limit the success of colonisation by arbuscular mycorrhizal fungi in a tomato wild relative

To understand whether domestication had an impact on susceptibility and responsiveness to arbuscular mycorrhizal fungi (AMF) in tomato (Solanum lycopersicum), we investigated two tomato cultivars ("M82" and "Moneymaker") and a panel of wild relatives including S. neorickii, S. habrochaites and S. pennellii encompassing the whole Lycopersicon clade. Most genotypes revealed good AM colonisation levels when inoculated with the AMF Funneliformis mosseae. By contrast, both S. pennellii accessions analysed showed a very low colonisation, but with normal arbuscule morphology, and a negative response in terms of root and shoot biomass. This behaviour was independent of fungal identity and environmental conditions. Genomic and transcriptomic analyses revealed in S. pennellii the lack of genes identified within QTLs for AM colonisation, a limited transcriptional reprogramming upon mycorrhization and a differential regulation of strigolactones and AM-related genes compared to tomato. Donor plants experiments indicated that the AMF could represent a cost for S. pennellii: F. mosseae could extensively colonise the root only when it was part of a mycorrhizal network, but a higher mycorrhization led to a higher inhibition of plant growth. These results suggest that genetics and functional traits of S. pennellii are responsible for the limited extent of AMF colonisation.

Modeling alternative translation initiation sites in plants reveals evolutionarily conserved cis-regulatory codes in eukaryotes

mRNA translation relies on identifying translation initiation sites (TISs) in mRNAs. Alternative TISs are prevalent across plant transcriptomes, but the mechanisms for their recognition are unclear. Using ribosome profiling and machine learning, we developed models for predicting alternative TISs in the tomato (Solanum lycopersicum). Distinct feature sets were predictive of AUG and nonAUG TISs in 5' untranslated regions and coding sequences, including a novel CU-rich sequence that promoted plant TIS activity, a translational enhancer found across dicots and monocots, and humans and viruses. Our results elucidate the mechanistic and evolutionary basis of TIS recognition, whereby cis-regulatory RNA signatures affect start site selection. The TIS prediction model provides global estimates of TISs to discover neglected protein-coding genes across plant genomes. The prevalence of cis-regulatory signatures across plant species, humans, and viruses suggests their broad and critical roles in reprogramming the translational landscape.

Predicting transcriptional responses to heat and drought stress from genomic features using a machine learning approach in rice

Plants have evolved various mechanisms to adapt to adverse environmental stresses, such as the modulation of gene expression. Expression of stress-responsive genes is controlled by specific regulators, including transcription factors (TFs), that bind to sequence-specific binding sites, representing key components of cis-regulatory elements and regulatory networks. Our understanding of the underlying regulatory code remains, however, incomplete. Recent studies have shown that, by training machine learning (ML) algorithms on genomic sequence features, it is possible to predict which genes will transcriptionally respond to a specific stress. By identifying the most important features for gene expression prediction, these trained ML models allow, in theory, to further elucidate the regulatory code underlying the transcriptional response to abiotic stress. Here, we trained random forest ML models to predict gene expression in rice (Oryza sativa) in response to heat or drought stress. Apart from thoroughly assessing model performance and robustness across various input training data, the importance of promoter and gene body sequence features to train ML models was evaluated. The use of enriched promoter oligomers, complementing known TF binding sites, allowed us to gain novel insights in DNA motifs contributing to the stress regulatory code. By comparing genomic feature importance scores for drought and heat stress over time, general and stress-specific genomic features contributing to the performance of the learned models and their temporal variation were identified. This study provides a solid foundation to build and interpret ML models accurately predicting transcriptional responses and enables novel insights in biological sequence features that are important for abiotic stress responses.

Plant domestication: setting biological clocks

Through domestication of wild species, humans have induced large changes in the developmental and circadian clocks of plants. As a result of these changes, modern crops are more productive and adaptive to contrasting environments from the center of origin of their wild ancestors, albeit with low genetic variability and abiotic stress tolerance. Likewise, a complete restructuring of plant metabolic timekeeping probably occurred during crop domestication. Here, we highlight that contrasting timings among organs in wild relatives of crops allowed them to recognize environmental adversities faster. We further propose that connections among biological clocks, which were established during plant domestication, may represent a fundamental source of genetic variation to improve crop resilience and yield.

Temporal regulation of cold transcriptional response in switchgrass

Switchgrass low-land ecotypes have significantly higher biomass but lower cold tolerance compared to up-land ecotypes. Understanding the molecular mechanisms underlying cold response, including the ones at transcriptional level, can contribute to improving tolerance of high-yield switchgrass under chilling and freezing environmental conditions. Here, by analyzing an existing switchgrass transcriptome dataset, the temporal cis-regulatory basis of switchgrass transcriptional response to cold is dissected computationally. We found that the number of cold-responsive genes and enriched Gene Ontology terms increased as duration of cold treatment increased from 30 min to 24 hours, suggesting an amplified response/cascading effect in cold-responsive gene expression. To identify genomic sequences likely important for regulating cold response, machine learning models predictive of cold response were established using k-mer sequences enriched in the genic and flanking regions of cold-responsive genes but not non-responsive genes. These k-mers, referred to as putative cis-regulatory elements (pCREs) are likely regulatory sequences of cold response in switchgrass. There are in total 655 pCREs where 54 are important in all cold treatment time points. Consistent with this, eight of 35 known cold-responsive CREs were similar to top-ranked pCREs in the models and only these eight were important for predicting temporal cold response. More importantly, most of the top-ranked pCREs were novel sequences in cold regulation. Our findings suggest additional sequence elements important for cold-responsive regulation previously not known that warrant further studies.

TDTHub, a web server tool for the analysis of transcription factor binding sites in plants

Transcriptional regulation underlies most developmental programs and physiological responses to environmental changes in plants. Transcription factors (TFs) play a key role in the regulation of gene expression by binding specifically to short DNA sequences in the regulatory regions of genes: the TF binding sites (TFBSs). In recent years, several bioinformatic tools have been developed to detect TFBSs in candidate genes, either by de novo prediction or by directly mapping experimentally known TFBSs. However, most of these tools contain information for only a few species or require multi-step procedures, and are not always intuitive for non-experienced researchers. Here we present TFBS-Discovery Tool Hub (TDTHub), a web server for quick and intuitive studies of transcriptional regulation in plants. TDTHub uses pre-computed TFBSs in 40 plant species and allows the choice of two mapping algorithms, providing a higher versatility. Besides the main TFBS enrichment tool, TDTHub includes additional tools to assist in the analysis and visualization of data. In order to demonstrate the effectiveness of TDTHub, we analyzed the transcriptional regulation of the anthocyanin biosynthesis pathway. We also analyzed the transcriptional cascades in response to jasmonate and wounding in Arabidopsis and tomato (Solanum lycopersicum), respectively. In these studies, TDTHub helped to verify the most relevant TF nodes and to propose new ones with a prominent role in these pathways. TDTHub is available at http://acrab.cnb.csic.es/TDTHub/, and it will be periodically upgraded and expanded for new species and gene annotations.

Modeling temporal and hormonal regulation of plant transcriptional response to wounding

Plants respond to wounding stress by changing gene expression patterns and inducing the production of hormones including jasmonic acid. This wounding transcriptional response activates specialized metabolism pathways such as the glucosinolate pathways in Arabidopsis thaliana. While the regulatory factors and sequences controlling a subset of wound-response genes are known, it remains unclear how wound response is regulated globally. Here, we how these responses are regulated by incorporating putative cis-regulatory elements, known transcription factor binding sites, in vitro DNA affinity purification sequencing, and DNase I hypersensitive sites to predict genes with different wound-response patterns using machine learning. We observed that regulatory sites and regions of open chromatin differed between genes upregulated at early and late wounding time-points as well as between genes induced by jasmonic acid and those not induced. Expanding on what we currently know, we identified cis-elements that improved model predictions of expression clusters over known binding sites. Using a combination of genome editing, in vitro DNA-binding assays, and transient expression assays using native and mutated cis-regulatory elements, we experimentally validated four of the predicted elements, three of which were not previously known to function in wound-response regulation. Our study provides a global model predictive of wound response and identifies new regulatory sequences important for wounding without requiring prior knowledge of the transcriptional regulators.

Intergenic spaces: a new frontier to improving plant health

To more sustainably mitigate the impact of crop diseases on plant health and productivity, there is a need for broader spectrum, long-lasting resistance traits. Defense response (DR) genes, located throughout the genome, participate in cellular and system-wide defense mechanisms to stave off infection by diverse pathogens. This multigenic resistance avoids rapid evolution of a pathogen to overcome host resistance. DR genes reside within resistance-associated quantitative trait loci (QTL), and alleles of DR genes in resistant varieties are more active during pathogen attack relative to susceptible haplotypes. Differential expression of DR genes results from polymorphisms in their regulatory regions, that includes cis-regulatory elements such as transcription factor binding sites as well as features that influence epigenetic structural changes to modulate chromatin accessibility during infection. Many of these elements are found in clusters, known as cis-regulatory modules (CRMs), which are distributed throughout the host genome. Regulatory regions involved in plant-pathogen interactions may also contain pathogen effector binding elements that regulate DR gene expression, and that, when mutated, result in a change in the plants' response. We posit that CRMs and the multiple regulatory elements that comprise them are potential targets for marker-assisted breeding for broad-spectrum, durable disease resistance.

A tomato receptor-like cytoplasmic kinase, SlZRK1, acts as a negative regulator in wound-induced jasmonic acid accumulation and insect resistance

Jasmonates accumulate rapidly and act as key regulators in response to mechanical wounding, but few studies have linked receptor-like cytoplasmic kinases (RLCKs) to wound-induced jasmonic acid (JA) signaling cascades. Here, we identified a novel wounding-induced RLCK-XII-2 subfamily member (SlZRK1) in tomato (Solanum lycopersicum) that was closely related to Arabidopsis HOPZ-ETI-DEFICIENT 1 (ZED1)-related kinases 1 based on phylogenetic analysis. SlZRK1 was targeted to the plasma membrane of tobacco mesophyll protoplasts as determined by transient co-expression with the plasma membrane marker mCherry-H+-ATPase. Catalytic residue sequence analysis and an in vitro kinase assay indicated that SlZRK1 may act as a pseudokinase. To further analyse the function of SlZRK1, we developed two stable knock-out mutants by CRISPR/Cas9. Loss of SlZRK1 significantly altered the expression of genes involved in JA biosynthesis, salicylic acid biosynthesis, and ethylene response. Furthermore, after mechanical wounding treatment, slzrk1 mutants increased transcription of early wound-inducible genes involved in JA biosynthesis and signaling. In addition, JA accumulation after wounding and plant resistance to herbivorous insects also were enhanced. Our findings expand plant regulatory networks in the wound-induced JA production by adding RLCKs as a new component in the wound signal transduction pathway.

Evolutionarily conserved hierarchical gene regulatory networks for plant salt stress response

Plant cells constantly alter their gene expression profiles to respond to environmental fluctuations. These continuous adjustments are regulated by multi-hierarchical networks of transcription factors. To understand how such gene regulatory networks (GRNs) have stabilized evolutionarily while allowing for species-specific responses, we compare the GRNs underlying salt response in the early-diverging and late-diverging plants Marchantia polymorpha and Arabidopsis thaliana. Salt-responsive GRNs, constructed on the basis of the temporal transcriptional patterns in the two species, share common trans-regulators but exhibit an evolutionary divergence in cis-regulatory sequences and in the overall network sizes. In both species, WRKY-family transcription factors and their feedback loops serve as central nodes in salt-responsive GRNs. The divergent cis-regulatory sequences of WRKY-target genes are probably associated with the expansion in network size, linking salt stress to tissue-specific developmental and physiological responses. The WRKY modules and highly linked WRKY feedback loops have been preserved widely in other plants, including rice, while keeping their binding-motif sequences mutable. Together, the conserved trans-regulators and the quickly evolving cis-regulatory sequences allow salt-responsive GRNs to adapt over a long evolutionary timescale while maintaining some consistent regulatory structure. This strategy may benefit plants as they adapt to changing environments.

Network-based approaches for understanding gene regulation and function in plants

Expression reprogramming directed by transcription factors is a primary gene regulation underlying most aspects of the biology of any organism. Our views of how gene regulation is coordinated are dramatically changing thanks to the advent and constant improvement of high-throughput profiling and transcriptional network inference methods: from activities of individual genes to functional interactions across genes. These technical and analytical advances can reveal the topology of transcriptional networks in which hundreds of genes are hierarchically regulated by multiple transcription factors at systems level. Here we review the state of the art of experimental and computational methods used in plant biology research to obtain large-scale datasets and model transcriptional networks. Examples of direct use of these network models and perspectives on their limitations and future directions are also discussed.

The cis-regulatory codes of response to combined heat and drought stress in Arabidopsis thaliana

Plants respond to their environment by dynamically modulating gene expression. A powerful approach for understanding how these responses are regulated is to integrate information about cis-regulatory elements (CREs) into models called cis-regulatory codes. Transcriptional response to combined stress is typically not the sum of the responses to the individual stresses. However, cis-regulatory codes underlying combined stress response have not been established. Here we modeled transcriptional response to single and combined heat and drought stress in Arabidopsis thaliana. We grouped genes by their pattern of response (independent, antagonistic and synergistic) and trained machine learning models to predict their response using putative CREs (pCREs) as features (median F-measure = 0.64). We then developed a deep learning approach to integrate additional omics information (sequence conservation, chromatin accessibility and histone modification) into our models, improving performance by 6.2%. While pCREs important for predicting independent and antagonistic responses tended to resemble binding motifs of transcription factors associated with heat and/or drought stress, important synergistic pCREs resembled binding motifs of transcription factors not known to be associated with stress. These findings demonstrate how in silico approaches can improve our understanding of the complex codes regulating response to combined stress and help us identify prime targets for future characterization.

Identification and evolution of gene regulatory networks: insights from comparative studies in plants

The availability of genome sequences, genome-wide assays of transcription factor binding, and accessible chromatin maps have unveiled gene regulatory landscapes in plants. This understanding has ushered in comparative gene regulatory network studies that assess network rewiring between species, across time, and between biological tissues. Comparisons of cis-regulatory elements across the plant kingdom have uncovered examples of conserved sequences, but also of divergence, indicating that selective pressures can vary in different plant families. Transcription factor duplication, followed by spatiotemporal expression divergence of the duplicates, also appears to be a key mechanism of network evolution. Here, we review recent literature describing the regulation of gene expression in plants, and how comparative studies provide insights into how these regulatory interactions change and lead to gene regulatory network rewiring.

Accumulation of Ascorbic Acid in Tomato Cell Culture: Influence of the Genotype, Source Explant and Time of In Vitro Cultivation

The production and commercialization of natural antioxidants is gaining increasing importance due to their wide range of biological effects and applications. In vitro cell culture is a valuable source of plant bioactive compounds, especially those highly dependent on environmental factors. Nonetheless, research on the accumulation in plant cultured cells of water-soluble antioxidant vitamins, such as the ascorbic acid (AsA), is very limited. Tomato fruits are a main dietary source of vitamin C and in this work, we explored the potential of in vitro cultured cells for AsA accumulation. Specifically, using a full factorial design, we examined the effect of the source explant, the time in tissue culture and the genetic difference present in two Introgression Line (IL7-3 and IL12-4) that harbor Quantitative Trait Loci (QTLs) for ascorbic acid in fruits. Moreover, we performed an expression analysis of genes involved in AsA metabolism to highlight the molecular mechanisms that can account for the difference between fruit explants and calli. Our work indicated that cultured tomato cells accumulate AsA well beyond the amount present in fruits and that the three factors under investigation and their interaction significantly influence AsA accumulation. The time in tissue culture is the main single factor and, different from the expectations for secondary metabolites, explants from unripe, mature green fruits provided the highest increase in AsA. Moreover, in controlled conditions the genetic differences between the ILs and the control genotype are less relevant for calli cultivated for longer time. Our work showed the potential of tomato cell culture to produce AsA and prompt further refinements towards its possible large-scale exploitation.

A rapid preparation procedure for laser microdissection-mediated harvest of plant tissues for gene expression analysis

BACKGROUND

Gene expression changes that govern essential biological processes can occur at the cell-specific level. To gain insight into such events, laser microdissection is applied to cut out specific cells or tissues from which RNA for gene expression analysis is isolated. However, the preparation of plant tissue sections for laser microdissection and subsequent RNA isolation usually involves fixation and embedding, processes that are often time-consuming and can lower the yield and quality of isolated RNA.

RESULTS

Infection sites of the parasitic plant Cuscuta reflexa growing on its compatible host plant Pelargonium zonale were sectioned using a vibratome and dried on glass slides at 4 °C before laser microdissection. High quality RNA (RQI > 7) was isolated from 1 mm2, 3 mm2 and 6 mm2 total surface areas of laser microdissection-harvested C. reflexa tissue, with the yield of RNA correlating to the amount of collected material (on average 7 ng total RNA/mm2). The expression levels of two parasite genes previously found to be highly expressed during host plant infection were shown to differ individually between specific regions of the infection site. By drying plant sections under low pressure to reduce the dehydration time, the induced expression of two wound-related genes during preparation was avoided.

CONCLUSIONS

Plants can be prepared quickly and easily for laser microdissection by direct sectioning of fresh tissue followed by dehydration on glass slides. We show that RNA isolated from material treated in this manner maintains high quality and enables the investigation of differential gene expression at a high morphological resolution.

Expression Divergence as an Evolutionary Alternative Mechanism Adopted by Two Rice Subspecies Against Rice Blast Infection

BACKGROUND

Rice (Oryza sativa L.) is one of the most important crops that serves as staple food for ~ 50% of the human population worldwide. Some important agronomic traits that allow rice to cope with numerous abiotic and biotic stresses have been selected and fixed during domestication. Knowledge on how expression divergence of genes gradually contributes to phenotypic differentiation in response to biotic stress and their contribution to rice population speciation is still limited.

RESULTS

Here, we explored gene expression divergence between a japonica rice cultivar Nipponbare and an indica rice cultivar 93-11 in response to invasion by the filamentous ascomycete fungus Magnaporthe oryzae (Pyricularia oryzae), a plant pathogen that causes significant loss to rice production worldwide. We investigated differentially expressed genes in the two cultivars and observed that evolutionarily conserved orthologous genes showed highly variable expression patterns under rice blast infection. Analysis of promoter region of these differentially expressed orthologous genes revealed the existence of cis-regulatory elements associated with the differentiated expression pattern of these genes in the two rice cultivars. Further comparison of these regions in global rice population indicated their fixation and close relationship with rice population divergence.

CONCLUSION

We proposed that variation in the expression patterns of these orthologous genes mediated by cis-regulatory elements in the two rice cultivars, may constitute an alternative evolutionary mechanism that distinguishes these two genetically and ecologically divergent rice cultivars in response to M. oryzae infection.