Seasonal pigment fluctuation in diploid and polyploid Arabidopsis revealed by machine learning-based phenotyping method PlantServation.

Long-term field monitoring of leaf pigment content is informative for understanding plant responses to environments distinct from regulated chambers but is impractical by conventional destructive measurements. We developed PlantServation, a method incorporating robust image-acquisition hardware and deep learning-based software that extracts leaf color by detecting plant individuals automatically. As a case study, we applied PlantServation to examine environmental and genotypic effects on the pigment anthocyanin content estimated from leaf color. We processed >4 million images of small individuals of four Arabidopsis species in the field, where the plant shape, color, and background vary over months. Past radiation, coldness, and precipitation significantly affected the anthocyanin content. The synthetic allopolyploid A. kamchatica recapitulated the fluctuations of natural polyploids by integrating diploid responses. The data support a long-standing hypothesis stating that allopolyploids can inherit and combine the traits of progenitors. PlantServation facilitates the study of plant responses to complex environments termed "in natura".

A low-coverage 3' RNA-seq to detect homeolog expression in polyploid wheat.

Although allopolyploid species are common among natural and crop species, it is not easy to distinguish duplicated genes, known as homeologs, during their genomic analysis. Yet, cost-efficient RNA sequencing (RNA-seq) is to be developed for large-scale transcriptomic studies such as time-series analysis and genome-wide association studies in allopolyploids. In this study, we employed a 3' RNA-seq utilizing 3' untranslated regions (UTRs) containing frequent mutations among homeologous genes, compared to coding sequence. Among the 3' RNA-seq protocols, we examined a low-cost method Lasy-Seq using an allohexaploid bread wheat, Triticum aestivum. HISAT2 showed the best performance for 3' RNA-seq with the least mapping errors and quick computational time. The number of detected homeologs was further improved by extending 1 kb of the 3' UTR annotation. Differentially expressed genes in response to mild cold treatment detected by the 3' RNA-seq were verified with high-coverage conventional RNA-seq, although the latter detected more differentially expressed genes. Finally, downsampling showed that even a 2 million sequencing depth can still detect more than half of expressed homeologs identifiable by the conventional 32 million reads. These data demonstrate that this low-cost 3' RNA-seq facilitates large-scale transcriptomic studies of allohexaploid wheat and indicate the potential application to other allopolyploid species.

WUSCHEL-RELATED HOMEOBOX 13 suppresses de novo shoot regeneration via cell fate control of pluripotent callus.

Plants can regenerate their bodies via de novo establishment of shoot apical meristems (SAMs) from pluripotent callus. Only a small fraction of callus cells is eventually specified into SAMs but the molecular mechanisms underlying fate specification remain obscure. The expression of WUSCHEL (WUS) is an early hallmark of SAM fate acquisition. Here, we show that a WUS paralog, WUSCHEL-RELATED HOMEOBOX 13 (WOX13), negatively regulates SAM formation from callus in Arabidopsis thaliana. WOX13 promotes non-meristematic cell fate via transcriptional repression of WUS and other SAM regulators and activation of cell wall modifiers. Our Quartz-Seq2-based single cell transcriptome revealed that WOX13 plays key roles in determining cellular identity of callus cell population. We propose that reciprocal inhibition between WUS and WOX13 mediates critical cell fate determination in pluripotent cell population, which has a major impact on regeneration efficiency.

Targeted single-cell gene induction by optimizing the dually regulated CRE/loxP system by a newly defined heat-shock promoter and the steroid hormone in Arabidopsis thaliana.

Multicellular organisms rely on intercellular communication systems to organize their cellular functions. In studies focusing on intercellular communication, the key experimental techniques include the generation of chimeric tissue using transgenic DNA recombination systems represented by the CRE/loxP system. If an experimental system enables the induction of chimeras at highly targeted cell(s), it will facilitate the reproducibility and precision of experiments. However, multiple technical limitations have made this challenging. The stochastic nature of DNA recombination events, especially, hampers reproducible generation of intended chimeric patterns. Infrared laser-evoked gene operator (IR-LEGO), a microscopic system that irradiates targeted cells using an IR laser, can induce heat shock-mediated expression of transgenes, for example, CRE recombinase gene, in the cells. In this study, we developed a method that induces CRE/loxP recombination in the target cell(s) of plant roots and leaves in a highly specific manner. We combined IR-LEGO, an improved heat-shock-specific promoter, and dexamethasone-dependent regulation of CRE. The optimal IR-laser power and irradiation duration were estimated via exhaustive irradiation trials and subsequent statistical modeling. Under optimized conditions, CRE/loxP recombination was efficiently induced without cellular damage. We also found that the induction efficiency varied among tissue types and cellular sizes. The developed method offers an experimental system to generate a precisely designed chimeric tissue, and thus, will be useful for analyzing intercellular communication at high resolution in roots and leaves.

Auxin-Induced WUSCHEL-RELATED HOMEOBOX13 Mediates Asymmetric Activity of Callus Formation upon Cutting.

Plants have the regenerative ability to reconnect cut organs, which is physiologically important to survive severe tissue damage. The ability to reconnect organs is utilized as grafting to combine two different individuals. Callus formation at the graft junction facilitates organ attachment and vascular reconnection. While it is well documented that local wounding signals provoke callus formation, how callus formation is differentially regulated at each cut end remains elusive. Here, we report that callus formation activity is asymmetrical between the top and bottom cut ends and is regulated by differential auxin accumulation. Gene expression analyses revealed that cellular auxin response is preferentially upregulated in the top part of the graft. Disruption of polar auxin transport inhibited callus formation from the top, while external application of auxin was sufficient to induce callus formation from the bottom, suggesting that asymmetric auxin accumulation is responsible for active callus formation from the top end. We further found that the expression of a key regulator of callus formation, WUSCHEL-RELATED HOMEOBOX 13 (WOX13), is induced by auxin. The ectopic callus formation from the bottom end, which is triggered by locally supplemented auxin, requires WOX13 function, demonstrating that WOX13 plays a pivotal role in auxin-dependent callus formation. The asymmetric WOX13 expression is observed both in grafted petioles and incised inflorescence stems, underscoring the generality of our findings. We propose that efficient organ reconnection is achieved by a combination of local wounding stimuli and disrupted long-distance signaling.

Development of a high-throughput field phenotyping rover optimized for size-limited breeding fields as open-source hardware.

Phenotyping is a critical process in plant breeding, especially when there is an increasing demand for streamlining a selection process in a breeding program. Since manual phenotyping has limited efficiency, high-throughput phenotyping methods are recently popularized owing to progress in sensor and image processing technologies. However, in a size-limited breeding field, which is common in Japan and other Asian countries, it is challenging to introduce large machinery in the field or fly unmanned aerial vehicles over the field. In this study, we developed a ground-based high-throughput field phenotyping rover that could be easily introduced to a field regardless of the scale and location of the field even without special facilities. We also made the field rover open-source hardware, making its system available to public for easy modification, so that anyone can build one for their own use at a low cost. The trial run of the field rover revealed that it allowed the collection of detailed remote-sensing images of plants and quantitative analyses based on the images. The results suggest that the field rover developed in this study could allow efficient phenotyping of plants especially in a small breeding field.

Wound-inducible WUSCHEL-RELATED HOMEOBOX 13 is required for callus growth and organ reconnection.

Highly efficient tissue repair is pivotal for surviving damage-associated stress. Plants generate callus upon injury to heal wound sites, yet regulatory mechanisms of tissue repair remain elusive. Here, we identified WUSCHEL-RELATED HOMEOBOX 13 (WOX13) as a key regulator of callus formation and organ adhesion in Arabidopsis (Arabidopsis thaliana). WOX13 belongs to an ancient subclade of the WOX family, and a previous study shows that WOX13 orthologs in the moss Physcomitrium patens (PpWOX13L) are involved in cellular reprogramming at wound sites. We found that the Arabidopsis wox13 mutant is totally defective in establishing organ reconnection upon grafting, suggesting that WOX13 is crucial for tissue repair in seed plants. WOX13 expression rapidly induced upon wounding, which was partly dependent on the activity of an AP2/ERF transcription factor, WOUND-INDUCED DEDIFFERENTIATION 1 (WIND1). WOX13 in turn directly upregulated WIND2 and WIND3 to further promote cellular reprogramming and organ regeneration. We also found that WOX13 orchestrates the transcriptional induction of cell wall-modifying enzyme genes, such as GLYCOSYL HYDROLASE 9Bs, PECTATE LYASE LIKEs and EXPANSINs. Furthermore, the chemical composition of cell wall monosaccharides was markedly different in the wox13 mutant. These data together suggest that WOX13 modifies cell wall properties, which may facilitate efficient callus formation and organ reconnection. Furthermore, we found that PpWOX13L complements the Arabidopsis wox13 mutant, suggesting that the molecular function of WOX13 is partly conserved between mosses and seed plants. This study provides key insights into the conservation and functional diversification of the WOX gene family during land plant evolution.

The boundary-expressed EPIDERMAL PATTERNING FACTOR-LIKE2 gene encoding a signaling peptide promotes cotyledon growth during Arabidopsis thaliana embryogenesis.

The shoot organ boundaries have important roles in plant growth and morphogenesis. It has been reported that a gene encoding a cysteine-rich secreted peptide of the EPIDERMAL PATTERNING FACTOR-LIKE (EPFL) family, EPFL2, is expressed in the boundary domain between the two cotyledon primordia of Arabidopsis thaliana embryo. However, its developmental functions remain unknown. This study aimed to analyze the role of EPFL2 during embryogenesis. We found that cotyledon growth was reduced in its loss-of-function mutants, and this phenotype was associated with the reduction of auxin response peaks at the tips of the primordia. The reduced cotyledon size of the mutant embryo recovered in germinating seedlings, indicating the presence of a factor that acted redundantly with EPFL2 to promote cotyledon growth in late embryogenesis. Our analysis suggests that the boundary domain between the cotyledon primordia acts as a signaling center that organizes auxin response peaks and promotes cotyledon growth.

Corrigendum: An Affordable Image-Analysis Platform to Accelerate Stomatal Phenotyping During Microscopic Observation.

[This corrects the article DOI: 10.3389/fpls.2021.715309.].

An Affordable Image-Analysis Platform to Accelerate Stomatal Phenotyping During Microscopic Observation.

Recent technical advances in the computer-vision domain have facilitated the development of various methods for achieving image-based quantification of stomata-related traits. However, the installation cost of such a system and the difficulties of operating it on-site have been hurdles for experimental biologists. Here, we present a platform that allows real-time stomata detection during microscopic observation. The proposed system consists of a deep neural network model-based stomata detector and an upright microscope connected to a USB camera and a graphics processing unit (GPU)-supported single-board computer. All the hardware components are commercially available at common electronic commerce stores at a reasonable price. Moreover, the machine-learning model is prepared based on freely available cloud services. This approach allows users to set up a phenotyping platform at low cost. As a proof of concept, we trained our model to detect dumbbell-shaped stomata from wheat leaf imprints. Using this platform, we collected a comprehensive range of stomatal phenotypes from wheat leaves. We confirmed notable differences in stomatal density (SD) between adaxial and abaxial surfaces and in stomatal size (SS) between wheat-related species of different ploidy. Utilizing such a platform is expected to accelerate research that involves all aspects of stomata phenotyping.

De Novo Genome Assembly of the Japanese Wheat Cultivar Norin 61 Highlights Functional Variation in Flowering Time and Fusarium-Resistant Genes in East Asian Genotypes.

Bread wheat is a major crop that has long been the focus of basic and breeding research. Assembly of its genome has been difficult because of its large size and allohexaploid nature (AABBDD genome). Following the first reported assembly of the genome of the experimental strain Chinese Spring (CS), the 10+ Wheat Genomes Project was launched to produce multiple assemblies of worldwide modern cultivars. The only Asian cultivar in the project is Norin 61, a representative Japanese cultivar adapted to grow across a broad latitudinal range, mostly characterized by a wet climate and a short growing season. Here, we characterize the key aspects of its chromosome-scale genome assembly spanning 15 Gb with a raw scaffold N50 of 22 Mb. Analysis of the repetitive elements identified chromosomal regions unique to Norin 61 that encompass a tandem array of the pathogenesis-related 13 family. We report novel copy-number variations in the B homeolog of the florigen gene FT1/VRN3, pseudogenization of its D homeolog and the association of its A homeologous alleles with the spring/winter growth habit. Furthermore, the Norin 61 genome carries typical East Asian functional variants different from CS, ranging from a single nucleotide to multi-Mb scale. Examples of such variation are the Fhb1 locus, which confers Fusarium head-blight resistance, Ppd-D1a, which confers early flowering, Glu-D1f for Asian noodle quality and Rht-D1b, which introduced semi-dwarfism during the green revolution. The adoption of Norin 61 as a reference assembly for functional and evolutionary studies will enable comprehensive characterization of the underexploited Asian bread wheat diversity.

Homeolog expression quantification methods for allopolyploids.

Genome duplication with hybridization, or allopolyploidization, occurs in animals, fungi and plants, and is especially common in crop plants. There is an increasing interest in the study of allopolyploids because of advances in polyploid genome assembly; however, the high level of sequence similarity in duplicated gene copies (homeologs) poses many challenges. Here we compared standard RNA-seq expression quantification approaches used currently for diploid species against subgenome-classification approaches which maps reads to each subgenome separately. We examined mapping error using our previous and new RNA-seq data in which a subgenome is experimentally added (synthetic allotetraploid Arabidopsis kamchatica) or reduced (allohexaploid wheat Triticum aestivum versus extracted allotetraploid) as ground truth. The error rates in the two species were very similar. The standard approaches showed higher error rates (>10% using pseudo-alignment with Kallisto) while subgenome-classification approaches showed much lower error rates (<1% using EAGLE-RC, <2% using HomeoRoq). Although downstream analysis may partly mitigate mapping errors, the difference in methods was substantial in hexaploid wheat, where Kallisto appeared to have systematic differences relative to other methods. Only approximately half of the differentially expressed homeologs detected using Kallisto overlapped with those by any other method in wheat. In general, disagreement in low-expression genes was responsible for most of the discordance between methods, which is consistent with known biases in Kallisto. We also observed that there exist uncertainties in genome sequences and annotation which can affect each method differently. Overall, subgenome-classification approaches tend to perform better than standard approaches with EAGLE-RC having the highest precision.

The functional balance between the WUSCHEL-RELATED HOMEOBOX1 gene and the phytohormone auxin is a key factor for cell proliferation in Arabidopsis seedlings.

The WUSCHEL-RELATED HOMEOBOX1 (WOX1) transcription factor and its homolog PRESSED FLOWER (PRS) are multifunctional regulators of leaf development that act as transcriptional repressors. These genes promote cell proliferation under certain conditions, but the related molecular mechanisms are not well understood. Here, we present a new function for WOX1 in cell proliferation. To identify the WOX1 downstream genes, we performed a microarray analysis of shoot apices of transgenic Arabidopsis thaliana lines harboring [35Sp::WOX1-glucocorticoid receptor (GR)] in which the WOX1 function was temporarily enhanced by dexamethasone. The downregulated genes were significantly enriched for the Gene Ontology term "response to auxin stimulus", whereas the significantly upregulated genes contained auxin transport-associated PIN1 and AUX1 and the auxin response factor MP, which are involved in formation of auxin response maxima. Simultaneous treatments of synthetic auxin and dexamethasone induced the formation of green compact calli and the unorganized proliferation of cells in the hypocotyl. A microarray analysis of 35Sp::WOX1-GR plants treated with indole-3-acetic acid and dexamethasone revealed that WOX1 and auxin additively influenced their common downstream genes. Furthermore, in the presence of an auxin-transport inhibitor, cell proliferation during leaf initiation was suppressed in the prs mutant but induced in a broad region of the peripheral zone of the shoot apical meristem in the ectopic WOX1-expressing line FILp::WOX1. Thus, our results clarify the additive effect of WOX1/PRS and auxin on their common downstream genes and highlight the importance of the balance between their functions in controlling cell proliferation.

Stem development through vascular tissues: EPFL-ERECTA family signaling that bounces in and out of phloem.

Plant cells communicate with each other using a variety of signaling molecules. Recent studies have revealed that various types of secreted peptides, as well as phytohormones known since long ago, mediate cell-cell communication in diverse contexts of plant life. These peptides affect cellular activities, such as proliferation and cell fate decisions, through their perception by cell surface receptors located on the plasma membrane of target cells. ERECTA (ER), an Arabidopsis thaliana receptor kinase gene, was first identified as a stem growth regulator, and since then an increasing number of studies have shown that ER is involved in a wide range of developmental and physiological processes. In particular, molecular functions of ER have been extensively studied in stomatal patterning. Furthermore, the importance of ER signaling in vascular tissues of inflorescence stems, especially in phloem cells, has recently been highlighted. In this review article, first we briefly summarize the history of ER research including studies on stomatal development, then introduce ER functions in vascular tissues, and discuss its interactions with phytohormones and other receptor kinase signaling pathways. Future questions and challenges will also be addressed.

Impact of erecta mutation on leaf serration differs between Arabidopsis accessions.

Serrations or teeth of plant leaves are a morphological trait regulated genetically and environmentally. Very recently, it has been reported that the receptor kinases encoded by three ERECTA (ER)-family genes, ER, ER-LIKE1 (ERL1) and ERL2, redundantly play a role in tooth growth in Arabidopsis thaliana. In the report, Columbia (Col) accession was used for analyses, where none of the signal mutant of the ER-family genes exhibited serration defects. The toothless, smooth leaf margin phenotype was evident only when two out of the three ER-family genes were lost. Interestingly, it has been widely recognized that the Arabidopsis accession Landsberg erecta (L.er), which carries a loss-of-function mutation in ER, develops round leaves with smaller leaf teeth. Here, we show that the functional ER transgene promotes the tooth growth in L.er to the level of Col, indicating that the er mutation in L.er is likely responsible for the reduced growth of leaf teeth. This suggests that er single mutation affects tooth growth in a different manner between Col and L.er backgrounds, though the molecular basis for this background-dependent effect remains to be addressed.

A Secreted Peptide and Its Receptors Shape the Auxin Response Pattern and Leaf Margin Morphogenesis.

Secreted peptides mediate intercellular communication [1, 2]. Several secreted peptides in the EPIDERMAL PATTERNING FACTOR-LIKE (EPFL) family regulate morphogenesis of tissues, such as stomata and inflorescences in plants [3-15]. The biological functions of other EPFL family members remain unknown. Here, we show that the EPFL2 gene is required for growth of leaf teeth. EPFL2 peptide physically interacts with ERECTA (ER) family receptor-kinases and, accordingly, the attenuation of ER family activities leads to formation of toothless leaves. During the tooth growth process, responses to the phytohormone auxin are maintained at tips of the teeth to promote their growth [16-19]. In the growing tooth tip of epfl2 and multiple er family mutants, the auxin response becomes broader. Conversely, overexpression of EPFL2 diminishes the auxin response, indicating that the EPFL2 signal restricts the auxin response to the tooth tip. Interestingly, the tip-specific auxin response in turn organizes characteristic expression patterns of ER family and EPFL2 by enhancing ER family expression at the tip while eliminating the EPFL2 expression from the tip. Our findings identify the novel ligand-receptor pairs promoting the tooth growth, and further reveal a feedback circuit between the peptide-receptor system and auxin response as a mechanism for maintaining proper auxin maxima during leaf margin morphogenesis.

Cell walls as a stage for intercellular communication regulating shoot meristem development.

Aboveground organs of plants are ultimately derived/generated from the shoot apical meristem (SAM), which is a proliferative tissue located at the apex of the stem. The SAM contains a population of stem cells that provide new cells for organ/tissue formation. The SAM is composed of distinct cell layers and zones with different properties. Primordia of lateral organs develop at the periphery of the SAM. The shoot apex is a dynamic and complex tissue, and as such intercellular communications among cells, layers and zones play significant roles in the coordination of cell proliferation, growth and differentiation to achieve elaborate morphogenesis. Recent findings have highlighted the importance of a number of signaling molecules acting in the cell wall space for the intercellular communication, including classic phytohormones and secretory peptides. Moreover, accumulating evidence has revealed that cell wall properties and their modifying enzymes modulate hormone actions. In this review, we outline how behaviors of signaling molecules and changes of cell wall properties are integrated for the shoot meristem regulation.

Pattern dynamics in adaxial-abaxial specific gene expression are modulated by a plastid retrograde signal during Arabidopsis thaliana leaf development.

The maintenance and reformation of gene expression domains are the basis for the morphogenic processes of multicellular systems. In a leaf primordium of Arabidopsis thaliana, the expression of FILAMENTOUS FLOWER (FIL) and the activity of the microRNA miR165/166 are specific to the abaxial side. This miR165/166 activity restricts the target gene expression to the adaxial side. The adaxial and abaxial specific gene expressions are crucial for the wide expansion of leaf lamina. The FIL-expression and the miR165/166-free domains are almost mutually exclusive, and they have been considered to be maintained during leaf development. However, we found here that the position of the boundary between the two domains gradually shifts from the adaxial side to the abaxial side. The cell lineage analysis revealed that this boundary shifting was associated with a sequential gene expression switch from the FIL-expressing (miR165/166 active) to the miR165/166-free (non-FIL-expressing) states. Our genetic analyses using the enlarged fil expression domain2 (enf2) mutant and chemical treatment experiments revealed that impairment in the plastid (chloroplast) gene expression machinery retards this boundary shifting and inhibits the lamina expansion. Furthermore, these developmental effects caused by the abnormal plastids were not observed in the genomes uncoupled1 (gun1) mutant background. This study characterizes the dynamic nature of the adaxial-abaxial specification process in leaf primordia and reveals that the dynamic process is affected by the GUN1-dependent retrograde signal in response to the failure of plastid gene expression. These findings advance our understanding on the molecular mechanism linking the plastid function to the leaf morphogenic processes.

Succinic semialdehyde dehydrogenase is involved in the robust patterning of Arabidopsis leaves along the adaxial-abaxial axis.

Polarity along the adaxial-abaxial axis of the leaf is essential for leaf development and morphogenesis. One of the genes that encodes a putative transcription factor regulating adaxial-abaxial polarity, FILAMENTOUS FLOWER (FIL), is expressed in the abaxial region of the leaf primordia. However, the molecular mechanisms controlling the polarized expression of FIL remain unclear. Here, we analyzed an enlarged fil expression domain1 (enf1) mutant of Arabidopsis, which forms both abaxialized leaves and adaxialized leaves. The ENF1 gene encodes SUCCINIC SEMIALDEHYDE DEHYDROGENASE (SSADH), which catalyzes the conversion of succinic semialdehyde (SSA) to succinate. The enf1 phenotype was suppressed by an additional mutation in GAMMA-AMINOBUTYRIC ACID AMINOTRANSFERASE1 (GABAT1), which encodes an SSA-producing enzyme, suggesting that SSA or its derivatives is the metabolite responsible for the defect in the adaxial-abaxial axis-dependent gene expression of enf1. In the shoot apical meristem, GABAT1 was expressed in the outermost layer but SSADH was not. Exogenous application of SSA induced adaxial characters on the abaxial side of the newly developed leaves. We suggest that a GABA shunt metabolite, SSA or its close derivatives, is involved in the robust leaf patterning and structure along the adaxial-abaxial axis.