PLAZA 4.0: an integrative resource for functional, evolutionary and comparative plant genomics

Comprehensive analysis of Alfin-like transcription factors associated with drought and salt stresses in wheat (Triticum aestivum L.)

BACKGROUND

Alfin-like proteins are a kind of plant-specific transcription factors, and play vital roles in plant growth, development and stress responses.

RESULTS

In this study, a total of 27 Alfin-like transcription factors were identified in wheat. TaAL genes were unevenly distributed on chromosome. Phylogenetic analysis showed TaAL genes were divided into AL-B and AL-C subfamilies, and TaALs with closer evolutionary relationships generally shared more similar exon-intron structures and conserved motifs. The cis-acting element analysis showed MBS, ABRE and CGTCA-motif were the most common in TaAL promoters. The interacting proteins and downstream target genes of TaAL genes were also investigated in wheat. The transcriptome data and real-time PCR results indicated TaAL genes were differentially expressed under drought and salt stresses, and TaAL1-B was significantly up-regulated in response to drought stress. In addition, association analysis revealed that TaAL1-B-Hap-I allelic variation had significantly higher survival rate compared to TaAL1-B-Hap-II under drought stress.

CONCLUSIONS

These results will provide vital information to increase our understanding of the Alfin-like gene family in wheat, and help us in breeding better wheat varieties in the future.

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.

Genome resources for three modern cotton lines guide future breeding efforts

Cotton (Gossypium hirsutum L.) is the key renewable fibre crop worldwide, yet its yield and fibre quality show high variability due to genotype-specific traits and complex interactions among cultivars, management practices and environmental factors. Modern breeding practices may limit future yield gains due to a narrow founding gene pool. Precision breeding and biotechnological approaches offer potential solutions, contingent on accurate cultivar-specific data. Here we address this need by generating high-quality reference genomes for three modern cotton cultivars ('UGA230', 'UA48' and 'CSX8308') and updating the 'TM-1' cotton genetic standard reference. Despite hypothesized genetic uniformity, considerable sequence and structural variation was observed among the four genomes, which overlap with ancient and ongoing genomic introgressions from 'Pima' cotton, gene regulatory mechanisms and phenotypic trait divergence. Differentially expressed genes across fibre development correlate with fibre production, potentially contributing to the distinctive fibre quality traits observed in modern cotton cultivars. These genomes and comparative analyses provide a valuable foundation for future genetic endeavours to enhance global cotton yield and sustainability.

Transcriptome Analysis of Stigmas of Vicia faba L. Flowers

Pollination in angiosperms depends on complex communication between pollen grains and stigmas, classified as wet or dry, depending on the presence or absence of secretions at the stigma surface, respectively. In species with wet stigma, the cuticle is disrupted and the presence of exudates is indicative of their receptivity. Most stigma studies are focused on a few species and families, many of them with self-incompatibility systems. However, there is scarce knowledge about the stigma composition in Fabaceae, the third angiosperm family, whose stigmas have been classified as semidry. Here we report the first transcriptome profiling and DEGs of Vicia faba L. styles and stigmas from autofertile (flowers able to self-fertilize in the absence of manipulation, whose exudate is released spontaneously) and autosterile (flowers that need to be manipulated to break the cuticle and release the exudates to be receptive) inbred lines. From the 76,269 contigs obtained from the de novo assembly, only 45.1% of the sequences were annotated with at least one GO term. A total of 115,920, 75,489, and 70,801 annotations were assigned to Biological Process (BP), Cellular Component (CC), and Molecular Function (MF) categories, respectively, and 5918 differentially expressed genes (DEGs) were identified between the autofertile and the autosterile lines. Among the most enriched metabolic pathways in the DEGs subset were those related with amino acid biosynthesis, terpenoid metabolism, or signal transduction. Some DEGs have been related with previous QTLs identified for autofertility traits, and their putative functions are discussed. The results derived from this work provide an important transcriptomic reference for style-stigma processes to aid our understanding of the molecular mechanisms involved in faba bean fertilization.

Exploring the potential role of four Rhizophagus irregularis nuclear effectors: opportunities and technical limitations

Arbuscular mycorrhizal fungi (AMF) are obligate symbionts that interact with the roots of most land plants. The genome of the AMF model species Rhizophagus irregularis contains hundreds of predicted small effector proteins that are secreted extracellularly but also into the plant cells to suppress plant immunity and modify plant physiology to establish a niche for growth. Here, we investigated the role of four nuclear-localized putative effectors, i.e., GLOIN707, GLOIN781, GLOIN261, and RiSP749, in mycorrhization and plant growth. We initially intended to execute the functional studies in Solanum lycopersicum, a host plant of economic interest not previously used for AMF effector biology, but extended our studies to the model host Medicago truncatula as well as the non-host Arabidopsis thaliana because of the technical advantages of working with these models. Furthermore, for three effectors, the implementation of reverse genetic tools, yeast two-hybrid screening and whole-genome transcriptome analysis revealed potential host plant nuclear targets and the downstream triggered transcriptional responses. We identified and validated a host protein interactors participating in mycorrhization in the host.S. lycopersicum and demonstrated by transcriptomics the effectors possible involvement in different molecular processes, i.e., the regulation of DNA replication, methylglyoxal detoxification, and RNA splicing. We conclude that R. irregularis nuclear-localized effector proteins may act on different pathways to modulate symbiosis and plant physiology and discuss the pros and cons of the tools used.

NKS1/ELMO4 is an integral protein of a pectin synthesis protein complex and maintains Golgi morphology and cell adhesion in Arabidopsis

Adjacent plant cells are connected by specialized cell wall regions, called middle lamellae, which influence critical agricultural characteristics, including fruit ripening and organ abscission. Middle lamellae are enriched in pectin polysaccharides, specifically homogalacturonan (HG). Here, we identify a plant-specific Arabidopsis DUF1068 protein, called NKS1/ELMO4, that is required for middle lamellae integrity and cell adhesion. NKS1 localizes to the Golgi apparatus and loss of NKS1 results in changes to Golgi structure and function. The nks1 mutants also display HG deficient phenotypes, including reduced seedling growth, changes to cell wall composition, and tissue integrity defects. These phenotypes are comparable to qua1 and qua2 mutants, which are defective in HG biosynthesis. Notably, genetic interactions indicate that NKS1 and the QUAs work in a common pathway. Protein interaction analyses and modeling corroborate that they work together in a stable protein complex with other pectin-related proteins. We propose that NKS1 is an integral part of a large pectin synthesis protein complex and that proper function of this complex is important to support Golgi structure and function.

D6PK plasma membrane polarity requires a repeated CXX(X)P motif and PDK1-dependent phosphorylation

D6 PROTEIN KINASE (D6PK) is a polarly localized plasma-membrane-associated kinase from Arabidopsis thaliana that activates polarly distributed PIN-FORMED auxin transporters. D6PK moves rapidly to and from the plasma membrane, independent of its PIN-FORMED targets. The middle D6PK domain, an insertion between kinase subdomains VII and VIII, is required and sufficient for association and polarity of the D6PK plasma membrane. How D6PK polarity is established and maintained remains to be shown. Here we show that cysteines from repeated middle domain CXX(X)P motifs are S-acylated and required for D6PK membrane association. While D6PK S-acylation is not detectably regulated during intracellular transport, phosphorylation of adjacent serine residues, in part in dependence on the upstream 3-PHOSPHOINOSITIDE-DEPENDENT PROTEIN KINASE, promotes D6PK transport, controls D6PK residence time at the plasma membrane and prevents its lateral diffusion. We thus identify new mechanisms for the regulation of D6PK plasma membrane interaction and polarity.

Sucrose synthase gene SUS3 could enhance cold tolerance in tomato

Tomatoes are susceptible to damage from cold temperatures in all stages of growth. Therefore, it is important to identify genetic resources and genes that can enhance tomato's ability to tolerate cold. In this study, a population of 223 tomato accessions was used to identify the sensitivity or tolerance of plants to cold stress. Transcriptome analysis of these accessions revealed that SUS3, a member of the sucrose synthase gene family, was induced by cold stress. We further investigated the role of SUS3 in cold stress by overexpression (OE) and RNA interference (RNAi). Compared with the wild type, SUS3-OE lines accumulated less MDA and electrolyte leakage and more proline and soluble sugar, maintained higher activities of SOD and CAT, reduced superoxide radicals, and suffered less membrane damage under cold. Thus, our findings indicate that SUS3 plays a crucial role in the response to cold stress. This study indicates that SUS3 may serve as a direct target for genetic engineering and improvement projects, which aim to augment the cold tolerance of tomato crops.

Phosphoproteome analyses pinpoint the F-box protein SLOW MOTION as a regulator of warm temperature-mediated hypocotyl growth in Arabidopsis

Hypocotyl elongation is controlled by several signals and is a major characteristic of plants growing in darkness or under warm temperature. While already several molecular mechanisms associated with this process are known, protein degradation and associated E3 ligases have hardly been studied in the context of warm temperature. In a time-course phosphoproteome analysis on Arabidopsis seedlings exposed to control or warm ambient temperature, we observed reduced levels of diverse proteins over time, which could be due to transcription, translation, and/or degradation. In addition, we observed differential phosphorylation of the LRR F-box protein SLOMO MOTION (SLOMO) at two serine residues. We demonstrate that SLOMO is a negative regulator of hypocotyl growth, also under warm temperature conditions, and protein-protein interaction studies revealed possible interactors of SLOMO, such as MKK5, DWF1, and NCED4. We identified DWF1 as a likely SLOMO substrate and a regulator of warm temperature-mediated hypocotyl growth. We propose that warm temperature-mediated regulation of SLOMO activity controls the abundance of hypocotyl growth regulators, such as DWF1, through ubiquitin-mediated degradation.

MINI-AC: inference of plant gene regulatory networks using bulk or single-cell accessible chromatin profiles

Gene regulatory networks (GRNs) represent the interactions between transcription factors (TF) and their target genes. Plant GRNs control transcriptional programs involved in growth, development, and stress responses, ultimately affecting diverse agricultural traits. While recent developments in accessible chromatin (AC) profiling technologies make it possible to identify context-specific regulatory DNA, learning the underlying GRNs remains a major challenge. We developed MINI-AC (Motif-Informed Network Inference based on Accessible Chromatin), a method that combines AC data from bulk or single-cell experiments with TF binding site (TFBS) information to learn GRNs in plants. We benchmarked MINI-AC using bulk AC datasets from different Arabidopsis thaliana tissues and showed that it outperforms other methods to identify correct TFBS. In maize, a crop with a complex genome and abundant distal AC regions, MINI-AC successfully inferred leaf GRNs with experimentally confirmed, both proximal and distal, TF-target gene interactions. Furthermore, we showed that both AC regions and footprints are valid alternatives to infer AC-based GRNs with MINI-AC. Finally, we combined MINI-AC predictions from bulk and single-cell AC datasets to identify general and cell-type specific maize leaf regulators. Focusing on C4 metabolism, we identified diverse regulatory interactions in specialized cell types for this photosynthetic pathway. MINI-AC represents a powerful tool for inferring accurate AC-derived GRNs in plants and identifying known and novel candidate regulators, improving our understanding of gene regulation in plants.

Evolution of a plant growth-regulatory protein interaction specificity

Specific protein-protein interactions (PPIs) enable biological regulation. However, the evolution of PPI specificity is little understood. Here we trace the evolution of the land-plant growth-regulatory DELLA-SLY1/GID2 PPI, revealing progressive increase in specificity of affinity of SLY1/GID2 for a particular DELLA form. While early-diverging SLY1s display relatively broad-range DELLA affinity, later-diverging SLY1s tend towards increasingly stringent affinity for a specific DELLA A' form generated by the growth-promoting phytohormone gibberellin (GA). Our novel mutational strategy reveals amino acid substitutions contributing to the evolution of Arabidopsis thaliana SLY1 A' specificity, also showing that routes permitting reversion to broader affinity became increasingly constrained over evolutionary time. We suggest that progressive affinity narrowing may be an important evolutionary driver of PPI specificity and that increase in SLY1/GID2-DELLA specificity enabled the enhanced flexibility of plant physiological environmental adaptation conferred by the GA-DELLA growth-regulatory mechanism.

Evolutionarily conserved protein motifs drive interactions between the plant nucleoskeleton and nuclear pores

The nucleoskeleton forms a filamentous meshwork under the nuclear envelope and contributes to the regulation of nuclear shape and gene expression. To understand how the Arabidopsis (Arabidopsis thaliana) nucleoskeleton physically connects to the nuclear periphery in plants, we investigated the Arabidopsis nucleoskeleton protein KAKU4 and sought for functional regions responsible for its localization at the nuclear periphery. We identified 3 conserved peptide motifs within the N-terminal region of KAKU4, which are required for intermolecular interactions of KAKU4 with itself, interaction with the nucleoskeleton protein CROWDED NUCLEI (CRWN), localization at the nuclear periphery, and nuclear elongation in differentiated tissues. Unexpectedly, we find these motifs to be present also in NUP82 and NUP136, 2 plant-specific nucleoporins from the nuclear pore basket. We further show that NUP82, NUP136, and KAKU4 have a common evolutionary history predating nonvascular land plants with KAKU4 mainly localizing outside the nuclear pore suggesting its divergence from an ancient nucleoporin into a new nucleoskeleton component. Finally, we demonstrate that both NUP82 and NUP136, through their shared N-terminal motifs, interact with CRWN and KAKU4 proteins revealing the existence of a physical continuum between the nuclear pore and the nucleoskeleton in plants.

Interspecies co-expression analysis of lateral root development using inducible systems in rice, Medicago, and Arabidopsis

Lateral roots are crucial for plant growth and development, making them an important target for research aiming to improve crop yields and food security. However, their endogenous ontogeny and, as it were, stochastic appearance challenge their study. Lateral Root Inducible Systems (LRIS) can be used to overcome these challenges by inducing lateral roots massively and synchronously. The combination of LRISs with transcriptomic approaches significantly advanced our insights in the molecular control of lateral root formation, in particular for Arabidopsis. Despite this success, LRISs have been underutilized for other plant species or for lateral root developmental stages later than the initiation. In this study, we developed and/or adapted LRISs in rice, Medicago, and Arabidopsis to perform RNA-sequencing during time courses that cover different developmental stages of lateral root formation and primordium development. As such, our study provides three extensive datasets of gene expression profiles during lateral root development in three different plant species. The three LRISs are highly effective but timing and spatial distribution of lateral root induction vary among the species. Detailed characterization of the stages in time and space in the respective species enabled an interspecies co-expression analysis to identify conserved players involved in lateral root development, as illustrated for the AUX/IAA and LBD gene families. Overall, our results provide a valuable resource to identify potentially conserved regulatory mechanisms in lateral root development, and as such will contribute to a better understanding of the complex regulatory network underlying lateral root development.

Dicer-like 3a mediates intergenerational resistance against root-knot nematodes in rice via hormone responses

In a continuously changing and challenging environment, passing down the memory of encountered stress factors to offspring could provide an evolutionary advantage. In this study, we demonstrate the existence of "intergenerational acquired resistance" in the progeny of rice (Oryza sativa) plants attacked by the belowground parasitic nematode Meloidogyne graminicola. Transcriptome analyses revealed that genes involved in defense pathways are generally downregulated in progeny of nematode-infected plants under uninfected conditions but show a stronger induction upon nematode infection. This phenomenon was termed "spring loading" and depends on initial downregulation by the 24-nucleotide (nt) siRNA biogenesis gene dicer-like 3a (dcl3a) involved in the RNA-directed DNA methylation pathway. Knockdown of dcl3a led to increased nematode susceptibility and abolished intergenerational acquired resistance, as well as jasmonic acid/ethylene spring loading in the offspring of infected plants. The importance of ethylene signaling in intergenerational resistance was confirmed by experiments on a knockdown line of ethylene insensitive 2 (ein2b), which lacks intergenerational acquired resistance. Taken together, these data indicate a role for DCL3a in regulating plant defense pathways during both within-generation and intergenerational resistance against nematodes in rice.

Complex origin, evolution, and diversification of non-canonically organized OVATE-OFP and OVATE-Like OFP gene pair across Embryophyta

Ovate Family Proteins (OFP) is a plant-specific gene family of negative transcriptional regulators. Till-date, a handful of in-silico studies have provided glimpses into family size, expansion patterns, and genic features across all major plant lineages. A major lacuna exists in understanding origin of organisation complexity of members such as those arranged in a head-to-head manner which may lead to transcriptional co-regulation via a common bi-directional promoter. To address this gap, we investigated the origin, organization and evolution of two head-to-head arranged gene pairs of homologs of AtOFP2-AtOFP17, and, AtOFP4-AtOFP20 across Archaeplastida. The ancestral forms of AtOFP2, AtOFP4, AtOFP17, and AtOFP20 are likely to have evolved in last common ancestors of Embryophyta (land plants) given their complete absence in Rhodophyta and Chlorophyta. The OFP gene family originated and expanded in Bryophyta, including protein variants with complete (OVATE-OFP) or partial (OVATE-Like OFP) OVATE domain; with head-to-head organization present only in Spermatophyta (gymnosperms and angiosperms). Ancestral State Reconstruction revealed the origin of head-to-head organized gene pair in gymnosperms, with both genes being OVATE-OFP (homologs of AtOFP2/4). Phylogenetic reconstruction and copy number analysis suggests the presence of a single copy of the head-to-head arranged pair of OFP2/4 (OVATE)-OFP17/20 (OVATE-Like) in all angiosperms except Brassicaceae, and a duplication event in last common ancestor of core Brassicaceae approximately 32-54 MYA leading to origin of AtOFP2-AtOFP17 and AtOFP4-AtOFP20 as paralogs. Synteny analysis of genomic regions harbouring homologs of AtOFP2-AtOFP17, AtOFP4-AtOFP20 and AtOFP2/4-AtOFP17/20 across angiosperms suggested ancestral nature of AtOFP2-AtOFP17 gene pair. The present study thus establishes the orthology and evolutionary history of two non-canonically organised gene pairs with variation in their OVATE domain. The non-canonical organisation, atleast in Brassicaceae, has the potential of generating complex transcriptional regulation mediated via a common bi-directional promoter. The study thus lays down a framework to understand evolution of gene and protein structure, transcriptional regulation and function across a phylogenetic lineage through comparative analyses.

Spatial transcriptomics of a lycophyte root sheds light on root evolution

Plant roots originated independently in lycophytes and euphyllophytes, whereas early vascular plants were rootless. The organization of the root apical meristem in euphyllophytes is well documented, especially in the model plant Arabidopsis. However, little is known about lycophyte roots and their molecular innovations during evolution. In this study, spatial transcriptomics was used to detect 97 root-related genes in the roots of the lycophyte Selaginella moellendorffii. A high number of genes showed expression patterns similar to what has been reported for seed plants, supporting the idea of a highly convergent evolution of mechanisms to control root development. Interaction and complementation data of SHORTROOT (SHR) and SCARECROW (SCR) homologs, furthermore, support a comparable regulation of the ground tissue (GT) between euphyllophytes and lycophytes. Root cap formation, in contrast, appears to be differently regulated. Several experiments indicated an important role of the WUSCHEL-RELATED HOMEOBOX13 gene SmWOX13a in Selaginella root cap formation. In contrast to multiple Arabidopsis WOX paralogs, SmWOX13a is able to induce root cap cells in Arabidopsis and has functionally conserved homologs in the fern Ceratopteris richardii. Lycophytes and a part of the euphyllophytes, therefore, may share a common mechanism regulating root cap formation, which was diversified or lost during seed plant evolution. In summary, we here provide a new spatial data resource for the Selaginella root, which in general advocates for conserved mechanisms to regulate root development but shows a clear divergence in the control of root cap formation, with a novel putative role of WOX genes in root cap formation in non-seed plants.

Evolution of the Tóxicos en Levadura 63 (TL63) gene family in plants and functional characterization of Arabidopsis thaliana TL63 under oxidative stress

MAIN CONCLUSION

TL63 orthologs were angiosperm specific and had undergone motifs loss and gain, and increased purifying selection. AtTL63 was involved in the response of yeast and Arabidopsis plants to oxidative stress. The Tóxicos en Levadura (TL) family, a class of E3 ubiquitin ligases with typical RING-H2 type zinc finger structure, plays a pivotal role in mediating physiological processes and responding to stress in plants. However, the evolution and function of TL63 remain unclear. In this study, TL63 homologs were dated roughly back to the origin of land plants and confirmed to have subjected to the gain and loss of motifs and increased purifying selection. Phylogenetic analysis displayed that 279 TL63s could be divided into four main clades (Clade A-D). Notably, the ancestral tandem TL40/41 cluster contributed to the expansion of modern Brassicaceae TL40/41. The substitution rate tests revealed that the TL63 lineage was evidently different from other lineages. The codon usage index exhibited that monocotyledons preferred to use not A3s and T3s, but C3s, G3s, CAI, CBI and Fop. Sequence analysis showed that the TL63 homologs had conserved TM and GLD motifs and RING-H2 domain whose key amino acid residues accounted for the high average abundance. Particularly, Arabidopsis thaliana TL63 (AtTL63) was located in the nuclei, cell membranes and peroxisomes and expressed universally and significantly throughout A. thaliana development. Under H2O2 treatment, low or moderate expression of the AtTL63 held beneficial effects on the growth and viability of yeast cells and the mutation or overexpression of the AtTL63 positively affected the growth of A. thaliana plants. In brief, this study could supply useful insight into the evolution of the plant TL63s and the AtTL63 functions under oxidative stress.

An mTRAN-mRNA interaction mediates mitochondrial translation initiation in plants

Plant mitochondria represent the largest group of respiring organelles on the planet. Plant mitochondrial messenger RNAs (mRNAs) lack Shine-Dalgarno-like ribosome-binding sites, so it is unknown how plant mitoribosomes recognize mRNA. We show that "mitochondrial translation factors" mTRAN1 and mTRAN2 are land plant-specific proteins, required for normal mitochondrial respiration chain biogenesis. Our studies suggest that mTRANs are noncanonical pentatricopeptide repeat (PPR)-like RNA binding proteins of the mitoribosomal "small" subunit. We identified conserved Adenosine (A)/Uridine (U)-rich motifs in the 5' regions of plant mitochondrial mRNAs. mTRAN1 binds this motif, suggesting that it is a mitoribosome homing factor to identify mRNAs. We demonstrate that mTRANs are likely required for translation of all plant mitochondrial mRNAs. Plant mitochondrial translation initiation thus appears to use a protein-mRNA interaction that is divergent from bacteria or mammalian mitochondria.

Functional annotation of proteins for signaling network inference in non-model species

Molecular biology aims to understand cellular responses and regulatory dynamics in complex biological systems. However, these studies remain challenging in non-model species due to poor functional annotation of regulatory proteins. To overcome this limitation, we develop a multi-layer neural network that determines protein functionality directly from the protein sequence. We annotate kinases and phosphatases in Glycine max. We use the functional annotations from our neural network, Bayesian inference principles, and high resolution phosphoproteomics to infer phosphorylation signaling cascades in soybean exposed to cold, and identify Glyma.10G173000 (TOI5) and Glyma.19G007300 (TOT3) as key temperature regulators. Importantly, the signaling cascade inference does not rely upon known kinase motifs or interaction data, enabling de novo identification of kinase-substrate interactions. Conclusively, our neural network shows generalization and scalability, as such we extend our predictions to Oryza sativa, Zea mays, Sorghum bicolor, and Triticum aestivum. Taken together, we develop a signaling inference approach for non-model species leveraging our predicted kinases and phosphatases.

CropGF: a comprehensive visual platform for crop gene family mining and analysis

A gene family refers to a group of genes that share a common ancestry and encode proteins or RNA molecules with similar functions or structural features. Gene families play a crucial role in determining the traits of plants and can be utilized to develop new crop varieties. Therefore, a comprehensive database of gene family is significant for gaining deep insight into crops. To address this need, we have developed CropGF (https://bis.zju.edu.cn/cropgf), a comprehensive visual platform that encompasses six important crops (rice, wheat, maize, barley, sorghum and foxtail millet) and one model plant (Arabidopsis), as well as genomics, transcriptomics and proteomics data for gene family mining and analysis, covering a total of 314 611 genes and 4399 types of domains. CropGF provides a versatile search system that allows for the identification of gene families and their members in a single crop or multiple crops. Users can customize their search based on gene family domains and/or homology using keywords or BLAST. To enhance usability, we have collected the corresponding ID information from various public databases for both genes and domains. Furthermore, CropGF comprises numerous downstream analysis modules, such as ka/ks analysis, phylogenetic tree construction, subcellular localization analysis and more. These visually-displayed modules provide intuitive insights into gene expression patterns, gene family expansion and functional relationships across different molecular levels and different species. We believe that CropGF will be a valuable resource for deep mining and analysis in future studies of crop gene families. Database URL https://bis.zju.edu.cn/cropgf.

The saponin bomb: a nucleolar-localized β-glucosidase hydrolyzes triterpene saponins in Medicago truncatula

Plants often protect themselves from their own bioactive defense metabolites by storing them in less active forms. Consequently, plants also need systems allowing correct spatiotemporal reactivation of such metabolites, for instance under pathogen or herbivore attack. Via co-expression analysis with public transcriptomes, we determined that the model legume Medicago truncatula has evolved a two-component system composed of a β-glucosidase, denominated G1, and triterpene saponins, which are physically separated from each other in intact cells. G1 expression is root-specific, stress-inducible, and coregulated with that of the genes encoding the triterpene saponin biosynthetic enzymes. However, the G1 protein is stored in the nucleolus and is released and united with its typically vacuolar-stored substrates only upon tissue damage, partly mediated by the surfactant action of the saponins themselves. Subsequently, enzymatic removal of carbohydrate groups from the saponins creates a pool of metabolites with an increased broad-spectrum antimicrobial activity. The evolution of this defense system benefited from both the intrinsic condensation abilities of the enzyme and the bioactivity properties of its substrates. We dub this two-component system the saponin bomb, in analogy with the mustard oil and cyanide bombs, commonly used to describe the renowned β-glucosidase-dependent defense systems for glucosinolates and cyanogenic glucosides.

Identification and characterization of CYP71 subclade cytochrome P450 enzymes involved in the biosynthesis of bitterness compounds in Cichorium intybus

Industrial chicory (Cichorium intybus var. sativum) and witloof (C. intybus var. foliosum) are crops with an important economic value, mainly cultivated for inulin production and as a leafy vegetable, respectively. Both crops are rich in nutritionally relevant specialized metabolites with beneficial effects for human health. However, their bitter taste, caused by the sesquiterpene lactones (SLs) produced in leaves and taproot, limits wider applications in the food industry. Changing the bitterness would thus create new opportunities with a great economic impact. Known genes encoding enzymes involved in the SL biosynthetic pathway are GERMACRENE A SYNTHASE (GAS), GERMACRENE A OXIDASE (GAO), COSTUNOLIDE SYNTHASE (COS) and KAUNIOLIDE SYNTHASE (KLS). In this study, we integrated genome and transcriptome mining to further unravel SL biosynthesis. We found that C. intybus SL biosynthesis is controlled by the phytohormone methyl jasmonate (MeJA). Gene family annotation and MeJA inducibility enabled the pinpointing of candidate genes related with the SL biosynthetic pathway. We specifically focused on members of subclade CYP71 of the cytochrome P450 family. We verified the biochemical activity of 14 C. intybus CYP71 enzymes transiently produced in Nicotiana benthamiana and identified several functional paralogs for each of the GAO, COS and KLS genes, pointing to redundancy in and robustness of the SL biosynthetic pathway. Gene functionality was further analyzed using CRISPR/Cas9 genome editing in C. intybus. Metabolite profiling of mutant C. intybus lines demonstrated a successful reduction in SL metabolite production. Together, this study increases our insights into the C. intybus SL biosynthetic pathway and paves the way for the engineering of C. intybus bitterness.

Functional and comparative analysis of THI1 gene in grasses with a focus on sugarcane

De novo synthesis of thiamine (vitamin B1) in plants depends on the action of thiamine thiazole synthase, which synthesizes the thiazole ring, and is encoded by the THI1 gene. Here, we investigated the evolution and diversity of THI1 in Poaceae, where C4 and C3 photosynthetic plants co-evolved. An ancestral duplication of THI1 is observed in Panicoideae that remains in many modern monocots, including sugarcane. In addition to the two sugarcane copies (ScTHI1-1 and ScTHI1-2), we identified ScTHI1-2 alleles showing differences in their sequence, indicating divergence between ScTHI1-2a and ScTHI1-2b. Such variations are observed only in the Saccharum complex, corroborating the phylogeny. At least five THI1 genomic environments were found in Poaceae, two in sugarcane, M. sinensis, and S. bicolor. The THI1 promoter in Poaceae is highly conserved at 300 bp upstream of the start codon ATG and has cis-regulatory elements that putatively bind to transcription factors associated with development, growth, development and biological rhythms. An experiment set to compare gene expression levels in different tissues across the sugarcane R570 life cycle showed that ScTHI1-1 was expressed mainly in leaves regardless of age. Furthermore, ScTHI1 displayed relatively high expression levels in meristem and culm, which varied with the plant age. Finally, yeast complementation studies with THI4-defective strain demonstrate that only ScTHI1-1 and ScTHI1-2b isoforms can partially restore thiamine auxotrophy, albeit at a low frequency. Taken together, the present work supports the existence of multiple origins of THI1 harboring genomic regions in Poaceae with predicted functional redundancy. In addition, it questions the contribution of the levels of the thiazole ring in C4 photosynthetic plant tissues or potentially the relevance of the THI1 protein activity.

6mA DNA Methylation on Genes in Plants Is Associated with Gene Complexity, Expression and Duplication

N6-methyladenine (6mA) DNA methylation has emerged as an important epigenetic modification in eukaryotes. Nevertheless, the evolution of the 6mA methylation of homologous genes after species and after gene duplications remains unclear in plants. To understand the evolution of 6mA methylation, we detected the genome-wide 6mA methylation patterns of four lotus plants (Nelumbo nucifera) from different geographic origins by nanopore sequencing and compared them to patterns in Arabidopsis and rice. Within lotus, the genomic distributions of 6mA sites are different from the widely studied 5mC methylation sites. Consistently, in lotus, Arabidopsis and rice, 6mA sites are enriched around transcriptional start sites, positively correlated with gene expression levels, and preferentially retained in highly and broadly expressed orthologs with longer gene lengths and more exons. Among different duplicate genes, 6mA methylation is significantly more enriched and conserved in whole-genome duplicates than in local duplicates. Overall, our study reveals the convergent patterns of 6mA methylation evolution based on both lineage and duplicate gene divergence, which underpin their potential role in gene regulatory evolution in plants.

Multispecies comparison of host responses to Fusarium circinatum challenge in tropical pines show consistency in resistance mechanisms

Fusarium circinatum poses a threat to both commercial and natural pine forests. Large variation in host resistance exists between species, with many economically important species being susceptible. Development of resistant genotypes could be expedited and optimised by investigating the molecular mechanisms underlying host resistance and susceptibility as well as increasing the available genetic resources. RNA-seq data, from F. circinatum inoculated and mock-inoculated ca. 6-month-old shoot tissue at 3- and 7-days postinoculation, was generated for three commercially important tropical pines, Pinus oocarpa, Pinus maximinoi and Pinus greggii. De novo transcriptomes were assembled and used to investigate the NLR and PR gene content within available pine references. Host responses to F. circinatum challenge were investigated in P. oocarpa (resistant) and P. greggii (susceptible), in comparison to previously generated expression profiles from Pinus tecunumanii (resistant) and Pinus patula (susceptible). Expression results indicated crosstalk between induced salicylate, jasmonate and ethylene signalling is involved in host resistance and compromised in susceptible hosts. Additionally, higher constitutive expression of sulfur metabolism and flavonoid biosynthesis in resistant hosts suggest involvement of these metabolites in resistance.

Puncta-localized TRAF domain protein TC1b contributes to the autoimmunity of snc1

Immune receptors play important roles in the perception of pathogens and initiation of immune responses in both plants and animals. Intracellular nucleotide-binding domain leucine-rich repeat (NLR)-type receptors constitute a major class of receptors in vascular plants. In the Arabidopsis thaliana mutant suppressor of npr1-1, constitutive 1 (snc1), a gain-of-function mutation in the NLR gene SNC1 leads to SNC1 overaccumulation and constitutive activation of defense responses. From a CRISPR/Cas9-based reverse genetics screen in the snc1 autoimmune background, we identified that mutations in TRAF CANDIDATE 1b (TC1b), a gene encoding a protein with four tumor necrosis factor receptor-associated factor (TRAF) domains, can suppress snc1 phenotypes. TC1b does not appear to be a general immune regulator as it is not required for defense mediated by other tested immune receptors. TC1b also does not physically associate with SNC1, affect SNC1 accumulation, or affect signaling of the downstream helper NLRs represented by ACTIVATED DISEASE RESISTANCE PROTEIN 1-L2 (ADR1-L2), suggesting that TC1b impacts snc1 autoimmunity in a unique way. TC1b can form oligomers and localizes to punctate structures of unknown function. The puncta localization of TC1b strictly requires its coiled-coil (CC) domain, whereas the functionality of TC1b requires the four TRAF domains in addition to the CC. Overall, we uncovered the TRAF domain protein TC1b as a novel positive contributor to plant immunity.

SMAP design: a multiplex PCR amplicon and gRNA design tool to screen for natural and CRISPR-induced genetic variation

Multiplex amplicon sequencing is a versatile method to identify genetic variation in natural or mutagenized populations through eco-tilling or multiplex CRISPR screens. Such genotyping screens require reliable and specific primer designs, combined with simultaneous gRNA design for CRISPR screens. Unfortunately, current tools are unable to combine multiplex gRNA and primer design in a high-throughput and easy-to-use manner with high design flexibility. Here, we report the development of a bioinformatics tool called SMAP design to overcome these limitations. We tested SMAP design on several plant and non-plant genomes and obtained designs for more than 80-90% of the target genes, depending on the genome and gene family. We validated the designs with Illumina multiplex amplicon sequencing and Sanger sequencing in Arabidopsis, soybean, and maize. We also used SMAP design to perform eco-tilling by tilling PCR amplicons across nine candidate genes putatively associated with haploid induction in Cichorium intybus. We screened 60 accessions of chicory and witloof and identified thirteen knockout haplotypes and their carriers. SMAP design is an easy-to-use command-line tool that generates highly specific gRNA and/or primer designs for any number of loci for CRISPR or natural variation screens and is compatible with other SMAP modules for seamless downstream analysis.

RAV1 family members function as transcriptional regulators and play a positive role in plant disease resistance

Phytopathogens pose a severe threat to agriculture and strengthening the plant defense response is an important strategy for disease control. Here, we report that AtRAV1, an AP2 and B3 domain-containing transcription factor, is required for basal plant defense in Arabidopsis thaliana. The atrav1 mutant lines demonstrate hyper-susceptibility against fungal pathogens (Rhizoctonia solani and Botrytis cinerea), whereas AtRAV1 overexpressing lines exhibit disease resistance against them. Enhanced expression of various defense genes and activation of mitogen-activated protein kinases (AtMPK3 and AtMPK6) are observed in the R. solani infected overexpressing lines, but not in the atrav1 mutant plants. An in vitro phosphorylation assay suggests AtRAV1 to be a novel phosphorylation target of AtMPK3. Bimolecular fluorescence complementation and yeast two-hybrid assays support physical interactions between AtRAV1 and AtMPK3. Overexpression of the native as well as phospho-mimic but not the phospho-defective variant of AtRAV1 imparts disease resistance in the atrav1 mutant A. thaliana lines. On the other hand, overexpression of AtRAV1 fails to impart disease resistance in the atmpk3 mutant. These analyses emphasize that AtMPK3-mediated phosphorylation of AtRAV1 is important for the elaboration of the defense response in A. thaliana. Considering that RAV1 homologs are conserved in diverse plant species, we propose that they can be gainfully deployed to impart disease resistance in agriculturally important crop plants. Indeed, overexpression of SlRAV1 (a member of the RAV1 family) imparts disease tolerance against not only fungal (R. solani and B. cinerea), but also against bacterial (Ralstonia solanacearum) pathogens in tomato, whereas silencing of the gene enhances disease susceptibility.

Plant lineage-specific PIKMIN1 drives APC/CCCS52A2 E3-ligase activity-dependent cell division

The anaphase-promoting complex/cyclosome (APC/C) marks key cell cycle proteins for proteasomal breakdown, thereby ensuring unidirectional progression through the cell cycle. Its target recognition is temporally regulated by activating subunits, one of which is called CELL CYCLE SWITCH 52 A2 (CCS52A2). We sought to expand the knowledge on the APC/C by using the severe growth phenotypes of CCS52A2-deficient Arabidopsis (Arabidopsis thaliana) plants as a readout in a suppressor mutagenesis screen, resulting in the identification of the previously undescribed gene called PIKMIN1 (PKN1). PKN1 deficiency rescues the disorganized root stem cell phenotype of the ccs52a2-1 mutant, whereas an excess of PKN1 inhibits the growth of ccs52a2-1 plants, indicating the need for control of PKN1 abundance for proper development. Accordingly, the lack of PKN1 in a wild-type background negatively impacts cell division, while its systemic overexpression promotes proliferation. PKN1 shows a cell cycle phase-dependent accumulation pattern, localizing to microtubular structures, including the preprophase band, the mitotic spindle, and the phragmoplast. PKN1 is conserved throughout the plant kingdom, with its function in cell division being evolutionarily conserved in the liverwort Marchantia polymorpha. Our data thus demonstrate that PKN1 represents a novel, plant-specific protein with a role in cell division that is likely proteolytically controlled by the CCS52A2-activated APC/C.

Identification of subfunctionalized aggregate-remodeling J-domain proteins in Arabidopsis thaliana

J-domain proteins (JDPs) are critical components of the cellular protein quality control machinery, playing crucial roles in preventing the formation and, solubilization of cytotoxic protein aggregates. Bacteria, yeast, and plants additionally have large, multimeric heat shock protein 100 (Hsp100)-class disaggregases that resolubilize protein aggregates. JDPs interact with aggregated proteins and specify the aggregate-remodeling activities of Hsp70s and Hsp100s. However, the aggregate-remodeling properties of plant JDPs are not well understood. Here we identify eight orthologs of Sis1 (an evolutionarily conserved Class II JDP of budding yeast) in Arabidopsis thaliana with distinct aggregate-remodeling functionalities. Six of these JDPs associate with heat-induced protein aggregates in vivo and co-localize with Hsp101 at heat-induced protein aggregate centers. Consistent with a role in solubilizing cytotoxic protein aggregates, an atDjB3 mutant had defects in both solubilizing heat-induced aggregates and acquired thermotolerance as compared with wild-type seedlings. Next, we used yeast prions as protein aggregate models to show that the six JDPs have distinct aggregate-remodeling properties. Results presented in this study, as well as findings from phylogenetic analysis, demonstrate that plants harbor multiple, evolutionarily conserved JDPs with capacity to process a variety of protein aggregate conformers induced by heat and other stressors.

A Survey of the Transcriptomic Resources in Durum Wheat: Stress Responses, Data Integration and Exploitation

Durum wheat (Triticum turgidum subsp. durum (Desf.) Husn.) is an allotetraploid cereal crop of worldwide importance, given its use for making pasta, couscous, and bulgur. Under climate change scenarios, abiotic (e.g., high and low temperatures, salinity, drought) and biotic (mainly exemplified by fungal pathogens) stresses represent a significant limit for durum cultivation because they can severely affect yield and grain quality. The advent of next-generation sequencing technologies has brought a huge development in transcriptomic resources with many relevant datasets now available for durum wheat, at various anatomical levels, also focusing on phenological phases and environmental conditions. In this review, we cover all the transcriptomic resources generated on durum wheat to date and focus on the corresponding scientific insights gained into abiotic and biotic stress responses. We describe relevant databases, tools and approaches, including connections with other "omics" that could assist data integration for candidate gene discovery for bio-agronomical traits. The biological knowledge summarized here will ultimately help in accelerating durum wheat breeding.

Extensive sequence duplication in Arabidopsis revealed by pseudo-heterozygosity

BACKGROUND

It is apparent that genomes harbor much structural variation that is largely undetected for technical reasons. Such variation can cause artifacts when short-read sequencing data are mapped to a reference genome. Spurious SNPs may result from mapping of reads to unrecognized duplicated regions. Calling SNP using the raw reads of the 1001 Arabidopsis Genomes Project we identified 3.3 million (44%) heterozygous SNPs. Given that Arabidopsis thaliana (A. thaliana) is highly selfing, and that extensively heterozygous individuals have been removed, we hypothesize that these SNPs reflected cryptic copy number variation.

RESULTS

The heterozygosity we observe consists of particular SNPs being heterozygous across individuals in a manner that strongly suggests it reflects shared segregating duplications rather than random tracts of residual heterozygosity due to occasional outcrossing. Focusing on such pseudo-heterozygosity in annotated genes, we use genome-wide association to map the position of the duplicates. We identify 2500 putatively duplicated genes and validate them using de novo genome assemblies from six lines. Specific examples included an annotated gene and nearby transposon that transpose together. We also demonstrate that cryptic structural variation produces highly inaccurate estimates of DNA methylation polymorphism.

CONCLUSIONS

Our study confirms that most heterozygous SNP calls in A. thaliana are artifacts and suggest that great caution is needed when analyzing SNP data from short-read sequencing. The finding that 10% of annotated genes exhibit copy-number variation, and the realization that neither gene- nor transposon-annotation necessarily tells us what is actually mobile in the genome suggests that future analyses based on independently assembled genomes will be very informative.

Genome-wide analysis of bromodomain gene family in Arabidopsis and rice

The bromodomain-containing proteins (BRD-proteins) belongs to family of 'epigenetic mark readers', integral to epigenetic regulation. The BRD-members contain a conserved 'bromodomain' (BRD/BRD-fold: interacts with acetylated-lysine in histones), and several additional domains, making them structurally/functionally diverse. Like animals, plants also contain multiple Brd-homologs, however the extent of their diversity and impact of molecular events (genomic duplications, alternative splicing, AS) therein, is relatively less explored. The present genome-wide analysis of Brd-gene families of Arabidopsis thaliana and Oryza sativa showed extensive diversity in structure of genes/proteins, regulatory elements, expression pattern, domains/motifs, and the bromodomain (w.r.t. length, sequence, location) among the Brd-members. Orthology analysis identified thirteen ortholog groups (OGs), three paralog groups (PGs) and four singleton members (STs). While more than 40% Brd-genes were affected by genomic duplication events in both plants, AS-events affected 60% A. thaliana and 41% O. sativa genes. These molecular events affected various regions (promoters, untranslated regions, exons) of different Brd-members with potential impact on expression and/or structure-function characteristics. RNA-Seq data analysis indicated differences in tissue-specificity and stress response of Brd-members. Analysis by RT-qPCR revealed differential abundance and salt stress response of duplicate A. thaliana and O. sativa Brd-genes. Further analysis of AtBrd gene, AtBrdPG1b showed salinity-induced modulation of splicing pattern. Bromodomain (BRD)-region based phylogenetic analysis placed the A. thaliana and O. sativa homologs into clusters/sub-clusters, mostly consistent with ortholog/paralog groups. The bromodomain-region displayed several conserved signatures in key BRD-fold elements (α-helices, loops), along with variations (1-20 sites) and indels among the BRD-duplicates. Homology modeling and superposition identified structural variations in BRD-folds of divergent and duplicate BRD-members, which might affect their interaction with the chromatin histones, and associated functions. The study also showed contribution of various duplication events in Brd-gene family expansion among diverse plants, including several monocot and dicot plant species.

A chromosome-level genome assembly for Erianthus fulvus provides insights into its biofuel potential and facilitates breeding for improvement of sugarcane

Erianthus produces substantial biomass, exhibits a good Brix value, and shows wide environmental adaptability, making it a potential biofuel plant. In contrast to closely related sorghum and sugarcane, Erianthus can grow in degraded soils, thus releasing pressure on agricultural lands used for biofuel production. However, the lack of genomic resources for Erianthus hinders its genetic improvement, thus limiting its potential for biofuel production. In the present study, we generated a chromosome-scale reference genome for Erianthus fulvus Nees. The genome size estimated by flow cytometry was 937 Mb, and the assembled genome size was 902 Mb, covering 96.26% of the estimated genome size. A total of 35 065 protein-coding genes were predicted, and 67.89% of the genome was found to be repetitive. A recent whole-genome duplication occurred approximately 74.10 million years ago in the E. fulvus genome. Phylogenetic analysis showed that E. fulvus is evolutionarily closer to S. spontaneum and diverged after S. bicolor. Three of the 10 chromosomes of E. fulvus formed through rearrangements of ancestral chromosomes. Phylogenetic reconstruction of the Saccharum complex revealed a polyphyletic origin of the complex and a sister relationship of E. fulvus with Saccharum sp., excluding S. arundinaceum. On the basis of the four amino acid residues that provide substrate specificity, the E. fulvus SWEET proteins were classified as mono- and disaccharide sugar transporters. Ortho-QTL genes identified for 10 biofuel-related traits may aid in the rapid screening of E. fulvus populations to enhance breeding programs for improved biofuel production. The results of this study provide valuable insights for breeding programs aimed at improving biofuel production in E. fulvus and enhancing sugarcane introgression programs.

Translational research from models to crops: comparative genomics for plant breeding

The concept of translational research, which originated in the medical field in the 1980s, consists in improving the efficient transfer of research results obtained in a species (which can be considered as a model or pivot) to all the species for which these results are of interest for its improvement in Agriculture. In this context, comparative genomics is an important tool for translational research, effectively identifying genes controlling common functions between species. Editing and phenotyping tools must thus allow the functional validation of the gene conserved within the species for which the knowledge has been extrapolated, that is to say transferred, and the identification of the best alleles and associated genotypes for exploitation in current breeding programs.

A long-read and short-read transcriptomics approach provides the first high-quality reference transcriptome and genome annotation for Pseudotsuga menziesii (Douglas-fir)

Douglas-fir (Pseudotsuga menziesii) is native to western North America. It grows in a wide range of environmental conditions and is an important timber tree. Although there are several studies on the gene expression responses of Douglas-fir to abiotic cues, the absence of high-quality transcriptome and genome data is a barrier to further investigation. Like for most conifers, the available transcriptome and genome reference dataset for Douglas-fir remains fragmented and requires refinement. We aimed to generate a highly accurate, and complete reference transcriptome and genome annotation. We deep-sequenced the transcriptome of Douglas-fir needles from seedlings that were grown under nonstress control conditions or a combination of heat and drought stress conditions using long-read (LR) and short-read (SR) sequencing platforms. We used 2 computational approaches, namely de novo and genome-guided LR transcriptome assembly. Using the LR de novo assembly, we identified 1.3X more high-quality transcripts, 1.85X more "complete" genes, and 2.7X more functionally annotated genes compared to the genome-guided assembly approach. We predicted 666 long noncoding RNAs and 12,778 unique protein-coding transcripts including 2,016 putative transcription factors. We leveraged the LR de novo assembled transcriptome with paired-end SR and a published single-end SR transcriptome to generate an improved genome annotation. This was conducted with BRAKER2 and refined based on functional annotation, repetitive content, and transcriptome alignment. This high-quality genome annotation has 51,419 unique gene models derived from 322,631 initial predictions. Overall, our informatics approach provides a new reference Douglas-fir transcriptome assembly and genome annotation with considerably improved completeness and functional annotation.

Evolution of the WRKY66 Gene Family and Its Mutations Generated by the CRISPR/Cas9 System Increase the Sensitivity to Salt Stress in Arabidopsis

Group Ⅲ WRKY transcription factors (TFs) play pivotal roles in responding to the diverse abiotic stress and secondary metabolism of plants. However, the evolution and function of WRKY66 remains unclear. Here, WRKY66 homologs were traced back to the origin of terrestrial plants and found to have been subjected to both motifs' gain and loss, and purifying selection. A phylogenetic analysis showed that 145 WRKY66 genes could be divided into three main clades (Clade A-C). The substitution rate tests indicated that the WRKY66 lineage was significantly different from others. A sequence analysis displayed that the WRKY66 homologs had conserved WRKY and C2HC motifs with higher proportions of crucial amino acid residues in the average abundance. The AtWRKY66 is a nuclear protein, salt- and ABA- inducible transcription activator. Simultaneously, under salt stress and ABA treatments, the superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT) activities, as well as the seed germination rates of Atwrky66-knockdown plants generated by the clustered, regularly interspaced, short palindromic repeats/CRISPR-associated 9 (CRISPR/Cas9) system, were all lower than those of wild type (WT) plants, but the relative electrolyte leakage (REL) was higher, indicating the increased sensitivities of the knockdown plants to the salt stress and ABA treatments. Moreover, RNA-seq and qRT-PCR analyses revealed that several regulatory genes in the ABA-mediated signaling pathway involved in stress response of the knockdown plants were significantly regulated, being evidenced by the more moderate expressions of the genes. Therefore, the AtWRKY66 likely acts as a positive regulator in the salt stress response, which may be involved in an ABA-mediated signaling pathway.

Auxin-Responsive (Phospho)proteome Analysis Reveals Key Biological Processes and Signaling Associated with Shoot-Borne Crown Root Development in Rice

The rice root system is primarily composed of shoot-borne adventitious/crown roots (ARs/CRs) that develop from the coleoptile base, and therefore, it is an excellent model system for studying shoot-to-root trans-differentiation process. We reveal global changes in protein and metabolite abundance and protein phosphorylation in response to an auxin stimulus during CR development. The liquid chromatography-tandem mass spectrometry (LC-MS/MS) and gas chromatography-mass spectrometry (GC-MS) analyses of developing crown root primordia (CRP) and emerged CRs identified 334 proteins and 12 amino acids, respectively, that were differentially regulated upon auxin treatment. Gene ontology enrichment analysis of global proteome data uncovered the biological processes associated with chromatin conformational change, gene expression and cell cycle that were regulated by auxin signaling. Spatial gene expression pattern analysis of differentially abundant proteins disclosed their stage-specific dynamic expression pattern during CRP development. Further, our tempo-spatial gene expression and functional analyses revealed that auxin creates a regulatory module during CRP development and activates ethylene biosynthesis exclusively during CRP initiation. Further, the phosphoproteome analysis identified 8,220 phosphosites, which could be mapped to 1,594 phosphoproteins and of which 66 phosphosites were differentially phosphorylated upon auxin treatment. Importantly, we observed differential phosphorylation of the cyclin-dependent kinase G-2 (OsCDKG;2) and cell wall proteins, in response to auxin signaling, suggesting that auxin-dependent phosphorylation may be required for cell cycle activation and cell wall synthesis during root organogenesis. Thus, our study provides evidence for the translational and post-translational regulation during CR development downstream of the auxin signaling pathway.

Rhizogenic Agrobacterium protein RolB interacts with the TOPLESS repressor proteins to reprogram plant immunity and development

Rhizogenic Agrobacterium strains comprise biotrophic pathogens that cause hairy root disease (HRD) on hydroponically grown Solanaceae and Cucurbitaceae crops, besides being widely explored agents for the creation of hairy root cultures for the sustainable production of plant-specialized metabolites. Hairy root formation is mediated through the expression of genes encoded on the T-DNA of the root-inducing (Ri) plasmid, of which several, including root oncogenic locus B (rolB), play a major role in hairy root development. Despite decades of research, the exact molecular function of the proteins encoded by the rol genes remains enigmatic. Here, by means of TurboID-mediated proximity labeling in tomato (Solanum lycopersicum) hairy roots, we identified the repressor proteins TOPLESS (TPL) and Novel Interactor of JAZ (NINJA) as direct interactors of RolB. Although these interactions allow RolB to act as a transcriptional repressor, our data hint at another in planta function of the RolB oncoprotein. Hence, by a series of plant bioassays, transcriptomic and DNA-binding site enrichment analyses, we conclude that RolB can mitigate the TPL functioning so that it leads to a specific and partial reprogramming of phytohormone signaling, immunity, growth, and developmental processes. Our data support a model in which RolB manipulates host transcription, at least in part, through interaction with TPL, to facilitate hairy root development. Thereby, we provide important mechanistic insights into this renowned oncoprotein in HRD.

Genome-Wide Analysis of Lipoxygenase (LOX) Genes in Angiosperms

Lipoxygenases (LOXs) are enzymes that catalyze the addition of an oxygen molecule to unsaturated fatty acids, thus forming hydroperoxides. In plants, these enzymes are encoded by a multigene family found in several organs with varying activity patterns, by which they are classified as LOX9 or LOX13. They are involved in several physiological functions, such as growth, fruit development, and plant defense. Despite several studies on genes of the LOX family in plants, most studies are restricted to a single species or a few closely related species. This study aimed to analyze the diversity, evolution, and expression of LOX genes in angiosperm species. We identified 247 LOX genes among 23 species of angiosperms and basal plants. Phylogenetic analyses identified clades supporting LOX13 and two main clades for LOX9: LOX9_A and LOX9_B. Eudicot species such as Tarenaya hassleriana, Capsella rubella, and Arabidopsis thaliana did not present LOX9_B genes; however, LOX9_B was present in all monocots used in this study. We identified that there were potential new subcellular localization patterns and conserved residues of oxidation for LOX9 and LOX13 yet unexplored. In summary, our study provides a basis for the further functional and evolutionary study of lipoxygenases in angiosperms.

The origin and evolution of salicylic acid signaling and biosynthesis in plants

Salicylic acid (SA) plays a pivotal role in plant response to biotic and abiotic stress. Several core SA signaling regulators and key proteins in SA biosynthesis have been well characterized. However, much remains unknown about the origin, evolution, and early diversification of core elements in plant SA signaling and biosynthesis. In this study, we identified 10 core protein families in SA signaling and biosynthesis across green plant lineages. We found that the key SA signaling receptors, the nonexpresser of pathogenesis-related (NPR) proteins, originated in the most recent common ancestor (MRCA) of land plants and formed divergent groups in the ancestor of seed plants. However, key transcription factors for SA signaling, TGACG motif-binding proteins (TGAs), originated in the MRCA of streptophytes, arguing for the stepwise evolution of core SA signaling in plants. Different from the assembly of the core SA signaling pathway in the ancestor of seed plants, SA exists extensively in green plants, including chlorophytes and streptophyte algae. However, the full isochorismate synthase (ICS)-based SA synthesis pathway was first assembled in the MRCA of land plants. We further revealed that the ancient abnormal inflorescence meristem 1 (AIM1)-based β-oxidation pathway is crucial for the biosynthesis of SA in chlorophyte algae, and this biosynthesis pathway may have facilitated the adaptation of early-diverging green algae to the high-light-intensity environment on land. Taken together, our findings provide significant insights into the early evolution and diversification of plant SA signaling and biosynthesis pathways, highlighting a crucial role of SA in stress tolerance during plant terrestrialization.

PEATmoss: A Gene Expression Atlas for Bryophytes

PEATmoss is an interactive gene expression atlas for bryophytes, which originally unified Physcomitrium patens RNA-seq and microarray expression data from multiple gene annotation versions. This atlas includes more than 100 experiments of P. patens, is expanding to host Anthoceros agrestis and Marchantia polymorpha, and aims to host data from more species in the future. PEATmoss has multiple visualization methods and tools for data downloading and is connected to the Physcomitrium patens Gene Model Lookup DB (PpGML DB), which links P. patens genes to annotations and resources from several databases and contains tools for gene version lookup and sequence and annotation extraction. Among the new features available in PEATmoss are dataset privacy control, multispecies menu, interactive color scale, co-expression network visualization, and replicate data downloading.

Genome-wide identification, phylogeny, and expression analysis of GRF transcription factors in pineapple (Ananas comosus)

Background

Pineapple is the only commercially grown fruit crop in the Bromeliaceae family and has significant agricultural, industrial, economic, and ornamental value. GRF (growth-regulating factor) proteins are important transcription factors that have evolved in seed plants (embryophytes). They contain two conserved domains, QLQ (Gln, Leu, Gln) and WRC (Trp, Arg, Cys), and regulate multiple aspects of plant growth and stress response, including floral organ development, leaf growth, and hormone responses. The GRF family has been characterized in a number of plant species, but little is known about this family in pineapple and other bromeliads.

Main discoveries

We identified eight GRF transcription factor genes in pineapple, and phylogenetic analysis placed them into five subfamilies (I, III, IV, V, VI). Segmental duplication appeared to be the major contributor to expansion of the AcGRF family, and the family has undergone strong purifying selection during evolution. Relative to that of other gene families, the gene structure of the GRF family showed less conservation. Analysis of promoter cis-elements suggested that AcGRF genes are widely involved in plant growth and development. Transcriptome data and qRT-PCR results showed that, with the exception of AcGRF5, the AcGRFs were preferentially expressed in the early stage of floral organ development and AcGRF2 was strongly expressed in ovules. Gibberellin treatment significantly induced AcGRF7/8 expression, suggesting that these two genes may be involved in the molecular regulatory pathway by which gibberellin promotes pineapple fruit expansion.

Conclusion

AcGRF proteins appear to play a role in the regulation of floral organ development and the response to gibberellin. The information reported here provides a foundation for further study of the functions of AcGRF genes and the traits they regulate.

Genetic control of tolerance to drought stress in soybean

BACKGROUND

Drought stress limits the production of soybean [Glycine max (L.) Merr.], which is the most grown high-value legume crop worldwide. Breeding for drought tolerance is a difficult endeavor and understanding the genetic basis of drought tolerance in soybean is therefore crucial for harnessing the genomic regions involved in the tolerance mechanisms. A genome-wide association study (GWAS) analysis was applied in a soybean germplasm collection (the EUCLEG collection) of 359 accessions relevant for breeding in Europe, to identify genomic regions and candidate genes involved in the response to short duration and long duration drought stress (SDS and LDS respectively) in soybean.

RESULTS

The phenotypic response to drought was stronger in the long duration drought (LDS) than in the short duration drought (SDS) experiment. Over the four traits considered (canopy wilting, leaf senescence, maximum absolute growth rate and maximum plant height) the variation was in the range of 8.4-25.2% in the SDS, and 14.7-29.7% in the LDS experiments. The GWAS analysis identified a total of 17 and 22 significant marker-trait associations for four traits in the SDS and LDS experiments, respectively. In the genomic regions delimited by these markers we identified a total of 12 and 16 genes with putative functions that are of particular relevance for drought stress responses including stomatal movement, root formation, photosynthesis, ABA signaling, cellular protection and cellular repair mechanisms. Some of these genomic regions co-localized with previously known QTLs for drought tolerance traits including water use efficiency, chlorophyll content and photosynthesis.

CONCLUSION

Our results indicate that the mechanism of slow wilting in the SDS might be associated with the characteristics of the root system, whereas in the LDS, slow wilting could be due to low stomatal conductance and transpiration rates enabling a high WUE. Drought-induced leaf senescence was found to be associated to ABA and ROS responses. The QTLs related to WUE contributed to growth rate and canopy height maintenance under drought stress. Co-localization of several previously known QTLs for multiple agronomic traits with the SNPs identified in this study, highlights the importance of the identified genomic regions for the improvement of agronomic performance in addition to drought tolerance in the EUCLEG collection.

syntenet: an R/Bioconductor package for the inference and analysis of synteny networks

SUMMARY

Interpreting and visualizing synteny relationships across several genomes is a challenging task. We previously proposed a network-based approach for better visualization and interpretation of large-scale microsynteny analyses. Here, we present syntenet, an R package to infer and analyze synteny networks from whole-genome protein sequence data. The package offers a simple and complete framework, including data preprocessing, synteny detection and network inference, network clustering and phylogenomic profiling, and microsynteny-based phylogeny inference. Graphical functions are also available to create publication-ready plots. Synteny networks inferred with syntenet can highlight taxon-specific gene clusters that likely contributed to the evolution of important traits, and microsynteny-based phylogenies can help resolve phylogenetic relationships under debate.

AVAILABILITY AND IMPLEMENTATION

syntenet is available on Bioconductor (https://bioconductor.org/packages/syntenet), and the source code is available on a GitHub repository (https://github.com/almeidasilvaf/syntenet).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Distinct chemical resistance-inducing stimuli result in common transcriptional, metabolic, and nematode community signatures in rice root and rhizosphere

Induced resistance (IR), a phenotypic state induced by an exogenous stimulus and characterized by enhanced resistance to future (a)biotic challenge, is an important component of plant immunity. Numerous IR-inducing stimuli have been described in various plant species, but relatively little is known about 'core' systemic responses shared by these distinct IR stimuli and the effects of IR on plant-associated microbiota. In this study, rice (Oryza sativa) leaves were treated with four distinct IR stimuli (β-aminobutyric acid, acibenzolar-S-methyl, dehydroascorbic acid, and piperonylic acid) capable of inducing systemic IR against the root-knot nematode Meloidogyne graminicola and evaluated their effect on the root transcriptome and exudome, and root-associated nematode communities. Our results reveal shared transcriptional responses-notably induction of jasmonic acid and phenylpropanoid metabolism-and shared alterations to the exudome that include increased amino acid, benzoate, and fatty acid exudation. In rice plants grown in soil from a rice field, IR stimuli significantly affected the composition of rhizosphere nematode communities 3 d after treatment, but by 14 d after treatment these changes had largely reverted. Notably, IR stimuli did not reduce nematode diversity, which suggests that IR might offer a sustainable option for managing plant-parasitic nematodes.

Orchid NAC Transcription Factors: A Focused Analysis of CUPULIFORMIS Genes

Plant transcription factors are involved in different developmental pathways. NAC transcription factors (No Apical Meristem, Arabidopsis thaliana Activating Factor, Cup-shaped Cotyledon) act in various processes, e.g., plant organ formation, response to stress, and defense mechanisms. In Antirrhinum majus, the NAC transcription factor CUPULIFORMIS (CUP) plays a role in determining organ boundaries and lip formation, and the CUP homologs of Arabidopsis and Petunia are involved in flower organ formation. Orchidaceae is one of the most species-rich families of angiosperms, known for its extraordinary diversification of flower morphology. We conducted a transcriptome and genome-wide analysis of orchid NACs, focusing on the No Apical Meristem (NAM) subfamily and CUP genes. To check whether the CUP homologs could be involved in the perianth formation of orchids, we performed an expression analysis on the flower organs of the orchid Phalaenopsis aphrodite at different developmental stages. The expression patterns of the CUP genes of P. aphrodite suggest their possible role in flower development and symmetry establishment. In addition, as observed in other species, the orchid CUP1 and CUP2 genes seem to be regulated by the microRNA, miR164. Our results represent a preliminary study of NAC transcription factors in orchids to understand the role of these genes during orchid flower formation.

Chloroplastic ascorbate peroxidases targeted to stroma or thylakoid membrane: The chicken or egg dilemma

Ascorbate peroxidases (APXs) are heme peroxidases that remove hydrogen peroxide in different subcellular compartments with concomitant ascorbate cycling. Here, we analysed and discussed phylogenetic and molecular features of the APX family. Ancient APX originated as a soluble stromal enzyme, and early during plant evolution, acquired both chloroplast-targeting and mitochondrion-targeting sequences and an alternative splicing mechanism whereby it could be expressed as a soluble or thylakoid membrane-bound enzyme. Later, independent duplication and neofunctionalization events in some angiosperm groups resulted in individual genes encoding stromal, thylakoidal and mitochondrial isoforms. These data reaffirm the complexity of plant antioxidant defenses that allow diverse plant species to acquire new means to adapt to changing environmental conditions.

Opposing effects of trans- and cis-cinnamic acid during rice coleoptile elongation

The phenylpropanoid cinnamic acid (CA) is a plant metabolite that can occur under a trans- or cis-form. In contrast to the proven bioactivity of the cis-form (c-CA), the activity of trans-CA (t-CA) is still a matter of debate. We tested both compounds using a submerged rice coleoptile assay and demonstrated that they have opposite effects on cell elongation. Notably, in the tip of rice coleoptile t-CA showed an inhibiting and c-CA a stimulating activity. By combining transcriptomics and (untargeted) metabolomics with activity assays and genetic and pharmacological experiments, we aimed to explain the underlying mechanistic processes. We propose a model in which c-CA treatment activates proton pumps and stimulates acidification of the apoplast, which in turn leads to the loosening of the cell wall, necessary for elongation. We hypothesize that c-CA also inactivates auxin efflux transporters, which might cause a local auxin accumulation in the tip of the coleoptile. For t-CA, the phenotype can partially be explained by a stimulation of cell wall polysaccharide feruloylation, leading to a more rigid cell wall. Metabolite profiling also demonstrated that salicylic acid (SA) derivatives are increased upon t-CA treatment. As SA is a known antagonist of auxin, the shift in SA homeostasis provides an additional explanation of the observed t-CA-mediated restriction on cell growth.

Cellular and gene expression patterns associated with root bifurcation in Selaginella

The roots of lycophytes branch through dichotomy or bifurcation, during which the root apex splits into two daughter roots. This is morphologically distinct from lateral root (LR) branching in the extant euphyllophytes, with LRs developing along the root axis at different distances from the apex. Although the process of root bifurcation is poorly understood, such knowledge can be important, because it may represent an evolutionarily ancient strategy that roots recruited to form new stem cells or meristems. In this study, we examined root bifurcation in the lycophyte Selaginella moellendorffii. We characterized an in vitro developmental time frame based on repetitive apex bifurcations, allowing us to sample different stages of dichotomous root branching and analyze the root meristem and root branching in S. moellendorffii at the microscopic and transcriptomic level. Our results showed that, in contrast to previous assumptions, initial cells (ICs) in the root meristem are mostly not tetrahedral but rather show an irregular shape. Tracking down the early stages of root branching argues for the occurrence of a symmetric division of the single IC, resulting in two apical stem cells that initiate root meristem bifurcation. Moreover, we generated a S. moellendorffii root branching transcriptome that resulted in the delineation of a subset of core meristem genes. The occurrence of multiple putative orthologs of meristem genes in this dataset suggests the presence of conserved pathways in the control of meristem and root stem cell establishment or maintenance.