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

Direct and indirect responses of the Arabidopsis transcriptome to an induced increase in trehalose 6-phosphate

Trehalose 6-phosphate (Tre6P) is an essential signal metabolite that regulates the level of sucrose, linking growth and development to the metabolic status. We hypothesized that Tre6P plays a role in mediating the regulation of gene expression by sucrose. To test this, we performed transcriptomic profiling on Arabidopsis (Arabidopsis thaliana) plants that expressed a bacterial TREHALOSE 6-PHOSPHATE SYNTHASE (TPS) under the control of an ethanol-inducible promoter. Induction led to a 4-fold rise in Tre6P levels, a concomitant decrease in sucrose, significant changes (FDR ≤ 0.05) of over 13,000 transcripts, and 2-fold or larger changes of over 5,000 transcripts. Comparison with nine published responses to sugar availability allowed some of these changes to be linked to the rise in Tre6P, while others were probably due to lower sucrose or other indirect effects. Changes linked to Tre6P included repression of photosynthesis-related gene expression and induction of many growth-related processes including ribosome biogenesis. About 500 starvation-related genes are known to be induced by SUCROSE-NON-FERMENTING-1-RELATED KINASE 1 (SnRK1). They were largely repressed by Tre6P in a manner consistent with SnRK1 inhibition by Tre6P. SnRK1 also represses many genes that are involved in biosynthesis and growth. These responded to Tre6P in a more complex manner, pointing toward Tre6P interacting with other C-signaling pathways. Additionally, elevated Tre6P modified the expression of genes encoding regulatory subunits of the SnRK1 complex and TPS class II and FCS-LIKE ZINC FINGER proteins that are thought to modulate SnRK1 function and genes involved in circadian, TARGET OF RAPAMYCIN, light, abscisic acid, and other hormone signaling.

Proteome reorganization and amino acid metabolism during germination and seedling establishment in Lupinus albus

During germination plants rely entirely on their seed storage compounds to provide energy and precursors for the synthesis of macromolecular structures until the seedling has emerged from the soil and photosynthesis can be established. Lupin seeds use proteins as their major storage compounds, accounting for up to 40% of the seed dry weight. Lupins are therefore a valuable complement to soy as a source of plant protein for human and animal nutrition. The aim of this study was to elucidate how storage protein metabolism is coordinated with other metabolic processes to meet the requirements of the growing seedling. In a quantitative approach, we analysed seedling growth, as well as alterations in biomass composition, the proteome, and metabolite profiles during germination and seedling establishment in Lupinus albus. The reallocation of nitrogen resources from seed storage proteins to functional seed proteins was mapped based on a manually curated functional protein annotation database. Although classified as a protein crop, Lupinus albus does not use amino acids as a primary substrate for energy metabolism during germination. However, fatty acid and amino acid metabolism may be integrated at the level of malate synthase to combine stored carbon from lipids and proteins into gluconeogenesis.

Coordinated metabolic adaptation of Arabidopsis thaliana to high light

In plants, exposure to high light irradiation induces various stress responses, which entail complex metabolic rearrangements. To explore these dynamics, we conducted time‐course experiments spanning 2 min to 72 h with Arabidopsis thaliana under high and control light. Comparative metabolomics, transcriptomics, redox proteomics, and stable isotope labeling on leaf rosettes identified a series of synchronous and successive responses that provide a deeper insight into well‐orchestrated mechanisms contributing to high‐light acclimation. We observed transient transcriptome downregulation related to light harvesting and electron flow before the profound remodeling of the photosynthetic apparatus. Throughout the entire time course, redox homeostasis is tightly balanced between downregulation of production and enhanced transformation of NADPH accompanied by redistribution of reducing equivalents across several subcellular compartments. In both light conditions, C4 acids such as malate and fumarate are produced via anaplerosis. In carbon units, their accumulation in vacuoles surpasses plastidic levels of starch and intensifies notably under high light. In parallel, citrate synthesis from pyruvate is significantly hindered diurnally. Isotopic labeling in 2‐oxoglutarate and glutamate suggests a moderate de novo synthesis of C5 acids from a vacuolar citrate reservoir during the light phase while they are largely renewed during the night. In the absence of a diurnal clockwise flow through the tricarboxylic acid (TCA) cycle, increased oxidation of photorespiratory glycine takes over as a source of reductants to fuel mitochondrial ATP production. These findings, along with previous research, contribute to a model integrating redox balance and linking increased carbon assimilation and nitrogen metabolism, especially in the context of an incomplete TCA cycle.

Spatiotemporal transcriptome atlas reveals gene regulatory patterns during the organogenesis of the rapid growing bamboo shoots

Bamboo with its remarkable growth rate and economic significance, offers an ideal system to investigate the molecular basis of organogenesis in rapidly growing plants, particular in monocots, where gene regulatory networks governing the maintenance and differentiation of shoot apical and intercalary meristems remain a subject of controversy. We employed both spatial and single-nucleus transcriptome sequencing on 10× platform to precisely dissect the gene functions in various tissues and early developmental stages of bamboo shoots. Our comprehensive analysis reveals distinct cell trajectories during shoot development, uncovering critical genes and pathways involved in procambium differentiation, intercalary meristem formation, and vascular tissue development. Spatial and temporal expression patterns of key regulatory genes, particularly those related to hormone signaling and lipid metabolism, strongly support the hypothesis that intercalary meristem origin from surrounded parenchyma cells. Specific gene expressions in intercalary meristem exhibit regular and dispersed distribution pattern, offering clues for understanding the intricate molecular mechanisms that drive the rapid growth of bamboo shoots. The single-nucleus and spatial transcriptome analysis reveal a comprehensive landscape of gene activity, enhancing the understanding of the molecular architecture of organogenesis and providing valuable resources for future genomic and genetic studies relying on identities of specific cell types.

Isoform-resolved genome annotation enables mapping of tissue-specific betalain regulation in amaranth

Betalains are coloring pigments produced in some families of the order Caryophyllales, where they replace anthocyanins as coloring pigments. While the betalain pathway itself is well studied, the tissue-specific regulation of the pathway remains mostly unknown. We enhance the high-quality Amaranthus hypochondriacus reference genome and produce a substantially more complete genome annotation, incorporating isoform details. We annotate betalain and anthocyanin pathway genes along with their regulators in amaranth and map the genetic control and tissue-specific regulation of the betalain pathway. Our improved genome annotation allowed us to identify causal mutations that lead to a knock-out of red betacyanins in natural accessions of amaranth. We reveal the tissue-specific regulation of flower color via a previously uncharacterized MYB transcription factor, AhMYB2. Downregulation of AhMYB2 in the flower leads to reduced expression of key betalain enzyme genes and loss of red flower color. Our improved amaranth reference genome represents the most complete genome of amaranth to date and is a valuable resource for betalain and amaranth research. High similarity of the flower betalain regulator AhMYB2 to anthocyanin regulators and a partially conserved interaction motif support the co-option of anthocyanin regulators for the betalain pathway as a possible reason for the mutual exclusiveness of the two pigments.

Widespread changes to the translational landscape in a maize microRNA biogenesis mutant

MicroRNAs are short, non-coding RNAs that repress gene expression in both plants and animals and have diverse functions related to growth, development, and stress responses. The ribonuclease, DICER-LIKE1 (DCL1) is required for two steps in plant miRNA biogenesis: cleavage of the primary miRNAs (pri-miRNAs) to release a hairpin structure, called the precursor miRNA (pre-miRNA) and cleavage of the pre-miRNA to generate the miRNA/miRNA* duplex. The mature miRNA guides the RNA-induced silencing complex to target RNAs with complementary sequences, resulting in translational repression and/or RNA cleavage of target mRNAs. However, the relative contribution of translational repression versus mRNA degradation by miRNAs remains unknown at the genome-level in crops, especially in maize. The maize fuzzy tassel (fzt) mutant contains a hypomorphic mutation in DCL1 resulting in broad developmental defects. While most miRNAs are reduced in fzt, the levels of miRNA-targeted mRNAs are not dramatically increased, suggesting that translational regulation by miRNAs may be common. To gain insight into the repression mechanism of plant miRNAs, we combined ribosome profiling and RNA-sequencing to globally survey miRNA activities in maize. Our data indicate that translational repression contributes significantly to regulation of most miRNA targets and that approximately one-third of miRNA targets are regulated primarily at the translational level. Surprisingly, ribosomes appear altered in fzt mutants suggesting that DCL1 may also have a role in ribosome biogenesis. Thus, DICER-LIKE1 shapes the translational landscape in plants through both miRNA-dependent and miRNA-independent mechanisms.

In-depth exploration of the genomic diversity in tea varieties based on a newly constructed pangenome of Camellia sinensis

Tea, one of the most widely consumed beverages globally, exhibits remarkable genomic diversity in its underlying flavour and health-related compounds. In this study, we present the construction and analysis of a tea pangenome comprising a total of 11 genomes, with a focus on three newly sequenced genomes comprising the purple-leaved assamica cultivar "Zijuan", the temperature-sensitive sinensis cultivar "Anjibaicha" and the wild accession "L618" whose assemblies exhibited excellent quality scores as they profited from latest sequencing technologies. Our analysis incorporates a detailed investigation of transposon complement across the tea pangenome, revealing shared patterns of transposon distribution among the studied genomes and improved transposon resolution with long read technologies, as shown by long terminal repeat (LTR) Assembly Index analysis. Furthermore, our study encompasses a gene-centric exploration of the pangenome, exploring the genomic landscape of the catechin pathway with our study, providing insights on copy number alterations and gene-centric variants, especially for Anthocyanidin synthases. We constructed a gene-centric pangenome by structurally and functionally annotating all available genomes using an identical pipeline, which both increased gene completeness and allowed for a high functional annotation rate. This improved and consistently annotated gene set will allow for a better comparison between tea genomes. We used this improved pangenome to capture the core and dispensable gene repertoire, elucidating the functional diversity present within the tea species. This pangenome resource might serve as a valuable resource for understanding the fundamental genetic basis of traits such as flavour, stress tolerance, and disease resistance, with implications for tea breeding programmes.

A stage-dependent seed defense response to explain efficient seed transmission of Xanthomonas citri pv. fuscans to common bean

Although seed represents an important means of plant pathogen dispersion, the seed-pathogen dialogue remains largely unexplored. A multiomic approach was performed at different seed developmental stages of common bean (Phaseolus vulgaris L.) during asymptomatic colonization by Xanthomonas citri pv. fuscans (Xcf), At the early seed developmental stages, we observed high transcriptional changes both in seeds with bacterial recognition and defense signal transduction genes, and in bacteria with up-regulation of the bacterial type 3 secretion system. This high transcriptional activity of defense genes in Xcf-colonized seeds during maturation refutes the widely diffused assumption considering seeds as passive carriers of microbes. At later seed maturation stages, few transcriptome changes indicated a less intense molecular dialogue between the host and the pathogen, but marked by changes in DNA methylation of plant defense genes, in response to Xcf colonization. We showed examples of pathogen-specific DNA methylations in colonized seeds acting as plant defense silencing to repress plant immune response during the germination process. Finally, we propose a novel plant-pathogen interaction model, specific to the seed tissues, highlighting the existence of distinct phases during seed-pathogen interaction with seeds being actively interacting with colonizing pathogens, then both belligerents switching to more passive mode at later stages.

C4-like Sesuvium sesuvioides (Aizoaceae) exhibits CAM in cotyledons and putative C4-like + CAM metabolism in adult leaves as revealed by transcriptome analysis

BACKGROUND

The co-occurrence of C4 and CAM photosynthesis in a single species seems to be unusual and rare. This is likely due to the difficulty in effectively co-regulating both pathways. Here, we conducted a comparative transcriptomic analysis of leaves and cotyledons of the C4-like species Sesuvium sesuvioides (Aizoaceae) using RNA-seq.

RESULTS

When compared to cotyledons, phosphoenolpyruvate carboxylase 4 (PEPC4) and some key C4 genes were found to be up-regulated in leaves. During the day, the expression of NADP-dependent malic enzyme (NADP-ME) was significantly higher in cotyledons than in leaves. The titratable acidity confirmed higher acidity in the morning than in the previous evening indicating the induction of weak CAM in cotyledons by environmental conditions. Comparison of the leaves of S. sesuvioides (C4-like) and S. portulacastrum (C3) revealed that PEPC1 was significantly higher in S. sesuvioides, while PEPC3 and PEPC4 were up-regulated in S. portulacastrum. Finally, potential key regulatory elements involved in the C4-like and CAM pathways were identified.

CONCLUSIONS

These findings provide a new species in which C4-like and CAM co-occur and raise the question if this phenomenon is indeed so rare or just hard to detect and probably more common in succulent C4 lineages.

Application of orthology and network biology to infer gene functions in non-model plants

Approximately 60% of the genes and gene products in the model species Arabidopsis thaliana have been functionally characterized. In non-model plant species, the functional annotation of the gene space is largely based on homology, with the assumption that genes with shared common ancestry have conserved functions. However, the wide variety in possible morphological, physiological, and ecological differences between plant species gives rise to many species- and clade-specific genes, for which this transfer of knowledge is not possible. Other complications, such as difficulties with genetic transformation, the absence of large-scale mutagenesis methods, and long generation times, further lead to the slow characterization of genes in non-model species. Here, we discuss different resources that integrate plant gene function information. Different approaches that support the functional annotation of gene products, based on orthology or network biology, are described. While sequence-based tools to characterize the functional landscape in non-model species are maturing and becoming more readily available, easy-to-use network-based methods inferring plant gene functions are not as prevalent and have limited functionality.

LSTrAP-denovo: Automated Generation of Transcriptome Atlases for Eukaryotic Species Without Genomes

Despite the abundance of species with transcriptomic data, a significant number of species still lack sequenced genomes, making it difficult to study gene function and expression in these organisms. While de novo transcriptome assembly can be used to assemble protein-coding transcripts from RNA-sequencing (RNA-seq) data, the datasets used often only feature samples of arbitrarily selected or similar experimental conditions, which might fail to capture condition-specific transcripts. We developed the Large-Scale Transcriptome Assembly Pipeline for de novo assembled transcripts (LSTrAP-denovo) to automatically generate transcriptome atlases of eukaryotic species. Specifically, given an NCBI TaxID, LSTrAP-denovo can (1) filter undesirable RNA-seq accessions based on read data, (2) select RNA-seq accessions via unsupervised machine learning to construct a sample-balanced dataset for download, (3) assemble transcripts via over-assembly, (4) functionally annotate coding sequences (CDS) from assembled transcripts and (5) generate transcriptome atlases in the form of expression matrices for downstream transcriptomic analyses. LSTrAP-denovo is easy to implement, written in Python, and is freely available at https://github.com/pengkenlim/LSTrAP-denovo/.

Using RNA sequencing to unravel molecular changes underlying the defense response in chickpea induced by Phytophthora medicaginis

Phytophthora root rot (PRR), caused by Phytophthora medicaginis, is a major soil-borne disease of chickpea in Australia. Breeding for PRR resistance is an effective approach to avoid significant yield loss. Genetic resistance has been identified in cultivated chickpea (Cicer arietinum) and in the wild relative C. echinospermum, with previous studies identifying independent genetic loci associated with each of these sources. However, the molecular mechanisms associated with PRR resistance are not known. RNA sequencing analysis employed in this study identified changes in gene expression in roots of three chickpea genotypes grown hydroponically, early post-infection with P. medicaginis zoospores. Analyses of differentially expressed genes (DEG) identified the activation of a higher number of non-specific R-genes in a PRR-susceptible variety than in the resistant genotypes, suggesting a whole plant resistance response occurring in chickpea against the pathogen. Contrasting molecular changes in signaling profiles, proteolysis and transcription factor pathways were observed in the cultivated and wild Cicer-derived resistant genotypes. DEG patterns supported a hypothesis that increased root elongation and reduced adventitious root formation limit the pathogen entry points in the genotype containing the wild Cicer source of PRR resistance. Candidate resistance genes, including an aquaporin and a maltose transporter in the wild Cicer source and GDSL esterases/lipases in the cultivated source of resistance, were oppositely regulated. Increased knowledge of these genes and pathways will improve our understanding of molecular mechanisms controlling PRR resistance in chickpea, and support the development of elite chickpea varieties through molecular breeding approaches.

Persulfidome of Sweet Pepper Fruits during Ripening: The Case Study of Leucine Aminopeptidase That Is Positively Modulated by H2S

Protein persulfidation is a thiol-based oxidative posttranslational modification (oxiPTM) that involves the modification of susceptible cysteine thiol groups present in peptides and proteins through hydrogen sulfide (H2S), thus affecting their function. Using sweet pepper (Capsicum annuum L.) fruits as a model material at different stages of ripening (immature green and ripe red), endogenous persulfidated proteins (persulfidome) were labeled using the dimedone switch method and identified using liquid chromatography and mass spectrometry analysis (LC-MS/MS). A total of 891 persulfidated proteins were found in pepper fruits, either immature green or ripe red. Among these, 370 proteins were exclusively present in green pepper, 237 proteins were exclusively present in red pepper, and 284 proteins were shared between both stages of ripening. A comparative analysis of the pepper persulfidome with that described in Arabidopsis leaves allowed the identification of 25% of common proteins. Among these proteins, glutathione reductase (GR) and leucine aminopeptidase (LAP) were selected to evaluate the effect of persulfidation using an in vitro approach. GR activity was unaffected, whereas LAP activity increased by 3-fold after persulfidation. Furthermore, this effect was reverted through treatment with dithiothreitol (DTT). To our knowledge, this is the first persulfidome described in fruits, which opens new avenues to study H2S metabolism. Additionally, the results obtained lead us to hypothesize that LAP could be involved in glutathione (GSH) recycling in pepper fruits.

Evidence-based unification of potato gene models with the UniTato collaborative genome browser

Potato (Solanum tuberosum) is the most popular tuber crop and a model organism. A variety of gene models for potato exist, and despite frequent updates, they are not unified. This hinders the comparison of gene models across versions, limits the ability to reuse experimental data without significant re-analysis, and leads to missing or wrongly annotated genes. Here, we unify the recent potato double monoploid v4 and v6 gene models by developing an automated merging protocol, resulting in a Unified poTato genome model (UniTato). We subsequently established an Apollo genome browser (unitato.nib.si) that enables public access to UniTato and further community-based curation. We demonstrate how the UniTato resource can help resolve problems with missing or misplaced genes and can be used to update or consolidate a wider set of gene models or genome information. The automated protocol, genome annotation files, and a comprehensive translation table are provided at github.com/NIB-SI/unitato.

Transcriptional changes during crown-root development and emergence in barley (Hordeum vulgare L.)

BACKGROUND

Roots play an important role during plant growth and development, ensuring water and nutrient uptake. Understanding the mechanisms regulating their initiation and development opens doors towards root system architecture engineering.

RESULTS

Here, we investigated by RNA-seq analysis the changes in gene expression in the barley stem base of 1 day-after-germination (DAG) and 10DAG seedlings when crown roots are formed. We identified 2,333 genes whose expression was lower in the stem base of 10DAG seedlings compared to 1DAG seedlings. Those genes were mostly related to basal cellular activity such as cell cycle organization, protein biosynthesis, chromatin organization, cytoskeleton organization or nucleotide metabolism. In opposite, 2,932 genes showed up-regulation in the stem base of 10DAG seedlings compared to 1DAG seedlings, and their function was related to phytohormone action, solute transport, redox homeostasis, protein modification, secondary metabolism. Our results highlighted genes that are likely involved in the different steps of crown root formation from initiation to primordia differentiation and emergence, and revealed the activation of different hormonal pathways during this process.

CONCLUSIONS

This whole transcriptomic study is the first study aiming at understanding the molecular mechanisms controlling crown root development in barley. The results shed light on crown root emergence that is likely associated with a strong cell wall modification, death of the cells covering the crown root primordium, and the production of defense molecules that might prevent pathogen infection at the site of root emergence.

SOS1 tonoplast neo-localization and the RGG protein SALTY are important in the extreme salinity tolerance of Salicornia bigelovii

The identification of genes involved in salinity tolerance has primarily focused on model plants and crops. However, plants naturally adapted to highly saline environments offer valuable insights into tolerance to extreme salinity. Salicornia plants grow in coastal salt marshes, stimulated by NaCl. To understand this tolerance, we generated genome sequences of two Salicornia species and analyzed the transcriptomic and proteomic responses of Salicornia bigelovii to NaCl. Subcellular membrane proteomes reveal that SbiSOS1, a homolog of the well-known SALT-OVERLY-SENSITIVE 1 (SOS1) protein, appears to localize to the tonoplast, consistent with subcellular localization assays in tobacco. This neo-localized protein can pump Na+ into the vacuole, preventing toxicity in the cytosol. We further identify 11 proteins of interest, of which SbiSALTY, substantially improves yeast growth on saline media. Structural characterization using NMR identified it as an intrinsically disordered protein, localizing to the endoplasmic reticulum in planta, where it can interact with ribosomes and RNA, stabilizing or protecting them during salt stress.

Persulfidation of plant and bacteroid proteins is involved in legume nodule development and senescence

Legumes establish symbiosis with rhizobia, forming nitrogen-fixing nodules. The central role of reactive oxygen species (ROS) and reactive nitrogen species (RNS) in nodule biology has been clearly established. Recently, hydrogen sulfide (H2S) and other reactive sulfur species (RSS) have emerged as novel signaling molecules in animals and plants. A major mechanism by which ROS, RNS, and RSS fulfil their signaling role is the post-translational modification of proteins. To identify possible functions of H2S in nodule development and senescence, we used the tag-switch method to quantify changes in the persulfidation profile of common bean (Phaseolus vulgaris) nodules at different developmental stages. Proteomic analyses indicate that persulfidation plays a regulatory role in plant and bacteroid metabolism and senescence. The effect of a H2S donor on nodule functioning and on several proteins involved in ROS and RNS homeostasis was also investigated. Our results using recombinant proteins and nodulated plants support a crosstalk among H2S, ROS, and RNS, a protective function of persulfidation on redox-sensitive enzymes, and a beneficial effect of H2S on symbiotic nitrogen fixation. We conclude that the general decrease of persulfidation levels observed in plant proteins of aging nodules is one of the mechanisms that disrupt redox homeostasis leading to senescence.

The road less taken: Dihydroflavonol 4-reductase inactivation and delphinidin anthocyanin loss underpins a natural intraspecific flower colour variation

Visual cues are of critical importance for the attraction of animal pollinators, however, little is known about the molecular mechanisms underpinning intraspecific floral colour variation. Here, we combined comparative spectral analysis, targeted metabolite profiling, multi-tissue transcriptomics, differential gene expression, sequence analysis and functional analysis to investigate a bee-pollinated orchid species, Glossodia major with common purple- and infrequent white-flowered morphs. We found uncommon and previously unreported delphinidin-based anthocyanins responsible for the conspicuous and pollinator-perceivable colour of the purple morph and three genetic changes underpinning the loss of colour in the white morph - (1) a loss-of-function (LOF; frameshift) mutation affecting dihydroflavonol 4-reductase (DFR1) coding sequence due to a unique 4-bp insertion, (2) specific downregulation of functional DFR1 expression and (3) the unexpected discovery of chimeric Gypsy transposable element (TE)-gene (DFR) transcripts with potential consequences to the genomic stability and post-transcriptional or epigenetic regulation of DFR. This is one of few known cases where regulatory changes and LOF mutation in an anthocyanin structural gene, rather than transcription factors, are important. Furthermore, if TEs prove to be a frequent source of mutation, the interplay between environmental stress-induced TE evolution and pollinator-mediated selection for adaptive colour variation may be an overlooked mechanism maintaining floral colour polymorphism in nature.

Multiomics analyses reveal the central role of the nucleolus and its machinery during heat stress acclimation in Pinus radiata

Global warming is causing rapid changes in mean annual temperature and more severe drought periods. These are major contributors of forest dieback, which is becoming more frequent and widespread. In this work, we investigated how the transcriptome of Pinus radiata changed during initial heat stress response and acclimation. To this end, we generated a high-density dataset employing Illumina technology. This approach allowed us to reconstruct a needle transcriptome, defining 12 164 and 13 590 transcripts as down- and up-regulated, respectively, during a time course stress acclimation experiment. Additionally, the combination of transcriptome data with other available omics layers allowed us to determine the complex inter-related processes involved in the heat stress response from the molecular to the physiological level. Nucleolus and nucleoid activities seem to be a central core in the acclimating process, producing specific RNA isoforms and other essential elements for anterograde-retrograde stress signaling such as NAC proteins (Pra_vml_051671_1 and Pra_vml_055001_5) or helicase RVB. These mechanisms are connected by elements already known in heat stress response (redox, heat-shock proteins, or abscisic acid-related) and with others whose involvement is not so well defined such as shikimate-related, brassinosteriods, or proline proteases together with their potential regulatory elements. This work provides a first in-depth overview about molecular mechanisms underlying the heat stress response and acclimation in P. radiata.

How does ionizing radiation affect amyloidogenesis in plants?

PURPOSE

Ionizing radiation is a harsh environmental factor that could induce plant senescence. We hypothesized that radiation-related senescence remodels proteome, particularly by triggering the accumulation of prion-like proteins in plant tissues. The object of this study, pea (Pisum sativum L.), is an agriculturally important legume. Research on the functional importance of amyloidogenic proteins was never performed on this species.

MATERIALS AND METHODS

Pea seeds were irradiated in the dose range 5-50 Gy of X-rays. Afterward, Fourier-transform infrared spectroscopy (FTIR) was used to investigate changes in the secondary structure of proteins in germinated 3-day-old seedlings. Specifically, we evaluated the ratio between the amide I and II peaks. Next, we performed protein staining with Congo red to compare the presence of amyloids in the samples. In parallel, we profiled the detergent-resistant proteome fraction by ultrahigh-performance liquid chromatography coupled with tandem mass spectrometry (UHPLC-MS). Differentially accumulated proteins were functionally analyzed in MapMan software, and the PLAAC tool was used to predict putative prion-like proteins.

RESULTS

We showed a reduced germination rate but higher plant height and faster appearance of reproductive organs in the irradiated at dose of 50 Gy group compared with the control; furthermore, we demonstrated more β-sheets and amyloid aggregates in the roots of stressed plants. We detected 531 proteins in detergent-resistant fraction extracted from roots, and 45 were annotated as putative prion-like proteins. Notably, 29 proteins were significantly differentially abundant between the irradiated and the control groups. These proteins belong to several functional categories: amino acid metabolism, carbohydrate metabolism, cytoskeleton organization, regulatory processes, protein biosynthesis, and RNA processing. Thus, the discovery proteomics provided deep data on novel aspects of plant stress biology.

CONCLUSION

Our data hinted that protein accumulation stimulated seedlings' growth as well as accelerated ontogenesis and, eventually, senescence, primarily through translation and RNA processing. The increased abundance of primary metabolism-related proteins indicates more intensive metabolic processes triggered in germinating pea seeds upon X-ray exposure. The functional role of detected putative amyloidogenic proteins should be validated in overexpression or knockout follow-up studies.

Chalkiness and premature controlled by energy homeostasis in OsNAC02 Ko-mutant during vegetative endosperm development

BACKGROUND

Chalkiness is a common phenotype induced by various reasons, such as abiotic stress or the imbalance of starch synthesis and metabolism during the development period. However, the reason mainly for one gene losing its function such as NAC (TFs has a large family in rice) which may cause premature is rarely known to us.

RESULTS

The Ko-Osnac02 mutant demonstrated an obviously early maturation stage compared to the wild type (WT) with 15 days earlier. The result showed that the mature endosperm of Ko-Osnac02 mutant exhibited chalkiness, characterized by white-core and white-belly in mature endosperm. As grain filling rate is a crucial factor in determining the yield and quality of rice (Oryza sativa, ssp. japonica), it's significant that mutant has a lower amylose content (AC) and higher soluble sugar content in the mature endosperm. Interestingly among the top DEGs in the RNA sequencing of N2 (3DAP) and WT seeds revealed that the OsBAM2 (LOC_Os10g32810) expressed significantly high in N2 mutant, which involved in Maltose up-regulated by the starch degradation. As Prediction of Protein interaction showed in the chalky endosperm formation in N2 seeds (3 DAP), seven genes were expressed at a lower-level which should be verified by a heatmap diagrams based on DEGs of N2 versus WT. The Tubulin genes controlling cell cycle are downregulated together with the MCM family genes MCM4 ( ↓), MCM7 ( ↑), which may cause white-core in the early endosperm development. In conclusion, the developing period drastically decreased in the Ko-Osnac02 mutants, which might cause the chalkiness in seeds during the early endosperm development.

CONCLUSIONS

The gene OsNAC02 which controls a great genetic co-network for cell cycle regulation in early development, and KO-Osnac02 mutant shows prematurity and white-core in endosperm.

The genome of Haberlea rhodopensis provides insights into the mechanisms for tolerance to multiple extreme environments

Haberlea rhodopensis, a resurrection species, is the only plant known to be able to survive multiple extreme environments, including desiccation, freezing temperatures, and long-term darkness. However, the molecular mechanisms underlying tolerance to these stresses are poorly studied. Here, we present a high-quality genome of Haberlea and found that ~ 23.55% of the 44,306 genes are orphan. Comparative genomics analysis identified 89 significantly expanded gene families, of which 25 were specific to Haberlea. Moreover, we demonstrated that Haberlea preserves its resurrection potential even in prolonged complete darkness. Transcriptome profiling of plants subjected to desiccation, darkness, and low temperatures revealed both common and specific footprints of these stresses, and their combinations. For example, PROTEIN PHOSPHATASE 2C (PP2C) genes were substantially induced in all stress combinations, while PHYTOCHROME INTERACTING FACTOR 1 (PIF1) and GROWTH RESPONSE FACTOR 4 (GRF4) were induced only in darkness. Additionally, 733 genes with unknown functions and three genes encoding transcription factors specific to Haberlea were specifically induced/repressed upon combination of stresses, rendering them attractive targets for future functional studies. The study provides a comprehensive understanding of the genomic architecture and reports details of the mechanisms of multi-stress tolerance of this resurrection species that will aid in developing strategies that allow crops to survive extreme and multiple abiotic stresses.

Transposable elements contribute to the establishment of the glycine shuttle in Brassicaceae species

C3 -C4 intermediate photosynthesis has evolved at least five times convergently in the Brassicaceae, despite this family lacking bona fide C4 species. The establishment of this carbon concentrating mechanism is known to require a complex suite of ultrastructural modifications, as well as changes in spatial expression patterns, which are both thought to be underpinned by a reconfiguration of existing gene-regulatory networks. However, to date, the mechanisms which underpin the reconfiguration of these gene networks are largely unknown. In this study, we used a pan-genomic association approach to identify genomic features that could confer differential gene expression towards the C3 -C4 intermediate state by analysing eight C3 species and seven C3 -C4 species from five independent origins in the Brassicaceae. We found a strong correlation between transposable element (TE) insertions in cis-regulatory regions and C3 -C4 intermediacy. Specifically, our study revealed 113 gene models in which the presence of a TE within a gene correlates with C3 -C4 intermediate photosynthesis. In this set, genes involved in the photorespiratory glycine shuttle are enriched, including the glycine decarboxylase P-protein whose expression domain undergoes a spatial shift during the transition to C3 -C4 photosynthesis. When further interrogating this gene, we discovered independent TE insertions in its upstream region which we conclude to be responsible for causing the spatial shift in GLDP1 gene expression. Our findings hint at a pivotal role of TEs in the evolution of C3 -C4 intermediacy, especially in mediating differential spatial gene expression.

Divergent roles of IREG/Ferroportin transporters from the nickel hyperaccumulator Leucocroton havanensis

In response to our ever-increasing demand for metals, phytotechnologies are being developed to limit the environmental impact of conventional metal mining. However, the development of these technologies, which rely on plant species able to tolerate and accumulate metals, is partly limited by our lack of knowledge of the underlying molecular mechanisms. In this work, we aimed to better understand the role of metal transporters of the IRON REGULATED 1/FERROPORTIN (IREG/FPN) family from the nickel hyperaccumulator Leucocroton havanensis from the Euphorbiaceae family. Using transcriptomic data, we identified two homologous genes, LhavIREG1 and LhavIREG2, encoding divalent metal transporters of the IREG/FPN family. Both genes are expressed at similar levels in shoots, but LhavIREG1 shows higher expression in roots. The heterologous expression of these transporters in A. thaliana revealed that LhavIREG1 is localized to the plasma membrane, whereas LhavIREG2 is located on the vacuole. In addition, the expression of each gene induced a significant increase in nickel tolerance. Taken together, our data suggest that LhavIREG2 is involved in nickel sequestration in vacuoles of leaf cells, whereas LhavIREG1 is mainly involved in nickel translocation from roots to shoots, but could also be involved in metal sequestration in cell walls. Our results suggest that paralogous IREG/FPN transporters may play complementary roles in nickel hyperaccumulation in plants.

PEO: Plant Expression Omnibus - a comparative transcriptomic database for 103 Archaeplastida

The Plant Expression Omnibus (PEO) is a web application that provides biologists with access to gene expression insights across over 100 plant species, ~60 000 manually annotated RNA-seq samples, and more than 4 million genes. The tool allows users to explore the expression patterns of genes across different organs, identify organ-specific genes, and discover top co-expressed genes for any gene of interest. PEO also provides functional annotations for each gene, allowing for the identification of genetic modules and pathways. PEO is designed to facilitate comparative kingdom-wide gene expression analysis and provide a valuable resource for plant biology research. We provide two case studies to demonstrate the utility of PEO in identifying candidate genes in pollen coat biosynthesis in Arabidopsis and investigating the biosynthetic pathway components of capsaicin in Capsicum annuum. The database is freely available at https://expression.plant.tools/.

Genome-wide profiling of histone (H3) lysine 4 (K4) tri-methylation (me3) under drought, heat, and combined stresses in switchgrass

BACKGROUND

Switchgrass (Panicum virgatum L.) is a warm-season perennial (C4) grass identified as an important biofuel crop in the United States. It is well adapted to the marginal environment where heat and moisture stresses predominantly affect crop growth. However, the underlying molecular mechanisms associated with heat and drought stress tolerance still need to be fully understood in switchgrass. The methylation of H3K4 is often associated with transcriptional activation of genes, including stress-responsive. Therefore, this study aimed to analyze genome-wide histone H3K4-tri-methylation in switchgrass under heat, drought, and combined stress.

RESULTS

In total, ~ 1.3 million H3K4me3 peaks were identified in this study using SICER. Among them, 7,342; 6,510; and 8,536 peaks responded under drought (DT), drought and heat (DTHT), and heat (HT) stresses, respectively. Most DT and DTHT peaks spanned 0 to + 2000 bases from the transcription start site [TSS]. By comparing differentially marked peaks with RNA-Seq data, we identified peaks associated with genes: 155 DT-responsive peaks with 118 DT-responsive genes, 121 DTHT-responsive peaks with 110 DTHT-responsive genes, and 175 HT-responsive peaks with 136 HT-responsive genes. We have identified various transcription factors involved in DT, DTHT, and HT stresses. Gene Ontology analysis using the AgriGO revealed that most genes belonged to biological processes. Most annotated peaks belonged to metabolite interconversion, RNA metabolism, transporter, protein modifying, defense/immunity, membrane traffic protein, transmembrane signal receptor, and transcriptional regulator protein families. Further, we identified significant peaks associated with TFs, hormones, signaling, fatty acid and carbohydrate metabolism, and secondary metabolites. qRT-PCR analysis revealed the relative expressions of six abiotic stress-responsive genes (transketolase, chromatin remodeling factor-CDH3, fatty-acid desaturase A, transmembrane protein 14C, beta-amylase 1, and integrase-type DNA binding protein genes) that were significantly (P < 0.05) marked during drought, heat, and combined stresses by comparing stress-induced against un-stressed and input controls.

CONCLUSION

Our study provides a comprehensive and reproducible epigenomic analysis of drought, heat, and combined stress responses in switchgrass. Significant enrichment of H3K4me3 peaks downstream of the TSS of protein-coding genes was observed. In addition, the cost-effective experimental design, modified ChIP-Seq approach, and analyses presented here can serve as a prototype for other non-model plant species for conducting stress studies.

Transcriptomics, proteomics, and metabolomics interventions prompt crop improvement against metal(loid) toxicity

The escalating challenges posed by metal(loid) toxicity in agricultural ecosystems, exacerbated by rapid climate change and anthropogenic pressures, demand urgent attention. Soil contamination is a critical issue because it significantly impacts crop productivity. The widespread threat of metal(loid) toxicity can jeopardize global food security due to contaminated food supplies and pose environmental risks, contributing to soil and water pollution and thus impacting the whole ecosystem. In this context, plants have evolved complex mechanisms to combat metal(loid) stress. Amid the array of innovative approaches, omics, notably transcriptomics, proteomics, and metabolomics, have emerged as transformative tools, shedding light on the genes, proteins, and key metabolites involved in metal(loid) stress responses and tolerance mechanisms. These identified candidates hold promise for developing high-yielding crops with desirable agronomic traits. Computational biology tools like bioinformatics, biological databases, and analytical pipelines support these omics approaches by harnessing diverse information and facilitating the mapping of genotype-to-phenotype relationships under stress conditions. This review explores: (1) the multifaceted strategies that plants use to adapt to metal(loid) toxicity in their environment; (2) the latest findings in metal(loid)-mediated transcriptomics, proteomics, and metabolomics studies across various plant species; (3) the integration of omics data with artificial intelligence and high-throughput phenotyping; (4) the latest bioinformatics databases, tools and pipelines for single and/or multi-omics data integration; (5) the latest insights into stress adaptations and tolerance mechanisms for future outlooks; and (6) the capacity of omics advances for creating sustainable and resilient crop plants that can thrive in metal(loid)-contaminated environments.

Comparative transcriptomics of the chilling stress response in two Asian mangrove species, Bruguiera gymnorhiza and Rhizophora apiculata

Low temperatures largely determine the geographic limits of plant species by reducing survival and growth. Inter-specific differences in the geographic distribution of mangrove species have been associated with cold tolerance, with exclusively tropical species being highly cold-sensitive and subtropical species being relatively cold-tolerant. To identify species-specific adaptations to low temperatures, we compared the chilling stress response of two widespread Indo-West Pacific mangrove species from Rhizophoraceae with differing latitudinal range limits-Bruguiera gymnorhiza (L.) Lam. ex Savigny (subtropical range limit) and Rhizophora apiculata Blume (tropical range limit). For both species, we measured the maximum photochemical efficiency of photosystem II (Fv/Fm) as a proxy for the physiological condition of the plants and examined gene expression profiles during chilling at 15 and 5 °C. At 15 °C, B. gymnorhiza maintained a significantly higher Fv/Fm than R. apiculata. However, at 5 °C, both species displayed equivalent Fv/Fm values. Thus, species-specific differences in chilling tolerance were only found at 15 °C, and both species were sensitive to chilling at 5 °C. At 15 °C, B. gymnorhiza downregulated genes related to the light reactions of photosynthesis and upregulated a gene involved in cyclic electron flow regulation, whereas R. apiculata downregulated more RuBisCo-related genes. At 5 °C, both species repressed genes related to CO2 assimilation. The downregulation of genes related to light absorption and upregulation of genes related to cyclic electron flow regulation are photoprotective mechanisms that likely contributed to the greater photosystem II photochemical efficiency of B. gymnorhiza at 15 °C. The results of this study provide evidence that the distributional range limits and potentially the expansion rates of plant species are associated with differences in the regulation of photosynthesis and photoprotective mechanisms under low temperatures.

The chromosome-scale genome and the genetic resistance machinery against insect herbivores of the Mexican toloache, Datura stramonium

Plant resistance refers to the heritable ability of plants to reduce damage caused by natural enemies, such as herbivores and pathogens, either through constitutive or induced traits like chemical compounds or trichomes. However, the genetic architecture-the number and genome location of genes that affect plant defense and the magnitude of their effects-of plant resistance to arthropod herbivores in natural populations remains poorly understood. In this study, we aimed to unveil the genetic architecture of plant resistance to insect herbivores in the annual herb Datura stramonium (Solanaceae) through quantitative trait loci mapping. We achieved this by assembling the species' genome and constructing a linkage map using an F2 progeny transplanted into natural habitats. Furthermore, we conducted differential gene expression analysis between undamaged and damaged plants caused by the primary folivore, Lema daturaphila larvae. Our genome assembly resulted in 6,109 scaffolds distributed across 12 haploid chromosomes. A single quantitative trait loci region on chromosome 3 was associated with plant resistance, spanning 0 to 5.17 cM. The explained variance by the quantitative trait loci was 8.44%. Our findings imply that the resistance mechanisms of D. stramonium are shaped by the complex interplay of multiple genes with minor effects. Protein-protein interaction networks involving genes within the quantitative trait loci region and overexpressed genes uncovered the key role of receptor-like cytoplasmic kinases in signaling and regulating tropane alkaloids and terpenoids, which serve as powerful chemical defenses against D. stramonium herbivores. The data generated in our study constitute important resources for delving into the evolution and ecology of secondary compounds mediating plant-insect interactions.

The gain-of-function mutation blf13 in the barley orthologue of the rice growth regulator NARROW LEAF1 is associated with increased leaf width

Canopy architecture in cereals plays an important role in determining yield. Leaf width represents one key aspect of this canopy architecture. However, our understanding of leaf width control in cereals remains incomplete. Classical mutagenesis studies in barely identified multiple morphological mutants, including those with differing leaf widths. Of these, we characterized the broad leaf13 (blf13) mutant in detail. Mutant plants form wider leaves due to increased post-initiation growth and cell proliferation. The mutant phenotype perfectly co-segregated with a missense mutation in the HvHNT1 gene which affected a highly conserved region of the encoded protein, orthologous to the rice NARROW LEAF1 (NAL1) protein. Causality of this mutation for the blf13 phenotype is further supported by correlative transcriptomic analyses and protein-protein interaction studies showing that the mutant HvNHT1 protein interacts more strongly with a known interactor than wild-type HvHNT1. The mutant HvHNT1 protein also showed stronger homodimerization compared with wild-type HvHNT1, and homology modelling suggested an additional interaction site between HvHNT1 monomers due to the blf13 mutation. Thus, the blf13 mutation parallels known gain-of-function NAL1 alleles in rice that increase leaf width and grain yield, suggesting that the blf13 mutation may have a similar agronomic potential in barley.

Transcriptome Analysis and Reactive Oxygen Species Detection Suggest Contrasting Molecular Mechanisms in Populus canadensis' Response to Different Formae Speciales of Marssonina brunnea

Revealing plant-pathogen interactions is important for resistance breeding, but it remains a complex process that presents many challenges. Marssonina leaf spot of poplars (MLSP) is the main disease in poplars; in China, its pathogens consist of two formae speciales, namely, Marssonina brunnea f. sp. Monogermtubi (MO) and M. brunnea f. sp. Multigermtubi (MU). However, the mechanism of the molecular interaction between poplars and the two formae speciales, especially for an incompatible system, remains unclear. In this study, we conducted transcriptome sequencing and reactive oxygen species (ROS) staining based on the interactions between Populus canadensis and the two formae speciales. The results show that the gene expression patterns of P. canadensis induced by MO and MU were significantly different, especially for the genes associated with biotic stress. Furthermore, MO and MU also triggered distinct ROS reactions of P. canadensis, and ROS (mainly H2O2) burst was only observed around the cells penetrated by MU. In conclusion, this study suggested that P. canadensis experienced different resistance reactions in response to the two formae speciales of M. brunnea, providing valuable insights for further understanding the host-pathogen interactions of MLSP.

Global transcriptomic analysis in avocado nursery trees reveals differential gene expression during asymptomatic infection by avocado sunblotch viroid (ASBVd)

Avocado sunblotch viroid (ASBVd) is the type species of the family Avsunviroidae and the causal agent of avocado sunblotch disease. The disease is characterised by the presence of chlorotic lesions on avocado fruit, leaves and/or stems. Infected trees may remain without chlorosis for extended periods of time, though distorted growth and reduced yield has been observed in these cases. The molecular effects of ASBVd on avocado, and members of the Avsunviroidae on their respective hosts in general, remain poorly understood. Host global transcriptomic studies within the family Pospiviroidae have identified several host pathways that are affected during these plant-pathogen interactions. In this study, we used RNA sequencing to investigate host gene expression in asymptomatic avocado nursery trees infected with ASBVd. Transcriptome data showed that 631 genes were differentially expressed, 63 % of which were upregulated during infection. Plant defence responses, phytohormone networks, gene expression pathways, secondary metabolism, cellular transport as well as protein modification and degradation were all significantly affected by ASBVd infection. This work represents the first global gene expression study of ASBVd-infected avocado, and the transcriptional reprogramming observed during this asymptomatic infection improves our understanding of the molecular interactions underlying broader avsunviroid-host interactions.

The nitrate transporter-sensor MtNPF6.8 regulates the branched chain amino acid/pantothenate metabolic pathway in barrel medic (Medicago truncatula Gaertn.) root tip

Nitrogen is the most limiting nutrient for plants, and it is preferentially absorbed in the form of nitrate by roots, which adapt to nitrate fluctuations by remodelling their architecture. Although core mechanisms of the response to nitrate availability are relatively well-known, signalling events controlling root growth and architecture have not all been identified, in particular in Legumes. However, the developmental effect of nitrate in Legumes is critical since external nitrate not only regulates root architecture but also N2-fixing nodule development. We have previously shown that in barrel medic (Medicago truncatula), the nitrate transporter MtNPF6.8 is required for nitrate sensitivity in root tip. However, uncertainty remains as to whether nitrogen metabolism itself is involved in the MtNPF6.8-mediated response. Here, we examine the metabolic effects of MtNPF6.8-dependent nitrate signalling using metabolomics and proteomics in WT and mtnpf6.8 root tips in presence or absence of nitrate. We found a reorchestration of metabolism due to the mutation, in favour of the branched chain amino acids/pantothenate metabolic pathway, and lipid catabolism via glyoxylate. That is, the mtnpf6.8 mutation was likely associated with a specific rerouting of acetyl-CoA production (glyoxylic cycle) and utilisation (pantothenate and branched chain amino acid synthesis). In agreement with our previous findings, class III peroxidases were confirmed as the main protein class responsive to nitrate, although in an MtNPF6.8-independent fashion. Our data rather suggest the involvement of other pathways within mtnpf6.8 root tips, such as Ca2+ signalling or cell wall methylation.

Development and analysis of de novo transcriptome assemblies of multiple genotypes of Cymbopogon spp. reveal candidate genes involved in the biosynthesis of aromatic monoterpenes

Despite the high economic value of the monoterpene-rich essential oils from different genotypes of Cymbopogon, the knowledge about the genes and metabolic route(s) involved in the biosynthesis of aromatic monoterpenes in this genus is limited. In the present study, a comprehensive transcriptome analysis of four genotypes of Cymbopogon, displaying diverse quantitative and qualitative profiles of volatile monoterpenes in their essential oils has been carried out. The comparative analysis of the deduced protein sequences corresponding to the transcriptomes of the four genotypes revealed 4609 genotype-specific orthogroups, which might contribute in defining genotype-specific phenotypes. The transcriptome data mining led to the identification of unigenes involved in the isoprenogenesis. The homology searches, combined with the phylogenetic and expression analyses provided information about candidate genes concerning the biosynthesis of monoterpene aldehyde, monoterpene alcohol, and monoterpene esters. In addition, the present study suggests a potential role of geranial reductase like enzyme in the biosynthesis of monoterpene aldehyde in Cymbopogon spp. The detailed analysis of the candidate pathway genes suggested that multiple enzymatic routes might be involved in the biosynthesis of aromatic monoterpenes in the genus Cymbopogon. The present study provides deeper insights into the biosynthesis of monoterpenes, which will be useful for the genetic improvement of these aromatic grasses.

Cuscuta campestris fine-tunes gene expression during haustoriogenesis as an adaptation to different hosts

The Cuscuta genus comprises obligate parasitic plants that have an unusually wide host range. Whether Cuscuta uses different infection strategies for different hosts or whether the infection strategy is mechanistically and enzymatically conserved remains unknown. To address this, we investigated molecular events during the interaction between field dodder (Cuscuta campestris) and two host species of the Solanum genus that are known to react differently to parasitic infection. We found that host gene induction, particularly of cell wall fortifying genes, coincided with a differential induction of genes for cell wall degradation in the parasite in the cultivated tomato (Solanum lycopersicum) but not in a wild relative (Solanum pennellii). This indicates that the parasite can adjust its gene expression in response to its host. This idea was supported by the increased expression of C. campestris genes encoding an endo-β-1,4-mannanase in response to exposure of the parasite to purified mono- and polysaccharides in a host-independent infection system. Our results suggest multiple key roles of the host cell wall in determining the outcome of an infection attempt.

Automatic annotation of the bHLH gene family in plants

BACKGROUND

The bHLH transcription factor family is named after the basic helix-loop-helix (bHLH) domain that is a characteristic element of their members. Understanding the function and characteristics of this family is important for the examination of a wide range of functions. As the availability of genome sequences and transcriptome assemblies has increased significantly, the need for automated solutions that provide reliable functional annotations is emphasised.

RESULTS

A phylogenetic approach was adapted for the automatic identification and functional annotation of the bHLH transcription factor family. The bHLH_annotator, designed for the automated functional annotation of bHLHs, was implemented in Python3. Sequences of bHLHs described in literature were collected to represent the full diversity of bHLH sequences. Previously described orthologs form the basis for the functional annotation assignment to candidates which are also screened for bHLH-specific motifs. The pipeline was successfully deployed on the two Arabidopsis thaliana accessions Col-0 and Nd-1, the monocot species Dioscorea dumetorum, and a transcriptome assembly of Croton tiglium. Depending on the applied search parameters for the initial candidates in the pipeline, species-specific candidates or members of the bHLH family which experienced domain loss can be identified.

CONCLUSIONS

The bHLH_annotator allows a detailed and systematic investigation of the bHLH family in land plant species and classifies candidates based on bHLH-specific characteristics, which distinguishes the pipeline from other established functional annotation tools. This provides the basis for the functional annotation of the bHLH family in land plants and the systematic examination of a wide range of functions regulated by this transcription factor family.

Identification of novel sources of partial and incomplete hypersensitive resistance to rust and associated genomic regions in common bean

Common bean (Phaseolus vulgaris) is one of the legume crops most consumed worldwide and bean rust is one of the most severe foliar biotrophic fungal diseases impacting its production. In this work, we searched for new sources of rust resistance (Uromyces appendiculatus) in a representative collection of the Portuguese germplasm, known to have accessions with an admixed genetic background between Mesoamerican and Andean gene pools. We identified six accessions with incomplete hypersensitive resistance and 20 partially resistant accessions of Andean, Mesoamerican, and admixed origin. We detected 11 disease severity-associated single-nucleotide polymorphisms (SNPs) using a genome-wide association approach. Six of the associations were related to partial (incomplete non-hypersensitive) resistance and five to incomplete hypersensitive resistance, and the proportion of variance explained by each association varied from 4.7 to 25.2%. Bean rust severity values ranged from 0.2 to 49.1% and all the infection types were identified, reflecting the diversity of resistance mechanisms deployed by the Portuguese germplasm.The associations with U. appendiculatus partial resistance were located in chromosome Pv08, and with incomplete hypersensitive resistance in chromosomes Pv06, Pv07, and Pv08, suggesting an oligogenic inheritance of both types of resistance. A resolution to the gene level was achieved for eight of the associations. The candidate genes proposed included several resistance-associated enzymes, namely β-amylase 7, acyl-CoA thioesterase, protein kinase, and aspartyl protease. Both SNPs and candidate genes here identified constitute promising genomics targets to develop functional molecular tools to support bean rust resistance precision breeding.

Data literacy in genome research

With an ever increasing amount of research data available, it becomes constantly more important to possess data literacy skills to benefit from this valuable resource. An integrative course was developed to teach students the fundamentals of data literacy through an engaging genome sequencing project. Each cohort of students performed planning of the experiment, DNA extraction, nanopore sequencing, genome sequence assembly, prediction of genes in the assembled sequence, and assignment of functional annotation terms to predicted genes. Students learned how to communicate science through writing a protocol in the form of a scientific paper, providing comments during a peer-review process, and presenting their findings as part of an international symposium. Many students enjoyed the opportunity to own a project and to work towards a meaningful objective.

Ontologies for increasing the FAIRness of plant research data

The importance of improving the FAIRness (findability, accessibility, interoperability, reusability) of research data is undeniable, especially in the face of large, complex datasets currently being produced by omics technologies. Facilitating the integration of a dataset with other types of data increases the likelihood of reuse, and the potential of answering novel research questions. Ontologies are a useful tool for semantically tagging datasets as adding relevant metadata increases the understanding of how data was produced and increases its interoperability. Ontologies provide concepts for a particular domain as well as the relationships between concepts. By tagging data with ontology terms, data becomes both human- and machine- interpretable, allowing for increased reuse and interoperability. However, the task of identifying ontologies relevant to a particular research domain or technology is challenging, especially within the diverse realm of fundamental plant research. In this review, we outline the ontologies most relevant to the fundamental plant sciences and how they can be used to annotate data related to plant-specific experiments within metadata frameworks, such as Investigation-Study-Assay (ISA). We also outline repositories and platforms most useful for identifying applicable ontologies or finding ontology terms.

KIPEs3: Automatic annotation of biosynthesis pathways

Flavonoids and carotenoids are pigments involved in stress mitigation and numerous other processes. Both pigment classes can contribute to flower and fruit coloration. Flavonoid aglycones and carotenoids are produced by a pathway that is largely conserved across land plants. Glycosylations, acylations, and methylations of the flavonoid aglycones can be species-specific and lead to a plethora of biochemically diverse flavonoids. We previously developed KIPEs for the automatic annotation of biosynthesis pathways and presented an application on the flavonoid aglycone biosynthesis. KIPEs3 is an improved version with additional features and the potential to identify not just the core biosynthesis players, but also candidates involved in the decoration steps and in the transport of flavonoids. Functionality of KIPEs3 is demonstrated through the analysis of the flavonoid biosynthesis in Arabidopsis thaliana Nd-1, Capsella grandiflora, and Dioscorea dumetorum. We demonstrate the applicability of KIPEs to other pathways by adding the carotenoid biosynthesis to the repertoire. As a technical proof of concept, the carotenoid biosynthesis was analyzed in the same species and Daucus carota. KIPEs3 is available as an online service to enable access without prior bioinformatics experience. KIPEs3 facilitates the automatic annotation and analysis of biosynthesis pathways with a consistent and high quality in a large number of plant species. Numerous genome sequencing projects are generating a huge amount of data sets that can be analyzed to identify evolutionary patterns and promising candidate genes for biotechnological and breeding applications.

Proteomic characterization of a lutein-hyperaccumulating Chlamydomonas reinhardtii mutant reveals photoprotection-related factors as targets for increasing cellular carotenoid content

BACKGROUND

Microalgae are emerging hosts for the sustainable production of lutein, a high-value carotenoid; however, to be commercially competitive with existing systems, their capacity for lutein sequestration must be augmented. Previous attempts to boost microalgal lutein production have focussed on upregulating carotenoid biosynthetic enzymes, in part due to a lack of metabolic engineering targets for expanding lutein storage.

RESULTS

Here, we isolated a lutein hyper-producing mutant of the model green microalga Chlamydomonas reinhardtii and characterized the metabolic mechanisms driving its enhanced lutein accumulation using label-free quantitative proteomics. Norflurazon- and high light-resistant C. reinhardtii mutants were screened to yield four mutant lines that produced significantly more lutein per cell compared to the CC-125 parental strain. Mutant 5 (Mut-5) exhibited a 5.4-fold increase in lutein content per cell, which to our knowledge is the highest fold increase of lutein in C. reinhardtii resulting from mutagenesis or metabolic engineering so far. Comparative proteomics of Mut-5 against its parental strain CC-125 revealed an increased abundance of light-harvesting complex-like proteins involved in photoprotection, among differences in pigment biosynthesis, central carbon metabolism, and translation. Further characterization of Mut-5 under varying light conditions revealed constitutive overexpression of the photoprotective proteins light-harvesting complex stress-related 1 (LHCSR1) and LHCSR3 and PSII subunit S regardless of light intensity, and increased accrual of total chlorophyll and carotenoids as light intensity increased. Although the photosynthetic efficiency of Mut-5 was comparatively lower than CC-125, the amplitude of non-photochemical quenching responses of Mut-5 was 4.5-fold higher than in CC-125 at low irradiance.

CONCLUSIONS

We used C. reinhardtii as a model green alga and identified light-harvesting complex-like proteins (among others) as potential metabolic engineering targets to enhance lutein accumulation in microalgae. These have the added value of imparting resistance to high light, although partially compromising photosynthetic efficiency. Further genetic characterization and engineering of Mut-5 could lead to the discovery of unknown players in photoprotective mechanisms and the development of a potent microalgal lutein production system.

Role of miRNAs in sucrose stress response, reactive oxygen species, and anthocyanin biosynthesis in Arabidopsis thaliana

In plants, sucrose is the main transported disaccharide that is the primary product of photosynthesis and controls a multitude of aspects of the plant life cycle including structure, growth, development, and stress response. Sucrose is a signaling molecule facilitating various stress adaptations by crosstalk with other hormones, but the molecular mechanisms are not well understood. Accumulation of high sucrose concentrations is a hallmark of many abiotic and biotic stresses, resulting in the accumulation of reactive oxygen species and secondary metabolite anthocyanins that have antioxidant properties. Previous studies have shown that several MYeloBlastosis family/MYB transcription factors are positive and negative regulators of sucrose-induced anthocyanin accumulation and subject to microRNA (miRNA)-mediated post-transcriptional silencing, consistent with the notion that miRNAs may be "nodes" in crosstalk signaling by virtue of their sequence-guided targeting of different homologous family members. In this study, we endeavored to uncover by deep sequencing small RNA and mRNA transcriptomes the effects of exogenous high sucrose stress on miRNA abundances and their validated target transcripts in Arabidopsis. We focused on genotype-by-treatment effects of high sucrose stress in Production of Anthocyanin Pigment 1-Dominant/pap1-D, an activation-tagged dominant allele of MYB75 transcription factor, a positive effector of secondary metabolite anthocyanin pathway. In the process, we discovered links to reactive oxygen species signaling through miR158/161/173-targeted Pentatrico Peptide Repeat genes and two novel non-canonical targets of high sucrose-induced miR408 and miR398b*(star), relevant to carbon metabolic fluxes: Flavonoid 3'-Hydroxlase (F3'H), an important enzyme in determining the B-ring hydroxylation pattern of flavonoids, and ORANGE a post-translational regulator of Phytoene Synthase expression, respectively. Taken together, our results contribute to understanding the molecular mechanisms of carbon flux shifts from primary to secondary metabolites in response to high sugar stress.

Functional analysis of Salix purpurea genes support roles for ARR17 and GATA15 as master regulators of sex determination

The Salicaceae family is of growing interest in the study of dioecy in plants because the sex determination region (SDR) has been shown to be highly dynamic, with differing locations and heterogametic systems between species. Without the ability to transform and regenerate Salix in tissue culture, previous studies investigating the mechanisms regulating sex in the genus Salix have been limited to genome resequencing and differential gene expression, which are mostly descriptive in nature, and functional validation of candidate sex determination genes has not yet been conducted. Here, we used Arabidopsis to functionally characterize a suite of previously identified candidate genes involved in sex determination and sex dimorphism in the bioenergy shrub willow Salix purpurea. Six candidate master regulator genes for sex determination were heterologously expressed in Arabidopsis, followed by floral proteome analysis. In addition, 11 transcription factors with predicted roles in mediating sex dimorphism downstream of the SDR were tested using DAP-Seq in both male and female S. purpurea DNA. The results of this study provide further evidence to support models for the roles of ARR17 and GATA15 as master regulator genes of sex determination in S. purpurea, contributing to a regulatory system that is notably different from that of its sister genus Populus. Evidence was also obtained for the roles of two transcription factors, an AP2/ERF family gene and a homeodomain-like transcription factor, in downstream regulation of sex dimorphism.

The Brassica napus boron deficient inflorescence transcriptome resembles a wounding and infection response

Oilseed rape and other crops of Brassica napus have a high demand for boron (B). Boron deficiencies result in the inhibition of root growth, and eventually premature flower abortion. Understanding the genetic mechanisms underlying flower abortion in B-limiting conditions could provide the basis to enhance B-efficiency and prevent B-deficiency-related yield losses. In this study, we assessed transcriptomic responses to B-deficiency in diverse inflorescence tissues at multiple time points of soil-grown plants that were phenotypically unaffected by B-deficiency until early flowering. Whilst transcript levels of known B transporters were higher in B-deficient samples, these remained remarkably stable as the duration of B-deficiency increased. Meanwhile, GO-term enrichment analysis indicated a growing response resembling that of a pathogen or pest attack, escalating to a huge transcriptome response in shoot heads at mid-flowering. Grouping differentially expressed genes within this tissue into MapMan functional bins indicated enrichment of genes related to wounding, jasmonic acid and WRKY transcription factors. Individual candidate genes for controlling the "flowering-without-seed-setting" phenotype from within MapMan biotic stress bins include those of the metacaspase family, which have been implicated in orchestrating programmed cell death. Overall temporal expression patterns observed here imply a dynamic response to B-deficiency, first increasing expression of B transporters before recruiting various biotic stress-related pathways to coordinate targeted cell death, likely in response to as yet unidentified B-deficiency induced damage-associated molecular patterns (DAMPs). This response indicates new pathways to target and dissect to control B-deficiency-induced flower abortion and to develop more B-efficient crops.

PLANTdataHUB: a collaborative platform for continuous FAIR data sharing in plant research

In modern reproducible, hypothesis-driven plant research, scientists are increasingly relying on research data management (RDM) services and infrastructures to streamline the processes of collecting, processing, sharing, and archiving research data. FAIR (i.e., findable, accessible, interoperable, and reusable) research data play a pivotal role in enabling the integration of interdisciplinary knowledge and facilitating the comparison and synthesis of a wide range of analytical findings. The PLANTdataHUB offers a solution that realizes RDM of scientific (meta)data as evolving collections of files in a directory - yielding FAIR digital objects called ARCs - with tools that enable scientists to plan, communicate, collaborate, publish, and reuse data on the same platform while gaining continuous quality control insights. The centralized platform is scalable from personal use to global communities and provides advanced federation capabilities for institutions that prefer to host their own satellite instances. This approach borrows many concepts from software development and adapts them to fit the challenges of the field of modern plant science undergoing digital transformation. The PLANTdataHUB supports researchers in each stage of a scientific project with adaptable continuous quality control insights, from the early planning phase to data publication. The central live instance of PLANTdataHUB is accessible at (https://git.nfdi4plants.org), and it will continue to evolve as a community-driven and dynamic resource that serves the needs of contemporary plant science.

The Effects of Rhizophagus irregularis Inoculation on Transcriptome of Medicago lupulina Leaves at Early Vegetative and Flowering Stages of Plant Development

The study is aimed at revealing the effects of Rhizophagus irregularis inoculation on the transcriptome of Medicago lupulina leaves at the early (second leaf formation) and later (flowering) stages of plant development. A pot experiment was conducted under conditions of low phosphorus (P) level in the substrate. M. lupulina plants were characterized by high mycorrhizal growth response and mycorrhization parameters. Library sequencing was performed on the Illumina HiseqXTen platform. Significant changes in the expression of 4863 (padj < 0.01) genes from 34049 functionally annotated genes were shown by Massive Analysis of cDNA Ends (MACE-Seq). GO enrichment analysis using the Kolmogorov-Smirnov test was performed, and 244 functional GO groups were identified, including genes contributing to the development of effective AM symbiosis. The Mercator online tool was used to assign functional classes of differentially expressed genes (DEGs). The early stage was characterized by the presence of six functional classes that included only upregulated GO groups, such as genes of carbohydrate metabolism, cellular respiration, nutrient uptake, photosynthesis, protein biosynthesis, and solute transport. At the later stage (flowering), the number of stimulated GO groups was reduced to photosynthesis and protein biosynthesis. All DEGs of the GO:0016036 group were downregulated because AM plants had higher resistance to phosphate starvation. For the first time, the upregulation of genes encoding thioredoxin in AM plant leaves was shown. It was supposed to reduce ROS level and thus, consequently, enhance the mechanisms of antioxidant protection in M. lupulina plants under conditions of low phosphorus level. Taken together, the obtained results indicate genes that are the most important for the effective symbiosis with M. lupulina and might be engaged in other plant species.

Global transcriptome analysis reveals fungal disease responsive core gene regulatory landscape in tea

Fungal infections are the inevitable limiting factor for productivity of tea. Transcriptome reprogramming recruits multiple regulatory pathways during pathogen infection. A comprehensive meta-analysis was performed utilizing previously reported, well-replicated transcriptomic datasets from seven fungal diseases of tea. The study identified a cumulative set of 18,517 differentially expressed genes (DEGs) in tea, implicated in several functional clusters, including the MAPK signaling pathway, transcriptional regulation, and the biosynthesis of phenylpropanoids. Gene set enrichment analyses under each pathogen stress elucidated that DEGs were involved in ethylene metabolism, secondary metabolism, receptor kinase activity, and various reactive oxygen species detoxification enzyme activities. Expressional fold change of combined datasets highlighting 2258 meta-DEGs shared a common transcriptomic response upon fungal stress in tea. Pervasive duplication events caused biotic stress-responsive core DEGs to appear in multiple copies throughout the tea genome. The co-expression network of meta-DEGs in multiple modules demonstrated the coordination of appropriate pathways, most of which involved cell wall organization. The functional coordination was controlled by a number of hub genes and miRNAs, leading to pathogenic resistance or susceptibility. This first-of-its-kind meta-analysis of host-pathogen interaction generated consensus candidate loci as molecular signatures, which can be associated with future resistance breeding programs in tea.

Differential metabolic reprogramming in developing soybean embryos in response to nutritional conditions and abscisic acid

Seed storage compound deposition is influenced by both maternal and filial tissues. Within this framework, we analyzed strategies that operate during the development and filling of soybean embryos, using in vitro culture systems combined with metabolomics and proteomics approaches. The carbon:nitrogen ratio (C:N) of the maternal supply and the hormone abscisic acid (ABA) are specific and interacting signals inducing differential metabolic reprogrammings linked to changes in the accumulation of storage macromolecules like proteins or oils. Differences in the abundance of sugars, amino acids, enzymes, transporters, transcription factors, and proteins involved in signaling were detected. Embryos adapted to the nutritional status by enhancing the metabolism of both carbon and nitrogen under lower C:N ratio condition or only carbon under higher C:N ratio condition. ABA turned off multiple pathways especially in high availability of amino acids, prioritizing the storage compounds biosynthesis. Common responses induced by ABA involved increased sucrose uptake (to increase the sink force) and oleosin (oil body structural component) accumulation. In turn, ABA differentially promoted protein degradation under lower nitrogen supply in order to sustain the metabolic demands. Further, the operation of a citrate shuttle was suggested by transcript quantification and enzymatic activity measurements. The results obtained are useful to help define biotechnological tools and technological approaches to improve oil and protein yields, with direct impact on human and animal nutrition as well as in green chemistry.

Breeding effects on durum wheat traits detected using GWAS and haplotype block analysis

Introduction

The recent boosting of genomic data in durum wheat (Triticum turgidum subsp. durum) offers the opportunity to better understand the effects of breeding on the genetic structures that regulate the expression of traits of agronomic interest. Furthermore, the identification of DNA markers useful for marker-assisted selection could also improve the reliability of technical protocols used for variety protection and registration.

Methods

Within this motivation context, 123 durum wheat accessions, classified into three groups: landraces (LR), ancient (OC) and modern cultivars (MC), were evaluated in two locations, for 34 agronomic traits, including UPOV descriptors, to assess the impact of changes that occurred during modern breeding.

Results

The association mapping analysis, performed with 4,241 SNP markers and six multi-locus-GWAS models, revealed 28 reliable Quantitative Trait Nucleotides (QTNs) related to plant morphology and kernel-related traits. Some important genes controlling flowering time and plant height were in linkage disequilibrium (LD) decay with QTNs identified in this study. A strong association for yellow berry was found on chromosome 6A (Q.Yb-6A) in a region containing the nadh-ubiquinone oxidoreductase subunit, a gene involved in starch metabolism. The Q.Kcp-2A harbored the PPO locus, with the associated marker (Ku_c13700_1196) in LD decay with Ppo-A1 and Ppo-A2. Interestingly, the Q.FGSGls-2B.1, identified by RAC875_c34512_685 for flag leaf glaucosity, mapped less than 1 Mb from the Epistatic inhibitors of glaucousness (Iw1), thus representing a good candidate for supporting the morphological DUS traits also with molecular markers. LD haplotype block approach revealed a higher diversity, richness and length of haploblocks in MC than OC and LR (580 in LR, 585 in OC and 612 in MC), suggesting a possible effect exerted by breeding programs on genomic regions associated with the agronomic traits.

Discussion

Our findings pave new ways to support the phenotypic characterization necessary for variety registration by using a panel of cost-effectiveness SNP markers associated also to the UPOV descriptors. Moreover, the panel of associated SNPs might represent a reservoir of favourable alleles to use in durum wheat breeding and genetics.

The role of microRNAs in responses to drought and heat stress in peanut (Arachis hypogaea)

MicroRNAs (miRNAs) are 21-24 nt small RNAs (sRNAs) that negatively regulate protein-coding genes and/or trigger phased small-interfering RNA (phasiRNA) production. Two thousand nine hundred miRNA families, of which ∼40 are deeply conserved, have been identified in ∼80 different plant species genomes. miRNA functions in response to abiotic stresses is less understood than their roles in development. Only seven peanut MIRNA families are documented in miRBase, yet a reference genome assembly is now published and over 480 plant-like MIRNA loci were predicted in the diploid peanut progenitor Arachis duranensis genome. We explored by computational analysis of a leaf sRNA library and publicly available sRNA, degradome, and transcriptome datasets the miRNA and phasiRNA space associated with drought and heat stresses in peanut. We characterized 33 novel candidate and 33 ancient conserved families of MIRNAs and present degradome evidence for their cleavage activities on mRNA targets, including several noncanonical targets and novel phasiRNA-producing noncoding and mRNA loci with validated novel targets such as miR1509 targeting serine/threonine-protein phosphatase7 and miRc20 and ahy-miR3514 targeting penta-tricopeptide repeats (PPRs), in contradistinction to other claims of miR1509/173/7122 superfamily miRNAs indirectly targeting PPRs via TAS-like noncoding RNA loci. We characterized the inverse correlations of significantly differentially expressed drought- and heat-regulated miRNAs, assayed by sRNA blots or transcriptome datasets, with target mRNA expressions in the same datasets. Meta-analysis of an expression atlas and over representation of miRNA target genes in co-expression networks suggest that miRNAs have functions in unique aspects of peanut gynophore development. Genome-wide MIRNA annotation of the published allopolyploid peanut genome can facilitate molecular breeding of value-added traits.