Regulators of early maize leaf development inferred from transcriptomes of laser capture microdissection (LCM)-isolated embryonic leaf cells.

The superior photosynthetic efficiency of C4 leaves over C3 leaves is owing to their unique Kranz anatomy, in which the vein is surrounded by one layer of bundle sheath (BS) cells and one layer of mesophyll (M) cells. Kranz anatomy development starts from three contiguous ground meristem (GM) cells, but its regulators and underlying molecular mechanism are largely unknown. To identify the regulators, we obtained the transcriptomes of 11 maize embryonic leaf cell types from five stages of pre-Kranz cells starting from median GM cells and six stages of pre-M cells starting from undifferentiated cells. Principal component and clustering analyses of transcriptomic data revealed rapid pre-Kranz cell differentiation in the first two stages but slow differentiation in the last three stages, suggesting early Kranz cell fate determination. In contrast, pre-M cells exhibit a more prolonged transcriptional differentiation process. Differential gene expression and coexpression analyses identified gene coexpression modules, one of which included 3 auxin transporter and 18 transcription factor (TF) genes, including known regulators of Kranz anatomy and/or vascular development. In situ hybridization of 11 TF genes validated their expression in early Kranz development. We determined the binding motifs of 15 TFs, predicted TF target gene relationships among the 18 TF and 3 auxin transporter genes, and validated 67 predictions by electrophoresis mobility shift assay. From these data, we constructed a gene regulatory network for Kranz development. Our study sheds light on the regulation of early maize leaf development and provides candidate leaf development regulators for future study.

Simultaneous detection of miRNA and mRNA at the single-cell level in plant tissues.

Detecting the simultaneous presence of a microRNA (miRNA) and a mRNA in a specific tissue can provide support for the prediction that the miRNA regulates the mRNA. Although two such methods have been developed for mammalian tissues, they have a low signal-noise ratio and/or poor resolution at the single-cell level. To overcome these drawbacks, we develop a method that uses sequence-specific miRNA-locked nucleic acid (LNA) and mRNA-LNA probes. Moreover, it augments the detection signal by rolling circle amplification, achieving a high signal-noise ratio at the single-cell level. Dot signals are counted for determining the expression levels of mRNA and miRNA molecules in specific cells. We show a high sequence specificity of our miRNA-LNA probe, revealing that it can discriminate single-base mismatches. Numerical quantification by our method is tested in transgenic rice lines with different gene expression levels. We conduct several applications. First, the spatial expression profiling of osa-miR156 and OsSPL12 in rice leaves reveals their specific expression in mesophyll cells. Second, studying rice and its mutant lines with our method reveals opposite expression patterns of miRNA and its target mRNA in tissues. Third, the dynamic expression profiles of ZmGRF8 and zma-miR396 during maize leaf development provide evidence that zma-miR396 regulates the preferential spatial expression of ZmGRF8 in bundle sheath cells. Finally, our method can be scaled up to simultaneously detect multiple miRNAs and mRNAs in a tissue. Thus, it is a sensitive and versatile technique for studying miRNA regulation of plant tissue development.

Maize Golden2-like transcription factors boost rice chloroplast development, photosynthesis, and grain yield.

Chloroplasts are the sites for photosynthesis, and two Golden2-like factors act as transcriptional activators of chloroplast development in rice (Oryza sativa L.) and maize (Zea mays L.). Rice OsGLK1 and OsGLK2 are orthologous to maize ZmGLK1 (ZmG1) and ZmGLK2 (ZmG2), respectively. However, while rice OsGLK1 and OsGLK2 act redundantly to regulate chloroplast development in mesophyll cells, maize ZmG1 and ZmG2 are functionally specialized and expressed in different cell-specific manners. To boost rice chloroplast development and photosynthesis, we generated transgenic rice plants overexpressing ZmG1 and ZmG2, individually or simultaneously, with constitutive promoters (pZmUbi::ZmG1 and p35S::ZmG2) or maize promoters (pZmG1::ZmG1, pZmG2::ZmG2, and pZmG1::ZmG1/pZmG2::ZmG2). Both ZmG1 and ZmG2 genes were highly expressed in transgenic rice leaves. Moreover, ZmG1 and ZmG2 showed coordinated expression in pZmG1::ZmG1/pZmG2::ZmG2 plants. All Golden2-like (GLK) transgenic plants had higher chlorophyll and protein contents, Rubisco activities and photosynthetic rates per unit leaf area in flag leaves. However, the highest grain yields occurred when maize promoters were used; pZmG1::ZmG1, pZmG2::ZmG2, and pZmG1::ZmG1/pZmG2::ZmG2 transgenic plants showed increases in grain yield by 51%, 47%, and 70%, respectively. In contrast, the pZmUbi::ZmG1 plant produced smaller seeds without yield increases. Transcriptome analysis indicated that maize GLKs act as master regulators promoting the expression of both photosynthesis-related and stress-responsive regulatory genes in both rice shoot and root. Thus, by promoting these important functions under the control of their own promoters, maize GLK1 and GLK2 genes together dramatically improved rice photosynthetic performance and productivity. A similar approach can potentially improve the productivity of many other crops.

Maize ANT1 modulates vascular development, chloroplast development, photosynthesis, and plant growth.

Arabidopsis AINTEGUMENTA (ANT), an AP2 transcription factor, is known to control plant growth and floral organogenesis. In this study, our transcriptome analysis and in situ hybridization assays of maize embryonic leaves suggested that maize ANT1 (ZmANT1) regulates vascular development. To better understand ANT1 functions, we determined the binding motif of ZmANT1 and then showed that ZmANT1 binds the promoters of millet SCR1, GNC, and AN3, which are key regulators of Kranz anatomy, chloroplast development, and plant growth, respectively. We generated a mutant with a single-codon deletion and two frameshift mutants of the ANT1 ortholog in the C4 millet Setaria viridis by the CRISPR/Cas9 technique. The two frameshift mutants displayed reduced photosynthesis efficiency and growth rate, smaller leaves, and lower grain yields than wild-type (WT) plants. Moreover, their leaves sporadically exhibited distorted Kranz anatomy and vein spacing. Conducting transcriptomic analysis of developing leaves in the WT and the three mutants we identified differentially expressed genes (DEGs) in the two frameshift mutant lines and found many down-regulated DEGs enriched in photosynthesis, heme, tetrapyrrole binding, and antioxidant activity. In addition, we predicted many target genes of ZmANT1 and chose 13 of them to confirm binding of ZmANT1 to their promoters. Based on the above observations, we proposed a model for ANT1 regulation of cell proliferation and leaf growth, vascular and vein development, chloroplast development, and photosynthesis through its target genes. Our study revealed biological roles of ANT1 in several developmental processes beyond its known roles in plant growth and floral organogenesis.

Whole-Genome Duplication Facilitated the Evolution of C4 Photosynthesis in Gynandropsis gynandra.

In higher plants, whole-genome duplication (WGD) is thought to facilitate the evolution of C4 photosynthesis from C3 photosynthesis. To understand this issue, we used new and existing leaf-development transcriptomes to construct two coding sequence databases for C4Gynandropsis gynandra and C3Tarenaya hassleriana, which shared a WGD before their divergence. We compared duplicated genes in the two species and found that the WGD contributed to four aspects of the evolution of C4 photosynthesis in G. gynandra. First, G. gynandra has retained the duplicates of ALAAT (alanine aminotransferase) and GOGAT (glutamine oxoglutarate aminotransferase) for nitrogen recycling to establish a photorespiratory CO2 pump in bundle sheath (BS) cells for increasing photosynthesis efficiency, suggesting that G. gynandra experienced a C3-C4 intermediate stage during the C4 evolution. Second, G. gynandra has retained almost all known vein-development-related paralogous genes derived from the WGD event, likely contributing to the high vein complexity of G. gynandra. Third, the WGD facilitated the evolution of C4 enzyme genes and their recruitment into the C4 pathway. Fourth, several genes encoding photosystem I proteins were derived from the WGD and are upregulated in G. gynandra, likely enabling the NADH dehydrogenase-like complex to produce extra ATPs for the C4 CO2 concentration mechanism. Thus, the WGD apparently played an enabler role in the evolution of C4 photosynthesis in G. gynandra. Importantly, an ALAAT duplicate became highly expressed in BS cells in G. gynandra, facilitating nitrogen recycling and transition to the C4 cycle. This study revealed how WDG may facilitate C4 photosynthesis evolution.

Comparative transcriptomics method to infer gene coexpression networks and its applications to maize and rice leaf transcriptomes.

Time-series transcriptomes of a biological process obtained under different conditions are useful for identifying the regulators of the process and their regulatory networks. However, such data are 3D (gene expression, time, and condition), and there is currently no method that can deal with their full complexity. Here, we developed a method that avoids time-point alignment and normalization between conditions. We applied it to analyze time-series transcriptomes of developing maize leaves under light-dark cycles and under total darkness and obtained eight time-ordered gene coexpression networks (TO-GCNs), which can be used to predict upstream regulators of any genes in the GCNs. One of the eight TO-GCNs is light-independent and likely includes all genes involved in the development of Kranz anatomy, which is a structure crucial for the high efficiency of photosynthesis in C4 plants. Using this TO-GCN, we predicted and experimentally validated a regulatory cascade upstream of SHORTROOT1, a key Kranz anatomy regulator. Moreover, we applied the method to compare transcriptomes from maize and rice leaf segments and identified regulators of maize C4 enzyme genes and RUBISCO SMALL SUBUNIT2 Our study provides not only a powerful method but also novel insights into the regulatory networks underlying Kranz anatomy development and C4 photosynthesis.

Chromosomal-level genome assembly of the semi-dwarf rice Taichung Native 1, an initiator of Green Revolution.

Here we report the 409.5 Mb chromosome-level assembly of the first bred semi-dwarf rice, the Taichung Native 1 (TN1), which served as the template for the development of the Green Revolution (GR) cultivar IR8 "miracle rice". We sequenced the TN1 genome utilizing multiple platforms and produced PacBio long reads, Illumina paired-end reads, Illumina mate-pair reads and 10x Genomics linked reads. We used a hybrid approach to assemble the 226× coverage of sequences by a combination of de novo and reference-guided approaches. The assembled TN1 genome has an N50 scaffold size of 33.1 Mb with the longest measuring 45.5 Mb. We annotated 37,526 genes, in which 24,102 (64.23%) were assigned Blast2GO annotations. The genome has 4672 or 95.4% complete BUSCOs and a repeat content of 51.52%. We developed our own method of creating a GR pangenome using the orthologous relationships of the proteins of TN1, IR8, MH63 and IR64, identifying 16,999 core orthologue groups of Green Revolution. From the pangenome, we identified a set of shared and unique gene ontology terms for the accessory clusters, characterizing TN1, IR8, MH63 and IR64. This TN1 genome assembly and GR pangenome will be a resource for new genomic discoveries about Green Revolution, and for improving the disease and insect resistances and the yield of rice.

Elevated auxin biosynthesis and transport underlie high vein density in C4 leaves.

High vein density, a distinctive trait of C4 leaves, is central to both C3-to-C4 evolution and conversion of C3 to C4-like crops. We tested the hypothesis that high vein density in C4 leaves is due to elevated auxin biosynthesis and transport in developing leaves. Up-regulation of genes in auxin biosynthesis pathways and higher auxin content were found in developing C4 leaves compared with developing C3 leaves. The same observation held for maize foliar (C4) and husk (C3) leaf primordia. Moreover, auxin content and vein density were increased in loss-of-function mutants of Arabidopsis MYC2, a suppressor of auxin biosynthesis. Treatment with an auxin biosynthesis inhibitor or an auxin transport inhibitor led to much fewer veins in new leaves. Finally, both Arabidopsis thaliana auxin efflux transporter pin1 and influx transporter lax2 mutants showed reduced vein numbers. Thus, development of high leaf vein density requires elevated auxin biosynthesis and transport.

Hydrogen cyanamide breaks grapevine bud dormancy in the summer through transient activation of gene expression and accumulation of reactive oxygen and nitrogen species.

BACKGROUND: Hydrogen cyanamide (HC) and pruning (P) have frequently been used to break dormancy in grapevine floral buds. However, the exact underlying mechanism remains elusive. This study aimed to address the early mode of action of these treatments on accumulation of reactive oxygen species (ROS) and reactive nitrogen species (RNS) and expression of related genes in the dormancy breaking buds of grapevine in the summer. RESULTS: The budbreak rates induced by pruning (P), hydrogen cyanamide (HC), pruning plus hydrogen cyanamide (PHC) and water (control) after 8 days were 33, 53, 95, and 0 %, respectively. Clearly, HC was more effective in stimulating grapevine budbreak and P further enhanced its potency. In situ staining of longitudinal bud sections after 12 h of treatments detected high levels of ROS and nitric oxide (NO) accumulated in the buds treated with PHC, compared with HC or P alone. The amounts of ROS and NO accumulated were highly correlated with the rates of budbreak among these treatments, highlighting the importance of a rapid, transient accumulation of sublethal levels of ROS and RNS in dormancy breaking. Microarray analysis revealed specific alterations in gene expression in dormancy breaking buds induced by P, HC and PHC after 24 h of treatment. Relative to control, PHC altered the expression of the largest number of genes, while P affected the expression of the least number of genes. PHC also exerted a greater intensity in transcriptional activation of these genes. Gene ontology (GO) analysis suggests that alteration in expression of ROS related genes is the major factor responsible for budbreak. qRT-PCR analysis revealed the transient expression dynamics of 12 specific genes related to ROS generation and scavenge during the 48 h treatment with PHC. CONCLUSION: Our results suggest that rapid accumulation of ROS and NO at early stage is important for dormancy release in grapevine in the summer, and the identification of the commonly expressed specific genes among the treatments allowed the construction of the signal transduction pathway related to ROS/RNS metabolism during dormancy release. The rapid accumulation of a sublethal level of ROS/RNS subsequently induces cell wall loosening and expansion for bud sprouting.