RESUMEN
Stomata play critical roles in gas exchange and immunity to pathogens. While many genes regulating early stomatal development up to the production of young guard cells (GCs) have been identified in Arabidopsis, much less is known about how young GCs develop into mature functional stomata. Here we perform a maturomics study on stomata, with "maturomics" defined as omics analysis of the maturation process of a tissue or organ. We develop an integrative scheme to analyze three public stomata-related single-cell RNA-seq datasets and identify a list of 586 genes that are specifically up-regulated in all three datasets during stomatal maturation and function formation. The list, termed sc_586, is enriched with known regulators of stomatal maturation and functions. To validate the reliability of the dataset, we selected two candidate G2-like transcription factor genes, MYS1 and MYS2, to investigate their roles in stomata. These two genes redundantly regulate the size and hoop rigidity of mature GCs, and the mys1 mys2 double mutants cause mature GCs with severe defects in regulating their stomatal apertures. Taken together, our results provide a valuable list of genes for studying GC maturation and function formation.
RESUMEN
Single-cell RNA-sequencing datasets of Arabidopsis roots have been generated, but related comprehensive gene co-expression network analyses are lacking. We conducted a single-cell gene co-expression network analysis with publicly available scRNA-seq datasets of Arabidopsis roots using a SingleCellGGM algorithm. The analysis identified 149 gene co-expression modules, which we considered to be gene expression programs (GEPs). By examining their spatiotemporal expression, we identified GEPs specifically expressed in major root cell types along their developmental trajectories. These GEPs define gene programs regulating root cell development at different stages and are enriched with relevant developmental regulators. As examples, a GEP specific for the quiescent center (QC) contains 20 genes regulating QC and stem cell niche homeostasis, and four GEPs are expressed in sieve elements (SEs) from early to late developmental stages, with the early-stage GEP containing 17 known SE developmental regulators. We also identified GEPs for metabolic pathways with cell-type-specific expression, suggesting the existence of cell-type-specific metabolism in roots. Using the GEPs, we discovered and verified a columella-specific gene, NRL27, as a regulator of the auxin-related root gravitropism response. Our analysis thus systematically reveals GEPs that regulate Arabidopsis root development and metabolism and provides ample resources for root biology studies.
RESUMEN
Identifying trait-associated genes is critical for rice (Oryza sativa) improvement, which usually relies on map-based cloning, quantitative trait locus analysis, or genome-wide association studies. Here we show that trait-associated genes tend to form modules within rice gene co-expression networks, a feature that can be exploited to discover additional trait-associated genes using reverse genetics. We constructed a rice gene co-expression network based on the graphical Gaussian model using 8,456 RNA-seq transcriptomes, which assembled into 1,286 gene co-expression modules functioning in diverse pathways. A number of the modules were enriched with genes associated with agronomic traits, such as grain size, grain number, tiller number, grain quality, leaf angle, stem strength, and anthocyanin content, and these modules are considered to be trait-associated gene modules. These trait-associated gene modules can be used to dissect the genetic basis of rice agronomic traits and to facilitate the identification of trait genes. As an example, we identified a candidate gene, OCTOPUS-LIKE 1 (OsOPL1), a homolog of the Arabidopsis (Arabidopsis thaliana) OCTOPUS gene, from a grain size module and verified it as a regulator of grain size via functional studies. Thus, our network represents a valuable resource for studying trait-associated genes in rice.