Conventional molecular analyses provide bulk genomic/transcriptomic data that are unable to reveal the cellular heterogeneity and to precisely define how gene networks orchestrate organ development. We will profile gene expression and identify open chromatin regions at the individual cells level, allowing to define cell-type specific regulatory elements, developmental trajectories and transcriptional networks orchestrating organ development and function. We will perform scRNA-seq and snATAC-seq on leaf/berry protoplasts and nuclei and combine them with the leaf/berry bulk tissues obtained results, where the analysis of transcripts, chromatin accessibility, histone modification and transcription factor binding sites showed that a large fraction of phenotypic variation appears to be determined by regulatory rather than coding variation and that many variants have an organ-specific effect.