A high resolution map of the Arabidopsis thaliana developmental transcriptome based on RNA-seq profiling

Uncovering genes involved in pollinator-driven mating system shifts and selfing syndrome evolution in Brassica rapa

Shifts in pollinator occurrence and their pollen transport effectiveness drive the evolution of mating systems in flowering plants. Understanding the genomic basis of these changes is essential for predicting the persistence of a species under environmental changes. We investigated the genomic changes in Brassica rapa over nine generations of pollination by hoverflies associated with rapid morphological evolution toward the selfing syndrome. We combined a genotyping-by-sequencing (GBS) approach with a genome-wide association study (GWAS) to identify candidate genes, and assessed their functional role in the observed morphological changes by studying mutations of orthologous genes in the model plant Arabidopsis thaliana. We found 31 candidate genes involved in a wide range of functions from DNA/RNA binding to transport. Our functional assessment of orthologous genes in A. thaliana revealed that two of the identified genes in B. rapa are involved in regulating the size of floral organs. We found a protein kinase superfamily protein involved in petal width, an important trait in plant attractiveness to pollinators. Moreover, we found a histone lysine methyltransferase (HKMT) associated with stamen length. Altogether, our study shows that hoverfly pollination leads to rapid evolution toward the selfing syndrome mediated by polygenic changes.

Phylogenetic insight into ABCE gene subfamily in plants

ATP-BINDING CASSETTE SUBFAMILY E MEMBER (ABCE) proteins are one of the most conserved proteins across eukaryotes and archaea. Yeast and most animals possess a single ABCE gene encoding the critical translational factor ABCE1. In several plant species, including Arabidopsis thaliana and Oryza sativa, two or more ABCE gene copies have been identified, however information related to plant ABCE gene family is still missing. In this study we retrieved ABCE gene sequences of 76 plant species from public genome databases and comprehensively analyzed them with the reference to A. thaliana ABCE2 gene (AtABCE2). Using bioinformatic approach we assessed the conservation and phylogeny of plant ABCEs. In addition, we performed haplotype analysis of AtABCE2 and its paralogue AtABCE1 using genomic sequences of 1,135 A. thaliana ecotypes. Plant ABCE proteins showed overall high sequence conservation, sharing at least 78% of amino acid sequence identity with AtABCE2. We found that over half of the selected species have two to eight ABCE genes, suggesting that in plants ABCE genes can be classified as a low-copy gene family, rather than a single-copy gene family. The phylogenetic trees of ABCE protein sequences and the corresponding coding sequences demonstrated that Brassicaceae and Poaceae families have independently undergone lineage-specific split of the ancestral ABCE gene. Other plant species have gained ABCE gene copies through more recent duplication events. We also noticed that ploidy level but not ancient whole genome duplications experienced by a species impacts ABCE gene family size. Deeper analysis of AtABCE2 and AtABCE1 from 1,135 A. thaliana ecotypes revealed four and 35 non-synonymous SNPs, respectively. The lower natural variation in AtABCE2 compared to AtABCE1 is in consistence with its crucial role for plant viability. Overall, while the sequence of the ABCE protein family is highly conserved in the plant kingdom, many plants have evolved to have more than one copy of this essential translational factor.

N2-Acetylornithine deacetylase functions as a Cys-Gly dipeptidase in the cytosolic glutathione degradation pathway in Arabidopsis thaliana

Glutathione (GSH) is required for various physiological processes in plants, including redox regulation and detoxification of harmful compounds. GSH also functions as a repository for assimilated sulfur and is actively catabolized in plants. In Arabidopsis, GSH is mainly degraded initially by cytosolic enzymes, γ-glutamyl cyclotransferase, and γ-glutamyl peptidase, which release cysteinylglycine (Cys-Gly). However, the subsequent enzyme responsible for catabolizing this dipeptide has not been identified to date. In the present study, we identified At4g17830 as a Cys-Gly dipeptidase, namely cysteinylglycine peptidase 1 (CGP1). CGP1 complemented the phenotype of the yeast mutant that cannot degrade Cys-Gly. The Arabidopsis cgp1 mutant had lower Cys-Gly degradation activity than the wild type and showed perturbed concentrations of thiol compounds. Recombinant CGP1 showed reasonable Cys-Gly degradation activity in vitro. Metabolomic analysis revealed that cgp1 exhibited signs of severe sulfur deficiency, such as elevated accumulation of O-acetylserine (OAS) and the decrease in sulfur-containing metabolites. Morphological changes observed in cgp1, including longer primary roots of germinating seeds, were also likely associated with sulfur starvation. Notably, At4g17830 has previously been reported to encode an N2-acetylornithine deacetylase (NAOD) that functions in the ornithine biosynthesis. The cgp1 mutant did not show a decrease in ornithine content, whereas the analysis of CGP1 structure did not rule out the possibility that CGP1 has Cys-Gly dipeptidase and NAOD activities. Therefore, we propose that CGP1 is a Cys-Gly dipeptidase that functions in the cytosolic GSH degradation pathway and may play dual roles in GSH and ornithine metabolism.

Defying gravity: WEEP promotes negative gravitropism in peach trees by establishing asymmetric auxin gradients

Trees with weeping shoot architectures are valued for their beauty and are a resource for understanding how plants regulate posture control. The peach (Prunus persica) weeping phenotype, which has elliptical downward arching branches, is caused by a homozygous mutation in the WEEP gene. Little is known about the function of WEEP despite its high conservation throughout Plantae. Here, we present the results of anatomical, biochemical, biomechanical, physiological, and molecular experiments that provide insight into WEEP function. Our data suggest that weeping peach trees do not have defects in branch structure. Rather, transcriptomes from the adaxial (upper) and abaxial (lower) sides of standard and weeping branch shoot tips revealed flipped expression patterns for genes associated with early auxin response, tissue patterning, cell elongation, and tension wood development. This suggests that WEEP promotes polar auxin transport toward the lower side during shoot gravitropic response, leading to cell elongation and tension wood development. In addition, weeping peach trees exhibited steeper root systems and faster lateral root gravitropic response. This suggests that WEEP moderates root gravitropism and is essential to establishing the set-point angle of lateral roots from the gravity vector. Additionally, size exclusion chromatography indicated that WEEP proteins self-oligomerize, like other proteins with sterile alpha motif domains. Collectively, our results from weeping peach provide insight into polar auxin transport mechanisms associated with gravitropism and lateral shoot and root orientation.

Arabidopsis AN3 and OLIGOCELLULA genes link telomere maintenance mechanisms with cell division and expansion control

Telomeres are conserved chromosomal structures necessary for continued cell division and proliferation. In addition to the classical telomerase pathway, multiple other genes including those involved in ribosome metabolism and chromatin modification contribute to telomere length maintenance. We previously reported that Arabidopsis thaliana ribosome biogenesis genes OLI2/NOP2A, OLI5/RPL5A and OLI7/RPL5B have critical roles in telomere length regulation. These three OLIGOCELLULA genes were also shown to function in cell proliferation and expansion control and to genetically interact with the transcriptional co-activator ANGUSTIFOLIA3 (AN3). Here we show that AN3-deficient plants progressively lose telomeric DNA in early homozygous mutant generations, but ultimately establish a new shorter telomere length setpoint by the fifth mutant generation with a telomere length similar to oli2/nop2a -deficient plants. Analysis of double an3 oli2 mutants indicates that the two genes are epistatic for telomere length control. Telomere shortening in an3 and oli mutants is not caused by telomerase inhibition; wild type levels of telomerase activity are detected in all analyzed mutants in vitro. Late generations of an3 and oli mutants are prone to stem cell damage in the root apical meristem, implying that genes regulating telomere length may have conserved functional roles in stem cell maintenance mechanisms. Multiple instances of anaphase fusions in late generations of oli5 and oli7 mutants were observed, highlighting an unexpected effect of ribosome biogenesis factors on chromosome integrity. Overall, our data implicate AN3 transcription coactivator and OLIGOCELLULA proteins in the establishment of telomere length set point in plants and further suggest that multiple regulators with pleiotropic functions can connect telomere biology with cell proliferation and cell expansion pathways.

Two adjacent NLR genes conferring quantitative resistance to clubroot disease in Arabidopsis are regulated by a stably inherited epiallelic variation

Clubroot caused by the protist Plasmodiophora brassicae is a major disease affecting cultivated Brassicaceae. Using a combination of quantitative trait locus (QTL) fine mapping, CRISPR-Cas9 validation, and extensive analyses of DNA sequence and methylation patterns, we revealed that the two adjacent neighboring NLR (nucleotide-binding and leucine-rich repeat) genes AT5G47260 and AT5G47280 cooperate in controlling broad-spectrum quantitative partial resistance to the root pathogen P. brassicae in Arabidopsis and that they are epigenetically regulated. The variation in DNA methylation is not associated with any nucleotide variation or any transposable element presence/absence variants and is stably inherited. Variations in DNA methylation at the Pb-At5.2 QTL are widespread across Arabidopsis accessions and correlate negatively with variations in expression of the two genes. Our study demonstrates that natural, stable, and transgenerationally inherited epigenetic variations can play an important role in shaping resistance to plant pathogens by modulating the expression of immune receptors.

Parasitic dodder expresses an arsenal of secreted cellulases with multi-substrate specificity during host invasion

Cuscuta campestris is a common and problematic parasitic plant which relies on haustoria to connect to and siphon nutrients from host plants. Glycoside hydrolase family 9 (GH9) cellulases (EC 3.2.1.4) play critical roles in plant cell wall biosynthesis and disassembly, but their roles during Cuscuta host invasion remains underexplored. In this study, we identified 22 full-length GH9 cellulase genes in C. campestris genome, which encoded fifteen secreted and seven membrane-anchored cellulases that showed distinct phylogenetic relationships. Expression profiles suggested that some of the genes are involved in biosynthesis and remodeling of the parasite's cell wall during haustoriogenesis, while other genes encoding secreted B- and C-type cellulases are tentatively associated with degrading host cell walls during invasion. Transcriptomic data in a host-free system and in the presence of susceptible or partially resistant tomato hosts, showed for especially GH9B7, GH9B11 and GH9B12 a shift in expression profiles in the presence of hosts, being more highly expressed during host attachment, indicating that Cuscuta can tune cellulase expression in response to a host. Functional analyses of recombinant B- and C-type cellulases showed endoglucanase activities over wide pH and temperature conditions, and activities towards multiple cellulose and hemicellulose substrates. These findings improve our understanding of host cell wall disassembly by Cuscuta, and cellulase activity towards broad substrate range potentially explain its wide host range. This is the first study to provide a broad biochemical insight into Cuscuta GH9 cellulases, which based on our study may have potential applications in industrial bioprocessing.

Natural variation in CYCLIC NUCLEOTIDE-GATED ION CHANNEL 4 reveals a novel role of calcium signaling in vegetative phase change in Arabidopsis

The timing of vegetative phase change (VPC) in plants is regulated by a temporal decline in the expression of miR156. Both exogenous cues and endogenous factors, such as temperature, light, sugar, nutrients, and epigenetic regulators, have been shown to affect VPC by altering miR156 expression. However, the genetic basis of natural variation in VPC remains largely unexplored. Here, we conducted a genome-wide association study on the variation of the timing of VPC in Arabidopsis. We identified CYCLIC NUCLEOTIDE-GATED ION CHANNEL 4 (CNGC4) as a significant locus associated with the diversity of VPC. Mutations in CNGC4 delayed VPC, accompanied by an increased expression level of miR156 and a corresponding decrease in SQUAMOSA PROMOTER BINDING-LIKE (SPL) gene expression. Furthermore, mutations in CNGC2 and CATION EXCHANGER 1/3 (CAX1/3) also led to a delay in VPC. Polymorphisms in the CNGC4 promoter contribute to the natural variation in CNGC4 expression and the diversity of VPC. Specifically, the early CNGC4 variant promotes VPC and enhances plant adaptation to local environments. In summary, our findings offer genetic insights into the natural variation in VPC in Arabidopsis, and reveal a previously unidentified role of calcium signaling in the regulation of VPC.

A pan-genome of 69 Arabidopsis thaliana accessions reveals a conserved genome structure throughout the global species range

Although originally primarily a system for functional biology, Arabidopsis thaliana has, owing to its broad geographical distribution and adaptation to diverse environments, developed into a powerful model in population genomics. Here we present chromosome-level genome assemblies of 69 accessions from a global species range. We found that genomic colinearity is very conserved, even among geographically and genetically distant accessions. Along chromosome arms, megabase-scale rearrangements are rare and typically present only in a single accession. This indicates that the karyotype is quasi-fixed and that rearrangements in chromosome arms are counter-selected. Centromeric regions display higher structural dynamics, and divergences in core centromeres account for most of the genome size variations. Pan-genome analyses uncovered 32,986 distinct gene families, 60% being present in all accessions and 40% appearing to be dispensable, including 18% private to a single accession, indicating unexplored genic diversity. These 69 new Arabidopsis thaliana genome assemblies will empower future genetic research.

Organ-delimited gene regulatory networks provide high accuracy in candidate transcription factor selection across diverse processes

Organ-specific gene expression datasets that include hundreds to thousands of experiments allow the reconstruction of organ-level gene regulatory networks (GRNs). However, creating such datasets is greatly hampered by the requirements of extensive and tedious manual curation. Here, we trained a supervised classification model that can accurately classify the organ-of-origin for a plant transcriptome. This K-Nearest Neighbor-based multiclass classifier was used to create organ-specific gene expression datasets for the leaf, root, shoot, flower, and seed in Arabidopsis thaliana. A GRN inference approach was used to determine the: i. influential transcription factors (TFs) in each organ and, ii. most influential TFs for specific biological processes in that organ. These genome-wide, organ-delimited GRNs (OD-GRNs), recalled many known regulators of organ development and processes operating in those organs. Importantly, many previously unknown TF regulators were uncovered as potential regulators of these processes. As a proof-of-concept, we focused on experimentally validating the predicted TF regulators of lipid biosynthesis in seeds, an important food and biofuel trait. Of the top 20 predicted TFs, eight are known regulators of seed oil content, e.g., WRI1, LEC1, FUS3. Importantly, we validated our prediction of MybS2, TGA4, SPL12, AGL18, and DiV2 as regulators of seed lipid biosynthesis. We elucidated the molecular mechanism of MybS2 and show that it induces purple acid phosphatase family genes and lipid synthesis genes to enhance seed lipid content. This general approach has the potential to be extended to any species with sufficiently large gene expression datasets to find unique regulators of any trait-of-interest.

Arabidopsis leucine-rich repeat malectin receptor-like kinases regulate pollen-stigma interactions

Flowering plants contain tightly controlled pollen-pistil interactions required for promoting intraspecific fertilization and preventing interspecific hybridizations. In Arabidopsis (Arabidopsis thaliana), several receptor kinases (RKs) are known to regulate the later stages of intraspecific pollen tube growth and ovular reception in the pistil, but less is known about RK regulation of the earlier stages. The Arabidopsis RECEPTOR-LIKE KINASE IN FLOWERS1 (RKF1)/RKF1-LIKE (RKFL) 1-3 cluster of 4 leucine-rich repeat malectin (LRR-MAL) RKs was previously found to function in the stigma to promote intraspecific pollen hydration. In this study, we tested additional combinations of up to 7 Arabidopsis LRR-MAL RK knockout mutants, including RKF1, RKFL1-3, LysM RLK1-INTERACTING KINASE1, REMORIN-INTERACTING RECEPTOR1, and NEMATODE-INDUCED LRR-RLK2. These LRR-MAL RKs were discovered to function in the female stigma to support intraspecific Arabidopsis pollen tube growth and to establish a prezygotic interspecific barrier against Capsella rubella pollen. Thus, this study uncovered additional biological functions for this poorly understood group of RKs in regulating the early stages of Arabidopsis sexual reproduction.

A FLOWERING LOCUS T ortholog is associated with photoperiod-insensitive flowering in hemp (Cannabis sativa L.)

Hemp (Cannabis sativa L.) is an extraordinarily versatile crop, with applications ranging from medicinal compounds to seed oil and fibre products. Cannabis sativa is a short-day plant, and its flowering is highly controlled by photoperiod. However, substantial genetic variation exists for photoperiod sensitivity in C. sativa, and photoperiod-insensitive ("autoflower") cultivars are available. Using a bi-parental mapping population and bulked segregant analysis, we identified Autoflower2, a 0.5 Mbp locus significantly associated with photoperiod-insensitive flowering in hemp. Autoflower2 contains an ortholog of the central flowering time regulator FLOWERING LOCUS T (FT) from Arabidopsis thaliana which we termed CsFT1. We identified extensive sequence divergence between alleles of CsFT1 from photoperiod-sensitive and insensitive cultivars of C. sativa, including a duplication of CsFT1 and sequence differences, especially in introns. Furthermore, we observed higher expression of one of the CsFT1 copies found in the photoperiod-insensitive cultivar. Genotyping of several mapping populations and a diversity panel confirmed a correlation between CsFT1 alleles and photoperiod response, affirming that at least two independent loci involved in the photoperiodic control of flowering, Autoflower1 and Autoflower2, exist in the C. sativa gene pool. This study reveals the multiple independent origins of photoperiod insensitivity in C. sativa, supporting the likelihood of a complex domestication history in this species. By integrating the genetic relaxation of photoperiod sensitivity into novel C. sativa cultivars, expansion to higher latitudes will be permitted, thus allowing the full potential of this versatile crop to be reached.

The small iron-deficiency-induced protein OLIVIA and its relation to the bHLH transcription factor POPEYE

Iron (Fe) is a crucial micronutrient needed in many metabolic processes. To balance needs and potential toxicity, plants control the amount of Fe they take up and allocate to leaves and seeds during their development. One important regulator of this process is POPEYE (PYE). PYE is a Fe deficiency-induced key bHLH transcription factor (TF) for allocation of internal Fe in plants. In the absence of PYE, there is altered Fe translocation and plants develop a leaf chlorosis. NICOTIANAMINE SYNTHASE4 (NAS4), FERRIC-REDUCTION OXIDASE3 (FRO3), and ZINC-INDUCED FACILITATOR1 (ZIF1) genes are expressed at higher level in pye-1 indicating that PYE represses these genes. PYE activity is controlled in a yet unknown manner. Here, we show that a small Fe deficiency-induced protein OLIVIA (OLV) can interact with PYE. OLV has a conserved C-terminal motif, that we named TGIYY. Through deletion mapping, we pinpointed that OLV TGIYY and several regions of PYE can be involved in the protein interaction. An OLV overexpressing (OX) mutant line exhibited an enhanced NAS4 gene expression. This was a mild Fe deficiency response phenotype that was related to PYE function. Leaf rosettes of olv mutants remained smaller than those of wild type, indicating that OLV promotes plant growth. Taken together, our study identified a small protein OLV as a candidate that may connect aspects of Fe homeostasis with regulation of leaf growth.

Single-cell and spatial RNA sequencing reveal the spatiotemporal trajectories of fruit senescence

The senescence of fruit is a complex physiological process, with various cell types within the pericarp, making it highly challenging to elucidate their individual roles in fruit senescence. In this study, a single-cell expression atlas of the pericarp of pitaya (Hylocereus undatus) is constructed, revealing exocarp and mesocarp cells undergoing the most significant changes during the fruit senescence process. Pseudotime analysis establishes cellular differentiation and gene expression trajectories during senescence. Early-stage oxidative stress imbalance is followed by the activation of resistance in exocarp cells, subsequently senescence-associated proteins accumulate in the mesocarp cells at late-stage senescence. The central role of the early response factor HuCMB1 is unveiled in the senescence regulatory network. This study provides a spatiotemporal perspective for a deeper understanding of the dynamic senescence process in plants.

Comparative secretome analysis of Striga and Cuscuta species identifies candidate virulence factors for two evolutionarily independent parasitic plant lineages

BACKGROUND

Many parasitic plants of the genera Striga and Cuscuta inflict huge agricultural damage worldwide. To form and maintain a connection with a host plant, parasitic plants deploy virulence factors (VFs) that interact with host biology. They possess a secretome that represents the complement of proteins secreted from cells and like other plant parasites such as fungi, bacteria or nematodes, some secreted proteins represent VFs crucial to successful host colonisation. Understanding the genome-wide complement of putative secreted proteins from parasitic plants, and their expression during host invasion, will advance understanding of virulence mechanisms used by parasitic plants to suppress/evade host immune responses and to establish and maintain a parasite-host interaction.

RESULTS

We conducted a comparative analysis of the secretomes of root (Striga spp.) and shoot (Cuscuta spp.) parasitic plants, to enable prediction of candidate VFs. Using orthogroup clustering and protein domain analyses we identified gene families/functional annotations common to both Striga and Cuscuta species that were not present in their closest non-parasitic relatives (e.g. strictosidine synthase like enzymes), or specific to either the Striga or Cuscuta secretomes. For example, Striga secretomes were strongly associated with 'PAR1' protein domains. These were rare in the Cuscuta secretomes but an abundance of 'GMC oxidoreductase' domains were found, that were not present in the Striga secretomes. We then conducted transcriptional profiling of genes encoding putatively secreted proteins for the most agriculturally damaging root parasitic weed of cereals, S. hermonthica. A significant portion of the Striga-specific secretome set was differentially expressed during parasitism, which we probed further to identify genes following a 'wave-like' expression pattern peaking in the early penetration stage of infection. We identified 39 genes encoding putative VFs with functions such as cell wall modification, immune suppression, protease, kinase, or peroxidase activities, that are excellent candidates for future functional studies.

CONCLUSIONS

Our study represents a comprehensive secretome analysis among parasitic plants and revealed both similarities and differences in candidate VFs between Striga and Cuscuta species. This knowledge is crucial for the development of new management strategies and delaying the evolution of virulence in parasitic weeds.

Comprehensive deciphering the alternative splicing patterns involved in leaf morphogenesis of Liriodendron chinense

Alternative splicing (AS), a pivotal post-transcriptional regulatory mechanism, profoundly amplifies diversity and complexity of transcriptome and proteome. Liriodendron chinense (Hemsl.) Sarg., an excellent ornamental tree species renowned for its distinctive leaf shape, which resembles the mandarin jacket. Despite the documented potential genes related to leaf development of L. chinense, the underlying post-transcriptional regulatory mechanisms remain veiled. Here, we conducted a comprehensive analysis of the transcriptome to clarify the genome-wide landscape of the AS pattern and the spectrum of spliced isoforms during leaf developmental stages in L. chinense. Our investigation unveiled 50,259 AS events, involving 10,685 genes (32.9%), with intron retention as the most prevalent events. Notably, the initial stage of leaf development witnessed the detection of 804 differentially AS events affiliated with 548 genes. Although both differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) were enriched into morphogenetic related pathways during the transition from fishhook (P2) to lobed (P7) leaves, there was only a modest degree of overlap between DASGs and DEGs. Furthermore, we conducted a comprehensively AS analysis on homologous genes involved in leaf morphogenesis, and most of which are subject to post-transcriptional regulation of AS. Among them, the AINTEGUMENTA-LIKE transcript factor LcAIL5 was characterization in detailed, which experiences skipping exon (SE), and two transcripts displayed disparate expression patterns across multiple stages. Overall, these findings yield a comprehensive understanding of leaf development regulation via AS, offering a novel perspective for further deciphering the mechanism of plant leaf morphogenesis.

A K homology (KH) domain protein identified by a forward genetic screen affects bundle sheath anatomy in Arabidopsis thaliana

Because of their photosynthetic capacity, leaves function as solar panels providing the basis for the growth of the entire plant. Although the molecular mechanisms of leaf development have been well studied in model dicot and monocot species, a lot of information is still needed about the interplay of the genes that regulate cell division and differentiation and thereby affect the photosynthetic performance of the leaf. We were specifically interested in understanding the differentiation of mesophyll and bundle sheath cells in Arabidopsis thaliana and aimed to identify genes that are involved in determining bundle sheath anatomy. To this end, we established a forward genetic screen by using ethyl methanesulfonate (EMS) for mutagenizing a reporter line expressing a chloroplast-targeted green fluorescent protein (sGFP) under the control of a bundle sheath-specific promoter. Based on the GFP fluorescence phenotype, numerous mutants were produced, and by pursuing a mapping-by-sequencing approach, the genomic segments containing mutated candidate genes were identified. One of the lines with an enhanced GFP fluorescence phenotype (named ELEVATED BUNDLE SHEATH CELLS SIGNAL 1 [ebss1]) was selected for further study, and the responsible gene was verified by CRISPR/Cas9-based mutagenesis of candidate genes located in the mapped genomic segment. The verified gene, At2g25970, encodes a K homology (KH) domain-containing protein.

Profiling of the gene expression and alternative splicing landscapes of Eucalyptus grandis

Eucalyptus is a widely planted hardwood tree species due to its fast growth, superior wood properties and adaptability. However, the post-transcriptional regulatory mechanisms controlling tissue development and stress responses in Eucalyptus remain poorly understood. In this study, we performed a comprehensive analysis of the gene expression profile and the alternative splicing (AS) landscape of E. grandis using strand-specific RNA-Seq, which encompassed 201 libraries including different organs, developmental stages, and environmental stresses. We identified 10 416 genes (33.49%) that underwent AS, and numerous differentially expressed and/or differential AS genes involved in critical biological processes, such as primary-to-secondary growth transition of stems, adventitious root formation, aging and responses to phosphorus- or boron-deficiency. Co-expression analysis of AS events and gene expression patterns highlighted the potential upstream regulatory role of AS events in multiple processes. Additionally, we highlighted the lignin biosynthetic pathway to showcase the potential regulatory functions of AS events in the KNAT3 and IRL3 genes within this pathway. Our high-quality expression atlas and AS landscape serve as valuable resources for unravelling the genetic control of woody plant development, long-term adaptation, and understanding transcriptional diversity in Eucalyptus. Researchers can conveniently access these resources through the interactive ePlant browser (https://bar.utoronto.ca/eplant_eucalyptus).

A novel tetratricopeptide-repeat protein, TTP1, forms complexes with glutamyl-tRNA reductase and protochlorophyllide oxidoreductase during tetrapyrrole biosynthesis

The biosynthesis of the tetrapyrrole end-products chlorophyll and heme depends on a multifaceted control mechanism that acts primarily at the post-translational level upon the rate-limiting step of 5-aminolevulinic acid synthesis and upon light-dependent protochlorophyllide oxidoreductase (POR). These regulatory processes require auxiliary factors that modulate the activity, stability, complex formation, and subplastidal localization of the relevant proteins. Together, they ensure optimal metabolic flow during the day and at night. As an Arabidopsis homolog of the POR-interacting tetratricopeptide-repeat protein (Pitt) first reported in Synechocystis, we characterize tetrapyrrole biosynthesis-regulating tetratricopeptide-repeat protein1 (TTP1). TTP1 is a plastid-localized, membrane-bound factor that interacts with POR, the Mg protoporphyrin monomethylester cyclase CHL27, glutamyl-tRNA reductase (GluTR), GluTR-binding protein, and FLUORESCENCE IN BLUE LIGHT. Lack of TTP1 leads to accumulation of GluTR, enhanced 5-aminolevulinic acid synthesis and lower levels of POR. Knockout mutants show enhanced sensitivity to reactive oxygen species and a slower greening of etiolated seedlings. Based on our studies, the interaction of TTP1 with GluTR and POR does not directly inhibit their enzymatic activity and contribute to the control of 5-aminolevulinic acid synthesis. Instead, we propose that TTP1 sequesters a fraction of these proteins on the thylakoid membrane, and contributes to their stability.

Mutation of negative regulatory gene CEHC1 encoding an FBXO3 protein results in normoxic expression of HYDA genes in Chlamydomonas reinhardtii

Oxygen is known to prevent hydrogen production in Chlamydomonas, both by inhibiting the hydrogenase enzyme and by preventing the accumulation of HYDA-encoding transcripts. We developed a screen for mutants showing constitutive accumulation of HYDA1 transcripts in the presence of oxygen. A reporter gene required for ciliary motility, placed under the control of the HYDA1 promoter, conferred motility only in hypoxic conditions. By selecting for mutants able to swim even in the presence of oxygen we obtained strains that express the reporter gene constitutively. One mutant identified a gene encoding an F-box only protein 3 (FBXO3), known to participate in ubiquitylation and proteasomal degradation pathways in other eukaryotes. Transcriptome profiles revealed that the mutation, termed cehc1-1 , leads to constitutive expression of HYDA1 and other genes regulated by hypoxia, and of many genes known to be targets of CRR1, a transcription factor in the nutritional copper signaling pathway. CRR1 was required for the constitutive expression of the HYDA1 reporter gene in cehc1-1 mutants. The CRR1 protein, which is normally degraded in Cu-supplemented cells, was stabilized in cehc1-1 cells, supporting the conclusion that CEHC1 acts to facilitate the degradation of CRR1. Our results reveal a novel negative regulator in the CRR1 pathway and possibly other pathways leading to complex metabolic changes associated with response to hypoxia.

Nucleotide-level distance metrics to quantify alternative splicing implemented in TranD

Advances in affordable transcriptome sequencing combined with better exon and gene prediction has motivated many to compare transcription across the tree of life. We develop a mathematical framework to calculate complexity and compare transcript models. Structural features, i.e. intron retention (IR), donor/acceptor site variation, alternative exon cassettes, alternative 5'/3' UTRs, are compared and the distance between transcript models is calculated with nucleotide level precision. All metrics are implemented in a PyPi package, TranD and output can be used to summarize splicing patterns for a transcriptome (1GTF) and between transcriptomes (2GTF). TranD output enables quantitative comparisons between: annotations augmented by empirical RNA-seq data and the original transcript models; transcript model prediction tools for longread RNA-seq (e.g. FLAIR versus Isoseq3); alternate annotations for a species (e.g. RefSeq vs Ensembl); and between closely related species. In C. elegans, Z. mays, D. melanogaster, D. simulans and H. sapiens, alternative exons were observed more frequently in combination with an alternative donor/acceptor than alone. Transcript models in RefSeq and Ensembl are linked and both have unique transcript models with empirical support. D. melanogaster and D. simulans, share many transcript models and long-read RNAseq data suggests that both species are under-annotated. We recommend combined references.

Dissection of quantitative trait nucleotides and candidate genes associated with agronomic and yield-related traits under drought stress in rapeseed varieties: integration of genome-wide association study and transcriptomic analysis

Introduction

An important strategy to combat yield loss challenge is the development of varieties with increased tolerance to drought to maintain production. Improvement of crop yield under drought stress is critical to global food security.

Methods

In this study, we performed multiomics analysis in a collection of 119 diverse rapeseed (Brassica napus L.) varieties to dissect the genetic control of agronomic traits in two watering regimes [well-watered (WW) and drought stress (DS)] for 3 years. In the DS treatment, irrigation continued till the 50% pod development stage, whereas in the WW condition, it was performed throughout the whole growing season.

Results

The results of the genome-wide association study (GWAS) using 52,157 single-nucleotide polymorphisms (SNPs) revealed 1,281 SNPs associated with traits. Six stable SNPs showed sequence variation for flowering time between the two irrigation conditions across years. Three novel SNPs on chromosome C04 for plant weight were located within drought tolerance-related gene ABCG16, and their pleiotropically effects on seed weight per plant and seed yield were characterized. We identified the C02 peak as a novel signal for flowering time, harboring 52.77% of the associated SNPs. The 288-kbps LD decay distance analysis revealed 2,232 candidate genes (CGs) associated with traits. The CGs BIG1-D, CAND1, DRG3, PUP10, and PUP21 were involved in phytohormone signaling and pollen development with significant effects on seed number, seed weight, and grain yield in drought conditions. By integrating GWAS and RNA-seq, 215 promising CGs were associated with developmental process, reproductive processes, cell wall organization, and response to stress. GWAS and differentially expressed genes (DEGs) of leaf and seed in the yield contrasting accessions identified BIG1-D, CAND1, and DRG3 genes for yield variation.

Discussion

The results of our study provide insights into the genetic control of drought tolerance and the improvement of marker-assisted selection (MAS) for breeding high-yield and drought-tolerant varieties.

Complex epistatic interactions between ELF3, PRR9, and PRR7 regulate the circadian clock and plant physiology

Circadian clocks are endogenous timekeeping mechanisms that coordinate internal physiological responses with the external environment. EARLY FLOWERING3 (ELF3), PSEUDO RESPONSE REGULATOR (PRR9), and PRR7 are essential components of the plant circadian clock and facilitate entrainment of the clock to internal and external stimuli. Previous studies have highlighted a critical role for ELF3 in repressing the expression of PRR9 and PRR7. However, the functional significance of activity in regulating circadian clock dynamics and plant development is unknown. To explore this regulatory dynamic further, we first employed mathematical modeling to simulate the effect of the prr9/prr7 mutation on the elf3 circadian phenotype. These simulations suggested that simultaneous mutations in prr9/prr7 could rescue the elf3 circadian arrhythmia. Following these simulations, we generated all Arabidopsis elf3/prr9/prr7 mutant combinations and investigated their circadian and developmental phenotypes. Although these assays could not replicate the results from the mathematical modeling, our results have revealed a complex epistatic relationship between ELF3 and PRR9/7 in regulating different aspects of plant development. ELF3 was essential for hypocotyl development under ambient and warm temperatures, while PRR9 was critical for root thermomorphogenesis. Finally, mutations in prr9 and prr7 rescued the photoperiod-insensitive flowering phenotype of the elf3 mutant. Together, our results highlight the importance of investigating the genetic relationship among plant circadian genes.

A review of artificial intelligence-assisted omics techniques in plant defense: current trends and future directions

Plants intricately deploy defense systems to counter diverse biotic and abiotic stresses. Omics technologies, spanning genomics, transcriptomics, proteomics, and metabolomics, have revolutionized the exploration of plant defense mechanisms, unraveling molecular intricacies in response to various stressors. However, the complexity and scale of omics data necessitate sophisticated analytical tools for meaningful insights. This review delves into the application of artificial intelligence algorithms, particularly machine learning and deep learning, as promising approaches for deciphering complex omics data in plant defense research. The overview encompasses key omics techniques and addresses the challenges and limitations inherent in current AI-assisted omics approaches. Moreover, it contemplates potential future directions in this dynamic field. In summary, AI-assisted omics techniques present a robust toolkit, enabling a profound understanding of the molecular foundations of plant defense and paving the way for more effective crop protection strategies amidst climate change and emerging diseases.

Transcriptomic and proteomic analyses of distinct Arabidopsis organs reveal high PSI-NDH complex accumulation in stems

In addition to leaves, the main site of photosynthetic reactions, active photosynthesis also takes place in stems, siliques and tree trunks. Although non-foliar photosynthesis has a marked effect on plant growth and yield, only limited information on the expression patterns of photosynthesis-related genes and the structure of photosynthetic machinery in different plant organs has been available. Here, we report the results of transcriptomic analysis of various organs of Arabidopsis thaliana and compare the gene expression profiles of young and mature leaves with a special focus on photosynthetic genes. Further, we analyzed the composition and organization of the photosynthetic electron transfer machinery in leaves, stems and green siliques at the protein level using BN-PAGE. RNA-Seq analysis revealed unique gene expression profiles in different plant organs and showed major differences in the expression of photosynthesis-related genes in young as compared to mature rosettes. Gel-based proteomic analysis of the thylakoid protein complex organization further showed that all studied plant organs contain the necessary components of the photosynthetic electron transfer chain. Intriguingly, stems accumulate high amounts of PSI-NDH complex, which has previously been implicated in cyclic electron transfer.

Proteome plasticity during Physcomitrium patens spore germination - from the desiccated phase to heterotrophic growth and reconstitution of photoautotrophy

The establishment of moss spores is considered a milestone in plant evolution. They harbor protein networks underpinning desiccation tolerance and accumulation of storage compounds that can be found already in algae and that are also utilized in seeds and pollen. Furthermore, germinating spores must produce proteins that drive the transition through heterotrophic growth to the autotrophic plant. To get insight into the plasticity of this proteome, we investigated it at five timepoints of moss (Physcomitrium patens) spore germination and in protonemata and gametophores. The comparison to previously published Arabidopsis proteome data of seedling establishment showed that not only the proteomes of spores and seeds are functionally related, but also the proteomes of germinating spores and young seedlings. We observed similarities with regard to desiccation tolerance, lipid droplet proteome composition, control of dormancy, and β-oxidation and the glyoxylate cycle. However, there were also striking differences. For example, spores lacked any obvious storage proteins. Furthermore, we did not detect homologs to the main triacylglycerol lipase in Arabidopsis seeds, SUGAR DEPENDENT1. Instead, we discovered a triacylglycerol lipase of the oil body lipase family and a lipoxygenase as being the overall most abundant proteins in spores. This finding indicates an alternative pathway for triacylglycerol degradation via oxylipin intermediates in the moss. The comparison of spores to Nicotiana tabacum pollen indicated similarities for example in regards to resistance to desiccation and hypoxia, but the overall developmental pattern did not align as in the case of seedling establishment and spore germination.

PSEUDO-ETIOLATION IN LIGHT proteins reduce greening by binding GLK transcription factors

Knocking out genes encoding proteins that downregulate the accumulation of pigments may lead to increases in crop quality and yield. PSEUDO-ETIOLATION IN LIGHT 1 (PEL1) downregulates the accumulation of carotenoids in carrot and chlorophyll in Arabidopsis and rice and may inhibit GOLDEN 2-LIKE (GLK) transcription factors. PEL1 belongs to a previously unstudied gene family found only in plants. We used CRISPR/Cas9 technology to knock out each member of the 4-member PEL gene family and both GLK genes in Arabidopsis. In pel mutants, chlorophyll levels were elevated in seedlings; after flowering, chloroplasts increased in size, and anthocyanin levels increased. Although the chlorophyll-deficient phenotype of glk1 glk2 was epistatic to pel1 pel2 pel3 pel4 in most of our experiments, glk1 glk2 was not epistatic to pel1 pel2 pel3 pel4 for the accumulation of anthocyanins in most of our experiments. The pel alleles attenuated growth, altered the accumulation of nutrients in seeds, disrupted an abscisic acid-inducible inhibition of seedling growth response that promotes drought tolerance, and affected the expression of genes associated with diverse biological functions, such as stress responses, cell wall metabolism hormone responses, signaling, growth, and the accumulation of phenylpropanoids and pigments. We found that PEL proteins specifically bind 6 transcription factors that influence the accumulation of anthocyanins, GLK2, and the carboxy termini of GLK1 and Arabidopsis thaliana myeloblastosis oncogene homolog 4 (AtMYB4). Our data indicate that the PEL proteins influence the accumulation of chlorophyll and many other processes, possibly by inhibiting GLK transcription factors and via other mechanisms, and that multiple mechanisms downregulate chlorophyll content.

Syntrichia ruralis: emerging model moss genome reveals a conserved and previously unknown regulator of desiccation in flowering plants

Water scarcity, resulting from climate change, poses a significant threat to ecosystems. Syntrichia ruralis, a dryland desiccation-tolerant moss, provides valuable insights into survival of water-limited conditions. We sequenced the genome of S. ruralis, conducted transcriptomic analyses, and performed comparative genomic and transcriptomic analyses with existing genomes and transcriptomes, including with the close relative S. caninervis. We took a genetic approach to characterize the role of an S. ruralis transcription factor, identified in transcriptomic analyses, in Arabidopsis thaliana. The genome was assembled into 12 chromosomes encompassing 21 169 protein-coding genes. Comparative analysis revealed copy number and transcript abundance differences in known desiccation-associated gene families, and highlighted genome-level variation among species that may reflect adaptation to different habitats. A significant number of abscisic acid (ABA)-responsive genes were found to be negatively regulated by a MYB transcription factor (MYB55) that was upstream of the S. ruralis ortholog of ABA-insensitive 3 (ABI3). We determined that this conserved MYB transcription factor, uncharacterized in Arabidopsis, acts as a negative regulator of an ABA-dependent stress response in Arabidopsis. The new genomic resources from this emerging model moss offer novel insights into how plants regulate their responses to water deprivation.

Petal abscission is promoted by jasmonic acid-induced autophagy at Arabidopsis petal bases

In angiosperms, the transition from floral-organ maintenance to abscission determines reproductive success and seed dispersion. For petal abscission, cell-fate decisions specifically at the petal-cell base are more important than organ-level senescence or cell death in petals. However, how this transition is regulated remains unclear. Here, we identify a jasmonic acid (JA)-regulated chromatin-state switch at the base of Arabidopsis petals that directs local cell-fate determination via autophagy. During petal maintenance, co-repressors of JA signaling accumulate at the base of petals to block MYC activity, leading to lower levels of ROS. JA acts as an airborne signaling molecule transmitted from stamens to petals, accumulating primarily in petal bases to trigger chromatin remodeling. This allows MYC transcription factors to promote chromatin accessibility for downstream targets, including NAC DOMAIN-CONTAINING PROTEIN102 (ANAC102). ANAC102 accumulates specifically at the petal base prior to abscission and triggers ROS accumulation and cell death via AUTOPHAGY-RELATED GENEs induction. Developmentally induced autophagy at the petal base causes maturation, vacuolar delivery, and breakdown of autophagosomes for terminal cell differentiation. Dynamic changes in vesicles and cytoplasmic components in the vacuole occur in many plants, suggesting JA-NAC-mediated local cell-fate determination by autophagy may be conserved in angiosperms.

Variability in drought gene expression datasets highlight the need for community standardization

Physiologically relevant drought stress is difficult to apply consistently, and the heterogeneity in experimental design, growth conditions, and sampling schemes make it challenging to compare water deficit studies in plants. Here, we re-analyzed hundreds of drought gene expression experiments across diverse model and crop species and quantified the variability across studies. We found that drought studies are surprisingly uncomparable, even when accounting for differences in genotype, environment, drought severity, and method of drying. Many studies, including most Arabidopsis work, lack high-quality phenotypic and physiological datasets to accompany gene expression, making it impossible to assess the severity or in some cases the occurrence of water deficit stress events. From these datasets, we developed supervised learning classifiers that can accurately predict if RNA-seq samples have experienced a physiologically relevant drought stress, and suggest this can be used as a quality control for future studies. Together, our analyses highlight the need for more community standardization, and the importance of paired physiology data to quantify stress severity for reproducibility and future data analyses.

A novel fluorescent protein pair facilitates FLIM-FRET analysis of plant immune receptor interaction under native conditions

Elucidating protein-protein interactions is crucial for our understanding of molecular processes within living organisms. Microscopy-based techniques can detect protein-protein interactions in vivo at the single-cell level and provide information on their subcellular location. Fluorescence lifetime imaging microscopy (FLIM)-Förster resonance energy transfer (FRET) is one of the most robust imaging approaches, but it is still very challenging to apply this method to proteins which are expressed under native conditions. Here we describe a novel combination of fluorescence proteins (FPs), mCitrine and mScarlet-I, which is ideally suited for FLIM-FRET studies of low abundance proteins expressed from their native promoters in stably transformed plants. The donor mCitrine displays excellent brightness in planta, near-mono-exponential fluorescence decay, and a comparatively long fluorescence lifetime. Moreover, the FRET pair has a good spectral overlap and a large Förster radius. This allowed us to detect constitutive as well as ligand-induced interaction of the Arabidopsis chitin receptor components CERK1 and LYK5 in a set of proof-of-principle experiments. Due to the good brightness of the acceptor mScarlet-I, the FP combination can be readily utilized for co-localization studies. The FP pair is also suitable for co-immunoprecipitation experiments and western blotting, facilitating a multi-method approach for studying and confirming protein-protein interactions.

Investigating low frequency somatic mutations in Arabidopsis with Duplex Sequencing

Mutations are the source of novel genetic diversity but can also lead to disease and maladaptation. The conventional view is that mutations occur randomly with respect to their environment-specific fitness consequences. However, intragenomic mutation rates can vary dramatically due to transcription coupled repair and based on local epigenomic modifications, which are non-uniformly distributed across genomes. One sequence feature associated with decreased mutation is higher expression level, which can vary depending on environmental cues. To understand whether the association between expression level and mutation rate creates a systematic relationship with environment-specific fitness effects, we perturbed expression through a heat treatment in Arabidopsis thaliana. We quantified gene expression to identify differentially expressed genes, which we then targeted for mutation detection using Duplex Sequencing. This approach provided a highly accurate measurement of the frequency of rare somatic mutations in vegetative plant tissues, which has been a recent source of uncertainty in plant mutation research. We included mutant lines lacking mismatch repair (MMR) and base excision repair (BER) capabilities to understand how repair mechanisms may drive biased mutation accumulation. We found wild type (WT) and BER mutant mutation frequencies to be very low (mean variant frequency 1.8×10-8 and 2.6×10-8, respectively), while MMR mutant frequencies were significantly elevated (1.13×10-6). These results show that somatic variant frequencies are extremely low in WT plants, indicating that larger datasets will be needed to address the fundamental evolutionary question as to whether environmental change leads to gene-specific changes in mutation rate.

ALLENE OXIDE SYNTHASE (AOS) induces petal senescence through a novel JA-associated regulatory pathway in Arabidopsis

Flowers are crucial for the reproduction of flowering plants and their senescence has drastic effects on plant-animal interactions as well as pollination. Petal senescence is the final phase of flower development which is regulated by hormones and genes. Among these, jasmonic acid (JA) has emerged as a major contributor to petal senescence, but its molecular mechanisms remain elusive. Here, the role of JA in petal senescence in Arabidopsis was investigated. We showed that petal senescence in aos mutant was significantly delayed, which also affected petal cell size and proliferation. Similar significant delays in petal senescence were observed in dad1 and coi1 mutants. However, MYB21/24 and MYC2/3/4, known downstream regulators of JA in flower development, played no role in petal senescence. This indicated that JA regulates petal senescence by modulating other unknown transcription factors. Transcriptomic analysis revealed that AOS altered the expression of 3681 genes associated, and identified groups of differentially expressed transcription factors, highlighting the potential involvement of AP-2, WRKY and NAC. Furthermore, bHLH13, bHLH17 and URH2 were identified as potential new regulators of JA-mediated petal senescence. In conclusion, our findings suggest a novel genetic pathway through which JA regulates petal senescence in Arabidopsis. This pathway operates independently of stamen development and leaf senescence, suggesting the evolution of specialized mechanisms for petal senescence.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01425-w.

Genome wide association study of Arabidopsis seed mucilage layers at a regional scale

The myxospermous species Arabidopsis thaliana extrudes a polysaccharidic mucilage from the seed coat epidermis during imbibition. The whole seed mucilage can be divided into a seed-adherent layer and a fully soluble layer, both layers presenting natural genetic variations. The adherent mucilage is variable in size and composition, while the soluble mucilage is variable in composition and physical properties. Studies reporting both the genetic architecture and the putative selective agents acting on this natural genetic variation are scarce. In this study, we set up a Genome Wide Association study (GWAS) based on 424 natural accessions collected from 166 natural populations of A. thaliana located south-west of France and previously characterized for a very important number of abiotic and biotic factors. We identified an extensive genetic variation for both mucilage layers. The adherent mucilage was mainly related to precipitation and temperature whereas the non-adherent mucilage was unrelated to any environmental factors. By combining a hierarchical Bayesian model with a local score approach, we identified 55 and 28 candidate genes, corresponding to 26 and 10 QTLs for the adherent and non-adherent mucilages, respectively. Putative or characterized function and expression data available in the literature were used to filter the candidate genes. Only one gene among our set of candidate genes was already described as a seed mucilage actor, leaving a large set of new candidates putatively implicated inseed mucilage synthesis or release. The present study lay out foundation to understand the influence of regional ecological factors acting on seed mucilage in A. thaliana.

Research progress on Brassicaceae plants: a bibliometrics analysis

The Brassicaceae is a worldwide family that produces ornamental flowers, edible vegetables, and oilseed plants, with high economic value in agriculture, horticulture, and landscaping. This study used the Web of Science core dataset and the CiteSpace bibliometric tool to quantitatively visualize the number of publications, authors, institutions, and countries of 3139 papers related to Brassicaceae plants from 2002 to 2022. The keywords and references were divided into two phases: Phase 1 (2002-2011) and Phase 2 (2012-2022) for quantitative and qualitative analysis. The results showed: An average annual publication volume of 149 articles, with an overall fluctuating upward trend; the research force was mainly led by Professor Ihsan A. Al-shehbaz from Missouri Botanical Garden; and the United States had the highest number of publications. In the first phase, research focused on the phylogeny of Brassicaceae plants, while the second phase delved into diverse research based on previous studies, research in areas such as polyploidy, molecular technique, physiology, and hyperaccumulator has been extended. Based on this research, we propounded some ideas for future studies on Brassicaceae plants and summarized the research gaps.

Plant mRNAs move into a fungal pathogen via extracellular vesicles to reduce infection

Cross-kingdom small RNA trafficking between hosts and microbes modulates gene expression in the interacting partners during infection. However, whether other RNAs are also transferred is unclear. Here, we discover that host plant Arabidopsis thaliana delivers mRNAs via extracellular vesicles (EVs) into the fungal pathogen Botrytis cinerea. A fluorescent RNA aptamer reporter Broccoli system reveals host mRNAs in EVs and recipient fungal cells. Using translating ribosome affinity purification profiling and polysome analysis, we observe that delivered host mRNAs are translated in fungal cells. Ectopic expression of two transferred host mRNAs in B. cinerea shows that their proteins are detrimental to infection. Arabidopsis knockout mutants of the genes corresponding to these transferred mRNAs are more susceptible. Thus, plants have a strategy to reduce infection by transporting mRNAs into fungal cells. mRNAs transferred from plants to pathogenic fungi are translated to compromise infection, providing knowledge that helps combat crop diseases.

BPG4 regulates chloroplast development and homeostasis by suppressing GLK transcription factors and involving light and brassinosteroid signaling

Chloroplast development adapts to the environment for performing suitable photosynthesis. Brassinosteroids (BRs), plant steroid hormones, have crucial effects on not only plant growth but also chloroplast development. However, the detailed molecular mechanisms of BR signaling in chloroplast development remain unclear. Here, we identify a regulator of chloroplast development, BPG4, involved in light and BR signaling. BPG4 interacts with GOLDEN2-LIKE (GLK) transcription factors that promote the expression of photosynthesis-associated nuclear genes (PhANGs), and suppresses their activities, thereby causing a decrease in the amounts of chlorophylls and the size of light-harvesting complexes. BPG4 expression is induced by BR deficiency and light, and is regulated by the circadian rhythm. BPG4 deficiency causes increased reactive oxygen species (ROS) generation and damage to photosynthetic activity under excessive high-light conditions. Our findings suggest that BPG4 acts as a chloroplast homeostasis factor by fine-tuning the expression of PhANGs, optimizing chloroplast development, and avoiding ROS generation.

Pectate lyase-like lubricates the male gametophyte's path toward its mating partner

The pollen tube is an extension of the male gametophyte in plants and mediates sexual reproduction by delivering the sperm cells to the female gametophyte. To accomplish this task, the elongating pollen tube must break through the thick wall of the pollen grain and penetrate multiple pistillar tissues. Both processes require the loosening of cell wall material-that of the pollen intine and that of the apoplast of the transmitting tract. The enzymatic toolbox for these cell wall modifying processes employed by the invading male gametophyte is elusive. We investigated the role of the pectin-digesting pectate lyase-like (PLL) by combining mutant analysis with microscopy observations, fluorescence recovery after photo-bleaching experiments, and immuno-detection. We show that in Arabidopsis (Arabidopsis thaliana), PLLs are required for intine loosening during the first steps of pollen tube germination. We provide evidence that during pollen tube elongation, PLLs are released by the pollen tube into the extracellular space, suggesting that they may be employed to soften the apoplast of the transmitting tissue. The synergistic enzymatic action of PLLs in the pollen grain, the pollen tube, and the transmitting track contribute to an effective fertilization process.

Old school, new rules: floral meristem development revealed by 3D gene expression atlases and high-resolution transcription factor-chromatin dynamics

The intricate morphology of the flower is primarily established within floral meristems in which floral organs will be defined and from where the developing flower will emerge. Floral meristem development involves multiscale-level regulation, including lineage and positional mechanisms for establishing cell-type identity, and transcriptional regulation mediated by changes in the chromatin environment. However, many key aspects of floral meristem development remain to be determined, such as: 1) the exact role of cellular location in connecting transcriptional inputs to morphological outcomes, and 2) the precise interactions between transcription factors and chromatin regulators underlying the transcriptional networks that regulate the transition from cell proliferation to differentiation during floral meristem development. Here, we highlight recent studies addressing these points through newly developed spatial reconstruction techniques and high-resolution transcription factor-chromatin environment interactions in the model plant Arabidopsis thaliana. Specifically, we feature studies that reconstructed 3D gene expression atlases of the floral meristem. We also discuss how the precise timing of floral meristem specification, floral organ patterning, and floral meristem termination is determined through temporally defined epigenetic dynamics for fine-tuning of gene expression. These studies offer fresh insights into the well-established principles of floral meristem development and outline the potential for further advances in this field in an age of integrated, powerful, multiscale resolution approaches.

A micro RNA mediates shoot control of root branching

Plants extract mineral nutrients from the soil, or from interactions with mutualistic soil microbes via their root systems. Adapting root architecture to nutrient availability enables efficient resource utilization, particularly in patchy and dynamic environments. Root growth responses to soil nitrogen levels are shoot-mediated, but the identity of shoot-derived mobile signals regulating root growth responses has remained enigmatic. Here we show that a shoot-derived micro RNA, miR2111, systemically steers lateral root initiation and nitrogen responsiveness through its root target TML (TOO MUCH LOVE) in the legume Lotus japonicus, where miR2111 and TML were previously shown to regulate symbiotic infections with nitrogen fixing bacteria. Intriguingly, systemic control of lateral root initiation by miR2111 and TML/HOLT (HOMOLOGUE OF LEGUME TML) was conserved in the nonsymbiotic ruderal Arabidopsis thaliana, which follows a distinct ecological strategy. Thus, the miR2111-TML/HOLT regulon emerges as an essential, conserved factor in adaptive shoot control of root architecture in dicots.

SCEP1 and SCEP2 are two new components of the synaptonemal complex central element

The synaptonemal complex (SC) is a proteinaceous structure that forms between homologous chromosomes during meiosis prophase. The SC is widely conserved across species, but its structure and roles during meiotic recombination are still debated. While the SC central region is made up of transverse filaments and central element proteins in mammals and fungi, few central element proteins have been identified in other species. Here we report the identification of two coiled-coil proteins, SCEP1 and SCEP2, that form a complex and localize at the centre of the Arabidopsis thaliana SC. In scep1 and scep2 mutants, chromosomes are aligned but not synapsed (the ZYP1 transverse filament protein is not loaded), crossovers are increased compared with the wild type, interference is lost and heterochiasmy is strongly reduced. We thus report the identification of two plant SC central elements, and homologues of these are found in all major angiosperm clades.

Subgenome dominance shapes novel gene evolution in the decaploid pitcher plant Nepenthes gracilis

Subgenome dominance after whole-genome duplication generates distinction in gene number and expression at the level of chromosome sets, but it remains unclear how this process may be involved in evolutionary novelty. Here we generated a chromosome-scale genome assembly of the Asian pitcher plant Nepenthes gracilis to analyse how its novel traits (dioecy and carnivorous pitcher leaves) are linked to genomic evolution. We found a decaploid karyotype and a clear indication of subgenome dominance. A male-linked and pericentromerically located region on the putative sex chromosome was identified in a recessive subgenome and was found to harbour three transcription factors involved in flower and pollen development, including a likely neofunctionalized LEAFY duplicate. Transcriptomic and syntenic analyses of carnivory-related genes suggested that the paleopolyploidization events seeded genes that subsequently formed tandem clusters in recessive subgenomes with specific expression in the digestive zone of the pitcher, where specialized cells digest prey and absorb derived nutrients. A genome-scale analysis suggested that subgenome dominance likely contributed to evolutionary innovation by permitting recessive subgenomes to diversify functions of novel tissue-specific duplicates. Our results provide insight into how polyploidy can give rise to novel traits in divergent and successful high-ploidy lineages.

Shoot Maturation Strengthens FLS2-Mediated Resistance to Pseudomonas syringae

Temporospatial regulation of immunity components is essential for properly activating plant defense response. Flagellin-sensing 2 (FLS2) is a surface-localized receptor that recognizes bacterial flagellin. The immune function of FLS2 is compromised in early stages of shoot development. However, the underlying mechanism for the age-dependent FLS2 signaling is not clear. Here, we show that the reduced basal immunity of juvenile leaves against Pseudomonas syringae pv. tomato DC3000 is independent of FLS2. The flg22-induced marker gene expression and reactive oxygen species activation were comparable in juvenile and adult stages, but callose deposition was more evident in the adult stage than the juvenile stage. We further demonstrated that microRNA156, a master regulator of plant aging, does not influence the expression of FLS2 and FRK1 (Flg22-induced receptor-like kinase 1) but mildly suppresses callose deposition in juvenile leaves. Our experiments revealed an intrinsic mechanism that regulates the amplitude of FLS2-mediated resistance during aging. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

The Arabidopsis SNARE complex genes regulate the early stages of pollen-stigma interactions

KEY MESSAGE

The VAMP721, VAMP722, SYP121, SYP122 and SNAP33 SNAREs are required in the Arabidopsis stigma for pollen hydration, further supporting a role for vesicle trafficking in the stigma's pollen responses. In the Brassicaceae, the process of accepting compatible pollen is a key step in successful reproduction and highly regulated following interactions between the pollen and the stigma. Central to this is the initiation of secretion in the stigma, which is proposed to provide resources to the pollen for hydration and germination and pollen tube growth. Previously, the eight exocyst subunit genes were shown to be required in the Arabidopsis stigma to support these pollen responses. One of the roles of the exocyst is to tether secretory vesicles at the plasma membrane for membrane fusion by the SNARE complex to enable vesicle cargo release. Here, we investigate the role of Arabidopsis SNARE genes in the stigma for pollen responses. Using a combination of different knockout and knockdown SNARE mutant lines, we show that VAMP721, VAMP722, SYP121, SYP122 and SNAP33 are involved in this process. Significant disruptions in pollen hydration were observed following pollination of wildtype pollen on the mutant SNARE stigmas. Overall, these results place the Arabidopsis SNARE complex as a contributor in the stigma for pollen responses and reaffirm the significance of secretion in the stigma to support the pollen-stigma interactions.

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.

The eINTACT method for studying nuclear changes in host plant cells targeted by bacterial effectors in native infection contexts

Type-III effector proteins are major virulence determinants that most gram-negative bacteria inject into host cells to manipulate cellular processes for infection. Because effector-targeted cells are embedded and underrepresented in infected plant tissues, it is technically challenging to isolate them for focused studies of effector-induced cellular changes. This protocol describes a novel technique, effector-inducible isolation of nuclei tagged in specific cell types (eINTACT), for isolating biotin-labeled nuclei from Arabidopsis plant cells that have received Xanthomonas bacterial effectors by using streptavidin-coated magnetic beads. This protocol is an extension of the existing Nature Protocols Protocol of the INTACT method for the affinity-based purification of nuclei of specific cell types in the context of developmental biology. In a phytopathology scenario, our protocol addresses how to obtain eINTACT transgenic lines and compatible bacterial mutants, verify the eINTACT system and purify nuclei of bacterial effector-recipient cells from infected tissues. Differential analyses of purified nuclei from plants infected by bacteria expressing the effector of interest and those from plants infected by effector-deletion bacterial mutants will reveal the effector-dependent nuclear changes in targeted host cells. Provided that the eINTACT system is available, the infection experiment takes 5 d, and the procedures, from collecting bacteria-infected leaves to obtaining nuclei of effector-targeted cells, can be completed in 4 h. eINTACT is a unique method for isolating high-quality nuclei from bacterial effector-targeted host cells in native infection contexts. This method is adaptable to study the functions of type-III effectors from numerous gram-negative bacteria in host plants that are amenable to transformation.

Epigenomic divergence correlates with sequence polymorphism in Arabidopsis paralogs

Processes affecting rates of sequence polymorphism are fundamental to the evolution of gene duplicates. The relationship between gene activity and sequence polymorphism can influence the likelihood that functionally redundant gene copies are co-maintained in stable evolutionary equilibria vs other outcomes such as neofunctionalization. Here, we investigate genic variation in epigenome-associated polymorphism rates in Arabidopsis thaliana and consider whether these affect the evolution of gene duplicates. We compared the frequency of sequence polymorphism and patterns of genetic differentiation between genes classified by exon methylation patterns: unmethylated (unM), gene-body methylated (gbM), and transposon-like methylated (teM) states, which reflect divergence in gene expression. We found that the frequency of polymorphism was higher in teM (transcriptionally repressed, tissue-specific) genes and lower in gbM (active, constitutively expressed) genes. Comparisons of gene duplicates were largely consistent with genome-wide patterns - gene copies that exhibit teM accumulate more variation, evolve faster, and are in chromatin states associated with reduced DNA repair. This relationship between expression, the epigenome, and polymorphism may lead to the breakdown of equilibrium states that would otherwise maintain genetic redundancies. Epigenome-mediated polymorphism rate variation may facilitate the evolution of novel gene functions in duplicate paralogs maintained over evolutionary time.

Methylome Response to Proteasome Inhibition by Pseudomonas syringae Virulence Factor Syringolin A

DNA methylation is an important epigenetic mark required for proper gene expression and silencing of transposable elements. DNA methylation patterns can be modified by environmental factors such as pathogen infection, in which modification of DNA methylation can be associated with plant resistance. To counter the plant defense pathways, pathogens produce effector molecules, several of which act as proteasome inhibitors. Here, we investigated the effect of proteasome inhibition by the bacterial virulence factor syringolin A (SylA) on genome-wide DNA methylation. We show that SylA treatment results in an increase of DNA methylation at centromeric and pericentromeric regions of Arabidopsis chromosomes. We identify several CHH differentially methylated regions (DMRs) that are enriched in the proximity of transcriptional start sites. SylA treatment does not result in significant changes in small RNA composition. However, significant changes in genome transcriptional activity can be observed, including a strong upregulation of resistance genes that are located on chromosomal arms. We hypothesize that DNA methylation changes could be linked to the upregulation of some atypical members of the de novo DNA methylation pathway, namely AGO3, AGO9, and DRM1. Our data suggests that modification of genome-wide DNA methylation resulting from an inhibition of the proteasome by bacterial effectors could be part of an epi-genomic arms race against pathogens. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.