Functional Redundancy and Divergence within the Arabidopsis RETICULATA-RELATED Gene Family

Proteomics- and metabolomics-based analysis of the regulation of germination in Norway maple and sycamore embryonic axes

Norway maple and sycamore belong to the Acer genus and produce desiccation-tolerant and desiccation-sensitive seeds, respectively. We investigated the seed germination process at the imbibed and germinated stages using metabolomic and proteomic approaches to determine why sycamore seeds germinate earlier and are more successful at establishing seedlings than Norway maple seeds under controlled conditions. Embryonic axes and embryonic axes with protruded radicles were analyzed at the imbibed and germinated stages, respectively. Among the 212 identified metabolites, 44 and 67 differentially abundant metabolites were found at the imbibed and germinated stages, respectively, in both Acer species. Higher levels of amines, growth and defense stimulants, including B vitamins, were found in sycamore. We identified 611 and 447 proteins specific to the imbibed and germinated stages, respectively, in addition to groups of proteins expressed at different levels. Functional analysis of significantly regulated proteins revealed that proteins with catalytic and binding activity were enriched during germination, and proteins possibly implicated in nitrogen metabolism and metabolite interconversion enzymes were the predominant classes. Proteins associated with the control of plant growth regulation and seed defense were observed in both species at both germination stages. Sycamore proteins possibly involved in abscisic acid signal transduction pathway, stress tolerance and alleviation, ion binding and oxygenase activities appeared to accompany germination in sycamore. We identified peptides containing methionine (Met) oxidized to methionine sulfoxide (MetO), and functional analyses of proteins with significantly regulated MetO sites revealed that translation, plant growth and development and metabolism of nitrogen compounds were the main processes under Met/MetO redox control. We propose that higher levels of storage proteins and amines, together with higher levels of B vitamins, supported more efficient nitrogen utilization in sycamore, resulting in faster seedling growth. In conclusion, omic signatures identified in sycamore seem to predispose germinated sycamore seeds to better postgerminative growth.

Functional conservation and divergence of arabidopsis VENOSA4 and human SAMHD1 in DNA repair

The human deoxyribonucleoside triphosphatase (dNTPase) Sterile alpha motif and histidine-aspartate domain containing protein 1 (SAMHD1) has a dNTPase-independent role in repairing DNA double-strand breaks (DSBs) by homologous recombination (HR). Here, we show that VENOSA4 (VEN4), the probable Arabidopsis thaliana ortholog of SAMHD1, also functions in DSB repair by HR. The ven4 loss-of-function mutants showed increased DNA ploidy and deregulated DNA repair genes, suggesting DNA damage accumulation. Hydroxyurea, which blocks DNA replication and generates DSBs, induced VEN4 expression. The ven4 mutants were hypersensitive to hydroxyurea, with decreased DSB repair by HR. Metabolomic analysis of the strong ven4-0 mutant revealed depletion of metabolites associated with DNA damage responses. In contrast to SAMHD1, VEN4 showed no evident involvement in preventing R-loop accumulation. Our study thus reveals functional conservation in DNA repair by VEN4 and SAMHD1.

Identification of the AHP family reveals their critical response to cytokinin regulation during adventitious root formation in apple rootstock

Adventitious root (AR) formation is a bottleneck for vegetative proliferation. In this study, 13 AHP genes (MdAHPs) were identified in the apple genome. Phylogenetic analysis grouped them into 3 clusters (I, II, III), with 4, 4, and 5 genes respectively. The 13 MdAHPs family members were named MdAHP1 to MdAHP13 by chromosome positions. The physicochemical properties, phylogenetic relationship, motifs, and elements of their proteins were also analyzed. The amino acid quantity varied from 60~189 aa, isoelectric point lay between 4.10 and 8.93, and there were 3~7 protein-conserving motifs. Excluding MdAHP6, other members' promoter sequences behaved 2-4 CTK response elements. Additionally, the expression characteristics of MdAHPs family members at key stages of AR formation and in different tissues were also examined with exogenous 6-BA and Lov treatments. The results showed that MdAHP3 might be a key member in AR formation. GUS staining indicated that the activity of the MdAHP3 promoter was also significantly enhanced by CTK treatment. The protein interactions of MdAHP3/MdAHP1 and MdAHP3/MdAHP6 were verified. Compared with WT, 35S::MdAHP3 transgenic poplars inhibited AR formation. The above experimental results suggested that MdAHP3, as a key family member, interacts with MdAHP1 and MdAHP6 proteins to jointly mediate AR formation in apple rootstocks.

A Comprehensive Analysis In Silico of KCS Genes in Maize Revealed Their Potential Role in Response to Abiotic Stress

β-ketoacyl-CoA synthase (KCS) enzymes play a pivotal role in plants by catalyzing the first step of very long-chain fatty acid (VLCFA) biosynthesis. This process is crucial for plant development and stress responses. However, the understanding of KCS genes in maize remains limited. In this study, we present a comprehensive analysis of ZmKCS genes, identifying 29 KCS genes that are unevenly distributed across nine maize chromosomes through bioinformatics approaches. These ZmKCS proteins varied in length and molecular weight, suggesting functional diversity. Phylogenetic analysis categorized 182 KCS proteins from seven species into six subgroups, with maize showing a closer evolutionary relationship to other monocots. Collinearity analysis revealed 102 gene pairs between maize and three other monocots, whereas only five gene pairs were identified between maize and three dicots, underscoring the evolutionary divergence of KCS genes between monocotyledonous and dicotyledonous plants. Structural analysis revealed that 20 out of 29 ZmKCS genes are intronless. Subcellular localization prediction and experimental validation suggest that most ZmKCS proteins are likely localized at the plasma membrane, with some also present in mitochondria and chloroplasts. Analysis of the cis-acting elements within the ZmKCS promoters suggested their potential involvement in abiotic stress responses. Notably, expression analysis under abiotic stresses highlighted ZmKCS17 as a potential key gene in the stress response of maize, which presented an over 10-fold decrease in expression under salt and drought stresses within 48 h. This study provides a fundamental understanding of ZmKCS genes, paving the way for further functional characterization and their potential application in maize breeding for enhanced stress tolerance.

Novel SNP markers and other stress-related genomic regions associated with nitrogen use efficiency in cassava

Cassava productivity is constrained by low soil nitrogen, which is predominant in most cassava-growing regions in the tropics and subtropical agroecology. Improving the low nitrogen tolerance of cassava has become an important breeding objective. The current study aimed to develop cassava varieties with improved nitrogen use efficiency by identifying genomic regions and candidate genes linked to nitrogen use efficiency in cassava. A genome-wide association study (GWAS) was performed using the Genome Association and Prediction Integrated Tool (GAPIT). A panel of 265 diverse cassava genotypes was phenotyped for 10 physiological and agronomic traits under optimum and low-nitrogen regimes. Whole-genome genotyping of these cassava cloneswas performed using the Diversity Arrays Technology (DArTseq) sequencing platform. A total of 68,814 single nucleotide polymorphisms (SNPs) were identified, which were spread across the entire 18 chromosomes of the cassava genome, of which 52 SNPs at various densities were found to be associated with nitrogen use efficiency in cassava and other yield-related traits. The putative genes identified through GWAS, especially those with significant associated SNP markers for NUE and related traits have the potential, if deployed appropriately, to develop cassava varieties with improved nitrogen use efficiency, which would translate to a reduction in the economic and environmental cost of cassava production.

Mapping the castor bean endosperm proteome revealed a metabolic interaction between plastid, mitochondria, and peroxisomes to optimize seedling growth

In this work, we studied castor-oil plant Ricinus communis as a classical system for endosperm reserve breakdown. The seeds of castor beans consist of a centrally located embryo with the two thin cotyledons surrounded by the endosperm. The endosperm functions as major storage tissue and is packed with nutritional reserves, such as oil, proteins, and starch. Upon germination, mobilization of the storage reserves requires inter-organellar interplay of plastids, mitochondria, and peroxisomes to optimize growth for the developing seedling. To understand their metabolic interactions, we performed a large-scale organellar proteomic study on castor bean endosperm. Organelles from endosperm of etiolated seedlings were isolated and subjected to liquid chromatography-tandem mass spectrometry (LC-MS/MS). Computer-assisted deconvolution algorithms were applied to reliably assign the identified proteins to their correct subcellular localization and to determine the abundance of the different organelles in the heterogeneous protein samples. The data obtained were used to build a comprehensive metabolic model for plastids, mitochondria, and peroxisomes during storage reserve mobilization in castor bean endosperm.

Analysis of the Arabidopsis venosa4-0 mutant supports the role of VENOSA4 in dNTP metabolism

Human Sterile alpha motif and histidine-aspartate domain containing protein 1 (SAMHD1) functions as a dNTPase to maintain dNTP pool balance. In eukaryotes, the limiting step in de novo dNTP biosynthesis is catalyzed by RIBONUCLEOTIDE REDUCTASE (RNR). In Arabidopsis, the RNR1 subunit of RNR is encoded by CRINKLED LEAVES 8 (CLS8), and RNR2 by three paralogous genes, including TSO MEANING 'UGLY' IN CHINESE 2 (TSO2). In plants, DIFFERENTIAL DEVELOPMENT OF VASCULAR ASSOCIATED CELLS 1 (DOV1) catalyzes the first step of the de novo biosynthesis of purines. Here, to explore the role of VENOSA4 (VEN4), the most likely Arabidopsis ortholog of human SAMHD1, we studied the ven4-0 point mutation, whose leaf phenotype was stronger than those of its insertional alleles. Structural predictions suggested that the E249L substitution in the mutated VEN4-0 protein rigidifies its 3D structure. The morphological phenotypes of the ven4, cls8, and dov1 single mutants were similar, and those of the ven4 tso2 and ven4 dov1 double mutants were synergistic. The ven4-0 mutant had reduced levels of four amino acids related to dNTP biosynthesis, including glutamine and glycine, which are precursors in the de novo purine biosynthesis. Our results reveal high functional conservation between VEN4 and SAMHD1 in dNTP metabolism.

Plant environmental sensing relies on specialized plastids

In plants, plastids are thought to interconvert to various forms that are specialized for photosynthesis, starch and oil storage, and diverse pigment accumulation. Post-endosymbiotic evolution has led to adaptations and specializations within plastid populations that align organellar functions with different cellular properties in primary and secondary metabolism, plant growth, organ development, and environmental sensing. Here, we review the plastid biology literature in light of recent reports supporting a class of 'sensory plastids' that are specialized for stress sensing and signaling. Abundant literature indicates that epidermal and vascular parenchyma plastids display shared features of dynamic morphology, proteome composition, and plastid-nuclear interaction that facilitate environmental sensing and signaling. These findings have the potential to reshape our understanding of plastid functional diversification.

Three Diverse Granule Preparation Methods for Proteomic Analysis of Mature Rice (Oryza sativa L.) Starch Grain

Starch is the primary form of reserve carbohydrate storage in plants. Rice (Oryza sativa L.) is a monocot whose reserve starch is organized into compounded structures within the amyloplast, rather than a simple starch grain (SG). The mechanism governing the assembly of the compound SG from polyhedral granules in apposition, however, remains unknown. To further characterize the proteome associated with these compounded structures, three distinct methods of starch granule preparation (dispersion, microsieve, and flotation) were performed. Phase separation of peptides (aqueous trypsin-shaving and isopropanol solubilization of residual peptides) isolated starch granule-associated proteins (SGAPs) from the distal proteome of the amyloplast and the proximal ‘amylome’ (the amyloplastic proteome), respectively. The term ‘distal proteome’ refers to SGAPs loosely tethered to the amyloplast, ones that can be rapidly proteolyzed, while proximal SGAPs are those found closer to the remnant amyloplast membrane fragments, perhaps embedded therein—ones that need isopropanol solvent to be removed from the mature organelle surface. These two rice starch-associated peptide samples were analyzed using nano-liquid chromatography–tandem mass spectrometry (Nano-HPLC-MS/MS). Known and novel proteins, as well as septum-like structure (SLS) proteins, in the mature rice SG were found. Data mining and gene ontology software were used to categorize these putative plastoskeletal components as a variety of structural elements, including actins, tubulins, tubulin-like proteins, and cementitious elements such as reticulata related-like (RER) proteins, tegument proteins, and lectins. Delineating the plastoskeletal proteome begins by understanding how each starch granule isolation procedure affects observed cytoplasmic and plastid proteins. The three methods described herein show how the technique used to isolate SGs differentially impacts the subsequent proteomic analysis and results obtained. It can thus be concluded that future investigations must make judicious decisions regarding the methodology used in extracting proteomic information from the compound starch granules being assessed, since different methods are shown to yield contrasting results herein. Data are available via ProteomeXchange with identifier PXD032314.

Haplotype analyses reveal novel insights into tomato history and domestication driven by long-distance migrations and latitudinal adaptations

A novel haplotype-based approach that uses Procrustes analysis and automatic classification was used to provide further insights into tomato history and domestication. Agrarian societies domesticated species of interest by introducing complex genetic modifications. For tomatoes, two species, one of which had two botanical varieties, are thought to be involved in its domestication: the fully wild Solanum pimpinellifolium (SP), the wild and semi-domesticated Solanum lycopersicum var. cerasiforme (SLC) and the cultivated S. l. var. lycopersicum (SLL). The Procrustes approach showed that SP evolved into SLC during a gradual migration from the Peruvian deserts to the Mexican rainforests and that Peruvian and Ecuadorian SLC populations were the result of more recent hybridizations. Our model was supported by independent evidence, including ecological data from the accession collection site and morphological data. Furthermore, we showed that photosynthesis-, and flowering time-related genes were selected during the latitudinal migrations.

Genetic Factors Underlying Single Fiber Quality in A-Genome Donor Asian Cotton (Gossypium arboreum)

The study of A-genome Asian cotton as a potential fiber donor in Gossypium species may offer an enhanced understanding of complex genetics and novel players related to fiber quality traits. Assessment of individual fibers providing classified fiber quality information to the textile industry is Advanced Fiber Information System (AFIS) in the recent technological era. Keeping the scenario, a diverse collection of 215 Asiatic cotton accessions were evaluated across three agro-ecological zones of China. Genome-Wide Association Studies (GWAS) was performed to detect association signals related to 17 AFIS fiber quality traits grouped into four categories viz: NEPs, fiber length, maturity, and fineness. Significant correlations were found within as well as among different categories of various traits related to fiber quality. Fiber fineness has shown a strong correlation to all other categories, whereas these categories are shown interrelationships via fiber-fineness. A total of 7,429 SNPs were found in association with 17 investigated traits, of which 177 were selected as lead SNPs. In the vicinity of these lead SNPs, 56 differentially expressed genes in various tissues/development stages were identified as candidate genes. This compendium connecting trait-SNP-genes may allow further prioritization of genes in GWAS loci to enable mechanistic studies. These identified quantitative trait nucleotides (QTNs) may prove helpful in fiber quality improvement in Asian cotton through marker-assisted breeding as well as in reviving eroded genetic factors of G. hirsutum via introgression breeding.

Riddled with holes: Understanding air space formation in plant leaves

Plants use energy from sunlight to transform carbon dioxide from the air into complex organic molecules, ultimately producing much of the food we eat. To make this complex chemistry more efficient, plant leaves are intricately constructed in 3 dimensions: They are flat to maximise light capture and contain extensive internal air spaces to increase gas exchange for photosynthesis. Many years of work has built up an understanding of how leaves form flat blades, but the molecular mechanisms that control air space formation are poorly understood. Here, I review our current understanding of air space formation and outline how recent advances can be harnessed to answer key questions and take the field forward. Increasing our understanding of plant air spaces will not only allow us to understand a fundamental aspect of plant development, but also unlock the potential to engineer the internal structure of crops to make them more efficient at photosynthesis with lower water requirements and more resilient in the face of a changing environment.

Leaves are interwoven with large air spaces to increase the efficiency of photosynthesis; however, how these air spaces form and how different patterns have evolved is almost unknown. This Unsolved Mystery article discusses the existing evidence and poses new avenues of research to answer this question.

Genome-wide identification and expression profiling of DREB genes in Saccharum spontaneum

Background

The dehydration-responsive element-binding proteins (DREBs) are important transcription factors that interact with a DRE/CRT (C-repeat) sequence and involve in response to multiple abiotic stresses in plants. Modern sugarcane are hybrids from the cross between Saccharum spontaneum and Saccharum officinarum, and the high sugar content is considered to the attribution of S. officinaurm, while the stress tolerance is attributed to S. spontaneum. To understand the molecular and evolutionary characterization and gene functions of the DREBs in sugarcane, based on the recent availability of the whole genome information, the present study performed a genome-wide in silico analysis of DREB genes and transcriptome analysis in the polyploidy S. spontaneum.

Results

Twelve DREB1 genes and six DREB2 genes were identified in S. spontaneum genome and all proteins contained a conserved AP2/ERF domain. Eleven SsDREB1 allele genes were assumed to be originated from tandem duplications, and two of them may be derived after the split of S. spontaneum and the proximal diploid species sorghum, suggesting tandem duplication contributed to the expansion of DREB1-type genes in sugarcane. Phylogenetic analysis revealed that one DREB2 gene was lost during the evolution of sugarcane. Expression profiling showed different SsDREB genes with variable expression levels in the different tissues, indicating seven SsDREB genes were likely involved in the development and photosynthesis of S. spontaneum. Furthermore, SsDREB1F, SsDREB1L, SsDREB2D, and SsDREB2F were up-regulated under drought and cold condition, suggesting that these four genes may be involved in both dehydration and cold response in sugarcane.

Conclusions

These findings demonstrated the important role of DREBs not only in the stress response, but also in the development and photosynthesis of S. spontaneum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07799-5.

Major histocompatibility complex genes and locus organization in the Komodo dragon (Varanus komodoensis)

We performed a meta-analysis of the newly assembled Komodo dragon (Varanus komodoensis) genome to characterize the major histocompatibility complex (MHC) of the species. The MHC gene clusters of the Komodo dragon are gene dense, complex, and contain counterparts of many genes of the human MHC. Our analysis identified 20 contigs encompassing ~ 6.9 Mbp of sequence with 223 annotated genes of which many are predicted orthologs to the genes of the human MHC. These MHC contigs range in size from 13.2 kb to 21.5 Mbp, contain an average of one gene per 30 kb, and are thought to occur on at least two chromosomes. Eight contigs, each > 100 kb, could be aligned to the human MHC based on gene content, and these represent gene clusters found in each of the recognized mammalian MHC subregions. The MHC of the Komodo dragon shares organizational features of other non-mammalian taxa. Multiple class Iα and class IIβ genes are indicated, with linkage between classical class I and immunoproteasome genes and between framework class I genes and genes associated with the mammalian class III subregion. These findings are supported in both Komodo genome assemblies and provide new insight into the MHC organization of these unique squamate reptiles.

Morpho-Physiological and Proteomic Response of Bt-Cotton and Non-Bt Cotton to Drought Stress

Drought stress impacts cotton plant growth and productivity across countries. Plants can initiate morphological, cellular, and proteomic changes to adapt to unfavorable conditions. However, our knowledge of how cotton plants respond to drought stress at the proteome level is limited. Herein, we elucidated the molecular coordination underlining the drought tolerance of two inbred cotton varieties, Bacillus thuringiensis-cotton [Bt-cotton + Cry1 Ac gene and Cry 2 Ab gene; NCS BG II BT (BTCS/BTDS)] and Hybrid cotton variety [Non-Bt-cotton; (HCS/HDS)]. Our morphological observations and biochemical experiments showed a different tolerance level between two inbred lines to drought stress. Our proteomic analysis using 2D-DIGE revealed that the changes among them were not obviously in respect to their controls apart from under drought stress, illustrating the differential expression of 509 and 337 proteins in BTDS and HDS compared to their controls. Among these, we identified eight sets of differentially expressed proteins (DEPs) and characterized them using MALDI-TOF/TOF mass spectrometry. Furthermore, the quantitative real-time PCR analysis was carried out with the identified drought-related proteins and confirmed differential expressions. In silico analysis of DEPs using Cytoscape network finds ATPB, NAT9, ERD, LEA, and EMB2001 to be functionally correlative to various drought-responsive genes LEA, AP2/ERF, WRKY, and NAC. These proteins play a vital role in transcriptomic regulation under stress conditions. The higher drought response in Bt cotton (BTCS/BTDS) attributed to the overexpression of photosynthetic proteins enhanced lipid metabolism, increased cellular detoxification and activation chaperones, and reduced synthesis of unwanted proteins. Thus, the Bt variety had enhanced photosynthesis, elevated water retention potential, balanced leaf stomata ultrastructure, and substantially increased antioxidant activity than the Non-Bt cotton. Our results may aid breeders and provide further insights into developing new drought-tolerant and high-yielding cotton hybrid varieties.

Nitric Oxide Alters the Pattern of Auxin Maxima and PIN-FORMED1 During Shoot Development

Hormone patterns tailor cell fate decisions during plant organ formation. Among them, auxins and cytokinins are critical phytohormones during early development. Nitric oxide (NO) modulates root architecture by the control of auxin spatial patterns. However, NO involvement during the coordination of shoot organogenesis remains unclear. Here, we explore the effect of NO during shoot development by using a phenotypic, cellular, and genetic analysis in Arabidopsis thaliana and get new insights into the characterization of NO-mediated leaf-related phenotypes. NO homeostasis mutants are impaired in several shoot architectural parameters, including phyllotactic patterns, inflorescence stem elongation, silique production, leaf number, and margin. Auxin distribution is a key feature for tissue differentiation and need to be controlled at different levels (i.e., synthesis, transport, and degradation mechanisms). The phenotypes resulting from the introduction of the cue1 mutation in the axr1 auxin resistant and pin1 backgrounds exacerbate the relationship between NO and auxins. Using the auxin reporter DR5:GUS, we observed an increase in auxin maxima under NO-deficient mutant backgrounds and NO scavenging, pointing to NO-ASSOCIATED 1 (NOA1) as the main player related to NO production in this process. Furthermore, polar auxin transport is mainly regulated by PIN-FORMED 1 (PIN1), which controls the flow along leaf margin and venations. Analysis of PIN1 protein levels shows that NO controls its accumulation during leaf development, impacting the auxin mediated mechanism of leaf building. With these findings, we also provide evidence for the NO opposite effects to determine root and shoot architecture, in terms of PIN1 accumulation under NO overproduction.

In Arabidopsis thaliana Substrate Recognition and Tissue- as Well as Plastid Type-Specific Expression Define the Roles of Distinct Small Subunits of Isopropylmalate Isomerase

In Arabidopsis thaliana, the heterodimeric isopropylmalate isomerase (IPMI) is composed of a single large (IPMI LSU1) and one of three different small subunits (IPMI SSU1 to 3). The function of IPMI is defined by the small subunits. IPMI SSU1 is required for Leu biosynthesis and has previously also been proposed to be involved in the first cycle of Met chain elongation, the first phase of the synthesis of Met-derived glucosinolates. IPMI SSU2 and IPMI SSU3 participate in the Met chain elongation pathway. Here, we investigate the role of the three IPMI SSUs through the analysis of the role of the substrate recognition region spanning five amino acids on the substrate specificity of IPMI SSU1. Furthermore, we analyze in detail the expression pattern of fluorophore-tagged IPMI SSUs throughout plant development. Our study shows that the substrate recognition region that differs between IPMI SSU1 and the other two IMPI SSUs determines the substrate preference of IPMI. Expression of IPMI SSU1 is spatially separated from the expression of IPMI SSU2 and IPMI SSU3, and IPMI SSU1 is found in small plastids, whereas IMPI SSU2 and SSU3 are found in chloroplasts. Our data show a distinct role for IMPI SSU1 in Leu biosynthesis and for IMPI SSU2 and SSU3 in the Met chain elongation pathway.

Evolution and Functional Divergence of the Fructokinase Gene Family in Populus

New kinase has emerged throughout evolution, but how new kinase evolve while maintaining their functions and acquiring new functions remains unclear. Fructokinase (FRK), the gateway kinase to fructose metabolism, plays essential roles in plant development, and stress tolerance. Here, we explored the evolution of FRK gene family in 20 plant species (from green algae to angiosperms) and their functional roles in Populus. We identified 125 putative FRK genes in the 20 plant species with an average of 6 members per species. Phylogenetic analysis separated these 125 genes into 8 clades including 3 conserved clades and 5 specific clades, the 5 of which only exist in green algae or angiosperms. Evolutionary analysis revealed that FRK genes in ancient land plants have the largest number of functional domains with the longest amino acid sequences, and the length of FRK genes became shorter during the transition to vascular plants. This was accompanied by loss, acquisition, and diversification of functional domains. In Populus, segmental duplication appears to be the main mechanism for the expansion of FRK genes. Specially, most FRK genes duplicated in salicoids are regulated by Populus-specific microRNAs. Furthermore, compared with common FRKs, Populus-specific FRKs have showed higher expression specificity and are associated with fewer growth and wood property traits, which suggests that these FRKs may have undergone functional divergence. Our study explores the specific roles of FRKs in the Populus genome and provides new insights for functional investigation of this gene family.

Intracellular phosphate homeostasis - A short way from metabolism to signaling

Phosphorus in plant cells occurs in inorganic form as both ortho- and pyrophosphate or bound to organic compounds, like e.g., nucleotides, phosphorylated metabolites, phospholipids, phosphorylated proteins, or phytate as P storage in the vacuoles of seeds. Individual compartments of the cell are surrounded by membranes that are selective barriers to avoid uncontrolled solute exchange. A controlled exchange of phosphate or phosphorylated metabolites is accomplished by specific phosphate transporters (PHTs) and the plastidial phosphate translocator family (PTs) of the inner envelope membrane. Plastids, in particular chloroplasts, are the site of various anabolic sequences of enzyme-catalyzed reactions. Apart from their role in metabolism PHTs and PTs are presumed to be also involved in communication between organelles and plant organs. Here we will focus on the integration of phosphate transport and homeostasis in signaling processes. Recent developments in this field will be critically assessed and potential future developments discussed. In particular, the occurrence of various plastid types in one organ (i.e. the leaf) with different functions with respect to metabolism or sensing, as has been documented recently following a tissue-specific proteomics approach (Beltran et al., 2018), will shed new light on functional aspects of phosphate homeostasis.

Imbalance of tyrosine by modulating TyrA arogenate dehydrogenases impacts growth and development of Arabidopsis thaliana

l-Tyrosine is an essential aromatic amino acid required for the synthesis of proteins and a diverse array of plant natural products; however, little is known on how the levels of tyrosine are controlled in planta and linked to overall growth and development. Most plants synthesize tyrosine by TyrA arogenate dehydrogenases, which are strongly feedback-inhibited by tyrosine and encoded by TyrA1 and TyrA2 genes in Arabidopsis thaliana. While TyrA enzymes have been extensively characterized at biochemical levels, their in planta functions remain uncertain. Here we found that TyrA1 suppression reduces seed yield due to impaired anther dehiscence, whereas TyrA2 knockout leads to slow growth with reticulate leaves. The tyra2 mutant phenotypes were exacerbated by TyrA1 suppression and rescued by the expression of TyrA2, TyrA1 or tyrosine feeding. Low-light conditions synchronized the tyra2 and wild-type growth, and ameliorated the tyra2 leaf reticulation. After shifting to normal light, tyra2 transiently decreased tyrosine and subsequently increased aspartate before the appearance of the leaf phenotypes. Overexpression of the deregulated TyrA enzymes led to hyper-accumulation of tyrosine, which was also accompanied by elevated aspartate and reticulate leaves. These results revealed that TyrA1 and TyrA2 have distinct and overlapping functions in flower and leaf development, respectively, and that imbalance of tyrosine, caused by altered TyrA activity and regulation, impacts growth and development of Arabidopsis. The findings provide critical bases for improving the production of tyrosine and its derived natural products, and further elucidating the coordinated metabolic and physiological processes to maintain tyrosine levels in plants.

Genome reorganization of the GmSHMT gene family in soybean showed a lack of functional redundancy in resistance to soybean cyst nematode

In soybeans, eighteen members constitute the serine hydroxymethyltransferase (GmSHMT) gene family, of which the cytosolic-targeted GmSHMT08c member has been reported to mediate resistance to soybean cyst nematode (SCN). This work presents a comprehensive study of the SHMT gene family members, including synteny, phylogeny, subcellular localizations, haplotypes, protein homology modeling, mutational, and expression analyses. Phylogenetic analysis showed that SHMT genes are divided into four classes reflecting their subcellular distribution (cytosol, nucleus, mitochondrion, and chloroplast). Subcellular localization of selected GmSHMT members supports their in-silico predictions and phylogenetic distribution. Expression and functional analyses showed that GmSHMT genes display many overlapping, but some divergent responses during SCN infection. Furthermore, mutational analysis reveals that all isolated EMS mutants that lose their resistance to SCN carry missense and nonsense mutations at the GmSHMT08c, but none of the Gmshmt08c mutants carried mutations in the other GmSHMT genes. Haplotype clustering analysis using the whole genome resequencing data from a collection of 106 diverse soybean germplams (15X) was performed to identify allelic variants and haplotypes within the GmSHMT gene family. Interestingly, only the cytosolic-localized GmSHMT08c presented SNP clusters that were associated with SCN resistance, supporting our mutational analysis. Although eight GmSHMT members respond to the nematode infestation, functional and mutational analysis has shown the absence of functional redundancy in resistance to SCN. Structural analysis and protein homology modeling showed the presence of spontaneous mutations at important residues within the GmSHMT proteins, suggesting the presence of altered enzyme activities based on substrate affinities. Due to the accumulation of mutations during the evolution of the soybean genome, the other GmSHMT members have undergone neofunctionalization and subfunctionalization events.

The Rice Rolled Fine Striped (RFS) CHD3/Mi-2 Chromatin Remodeling Factor Epigenetically Regulates Genes Involved in Oxidative Stress Responses During Leaf Development

In rice (Oryza sativa), moderate leaf rolling increases photosynthetic competence and raises grain yield; therefore, this important agronomic trait has attracted much attention from plant biologists and breeders. However, the relevant molecular mechanism remains unclear. Here, we isolated and characterized Rolled Fine Striped (RFS), a key gene affecting rice leaf rolling, chloroplast development, and reactive oxygen species (ROS) scavenging. The rfs-1 gamma-ray allele and the rfs-2 T-DNA insertion allele of RFS failed to complement each other and their mutants had similar phenotypes, producing extremely incurved leaves due to defective development of vascular cells on the adaxial side. Map-based cloning showed that the rfs-1 mutant harbors a 9-bp deletion in a gene encoding a predicted CHD3/Mi-2 chromatin remodeling factor belonging to the SNF2-ATP-dependent chromatin remodeling family. RFS was expressed in various tissues and accumulated mainly in the vascular cells throughout leaf development. Furthermore, RFS deficiency resulted in a cell death phenotype that was caused by ROS accumulation in developing leaves. We found that expression of five ROS-scavenging genes [encoding catalase C, ascorbate peroxidase 8, a putative copper/zinc superoxide dismutase (SOD), a putative SOD, and peroxiredoxin IIE2] decreased in rfs-2 mutants. Western-blot and chromatin immunoprecipitation (ChIP) assays demonstrated that rfs-2 mutants have reduced H3K4me3 levels in ROS-related genes. Loss-of-function in RFS also led to multiple developmental defects, affecting pollen development, grain filling, and root development. Our results suggest that RFS is required for many aspects of plant development and its function is closely associated with epigenetic regulation of genes that modulate ROS homeostasis.

Non-specific activities of the major herbicide-resistance gene BAR

Bialaphos resistance (BAR) and phosphinothricin acetyltransferase (PAT) genes, which convey resistance to the broad-spectrum herbicide phosphinothricin (also known as glufosinate) via N-acetylation, have been globally used in basic plant research and genetically engineered crops 1-4 . Although early in vitro enzyme assays showed that recombinant BAR and PAT exhibit substrate preference toward phosphinothricin over the 20 proteinogenic amino acids 1 , indirect effects of BAR-containing transgenes in planta, including modified amino acid levels, have been seen but without the identification of their direct causes 5,6 . Combining metabolomics, plant genetics and biochemical approaches, we show that transgenic BAR indeed converts two plant endogenous amino acids, aminoadipate and tryptophan, to their respective N-acetylated products in several plant species. We report the crystal structures of BAR, and further delineate structural basis for its substrate selectivity and catalytic mechanism. Through structure-guided protein engineering, we generated several BAR variants that display significantly reduced non-specific activities compared with its wild-type counterpart in vivo. The transgenic expression of enzymes can result in unintended off-target metabolism arising from enzyme promiscuity. Understanding such phenomena at the mechanistic level can facilitate the design of maximally insulated systems featuring heterologously expressed enzymes.

Surveying the Oligomeric State of Arabidopsis thaliana Chloroplasts

Blue native-PAGE (BN-PAGE) resolves protein complexes in their native state. When combined with immunoblotting, it can be used to identify the presence of high molecular weight complexes at high resolution for any protein, given a suitable antibody. To identify proteins in high molecular weight complexes on a large scale and to bypass the requirement for specific antibodies, we applied a tandem mass spectrometry (MS/MS) approach to BN-PAGE-resolved chloroplasts. Fractionation of the gel into six bands allowed identification and label-free quantification of 1000 chloroplast proteins with native molecular weight separation. Significantly, this approach achieves a depth of identification comparable with traditional shotgun proteomic analyses of chloroplasts, indicating much of the known chloroplast proteome is amenable to MS/MS identification under our fractionation scheme. By limiting the number of fractionation bands to six, we facilitate scaled-up comparative analyses, as we demonstrate with the reticulata chloroplast mutant displaying a reticulated leaf phenotype. Our comparative proteomics approach identified a candidate interacting protein of RETICULATA as well as effects on lipid remodeling proteins, amino acid metabolic enzymes, and plastid division machinery. We additionally highlight selected proteins from each sub-compartment of the chloroplast that provide novel insight on known or hypothesized protein complexes to further illustrate the utility of this approach. Our results demonstrate the high sensitivity and reproducibility of this technique, which is anticipated to be widely adaptable to other sub-cellular compartments.

PP2A Phosphatase as a Regulator of ROS Signaling in Plants

Reactive oxygen species (ROS) carry out vital functions in determining appropriate stress reactions in plants, but the molecular mechanisms underlying the sensing, signaling and response to ROS as signaling molecules are not yet fully understood. Recent studies have underscored the role of Protein Phosphatase 2A (PP2A) in ROS-dependent responses involved in light acclimation and pathogenesis responses in Arabidopsis thaliana. Genetic, proteomic and metabolomic studies have demonstrated that trimeric PP2A phosphatases control metabolic changes and cell death elicited by intracellular and extracellular ROS signals. Associated with this, PP2A subunits contribute to transcriptional and post-translational regulation of pro-oxidant and antioxidant enzymes. This review highlights the emerging role of PP2A phosphatases in the regulatory ROS signaling networks in plants.

GUN1, a Jack-Of-All-Trades in Chloroplast Protein Homeostasis and Signaling

The GENOMES UNCOUPLED 1 (GUN1) gene has been reported to encode a chloroplast-localized pentatricopeptide-repeat protein, which acts to integrate multiple indicators of plastid developmental stage and altered plastid function, as part of chloroplast-to-nucleus retrograde communication. However, the molecular mechanisms underlying signal integration by GUN1 have remained elusive, up until the recent identification of a set of GUN1-interacting proteins, by co-immunoprecipitation and mass-spectrometric analyses, as well as protein-protein interaction assays. Here, we review the molecular functions of the different GUN1 partners and propose a major role for GUN1 as coordinator of chloroplast translation, protein import, and protein degradation. This regulatory role is implemented through proteins that, in most cases, are part of multimeric protein complexes and whose precise functions vary depending on their association states. Within this framework, GUN1 may act as a platform to promote specific functions by bringing the interacting enzymes into close proximity with their substrates, or may inhibit processes by sequestering particular pools of specific interactors. Furthermore, the interactions of GUN1 with enzymes of the tetrapyrrole biosynthesis (TPB) pathway support the involvement of tetrapyrroles as signaling molecules in retrograde communication.

Protein phosphatase 2A (PP2A) regulatory subunit B'γ interacts with cytoplasmic ACONITASE 3 and modulates the abundance of AOX1A and AOX1D in Arabidopsis thaliana

Organellar reactive oxygen species (ROS) signalling is a key mechanism that promotes the onset of defensive measures in stress-exposed plants. The underlying molecular mechanisms and feedback regulation loops, however, still remain poorly understood. Our previous work has shown that a specific regulatory B'γ subunit of protein phosphatase 2A (PP2A) is required to control organellar ROS signalling and associated metabolic adjustments in Arabidopsis thaliana. Here, we addressed the mechanisms through which PP2A-B'γ impacts on organellar metabolic crosstalk and ROS homeostasis in leaves. Genetic, biochemical and pharmacological approaches, together with a combination of data-dependent acquisition (DDA) and selected reaction monitoring (SRM) MS techniques, were utilized to assess PP2A-B'γ-dependent adjustments in Arabidopsis thaliana. We show that PP2A-B'γ physically interacts with the cytoplasmic form of aconitase, a central metabolic enzyme functionally connected with mitochondrial respiration, oxidative stress responses and regulation of cell death in plants. Furthermore, PP2A-B'γ impacts ROS homeostasis by controlling the abundance of specific alternative oxidase isoforms, AOX1A and AOX1D, in leaf mitochondria. We conclude that PP2A-B'γ-dependent regulatory actions modulate the functional status of metabolic enzymes that essentially contribute to intracellular ROS signalling and metabolic homeostasis in plants.

Rapid identification of angulata leaf mutations using next-generation sequencing

Map-based (positional) cloning has traditionally been the preferred strategy for identifying the causal genes underlying the phenotypes of mutants isolated in forward genetic screens. Massively parallel sequencing technologies are enabling the rapid cloning of genes identified in such screens. We have used a combination of linkage mapping and whole-genome re-sequencing to identify the causal mutations in four loss-of-function angulata (anu) mutants. These mutants were isolated in a screen for mutants with defects in leaf shape and leaf pigmentation. Our results show that the anu1-1, anu4-1, anu9-1 and anu12-1 mutants carry new alleles of the previously characterized SECA2, TRANSLOCON AT THE OUTER MEMBRANE OF CHLOROPLASTS 33 (TOC33), NON-INTRINSIC ABC PROTEIN 14 (NAP14) and CLP PROTEASE PROTEOLYTIC SUBUNIT 1 (CLPR1) genes. Re-sequencing the genomes of fine mapped mutants is a feasible approach that has allowed us to identify a moderate number of candidate mutations, including the one that causes the mutant phenotype, in a nonstandard genetic background. Our results indicate that anu mutations specifically affect plastid-localized proteins involved in diverse processes, such as the movement of peptides through chloroplast membranes (ANU1 and ANU4), metal homeostasis (ANU9) and protein degradation (ANU12).

DUF581 is plant specific FCS-like zinc finger involved in protein-protein interaction

Zinc fingers are a ubiquitous class of protein domain with considerable variation in structure and function. Zf-FCS is a highly diverged group of C₂-C₂ zinc finger which is present in animals, prokaryotes and viruses, but not in plants. In this study we identified that a plant specific domain of unknown function, DUF581 is a zf-FCS type zinc finger. Based on HMM-HMM comparison and signature motif similarity we named this domain as FCS-Like Zinc finger (FLZ) domain. A genome wide survey identified that FLZ domain containing genes are bryophytic in origin and this gene family is expanded in spermatophytes. Expression analysis of selected FLZ gene family members of A. thaliana identified an overlapping expression pattern suggesting a possible redundancy in their function. Unlike the zf-FCS domain, the FLZ domain found to be highly conserved in sequence and structure. Using a combination of bioinformatic and protein-protein interaction tools, we identified that FLZ domain is involved in protein-protein interaction.

Arabidopsis ANGULATA10 is required for thylakoid biogenesis and mesophyll development

The chloroplasts of land plants contain internal membrane systems, the thylakoids, which are arranged in stacks called grana. Because grana have not been found in Cyanobacteria, the evolutionary origin of genes controlling the structural and functional diversification of thylakoidal membranes in land plants remains unclear. The angulata10-1 (anu10-1) mutant, which exhibits pale-green rosettes, reduced growth, and deficient leaf lateral expansion, resulting in the presence of prominent marginal teeth, was isolated. Palisade cells in anu10-1 are larger and less packed than in the wild type, giving rise to large intercellular spaces. The ANU10 gene encodes a protein of unknown function that localizes to both chloroplasts and amyloplasts. In chloroplasts, ANU10 associates with thylakoidal membranes. Mutant anu10-1 chloroplasts accumulate H2O2, and have reduced levels of chlorophyll and carotenoids. Moreover, these chloroplasts are small and abnormally shaped, thylakoidal membranes are less abundant, and their grana are absent due to impaired thylakoid stacking in the anu10-1 mutant. Because the trimeric light-harvesting complex II (LHCII) has been reported to be required for thylakoid stacking, its levels were determined in anu10-1 thylakoids and they were found to be reduced. Together, the data point to a requirement for ANU10 for chloroplast and mesophyll development.

Towards Understanding Extracellular ROS Sensory and Signaling Systems in Plants

Reactive Oxygen Species (ROS) are ubiquitous metabolites in all aerobic organisms. Traditionally ROS have been considered as harmful, accidental byproducts of cellular functions involving electron transport chains or electron transfer. However, it is now recognized that controlled production of ROS has significant signaling functions, for example, in pathogen defense, in the regulation of stomatal closure, or in cell-to-cell signaling. ROS formation in subcellular compartments is critical to act as “alarm” signal in the response to stress, and the concept of ROS as primarily signaling substances has emerged. The involvement of ROS in several developmental and inducible processes implies that there must be coordinated function of signaling network(s) that govern ROS responses and subsequent processes. The air pollutant ozone can be used as a useful tool to elucidate the function of apoplastic ROS: O3 degrades in cell wall into various ROS which are interpreted as ROS with signaling function inducing downstream responses. We have used ozone as a tool in mutant screens and transcript profiling-reverse genetics to identify genes involved in processes related to the signaling function of ROS. We review here our recent findings in the elucidation of apoplastic ROS sensing, signaling, and interaction with various symplastic components.

The protein phosphatase subunit PP2A-B'γ is required to suppress day length-dependent pathogenesis responses triggered by intracellular oxidative stress

Oxidative stress responses are influenced by growth day length, but little is known about how this occurs. A combined reverse genetics, metabolomics and proteomics approach was used to address this question in Arabidopsis thaliana. A catalase-deficient mutant (cat2), in which intracellular oxidative stress drives pathogenesis-related responses in a day length-dependent manner, was crossed with a knockdown mutant for a specific type 2A protein phosphatase subunit (pp2a-b'γ). In long days (LD), the pp2a-b'γ mutation reinforced cat2-triggered pathogenesis responses. In short days (SD), conditions in which pathogenesis-related responses were not activated in cat2, the additional presence of the pp2a-b'γ mutation allowed lesion formation, PATHOGENESIS-RELATED GENE1 (PR1) induction, salicylic acid (SA) and phytoalexin accumulation and the establishment of metabolite profiles that were otherwise observed in cat2 only in LD. Lesion formation in cat2 pp2a-b'γ in SD was genetically dependent on SA synthesis, and was associated with decreased PHYTOCHROME A transcripts. Phosphoproteomic analyses revealed that several potential protein targets accumulated in the double mutant, including recognized players in pathogenesis and key enzymes of primary metabolism. We conclude that the cat2 and pp2a-b'γ mutations interact synergistically, and that PP2A-B'γ is an important player in controlling day length-dependent responses to intracellular oxidative stress, possibly through phytochrome-linked pathways.

Signalling crosstalk in light stress and immune reactions in plants

The evolutionary history of plants is tightly connected with the evolution of microbial pathogens and herbivores, which use photosynthetic end products as a source of life. In these interactions, plants, as the stationary party, have evolved sophisticated mechanisms to sense, signal and respond to the presence of external stress agents. Chloroplasts are metabolically versatile organelles that carry out fundamental functions in determining appropriate immune reactions in plants. Besides photosynthesis, chloroplasts host key steps in the biosynthesis of amino acids, stress hormones and secondary metabolites, which have a great impact on resistance against pathogens and insect herbivores. Changes in chloroplast redox signalling pathways and reactive oxygen species metabolism also mediate local and systemic signals, which modulate plant resistance to light stress and disease. Moreover, interplay among chloroplastic signalling networks and plasma membrane receptor kinases is emerging as a key mechanism that modulates stress responses in plants. This review highlights the central role of chloroplasts in the signalling crosstalk that essentially determines the outcome of plant-pathogen interactions in plants.

Plastid signals and the bundle sheath: mesophyll development in reticulate mutants

The development of a plant leaf is a meticulously orchestrated sequence of events producing a complex organ comprising diverse cell types. The reticulate class of leaf variegation mutants displays contrasting pigmentation between veins and interveinal regions due to specific aberrations in the development of mesophyll cells. Thus, the reticulate mutants offer a potent tool to investigate cell-type-specific developmental processes. The discovery that most mutants are affected in plastid-localized, metabolic pathways that are strongly expressed in vasculature-associated tissues implicates a crucial role for the bundle sheath and their chloroplasts in proper development of the mesophyll cells. Here, we review the reticulate mutants and their phenotypic characteristics, with a focus on those in Arabidopsis thaliana. Two alternative models have been put forward to explain the relationship between plastid metabolism and mesophyll cell development, which we call here the supply and the signaling hypotheses. We critically assess these proposed models and discuss their implications for leaf development and bundle sheath function in C3 species. The characterization of the reticulate mutants supports the significance of plastid retrograde signaling in cell development and highlights the significance of the bundle sheath in C3 photosynthesis.

Mutation of an Arabidopsis NatB N-alpha-terminal acetylation complex component causes pleiotropic developmental defects

N-α-terminal acetylation is one of the most common, but least understood modifications of eukaryotic proteins. Although a high degree of conservation exists between the N-α-terminal acetylomes of plants and animals, very little information is available on this modification in plants. In yeast and humans, N-α-acetyltransferase complexes include a single catalytic subunit and one or two auxiliary subunits. Here, we report the positional cloning of TRANSCURVATA2 (TCU2), which encodes the auxiliary subunit of the NatB N-α-acetyltransferase complex in Arabidopsis. The phenotypes of loss-of-function tcu2 alleles indicate that NatB complex activity is required for flowering time regulation and for leaf, inflorescence, flower, fruit and embryonic development. In double mutants, tcu2 alleles synergistically interact with alleles of ARGONAUTE10, which encodes a component of the microRNA machinery. In summary, NatB-mediated N-α-terminal acetylation of proteins is pleiotropically required for Arabidopsis development and seems to be functionally related to the microRNA pathway.