MapMan4: A Refined Protein Classification and Annotation Framework Applicable to Multi-Omics Data Analysis

Genomic signatures of the evolution of defence against its natural enemies in the poisonous and medicinal plant Datura stramonium (Solanaceae)

Tropane alkaloids and terpenoids are widely used in the medicine and pharmaceutic industry and evolved as chemical defenses against herbivores and pathogens in the annual herb Datura stramonium (Solanaceae). Here, we present the first draft genomes of two plants from contrasting environments of D. stramonium. Using these de novo assemblies, along with other previously published genomes from 11 Solanaceae species, we carried out comparative genomic analyses to provide insights on the genome evolution of D. stramonium within the Solanaceae family, and to elucidate adaptive genomic signatures to biotic and abiotic stresses in this plant. We also studied, in detail, the evolution of four genes of D. stramonium-Putrescine N-methyltransferase, Tropinone reductase I, Tropinone reductase II and Hyoscyamine-6S-dioxygenase-involved in the tropane alkaloid biosynthesis. Our analyses revealed that the genomes of D. stramonium show signatures of expansion, physicochemical divergence and/or positive selection on proteins related to the production of tropane alkaloids, terpenoids, and glycoalkaloids as well as on R defensive genes and other important proteins related with biotic and abiotic pressures such as defense against natural enemies and drought.

Barley Anther and Meiocyte Transcriptome Dynamics in Meiotic Prophase I

In flowering plants, successful germinal cell development and meiotic recombination depend upon a combination of environmental and genetic factors. To gain insights into this specialized reproductive development program we used short- and long-read RNA-sequencing (RNA-seq) to study the temporal dynamics of transcript abundance in immuno-cytologically staged barley (Hordeum vulgare) anthers and meiocytes. We show that the most significant transcriptional changes in anthers occur at the transition from pre-meiosis to leptotene-zygotene, which is followed by increasingly stable transcript abundance throughout prophase I into metaphase I-tetrad. Our analysis reveals that the pre-meiotic anthers are enriched in long non-coding RNAs (lncRNAs) and that entry to meiosis is characterized by their robust and significant down regulation. Intriguingly, only 24% of a collection of putative meiotic gene orthologs showed differential transcript abundance in at least one stage or tissue comparison. Argonautes, E3 ubiquitin ligases, and lys48 specific de-ubiquitinating enzymes were enriched in prophase I meiocyte samples. These developmental, time-resolved transcriptomes demonstrate remarkable stability in transcript abundance in meiocytes throughout prophase I after the initial and substantial reprogramming at meiosis entry and the complexity of the regulatory networks involved in early meiotic processes.

Dataset for transcriptome and physiological response of mature tomato seed tissues to light and heat during fruit ripening

Seed vigor is an estimate of how successfully a seed lot will establish seedlings under a wide range of environmental conditions, with both the embryo and the surrounding endosperm playing distinct roles in the germination behaviour. Germination and seedling establishment are essential for crop production to be both sustainable and profitable. Seed vigor traits are sequentially acquired during development via genetic programs that are poorly understood, but known to be under the strong influence of environmental conditions. To investigate how light and temperature have an impact on the molecular mechanisms governing seed vigor at harvest, RNA sequencing was performed on Solanum lycopersicum cv. Moneymaker seed tissues (i.e. embryo and endosperm) that were dissected from fruits that were submitted to standard or high temperature and/or standard or dim light. The dataset encompassed a total of 26.5 Gb raw data from mature embryo and endosperm tissues transcriptomes. The raw and mapped reads data on build SL4.0 and annotation ITAG4.0 are available under accession GSE158641 at NCBI Gene Expression Omnibus (GEO) database. Data on seed vigor characteristics are presented together with the differentially expressed gene transcripts. GO and Mapman annotations were generated on ITAG4.0 to analyse this dataset and are provided for datamining future datasets.

Zucchini Plants Alter Gene Expression and Emission of (E)-β-Caryophyllene Following Aphis gossypii Infestation

Zucchini (Cucurbita pepo L.) is widely cultivated in temperate regions. One of the major production challenges is the damage caused by Aphis gossypii (Homoptera: Aphididae), a polyphagous aphid, which can negatively affect its host plant, both directly by feeding and indirectly by vectoring viruses. To gain insights into the transcriptome events that occur during the zucchini-aphid interaction and to understand the early-to-late defense response through gene expression profiles, we performed RNA-sequencing (RNA-Seq) on zucchini leaves challenged by A. gossypii (24, 48, and 96 h post-infestation; hpi). Data analysis indicated a complex and dynamic pattern of gene expression and a transient transcriptional reconfiguration that involved more than 700 differentially expressed genes (DEGs), including a large number of defense-related genes. The down-regulation of key genes of plant immunity, such as leucine-rich repeat (LRR) protein kinases, transcription factors, and genes associated with direct (i.e., protease inhibitors, cysteine peptidases, etc.) and indirect (i.e., terpene synthase) defense responses, suggests the aphid ability to manipulate plant immune responses. We also investigated the emission of volatile organic compounds (VOCs) from infested plants and observed a reduced emission of (E)-β-caryophyllene at 48 hpi, likely the result of aphid effectors, which reflects the down-regulation of two genes involved in the biosynthesis of terpenoids. We showed that (E)-β-caryophyllene emission was modified by the duration of plant infestation and by aphid density and that this molecule highly attracts Aphidius colemani, a parasitic wasp of A. gossypii. With our results we contributed to the identification of genes involved in cucurbit plant interactions with phloem feeders. Our findings may also help pave the way toward developing tolerant zucchini varieties and to identify molecules for sustainable management of harmful insect populations.

Comparative RNA-Seq analysis reveals genes associated with masculinization in female Cannabis sativa

Using RNA profiling, we identified several silver thiosulfate-induced genes that potentially control the masculinization of female Cannabis sativa plants. Genetically female Cannabis sativa plants normally bear female flowers, but can develop male flowers in response to environmental and developmental cues. In an attempt to elucidate the molecular elements responsible for sex expression in C. sativa plants, we developed genetically female lines producing both female and chemically-induced male flowers. Furthermore, we carried out RNA-Seq assays aimed at identifying differentially expressed genes responsible for male flower development in female plants. The results revealed over 10,500 differentially expressed genes, of which around 200 potentially control masculinization of female cannabis plants. These genes include transcription factors and other genes involved in male organ (i.e., anther and pollen) development, as well as genes involved in phytohormone signalling and male-biased phenotypes. The expressions of 15 of these genes were further validated by qPCR assay confirming similar expression patterns to that of RNA-Seq data. These genes would be useful for understanding predisposed plants producing flowers of both sex types in the same plant, and help breeders to regulate the masculinization of female plants through targeted breeding and plant biotechnology.

Salinity Effects on Guard Cell Proteome in Chenopodium quinoa

Epidermal fragments enriched in guard cells (GCs) were isolated from the halophyte quinoa (Chenopodium quinoa Wild.) species, and the response at the proteome level was studied after salinity treatment of 300 mM NaCl for 3 weeks. In total, 2147 proteins were identified, of which 36% were differentially expressed in response to salinity stress in GCs. Up and downregulated proteins included signaling molecules, enzyme modulators, transcription factors and oxidoreductases. The most abundant proteins induced by salt treatment were desiccation-responsive protein 29B (50-fold), osmotin-like protein OSML13 (13-fold), polycystin-1, lipoxygenase, alpha-toxin, and triacylglycerol lipase (PLAT) domain-containing protein 3-like (eight-fold), and dehydrin early responsive to dehydration (ERD14) (eight-fold). Ten proteins related to the gene ontology term “response to ABA” were upregulated in quinoa GC; this included aspartic protease, phospholipase D and plastid-lipid-associated protein. Additionally, seven proteins in the sucrose–starch pathway were upregulated in the GC in response to salinity stress, and accumulation of tryptophan synthase and L-methionine synthase (enzymes involved in the amino acid biosynthesis) was observed. Exogenous application of sucrose and tryptophan, L-methionine resulted in reduction in stomatal aperture and conductance, which could be advantageous for plants under salt stress. Eight aspartic proteinase proteins were highly upregulated in GCs of quinoa, and exogenous application of pepstatin A (an inhibitor of aspartic proteinase) was accompanied by higher oxidative stress and extremely low stomatal aperture and conductance, suggesting a possible role of aspartic proteinase in mitigating oxidative stress induced by saline conditions.

Phosphorylation Site Motifs in Plant Protein Kinases and Their Substrates

Protein phosphorylation is an important cellular regulatory mechanism affecting the activity, localization, conformation, and interaction of proteins. Protein phosphorylation is catalyzed by kinases, and thus kinases are the enzymes regulating cellular signaling cascades. In the model plant Arabidopsis, 940 genes encode for kinases. The substrate proteins of kinases are phosphorylated at defined sites, which consist of common patterns around the phosphorylation site, known as phosphorylation motifs. The discovery of kinase specificity with a preference of phosphorylation of certain motifs and application of such motifs in deducing signaling cascades helped to reveal underlying regulation mechanisms, and facilitated the prediction of kinase-target pairs. In this mini-review, we took advantage of retrieved data as examples to present the functions of kinase families along with their commonly found phosphorylation motifs from their substrates.

Phosphoproteomics Profiling of Receptor Kinase Mutants

The transmembrane receptor kinase family is the largest protein kinase family in Arabidopsis. Many members of this family play critical roles in plant signaling pathways. However, many of these kinases have yet uncharacterized functions and very little is known about the direct substrates of these kinases. We have developed the "ShortPhos" method, an efficient and simple mass spectrometry (MS)-based phosphoproteomics protocol to perform comparative phosphopeptide profiling of knockout mutants of receptor-like kinases. Through this method, we are able to better understand the functional roles of plant kinases in the context of their signaling networks.

A phased genome based on single sperm sequencing reveals crossover pattern and complex relatedness in tea plants

For diploid organisms that are highly heterozygous, a phased haploid genome can greatly aid in functional genomic, population genetic and breeding studies. Based on the genome sequencing of 135 single sperm cells of the elite tea cultivar 'Fudingdabai', we herein phased the genome of Camellia sinensis, one of the most popular beverage crops worldwide. High-resolution genetic and recombination maps of Fudingdabai were constructed, which revealed that crossover (CO) positions were frequently located in the 5' and 3' ends of annotated genes, while CO distributions across the genome were random. The low CO frequency in tea can be explained by strong CO interference, and CO simulation revealed the proportion of interference insensitive CO ranged from 5.2% to 11.7%. We furthermore developed a method to infer the relatedness between tea accessions and detected complex kinship and genetic signatures of 106 tea accessions. Among them, 59 accessions were closely related with Fudingdabai and 31 of them were first-degree relatives. We additionally identified genes displaying allele specific expression patterns between the two haplotypes of Fudingdabai and genes displaying significantly differential expression levels between Fudingdabai and other haplotypes. These results lay the foundation for further investigation of genetic and epigenetic factors underpinning the regulation of gene expression and provide insights into the evolution of tea plants as well as a valuable genetic resource for future breeding efforts.

The Genomic Impact of Mycoheterotrophy in Orchids

Mycoheterotrophic plants have lost the ability to photosynthesize and obtain essential mineral and organic nutrients from associated soil fungi. Despite involving radical changes in life history traits and ecological requirements, the transition from autotrophy to mycoheterotrophy has occurred independently in many major lineages of land plants, most frequently in Orchidaceae. Yet the molecular mechanisms underlying this shift are still poorly understood. A comparison of the transcriptomes of Epipogium aphyllum and Neottia nidus-avis, two completely mycoheterotrophic orchids, to other autotrophic and mycoheterotrophic orchids showed the unexpected retention of several genes associated with photosynthetic activities. In addition to these selected retentions, the analysis of their expression profiles showed that many orthologs had inverted underground/aboveground expression ratios compared to autotrophic species. Fatty acid and amino acid biosynthesis as well as primary cell wall metabolism were among the pathways most impacted by this expression reprogramming. Our study suggests that the shift in nutritional mode from autotrophy to mycoheterotrophy remodeled the architecture of the plant metabolism but was associated primarily with function losses rather than metabolic innovations.

The Regulation of Adaptation to Cold and Drought Stresses in Poa crymophila Keng Revealed by Integrative Transcriptomics and Metabolomics Analysis

Poa crymophila Keng is highly adaptable to long-term low temperature and drought conditions, making it a desirable foraging grass of the Qinghai-Tibet Plateau. Here, the widely targeted metabolomics and comparative transcriptome analyses were utilized for the discovery of metabolites and genes in P. crymophila in response to cold and drought stresses. P. crymophila were exposed to -5°C for 24 h and recovered to 22°C for 48 h, as well as drought for 10 days followed by re-watering for 1 day. In total, 779 metabolic features were assigned to metabolites and 167,845 unigenes were generated. Seventeen compounds showed significant up-regulation (variable importance in project >1) under both stresses in the metabolic profiling, mainly annotated as carbohydrates, flavones, and phenylpropanoids. The genes which were positively correlated with these metabolites were assigned to pathways (sucrose-starch, raffinose, phenylpropanoid, and flavone metabolism) using the Mapman software package. Alpha-amylase, beta-fructofuranosidase, and sugar transport genes degraded the glucose and starch to small molecule sugars for the purpose of osmotic adjustment and to provide more energy for the growth of P. crymophila in an adverse environment. The induction of cinnamoyl-CoA reductase (CCR) and the MYB gene as well as the sharp increase in schizandrin, a kind of lignan, showed that this likely has the closest connection with the tolerance to both stresses. Four significantly induced flavone compounds are probably involved in reducing oxidative damage. Our results indicated that activation of the phenlypropanoid pathway plays the primary role in P. crymophila adapting to harsh environments. This study showed the mechanism of P. crymophila responding to both cold and drought stresses and showed the discovery of a new biological regulator against stresses.

MapMan Visualization of RNA-Seq Data Using Mercator4 Functional Annotations

Plant omics research has advanced to the stage where it is feasible to generate data from multiple samples and multiple time points to gain insight into biological processes. This impressive array of data can prove challenging to interpret. In this chapter, we describe a solution to this problem, consisting of the MapMan transcript visualization application and the associated MapMan4 ontology and Mercator4 online annotation process.

RNA Sequencing Analyses for Deciphering Potato Molecular Responses

Understanding the molecular mechanisms of potato development and responses to environmental stressors is of utmost importance for achieving stable crop yields. RNA sequencing (RNA-Seq) provides an insight into responses of all of the organism genes to the environmental and developmental cues and thus provides insights into underlying modes of action. In this chapter, we guide a researcher through some of the most important steps in the analysis of transcriptomics data. The initial topic of experimental design is followed by a more wet-lab-oriented section on RNA-Seq sample preparation. Next, we present intermediate steps of data retrieval, quality control, mapping, and differential expression of the dataset and a section on how to expose your data to the public (i.e., public repositories) and make it findable, accessible, interoperable, and reusable (FAIR). In the last four sections, we describe specific tools or Web applications, which ease the exploration of generated results in the context of their gene function and network-based visualizations, specifically GoMapMan, GSEA, DiNAR, and Biomine Explorer. All sections are accompanied by potato dataset examples and include general hints and tricks, as well as potato specificities that one should be aware of.

Transcriptomic Analysis of Changes in Gene Expression During Flowering Induction in Sugarcane Under Controlled Photoperiodic Conditions

Flowering is of utmost relevance for the agricultural productivity of the sugarcane bioeconomy, but data and knowledge of the genetic mechanisms underlying its photoperiodic induction are still scarce. An understanding of the molecular mechanisms that regulate the transition from vegetative to reproductive growth in sugarcane could provide better control of flowering for breeding. This study aimed to investigate the transcriptome of +1 mature leaves of a sugarcane cultivar subjected to florally inductive and non-inductive photoperiodic treatments to identify gene expression patterns and molecular regulatory modules. We identified 7,083 differentially expressed (DE) genes, of which 5,623 showed significant identity to other plant genes. Functional group analysis showed differential regulation of important metabolic pathways involved in plant development, such as plant hormones (i.e., cytokinin, gibberellin, and abscisic acid), light reactions, and photorespiration. Gene ontology enrichment analysis revealed evidence of upregulated processes and functions related to the response to abiotic stress, photoprotection, photosynthesis, light harvesting, and pigment biosynthesis, whereas important categories related to growth and vegetative development of plants, such as plant organ morphogenesis, shoot system development, macromolecule metabolic process, and lignin biosynthesis, were downregulated. Also, out of 76 sugarcane transcripts considered putative orthologs to flowering genes from other plants (such as Arabidopsis thaliana, Oryza sativa, and Sorghum bicolor), 21 transcripts were DE. Nine DE genes related to flowering and response to photoperiod were analyzed either at mature or spindle leaves at two development stages corresponding to the early stage of induction and inflorescence primordia formation. Finally, we report a set of flowering-induced long non-coding RNAs and describe their level of conservation to other crops, many of which showed expression patterns correlated against those in the functionally grouped gene network.

PhosPhAt 4.0: An Updated Arabidopsis Database for Searching Phosphorylation Sites and Kinase-Target Interactions

The PhosPhAt 4.0 database contains information on Arabidopsis phosphorylation sites identified by mass spectrometry in large-scale experiments from different research groups. So far PhosPhAt 4.0 has been one of the most significant large-scale data resources for plant phosphorylation studies. Functionalities of the web application, besides display of phosphorylation sites, include phosphorylation site prediction and kinase-target relationships retrieval. Here, we present an overview and user instructions for the PhosPhAt 4.0 database, with strong emphasis on recent renewals regarding protein annotation by SUBA4.0 and Mapman4, and additional phosphorylation site information imported from other databases, such as UniProt. Here, we provide a user guide for the retrieval of phosphorylation motifs from the kinase-target database and how to visualize these results. The improvements incorporated into the PhosPhAt 4.0 database have produced much more functionality and user flexibility for phosphoproteomic analysis.

Far-Red Light-Induced Azolla filiculoides Symbiosis Sexual Reproduction: Responsive Transcripts of Symbiont Nostoc azollae Encode Transporters Whilst Those of the Fern Relate to the Angiosperm Floral Transition

Water ferns of the genus Azolla and the filamentous cyanobacteria Nostoc azollae constitute a model symbiosis that enabled the colonization of the water surface with traits highly desirable for the development of more sustainable crops: their floating mats capture CO₂ and fix N₂ at high rates using light energy. Their mode of sexual reproduction is heterosporous. The regulation of the transition from the vegetative phase to the spore forming phase in ferns is largely unknown, yet a prerequisite for Azolla domestication, and of particular interest as ferns represent the sister lineage of seed plants. Sporocarps induced with far red light could be crossed so as to verify species attribution of strains from the Netherlands but not of the strain from the Anzali lagoon in Iran; the latter strain was assigned to a novel species cluster from South America. Red-dominated light suppresses the formation of dissemination stages in both gametophyte- and sporophyte-dominated lineages of plants, the response likely is a convergent ecological strategy to open fields. FR-responsive transcripts included those from MIKC^C homologues of CMADS1 and miR319-controlled GAMYB transcription factors in the fern, transporters in N. azollae, and ycf2 in chloroplasts. Loci of conserved microRNA (miRNA) in the fern lineage included miR172, yet FR only induced miR529 and miR535, and reduced miR319 and miR159. Phylogenomic analyses of MIKC^C TFs suggested that the control of flowering and flower organ specification may have originated from the diploid to haploid phase transition in the homosporous common ancestor of ferns and seed plants.

Common bean SNP alleles and candidate genes affecting photosynthesis under contrasting water regimes

Water deficit is a major worldwide constraint to common bean (Phaseolus vulgaris L.) production, being photosynthesis one of the most affected physiological processes. To gain insights into the genetic basis of the photosynthetic response of common bean under water-limited conditions, a collection of 158 Portuguese accessions was grown under both well-watered and water-deficit regimes. Leaf gas-exchange parameters were measured and photosynthetic pigments quantified. The same collection was genotyped using SNP arrays, and SNP-trait associations tested considering a linear mixed model accounting for the genetic relatedness among accessions. A total of 133 SNP-trait associations were identified for net CO2 assimilation rate, transpiration rate, stomatal conductance, and chlorophylls a and b, carotenes, and xanthophyll contents. Ninety of these associations were detected under water-deficit and 43 under well-watered conditions, with only two associations common to both treatments. Identified candidate genes revealed that stomatal regulation, protein translocation across membranes, redox mechanisms, hormone, and osmotic stress signaling were the most relevant processes involved in common bean response to water-limited conditions. These candidates are now preferential targets for common bean water-deficit-tolerance breeding. Additionally, new sources of water-deficit tolerance of Andean, Mesoamerican, and admixed origin were detected as accessions valuable for breeding, and not yet explored.

Differentiation of Tracheary Elements in Sugarcane Suspension Cells Involves Changes in Secondary Wall Deposition and Extensive Transcriptional Reprogramming

Plant lignocellulosic biomass, mostly composed of polysaccharide-rich secondary cell walls (SCWs), provides fermentable sugars that may be used to produce biofuels and biomaterials. However, the complex chemical composition and physical structure of SCWs hinder efficient processing of plant biomass. Understanding the molecular mechanisms underlying SCW deposition is, thus, essential to optimize bioenergy feedstocks. Here, we establish a xylogenic culture as a model system to study SCW deposition in sugarcane; the first of its kind in a C4 grass species. We used auxin and brassinolide to differentiate sugarcane suspension cells into tracheary elements, which showed metaxylem-like reticulate or pitted SCW patterning. The differentiation led to increased lignin levels, mainly caused by S-lignin units, and a rise in p-coumarate, leading to increased p-coumarate:ferulate ratios. RNAseq analysis revealed massive transcriptional reprogramming during differentiation, with upregulation of genes associated with cell wall biogenesis and phenylpropanoid metabolism and downregulation of genes related to cell division and primary metabolism. To better understand the differentiation process, we constructed regulatory networks of transcription factors and SCW-related genes based on co-expression analyses. Accordingly, we found multiple regulatory modules that may underpin SCW deposition in sugarcane. Our results provide important insights and resources to identify biotechnological strategies for sugarcane biomass optimization.

Bringing New Methods to the Seed Proteomics Platform: Challenges and Perspectives

For centuries, crop plants have represented the basis of the daily human diet. Among them, cereals and legumes, accumulating oils, proteins, and carbohydrates in their seeds, distinctly dominate modern agriculture, thus play an essential role in food industry and fuel production. Therefore, seeds of crop plants are intensively studied by food chemists, biologists, biochemists, and nutritional physiologists. Accordingly, seed development and germination as well as age- and stress-related alterations in seed vigor, longevity, nutritional value, and safety can be addressed by a broad panel of analytical, biochemical, and physiological methods. Currently, functional genomics is one of the most powerful tools, giving direct access to characteristic metabolic changes accompanying plant development, senescence, and response to biotic or abiotic stress. Among individual post-genomic methodological platforms, proteomics represents one of the most effective ones, giving access to cellular metabolism at the level of proteins. During the recent decades, multiple methodological advances were introduced in different branches of life science, although only some of them were established in seed proteomics so far. Therefore, here we discuss main methodological approaches already employed in seed proteomics, as well as those still waiting for implementation in this field of plant research, with a special emphasis on sample preparation, data acquisition, processing, and post-processing. Thereby, the overall goal of this review is to bring new methodologies emerging in different areas of proteomics research (clinical, food, ecological, microbial, and plant proteomics) to the broad society of seed biologists.

High Resolution Proteome of Lipid Droplets Isolated from the Pennate Diatom Phaeodactylum tricornutum (Bacillariophyceae) Strain pt4 provides mechanistic insights into complex intracellular coordination during nitrogen deprivation

Lipid droplets (LDs) are an organelle conserved amongst all eukaryotes, consisting of a neutral lipid core surrounded by a polar lipid monolayer. Many species of microalgae accumulate LDs in response to stress conditions, such as nitrogen starvation. Here, we report the isolation and proteomic profiling of LD proteins from the model oleaginous pennate diatom Phaeodactylum tricornutum, strain Pt4 (UTEX 646). We also provide a quantitative description of LD morphological ontogeny, and fatty acid content. Novel cell disruption and LD isolation methods, combined with suspension-trapping and nanoflow liquid chromatography coupled to high resolution mass spectrometry, yielded an unprecedented number of LD proteins. Predictive annotation of the LD proteome suggests a broad assemblage of proteins with diverse functions, including lipid metabolism and vesicle trafficking, as well as ribosomal and proteasomal machinery. These proteins provide mechanistic insights into LD processes, and evidence for interactions between LDs and other organelles. We identify for the first time several key steps in diatom LD-associated triacylglycerol biosynthesis. Bioinformatic analyses of the LD proteome suggests multiple protein targeting mechanisms, including amphipathic helices, post-translational modifications, and translocation machinery. This work corroborates recent findings from other strains of P. tricornutum, other diatoms, and other eukaryotic organisms, suggesting that the fundamental proteins orchestrating LDs are conserved, and represent an ancient component of the eukaryotic endomembrane system. We postulate a comprehensive model of nitrogen starvation-induced diatom LDs on a molecular scale, and provide a wealth of candidates for metabolic engineering, with the potential to eventually customize LD contents.

De novo RNA-Seq analysis in sensitive rice cultivar and comparative transcript profiling in contrasting genotypes reveal genetic biomarkers for fluoride-stress response

The fluoride-sensitive indica rice cultivar, IR-64 was subjected to NaF-treatment for 25 days, following which RNA-Seq analysis identified significant up and down regulation of 1,303 and 93 transcripts respectively. Gene ontology (GO) enrichment analysis classified transcripts into groups related to 'cellular part', 'membrane', 'catalytic activity', 'transporter activity', 'binding', 'metabolic processes' and 'cellular processes'. Analysis of differentially expressed genes (DEGs) revealed fluoride-mediated suppression of abscisic acid (ABA) biosynthesis and signaling. Instead, the gibberellin-dependent pathway and signaling via ABA-independent transcription factors (TFs) was activated. Comparative profiling of selected DEGs in IR-64 and fluoride-tolerant variety, Khitish revealed significant cytoskeletal and nucleosomal remodelling, accompanied with escalated levels of autophagy in stressed IR-64 (unlike that in stressed Khitish). Genes associated with ion, solute and xenobiotic transport were strongly up regulated in stressed IR-64, indicating potential fluoride entry through these channels. On the contrary, genes associated with xenobiotic mobility were suppressed in the tolerant cultivar, which restricted bioaccumulation and translocation of fluoride. Pairwise expression profile analysis between stressed IR-64 and Khitish, supported by extensive statistical modelling predicted that fluoride susceptibility was associated with high expression of genes like amino acid transporter, ABC transporter2, CLCd, MFS monosaccharide transporter, SulfT2.1 and PotT2 while fluoride tolerance with high expression of Sweet11.

The plastid proteome of the nonphotosynthetic chlorophycean alga Polytomella parva

The unicellular, free-living, nonphotosynthetic chlorophycean alga Polytomella parva, closely related to Chlamydomonas reinhardtii and Volvox carteri, contains colorless, starch-storing plastids. The P. parva plastids lack all light-dependent processes but maintain crucial metabolic pathways. The colorless alga also lacks a plastid genome, meaning no transcription or translation should occur inside the organelle. Here, using an algal fraction enriched in plastids as well as publicly available transcriptome data, we provide a morphological and proteomic characterization of the P. parva plastid, ultimately identifying several plastid proteins, both by mass spectrometry and bioinformatic analyses. Data are available via ProteomeXchange with identifier PXD022051. Altogether these results led us to propose a plastid proteome for P. parva, i.e., a set of proteins that participate in carbohydrate metabolism; in the synthesis and degradation of starch, amino acids and lipids; in the biosynthesis of terpenoids and tetrapyrroles; in solute transport and protein translocation; and in redox homeostasis. This is the first detailed plastid proteome from a unicellular, free-living colorless alga.

Erythromycin Treatment of Brassica campestris Seedlings Impacts the Photosynthetic and Protein Synthesis Pathways

Erythromycin (Ery) is a commonly used veterinary drug that prevents infections and promotes the growth of farm animals. Ery is often detected in agricultural fields due to the effects of manure application in the ecosystem. However, there is a lack of information on Ery toxicity in crops. In this study, we performed a comparative proteomic analysis to identify the molecular mechanisms of Ery toxicity during seedling growth based on our observation of a decrease in chlorophyll (Chl) contents using Brassica campestris. A total of 452 differentially abundant proteins (DAPs) were identified including a ribulose-1,5-bisphosphate carboxylase (RuBisCO). The proteomic analysis according to gene ontology (GO) classification revealed that many of these DAPs responding to Ery treatment functioned in a cellular process and a metabolic process. The molecular function analysis showed that DAPs classified within catalytic activity were predominantly changed by Ery, including metabolite interconversion enzyme and protein modifying enzyme. An analysis of functional pathways using MapMan revealed that many photosynthesis components were downregulated, whereas many protein biosynthesis components were upregulated. A good relationship was observed between protein and transcript abundance in a photosynthetic pathway, as determined by qPCR analysis. These combined results suggest that Ery affects plant physiological activity by downregulating protein abundance in the photosynthetic pathway.

TOA: A software package for automated functional annotation in non-model plant species

The increase of sequencing capacity provided by high-throughput platforms has made it possible to routinely obtain large sets of genomic and transcriptomic sequences from model and non-model organisms. Subsequent genomic analysis and gene discovery in next-generation sequencing experiments are, however, bottlenecked by functional annotation. One common way to perform functional annotation of sets of sequences obtained from next-generation sequencing experiments, is by searching for homologous sequences and accessing the related functional information deposited in genomic databases. Functional annotation is especially challenging for non-model organisms, like many plant species. In such cases, existing free and commercial general-purpose applications may not offer complete and accurate results. We present TOA (Taxonomy-oriented annotation), a Python-based user-friendly open source application designed to establish functional annotation pipelines geared towards non-model plant species that can run in Linux/Mac computers, HPCs and cloud servers. TOA performs homology searches against proteins stored in the PLAZA databases, NCBI RefSeq Plant, Nucleotide Database and Non-Redundant Protein Sequence Database, and outputs functional information from several ontology systems: Gene Ontology, InterPro, EC, KEGG, Mapman and MetaCyc. The software performance was validated by comparing the runtimes, total number of annotated sequences and accuracy of the functional information obtained for several plant benchmark data sets with TOA and other functional annotation solutions. TOA outperformed the other software in terms of number of annotated sequences and accuracy of the annotation and constitutes a good alternative to improve functional annotation in plants. TOA is especially recommended for gymnosperms or for low quality sequence data sets of non-model plants.

CropPAL for discovering divergence in protein subcellular location in crops to support strategies for molecular crop breeding

Agriculture faces increasing demand for yield, higher plant-derived protein content and diversity while facing pressure to achieve sustainability. Although the genomes of many of the important crops have been sequenced, the subcellular locations of most of the encoded proteins remain unknown or are only predicted. Protein subcellular location is crucial in determining protein function and accumulation patterns in plants, and is critical for targeted improvements in yield and resilience. Integrating location data from over 800 studies for 12 major crop species into the cropPAL2020 data collection showed that while >80% of proteins in most species are not localised by experimental data, combining species data or integrating predictions can help bridge gaps at similar accuracy. The collation and integration of over 61 505 experimental localisations and more than 6 million predictions showed that the relative sizes of the protein catalogues located in different subcellular compartments are comparable between crops and Arabidopsis. A comprehensive cross-species comparison showed that between 50% and 80% of the subcellulomes are conserved across species and that conservation only depends to some degree on the phylogenetic relationship of the species. Protein subcellular locations in major biosynthesis pathways are more often conserved than in metabolic pathways. Underlying this conservation is a clear potential for subcellular diversity in protein location between species by means of gene duplication and alternative splicing. Our cropPAL data set and search platform (https://crop-pal.org) provide a comprehensive subcellular proteomics resource to drive compartmentation-based approaches for improving yield, protein composition and resilience in future crop varieties.

Utilizing PacBio Iso-Seq for Novel Transcript and Gene Discovery of Abiotic Stress Responses in Oryza sativa L

The wide natural variation present in rice is an important source of genes to facilitate stress tolerance breeding. However, identification of candidate genes from RNA-Seq studies is hampered by the lack of high-quality genome assemblies for the most stress tolerant cultivars. A more targeted solution is the reconstruction of transcriptomes to provide templates to map RNA-seq reads. Here, we sequenced transcriptomes of ten rice cultivars of three subspecies on the PacBio Sequel platform. RNA was isolated from different organs of plants grown under control and abiotic stress conditions in different environments. Reconstructed de novo reference transcriptomes resulted in 37,500 to 54,600 plant-specific high-quality isoforms per cultivar. Isoforms were collapsed to reduce sequence redundancy and evaluated, e.g., for protein completeness (BUSCO). About 40% of all identified transcripts were novel isoforms compared to the Nipponbare reference transcriptome. For the drought/heat tolerant aus cultivar N22, 56 differentially expressed genes in developing seeds were identified at combined heat and drought in the field. The newly generated rice transcriptomes are useful to identify candidate genes for stress tolerance breeding not present in the reference transcriptomes/genomes. In addition, our approach provides a cost-effective alternative to genome sequencing for identification of candidate genes in highly stress tolerant genotypes.

High light intensity aggravates latent manganese deficiency in maize

Manganese (Mn) plays an important role in the oxygen-evolving complex, where energy from light absorption is used for water splitting. Although changes in light intensity and Mn status can interfere with the functionality of the photosynthetic apparatus, the interaction between these two factors and the underlying mechanisms remain largely unknown. Here, maize seedlings were grown hydroponically and exposed to two different light intensities under Mn-sufficient or -deficient conditions. No visual Mn deficiency symptoms appeared even though the foliar Mn concentration in the Mn-deficient treatments was reduced to 2 µg g-1. However, the maximum quantum yield efficiency of PSII and the net photosynthetic rate declined significantly, indicating latent Mn deficiency. The reduction in photosynthetic performance by Mn depletion was further aggravated when plants were exposed to high light intensity. Integrated transcriptomic and proteomic analyses showed that a considerable number of genes encoding proteins in the photosynthetic apparatus were only suppressed by a combination of Mn deficiency and high light, thus indicating interactions between changes in Mn nutritional status and light intensity. We conclude that high light intensity aggravates latent Mn deficiency in maize by interfering with the abundance of PSII proteins.

Differential expression in leaves of Saccharum genotypes contrasting in biomass production provides evidence of genes involved in carbon partitioning

BACKGROUND

The development of biomass crops aims to meet industrial yield demands, in order to optimize profitability and sustainability. Achieving these goals in an energy crop like sugarcane relies on breeding for sucrose accumulation, fiber content and stalk number. To expand the understanding of the biological pathways related to these traits, we evaluated gene expression of two groups of genotypes contrasting in biomass composition.

RESULTS

First visible dewlap leaves were collected from 12 genotypes, six per group, to perform RNA-Seq. We found a high number of differentially expressed genes, showing how hybridization in a complex polyploid system caused extensive modifications in genome functioning. We found evidence that differences in transposition and defense related genes may arise due to the complex nature of the polyploid Saccharum genomes. Genotypes within both biomass groups showed substantial variability in genes involved in photosynthesis. However, most genes coding for photosystem components or those coding for phosphoenolpyruvate carboxylases (PEPCs) were upregulated in the high biomass group. Sucrose synthase (SuSy) coding genes were upregulated in the low biomass group, showing that this enzyme class can be involved with sucrose synthesis in leaves, similarly to sucrose phosphate synthase (SPS) and sucrose phosphate phosphatase (SPP). Genes in pathways related to biosynthesis of cell wall components and expansins coding genes showed low average expression levels and were mostly upregulated in the high biomass group.

CONCLUSIONS

Together, these results show differences in carbohydrate synthesis and carbon partitioning in the source tissue of distinct phenotypic groups. Our data from sugarcane leaves revealed how hybridization in a complex polyploid system resulted in noticeably different transcriptomic profiles between contrasting genotypes.

Oxford Nanopore sequencing: new opportunities for plant genomics?

DNA sequencing was dominated by Sanger's chain termination method until the mid-2000s, when it was progressively supplanted by new sequencing technologies that can generate much larger quantities of data in a shorter time. At the forefront of these developments, long-read sequencing technologies (third-generation sequencing) can produce reads that are several kilobases in length. This greatly improves the accuracy of genome assemblies by spanning the highly repetitive segments that cause difficulty for second-generation short-read technologies. Third-generation sequencing is especially appealing for plant genomes, which can be extremely large with long stretches of highly repetitive DNA. Until recently, the low basecalling accuracy of third-generation technologies meant that accurate genome assembly required expensive, high-coverage sequencing followed by computational analysis to correct for errors. However, today's long-read technologies are more accurate and less expensive, making them the method of choice for the assembly of complex genomes. Oxford Nanopore Technologies (ONT), a third-generation platform for the sequencing of native DNA strands, is particularly suitable for the generation of high-quality assemblies of highly repetitive plant genomes. Here we discuss the benefits of ONT, especially for the plant science community, and describe the issues that remain to be addressed when using ONT for plant genome sequencing.

ddRAD sequencing-based genotyping for population structure analysis in cultivated tomato provides new insights into the genomic diversity of Mediterranean 'da serbo' type long shelf-life germplasm

Double digest restriction-site associated sequencing (ddRAD-seq) is a flexible and cost-effective strategy for providing in-depth insights into the genetic architecture of germplasm collections. Using this methodology, we investigated the genomic diversity of a panel of 288 diverse tomato (Solanum lycopersicum L.) accessions enriched in 'da serbo' (called 'de penjar' in Spain) long shelf life (LSL) materials (152 accessions) mostly originating from Italy and Spain. The rest of the materials originate from different countries and include landraces for fresh consumption, elite cultivars, heirlooms, and breeding lines. Apart from their LSL trait, 'da serbo' landraces are of remarkable interest for their resilience. We identified 32,799 high-quality SNPs, which were used for model ancestry population structure and non-parametric hierarchical clustering. Six genetic subgroups were revealed, clearly separating most 'da serbo' landraces, but also the Spanish germplasm, suggesting a subdivision of the population based on type and geographical provenance. Linkage disequilibrium (LD) in the collection decayed very rapidly within <5 kb. We then investigated SNPs showing contrasted minor frequency allele (MAF) in 'da serbo' materials, resulting in the identification of high frequencies in this germplasm of several mutations in genes related to stress tolerance and fruit maturation such as CTR1 and JAR1. Finally, a mini-core collection of 58 accessions encompassing most of the diversity was selected for further exploitation of key traits. Our findings suggest the presence of a genetic footprint of the 'da serbo' germplasm selected in the Mediterranean basin. Moreover, we provide novel insights on LSL 'da serbo' germplasm as a promising source of alleles for tolerance to stresses.

Molecular Remodeling in Populus PdKOR RNAi Roots Profiled Using LC-MS/MS Proteomics

Plant endo-β-1,4-glucanases belonging to the Glycoside Hydrolase Family 9 have functional roles in cell wall biosynthesis and remodeling via endohydrolysis of (1→4)-β-d-glucosidic linkages. Modification of cell wall chemistry via RNA interference (RNAi)-mediated downregulation of Populus deltoides KORRIGAN (PdKOR), an endo-β-1,4-glucanase familygene was shown to have functional consequences on the composition of secondary metabolome and the ability of modified roots to interact with beneficial microbes. The molecular remodeling that underlies the observed differences at metabolic, physiological, and morphological levels in roots is not well understood. Here a liquid chromatography (LC)-tandem mass spectrometry (MS/MS)-based proteome profiling approach is used to survey the molecular remodeling in root tissues of PdKOR and control plants. A total of 14316 peptides are identified and these mapped to 7139 P. deltoides proteins. Based on 90% sequence identity, the measured protein accessions represent 1187 functional protein groups. Analysis of Gene Ontology (GO) categories and specific individual proteins show differential expression of proteins relevant to plant-microbe interactions, cell wall chemistry, and metabolism. The new proteome dataset serves as a useful resource for deriving new hypotheses and empirical testing pertaining to functional roles of proteins and pathways in differential priming of plant roots to interactions with microbes.

Transcriptomic analysis at organ and time scale reveals gene regulatory networks controlling the sulfate starvation response of Solanum lycopersicum

BACKGROUND

Sulfur is a major component of biological molecules and thus an essential element for plants. Deficiency of sulfate, the main source of sulfur in soils, negatively influences plant growth and crop yield. The effect of sulfate deficiency on plants has been well characterized at the physiological, transcriptomic and metabolomic levels in Arabidopsis thaliana and a limited number of crop plants. However, we still lack a thorough understanding of the molecular mechanisms and regulatory networks underlying sulfate deficiency in most plants. In this work we analyzed the impact of sulfate starvation on the transcriptome of tomato plants to identify regulatory networks and key transcriptional regulators at a temporal and organ scale.

RESULTS

Sulfate starvation reduces the growth of roots and leaves which is accompanied by major changes in the organ transcriptome, with the response being temporally earlier in roots than leaves. Comparative analysis showed that a major part of the Arabidopsis and tomato transcriptomic response to sulfate starvation is conserved between these plants and allowed for the identification of processes specifically regulated in tomato at the transcript level, including the control of internal phosphate levels. Integrative gene network analysis uncovered key transcription factors controlling the temporal expression of genes involved in sulfate assimilation, as well as cell cycle, cell division and photosynthesis during sulfate starvation in tomato roots and leaves. Interestingly, one of these transcription factors presents a high identity with SULFUR LIMITATION1, a central component of the sulfate starvation response in Arabidopsis.

CONCLUSIONS

Together, our results provide the first comprehensive catalog of sulfate-responsive genes in tomato, as well as novel regulatory targets for future functional analyses in tomato and other crops.

Characterization of In Vivo Function(s) of Members of the Plant Mitochondrial Carrier Family

Although structurally related, mitochondrial carrier family (MCF) proteins catalyze the specific transport of a range of diverse substrates including nucleotides, amino acids, dicarboxylates, tricarboxylates, cofactors, vitamins, phosphate and H⁺. Despite their name, they do not, however, always localize to the mitochondria, with plasma membrane, peroxisomal, chloroplast and thylakoid and endoplasmic reticulum localizations also being reported. The existence of plastid-specific MCF proteins is suggestive that the evolution of these proteins occurred after the separation of the green lineage. That said, plant-specific MCF proteins are not all plastid-localized, with members also situated at the endoplasmic reticulum and plasma membrane. While by no means yet comprehensive, the in vivo function of a wide range of these transporters is carried out here, and we discuss the employment of genetic variants of the MCF as a means to provide insight into their in vivo function complementary to that obtained from studies following their reconstitution into liposomes.

Systematic analysis of 1298 RNA-Seq samples and construction of a comprehensive soybean (Glycine max) expression atlas

Soybean (Glycine max [L.] Merr.) is a major crop in animal feed and human nutrition, mainly for its rich protein and oil contents. The remarkable rise in soybean transcriptome studies over the past 5 years generated an enormous amount of RNA-seq data, encompassing various tissues, developmental conditions and genotypes. In this study, we have collected data from 1298 publicly available soybean transcriptome samples, processed the raw sequencing reads and mapped them to the soybean reference genome in a systematic fashion. We found that 94% of the annotated genes (52 737/56 044) had detectable expression in at least one sample. Unsupervised clustering revealed three major groups, comprising samples from aerial, underground and seed/seed-related parts. We found 452 genes with uniform and constant expression levels, supporting their roles as housekeeping genes. On the other hand, 1349 genes showed heavily biased expression patterns towards particular tissues. A transcript-level analysis revealed that 95% (70 963 of 74 490) of the assembled transcripts have intron chains exactly matching those from known transcripts, whereas 3256 assembled transcripts represent potentially novel splicing isoforms. The dataset compiled here constitute a new resource for the community, which can be downloaded or accessed through a user-friendly web interface at http://venanciogroup.uenf.br/resources/. This comprehensive transcriptome atlas will likely accelerate research on soybean genetics and genomics.

Slower development of lower canopy beans produces better coffee

The production of high-quality coffee is being challenged by changing climates in coffee-growing regions. The coffee beans from the upper and lower canopy at different development stages of the same plants were analyzed to investigate the impact of the microenvironment on gene expression and coffee quality. Compared with coffee beans from the upper canopy, lower canopy beans displayed more intense aroma with higher caffeine, trigonelline, and sucrose contents, associated with greater gene expression in the representative metabolic pathways. Global gene expression indicated a longer ripening in the lower canopy, resulting from higher expression of genes relating to growth inhibition and suppression of chlorophyll degradation during early bean ripening. Selection of genotypes or environments that enhance expression of the genes slowing bean development may produce higher quality coffee beans, allowing coffee production in a broader range of available future environments.

Specialized phenolic compounds in seeds: structures, functions, and regulations

Plants produce a huge diversity of specialized metabolites (SM) throughout their life cycle that play important physiological and ecological functions. SM can protect plants and seeds against diseases, predators, and abiotic stresses, or support their interactions with beneficial or symbiotic organisms. They also have strong impacts on human nutrition and health. Despite this importance, the biosynthesis and biological functions of most of the SM remain elusive and their diversity and/or quantity have been reduced in most crops during domestication. Seeds present a large number of SM that are important for their physiological, agronomic, nutritional or industrial qualities and hence, provide interesting models for both studying biosynthesis and producing large amounts of specialized metabolites. For instance, phenolics are abundant and widely distributed in seeds. More specifically, flavonoid pathway has been instrumental for understanding environmental or developmental regulations of specialized metabolic pathways, at the molecular and cellular levels. Here, we summarize current knowledge on seed phenolics as model, and discuss how recent progresses in omics approaches could help to further characterize their diversity, regulations, and the underlying molecular mechanisms involved.

Hayai-Annotation Plants: an ultra-fast and comprehensive functional gene annotation system in plants

Summary

Hayai-Annotation Plants is a browser-based interface for an ultra-fast and accurate functional gene annotation system for plant species using R. The pipeline combines the sequence-similarity searches, using USEARCH against UniProtKB (taxonomy Embryophyta), with a functional annotation step. Hayai-Annotation Plants provides five layers of annotation: i) protein name; ii) gene ontology terms consisting of its three main domains (Biological Process, Molecular Function and Cellular Component); iii) enzyme commission number; iv) protein existence level; and v) evidence type. It implements a new algorithm that gives priority to protein existence level to propagate GO and EC information and annotated Arabidopsis thaliana representative peptide sequences (Araport11) within 5 min at the PC level.

Availability and implementation

The software is implemented in R and runs on Macintosh and Linux systems. It is freely available at https://github.com/kdri-genomics/Hayai-Annotation-Plants under the GPLv3 license.

Supplementary information

Supplementary data are available at Bioinformatics online.

Systematic Multi-Omics Integration (MOI) Approach in Plant Systems Biology

Across all facets of biology, the rapid progress in high-throughput data generation has enabled us to perform multi-omics systems biology research. Transcriptomics, proteomics, and metabolomics data can answer targeted biological questions regarding the expression of transcripts, proteins, and metabolites, independently, but a systematic multi-omics integration (MOI) can comprehensively assimilate, annotate, and model these large data sets. Previous MOI studies and reviews have detailed its usage and practicality on various organisms including human, animals, microbes, and plants. Plants are especially challenging due to large poorly annotated genomes, multi-organelles, and diverse secondary metabolites. Hence, constructive and methodological guidelines on how to perform MOI for plants are needed, particularly for researchers newly embarking on this topic. In this review, we thoroughly classify multi-omics studies on plants and verify workflows to ensure successful omics integration with accurate data representation. We also propose three levels of MOI, namely element-based (level 1), pathway-based (level 2), and mathematical-based integration (level 3). These MOI levels are described in relation to recent publications and tools, to highlight their practicality and function. The drawbacks and limitations of these MOI are also discussed for future improvement toward more amenable strategies in plant systems biology.

Parallel Metabolomic and Transcriptomic Analysis Reveals Key Factors for Quality Improvement of Tea Plants

As one of the most popular beverages globally, tea has enormous economic, cultural, and medicinal importance that necessitates a comprehensive metabolomics study of this species. In this study, a large-scale targeted metabolomics analysis on two types of leaf tissues of nine tea cultivars from five representative geographical origins within China was carried out using the liquid chromatography-mass spectrometry technique. RNA-seq-based transcriptomic analysis was in parallel conducted on the same samples, and gene expression and metabolic differentiation between tissues as well as between the multiple tea cultivars were investigated. The data obtained provide an accessible resource for further studies of naturally occurring metabolic variation of tea plants, which will aid in thoroughly interpreting the underlying genetic and molecular mechanisms of biosynthesis of specialized metabolites in this critical species. Candidate genes including a transcription factor (CsMYB5-like), which were highly correlated with both the content of flavonoids and the expression level of genes participating in the phenylpropanoid and flavonoid biosynthesis pathway, were identified as potential targets for quality improvement of tea.

Phenotyping: New Windows into the Plant for Breeders

Plant phenotyping enables noninvasive quantification of plant structure and function and interactions with environments. High-capacity phenotyping reaches hitherto inaccessible phenotypic characteristics. Diverse, challenging, and valuable applications of phenotyping have originated among scientists, prebreeders, and breeders as they study the phenotypic diversity of genetic resources and apply increasingly complex traits to crop improvement. Noninvasive technologies are used to analyze experimental and breeding populations. We cover the most recent research in controlled-environment and field phenotyping for seed, shoot, and root traits. Select field phenotyping technologies have become state of the art and show promise for speeding up the breeding process in early generations. We highlight the technologies behind the rapid advances in proximal and remote sensing of plants in fields. We conclude by discussing the new disciplines working with the phenotyping community: data science, to address the challenge of generating FAIR (findable, accessible, interoperable, and reusable) data, and robotics, to apply phenotyping directly on farms.

LSTrAP-Cloud: A User-Friendly Cloud Computing Pipeline to Infer Coexpression Networks

As genomes become more and more available, gene function prediction presents itself as one of the major hurdles in our quest to extract meaningful information on the biological processes genes participate in. In order to facilitate gene function prediction, we show how our user-friendly pipeline, the Large-Scale Transcriptomic Analysis Pipeline in Cloud (LSTrAP-Cloud), can be useful in helping biologists make a shortlist of genes involved in a biological process that they might be interested in, by using a single gene of interest as bait. The LSTrAP-Cloud is based on Google Colaboratory, and provides user-friendly tools that process quality-control RNA sequencing data streamed from the European Nucleotide Archive. The LSTRAP-Cloud outputs a gene coexpression network that can be used to identify functionally related genes for any organism with a sequenced genome and publicly available RNA sequencing data. Here, we used the biosynthesis pathway of Nicotiana tabacum as a case study to demonstrate how enzymes, transporters, and transcription factors involved in the synthesis, transport, and regulation of nicotine can be identified using our pipeline.

Sugar Beet (Beta vulgaris) Guard Cells Responses to Salinity Stress: A Proteomic Analysis

Soil salinity is a major environmental constraint affecting crop growth and threatening global food security. Plants adapt to salinity by optimizing the performance of stomata. Stomata are formed by two guard cells (GCs) that are morphologically and functionally distinct from the other leaf cells. These microscopic sphincters inserted into the wax-covered epidermis of the shoot balance CO₂ intake for photosynthetic carbon gain and concomitant water loss. In order to better understand the molecular mechanisms underlying stomatal function under saline conditions, we used proteomics approach to study isolated GCs from the salt-tolerant sugar beet species. Of the 2088 proteins identified in sugar beet GCs, 82 were differentially regulated by salt treatment. According to bioinformatics analysis (GO enrichment analysis and protein classification), these proteins were involved in lipid metabolism, cell wall modification, ATP biosynthesis, and signaling. Among the significant differentially abundant proteins, several proteins classified as "stress proteins" were upregulated, including non-specific lipid transfer protein, chaperone proteins, heat shock proteins, inorganic pyrophosphatase 2, responsible for energized vacuole membrane for ion transportation. Moreover, several antioxidant enzymes (peroxide, superoxidase dismutase) were highly upregulated. Furthermore, cell wall proteins detected in GCs provided some evidence that GC walls were more flexible in response to salt stress. Proteins such as L-ascorbate oxidase that were constitutively high under both control and high salinity conditions may contribute to the ability of sugar beet GCs to adapt to salinity by mitigating salinity-induced oxidative stress.

Transcriptome Analysis of Apple Leaves Infected by the Rust Fungus Gymnosporangium yamadae at Two Sporulation Stages

Apple rust disease caused by Gymnosporangium yamadae is one of the major threats to apple orchards. In this study, dual RNA-seq analysis was conducted to simultaneously monitor gene expression profiles of G. yamadae and infected apple leaves during the formation of rust spermogonia and aecia. The molecular mechanisms underlying this compatible interaction at 10 and 30 days postinoculation (dpi) indicate a significant reaction from the host plant and comprise detoxication pathways at the earliest stage and the induction of secondary metabolism pathways at 30 dpi. Such host reactions have been previously reported in other rust pathosystems and may represent a general reaction to rust infection. G. yamadae transcript profiling indicates a conserved genetic program in spermogonia and aecia that is shared with other rust fungi, whereas secretome prediction reveals the presence of specific secreted candidate effector proteins expressed during apple infection. Unexpectedly, the survey of fungal unigenes in the transcriptome assemblies of inoculated and mock-inoculated apple leaves reveals that G. yamadae infection may modify the fungal community composition in the apple phyllosphere at 30 dpi. Collectively, our results provide novel insights into the compatible apple-G. yamadae interaction and advance the knowledge of this heteroecious demicyclic rust fungus.

Gene expression analysis of Cyanophora paradoxa reveals conserved abiotic stress responses between basal algae and flowering plants

The glaucophyte Cyanophora paradoxa represents the most basal member of the kingdom Archaeplastida, but the function and expression of most of its genes are unknown. This information is needed to uncover how functional gene modules, that is groups of genes performing a given function, evolved in the plant kingdom. We have generated a gene expression atlas capturing responses of Cyanophora to various abiotic stresses. The data were included in the CoNekT-Plants database, enabling comparative transcriptomic analyses across two algae and six land plants. We demonstrate how the database can be used to study gene expression, co-expression networks and gene function in Cyanophora, and how conserved transcriptional programs can be identified. We identified gene modules involved in phycobilisome biosynthesis, response to high light and cell division. While we observed no correlation between the number of differentially expressed genes and the impact on growth of Cyanophora, we found that the response to stress involves a conserved, kingdom-wide transcriptional reprogramming, which is activated upon most stresses in algae and land plants. The Cyanophora stress gene expression atlas and the tools found in the https://conekt.plant.tools/ database thus provide a useful resource to reveal functionally related genes and stress responses in the plant kingdom.

Comparative transcriptome analysis reveals the molecular regulation underlying the adaptive mechanism of cherry (Cerasus pseudocerasus Lindl.) to shelter covering

BACKGROUND

Rain-shelter covering is widely applied during cherry fruit development in subtropical monsoon climates with the aim of decreasing the dropping and cracking of fruit caused by excessive rainfall. Under rain-shelter covering, the characteristics of the leaves and fruit of the cherry plant may adapt to the changes in the microclimate. However, the molecular mechanism underlying such adaptation remains unclear, although clarifying it may be helpful for improving the yield and quality of cherry under rain-shelter covering.

RESULTS

To better understand the regulation and adaptive mechanism of cherry under rain-shelter covering, 38,621 and 3584 differentially expressed genes were identified with a combination of Illumina HiSeq and single-molecule real-time sequencing in leaves and fruits, respectively, at three developmental stages. Among these, key genes, such as those encoding photosynthetic-antenna proteins (Lhca and Lhcb) and photosynthetic electron transporters (PsbP, PsbR, PsbY, and PetF), were up-regulated following the application of rain-shelter covering, leading to increased efficiency of light utilization. The mRNA levels of genes involved in carbon fixation, namely, rbcL and rbcS, were clearly increased compared with those under shelter-free conditions, resulting in improved CO₂ utilization. Furthermore, the transcription levels of genes involved in chlorophyll (hemA, hemN, and chlH) and carotenoid synthesis (crtB, PDS, crtISO, and lcyB) in the sheltered leaves peaked earlier than those in the unsheltered leaves, thereby promoting organic matter accumulation in leaves. Remarkably, the expression levels of key genes involved in the metabolic pathways of phenylpropanoid (PAL, C4H, and 4CL) and flavonoid (CHS, CHI, F3'H, DFR, and ANS) in the sheltered fruits were also up-regulated earlier than of those in the unsheltered fruits, conducive to an increase in anthocyanin content in the fruits.

CONCLUSIONS

According to the physiological indicators and transcriptional expression levels of the related genes, the adaptive regulation mechanism of cherry plants was systematically revealed. These findings can help understand the effect of rain-shelter covering on Chinese cherry cultivation in rainy regions.

Plant Reactome: a knowledgebase and resource for comparative pathway analysis

Plant Reactome (https://plantreactome.gramene.org) is an open-source, comparative plant pathway knowledgebase of the Gramene project. It uses Oryza sativa (rice) as a reference species for manual curation of pathways and extends pathway knowledge to another 82 plant species via gene-orthology projection using the Reactome data model and framework. It currently hosts 298 reference pathways, including metabolic and transport pathways, transcriptional networks, hormone signaling pathways, and plant developmental processes. In addition to browsing plant pathways, users can upload and analyze their omics data, such as the gene-expression data, and overlay curated or experimental gene-gene interaction data to extend pathway knowledge. The curation team actively engages researchers and students on gene and pathway curation by offering workshops and online tutorials. The Plant Reactome supports, implements and collaborates with the wider community to make data and tools related to genes, genomes, and pathways Findable, Accessible, Interoperable and Re-usable (FAIR).