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

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.

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.

Localization of the MTP4 transporter to trans-Golgi network in pollen tubes of Arabidopsis thaliana

Zinc (Zn) is an essential element for plants. Numerous proteins in different cellular compartments require Zn for their structure and function. Zn can be toxic when it accumulates in high levels in the cytoplasm. Therefore, Zn homeostasis at tissue, cell, and organelle levels is vital for plant growth. A part of the metal tolerance protein (MTP) / Cation Diffusion Facilitator (CDF) transporters functions as Zn transporters, exporting Zn from the cytosol to various membrane compartments. In Arabidopsis thaliana, MTP1, MTP2, MTP3, MTP4, MTP5, and MTP12 are classified as Zn transporters (Zn-CDF). In this study, we systematically analyzed the localization of GFP-fused Zn-CDFs in the leaf epidermal cells of Nicotiana benthamiana. As previously reported, MTP1 and MTP3 were localized to tonoplast, MTP2 to endoplasmic reticulum, and MTP5 to Golgi. In addition, we identified the localization of MTP4 to trans-Golgi Network (TGN). Since MTP4 is specifically expressed in pollen, we analyzed the localization of MTP4-GFP in the Arabidopsis pollen tubes and confirmed that it is in the TGN. We also showed the Zn transport capability of MTP4 in yeast cells. We then analyzed the phenotype of an mtp4 T-DNA insertion mutant under both limited and excess Zn conditions. We found that their growth and fertility were not largely different from the wild-type. Our study has paved the way for investigating the possible roles of MTP4 in metallating proteins in the secretory pathway or in exporting excess Zn through exocytosis. In addition, our system of GFP-fused MTPs will help study the mechanisms for targeting transporters to specific membrane compartments.

Impact of protein domains on the MEL2 granule, a cytoplasmic ribonucleoprotein complex maintaining faithful meiosis progression in rice

Cytoplasmic ribonucleoprotein (RNP) granules are membraneless structures composed of various RNAs and proteins that play important roles in post-transcriptional regulation. While RNP granules are known to regulate the meiotic entry in some organisms, little is known about their roles in plants. In this study, we observed the cytoplasmic granular structures of rice RNA-binding protein MEIOSIS ARRESTED AT LEPTOTENE2 (MEL2), which contributes to the control of meiotic entry timing, in leaf protoplasts and spore mother cells. We performed colocalization analysis with known cytoplasmic RNP factors, and domain deletion analysis to assess their impact on granule formation and meiosis progression. Conservation of MEL2 domains across plant species was also explored. Our results indicated that MEL2 granules colocalized with processing body and stress granule factors. The maintenance of granule properties modulated by LOTUS domain and the intrinsically disordered region (IDR) is essential for proper MEL2 function in meiosis progression. MEL2-like proteins widely found in plant kingdom conserved LOTUS domain followed by the IDR despite their diverse domain structures, suggesting the functional conservation of these domains among plant species. This study highlights the role of MEL2 granule dynamics and its impact on meiotic transition and progression.

Identification of QTLs associated with yield-related traits and superior genotype prediction using recombinant inbred line population in tobacco

Tobacco is an economically significant industrial crop and model plant for genetic research, yet little is known about its genetic architecture. Quantitative trait loci (QTL) analysis was performed for six agronomic traits on an F_7 population of 341 genotypes, parents, and F1 plants using 1974 SSR markers across two environments. 31 QTLs contributing single-locus additive effects on 13 linkage groups (LGs) and 6 QTL pairs contributing epistatic effects on 6 LGs, were detected by the QTLNetwork 2.0 which was developed for the mixed-linear-model-based composite interval mapping (MCIM). Notably, 5 QTLs and 1 epistatic QTL pair were found to have pleiotropic effects on some genetically related traits. Moreover, the Broad sense heritability of the detected QTLs ranged from 1.05% to 43.33%, while genotype-by-environment interaction heritability spanned from 27.09% to 56.25%. Based on the results of QTL mapping, the potential superior lines for all or specific environments were designed and evaluated. Five major QTLs were finely dissected based on the tobacco reference genome of K326, and 31 candidate genes were predicted. This study offered new insights into the complicated genetic architecture and QTL resources for efficient breeding design for genetic improvement of agronomic traits in tobacco.

The role of D3-type cyclins is related to cytokinin and the bHLH transcription factor SPATULA in Arabidopsis gynoecium development

MAIN CONCLUSION

We studied the D3-type cyclin function during gynoecium development in Arabidopsis and how they are related to the hormone cytokinin and the transcription factor SPATULA. Growth throughout the life of plants is sustained by cell division and differentiation processes in meristematic tissues. In Arabidopsis, gynoecium development implies a multiphasic process where the tissues required for pollination, fertilization, and seed development form. The Carpel Margin Meristem (CMM) is a mass of undifferentiated cells that gives rise to the gynoecium internal tissues, such as septum, ovules, placenta, funiculus, transmitting tract, style, and stigma. Different genetic and hormonal factors, including cytokinin, control the CMM function. Cytokinin regulates the cell cycle transitions through the activation of cell cycle regulators as cyclin genes. D3-type cyclins are expressed in proliferative tissues, favoring the mitotic cell cycle over the endoreduplication. Though the role of cytokinin in CMM and gynoecium development is highly studied, its specific role in regulating the cell cycle in this tissue remains unclear. Additionally, despite extensive research on the relationship between CYCD3 genes and cytokinin, the regulatory mechanism that connects them remains elusive. Here, we found that D3-type cyclins are expressed in proliferative medial and lateral tissues. Conversely, the depletion of the three CYCD3 genes showed that they are not essential for gynoecium development. However, the addition of exogenous cytokinin showed that they could control the division/differentiation balance in gynoecium internal tissues and outgrowths. Finally, we found that SPATULA can be a mechanistic link between cytokinin and the D3-type cyclins. The data suggest that the role of D3-type cyclins in gynoecium development is related to the cytokinin response, and they might be activated by the transcription factor SPATULA.

LAFL Factors in Seed Development and Phase Transitions

Development is a chain reaction in which one event leads to another until the completion of a life cycle. Phase transitions are milestone events in the cycle of life. LEAFY COTYLEDON1 (LEC1), ABA INSENSITIVE3 (ABI3), FUSCA3 (FUS3), and LEC2 proteins, collectively known as LAFL, are master transcription factors (TFs) regulating seed and other developmental processes. Since the initial characterization of the LAFL genes, more than three decades of active research has generated tremendous amounts of knowledge about these TFs, whose roles in seed development and germination have been comprehensively reviewed. Recent advances in cell biology with genetic and genomic tools have allowed the characterization of the LAFL regulatory networks in previously challenging tissues at a higher throughput and resolution in reference species and crops. In this review, we provide a holistic perspective by integrating advances at the epigenetic, transcriptional, posttranscriptional, and protein levels to exemplify the spatiotemporal regulation of the LAFL networks in Arabidopsis seed development and phase transitions, and we briefly discuss the evolution of these TF networks.

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.

Comparative transcriptome profiling reveals distinct regulatory responses of secondary defensive metabolism in Datura species (Solanaceae) under plant development and herbivory-mediated stress

Differential expression of genes is key to mediating developmental and stress-related plant responses. Here, we addressed the regulation of plant metabolic responses to biotic stress and the developmental variation of defense-related genes in four species of the genus Datura with variable patterns of metabolite accumulation and development. We combine transcriptome profiling with phylogenomic techniques to analyze gene expression and coexpression in plants subjected to damage by a specialist folivore insect. We found (1) common overall gene expression in species of similar chemical profiles, (2) species-specific responses of proteins involved in specialized metabolism, characterized by constant levels of gene expression coupled with transcriptional rearrangement, and (3) induction of transcriptional rearrangement of major terpene and tropane alkaloid genes upon herbivory. Our results indicate differential modulation of terpene and tropane metabolism linked to jasmonate signaling and specific transcription factors to regulate developmental variation and stress programs, and suggest plastic adaptive responses to cope with herbivory. The transcriptional profiles of specialized metabolism shown here reveal complex genetic control of plant metabolism and contribute to understanding the molecular basis of adaptations and the physiological variation of significant ecological traits.

Sporophytic control of tapetal development and pollen fertility by a mitogen-activated protein kinase cascade in rice

Tapetum, the innermost layer of the anther wall, provides essential nutrients and materials for pollen development. Timely degradation of anther tapetal cells is a prerequisite for normal pollen development in flowering plants. Tapetal cells facilitate male gametogenesis by providing cellular contents after highly coordinated programmed cell death (PCD). Tapetal development is regulated by a transcriptional network. However, the signaling pathway(s) involved in this process are poorly understood. In this study, we report that a mitogen-activated protein kinase (MAPK) cascade composed of OsYDA1/OsYDA2-OsMKK4-OsMPK6 plays an important role in tapetal development and male gametophyte fertility. Loss of function of this MAPK cascade leads to anther indehiscence, enlarged tapetum, and aborted pollen grains. Tapetal cells in osmkk4 and osmpk6 mutants exhibit an increased presence of lipid body-like structures within the cytoplasm, which is accompanied by a delayed occurrence of PCD. Expression of a constitutively active version of OsMPK6 (CA-OsMPK6) can rescue the pollen defects in osmkk4 mutants, confirming that OsMPK6 functions downstream of OsMKK4 in this pathway. Genetic crosses also demonstrated that the MAPK cascade sporophyticly regulates pollen development. Our study reveals a novel function of rice MAPK cascade in plant male reproductive biology.

Investigating a role for PUB17 and PUB16 in the self-incompatibility signaling pathway in transgenic Arabidopsis thaliana

In Brassicaceae self-incompatibility (SI), self-pollen rejection is initiated by the S-haplotype specific interactions between the pollen S cysteine-rich/S-locus protein 11 (SCR/SP11) ligands and the stigma S receptor kinases (SRK). In Brassica SI, a member of the Plant U-Box (PUB) E3 ubiquitin ligases, ARM-repeat containing 1 (ARC1), is then activated by SRK in this stigma and cellular events downstream of this cause SI pollen rejection by inhibiting pollen hydration and pollen tube growth. During the transition to selfing, Arabidopsis thaliana lost the SI components, SCR, SRK, and ARC1. However, this trait can be reintroduced into A. thaliana by adding back functional copies of these genes from closely related SI species. Both SCR and SRK are required for this, though the degree of SI pollen rejection varies between A. thaliana accessions, and ARC1 is not always needed to produce a strong SI response. For the A. thaliana C24 accession, only transforming with Arabidopsis lyrata SCR and SRK confers a strong SI trait (SI-C24), and so here, we investigated if ARC1-related PUBs were involved in the SI pathway in the transgenic A. thaliana SI-C24 line. Two close ARC1 homologs, PUB17 and PUB16, were selected, and (Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) technology was used to generate pub17 and pub16 mutations in the C24 accession. These mutants were then crossed into the transgenic A. thaliana SI-C24 line and their potential impact on SI pollen rejection was investigated. Overall, we did not observe any significant differences in SI responses to implicate PUB17 and PUB16 functioning in the transgenic A. thaliana SI-C24 stigma to reject SI pollen.

Debottlenecking the L-DOPA 4,5-dioxygenase step with enhanced tyrosine supply boosts betalain production in Nicotiana benthamiana

Synthetic biology provides emerging tools to produce valuable compounds in plant hosts as sustainable chemical production platforms. However, little is known about how supply and utilization of precursors is coordinated at the interface of plant primary and specialized metabolism, limiting our ability to efficiently produce high levels of target specialized metabolites in plants. L-Tyrosine is an aromatic amino acid precursor of diverse plant natural products including betalain pigments, which are used as the major natural food red colorants and more recently a visual marker for plant transformation. Here, we studied the impact of enhanced L-tyrosine supply on the production of betalain pigments by expressing arogenate dehydrogenase (TyrA) from table beet (Beta vulgaris, BvTyrAα), which has relaxed feedback inhibition by L-tyrosine. Unexpectedly, betalain levels were reduced when BvTyrAα was coexpressed with the betalain pathway genes in Nicotiana benthamiana leaves; L-tyrosine and 3,4-dihydroxy-L-phenylalanine (L-DOPA) levels were drastically elevated but not efficiently converted to betalains. An additional expression of L-DOPA 4,5-dioxygenase (DODA), but not CYP76AD1 or cyclo-DOPA 5-O-glucosyltransferase, together with BvTyrAα and the betalain pathway, drastically enhanced betalain production, indicating that DODA is a major rate-limiting step of betalain biosynthesis in this system. Learning from this initial test and further debottlenecking the DODA step maximized betalain yield to an equivalent or higher level than that in table beet. Our data suggest that balancing between enhanced supply ("push") and effective utilization ("pull") of precursor by alleviating a bottleneck step is critical in successful plant synthetic biology to produce high levels of target compounds.

Evolution of Whirly1 in the angiosperms: sequence, splicing, and expression in a clade of early transitional mycoheterotrophic orchids

The plastid-targeted transcription factor Whirly1 (WHY1) has been implicated in chloroplast biogenesis, plastid genome stability, and fungal defense response, which together represent characteristics of interest for the study of autotrophic losses across the angiosperms. While gene loss in the plastid and nuclear genomes has been well studied in mycoheterotrophic plants, the evolution of the molecular mechanisms impacting genome stability is completely unknown. Here, we characterize the evolution of WHY1 in four early transitional mycoheterotrophic orchid species in the genus Corallorhiza by synthesizing the results of phylogenetic, transcriptomic, and comparative genomic analyses with WHY1 genomic sequences sampled from 21 orders of angiosperms. We found an increased number of non-canonical WHY1 isoforms assembled from all but the greenest Corallorhiza species, including intron retention in some isoforms. Within Corallorhiza, phylotranscriptomic analyses revealed the presence of tissue-specific differential expression of WHY1 in only the most photosynthetically capable species and a coincident increase in the number of non-canonical WHY1 isoforms assembled from fully mycoheterotrophic species. Gene- and codon-level tests of WHY1 selective regimes did not infer significant signal of either relaxed selection or episodic diversifying selection in Corallorhiza but did so for relaxed selection in the late-stage full mycoheterotrophic orchids Epipogium aphyllum and Gastrodia elata. Additionally, nucleotide substitutions that most likely impact the function of WHY1, such as nonsense mutations, were only observed in late-stage mycoheterotrophs. We propose that our findings suggest that splicing and expression changes may precede the selective shifts we inferred for late-stage mycoheterotrophic species, which therefore does not support a primary role for WHY1 in the transition to mycoheterotrophy in the Orchidaceae. Taken together, this study provides the most comprehensive view of WHY1 evolution across the angiosperms to date.

The cacao gene atlas: a transcriptome developmental atlas reveals highly tissue-specific and dynamically-regulated gene networks in Theobroma cacao L

BACKGROUND

Theobroma cacao, the cocoa tree, is a tropical crop grown for its highly valuable cocoa solids and fat which are the basis of a 200-billion-dollar annual chocolate industry. However, the long generation time and difficulties associated with breeding a tropical tree crop have limited the progress of breeders to develop high-yielding disease-resistant varieties. Development of marker-assisted breeding methods for cacao requires discovery of genomic regions and specific alleles of genes encoding important traits of interest. To accelerate gene discovery, we developed a gene atlas composed of a large dataset of replicated transcriptomes with the long-term goal of progressing breeding towards developing high-yielding elite varieties of cacao.

RESULTS

We describe the creation of the Cacao Transcriptome Atlas, its global characterization and define sets of genes co-regulated in highly organ- and temporally-specific manners. RNAs were extracted and transcriptomes sequenced from 123 different tissues and stages of development representing major organs and developmental stages of the cacao lifecycle. In addition, several experimental treatments and time courses were performed to measure gene expression in tissues responding to biotic and abiotic stressors. Samples were collected in replicates (3-5) to enable statistical analysis of gene expression levels for a total of 390 transcriptomes. To promote wide use of these data, all raw sequencing data, expression read mapping matrices, scripts, and other information used to create the resource are freely available online. We verified our atlas by analyzing the expression of genes with known functions and expression patterns in Arabidopsis (ACT7, LEA19, AGL16, TIP13, LHY, MYB2) and found their expression profiles to be generally similar between both species. We also successfully identified tissue-specific genes at two thresholds in many tissue types represented and a set of genes highly conserved across all tissues.

CONCLUSION

The Cacao Gene Atlas consists of a gene expression browser with graphical user interface and open access to raw sequencing data files as well as the unnormalized and CPM normalized read count data mapped to several cacao genomes. The gene atlas is a publicly available resource to allow rapid mining of cacao gene expression profiles. We hope this resource will be used to help accelerate the discovery of important genes for key cacao traits such as disease resistance and contribute to the breeding of elite varieties to help farmers increase yields.

The Arabidopsis AtSWEET13 transporter discriminates sugars by selective facial and positional substrate recognition

Transporters are targeted by endogenous metabolites and exogenous molecules to reach cellular destinations, but it is generally not understood how different substrate classes exploit the same transporter's mechanism. Any disclosure of plasticity in transporter mechanism when treated with different substrates becomes critical for developing general selectivity principles in membrane transport catalysis. Using extensive molecular dynamics simulations with an enhanced sampling approach, we select the Arabidopsis sugar transporter AtSWEET13 as a model system to identify the basis for glucose versus sucrose molecular recognition and transport. Here we find that AtSWEET13 chemical selectivity originates from a conserved substrate facial selectivity demonstrated when committing alternate access, despite mono-/di-saccharides experiencing differing degrees of conformational and positional freedom throughout other stages of transport. However, substrate interactions with structural hallmarks associated with known functional annotations can help reinforce selective preferences in molecular transport.

Interplay between phosphorylation and oligomerization tunes the conformational ensemble of SWEET transporters

SWEET sugar transporters are desirable biotechnological targets for improving plant growth. One engineering strategy includes modulating how SWEET transporters are regulated. Phosphorylation and oligomerization have been shown to positively regulate SWEET function, leading to increased sugar transport activity. However, constitutive phosphorylation may not be beneficial to plant health under basal conditions. Structural and mechanistic understanding of the interplay between phosphorylation and oligomerization in functional regulation of SWEETs remains limited. Using extensive molecular dynamics simulations coupled with Markov state models, we demonstrate the thermodynamic and kinetic effects of SWEET phosphorylation and oligomerization using OsSWEET2b as a model. We report that the beneficial effects of these SWEET regulatory mechanisms bias outward-facing states and improved extracellular gating, which complement published experimental findings. Our results offer molecular insights to SWEET regulation and may guide engineering strategies throughout the SWEET transport family.

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.

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.

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.