ePlant: Visualizing and Exploring Multiple Levels of Data for Hypothesis Generation in Plant Biology

Advances and challenges of immunocheckpoint inhibitors in the treatment of primary liver cancer

Primary liver cancer (PLC) is one of the most common malignant tumors, which clinically characterized by occult onset, rapid development, easy recurrence and poor prognosis. With the rapid development of tumor immunotherapy research, tumor immunotherapy has also achieved remarkable clinical efficacy, and jointly promoted the overall improvement of tumor immunology from mechanism research to clinical transformation, from single discipline to multi-disciplinary integration. Immunotherapy has obvious advantages in treatment-related toxicity and efficacy compared with traditional therapy. In hepatocellular carcinoma (HCC), immunotherapy alone or in combination with other therapies may help to control tumor progression, and there are many immune checkpoint inhibitors (ICIs) widely used in clinical or ongoing clinical trials. However, tumor immunology research is still facing many challenges. How to effectively evaluate the efficacy, whether there are related biomarkers, the generation of immune tolerance and the lack of clinical trials to objectively evaluate the efficacy are still urgent problems to be solved, but it also brings new research opportunities for basic and clinical immunology researchers. The study of treatment of ICIs of PLC has become a hot spot in clinical research field. This paper summarizes and prospects the research progress and challenges of ICIs for PLC.

Core circadian clock and light signaling genes brought into genetic linkage across the green lineage

The circadian clock is conserved at both the level of transcriptional networks as well as core genes in plants, ensuring that biological processes are phased to the correct time of day. In the model plant Arabidopsis (Arabidopsis thaliana), the core circadian SHAQKYF-type-MYB (sMYB) genes CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) and REVEILLE (RVE4) show genetic linkage with PSEUDO-RESPONSE REGULATOR 9 (PRR9) and PRR7, respectively. Leveraging chromosome-resolved plant genomes and syntenic ortholog analysis enabled tracing this genetic linkage back to Amborella trichopoda, a sister lineage to the angiosperm, and identifying an additional evolutionarily conserved genetic linkage in light signaling genes. The LHY/CCA1-PRR5/9, RVE4/8-PRR3/7, and PIF3-PHYA genetic linkages emerged in the bryophyte lineage and progressively moved within several genes of each other across an array of angiosperm families representing distinct whole-genome duplication and fractionation events. Soybean (Glycine max) maintained all but two genetic linkages, and expression analysis revealed the PIF3-PHYA linkage overlapping with the E4 maturity group locus was the only pair to robustly cycle with an evening phase, in contrast to the sMYB-PRR morning and midday phase. While most monocots maintain the genetic linkages, they have been lost in the economically important grasses (Poaceae), such as maize (Zea mays), where the genes have been fractionated to separate chromosomes and presence/absence variation results in the segregation of PRR7 paralogs across heterotic groups. The environmental robustness model is put forward, suggesting that evolutionarily conserved genetic linkages ensure superior microhabitat pollinator synchrony, while wide-hybrids or unlinking the genes, as seen in the grasses, result in heterosis, adaptation, and colonization of new ecological niches.

Diversification of heat shock transcription factors expanded thermal stress responses during early plant evolution

The copy numbers of many plant transcription factor (TF) genes substantially increased during terrestrialization. This allowed TFs to acquire new specificities and thus create gene regulatory networks (GRNs) with new biological functions to help plants adapt to terrestrial environments. Through characterizing heat shock factor (HSF) genes MpHSFA1 and MpHSFB1 in the liverwort Marchantia polymorpha, we explored how heat-responsive GRNs widened their functions in M. polymorpha and Arabidopsis thaliana. An interspecies comparison of heat-induced transcriptomes and the evolutionary rates of HSFs demonstrated the emergence and subsequent rapid evolution of HSFB prior to terrestrialization. Transcriptome and metabolome analyses of M. polymorpha HSF-null mutants revealed that MpHSFA1 controls canonical heat responses such as thermotolerance and metabolic changes. MpHSFB1 also plays essential roles in heat responses, as well as regulating developmental processes including meristem branching and antheridiophore formation. Analysis of cis-regulatory elements revealed development- and stress-related TFs that function directly or indirectly downstream of HSFB. Male gametophytes of M. polymorpha showed higher levels of thermotolerance than female gametophytes, which could be explained by different expression levels of MpHSFA1U and MpHSFA1V on sex chromosome. We propose that the diversification of HSFs is linked to the expansion of HS responses, which enabled coordinated multicellular reactions in land plants.

Identification and characterization of putative targets of VEGETATIVE1/FULc, a key regulator of development of the compound inflorescence in pea and related legumes

Inflorescence architecture contributes to essential plant traits. It determines plant shape, contributing to morphological diversity, and also determines the position and number of flowers and fruits produced by the plant, thus influencing seed yield. Most legumes have compound inflorescences, where flowers are produced in secondary inflorescences (I2), formed at the flanks of the main primary inflorescence (I1), in contrast to simple inflorescences of plants like Arabidopsis, in which flowers are directly formed on the I1. The pea VEGETATIVE1/FULc (VEG1) gene, and its homologs in other legumes, specify the formation of the I2 meristem, a function apparently restricted to legumes. To understand the control of I2 development, it is important to identify the genes working downstream of VEG1. In this study, we adopted a novel strategy to identify genes expressed in the I2 meristem, as potential regulatory targets of VEG1. To identify pea I2-meristem genes, we compared the transcriptomes of inflorescence apices from wild-type and mutants affected in I2 development, such as proliferating inflorescence meristems (pim, with more I2 meristems), and veg1 and vegetative2 (both without I2 meristems). Analysis of the differentially expressed genes using Arabidopsis genome databases combined with RT-qPCR expression analysis in pea allowed the selection of genes expressed in the pea inflorescence apex. In situ hybridization of four of these genes showed that all four genes are expressed in the I2 meristem, proving our approach to identify I2-meristem genes was successful. Finally, analysis by VIGS (virus-induced gene silencing) in pea identified one gene, PsDAO1, whose silencing leads to small plants, and another gene, PsHUP54, whose silencing leads to plants with very large stubs, meaning that this gene controls the activity of the I2 meristem. PsHUP54-VIGS plants are also large and, more importantly, produce large pods with almost double the seeds as the control. Our study shows a new useful strategy to isolate I2-meristem genes and identifies a novel gene, PsHUP54, which seems to be a promising tool to improve yield in pea and in other legumes.

Uncovering natural variation in root system architecture and growth dynamics using a robotics-assisted phenomics platform

The plant kingdom contains a stunning array of complex morphologies easily observed above-ground, but more challenging to visualize below-ground. Understanding the magnitude of diversity in root distribution within the soil, termed root system architecture (RSA), is fundamental in determining how this trait contributes to species adaptation in local environments. Roots are the interface between the soil environment and the shoot system and therefore play a key role in anchorage, resource uptake, and stress resilience. Previously, we presented the GLO-Roots (Growth and Luminescence Observatory for Roots) system to study the RSA of soil-grown Arabidopsis thaliana plants from germination to maturity (Rellán-Álvarez et al., 2015). In this study, we present the automation of GLO-Roots using robotics and the development of image analysis pipelines in order to examine the temporal dynamic regulation of RSA and the broader natural variation of RSA in Arabidopsis, over time. These datasets describe the developmental dynamics of two independent panels of accessions and reveal highly complex and polygenic RSA traits that show significant correlation with climate variables of the accessions' respective origins.

Ethylene Insensitive 3-Like 2 is a Brassicaceae-specific transcriptional regulator involved in fine-tuning ethylene responses in Arabidopsis thaliana

Ethylene signaling directs a pleiotropy of developmental processes in plants. In Arabidopsis, ethylene signaling converges at the master transcription factor Ethylene Insensitive 3 (EIN3), which has five homologs, EIN3-like 1-5 (EIL1-EIL5). EIL1 is most fully characterized and operates similarly to EIN3, while EIL3-5 are not involved in ethylene signaling. EIL2 remains less investigated. Our phylogenetic analysis revealed that EIL2 homologs have only been retrieved in the Brassicaceae family, suggesting that EIL2 diverged to have specific functions in the mustard family. By characterizing eil2 mutants, we found that EIL2 is involved in regulating ethylene-specific developmental processes in Arabidopsis thaliana, albeit in a more subtle way compared with EIN3/EIL1. EIL2 steers ethylene-triggered hypocotyl elongation in light-grown seedlings and is involved in lateral root formation. Furthermore, EIL2 takes part in regulating flowering time as eil2 mutants flower on average 1 d earlier and have fewer leaves. A pEIL2:EIL2:GFP translational reporter line revealed that EIL2 protein abundance is restricted to the stele of young developing roots. EIL2 expression, and not EIL2 protein stability, is regulated by ethylene in an EIN3/EIL1-dependent way. Despite EIL2 taking part in several developmental processes, the precise upstream and downstream regulation of this ethylene- and Brassicaceae-specific transcription factor remains to be elucidated.

NetREx: Network-based Rice Expression Analysis Server for abiotic stress conditions

Recent focus on transcriptomic studies in food crops like rice, wheat and maize provide new opportunities to address issues related to agriculture and climate change. Re-analysis of such data available in public domain supplemented with annotations across molecular hierarchy can be of immense help to the plant research community, particularly co-expression networks representing transcriptionally coordinated genes that are often part of the same biological process. With this objective, we have developed NetREx, a Network-based Rice Expression Analysis Server, that hosts ranked co-expression networks of Oryza sativa using publicly available messenger RNA sequencing data across uniform experimental conditions. It provides a range of interactable data viewers and modules for analysing user-queried genes across different stress conditions (drought, flood, cold and osmosis) and hormonal treatments (abscisic and jasmonic acid) and tissues (root and shoot). Subnetworks of user-defined genes can be queried in pre-constructed tissue-specific networks, allowing users to view the fold change, module memberships, gene annotations and analysis of their neighbourhood genes and associated pathways. The web server also allows querying of orthologous genes from Arabidopsis, wheat, maize, barley and sorghum. Here, we demonstrate that NetREx can be used to identify novel candidate genes and tissue-specific interactions under stress conditions and can aid in the analysis and understanding of complex phenotypes linked to stress response in rice. Database URL: https://bioinf.iiit.ac.in/netrex/index.html.

Identification of genetic loci conferring seed coat color based on a high-density map in soybean

Seed coat color is a typical evolutionary trait. Identification of the genetic loci that control seed coat color during the domestication of wild soybean could clarify the genetic variations between cultivated and wild soybean. We used 276 F10 recombinant inbred lines (RILs) from the cross between a cultivated soybean (JY47) and a wild soybean (ZYD00321) as the materials to identify the quantitative trait loci (QTLs) for seed coat color. We constructed a high-density genetic map using re-sequencing technology. The average distance between adjacent markers was 0.31 cM on this map, comprising 9,083 bin markers. We identified two stable QTLs (qSC08 and qSC11) for seed coat color using this map, which, respectively, explained 21.933 and 26.934% of the phenotypic variation. Two candidate genes (CHS3C and CHS4A) in qSC08 were identified according to the parental re-sequencing data and gene function annotations. Five genes (LOC100786658, LOC100801691, LOC100806824, LOC100795475, and LOC100787559) were predicted in the novel QTL qSC11, which, according to gene function annotations, might control seed coat color. This result could facilitate the identification of beneficial genes from wild soybean and provide useful information to clarify the genetic variations for seed coat color in cultivated and wild soybean.

Adenylates regulate Arabidopsis plastidial thioredoxin activities through the binding of a CBS domain protein

Cystathionine-β-synthase (CBS) domains are found in proteins of all living organisms and have been proposed to play a role as energy sensors regulating protein activities through their adenosyl ligand binding capacity. In plants, members of the CBSX protein family carry a stand-alone pair of CBS domains. In Arabidopsis (Arabidopsis thaliana), CBSX1 and CBSX2 are targeted to plastids where they have been proposed to regulate thioredoxins (TRXs). TRXs are ubiquitous cysteine thiol oxido-reductases involved in the redox-based regulation of numerous enzymatic activities as well as in the regeneration of thiol-dependent peroxidases. In Arabidopsis, 10 TRX isoforms have been identified in plastids and divided into five sub-types. Here, we show that CBSX2 specifically inhibits the activities of m-type TRXs toward two chloroplast TRX-related targets. By testing activation of NADP-malate dehydrogenase and reduction of 2-Cys peroxiredoxin, we found that TRXm1/2 inhibition by CBSX2 was alleviated in the presence of AMP or ATP. We also determined, by pull-down assays, a direct interaction of CBSX2 with reduced TRXm1 and m2 that was abolished in the presence of adenosyl ligands. In addition, we report that, compared with wild-type plants, the Arabidopsis T-DNA double mutant cbsx1 cbsx2 exhibits growth and chlorophyll accumulation defects in cold conditions, suggesting a function of plastidial CBSX proteins in plant stress adaptation. Together, our results show an energy-sensing regulation of plastid TRX m activities by CBSX, possibly allowing a feedback regulation of ATP homeostasis via activation of cyclic electron flow in the chloroplast, to maintain a high energy level for optimal growth.

De-etiolation-induced protein 1 (DEIP1) mediates assembly of the cytochrome b6f complex in Arabidopsis

The conversion of light energy to chemical energy by photosynthesis requires the concerted action of large protein complexes in the thylakoid membrane. Recent work has provided fundamental insights into the three-dimensional structure of these complexes, but how they are assembled from hundreds of parts remains poorly understood. Particularly little is known about the biogenesis of the cytochrome b₆f complex (Cytb₆f), the redox-coupling complex that interconnects the two photosystems. Here we report the identification of a factor that guides the assembly of Cytb₆f in thylakoids of chloroplasts. The protein, DE-ETIOLATION-INDUCED PROTEIN 1 (DEIP1), resides in the thylakoid membrane and is essential for photoautotrophic growth. Knock-out mutants show a specific loss of Cytb₆f, and are defective in complex assembly. We demonstrate that DEIP1 interacts with the two cytochrome subunits of the complex, PetA and PetB, and mediates the assembly of intermediates in Cytb₆f biogenesis. The identification of DEIP1 provides an entry point into the study of the assembly pathway of a crucial complex in photosynthetic electron transfer.

The Cytb6f complex is a multi-subunit enzyme that couples the two photosystems during the light reactions of photosynthesis. Here the authors show that the thylakoid-localized DEIP1 protein interacts with the PetA and PetB subunits, and is essential for Cytb6f complex assembly in Arabidopsis.

Exploiting plant transcriptomic databases: Resources, tools, and approaches

There are now more than 300 000 RNA sequencing samples available, stemming from thousands of experiments capturing gene expression in organs, tissues, developmental stages, and experimental treatments for hundreds of plant species. The expression data have great value, as they can be re-analyzed by others to ask and answer questions that go beyond the aims of the study that generated the data. Because gene expression provides essential clues to where and when a gene is active, the data provide powerful tools for predicting gene function, and comparative analyses allow us to study plant evolution from a new perspective. This review describes how we can gain new knowledge from gene expression profiles, expression specificities, co-expression networks, differential gene expression, and experiment correlation. We also introduce and demonstrate databases that provide user-friendly access to these tools.

Prediction and expression analysis of deleterious nonsynonymous SNPs of Arabidopsis ACD11 gene by combining computational algorithms and molecular docking approach

Accelerated cell death 11 (ACD11) is an autoimmune gene that suppresses pathogen infection in plants by preventing plant cells from becoming infected by any pathogen. This gene is widely known for growth inhibition, premature leaf chlorosis, and defense-related programmed cell death (PCD) in seedlings before flowering in Arabidopsis plant. Specific amino acid changes in the ACD11 protein’s highly conserved domains are linked to autoimmune symptoms including constitutive defensive responses and necrosis without pathogen awareness. The molecular aspect of the aberrant activity of the ACD11 protein is difficult to ascertain. The purpose of our study was to find the most deleterious mutation position in the ACD11 protein and correlate them with their abnormal expression pattern. Using several computational methods, we discovered PCD vulnerable single nucleotide polymorphisms (SNPs) in ACD11. We analysed the RNA-Seq data, identified the detrimental nonsynonymous SNPs (nsSNP), built genetically mutated protein structures and used molecular docking to assess the impact of mutation. Our results demonstrated that the A15T and A39D mutations in the GLTP domain were likely to be extremely detrimental mutations that inhibit the expression of the ACD11 protein domain by destabilizing its composition, as well as disrupt its catalytic effectiveness. When compared to the A15T mutant, the A39D mutant was more likely to destabilize the protein structure. In conclusion, these mutants can aid in the better understanding of the vast pool of PCD susceptibilities connected to ACD11 gene GLTP domain activation.

Non synonymous single nucleotide polymorphism (nsSNP) is a process in which amino acid sequence of a protein is altered as a result of single nucleotide alteration in the coding region (mRNA) of any living organism. Therefore, the entire protein structure, interactions and stability are altered, which may have a negative impact on living organisms. Hence, to completely comprehend this biological process, we must first solve the unresolved mutational protein structure and mutated protein interactions. The major goal of our research is to identify the most harmful mutation in our target protein structure and how it interacts within cells. However, it was discovered that only a few alterations in residues had the largest negative impact on the protein’s internal structure and also on the protein-ligand interactions. We show that based on the amino acid sequence of a protein computationally, it is feasible to discover mutational positions in the sequence, generate mutation protein structure and interactions with related ligands. Our findings show that the essential mechanisms underlying protein mutations generated by this process are identical. The capacity to correctly detect mutations from sequence allows the annotation and study of protein-ligand interactions throughout a whole organism, which might aid function prediction and gene expression.

REGENOMICS: A web-based application for plant REGENeration-associated transcriptOMICS analyses

In plants, differentiated somatic cells exhibit an exceptional ability to regenerate new tissues, organs, or whole plants. Recent studies have unveiled core genetic components and pathways underlying cellular reprogramming and de novo tissue regeneration in plants. Although high-throughput analyses have led to key discoveries in plant regeneration, a comprehensive organization of large-scale data is needed to further enhance our understanding of plant regeneration. Here, we collected all currently available transcriptome datasets related to wounding responses, callus formation, de novo organogenesis, somatic embryogenesis, and protoplast regeneration to construct REGENOMICS, a web-based application for plant REGENeration-associated transcriptOMICS analyses. REGENOMICS supports single- and multi-query analyses of plant regeneration-related gene-expression dynamics, co-expression networks, gene-regulatory networks, and single-cell expression profiles. Furthermore, it enables user-friendly transcriptome-level analysis of REGENOMICS-deposited and user-submitted RNA-seq datasets. Overall, we demonstrate that REGENOMICS can serve as a key hub of plant regeneration transcriptome analysis and greatly enhance our understanding on gene-expression networks, new molecular interactions, and the crosstalk between genetic pathways underlying each mode of plant regeneration. The REGENOMICS web-based application is available at http://plantregeneration.snu.ac.kr.

Hydroxyproline-O-Galactosyltransferases Synthesizing Type II Arabinogalactans Are Essential for Male Gametophytic Development in Arabidopsis

In flowering plants, male reproductive function is determined by successful development and performance of stamens, pollen grains, and pollen tubes. Despite the crucial role of highly glycosylated arabinogalactan-proteins (AGPs) in male gamete formation, pollen grain, and pollen tube cell walls, the underlying mechanisms defining these functions of AGPs have remained elusive. Eight partially redundant Hyp-galactosyltransferases (named GALT2-GALT9) genes/enzymes are known to initiate Hyp-O-galactosylation for Hyp-arabinogalactan (AG) production in Arabidopsis thaliana. To assess the contributions of these Hyp-AGs to male reproductive function, we used a galt2galt5galt7galt8galt9 quintuple Hyp-GALT mutant for this study. Both anther size and pollen viability were compromised in the quintuple mutants. Defects in male gametogenesis were observed in later stages of maturing microspores after meiosis, accompanied by membrane blebbing and numerous lytic vacuoles. Cytological and ultramicroscopic observations revealed that pollen exine reticulate architecture and intine layer development were affected such that non-viable collapsed mature pollen grains were produced, which were devoid of cell content and nuclei, with virtually no intine. AGP immunolabeling demonstrated alterations in cell wall architecture of the anther, pollen grains, and pollen tube. Specifically, the LM2 monoclonal antibody (which recognized β-GlcA epitopes on AGPs) showed a weak signal for the endothecium, microspores, and pollen tube apex. Pollen tube tips also displayed excessive callose deposition. Interestingly, expression patterns of pollen-specific AGPs, namely AGP6, AGP11, AGP23, and AGP40, were determined to be higher in the quintuple mutants. Taken together, our data illustrate the importance of type-II AGs in male reproductive function for successful fertilization.

Evolutionary and Gene Expression Analyses Reveal New Insights into the Role of LSU Gene-Family in Plant Responses to Sulfate-Deficiency

LSU proteins belong to a plant-specific gene family initially characterized by their strong induction in response to sulfate (S) deficiency. In the last few years, LSUs have arisen as relevant hubs in protein–protein interaction networks, in which they play relevant roles in the response to abiotic and biotic stresses. Most of our knowledge on LSU genomic organization, expression and function comes from studies in Arabidopsis and tobacco, while little is known about the LSU gene repertoire and evolution of this family in land plants. In this work, a total of 270 LSU family members were identified using 134 land plant species with whole-genome sequences available. Phylogenetic analysis revealed that LSU genes belong to a Spermatophyta-specific gene family, and their homologs are distributed in three major groups, two for dicotyledons and one group for monocotyledons. Protein sequence analyses showed four new motifs that further support the subgroup classification by phylogenetic analyses. Moreover, we analyzed the expression of LSU genes in one representative species of each phylogenetic group (wheat, tomato and Arabidopsis) and found a conserved response to S deficiency, suggesting that these genes might play a key role in S stress responses. In summary, our results indicate that LSU genes belong to the Spermatophyta-specific gene family and their response to S deficiency is conserved in angiosperms.

The evolutionary history of small RNAs in Solanaceae

The Solanaceae or "nightshade" family is an economically important group with remarkable diversity. To gain a better understanding of how the unique biology of the Solanaceae relates to the family's small RNA (sRNA) genomic landscape, we downloaded over 255 publicly available sRNA data sets that comprise over 2.6 billion reads of sequence data. We applied a suite of computational tools to predict and annotate two major sRNA classes: (1) microRNAs (miRNAs), typically 20- to 22-nucleotide (nt) RNAs generated from a hairpin precursor and functioning in gene silencing and (2) short interfering RNAs (siRNAs), including 24-nt heterochromatic siRNAs typically functioning to repress repetitive regions of the genome via RNA-directed DNA methylation, as well as secondary phased siRNAs and trans-acting siRNAs generated via miRNA-directed cleavage of a polymerase II-derived RNA precursor. Our analyses described thousands of sRNA loci, including poorly understood clusters of 22-nt siRNAs that accumulate during viral infection. The birth, death, expansion, and contraction of these sRNA loci are dynamic evolutionary processes that characterize the Solanaceae family. These analyses indicate that individuals within the same genus share similar sRNA landscapes, whereas comparisons between distinct genera within the Solanaceae reveal relatively few commonalities.

The RNA binding protein OsLa influences grain and anther development in rice

La proteins are found widely in eukaryotes and play a variety of vital roles. AtLa1 has been identified as an La protein that is necessary for embryogenesis in Arabidopsis; however, the existence and biological functions of La proteins in rice (Oryza sativa L.) remain unclear. In this study, we identified and characterized two La proteins in rice that are homologous to AtLa1 and named them OsLa1 and OsLa2. Both the OsLa1 and OsLa2 genes encode RNA-binding proteins with an La domain and two RNA-binding domains. Mutant OsLa1 reduced grain length and pollen fertility, whereas OsLa1 overexpression caused the opposite phenotypes. Further experiments indicated that OsLa1 modulates grain size by influencing cell expansion. Interestingly, mutant OsLa2 resulted in thin grains with decreased weight and a low seed-setting rate. We also found that OsLa1 interacted with OsLa2 and that both OsLa1 and OsLa2 interacted with OseIF6.1, a eukaryotic translation initiation factor involved in ribosome biogenesis. In addition, OsLa1 was able to bind to OseIF6.1 mRNA to modulate its expression. Complete OseIF6.1 knockout caused lethality and OseIF6.1/oseif6.1 heterozygous plants displayed low fertility and low seed setting. Together, our results enrich our knowledge of the role of La proteins in rice growth and development, as well as the relationship between La and eIF6 in rice.

Arabidopsis CAMTA3/SR1 is involved in drought stress tolerance and ABA signaling

Calcium/calmodulin signals are important for various cellular and physiological activities in plants. Calmodulin binding transcription activators also named Signal Responsive (SR) proteins belong to an important calcium/calmodulin-dependent transcription factor family that plays critical roles in stress responses. However, the role of SRs in abscisic acid (ABA) regulated plant responses to drought stress is largely unknown. Here, we characterized the role of Arabidopsis SR1 in drought stress tolerance and ABA response by analyzing the phenotypes of SR1 knockout and SR1-overexpression plants. sr1 mutants which accumulate salicylic acid (SA) were found more sensitive to drought stress and showed a higher water loss rate as compared with wild-type. By contrast, SR1-overexpression lines exhibited increased drought tolerance and less water loss than wild-type. Furthermore, sr1 mutants showed reduced ABA response in seed germination, root elongation, and stomatal closure, while SR1-overexpression lines displayed more sensitive to ABA than wild-type. In addition, the drought-sensitive and ABA-insensitive phenotypes of sr1 mutants were recovered by diminishing SA accumulation via knockouts of SA synthesizer ICS1 or activator PAD4, or through expression of SA-degrading enzyme NahG. Some drought/ABA-responsive genes exhibited differentially expressed in sr1 mutants and SR1-overexpression plants. These results suggest that SR1 plays a positive role in drought stress tolerance and ABA response, and drought/ABA responses are antagonized by SA accumulation that is negatively regulated by SR1.

SEED LIPID DROPLET PROTEIN1, SEED LIPID DROPLET PROTEIN2, and LIPID DROPLET PLASMA MEMBRANE ADAPTOR mediate lipid droplet-plasma membrane tethering

Membrane contact sites (MCSs) are interorganellar connections that allow for the direct exchange of molecules, such as lipids or Ca2+ between organelles, but can also serve to tether organelles at specific locations within cells. Here, we identified and characterized three proteins of Arabidopsis thaliana that form a lipid droplet (LD)-plasma membrane (PM) tethering complex in plant cells, namely LD-localized SEED LD PROTEIN (SLDP) 1 and SLDP2 and PM-localized LD-PLASMA MEMBRANE ADAPTOR (LIPA). Using proteomics and different protein-protein interaction assays, we show that both SLDPs associate with LIPA. Disruption of either SLDP1 and SLDP2 expression, or that of LIPA, leads to an aberrant clustering of LDs in Arabidopsis seedlings. Ectopic co-expression of one of the SLDPs with LIPA is sufficient to reconstitute LD-PM tethering in Nicotiana tabacum pollen tubes, a cell type characterized by dynamically moving LDs in the cytosolic streaming. Furthermore, confocal laser scanning microscopy revealed both SLDP2.1 and LIPA to be enriched at LD-PM contact sites in seedlings. These and other results suggest that SLDP and LIPA interact to form a tethering complex that anchors a subset of LDs to the PM during post-germinative seedling growth in Arabidopsis.

KNOX II transcription factor HOS59 functions in regulating rice grain size

Plant Knotted1-like homeobox (KNOX) genes encode homeodomain-containing transcription factors. In rice (Oryza sativa L.), little is known about the downstream target genes of KNOX Class II subfamily proteins. Here we generated chromatin immunoprecipitation (ChIP)-sequencing datasets for HOS59, a member of the rice KNOX Class II subfamily, and characterized the genome-wide binding sites of HOS59. We conducted trait ontology (TO) analysis of 9705 identified downstream target genes, and found that multiple TO terms are related to plant structure morphology and stress traits. ChIP-quantitative PCR (qPCR) was conducted to validate some key target genes. Meanwhile, our IP-MS datasets showed that HOS59 was closely associated with BELL family proteins, some grain size regulators (OsSPL13, OsSPL16, OsSPL18, SLG, etc.), and some epigenetic modification factors such as OsAGO4α and OsAGO4β, proteins involved in small interfering RNA-mediated gene silencing. Furthermore, we employed CRISPR/Cas9 editing and transgenic approaches to generate hos59 mutants and overexpression lines, respectively. Compared with wild-type plants, the hos59 mutants have longer grains and increased glume cell length, a loose plant architecture, and drooping leaves, while the overexpression lines showed smaller grain size, erect leaves, and lower plant height. The qRT-PCR results showed that mutation of the HOS59 gene led to upregulation of some grain size-related genes such as OsSPL13, OsSPL18, and PGL2. In summary, our results indicate that HOS59 may be a repressor of the downstream target genes, negatively regulating glume cell length, rice grain size, plant architecture, etc. The identified downstream target genes and possible interaction proteins of HOS59 improve our understanding of the KNOX regulatory networks.

Sucrose synthases are not involved in starch synthesis in Arabidopsis leaves

Many plants accumulate transitory starch reserves in their leaves during the day to buffer their carbohydrate supply against fluctuating light conditions, and to provide carbon and energy for survival at night. It is universally accepted that transitory starch is synthesized from ADP-glucose (ADPG) in the chloroplasts. However, the consensus that ADPG is made in the chloroplasts by ADPG pyrophosphorylase has been challenged by a controversial proposal that ADPG is made primarily in the cytosol, probably by sucrose synthase (SUS), and then imported into the chloroplasts. To resolve this long-standing controversy, we critically re-examined the experimental evidence that appears to conflict with the consensus pathway. We show that when precautions are taken to avoid artefactual changes during leaf sampling, Arabidopsis thaliana mutants that lack SUS activity in mesophyll cells (quadruple sus1234) or have no SUS activity (sextuple sus123456) have wild-type levels of ADPG and starch, while ADPG is 20 times lower in the pgm and adg1 mutants that are blocked in the consensus chloroplastic pathway of starch synthesis. We conclude that the ADPG needed for starch synthesis in leaves is synthesized primarily by ADPG pyrophosphorylase in the chloroplasts.

CYTOKININ RESPONSE FACTOR 2 is involved in modulating the salt stress response

Cytokinin has strong connections to development and a growing role in the abiotic stress response. Here we show that CYTOKININ RESPONSE FACTOR 2 (CRF2) is additionally involved in the salt (NaCl) stress response. CRF2 promoter-GUS expression indicates CRF2 involvement in the response to salt stress as well as the previously known cytokinin response. Interestingly, CRF2 mutant seedlings are quite similar to the wild type (WT) under non-stressed conditions yet have many distinct changes in response to salt stress. Cytokinin levels measured by liquid chromatography-tandem mass spectrometry (LC-MS/MS) that increased in the WT after salt stress are decreased in crf2, potentially from CRF2 regulation of cytokinin biosynthesis genes. Ion content measured by inductively coupled plasma optical emission spectrometry (ICP-OES) was increased in the WT for Na, K, Mn, Ca and Mg after salt stress, whereas the corresponding Ca and Mg increases are lacking in crf2. Many genes examined by RNA-seq analysis were altered transcriptionally by salt stress in both the WT and crf2, yet interestingly approximately one-third of salt-modified crf2 transcripts (2655) showed unique regulation. Different transcript profiles for salt stress in crf2 compared with the WT background was further supported through an examination of co-expressed genes by weighted gene correlation network analysis (WGCMA) and principal component analysis (PCA). Additionally, Gene Ontology (GO) enrichment terms found from salt-treated transcripts revealed most photosynthesis-related terms as only being affected in crf2, leading to an examination of chlorophyll levels and the efficiency of photosystem II (via the ratio of variable fluorescence to maximum fluorescence, F_v /F_m ) as well as physiology after salt treatment. Salt stress-treated crf2 plants had both reduced chlorophyll levels and lower F_v /F_m values compared with the WT, suggesting that CRF2 plays a role in the modulation of salt stress responses linked to photosynthesis.

A native promoter–gene fusion created by CRISPR/Cas9‐mediated genomic deletion offers a transgene‐free method to drive oil accumulation in leaves

Achieving gain‐of‐function phenotypes without inserting foreign DNA is an important challenge for plant biotechnologists. Here, we show that a gene can be brought under the control of a promoter from an upstream gene by deleting the intervening genomic sequence using dual‐guide CRISPR/Cas9. We fuse the promoter of a nonessential photosynthesis‐related gene to DIACYLGLYCEROL ACYLTRANSFERASE 2 (DGAT2) in the lipase‐deficient sugar‐dependent 1 mutant of Arabidopsis thaliana to drive ectopic oil accumulation in leaves. DGAT2 expression is enhanced more than 20‐fold and the triacylglycerol content increases by around 30‐fold. This deletion strategy offers a transgene‐free route to engineering traits that rely on transcriptional gain‐of‐function, such as producing high lipid forage to increase the productivity and sustainability of ruminant farming.

Identification and Characterization of SOG1 (Suppressor of Gamma Response 1) Homologues in Plants Using Data Mining Resources and Gene Expression Profiling

SOG1 (Suppressor of the Gamma response 1) is the master-regulator of plant DNA damage response (DDR), a highly coordinated network of DNA damage sensors, transducers, mediators, and effectors, with highly coordinated activities. SOG1 transcription factor belongs to the NAC/NAM protein family, containing the well-conserved NAC domain and five serine-glutamine (SQ) motifs, preferential targets for phosphorylation by ATM and ATR. So far, the information gathered for the SOG1 function comes from studies on the model plant Arabidopsis thaliana. To expand the knowledge on plant-specific DDR, it is opportune to gather information on other SOG1 orthologues. The current study identified plants where multiple SOG1 homologues are present and evaluated their functions by leveraging the information contained in publicly available transcriptomics databases. This analysis revealed the presence of multiple SOG1 sequences in thirteen plant species, and four (Medicago truncatula, Glycine max, Kalankoe fedtschenkoi, Populus trichocarpa) were selected for gene expression data mining based on database availability. Additionally, M. truncatula seeds and seedlings exposed to treatments known to activate DDR pathways were used to evaluate the expression profiles of MtSOG1a and MtSOG1b. The experimental workflow confirmed the data retrieved from transcriptomics datasets, suggesting that the SOG1 homologues have redundant functions in different plant species.

Deep roots and many branches: Origins of plant-specialized metabolic enzymes in general metabolism

Collectively, plants produce hundreds of thousands of specialized metabolites from simple building blocks such as amino acids, fatty acids, and isoprenoids. As additional specialized metabolic enzymes are described, there is increasing recognition of the importance of cooption of general metabolic enzymes to specialized metabolism by gene duplication, narrowing of expression, and alteration of enzymatic activities. Here, we examine how several classes of enzymes were each coopted multiple times. We demonstrate the simplicity of achieving the synthesis of analogous chemicals by coopting existing enzymes and summarize emerging insights that could inform rational metabolic engineering of both general and specialized metabolic enzymes.

Research progress of aldehyde oxidases in plants

Plant aldehyde oxidases (AOs) are multi-functional enzymes, and they could oxidize abscisic aldehyde into ABA (abscisic acid) or indole acetaldehyde into IAA (indoleacetic acid) as the last step, respectively. AOs can be divided into four groups based on their biochemical and physiological functions. In this review, we summarized the recent studies about AOs in plants including the motif information, biochemical, and physiological functions. Besides their role in phytohormones biosynthesis and stress response, AOs could also involve in reactive oxygen species homeostasis, aldehyde detoxification and stress tolerance.

Apoplastic class III peroxidases PRX62 and PRX69 promote Arabidopsis root hair growth at low temperature

Root Hairs (RHs) growth is influenced by endogenous and by external environmental signals that coordinately regulate its final cell size. We have recently determined that RH growth was unexpectedly boosted when Arabidopsis thaliana seedlings are cultivated at low temperatures. It was proposed that RH growth plasticity in response to low temperature was linked to a reduced nutrient availability in the media. Here, we explore the molecular basis of this RH growth response by using a Genome Wide Association Study (GWAS) approach using Arabidopsis thaliana natural accessions. We identify the poorly characterized PEROXIDASE 62 (PRX62) and a related protein PRX69 as key proteins under moderate low temperature stress. Strikingly, a cell wall protein extensin (EXT) reporter reveals the effect of peroxidase activity on EXT cell wall association at 10 °C in the RH apical zone. Collectively, our results indicate that PRX62, and to a lesser extent PRX69, are key apoplastic PRXs that modulate ROS-homeostasis and cell wall EXT-insolubilization linked to RH elongation at low temperature.

GXP: Analyze and Plot Plant Omics Data in Web Browsers

Next-generation sequencing and metabolomics have become very cost and work efficient and are integrated into an ever-growing number of life science research projects. Typically, established software pipelines analyze raw data and produce quantitative data informing about gene expression or concentrations of metabolites. These results need to be visualized and further analyzed in order to support scientific hypothesis building and identification of underlying biological patterns. Some of these tools already exist, but require installation or manual programming. We developed “Gene Expression Plotter” (GXP), an RNAseq and Metabolomics data visualization and analysis tool entirely running in the user’s web browser, thus not needing any custom installation, manual programming or uploading of confidential data to third party servers. Consequently, upon receiving the bioinformatic raw data analysis of RNAseq or other omics results, GXP immediately enables the user to interact with the data according to biological questions by performing knowledge-driven, in-depth data analyses and candidate identification via visualization and data exploration. Thereby, GXP can support and accelerate complex interdisciplinary omics projects and downstream analyses. GXP offers an easy way to publish data, plots, and analysis results either as a simple exported file or as a custom website. GXP is freely available on GitHub (see introduction)

Chlorophyll and Pheophytin Dephytylating Enzymes Required for Efficient Repair of PSII in Synechococcus elongatus PCC 7942

The Chlorophyll Dephytylase1 (CLD1) and pheophytinase (PPH) proteins of Arabidopsis thaliana are homologous proteins characterized respectively as a dephytylase for chlorophylls (Chls) and pheophytin a (Phein a) and a Phein a-specific dephytylase. Three genes encoding CLD1/PPH homologs (dphA1, dphA2 and dphA3) were found in the genome of the cyanobacterium Synechococcus elongatus PCC 7942 and shown to be conserved in most cyanobacteria. His6-tagged DphA1, DphA2 and DphA3 proteins were expressed in Escherichia coli, purified to near homogeneity, and shown to exhibit significant levels of dephytylase activity for Chl a and Phein a. Each DphA protein showed similar dephytylase activities for Chl a and Phein a, but the three proteins were distinct in their kinetic properties, with DphA3 showing the highest and lowest Vmax and Km values, respectively, among the three. Transcription of dphA1 and dphA3 was enhanced under high-light conditions, whereas that of dphA2 was not affected by the light conditions. None of the dphA single mutants of S. elongatus showed profound growth defects under low (50 µmol photons m-2 s-1) or high (400 µmol photons m-2 s-1) light conditions. The triple dphA mutant did not show obvious growth defects under these conditions, either, but under illumination of 1,000 µmol photons m-2 s-1, the mutant showed more profound growth retardation compared with wild type (WT). The repair of photodamaged photosystem II (PSII) was much slower in the triple mutant than in WT. These results revealed that dephytylation of Chl a or Phein a or of both is required for efficient repair of photodamaged PSII.

Differential Involvement of Arabidopsis β’-COP Isoforms in Plant Development

Coat protein I (COPI) is necessary for intra-Golgi transport and retrograde transport from the Golgi apparatus back to the endoplasmic reticulum. The key component of the COPI coat is the coatomer complex, which is composed of seven subunits (α/β/β’/γ/δ/ε/ζ) and is recruited en bloc from the cytosol onto Golgi membranes. In mammals and yeast, α- and β’-COP WD40 domains mediate cargo-selective interactions with dilysine motifs present in canonical cargoes of COPI vesicles. In contrast to mammals and yeast, three isoforms of β’-COP (β’1-3-COP) have been identified in Arabidopsis. To understand the role of Arabidopsis β’-COP isoforms in plant biology, we have identified and characterized loss-of-function mutants of the three isoforms, and double mutants were also generated. We have found that the trafficking of a canonical dilysine cargo (the p24 family protein p24δ5) is affected in β’-COP double mutants. By western blot analysis, it is also shown that protein levels of α-COP are reduced in the β’-COP double mutants. Although none of the single mutants showed an obvious growth defect, double mutants showed different growth phenotypes. The double mutant analysis suggests that, under standard growth conditions, β’1-COP can compensate for the loss of both β’2-COP and β’3-COP and may have a prominent role during seedling development.

Subcellular Proteomics as a Unified Approach of Experimental Localizations and Computed Prediction Data for Arabidopsis and Crop Plants

In eukaryotic organisms, subcellular protein location is critical in defining protein function and understanding sub-functionalization of gene families. Some proteins have defined locations, whereas others have low specificity targeting and complex accumulation patterns. There is no single approach that can be considered entirely adequate for defining the in vivo location of all proteins. By combining evidence from different approaches, the strengths and weaknesses of different technologies can be estimated, and a location consensus can be built. The Subcellular Location of Proteins in Arabidopsis database ( http://suba.live/ ) combines experimental data sets that have been reported in the literature and is analyzing these data to provide useful tools for biologists to interpret their own data. Foremost among these tools is a consensus classifier (SUBAcon) that computes a proposed location for all proteins based on balancing the experimental evidence and predictions. Further tools analyze sets of proteins to define the abundance of cellular structures. Extending these types of resources to plant crop species has been complex due to polyploidy, gene family expansion and contraction, and the movement of pathways and processes within cells across the plant kingdom. The Crop Proteins of Annotated Location database ( http://crop-pal.org/ ) has developed a range of subcellular location resources including a species-specific voting consensus for 12 plant crop species that offers collated evidence and filters for current crop proteomes akin to SUBA. Comprehensive cross-species comparison of these data shows that the sub-cellular proteomes (subcellulomes) depend only to some degree on phylogenetic relationship and are more conserved in major biosynthesis than in metabolic pathways. Together SUBA and cropPAL created reference subcellulomes for plants as well as species-specific subcellulomes for cross-species data mining. These data collections are increasingly used by the research community to provide a subcellular protein location layer, inform models of compartmented cell function and protein-protein interaction network, guide future molecular crop breeding strategies, or simply answer a specific question-where is my protein of interest inside the cell?

ZmERF21 directly regulates hormone signaling and stress-responsive gene expression to influence drought tolerance in maize seedlings

Drought stress adversely impacts crop development and yield. Maize frequently encounters drought stress during its life cycle. Improvement of drought tolerance is a priority of maize breeding programs. Here, we identified a novel transcription factor encoding gene, APETALA2 (AP2)/Ethylene response factor (ERF), which is tightly associated with drought tolerance in maize seedlings. ZmERF21 is mainly expressed in the root and leaf and it can be highly induced by polyethylene glycol treatment. Genetic analysis showed that the zmerf21 mutant plants displayed a reduced drought tolerance phenotype, accompanied by phenotypical and physiological changes that are commonly observed in drought conditions. Overexpression of ZmERF21 in maize significantly increased the chlorophyll content and activities of antioxidant enzymes under drought conditions. RNA-Seq and DNA affinity purification sequencing analysis further revealed that ZmERF21 may directly regulate the expression of genes related to hormone (ethylene, abscisic acid) and Ca signaling as well as other stress-response genes through binding to the promoters of potential target genes. Our results thereby provided molecular evidence of ZmERF21 is involved in the drought stress response of maize.

Elevated tyrosine results in the cytosolic retention of 3-deoxy-d-arabino-heptulosonate 7-phosphate synthase in Arabidopsis thaliana

The shikimate pathway plays a central role in the biosynthesis of aromatic amino acids and specialized metabolites in plants. The first enzyme, 3-deoxy-d-arabino-heptulosonate 7-phosphate synthase (DAHPS) serves as a key regulatory point for the pathway in various organisms. These enzymes are important in regulating the shikimate pathway in multiple microbial systems. The mechanism of regulation of DAHPS is poorly understood in plants, and the role of tyrosine (Tyr) with respect to the three DAHPS isozymes from Arabidopsis thaliana was investigated. In vitro enzymatic analyses established that Tyr does not function as an allosteric regulator for the A. thaliana DAHPS isozymes. In contrast, Arabidopsis T-DNA insertional mutants for the DAHPS1 locus, dahps1, are hypersensitive to elevated Tyr. Tyr hypersensitivity can be reversed with tryptophan and phenylalanine supplementation, indicating that Tyr is affecting the shikimate pathway flux in the dahps1 mutant. Tyr treatment of Arabidopsis seedlings showed reduced accumulation of overexpressed DAHPS2 in the chloroplast. Further, bimolecular fluorescence complementation studies revealed that DAHPS2 interacts with a 14-3-3 protein in the cytosol, and this interaction is enhanced with Tyr treatment. This interaction with 14-3-3 may retain DAHPS2 in the cytosol, which prevents its ability to function in the chloroplast with elevated Tyr.

Bibenzyl synthesis in Cannabis sativa L

This study focuses on the biosynthesis of a suite of specialized metabolites from Cannabis that are known as the 'bibenzyls'. In planta, bibenzyls accumulate in response to fungal infection and various other biotic stressors; however, it is their widely recognized anti-inflammatory properties in various animal cell models that have garnered recent therapeutic interest. We propose that these compounds are synthesized via a branch point from the core phenylpropanoid pathway in Cannabis, in a three-step sequence. First, various hydroxycinnamic acids are esterified to acyl-coenzyme A (CoA) by a member of the 4-coumarate-CoA ligase family (Cs4CL4). Next, these CoA esters are reduced by two double-bond reductases (CsDBR2 and CsDBR3) that form their corresponding dihydro-CoA derivatives from preferred substrates. Finally, the bibenzyl backbone is completed by a polyketide synthase that specifically condenses malonyl-CoA with these dihydro-hydroxycinnamoyl-CoA derivatives to form two bibenzyl scaffolds: dihydropiceatannol and dihydroresveratrol. Structural determination of this 'bibenzyl synthase' enzyme (CsBBS2) indicates that a narrowing of the hydrophobic pocket surrounding the active site evolved to sterically favor the non-canonical and more flexible dihydro-hydroxycinnamoyl-CoA substrates in comparison with their oxidized relatives. Accordingly, three point mutations that were introduced into CsBBS2 proved sufficient to restore some enzymatic activity with an oxidized substrate, in vitro. Together, the identification of this set of Cannabis enzymes provides a valuable contribution to the growing 'parts prospecting' inventory that supports the rational metabolic engineering of natural product therapeutics.

Multiple haplotype-based analyses provide genetic and evolutionary insights into tomato fruit weight and composition

Improving fruit quality traits such as metabolic composition remains a challenge for tomato breeders. To better understand the genetic architecture of these traits and decipher the demographic history of the loci controlling tomato quality traits, we applied an innovative approach using multiple haplotype-based analyses, aiming to test the potentials of haplotype based study in association and genomic prediction studies. We performed and compared haplotype vs SNP-based associations (hapQTL) with multi-locus mixed model (MLMM), focusing on tomato fruit weight and metabolite contents (i.e. sugars, organic acids and amino acids). Using a panel of 163 tomato accessions genotyped with 5995 SNPs, we detected a total of 784 haplotype blocks, with an average size of haplotype blocks ~58 kb. A total of 108 significant associations for 26 traits were detected thanks to Haplotype/SNP-based Bayes models. Haplotype-based Bayes model (97 associations) outperformed SNP-based Bayes model (50 associations) and MLMM (53 associations) in identifying marker-trait associations as well as in genomic prediction (especially for those traits with moderate to low heritability). To decipher the demographic history, we identified 24 positive selective sweeps using the integrated haplotype score (iHS). Most of the significant associations for tomato quality traits were located within selective sweeps (54.63% and 71.7% in hapQTL and MLMM models, respectively). Promising candidate genes were identified controlling tomato fruit weight and metabolite contents. We thus demonstrated the benefits of using haplotypes for evolutionary and genetic studies, providing novel insights into tomato quality improvement and breeding history.

The CLASSY family controls tissue-specific DNA methylation patterns in Arabidopsis

DNA methylation shapes the epigenetic landscape of the genome, plays critical roles in regulating gene expression, and ensures transposon silencing. As is evidenced by the numerous defects associated with aberrant DNA methylation landscapes, establishing proper tissue-specific methylation patterns is critical. Yet, how such differences arise remains a largely open question in both plants and animals. Here we demonstrate that CLASSY1-4 (CLSY1-4), four locus-specific regulators of DNA methylation, also control tissue-specific methylation patterns, with the most striking pattern observed in ovules where CLSY3 and CLSY4 control DNA methylation at loci with a highly conserved DNA motif. On a more global scale, we demonstrate that specific clsy mutants are sufficient to shift the epigenetic landscape between tissues. Together, these findings reveal substantial epigenetic diversity between tissues and assign these changes to specific CLSY proteins, elucidating how locus-specific targeting combined with tissue-specific expression enables the CLSYs to generate epigenetic diversity during plant development.

CLASSY (CLSY) proteins regulate DNA methylation at specific loci in the Arabidopsis genome. Here the authors show that the CLSYs also control tissue-specific DNA methylation, including at siren loci in ovules, and that the lack of an individual CLSYs can shift the epigenetic landscape between tissues.

A putative SUBTILISIN-LIKE SERINE PROTEASE 1 (SUBSrP1) regulates anther cuticle biosynthesis and panicle development in rice

INTRODUCTION

Panicle abortion is a severe physiological defect and causes a reduction in grain yield.

OBJECTIVES

In this study, we aim to provide the characterization and functional analysis of a mutant apa1331 (apical panicle abortion1331).

METHODS

The isolated mutant from an EMS-mutagenized population was subjected to SSR analysis and Mutmap assay for candidate gene mapping. We performed phenotypic analysis, anthers cross-sections morphology, wax and cutin profiling, biochemical assays and phylogenetic analysis for characterization and evaluation of apa1331. We used CRISPR/Cas9 disruption for functional validation of its candidate gene. Furthermore, comparative RNA-seq and relative expression analysis were performed to get further insights into mechanistic role of the candidate gene.

RESULTS

The anthers from the apical spikelets of apa1331 were degenerated, pollen-less and showed defects in cuticle formation. Transverse sections of apa1331 anthers showed defects in post-meiotic microspore development at stage 8-9. Gas Chromatography showed a significant reduction of wax and cutin in anthers of apa1331 compared to Wildtype (WT). Quantification of H₂O₂ and MDA has indicated the excessive ROS (reactive oxygen species) in apa1331. Trypan blue staining and TUNEL assay revealed cell death and excessive DNA fragmentation in apa1331. Map-based cloning and Mutmap analysis revealed that LOC_Os04g40720, encoding a putative SUBTILISIN-LIKE SERINE PROTEASE (OsSUBSrP1), harbored an SNP (A > G) in apa1331. Phenotypic defects were only seen in apical spikelets due to highest expression of OsSUBSrP1 in upper panicle portion. CRISPR-mediated knock-out lines of OsSUBSrP1 displayed spikelet abortion comparable to apa1331. Global gene expression analysis revealed a significant downregulation of wax and cutin biosynthesis genes.

CONCLUSIONS

Our study reports the novel role of SUBSrP1 in anther cuticle biosynthesis by ROS-mediated programmed cell death in rice.

Genome-Wide Characterization and Expression Analysis of KH Family Genes Response to ABA and SA in Arabidopsis thaliana

K-homologous (KH) family is a type of nucleic acid-binding protein containing the KH domain and has been found to affect splicing and transcriptional regulation. However, KH family genes haven’t been investigated in plant species systematically. In this study, we identified 30 genes that belonged to the KH family based on HMM of the KH domain in Arabidopsis thaliana. Phylogenetic tree analysis showed that the KH family is grouped into three subgroups. Synteny analysis showed that AtKH9 and AtKH29 have the conserved synteny relationship between A. thaliana and the other five species. The AtKH9 and AtKH29 were located in the cytoplasm and nucleus. The seed germination rates of the mutants atkh9 and atkh29 were higher than wild-type after abscisic acid (ABA) and salicylic acid (SA) treatments. In addition, the expression of ABA-related genes, such as ABRE-binding factor 2 (ABF2), ABRE-binding factor 4 (ABF4), and delta 1-pyrroline-5-carboxylate synthase (P5CS), and an SA-related gene pathogenesis-related proteins b (PR1b) were downregulated after ABA and SA treatments, respectively. These results suggested that atkh9 and atkh29 mutants inhibit the effect of ABA and SA on seed germination. In conclusion, our results provide valuable information for further exploration of the function of KH family genes and propose directions and ideas for the identification and characterization of KH family genes in other plants.

Identification of new proteins in mature sieve elements

The phloem enables vascular plants to transport photoassimilates from source tissues to heterotrophic sink tissues. In the phloem, unbroken strings of enucleated sieve elements, which lose the majority of their cellular contents upon maturation, provide a low resistance path for mass flow. The protein machinery in mature sieve elements performs vital functions to maintain the flow, transmit systemic signals and defend the sugar stream against pests. However, our knowledge of this particular protein population is very limited since mature sieve elements are difficult to isolate and not amenable to transcriptomic analysis due to their enucleate nature. Here, we used co-expression analysis and published gene clusters from transcriptomic studies to generate a list of sieve element proteins that potentially survive the enucleation process to reside in mature sieve elements. We selected seven candidates and show that they all localize in sieve elements in Arabidopsis roots and six of them in bolting stems. Our results support the idea that nascent sieve elements prior to enucleation translate part of the protein machinery found in mature sieve elements. Our co-expression list and the publicly available gene clusters expressed in late proto- and meta-phloem sieve elements are valuable resources for uncharacterized genes that may function in mature sieve elements.

A Comprehensive Analysis of Calmodulin-Like Proteins of Glycine max Indicates Their Role in Calcium Signaling and Plant Defense Against Insect Attack

The calcium (Ca²⁺) signaling is a crucial event during plant-herbivore interaction, which involves a transient change in cytosolic Ca²⁺ concentration, which is sensed by Ca²⁺-sensors, and the received message is transduced to downstream target proteins leading to appropriate defense response. Calmodulin-like proteins (CMLs) are calcium-sensing plant-specific proteins. Although CMLs have been identified in a few plants, they remained uncharacterized in leguminous crop plants. Therefore, a wide-range analysis of CMLs of soybean was performed, which identified 41 true CMLs with greater than 50% similarity with Arabidopsis CMLs. The phylogenetic study revealed their evolutionary relatedness with known CMLs. Further, the identification of conserved motifs, gene structure analysis, and identification of cis-acting elements strongly supported their identity as members of this family and their involvement in stress responses. Only a few Glycine max CMLs (GmCMLs) exhibited differential expression in different tissue types, and rest of them had minimal expression. Additionally, differential expression patterns of GmCMLs were observed during Spodoptera litura-feeding, wounding, and signaling compound treatments, indicating their role in plant defense. The three-dimensional structure prediction, identification of interactive domains, and docking with Ca²⁺ ions of S. litura-inducible GmCMLs, indicated their identity as calcium sensors. This study on the characterization of GmCMLs provided insights into their roles in calcium signaling and plant defense during herbivory.

CAND1 is required for pollen viability in Arabidopsis thaliana-a test of the adaptive exchange hypothesis

The dynamic assembly of SKP1•CUL1•F-box protein (SCF) ubiquitin ligases is important for protein ubiquitination and degradation. This process is enabled by CAND1, which exchanges F-box proteins associated with the common CUL1 scaffold, and thereby, recycles the limited CUL1 core and allows diverse F-box proteins to assemble active SCFs. Previous human cell biological and computational studies have led to the adaptive exchange hypothesis, which suggests that the CAND1-mediated exchange confers plasticity on the SCF system, allowing cells to tolerate large variations in F-box protein expression. Here, we tested this hypothesis using Arabidopsis thaliana, a multicellular organism expressing hundreds of F-box protein genes at variable levels in different tissues. The cand1 null mutant in Arabidopsis is viable but produce almost no seeds. Bioinformatic, cell biological, and developmental analyses revealed that the low fertility in the cand1 mutant is associated with cell death in pollen, where the net expression of F-box protein genes is significantly higher than any other Arabidopsis tissue. In addition, we show that the transmission efficiency of the cand1 null allele was reduced through the male but not the female gametophyte. Our results suggest that CAND1 activity is essential in cells or tissues expressing high levels of F-box proteins. This finding is consistent with the proposed adaptive exchange hypothesis, demonstrating the necessity of the evolutionarily conserved CAND1-mediated exchange system in the development of a multicellular organism.

BUMPY STEM Is an Arabidopsis Choline/Ethanolamine Kinase Required for Normal Development and Chilling Responses

Phospholipid biosynthesis is a core metabolic pathway that affects all aspects of plant growth and development. One of the earliest step in this pathway is mediated by choline/ethanolamine kinases (CEKs), enzymes in the Kennedy pathway that catalyze the synthesis of the polar head groups found on the most abundant plant phospholipids. The Arabidopsis genome encodes four CEKs. CEK1-3 have been well characterized using viable mutants while CEK4 encodes an essential gene, making it difficult to characterize its effects on plant development and responses to the environment. We have isolated an EMS-induced allele of CEK4 called bumpy stem (bst). bst plants are viable, allowing the effects of decreased CEK4 function to be characterized throughout the Arabidopsis life cycle. bst mutants have a range of developmental defects including ectopic stem growths at the base of their flowers, reduced fertility, and short roots and stems. They are also sensitive to cold temperatures. Supplementation with choline, phosphocholine, ethanolamine, and phosphoethanolamine rescues bst root phenotypes, highlighting the flow of metabolites between the choline and ethanolamine branches of the Kennedy pathway. The identification of bst and characterization of its phenotypes defines new roles for CEK4 that go beyond its established biochemical function as an ethanolamine kinase.

Advances in Cell Wall Matrix Research with a Focus on Mixed-Linkage Glucan

Mixed β(1,3;1,4)-linkage glucan (MLG) is commonly found in the monocot lineage, at particularly high levels in the Poaceae family, but also in the evolutionally distant genus, Equisetum. MLG has several properties that make it unique from other plant cell wall polysaccharides. It consists of β1,4-linked polymers of glucose interspersed with β1,3-linkages, but the presence of β1,3-linkages provides quite different physical properties compared to its closest form of the cell wall component, cellulose. The mechanisms of MLG biosynthesis have been investigated to understand whether single or multiple enzymes are required to build mixed linkages in the glucan chain. Currently, MLG synthesis by a single enzyme is supported by mutagenesis analyses of cellulose synthase-like F6, the major MLG synthase, but further investigation is needed to gather mechanistic insights. Because of transient accumulation of MLG in elongating cells and vegetative tissues, several hypotheses have been proposed to explain the role of MLG in the plant cell wall. Studies have been carried out to identify gene expression regulators during development and light cycles as well as enzymes involved in MLG organization in the cell wall. A role of MLG as a storage molecule in grains is evident, but the role of MLG in vegetative tissues is still not well understood. Characterization of a cell wall component is difficult due to the complex heterogeneity of the plant cell wall. However, as detailed in this review, recent exciting research has made significant impacts in the understanding of MLG biology in plants.

Topology of the redox network during induction of photosynthesis as revealed by time-resolved proteomics in tobacco

Photosynthetically produced electrons provide energy for various metabolic pathways, including carbon reduction. Four Calvin-Benson cycle enzymes and several other plastid proteins are activated in the light by reduction of specific cysteines via thioredoxins, a family of electron transporters operating in redox regulation networks. How does this network link the photosynthetic chain with cellular metabolism? Using a time-resolved redox proteomic method, we have investigated the redox network in vivo during the dark–to–low light transition. We show that redox states of some thioredoxins follow the photosynthetic linear electron transport rate. While some redox targets have kinetics compatible with an equilibrium with one thioredoxin (TRXf), reduction of other proteins shows specific kinetic limitations, allowing fine-tuning of each redox-regulated step of chloroplast metabolism. We identified five new redox-regulated proteins, including proteins involved in Mg2+ transport and 1O2 signaling. Our results provide a system-level functional view of the photosynthetic redox regulation network.

Temporal Gene Expression in Apical Culms Shows Early Changes in Cell Wall Biosynthesis Genes in Sugarcane

Multiple genes in sugarcane control sucrose accumulation and the biosynthesis of cell wall components; however, it is unclear how these genes are expressed in its apical culms. To better understand this process, we sequenced mRNA from +1 stem internodes collected from four genotypes with different concentrations of soluble solids. Culms were collected at four different time points, ranging from six to 12-month-old plants. Here we show differentially expressed genes related to sucrose metabolism and cell wall biosynthesis, including genes encoding invertases, sucrose synthase and cellulose synthase. Our results showed increased expression of invertases in IN84-58, the genotype with lower sugar and higher fiber content, as well as delayed expression of secondary cell wall-related cellulose synthase for the other genotypes. Interestingly, genes involved with hormone metabolism were differentially expressed across time points in the three genotypes with higher soluble solids content. A similar result was observed for genes controlling maturation and transition to reproductive stages, possibly a result of selection against flowering in sugarcane breeding programs. These results indicate that carbon partitioning in apical culms of contrasting genotypes is mainly associated with differential cell wall biosynthesis, and may include early modifications for subsequent sucrose accumulation. Co-expression network analysis identified transcription factors related to growth and development, showing a probable time shift for carbon partitioning occurred in 10-month-old plants.

Arabidopsis bioinformatics: tools and strategies

The sequencing of the Arabidopsis thaliana genome 21 years ago ushered in the genomics era for plant research. Since then, an incredible variety of bioinformatic tools permit easy access to large repositories of genomic, transcriptomic, proteomic, epigenomic and other '-omic' data. In this review, we cover some more recent tools (and highlight the 'classics') for exploring such data in order to help formulate quality, testable hypotheses, often without having to generate new experimental data. We cover tools for examining gene expression and co-expression patterns, undertaking promoter analyses and gene set enrichment analyses, and exploring protein-protein and protein-DNA interactions. We will touch on tools that integrate different data sets at the end of the article.

The Story of Auxin-Binding Protein 1 (ABP1)

The auxin-binding protein 1 (ABP1) has endured a history of undulating prominence as a candidate receptor for this important phytohormone. Its capacity for binding auxin has not been in doubt, a feature adequately explained by its crystal structure, but any relevance of this to auxin signaling and plant development has been far more demanding to define. Over its research lifetime, it has been associated with many auxin-induced activities, including ion fluxes across the plasma membrane, rearrangement of the cytoskeleton and cell shape, and the abundance of PIN proteins at the plasma membrane via control of endocytosis, all of which required its presence in the apoplast. Yet, ABP1 has a KDEL sequence that targets it to the endoplasmic reticulum, where most of it remains. This mismatch has been more than adequately compensated for by the need for an auxin receptor to account for responses far too rapid to be executed through transcription and translation and the TIR1/AuxIAA coreceptor system. However, discoveries showing that abp1-null mutants are not compromised for auxin signaling or development, that TIR1 or AFB1 are necessarily involved with very rapid responses at the plasma membrane, and that these rapid responses are mediated with intracellular auxin all suggest that ABP1's auxin-binding capacity is not physiologically relevant. Nevertheless, ABP1 is ubiquitous in higher plants and throughout plant tissues. We need to complete its history by defining its function inside plant cells.

MtExpress, a Comprehensive and Curated RNAseq-based Gene Expression Atlas for the Model Legume Medicago truncatula

Although RNA sequencing (RNAseq) has been becoming the main transcriptomic approach in the model legume Medicago truncatula, there is currently no genome-wide gene expression atlas covering the whole set of RNAseq data published for this species. Nowadays, such a tool is highly valuable to provide a global view of gene expression in a wide range of conditions and tissues/organs. Here, we present MtExpress, a gene expression atlas that compiles an exhaustive set of published M. truncatula RNAseq data (https://medicago.toulouse.inrae.fr/MtExpress). MtExpress makes use of recent releases of M. truncatula genome sequence and annotation, as well as up-to-date tools to perform mapping, quality control, statistical analysis and normalization of RNAseq data. MtExpress combines semi-automated pipelines with manual re-labeling and organization of samples to produce an attractive and user-friendly interface, fully integrated with other available Medicago genomic resources. Importantly, MtExpress is highly flexible, in terms of both queries, e.g. allowing searches with gene names and orthologous gene IDs from Arabidopsis and other legume species, and outputs, to customize visualization and redirect gene study to relevant Medicago webservers. Thanks to its semi-automated pipeline, MtExpress will be frequently updated to follow the rapid pace of M. truncatula RNAseq data publications, as well as the constant improvement of genome annotation. MtExpress also hosts legacy GeneChip expression data originally stored in the Medicago Gene Expression Atlas, as a very valuable and complementary resource.

GWAS identifies an ortholog of the rice D11 gene as a candidate gene for grain size in an international collection of hexaploid wheat

Grain size is a key agronomic trait that contributes to grain yield in hexaploid wheat. Grain length and width were evaluated in an international collection of 157 wheat accessions. These accessions were genetically characterized using a genotyping-by-sequencing (GBS) protocol that produced 73,784 single nucleotide polymorphism (SNP) markers. GBS-derived genotype calls obtained on Chinese Spring proved extremely accurate when compared to the reference (> 99.9%) and showed > 95% agreement with calls made at SNP loci shared with the 90 K SNP array on a subset of 71 Canadian wheat accessions for which both types of data were available. This indicates that GBS can yield a large amount of highly accurate SNP data in hexaploid wheat. The genetic diversity analysis performed using this set of SNP markers revealed the presence of six distinct groups within this collection. A GWAS was conducted to uncover genomic regions controlling variation for grain length and width. In total, seven SNPs were found to be associated with one or both traits, identifying three quantitative trait loci (QTLs) located on chromosomes 1D, 2D and 4A. In the vicinity of the peak SNP on chromosome 2D, we found a promising candidate gene (TraesCS2D01G331100), whose rice ortholog (D11) had previously been reported to be involved in the regulation of grain size. These markers will be useful in breeding for enhanced wheat productivity.