Investigating the cis-Regulatory Basis of C3 and C4 Photosynthesis in Grasses at Single-Cell Resolution.

While considerable knowledge exists about the enzymes pivotal for C4 photosynthesis, much less is known about the cis-regulation important for specifying their expression in distinct cell types. Here, we use single-cell-indexed ATAC-seq to identify cell-type-specific accessible chromatin regions (ACRs) associated with C4 enzymes for five different grass species. This study spans four C4 species, covering three distinct photosynthetic subtypes: Zea mays and Sorghum bicolor (NADP-ME), Panicum miliaceum (NAD-ME), Urochloa fusca (PEPCK), along with the C3 outgroup Oryza sativa. We studied the cis-regulatory landscape of enzymes essential across all C4 species and those unique to C4 subtypes, measuring cell-type-specific biases for C4 enzymes using chromatin accessibility data. Integrating these data with phylogenetics revealed diverse co-option of gene family members between species, showcasing the various paths of C4 evolution. Besides promoter proximal ACRs, we found that, on average, C4 genes have two to three distal cell-type-specific ACRs, highlighting the complexity and divergent nature of C4 evolution. Examining the evolutionary history of these cell-type-specific ACRs revealed a spectrum of conserved and novel ACRs, even among closely related species, indicating ongoing evolution of cis-regulation at these C4 loci. This study illuminates the dynamic and complex nature of CRE evolution in C4 photosynthesis, particularly highlighting the intricate cis-regulatory evolution of key loci. Our findings offer a valuable resource for future investigations, potentially aiding in the optimization of C3 crop performance under changing climatic conditions.

Evolution of plant cell-type-specific cis -regulatory elements.

Cis -regulatory elements (CREs) are critical in regulating gene expression, and yet our understanding of CRE evolution remains a challenge. Here, we constructed a comprehensive single-cell atlas of chromatin accessibility in Oryza sativa , integrating data from 104,029 nuclei representing 128 discrete cell states across nine distinct organs. We used comparative genomics to compare cell-type resolved chromatin accessibility between O. sativa and 57,552 nuclei from four additional grass species ( Zea mays, Sorghum bicolor, Panicum miliaceum , and Urochloa fusca ). Accessible chromatin regions (ACRs) had different levels of conservation depending on the degree of cell-type specificity. We found a complex relationship between ACRs with conserved noncoding sequences, cell-type specificity, conservation, and tissue-specific switching. Additionally, we found that epidermal ACRs were less conserved compared to other cell types, potentially indicating that more rapid regulatory evolution has occurred in the L1 epidermal layer of these species. Finally, we identified and characterized a conserved subset of ACRs that overlapped the repressive histone modification H3K27me3, implicating them as potentially critical silencer CREs maintained by evolution. Collectively, this comparative genomics approach highlights the dynamics of cell-type-specific CRE evolution in plants.

Comparative genomic analyses reveal the genetic basis of the yellow-seed trait in Brassica napus.

Yellow-seed trait is a desirable breeding characteristic of rapeseed (Brassica napus) that could greatly improve seed oil yield and quality. However, the underlying mechanisms controlling this phenotype in B. napus plants are difficult to discern because of their complexity. Here, we assemble high-quality genomes of yellow-seeded (GH06) and black-seeded (ZY821). Combining in-depth fine mapping of a quantitative trait locus (QTL) for seed color with other omics data reveal BnA09MYB47a, encoding an R2R3-MYB-type transcription factor, as the causal gene of a major QTL controlling the yellow-seed trait. Functional studies show that sequence variation of BnA09MYB47a underlies the functional divergence between the yellow- and black-seeded B. napus. The black-seed allele BnA09MYB47aZY821, but not the yellow-seed allele BnA09MYB47aGH06, promotes flavonoid biosynthesis by directly activating the expression of BnTT18. Our discovery suggests a possible approach to breeding B. napus for improved commercial value and facilitates flavonoid biosynthesis studies in Brassica crops.

A multi-omic Nicotiana benthamiana resource for fundamental research and biotechnology.

Nicotiana benthamiana is an invaluable model plant and biotechnology platform with a ~3 Gb allotetraploid genome. To further improve its usefulness and versatility, we have produced high-quality chromosome-level genome assemblies, coupled with transcriptome, epigenome, microRNA and transposable element datasets, for the ubiquitously used LAB strain and a related wild accession, QLD. In addition, single nucleotide polymorphism maps have been produced for a further two laboratory strains and four wild accessions. Despite the loss of five chromosomes from the ancestral tetraploid, expansion of intergenic regions, widespread segmental allopolyploidy, advanced diploidization and evidence of recent bursts of Copia pseudovirus (Copia) mobility not seen in other Nicotiana genomes, the two subgenomes of N. benthamiana show large regions of synteny across the Solanaceae. LAB and QLD have many genetic, metabolic and phenotypic differences, including disparate RNA interference responses, but are highly interfertile and amenable to genome editing and both transient and stable transformation. The LAB/QLD combination has the potential to be as useful as the Columbia-0/Landsberg errecta partnership, utilized from the early pioneering days of Arabidopsis genomics to today.

Phytohormone ethylene mediates oligogalacturonic acid-induced growth inhibition in tomato etiolated seedlings.

Plant growth and immunity are tightly interconnected. Oligogalacturonic acids (OGs) are pectic fragments and have been well investigated in plant immunity as a damage-associated molecular pattern. However, little is known regarding how OGs affect plant growth. Here, we reveal that OGs inhibit the growth of intact etiolated seedling by using the horticultural crop tomato as a model. This inhibitory effect is partially suppressed by the action of ethylene biosynthesis inhibitors, or the gene silencing of SlACS2, an essential rate-limiting enzyme for ethylene biosynthesis, suggesting that SlACS2-mediated ethylene production promotes OG-induced growth inhibition. Furthermore, OGs treatment elevates the SlACS2 protein phosphorylation, and its decrease by the kinase inhibitor K252a partially rescue OG-induced growth inhibition, indicating that SlACS2 phosphorylation involves in OG-induced growth inhibition. Moreover, the mitogen-activated protein kinase SlMPK3 could be activated by OGs treatment and can directly phosphorylate SlACS2 in vitro, and the bimolecular fluorescence complementation combining with the yeast two-hybrid assay shows that SlMPK3 interacts with SlACS2, indicating that SlMPK3 may participate in modulating the OG-induced SlACS2 phosphorylation and growth inhibition. Our results reveal a regulatory mechanism at both the transcriptional and post-transcriptional levels by which OGs inhibit the growth of intact plant seedlings.

Tomato MED25 regulates fruit ripening by interacting with EIN3-like transcription factors.

Fruit ripening relies on the precise spatiotemporal control of RNA polymerase II (Pol II)-dependent gene transcription, and the evolutionarily conserved Mediator (MED) coactivator complex plays an essential role in this process. In tomato (Solanum lycopersicum), a model climacteric fruit, ripening is tightly coordinated by ethylene and several key transcription factors. However, the mechanism underlying the transmission of context-specific regulatory signals from these ripening-related transcription factors to the Pol II transcription machinery remains unknown. Here, we report the mechanistic function of MED25, a subunit of the plant Mediator transcriptional coactivator complex, in controlling the ethylene-mediated transcriptional program during fruit ripening. Multiple lines of evidence indicate that MED25 physically interacts with the master transcription factors of the ETHYLENE-INSENSITIVE 3 (EIN3)/EIN3-LIKE (EIL) family, thereby playing an essential role in pre-initiation complex formation during ethylene-induced gene transcription. We also show that MED25 forms a transcriptional module with EIL1 to regulate the expression of ripening-related regulatory as well as structural genes through promoter binding. Furthermore, the EIL1-MED25 module orchestrates both positive and negative feedback transcriptional circuits, along with its downstream regulators, to fine-tune ethylene homeostasis during fruit ripening.

Limited conservation in cross-species comparison of GLK transcription factor binding suggested wide-spread cistrome divergence.

Non-coding cis-regulatory variants in animal genomes are an important driving force in the evolution of transcription regulation and phenotype diversity. However, cistrome dynamics in plants remain largely underexplored. Here, we compare the binding of GOLDEN2-LIKE (GLK) transcription factors in tomato, tobacco, Arabidopsis, maize and rice. Although the function of GLKs is conserved, most of their binding sites are species-specific. Conserved binding sites are often found near photosynthetic genes dependent on GLK for expression, but sites near non-differentially expressed genes in the glk mutant are nevertheless under purifying selection. The binding sites' regulatory potential can be predicted by machine learning model using quantitative genome features and TF co-binding information. Our study show that genome cis-variation caused wide-spread TF binding divergence, and most of the TF binding sites are genetically redundant. This poses a major challenge for interpreting the effect of individual sites and highlights the importance of quantitatively measuring TF occupancy.

Characterization of regulatory modules controlling leaf angle in maize.

Leaf angle is an important agronomic trait determining maize (Zea mays) planting density and light penetration into the canopy and contributes to the yield gain in modern maize hybrids. However, little is known about the molecular mechanisms underlying leaf angle beyond the ZmLG1 (liguleless1) and ZmLG2 (Liguleless2) genes. In this study, we found that the transcription factor (TF) ZmBEH1 (BZR1/BES1 homolog gene 1) is targeted by ZmLG2 and regulates leaf angle formation by influencing sclerenchyma cell layers on the adaxial side. ZmBEH1 interacted with the TF ZmBZR1 (Brassinazole Resistant 1), whose gene expression was also directly activated by ZmLG2. Both ZmBEH1 and ZmBZR1 are bound to the promoter of ZmSCL28 (SCARECROW-LIKE 28), a third TF that influences leaf angle. Our study demonstrates regulatory modules controlling leaf angle and provides gene editing targets for creating optimal maize architecture suitable for dense planting.

A combinatorial indexing strategy for low-cost epigenomic profiling of plant single cells.

Understanding how cis-regulatory elements facilitate gene expression is a key question in biology. Recent advances in single-cell genomics have led to the discovery of cell-specific chromatin landscapes that underlie transcription programs in animal models. However, the high equipment and reagent costs of commercial systems limit their applications for many laboratories. In this study, we developed a combinatorial index and dual PCR barcode strategy to profile the Arabidopsis thaliana root single-cell epigenome without any specialized equipment. We generated chromatin accessibility profiles for 13 576 root nuclei with an average of 12 784 unique Tn5 integrations per cell. Integration of the single-cell assay for transposase-accessible chromatin sequencing and RNA sequencing data sets enabled the identification of 24 cell clusters with unique transcription, chromatin, and cis-regulatory signatures. Comparison with single-cell data generated using the commercial microfluidic platform from 10X Genomics revealed that this low-cost combinatorial index method is capable of unbiased identification of cell-type-specific chromatin accessibility. We anticipate that, by removing cost, instrumentation, and other technical obstacles, this method will be a valuable tool for routine investigation of single-cell epigenomes and provide new insights into plant growth and development and plant interactions with the environment.

Two zinc finger proteins with functions in m6A writing interact with HAKAI.

The methyltransferase complex (m6A writer), which catalyzes the deposition of N6-methyladenosine (m6A) in mRNAs, is highly conserved across most eukaryotic organisms, but its components and interactions between them are still far from fully understood. Here, using in vivo interaction proteomics, two HAKAI-interacting zinc finger proteins, HIZ1 and HIZ2, are discovered as components of the Arabidopsis m6A writer complex. HAKAI is required for the interaction between HIZ1 and MTA (mRNA adenosine methylase A). Whilst HIZ1 knockout plants have normal levels of m6A, plants in which it is overexpressed show reduced methylation and decreased lateral root formation. Mutant plants lacking HIZ2 are viable but have an 85% reduction in m6A abundance and show severe developmental defects. Our findings suggest that HIZ2 is likely the plant equivalent of ZC3H13 (Flacc) of the metazoan m6A-METTL Associated Complex.

Identification of Chilling Accumulation-Associated Genes for Litchi Flowering by Transcriptome-Based Genome-Wide Association Studies.

Litchi is an important Sapindaceae fruit tree. Flowering in litchi is triggered by low temperatures in autumn and winter. It can be divided into early-, medium-, and late-flowering phenotypes according to the time for floral induction. Early-flowering varieties need low chilling accumulation level for floral induction, whereas the late-flowering varieties require high chilling accumulation level. In the present study, RNA-Seq of 87 accessions was performed and transcriptome-based genome-wide association studies (GWAS) was used to identify candidate genes involved in chilling accumulation underlying the time for floral induction. A total of 98,155 high-quality single-nucleotide polymorphism (SNP) sites were obtained. A total of 1,411 significantly associated SNPs and 1,115 associated genes (AGs) were identified, of which 31 were flowering-related, 23 were hormone synthesis-related, and 27 were hormone signal transduction-related. Association analysis between the gene expression of the AGs and the flowering phenotypic data was carried out, and differentially expressed genes (DEGs) in a temperature-controlled experiment were obtained. As a result, 15 flowering-related candidate AGs (CAGs), 13 hormone synthesis-related CAGs, and 11 hormone signal transduction-related CAGs were further screened. The expression levels of the CAGs in the early-flowering accessions were different from those in the late-flowering ones, and also between the flowering trees and non-flowering trees. In a gradient chilling treatment, flowering rates of the trees and the CAGs expression were affected by the treatment. Our present work for the first time provided candidate genes for genetic regulation of flowering in litchi using transcriptome-based GWAS.

SnRK2 subfamily I protein kinases regulate ethylene biosynthesis by phosphorylating HB transcription factors to induce ACO1 expression in apple.

Ethylene (ETH) controls climacteric fruit ripening and can be triggered by osmotic stress. However, the mechanism regulating ETH biosynthesis during fruit ripening and under osmotic stress is largely unknown in apple (Malus domestica). Here, we explored the roles of SnRK2 protein kinases in ETH biosynthesis related to fruit ripening and osmoregulation. We identified the substrates of MdSnRK2-I using phosphorylation analysis techniques. Finally, we identified the MdSnRK2-I-mediated signaling pathway for ETH biosynthesis related to fruit ripening and osmoregulation. The activity of two MdSnRK2-I members, MdSnRK2.4 and MdSnRK2.9, was significantly upregulated during ripening or following mannitol treatment. Overexpression of MdSnRK2-I increased ETH biosynthesis under normal and osmotic conditions in apple fruit. MdSnRK2-I phosphorylated the transcription factors MdHB1 and MdHB2 to enhance their protein stability and transcriptional activity on MdACO1. MdSnRK2-I also interacted with MdACS1 and increased its protein stability through two phosphorylation sites. The increased MdACO1 expression and MdACS1 protein stability resulted in higher ETH production in apple fruit. In addition, heterologous expression of MdSnRK2-I or manipulation of SlSnRK2-I expression in tomato (Solanum lycopersicum) fruit altered fruit ripening and ETH biosynthesis. We established that MdSnRK2-I functions in fruit ripening and osmoregulation, and identified the MdSnRK2-I-mediated signaling pathway controlling ETH biosynthesis.

Chromatin and regulatory differentiation between bundle sheath and mesophyll cells in maize.

C4 plants partition photosynthesis enzymes between the bundle sheath (BS) and the mesophyll (M) cells for the better delivery of CO2 to RuBisCO and to reduce photorespiration. To better understand how C4 photosynthesis is regulated at the transcriptional level, we performed RNA-seq, ATAC-seq, ChIP-seq and Bisulfite-seq (BS-seq) on BS and M cells isolated from maize leaves. By integrating differentially expressed genes with chromatin features, we found that chromatin accessibility coordinates with epigenetic features, especially H3K27me3 modification and CHH methylation, to regulate cell type-preferentially enriched gene expression. Not only the chromatin-accessible regions (ACRs) proximal to the genes (pACRs) but also the distal ACRs (dACRs) are determinants of cell type-preferentially enriched expression. We further identified cell type-preferentially enriched motifs, e.g. AAAG for BS cells and TGACC/T for M cells, and determined their corresponding transcription factors: DOFs and WRKYs. The complex interaction between cis and trans factors in the preferential expression of C4 genes was also observed. Interestingly, cell type-preferentially enriched gene expression can be fine-tuned by the coordination of multiple chromatin features. Such coordination may be critical in ensuring the cell type-specific function of key C4 genes. Based on the observed cell type-preferentially enriched expression pattern and coordinated chromatin features, we predicted a set of functionally unknown genes, e.g. Zm00001d042050 and Zm00001d040659, to be potential key C4 genes. Our findings provide deep insight into the architectures associated with C4 gene expression and could serve as a valuable resource to further identify the regulatory mechanisms present in C4 species.

McMYB4 improves temperature adaptation by regulating phenylpropanoid metabolism and hormone signaling in apple.

Temperature changes affect apple development and production. Phenylpropanoid metabolism and hormone signaling play a crucial role in regulating apple growth and development in response to temperature changes. Here, we found that McMYB4 is induced by treatment at 28 °C and 18 °C, and McMYB4 overexpression results in flavonol and lignin accumulation in apple leaves. Yeast one-hybrid (Y1H) assays and electrophoretic mobility shift assays (EMSAs) further revealed that McMYB4 targets the promoters of the flavonol biosynthesis genes CHS and FLS and the lignin biosynthesis genes CAD and F5H. McMYB4 expression resulted in higher levels of flavonol and lignin biosynthesis in apple during growth at 28 °C and 18 °C than during growth at 23 °C. At 28 °C and 18 °C, McMYB4 also binds to the AUX/ARF and BRI/BIN promoters to activate gene expression, resulting in acceleration of the auxin and brassinolide signaling pathways. Taken together, our results demonstrate that McMYB4 promotes flavonol biosynthesis and brassinolide signaling, which decreases ROS contents to improve plant resistance and promotes lignin biosynthesis and auxin signaling to regulate plant growth. This study suggests that McMYB4 participates in the abiotic resistance and growth of apple in response to temperature changes by regulating phenylpropanoid metabolism and hormone signaling.

Missing Nurse Bees-Early Transcriptomic Switch From Nurse Bee to Forager Induced by Sublethal Imidacloprid.

The environmental residue/sublethal doses of neonicotinoid insecticides are believed to generate a negative impact on pollinators, including honey bees. Here we report our recent investigation on how imidacloprid, one of the major neonicotinoids, affects worker bees by profiling the transcriptomes of various ages of bees exposed to different doses of imidacloprid during the larval stage. The results show that imidacloprid treatments during the larval stage severely altered the gene expression profiles and may induce precocious foraging. Differential expression of foraging regulators was found in 14-day-old treated adults. A high transcriptome similarity between larvae-treated 14-day-old adults and 20-day-old controls was also observed, and the similarity was positively correlated with the dose of imidacloprid. One parts per billion (ppb) of imidacloprid was sufficient to generate a long-term impact on the bee's gene expression as severe as with 50 ppb imidacloprid. The disappearance of nurse bees may be driven not only by the hive member constitution but also by the neonicotinoid-induced precocious foraging behavior.

Effects of artificial light at night (ALAN) on gene expression of Aquatica ficta firefly larvae.

Artificial light at night (ALAN) is a major driver of firefly population declines, but its physiological effects are not well understood. To investigate the impact of ALAN on firefly development, we exposed larval Aquatica ficta fireflies to ALAN for two weeks. High larval mortality was observed in the periods of 1-68 days and 106-134 days post-treatment, which may represent the short- and long-term impacts of ALAN. We then profiled the transcriptome of larval Aquatica ficta fireflies following two weeks of ALAN exposure. A total of 1262 (1.67% out of 75777 unigenes) were differentially expressed in the treatment group: 1157 were down-regulated, and 105 were up-regulated. Up-regulated unigenes were related to regulation of hormone levels, ecdysteroid metabolic process, and response to stimulus; down-regulated unigenes were related to negative regulation of insulin receptor signaling, germ cell development, oogenesis, spermatid development, and regulation of neuron differentiation. Transcriptome results suggest that the endocrine, reproductive, and neural development of firefly larvae could be impaired by even relatively brief period of ALAN exposure. This report contributes a much-needed molecular perspective to the growing body of research documenting the fitness impacts of ALAN on bioluminescent fireflies.

Non-homology-based prediction of gene functions in maize (Zea mays ssp. mays).

Advances in genome sequencing and annotation have eased the difficulty of identifying new gene sequences. Predicting the functions of these newly identified genes remains challenging. Genes descended from a common ancestral sequence are likely to have common functions. As a result, homology is widely used for gene function prediction. This means functional annotation errors also propagate from one species to another. Several approaches based on machine learning classification algorithms were evaluated for their ability to accurately predict gene function from non-homology gene features. Among the eight supervised classification algorithms evaluated, random-forest-based prediction consistently provided the most accurate gene function prediction. Non-homology-based functional annotation provides complementary strengths to homology-based annotation, with higher average performance in Biological Process GO terms, the domain where homology-based functional annotation performs the worst, and weaker performance in Molecular Function GO terms, the domain where the accuracy of homology-based functional annotation is highest. GO prediction models trained with homology-based annotations were able to successfully predict annotations from a manually curated "gold standard" GO annotation set. Non-homology-based functional annotation based on machine learning may ultimately prove useful both as a method to assign predicted functions to orphan genes which lack functionally characterized homologs, and to identify and correct functional annotation errors which were propagated through homology-based functional annotations.

Reconstructing the maize leaf regulatory network using ChIP-seq data of 104 transcription factors.

The transcription regulatory network inside a eukaryotic cell is defined by the combinatorial actions of transcription factors (TFs). However, TF binding studies in plants are too few in number to produce a general picture of this complex network. In this study, we use large-scale ChIP-seq to reconstruct it in the maize leaf, and train machine-learning models to predict TF binding and co-localization. The resulting network covers 77% of the expressed genes, and shows a scale-free topology and functional modularity like a real-world network. TF binding sequence preferences are conserved within family, while co-binding could be key for their binding specificity. Cross-species comparison shows that core network nodes at the top of the transmission of information being more conserved than those at the bottom. This study reveals the complex and redundant nature of the plant transcription regulatory network, and sheds light on its architecture, organizing principle and evolutionary trajectory.