Microbial-based approaches offer a promising strategy to decrease the use of chemical fertilizers in agriculture. Among them, arbuscular mycorrhizal fungi (AMF), which extend root surface area and enhance phosphate uptake, and phosphate-solubilizing bacteria (PSB) are particularly relevant. However, their effectiveness depends strongly on plant genetic diversity. To identify genetic markers underlying plant responses to beneficial soil microbes, we studied a panel of 128 fully sequenced Lactuca sativa varieties under controlled phosphate-starvation conditions and treated with AMF and PSB. Lettuce genetic variation showed a strong effect on physiological and morphological responses to microbial inoculation. Genome-wide association studies identified specific genomic regions associated with changes in leaf phosphate content and shoot biomass following treatment. Beyond genetic factors, we observed shifts in fungal β-diversity and increased bacterial α-diversity associated with phenotypic variation. We also identified 44 amplicon sequence variants associated with agriculturally relevant traits. Among these, six bacterial strains were experimentally validated through in vitro and pot experiments for their effects on leaf phosphate concentration and shoot biomass. Overall, we highlighted key genetic, microbial, and physiological mechanisms that may enhance microbial treatments for improved plant phosphate management in lettuce.

Plasmodiophora brassicae, the causal agent of clubroot disease, actively manipulates host hormone pathways to promote infection. Previous studies identified PbGH3 as a putative GH3-like effector with in vitro auxin-conjugating activity, although its biological function in planta was investigated but remained unresolved. Here, we studied the role of PbGH3 during host colonization using PbGH3-overexpressing Arabidopsis thaliana and Brassica napus lines. PbGH3 overexpression induced conserved developmental phenotypes in both species, including epinastic leaves, reduced apical dominance, enhanced lateral branching, and increased root hair formation, consistent with altered hormone balance. However, hormone profiling and exogenous auxin assays did not support a role for PbGH3 as a canonical auxin-conjugating enzyme in planta. Instead, structural analyses revealed a strong similarity between PbGH3 and the Arabidopsis clade III GH3 proteins GH3.12/PBS3 and GH3.7, which are associated with salicylic acid (SA) metabolism. Consistent with this observation, loss of GH3.12/PBS3 increased susceptibility to clubroot, whereas expression of PbGH3 in the gh3.12 mutant background partially restored SA accumulation and reduced disease susceptibility. PbGH3 expression peaked during the early stages of infection, and overexpressing lines displayed enhanced root hair colonization by P. brassicae, suggesting a role during primary infection rather than later gall development. Together, our results support a model in which PbGH3 modulates SA-associated hormonal crosstalk and root morphology that could be harnessed by the clubroot pathogen to ensure host colonization. These findings provide new insight into how P. brassicae manipulates host hormone homeostasis during infection. [Formula: see text] Copyright © 2026 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Thylakoid membranes are indispensable for oxygenic photosynthesis, yet the mechanisms that protect these membranes from photooxidative damage remain poorly understood. By screening poorly characterized proteins induced during the chloroplast unfolded protein response, we identify VIA1 as an essential factor for preserving thylakoid integrity under high light in the model green alga Chlamydomonas reinhardtii. Loss of VIA1 causes hypersensitivity to photo-oxidative stress and rapid thylakoid swelling. VIA1 localizes to thylakoid membranes and directly binds Vesicle-Inducing Protein in Plastids 1 (VIPP1), an ESCRT-III-like protein essential for thylakoid biogenesis and remodeling. Structure-guided mutagenesis shows that this interaction is required for VIA1 function and is mediated by a winged-helix domain interface reminiscent of ESCRT-II/ESCRT-III binding mode. VIA1 orthologs from cyanobacteria and land plants rescue the Chlamydomonas via1 mutant phenotype, and disruption of VIA1 in Synechocystis sp. PCC 6803 impairs growth, especially under light stress. Together, these findings establish VIA1 as an evolutionarily conserved protein that contributes to thylakoid membrane homeostasis via its interaction with VIPP1.

MicroRNAs (miRNAs) are small non-coding RNA molecules essential for growth and development in eukaryotes. In plants, the master gene DICER-LIKE 1 (DCL1) catalyzes the biogenesis of miRNAs by processing double-stranded precursors that give rise to mature miRNAs. We sought to understand the function and evolution of microRNAs using Marchantia polymorpha, a model bryophyte that allows comparative approaches to infer characteristics of the ancestral land plant. We functionally characterized loss-of-function mutants of MpDCL1 a generated by means of CRISPR-Cas9-mediated genomic edition and the miR166/CLASS III HD-ZIP regulatory circuit in Marchantia polymorpha. We report that MpDCL1a is required for the biogenesis of miRNAs and uncovered a central role for miR166/Homeodomain Zipper Class III-regulated auxin synthesis in the specification of cell identity, patterning, meristem function, laminar expansion, and the development of the body in the last common ancestor of the bryophytes and vascular plants. Our findings indicate that DCL1, Class III HD-ZIP, miR166, and auxin functioned in the development of the body of the last common ancestor of extant land plants and provide a novel working framework to interrogate the basic principles of cell specification and patterning in plants.

Microbial host-parasite interactions shape global biogeochemical cycles, yet their cellular and evolutionary mechanisms remain poorly understood. Here we developed a model pathosystem to explore the interaction between the eukaryotic microparasite, Pirsonia diadema, and its host, the bloom-forming diatom, Coscinodiscus. Culture conditions recapitulated Pirsonia's rapid life cycle in vitro, and live-cell imaging enabled the quantitative analysis of infection dynamics and host mortality. Genome sequencing and time-resolved dual RNA sequencing of Pirsonia revealed unique genetic innovations relative to their oomycete relatives and an upregulation of an expanded multi-gene family of integrin-like proteins in Pirsonia zoospores during infection. Targeted drug perturbations combined with ultrastructure expansion microscopy confirm stage specific roles for Pirsonia's cytoskeleton, including actin-dependent host attachment and formation of intracellular parasite-induced feeding structures and microtubule-associated processes in zoospore differentiation. Together, these findings reveal mechanisms underpinning Pirsonia infection and provide insights into the evolution of strategies that underpin parasitic shunts in marine trophic networks.

Cytosine methylation at CG sites is proposed to be maintained during DNA replication through a semiconservative mechanism. In mammals, this maintenance involves replication-coupled and replication-uncoupled steps; however, the underlying mechanisms and their possible conservation across eukaryotes remain unclear. Using live-cell sensors developed in this study to track CG methylation (mCG) dynamics, we uncovered that mCG maintenance extends far beyond DNA replication in the model plant Arabidopsis, with the major CG methyltransferase 1 (MET1) remaining associated with chromatin from S phase till G2 phase. Genetic dissection demonstrated that MET1, like its mammalian homolog DNA methyltransferase 1 (DNMT1), methylates nucleosome-associated DNA. Single-molecule duplex sequencing revealed that MET1/DNMT1-dependent methylation follows the rotational accessibility of nucleosomal DNA and is only partly templated by hemimethylated sites. This nucleosome-constrained de novo activity not only promotes mCG recovery after accidental loss but also shapes short- and long-term DNA methylation variation in plants and mammals.

We discovered a previously uncharacterized moss-specific protein, Rosetta NATD-MAPK 1 (RAK1) in Physcomitrium patens, which uniquely integrates MAP kinase (MAPK)-dependent signaling with N-acetyltransferase activity. Through phenotypical and biochemical analyses, we characterized RAK1 function in the regulation of the 2D-to-3D growth transition. We identified differential acetylation events associated with metabolic reprogramming in rak1 mutants and demonstrated that RAK1 has acetyltransferase activity enhanced by MAPK domain activation. Consistently, RAK1 was found to interact with proteins involved in metabolic processes. Collectively, this study uncovers a previously unknown multidomain protein and provides mechanistic insights into the interplay of post-translational modifications during developmental reprogramming.

We present a novel diversity-focused haplotype map (HapMap) that characterizes over 64.5 million maize (Zea mays ssp. mays) single nucleotide polymorphisms (SNPs) genotyped across 818 individuals from diverse backgrounds. This HapMap aims to balance the variation obtained from domesticated landraces and inbred lines, outgroup Zea spp. and more distant Tripsacum spp. in order to minimize ascertainment bias for diversity studies. Included individuals derive from public data from various experimental setups and coverages, which is challenging for standard SNP callers to accommodate. We provide evidence of coverage biases associated with standard callers that influence resulting variation and introduce a novel approach called Unified Multi-Caller Ensemble (UME), which enhances variant calling accuracy in low-coverage and mixed-coverage genomic datasets. UME corrects for coverage bias resulting from inter-sample coverage heterogeneity by leveraging evidence from variant callers with orthogonal strategies, re-calibrating the error probabilities across callers to minimize the impact of error biases inherent to a given caller. It outperforms individual strategies and excels in de novo variant calling, taking advantage of instances of higher depth reads, even in low coverage individuals, while preserving biologically informative variant relationships across coverage levels. An important feature of UME is the independence from population allele frequencies in the discovery panel, thus avoiding ascertainment bias resulting from unbalanced input genetic diversity. Discovered variants are less affected by ascertainment bias because no population filtering is used, and the full diversity of SNPs is retained in the final variant call set to maximize the utility of the dataset for production calling newly sequenced samples. We present a strategy for filtering the recalibrated error profiles that relies on maximizing demographic signals to retain genetic relationships within the population while reducing sequencing error. After the variant discovery phase, we employ the UME production stage, which enriches genotype calling across all coverage levels, benefiting low-coverage samples. Error introduced in this process is removed through subsequent filtering. Using this approach, we generated a coverage bias-controlled maize HapMap database, providing a comprehensive representation of maize accessions and emphasizing landrace diversity. This diverse panel of domesticated maize and outgroups from across the Americas enables accurate genotyping in low-coverage samples while offering crucial context for interpreting diversity, particularly for natural diversity and paleogenomic analyses.

Heterochromatin formation is pivotal in many eukaryotes with repetitive sequences, such as transposable elements (TEs). However, in plants, where the known de novo DNA methylation mechanism (RdDM) targets euchromatin, how heterochromatin is formed in a region-specific manner remains unclear. We previously reported an RdDM-independent de novo establishment of H3K9me and non-CpG methylation, both of which localize in heterochromatin. Here we show that the mutually exclusive histone H2A variants, H2A.W and H2A.Z, function as guides to initiate heterochromatin formation; H2A.W and H2A.Z promotes and inhibits heterochromatin establishment, respectively, especially in chromosomal arm regions with dispersed TEs. In contrast, pericentromeric TEs demonstrate autonomous heterochromatin formation, less dependently on these H2A variants. Furthermore, H2A.Z protects protein-coding genes from ectopic heterochromatin formation, possibly by preventing its spreading. We propose that the genome indexing mechanism driven by H2A variants, as well as the autonomous formation of pericentromeric heterochromatin, shapes proper epigenomic patterns in Arabidopsis.

Although research on the nucleus has historically focused on proteins and nucleic acids, the crucial roles of lipids in nuclear biology have remained underappreciated. Lipids are not merely structural components but are now understood to be active regulators of diverse processes vital to nuclear homeostasis and function, including nuclear pore complex assembly, RNA processing, the DNA damage response and gene expression control. Furthermore, the nucleus is emerging as a critical hub for both lipid metabolism and mechanosensing, contributing to cellular lipid and mechanical homeostasis. This Perspective synthesizes the current understanding of lipids in nuclear function and genome integrity, highlighting the key findings and open questions that were the focus of The Company of Biologists Workshop 'Roles of Lipids in Nuclear Homeostasis and Genome Stability'.

Ultraviolet-B radiation (UV-B) poses a major challenge to all forms of plant life. The liverwort Marchantia (Marchantia polymorpha) serves as a key model organism for studying signaling pathways and inferring their evolution throughout the green lineage. Marchantia expresses key components of UV-B signaling, including the photoreceptor UV RESISTANCE LOCUS 8 (MpUVR8), the WD40-repeat protein REPRESSOR OF UV-B PHOTOMORPHOGENESIS (MpRUP), the E3 ubiquitin ligase complex CONSTITUTIVELY PHOTOMORPHOGENIC 1/SUPPRESSOR OF phyA-105 (MpCOP1/MpSPA), and the transcriptional regulator ELONGATED HYPOCOTYL 5 (MpHY5). Here, we show that whereas the core UVR8 photocycle is well conserved, regulatory roles of SPA and HY5 diverged during plant evolution. Unlike the case in Arabidopsis (Arabidopsis thaliana) and in contrast to the strong Mpcop1 mutant phenotype, Mpspa mutants developed only a very weak constitutive photomorphogenesis phenotype, indicating that COP1 function is much more independent of SPA in Marchantia than in Arabidopsis. Moreover, in contrast to Arabidopsis SPAs, Mpspa showed enhanced UV-B acclimation and UV stress tolerance, indicating that MpSPA is a negative regulator of MpUVR8 signaling. Similar to Arabidopsis HY5/HYH, MpHY5 functioned antagonistically to MpCOP1, but its role in UV-B-mediated gene expression changes was more limited. Our findings demonstrate that although core components of UV-B signaling existed in the last common ancestor of extant land plants, regulatory interactions have diversified in different lineages since their divergence more than 400 million years ago.

RNA viruses co-opt host endomembranes to form replication complexes, often triggering cellular stress and immune responses. Here, we show that Arabidopsis thaliana activates selective autophagy to respond to viruses targeting mitochondria, chloroplasts, and the endoplasmic reticulum. Rather than degrading viral components, autophagy selectively removes the immune regulator Enhanced Disease Susceptibility 1 (EDS1) to prevent cell death. This targeted mechanism is mediated by oligomeric metabolic enzymes that moonlight as selective autophagy receptors, linking organelle stress to immune homeostasis. Our findings establish selective autophagy as an essential immune rheostat that fine-tunes defense responses and safeguards cellular integrity to promote host survival during viral infections.

The unicellular red algae, Cyanidiophyceae, that diverged early during Archaeplastida (algal and plant) evolution, occupy a variety of extreme habitats that are inhospitable for most other eukaryotes. With the use of modern genomics and genetics methods, Cyanidiophyceae show a remarkable taxonomic diversity, share haplodiplophasic life cycles, and are engaged in complex trophic interactions with microbes that occupy geothermal niches. Amenable to molecular engineering, Cyanidiophyceae are excellent models for understanding evolutionary mechanisms that underpin their extremophilic lifestyles. Their unique growth conditions make these choice red algae of high interest for biotechnological exploitation in environments unsuitable for crops.

Autophagy is a conserved intracellular catabolic process, critical for plant stress tolerance. Upon their delivery in the vacuole, how autophagic bodies containing cargo are hydrolyzed to warrant autophagy degradation remains unclear in multicellular organisms. Here, we found that two Arabidopsis phospholipases, LCAT4 and LCAT3, traffic to the vacuolar lumen and converge on autophagic bodies through fundamentally different routes. While LCAT4 directly binds ATG8 and uses autophagy as a transport system to reach the vacuole prepackaged within autophagosomes, LCAT3 traffics to the lytic compartment independently of autophagosome formation. Knocking out both genes causes an accumulation of autophagic bodies accompanied with a reduction in autophagy degradation. In vivo reconstitution demonstrated that LCAT3 can hydrolyse the membrane of autophagic bodies, enabling the activity of LCAT4 to enhance this process. Together, this work sheds light on the vacuolar stages of autophagy, showing that plants have evolved a multi-component pathway for the efficient disruption of autophagosomal membranes as a critical step for the completion of the autophagy pathway.

Autophagy is a vital cellular quality control pathway that maintains cellular homoeostasis under changing environments. While the organismal phenotypes of autophagy-deficient plants under stress are well characterized, the contribution of cell-type-specific autophagy responses to whole-plant homoeostasis remains poorly understood. Here we show that root-hair-forming cells (trichoblasts) in Arabidopsis thaliana exhibit higher autophagic flux than adjacent non-hair cells (atrichoblasts). This differential autophagy is genetically linked to cell fate determination during early development. Mutants with disrupted trichoblast or atrichoblast identity lose the autophagy distinction between these cell types. Functional analyses reveal that elevated autophagy in trichoblasts is essential for sodium ion sequestration in vacuoles-a key mechanism for salt stress tolerance. Disrupting autophagy specifically in trichoblasts impairs sodium accumulation and reduces plant survival under salt stress. Conversely, cell-type-specific complementation restores both sodium sequestration and stress tolerance. Our findings uncover a cell-type-specific autophagy program in root hairs and demonstrate how developmental cues shape autophagy to enhance stress resilience. This work establishes a direct link between cell identity, autophagy and environmental adaptation in plants.

Root hairs are tubular tip-growing extensions of root epidermal cells that enhance root surface area for water and nutrient acquisition. While the mechanisms governing root hair fate, polarity and tip growth are well understood, the regulation of root hair longevity remains largely unknown. Here we show that root hair cells employ high levels of autophagy to maximize their lifespan in Arabidopsis thaliana. Loss-of-function mutations in the genes encoding the autophagy regulators ATG2, ATG5 or ATG7 induce a premature, cell-autonomous cell death program. This cell death is activated via a gene regulatory network promoted by the NAC transcription factors ANAC046 and ANAC087. Our findings uncover an antagonistic relationship between autophagy and senescence-controlled programmed cell death in root hair lifespan regulation, with potential implications for optimizing plant nutrient and water uptake in crop breeding.

Pathogens use sophisticated strategies to modulate host protein homeostasis by targeting proteolytic pathways, but their impact on protein synthesis remains elusive. We report that pathogenic bacteria Pseudomonas syringae (Pst) targets ribonucleoprotein condensates, known as processing bodies (P-bodies), to attenuate host translation through two effectors with liquid-like properties. We uncovered a previously unknown link that Pst-mediated repression of the endoplasmic reticulum stress response is required for P-body assembly. Furthermore, we identify a functional link between P-bodies and autophagy, demonstrating that autophagic clearance of P-bodies is crucial for maintaining the balance between translationally active and inactive messenger RNAs. Together, our findings provide insights on how host translation is attenuated by bacteria to dampen plant immunity and uncover unknown connections between ER stress responses and autophagy with P-body dynamics.

Histone variants are key regulators of chromatin function. The H3 variants H3.1 and H3.3 evolved independently in animals and plants and differ at amino acid position 31, where H3.1 contains alanine (A), whereas H3.3 carries serine (S) in animals or threonine (T) in plants. Although S and T can both be phosphorylated, the biological significance of plants having selectively adopted T over S remains unclear. Here, we report that H3.3T31 plays a critical role in plant development and stress responses by promoting H3K36me3 on H3.3. T31 prevents plant-specific H3K27 methyltransferases ATXR5 and ATXR6 from depositing H3K27me1, which otherwise inhibits the H3K36 methyltransferase EFS. The substitution of H3.3T31 with S or A increases ATXR5/6 activity and elevates H3K27me1, leading to reduced H3K36me3. Together, these findings suggest co-selection of the plant-specific H3.3T31 residue and ATXR5/6 to ensure the preferential accumulation of H3K27me1 on H3.1 and H3K36me3 on H3.3, thereby supporting chromatin function in plants.

Processing bodies (P-bodies) are conserved ribonucleoprotein granules central to RNA metabolism across eukaryotes. Although the mechanisms underlying their assembly are well understood, the pathways governing their selective turnover remain unclear. Here, we identify the conserved decapping proteins Enhancer of mRNA decapping 4 (EDC4) and decapping protein 1 (DCP1) as a selective autophagy receptor pair responsible for P-body turnover in the model plant Marchantia polymorpha. MpEDC4 engages ATG8 via a canonical ATG8-interacting motif, while MpDCP1 contains a previously unrecognized reverse ATG8-interacting motif within its intrinsically disordered region. Mutations disrupting these motifs impair the autophagic degradation of P-bodies, demonstrating a cooperative receptor mechanism. Notably, this autophagic function is lineage-specific, as orthologs in Arabidopsis and humans lack ATG8-binding capacity. Strikingly, the heterologous expression of MpEDC4 in human cells promotes the degradation of α-synuclein, a protein linked to Parkinson's disease etiology. Our findings uncover an evolutionary innovation that links RNA metabolism to selective autophagy and open avenues for the cross-kingdom engineering of targeted protein degradation pathways.

The long-standing notion that genotypes map to phenotypes through simple one gene-one trait relationships continues to shape both research in the life sciences and public understanding, with implications for policy and funding priorities. Yet this paradigm is increasingly recognized as inadequate for explaining continuous phenotypic variation and the complex genetic architectures of the genotype-phenotype map. Modern genetics emerged from the early 20th-century synthesis of Mendelian and biometric schools of heredity, with R.A. Fisher demonstrating early on how multiple discrete loci could collectively produce continuous variation. Despite this fundamental insight, Mendelism-with its focus on single genes and standardized genetic backgrounds-became the dominant framework, shaping current genetics research and molecular biology as well as science education. The advent of large-scale genomic data has revealed yet again the limitations of this reductionist approach. Evidence from quantitative genetics now shows that most phenotypes arise from complex networks of many interdependent genes and their dynamic responses to environmental perturbations. Here we trace the historical roots of how Mendelian classical genetics departed from the biometric school to create the current predominant paradigm in genetics, despite fundamentally unresolved issues. Moving on from this one-sided paradigm will require systematic development of integrative, evolutionarily grounded experimental approaches that better capture the multigenic and context-dependent nature of inheritance. Achieving such an extended perspective will require methodological innovation, including advances in large-scale (e.g. automated) phenotyping. Dedicated research programs will be necessary to advance a new era of genetic research into the complex mechanisms underlying phenotypic variation.

The growth orientation of the Marchantia polymorpha thallus-a system of dorsiventralized, indeterminate axes-is modulated by light. We show that red and blue light act antagonistically to control thallus growth orientation, with red-light signaling promoting epinasty and blue-light signaling promoting hyponasty. We found that loss-of-function mutations in the blue-light receptor MpPHOT led to epinasty, while loss-of-function mutations in the red-light receptor MpPHY resulted in hyponasty. We hypothesize that these antagonistic activities of blue- and red-light signaling are balanced in white light, resulting in the development of flat thalli. Using time-resolved transcriptomics, we identified genes that were rapidly induced upon light exposure. Among these genes were all six members of the M. polymorpha BBX gene family. Mutants harboring loss-of-function mutations in two of the six MpBBX transcription factors developed defective thallus growth orientation. Mpbbx1 loss-of-function mutants formed hyponastic thalli, while Mpbbx5 loss-of-function mutants developed epinastic thalli. Double mutants Mpbbx1 and Mpbbx5 grew flat, supporting the hypothesis that they function antagonistically. Together, these data indicate that phototropin-mediated blue-light and phytochrome-mediated red-light signaling antagonistically modulate thallus growth orientation and that BBX transcription factors also act antagonistically to regulate thallus flatness.

Clearance of arrested nascent polypeptides resulting from ribosomal stalling is essential for proteostasis. Stalled endoplasmic reticulum (ER)-bound ribosomes are marked by ubiquitin-fold modifier 1 (UFM1) on the large ribosomal subunit protein RPL26, but the precise role of this modification in ribosome-associated quality control (RQC) remains poorly understood. Here, we define the interplay between the UFMylation machinery and the RQC in clearing arrested polypeptides upon ribosome stalling at the ER. Proteomic analysis shows that RQC factors associate with UFMylated ribosomes. Functional assays demonstrate that ribosome rescue factors ZNF598 and ASC-1 recognize and split stalled ribosomes at the ER, a prerequisite for RPL26 UFMylation. The UFM1 E3 ligase complex then binds and UFMylates the post-split 60S-peptidyl-tRNA complex, facilitating access of RQC factors. Depletion of the NEMF/LTN1 complex leads to accumulation of UFMylated ribosomes, whereas impaired UFMylation weakens NEMF/LTN1 binding to ER-stalled ribosomes, supporting a physical link between these pathways. These findings demonstrate that RQC cooperates with the UFMylation machinery to overcome the topological constraints of clearing the arrested polypeptides at the ER.

Field Asymmetric Ion Mobility Spectrometry (FAIMS) is an indispensable tool in single-cell and low-input proteomics. Here, we show that tuning FAIMS resolution via electrode temperature modulation improves sensitivity. We demonstrate that lowering FAIMS resolution broadens the compensation voltage window, thereby increasing ion transmission. This significantly improves ion counts, quantitative precision and peptide identifications by 25-34%. These results were mirrored when analyzing single cell samples.

Maintaining transcriptional fidelity is essential for precise gene regulation and genome stability. Despite this, cryptic antisense transcription, occurring opposite to canonical coding sequences, is a pervasive feature across all domains of life. How such potentially harmful cryptic sites are regulated remains incompletely understood. Here, we show that nucleosome arrays within gene bodies play a key role in suppressing cryptic transcription. Using the fission yeast Schizosaccharomyces pombe as a model, we demonstrate that the CHD-family chromatin remodeler Hrp3 coordinates with the transcription elongation machinery, via the transcriptional regulator Prf1/RTF1, to position nucleosomes at sites of cryptic transcription initiation within gene bodies. In the absence of Hrp3, AT-rich sequences within gene bodies lose nucleosome occupancy, exposing promoter-like sequences that drive cryptic initiation. While cryptic transcription is generally detrimental, we identify a subset of antisense transcripts that encode critical meiotic genes, suggesting that cryptic transcription can also serve as a source of regulatory innovation. These findings define an elongation‑coupled chromatin pathway that preserves transcriptional fidelity and reveal how nucleosome remodeling shapes antisense transcription, cellular homeostasis, and adaptive potential.

Injury-induced disruption of cellular homeostasis leads to accumulation of stress at sites adjacent to a wound. How do these cells mitigate wound-induced stress and restore cellular homeostasis to promote regeneration? To address this question, we examined the role of autophagy-a conserved quality control pathway that recycles defective cellular components to sustain cellular homeostasis-in facilitating wound repair in plants. We demonstrate that transcriptional activation of autophagy-related gene 8 (ATG8) genes is essential for de novo root regeneration, but dispensable for wound induced callus formation from an excised leaf. Plant-specific transcription factors PLETHORA (PLT) activate the transcription of a subset of ATG8 genes and function nonredundantly in this process. Disrupting the PLT-ATG8 regulatory axis severely impairs organelle turnover, resulting in increased intracellular stress and ectopic accumulation of reactive oxygen species (ROS). PLT-ATG8 mediated positioning of optimal ROS levels promotes the expression of stem-cell regulators for successful de novo root regeneration. Altogether, our findings illustrate how plants utilize kingdom-specific developmental regulators, such as PLTs, to activate evolutionarily conserved pathways, effectively managing wound-induced cellular stress and facilitating organ regeneration.

Somatic embryogenesis (SE), the process by which differentiated somatic cells are reprogrammed to form embryos, represents a unique manifestation of plant cell totipotency. Despite its fundamental and applied importance, the molecular mechanisms that initiate embryogenic reprogramming remain largely unknown, mainly because explant tissues are cell-type heterogeneous and contain only a small fraction of SE-competent cells. Here, we applied fluorescence-activated nuclei sorting (FANS) in Arabidopsis thaliana to isolate cells expressing the WOX2 gene, a marker for early embryogenesis. Comparative transcriptomic analysis of WOX2(+)-positive and WOX2(-)-negative nuclei revealed that genes up-regulated in the embryogenic nuclei were strongly enriched in biological processes related to embryo and tissue development, while down-regulated transcripts were linked to primary metabolism, suggesting a transcriptional switch from differentiated to totipotent identity. Several transcription factor genes showed strong induction, including MYB46, ZAT14, MYB98, GRF7, MYR2, and TCP19. Functional analyses using β-estradiol-inducible overexpression and loss-of-function mutants confirmed that these genes modulate embryogenic competence. In particular, MYB46 and ZAT14 emerged as candidate regulators acting downstream of WOX2 to coordinate complementary aspects of cell wall remodeling required for cellular isolation and the acquisition of embryogenic identity. These findings reveal an early, WOX2-centered regulatory network that precedes the activation of canonical SE regulators such as LEC2 and BBM, uncovering a new transcriptional layer of totipotency control.

Duckweeds (Lemnaceae) have excellent potential for fundamental and applied research due to ease of cultivation, small size, and continuous fast clonal growth. However, their usage as model organisms and platforms for biotechnological applications is often limited by the lack of universal genetic manipulation methods necessary for transgene expression, gene editing, and other methods to modify gene expression. To identify suitable strains for genetic manipulation of the giant duckweed, Spirodela polyrhiza, we screened several genotypes for callus induction and regeneration and established genetic transformation. We identified SP162 to be amenable to Agrobacterium-mediated transformation via tissue culture. The procedure is robust and reproducible across laboratories, allowing stable expression of different reporter genes and selectable markers, enabling CRISPR/Cas9-mediated genome editing. In addition, due to a weak small RNA-based silencing response, S. polyrhiza sustains prolonged periods of transgene activity in transient expression assays. To promote duckweed research and encourage the adoption of S. polyrhiza, we have made SP162 (ID#: 5676) and its genome publicly available and provide here detailed procedures for its cultivation and transformation. Furthermore, we created a web server to explore its genome, retrieve gene sequences, and implement orthologous gene search and a gRNA design function for diverse CRISPR/Cas-based applications (https://agxu.uni-mainz.de/SP162/).

The complexity of varied modifications of chromatin composition is integrated in archetypal combinations called chromatin states that predict the local potential for transcription. The degree of conservation of chromatin states has not been established amongst plants, and how they interact with transcription factors is unknown. Here we identify and characterize chromatin states in the flowering plant Arabidopsis thaliana and the bryophyte Marchantia polymorpha, showing a large degree of functional conservation over more than 450 million years of land plant evolution. We used this new resource of conserved plant chromatin states to understand the influence of chromatin states on gene regulation. We established the preferential association of chromatin states with binding sites and activity of transcription factors. These associations define three main groups of transcription factors that bind upstream of the transcription start site, at the + 1 nucleosome or further downstream of the transcription start site and broadly associate with distinct biological functions including a list of potential candidate pioneer factors we know little about in plants, compared to their important roles in animal stem cells and early development.

Land plant evolution has been marked by bursts of novelty, often underpinned by extensive genomic innovation. A key mechanism driving these changes is horizontal gene transfer (HGT), the process by which genes move between species and even across taxonomic kingdoms. HGT can accelerate evolutionary change through the rapid introduction of new genes yet its importance in plant biology is only beginning to be understood. Here, we review the functional contributions of HGT during the origin and diversification of land plants. We discuss the occurrence of HGT throughout plant evolution and its impact on the origin of defining traits from cell walls to developmental programs. Beyond ancient contributions, HGT continues to drive the emergence of lineage-specific innovations. Recently acquired bacterial and fungal genes make complex functional contributions to processes including stress response, pathogen defence, and development across plant phylogeny. These observations suggest that HGT was, and continues to be, a major force shaping plant evolution, exemplifying the potential significance of HGT in eukaryotic biology more broadly.

Meristems are generative centers from which the bodies of land plants develop. Marchantia polymorpha spores are single cells formed during meiosis that-upon germination-form an early cell mass on which a flat prothallus develops. A single stem cell niche (meristem) develops de novo at the margin of the prothallus to generate the thallus plant body. Here we show that the prothallus forms at the apical pole of the early cell mass and represses the formation of other prothalli in a light- and auxin-dependent manner. Light is required for prothallus development and the formation of an apical auxin signaling minimum, marked by expression of LOW AUXIN RESPONSIVE (MpLAXR). MpLAXR marks the sporeling apex-the site of organogenesis-and predicts the positioning of the future apical stem cell on the prothallus. Ablation of the dominant prothallus induces the development of a suppressed prothallus in a process comparable to apical dominance. A similar molecular program operates during plant regeneration from a single differentiated thallus cell, which regains stemness (pluripotency) upon surgical isolation from its tissue context, and the isolated cell divides, forming an early cell mass with an MpLAXR-expressing apex where a flat prothallus with a single meristem forms. We conclude that a light-dependent, apical auxin signaling minimum is required for the formation of the prothallus and the de novo development of the first meristem in Marchantia polymorpha.

Selective autophagy is a fundamental protein quality control pathway that safeguards proteostasis by degrading damaged or surplus cellular components, particularly under stress. This process is orchestrated by selective autophagy receptors (SARs) that recruit specific cargo for degradation. Although significant strides have been made in understanding the molecular framework of selective autophagy, the diversity of SAR repertoires across species remains largely unexplored. Through a comparative interactome analysis across five model organisms, we identified a suite of conserved and lineage-specific SAR candidates. Among these, we validated coupling of ubiquitin to endoplasmic reticulum degradation- and protein rich in the amino acids E, L, K, and S-domain-containing SAR (CESAR) as a conserved SAR critical for proteostasis under heat stress. CESAR specifically facilitates the degradation of ubiquitinated protein aggregates and is indispensable for heat stress tolerance. Altogether, our study establishes a robust pipeline and a rich resource for SAR discovery. It also positions CESAR as a pivotal regulator of proteostasis, with broad implications for improving stress resilience in plants.

The liverwort Marchantia polymorpha is a key model organism for understanding land plant evolution, development, and gene regulation. To support the growing demand for high-quality genomic resources, we present MarpolBase, a comprehensive and integrated genome database that hosts newly assembled, high-accuracy reference genomes for both the male Tak-1 and female Tak-2 accessions, designated as ver. 7.1 reference genomes. These new assemblies, generated using PacBio HiFi long-read sequencing, represent nearly telomere-to-telomere chromosome-level genomes, with improvements in assembly continuity, annotation accuracy, and structural resolution-especially for repeat-rich regions and sex chromosomes. MarpolBase offers not only access to genome sequences and gene annotations but also provides a unified platform for data exploration, comparative analysis, and community-driven gene nomenclature for Marchantia polymorpha. It includes keyword-searchable gene pages with structural and functional annotations, expression data integration, genome browser visualization, and online analytical and utility tools. By unifying genome assembly, annotation, nomenclature, and analysis tools in a single platform, MarpolBase serves as a central resource for functional genomics and evolutionary studies in M. polymorpha, and a model for future plant genome databases. The genomic resources of MarpolBase are freely available at https://marchantia.info.

Endophytic colonization of Arabidopsis (Arabidopsis thaliana) by the beneficial root endophyte Serendipita indica is characterized by an initial biotrophic phase, followed by a confined host cell death phase that facilitates fungal accommodation. However, the host molecular pathways that restrict S. indica proliferation and regulate symbiosis-associated cell death remain largely unknown. Our study demonstrates that autophagy, a key cellular degradation pathway that maintains homeostasis, is locally activated during colonization and is required to limit fungal proliferation and immunometabolic stress. Autophagy-deficient mutants exhibit elevated basal root cell death, increased colonization, and hypersensitivity to the fungal-derived purine metabolite 2'-deoxyadenosine (dAdo), an immunometabolic signal that modulates host cell viability and reprograms immune and metabolic responses via ENT3 (equilibrative nucleoside transporter 3)-mediated uptake. In ent3 and atg5 ent3 mutants, suppression of dAdo import reduces S. indica-induced cell death, confirming the central role of ENT3-mediated uptake. Despite increased colonization and stress sensitivity, autophagy-deficient plants retain S. indica-mediated root growth promotion, indicating that mutualistic benefits can occur independently of immunometabolic stress resilience. Based on these findings, we propose that autophagy-mediated pro-survival responses are essential for maintaining symbiotic homeostasis by integrating immunometabolic signals and preserving host cell viability.

Protein neofunctionalization is a key driver of cellular complexity. However, subunits of multimeric protein complexes are often thought to be evolutionarily constrained, limiting their capacity for functional divergence. This presents a paradox in plants, where the Exo70 subunit of the exocyst-an octameric complex essential for exocytosis-has undergone striking expansion and diversification. Here we show that electrostatic changes in the N-terminal helix of Exo70 facilitated its physical and functional dissociation from the exocyst, relieving constraints imposed by complex integration. Using Marchantia polymorpha and Arabidopsis thaliana, we demonstrate that this 'complex escape' enables Exo70 paralogues to acquire distinct localizations, interactomes and functions independent of canonical exocytosis. Ancestral reconstructions across land plants reveal that this electrostatic shift predates the extensive radiation of the plant Exo70 protein family, with some lineages later reassociating with the complex. Our findings reveal a reversible mechanism that enabled Exo70 to circumvent the evolutionary and biophysical constraints imposed by complex integration and diversify-a mechanism that could represent a generalizable route to protein neofunctionalization and cellular innovation.

DNA methylation is a key epigenetic mark that impacts gene expression and represses transposable elements in eukaryotes. Numerous examples of cis elements targeted by DNA methylation, particularly at CG sites (mCG), have been reported to be under selective pressure in animals and plants. By contrast, there is limited knowledge of trans regulators of mCG leading to adaptation. Here, a genome-wide association study identifies CELL DIVISION CYCLE-ASSOCIATED PROTEIN 7 (CDCA7) as a major trans determinant of mCG in natural populations of Arabidopsis thaliana. CDCA7 or its paralogue physically binds the chromatin remodeller DECREASE IN DNA METHYLATION 1 (DDM1), which facilitates access of methyltransferases to DNA. Epigenomic analysis shows that while CDCA7 proteins control all DDM1-dependent processes, their predominant function is the maintenance of mCG. We identify a 26-bp promoter indel modulating CDCA7 expression in natural populations and determining the degree of mCG and transposable element silencing. The geographic distribution of CDCA7 alleles suggests that new alleles have repeatedly expanded to novel ecological niches, indicating a potential role in local adaptation. Our findings establish CDCA7 proteins as dedicated regulators of mCG and imply that DDM1 relies on alternative partners to regulate other chromatin features. Broadly, they illustrate how changes in global DNA methylation levels through transcriptional regulation of the epigenetic machinery have the capacity to facilitate local adaptation.

BACKGROUND AND AIMS: Isoetales is a clade of lycopsids that evolved colossal arborescent forms during their Palaeozoic prime but today are represented solely by the small herbaceous monogeneric Isoetes. Despite the differences in scale of taxa in the clade, the rooting system of all members consists of two-parts; rootlets develop from a rhizomorph in a regular pattern termed rhizotaxy. Rhizomorphs are highly diverse in morphology, leading to different terms being used to describe aspects of rhizotaxy in contrasting lineages. Here we set out to investigate the degree to which rhizotaxy was conserved among taxa, aiming to provide a standard geometric definition and developmental interpretation of rhizotaxy. METHODS: We developed a pipeline to quantitatively describe rhizotaxy. This pipeline allowed rootlet arrangement to be captured in 3D, before being visualised on a 2D lattice to which Delaunay triangulation could be applied. This approach offers a standard quantitative method of comparing rhizotaxy across disparate rhizomorphs. Next, to investigate the evolution and development of rhizotaxy we applied our pipeline to 3D reconstructions we generated of the rooting system of the extinct Carboniferous lycopsid, Oxroadia. Finally, we made direct observations of rootlet development in Isoetes using time-course imaging. KEY RESULTS: We demonstrate that rhizotaxy can be described as an equilateral triangular lattice for all members of the Isoetales, including Oxroadia. By combining evidence from direct observation of rootlet development in Isoetes with inferences of rootlet development and the early stages of sporophyte ontogeny of Oxroadia we conclude that the conserved rhizotaxy developed through the process of rootlet intercalation. CONCLUSIONS: We provide a single geometric definition and predicted developmental mechanism for rhizotaxy that applies to all Isoetales. Our findings call into question the literal interpretation that the rhizomorph is a modified shoot.

The prevalence of synthetic associations in GWAS, where non-causal variants become significant by tagging multiple undetected causal variants and not necessarily in strong linkage disequilibrium with any single one, remains unexplored. We introduce a novel machine-learning approach using only genotype data to infer such associations in human GWAS. Our analysis reveals that 3-5% of GWAS Catalog peaks may represent potential synthetic associations, often arising from epistatic interactions between common variants rather than multiple rare variants acting independently. Our findings highlight the need for multi-locus models and emphasize careful GWAS interpretation and follow-up analyses like fine-mapping and trait prediction.

DNA methylation loss at transposable elements (TEs) can affect neighboring genes and be epigenetically inherited in plants, yet the determinants and significance of this additional system of inheritance are unknown. Here, we demonstrate in Arabidopsis thaliana that transgenerational stability of experimentally-induced hypomethylation at TE loci is constrained by small RNAs derived from related copies. Using data from >700 strains collected worldwide, we uncover similar and recurrent hypomethylation at hundreds of these TE loci, often near genes. Most natural epivariants we tested can be inherited without DNA sequence changes and are therefore bona fide epialleles, although genetic factors modulate their recurrence or persistence. Epiallelic variants often cause gene expression changes and may be targets of selection, thus revealing their contribution to heritable phenotypic variation in nature.

The germination of most seeds is influenced by the duration, intensity, and quality of light. The seeds of the model plant Arabidopsis are positive photoblastic and require light to germinate. The germination of negative photoblastic seeds is inhibited by white light. The molecular mechanisms that regulate negative photoblastic germination are unknown due to the lack of a suitable model plant. We identified an accession with negative photoblastic germination in Aethionema arabicum that grows in semi-arid natural habitats. In a forward genetic screen, we identified a mutant - koyash2 (koy2) - that is defective in negative photoblastic germination. There is a nonsense mutation in the gene encoding the phytochrome A photoreceptor in the koy2 mutant. Here, we show that phytochrome A is required for negative photoblastic germination. The defective negative photoblastic phenotype of the koy2 mutant is the result of defective inhibition of germination by the phytochrome A-mediated high irradiance response. This is the first example of phytochrome A-mediated response controlling negative photoblastic seed germination in white, red, far-red, and blue light. We speculate that genetically encoded variation in phytochrome A-mediated germination responses is responsible for local adaptation of Ae. arabicum throughout the Irano-Turanian region.

Making sense of whole-genome polymorphism data is challenging, but it is essential for overcoming the biases in SNP data. Here we analyze 27 genomes of Arabidopsis thaliana to illustrate these issues. Genome size variation is mostly due to tandem repeat regions that are difficult to assemble. However, while the rest of the genome varies little in length, it is full of structural variants, mostly due to transposon insertions. Because of this, the pangenome coordinate system grows rapidly with sample size and ultimately becomes 70% larger than the size of any single genome, even for n = 27. Finally, we show how short-read data are biased by read mapping. SNP calling is biased by the choice of reference genome, and both transcriptome and methylome profiling results are affected by mapping reads to a reference genome rather than to the genome of the assayed individual.

Galls, especially those induced by insects, represent one of the most dramatic examples of plant developmental reprogramming, combining complex de novo organogenesis with compromised defence. Insect-induced galls are not just a fascinating natural phenomenon but a unique system for future discoveries in developmental biology, plant defence, and evolutionary ecology. Gall development is under the control of their insect manipulators and in sync with insect growth to provide tailored nutritive and protective environments. But this alone does not explain the huge diversity in their morphology which evolved within complex ecological niches. In this review, we summarise recent findings in this underexplored field and examine the defining features of insect-induced galls compared to non-gall herbivores, microbial gall inducers, and symbionts. By exploring commonalities and differences in developmental reprogramming, defence and nutrition, we highlight the uniqueness of insect-induced galls and their potential for discoveries in plant biology.

The plant vacuole, traditionally viewed as a static storage organelle, has recently emerged as a dynamic hub orchestrating signaling, metabolic integration, and stress responses. This review synthesizes recent advances that position the vacuole as a pivotal regulator of plant development and environmental adaptation. We discuss the vacuole's multifaceted roles in ion sequestration, lipid trafficking, mechanosensing, and signal transduction, highlighting its central role in preserving cellular homeostasis. We summarize recent data supporting two distinct forms of vacuolar biogenesis: inheritance from existing organelles and de novo formation. Lastly, we discuss our recent findings that define a vacuolar quality control (VQC) pathway, safeguarding tonoplast integrity during stress. Collectively, these insights redefine our understanding of the vacuole's essential contributions to plant physiology and resilience, advocating for an updated conceptual framework that recognizes the vacuole as a central hub for developmental processes and environmental adaptation.

Autophagy sustains cellular health by recycling damaged or excess components through autophagosomes. It is mediated by conserved ATG proteins, among which the ubiquitin-like ATG8 protein plays a central role by linking cargo to the growing autophagosomes. Unlike most ATG proteins, the ATG8 gene family is significantly expanded in vascular plants, but its functional specialization remains poorly understood. Using transcriptional and translational reporters in Arabidopsis thaliana, we revealed that ATG8 isoforms are differentially expressed across tissues and form distinct autophagosomes. To explore ATG8 specialization, we generated the nonuple Δatg8 mutant lacking all nine ATG8 isoforms. The mutant displayed hypersensitivity to carbon and nitrogen starvation, coupled with defects in bulk and selective autophagy as shown by biochemical and ultrastructural analyses. Complementation experiments demonstrated that ATG8A could rescue both carbon and nitrogen starvation phenotypes, whereas ATG8H could only complement carbon starvation. Proximity labeling proteomics further identified isoform-specific interactors under nitrogen starvation, underscoring their functional divergence. These findings provide genetic evidence for functional specialization of ATG8 isoforms in plants and lay the foundation for investigating their roles in diverse cell types and stress conditions.

Plant galls are abnormal growing tissues induced by various parasitic organisms, exhibiting diverse and complex morphologies. Typically, these galls differ significantly in appearance from their host plants. Here, we report that larvae of a parasitic fly generate unique, rosette galls on Aster scaber, a perennial herb. These galls develop from vegetative organs after the larvae reprogram floral gene expression. To investigate the underlying mechanisms, we conducted whole-genome sequencing and transcriptome analysis. Our findings reveal that the larvae induce host organ dedifferentiation into an amorphous callus, activate floral genes, and selectively suppress genes associated with carpel development. As a result, the pseudoflowers consist solely of tepal-like leaflets and a specialized chamber, and the larvae influence pigment biosynthesis. Hijacking plants developmental gene networks by insects to sequentially mediate dedifferentiation, cytokinin regulation, and tepal-like leaflets formation provides a framework to study highly elaborate forms of parasitism and symbiosis between plants and insects.

Evaluating population genetic inference methods is challenging due to the complexity of evolutionary histories, potential model misspecification, and unconscious biases in self-assessment. The Genomic History Inference Strategies Tournament (GHIST) is a community-driven competition designed to evaluate methods for inferring evolutionary history from population genomic data. The inaugural GHIST competition ran from July to November 2024 and featured four demographic history inference challenges of varying complexity: a bottleneck model, a split with isolation model, a secondary contact model with demographic complexity, and an archaic admixture model. Data were provided as error-free VCF files, and participants submitted numerical parameter estimates that were scored by relative root mean squared error. Approximately 60 participants competed, using diverse approaches. Results revealed the current dominance of methods based on site frequency spectra, while highlighting the advantages of flexible model-building approaches for complex demographic histories. We discuss insights regarding the competition and outline the next iteration, which is ongoing with expanded challenge diversity. By providing standardized benchmarks and highlighting areas for improvement, GHIST represents a substantial step toward more reliable inference of evolutionary history from genomic data.

The ability of a cell to polarise, and direct cell growth or orient cell division, for example, is fundamental for the morphogenesis of multicellular organisms. A key molecular system for signalling cell polarity in diverse eukaryotes involves the RHO family of small GTPases. Since its origin in early eukaryotes, the RHO family has evolved independently in different lineages, and the plant-specific subfamily of RHO - RHO of plants (ROP) - was established in the streptophyte algal ancestors of land plants. Insights from both bryophytes and vascular plants have revealed conserved roles for ROP signalling in land plant morphogenesis. Synthesising our current understanding of how ROP signalling evolved and how it regulates morphogenesis in extant land plants, we propose that this polarity signalling system was co-opted to spatially coordinate morphogenetic mechanisms that evolved in the algal ancestors of land plants.

Plant shoot branches are formed by the initiation and activation of generative centers known as meristems. In dichotomously branching plants, such as many bryophytes and lycophytes, new meristems are formed when a pre-existing meristem splits into two daughter meristems. These meristems may be active and produce shoot branches or may be inactive. Here, we show that in conditions where meristem inactivation occurs, such as simulated shade, the position of the inactive meristem alternated between either side of the plant body in the liverwort Marchantia polymorpha. Using this predictable pattern, we generated transcriptomes of active and inactive meristems and identified the cytochrome P450 monooxygenase, MpCYP78E1, as a novel regulator of meristem activity. MpCYP78E1 reporter expression was higher in active meristems than in inactive meristems. More meristems were active in loss-of-function mutants than in wild type, while fewer meristems were active in gain-of-function mutants, indicating that MpCYP78E1 inhibits meristem activity. Furthermore, unlike wild type, Mpcyp78e1 loss-of-function mutants produced supernumerary meristems from the center of the mature plant body. We conclude that MpCYP78E1 inhibits both meristem initiation and activity to modulate shoot branching architecture.

Inferring genealogical relationships of wild populations is useful because it gives direct estimates of mating patterns and variance in reproductive success. Inference can be improved by including information about parentage shared between siblings, or by modelling phenotypes or population data related to mating. However, we currently lack a framework to infer parent-offspring relationships, sibships and population parameters in a single analysis. To address this, we here extend a previous method, Fractional Analysis of Paternity and Sibships, to include population data for the case where one parent is known. We illustrate this with the example of pollen dispersal in a natural hybrid zone population of the snapdragon Antirrhinum majus. Pollen dispersal is leptokurtic, with half of mating events occurring within 30 m, but with a long tail of mating events up to 859 m. Using simulations, we find that both sibship and population information substantially improve pedigree reconstruction, and that we can expect to resolve median dispersal distances with high accuracy.

Nucleosomes are thought to be structural barriers to transcription, establishing a restrictive ground state that must be destabilized for gene expression. However, structural insights have revealed that transcription can proceed in the presence of nucleosomes, suggesting that this model is incomplete. Here, we reconstituted H2A.Z sequences resulting from more than a billion years of eukaryotic evolution in a single synthetic host system, interrogating their impact on transcription. We identified single-residue substitutions within the ultra-conserved core domain loop 2 (L2) of H2A.Z as sufficient to confer emergent properties and drive neofunctionalization. Such L2 neomorphs acquired a direct interaction with transcription elongation factor Spt6, rewiring gene expression by tuning polymerase processivity. We conclude that even minimal changes in histone sequences can transform their function, underscoring the evolutionary potential of the histone core domain to drive regulatory innovation and highlighting a previously unappreciated role of the histone core domain in transcriptional regulation.

Extrachromosomal circular DNA (eccDNA) has been described in several eukaryotic species and has been shown to impact phenomena as diverse as cancer and herbicide tolerance. EccDNA is thought to arise mainly through transposable element (TE) mobilization. Because studies based on short-read sequencing cannot efficiently identify full-length eccDNA forms generated from TEs, we employed the CIDER-Seq pipeline based on long-read sequencing, to obtain full-length eccDNAs from Arabidopsis. The generated eccDNA datasets identified centromeric/pericentromeric regions as hotspots of eccDNAs with several eccDNA molecules originating from Helitron and LTR TEs. To investigate the role of epigenetic marks on TE-derived eccDNA biogenesis, we studied Arabidopsis methylation mutants dcl3, rdr6, ros1, and ddm1. Contrasting the TE-suppression previously reported in the hypermethylated ros1 mutants, we identified activation of TEs in ros1, specifically of LTR/Gypsy TEs. An enrichment of LTR/Copia elements was identified in actively dividing calli and the shoot apical meristem (SAM). Uncharacterized "variable TEs" with high eccDNA and expression were identified in the SAM, including ATCOPIA58. Together, our study reveals the genomic origins of eccDNAs and delineates the link between epigenetic regulation, transposon mobilization, and eccDNA biogenesis.

Wheat genetic resources hold the diversity required to mitigate agricultural challenges from climate change and reduced inputs. Using DArTseq, we genotype 461 wheat landraces and cultivars and evaluate them for powdery mildew resistance. By developing a k-mer-based GWAS approach with fully assembled genomes of Triticum aestivum and its progenitors, we uncover 25% more resistance-associated k-mers than single-reference methods, outperforming SNP-based GWAS in both loci detection and mapping precision. In total, we detect 34 powdery mildew resistance loci, including 27 potentially novel regions. Our approach underscores the importance of integrating multiple reference genomes to unlock the potential of wheat germplasm.

In this study, we characterized a panel of 1264 maize near-isogenic lines (NILs), developed from crosses between 18 diverse inbred lines and the recurrent parent B73, referred to as nested NILs (nNILs). In this study, 888 of the nNILs were genotyped using genotyping-by-sequencing (GBS). Subsequently, 24 of these nNILs, and all the parental lines, were re-genotyped using a high-density single nucleotide polymorphism (SNP) chip. A novel pipeline for calling introgressions, which does not rely on knowing the donor parent of each nNIL, was developed based on a hidden Markov model (HMM) algorithm. By comparing the introgressions detected using GBS data with those identified using chip data, we optimized the HMM parameters for analyzing the entire nNIL population. A total of 2969 introgressions were identified across the 888 nNILs. Individual introgression blocks ranged from 21 bp to 204 Mbp, with an average size of 17 Mbp. By comparing SNP genotypes within introgressed segments to the known genotypes of the donor lines, we determined that in about one third of the lines, the identity of the donors did not match expectation based on their pedigrees. We characterized the entire nNIL population for three foliar diseases. Using these data, we mapped a number of quantitative trait loci (QTL) for disease resistance in the nNIL population and observed extensive variation in effects among the alleles from different donor parents at most QTL identified. This population will be of significant utility for dissecting complex agronomic traits and allelic series in maize.

Roots are essential for plant adaptation to changing environments, yet the role of roots in crop domestication remains unclear. This study examined the evolution of root phenotypes from teosinte to maize, a transition resulting in reduced nodal root number (NRN), multiseriate cortical sclerenchyma (MCS), and increased seminal root number (SRN). We reconstructed the root phenotypes of maize and teosinte, as well as the environments of the Tehuacán Valley over the last 18 000 yr using a combination of ancient DNA, paleobotany, and functional-structural modeling. Our models reveal that increasing Holocene atmospheric CO2 concentrations favored the appearance of reduced NRN and MCS between 12 000 and 8000 yr before present (yrBP), promoting deeper root systems. The advent of irrigation by 6000 yrBP switched nitrogen distribution from topsoil to subsoil domains, a change that increased the utility of reduced NRN and MCS. Comparison of allelic frequencies among ancient samples ranging from 5500 to 500 yrBP suggests that increased SRN may have appeared c. 3500 yrBP, coinciding with a period of increased human population, agricultural intensification, and soil degradation. Our results suggest that root phenotypes that enhance plant performance under nitrogen stress are important for maize adaptation to changing agricultural practices in the Tehuacán Valley.

N4-methylcytosine (4mC) is an important DNA modification in prokaryotes, but its relevance and even its presence in eukaryotes have been mysterious. Here we show that spermatogenesis in the liverwort Marchantia polymorpha involves two waves of extensive DNA methylation reprogramming. First, 5-methylcytosine (5mC) expands from transposons to the entire genome. Notably, the second wave installs 4mC throughout genic regions, covering over 50% of CG sites in sperm. 4mC requires a methyltransferase (MpDN4MT1a) that is specifically expressed during late spermiogenesis. Deletion of MpDN4MT1a alters the sperm transcriptome, causes sperm swimming and fertility defects, and impairs post-fertilization development. Our results reveal extensive 4mC in a eukaryote, identify a family of eukaryotic methyltransferases, and elucidate the biological functions of 4mC in reproductive development, thereby expanding the repertoire of functional eukaryotic DNA modifications.

Local adaptation has been demonstrated in many organisms, but the traits involved, and the temporal and spatial scales at which selection acts are generally unknown. We carried out a multi-year study of 200 accessions (natural inbred lines) of Swedish Arabidopsis thaliana using local field sites and a combination of common-garden experiments that measured adult survival and fecundity, and selection experiments that measured fitness over the full life cycle. We found evidence of strong and variable selection, with particular genotypes favored more than five-fold in certain years and locations. Fecundity showed evidence of classical local adaptation, with accessions generally performing better close to their home. However, southern accessions usually had the highest fecundity-but were far more sensitive to harsh winters and slug herbivory, which strongly decreased both survival and fecundity. Accessions originally sampled on beaches had low fecundity in all environments, but massively outperformed all other accessions in the selection experiments, presumably due to an advantage during seedling establishment associated with their very large seeds. We conclude that local adaptation in A. thaliana reflects strong temporally and spatially varying selection on multiple traits, generally involving trade-offs and different life-history strategies, making fitness difficult to predict and measure.

'Tagmentation' approaches to bisulphite sequencing use a transposase to simultaneously make double-stranded breaks and ligate adaptors to the resulting fragments, allowing for higher throughput with less starting material. However, it has also been noted that certain tagmentation protocols have an unusually high number unmethylated cytosines that are not converted to thymine. Here we describe this phenomenon in detail, and find that results are consistent with single strand nicks by the transposase, followed by strand displacement of part or all of the DNA fragment, leading to erroneous incorporation of methylated cytosines. Nevertheless we show that these errors can be accounted for in downstream analysis and need not impede biological conclusions. We provide a Python package to allow users to implement this framework. Ultimately the additional effort of accounting for errors must be traded off against the scalability of the protocol in planning experiments.

Autophagy has emerged as an essential quality control pathway in plants that selectively and rapidly removes damaged or unwanted cellular components to maintain cellular homeostasis. It can recycle a broad range of cargoes, including entire organelles, protein aggregates, and even invading microbes. It involves the de novo biogenesis of a new cellular compartment, making it intimately linked to endomembrane trafficking pathways. Autophagy is induced by a wide range of biotic and abiotic stress factors, and autophagy mutant plants are highly sensitive to stress, making it an attractive target for improving plant stress resilience. Here, we critically discuss recent discoveries related to plant autophagy and highlight open questions and future research areas.

To identify novel genes engaged in plant epidermal development, we characterized the phenotypic variability of rosette leaf epidermis of 310 sequenced Arabidopsis thaliana accessions, focusing on trichome shape and distribution, compositional characteristics of the trichome cell wall, and histologically detectable metal ion distribution. Some of these traits correlated with cLimate parameters of our accession's locations of origin, suggesting environmental selection. A novel metal deposition pattern in stomatal guard cells was observed in some accessions. Subsequent GWAS analysis identified 1546 loci with protein sequence-altering SNPs associated with one or more traits, including 5 genes with previously reported relevant mutant phenotypes and 80 additional genes with known or predicted roles in relevant developmental and cellular processes. Some candidates, including GFS9/TT9, exhibited environmentally correlated allele distribution. Several large gene famiLies, namely DUF674, DUF784, DUF1262, DUF1985, DUF3741, cytochrome P450, receptor-Like kinases, Cys/His-rich C1 domain proteins and formins were overrepresented among the candidates for various traits, suggesting epidermal development-related functions. A possible participation of formins in guard cell metal deposition was supported by observations in available loss of function mutants. Screening of candidate gene lists against the STRING interactome database uncovered several predominantly nuclear protein interaction networks with possible novel roles in epidermal development.

Paralogous variants of canonical histones guide accessibility to DNA and function as additional layers of genome regulation. Across eukaryotes, the mechanism of action and functional significance of several variants of core histones are well known except those of histone H4. Here we show that a variant of H4 (H4.V) expressing tissue-specifically among Oryza members mediated specific epigenetic changes contributing to salt tolerance. H4.V was incorporated into specific heterochromatic sites, where it blocked the deposition of active histone marks. Stress-dependent redistribution of H4.V enabled the incorporation of acetylated H4 lysine 5 (H4K5ac) in the gene bodies. The misexpression of H4.V led to defects in reproductive development and in mounting salt stress responses. H4.V formed homotypic nucleosomes and mediated these alterations by conferring distinct molecular properties to the nucleosomes, as seen with cryo electron microscopy structures and biochemical assays. These results reveal not only an H4 variant among plants but also a chromatin regulation that might have contributed to the adaptation of semi-aquatic Oryza members.

The discovery of histone H3 copper reductase activity provides a novel metabolic framework for understanding the functions of core histone residues, which, unlike N-terminal residues, have remained largely unexplored. We previously demonstrated that histone H3 cysteine 110 (H3C110) contributes to cupric (Cu2+) ion binding and its reduction to the cuprous (Cu1+) form. However, this residue is absent in Saccharomyces cerevisiae, raising questions about its evolutionary and functional significance. Here, we report that H3C110 has been lost in many fungal lineages despite near-universal conservation across eukaryotes. Introduction of H3C110 into S. cerevisiae increased intracellular Cu1+ levels and ameliorated the iron homeostasis defects caused by inactivation of the Cup1 metallothionein or glutathione depletion. Enhanced histone copper reductase activity also extended replicative life span under oxidative growth conditions but reduced it under fermentative conditions. Our findings suggest that a trade-off between histone copper reductase activity, iron metabolism, and life span may underlie the loss or retention of H3C110 across eukaryotes.

Meiotic crossovers rearrange allele combinations and create offspring diversity. Crossovers occur nonrandomly along chromosomes, predominantly in distal euchromatin and less in pericentromeric heterochromatin marked with histone H3 lysine 9 dimethylation (H3K9me2) and the H2A variant H2A.W in Arabidopsis thaliana. Loss of H3K9me2 increases heterochromatic crossovers, but how H2A.W affects crossover formation in pericentromeric regions is unknown. Here, we report that H2A.W is required to restrict heterochromatic crossovers in Arabidopsis. Using meiosis-specific microRNA-induced gene silencing (meiMIGS) and fluorescence-tagged recombination reporters, we show that meiotic knockdown of H2A.W.6, H2A.W.7, and H2A.W.12 (meiMIGS-H2A.W.6/7/12) increases pericentromeric crossovers. High-resolution genomic maps of crossovers show that meiMIGS-H2A.W.6/7/12 enhances heterochromatic crossovers, similar to meiMIGS plants silencing the H3K9me2 pathway. Consistently, genome-wide crossover maps show that the mutants h2a.w.6, h2a.w.7, h2a.w.6 h2a.w.7, and h2a.w.6 h2a.w.7 h2a.w.12, but not h2a.w.12, exhibit a similar increase in heterochromatic crossovers to meiMIGS-H2A.W.6/7/12, demonstrating that H2A.W.6 and H2A.W.7 limit heterochromatic crossovers. Profiling of genome-wide nucleosome density using micrococcal nuclease sequencing reveals that h2a.w mutants with increased heterochromatic crossovers have increased heterochromatin accessibility, with lower H3K9me2 levels during meiosis. Our findings shed light on the role of H2A.W variants as heterochromatin compaction factors that suppress meiotic crossovers within the pericentromeric regions.

Epigenetic silencing is essential for regulating gene expression and cellular diversity in eukaryotes. While DNA and H3K9 methylation silence transposable elements (TEs), H3K27me3 marks deposited by the Polycomb repressive complex 2 (PRC2) silence varying proportions of TEs and genes across different lineages. Despite the major development role epigenetic silencing plays in multicellular eukaryotes, little is known about how epigenetic regulatory networks were shaped over evolutionary time. Here, we analyze epigenomes from diverse species across the green lineage to infer the chronological epigenetic recruitment of genes during land plant evolution. We first reveal the nature of plant heterochromatin in the unicellular chlorophyte microalga Chlorella sorokiniana and identify several genes marked with H3K27me3, highlighting the deep origin of PRC2-regulated genes in the green lineage. By incorporating genomic phylostratigraphy, we show how genes of differing evolutionary age occupy distinct epigenetic states in plants. While young genes tend to be silenced by H3K9 methylation, genes that emerged in land plants are preferentially marked with H3K27me3, some of which form part of a common network of PRC2-repressed genes across distantly related species. Finally, we analyze the potential recruitment of PRC2 to plant H3K27me3 domains and identify conserved DNA-binding sites of ancient transcription factor families known to interact with PRC2. Our findings shed light on the conservation and potential origin of epigenetic regulatory networks in the green lineage, while also providing insight into the evolutionary dynamics and molecular triggers that underlie the adaptation and elaboration of epigenetic regulation, laying the groundwork for its future consideration in other eukaryotic lineages.

Plants have an impressive repertoire to react to stress conditions that limit regular growth. Elevated temperatures beyond the optimal range cause rapid and specific transcriptional responses, resulting in developmental alterations and plasticity. Heat stress also causes chromatin decondensation and activation of some transposable elements (TEs), endangering genomic integrity. This is especially risky for stem cells in the shoot apical meristem (SAM) that potentially contribute to the next generation. We examined how the heat stress response in SAM stem cells of Arabidopsis (Arabidopsis thaliana) is different from that in other tissues and whether the elements of epigenetic TE control active in the meristem are involved in specific heat protection of stem cells. Using fluorescence-activated nuclear sorting to isolate and characterize SAM stem cells after exposure to conditions that activate a heat-responsive TE, we found that SAM stem cells maintain their developmental program and suppress the heat-response pathways dominating in surrounding somatic cells. Furthermore, mutants defective in DNA methylation recovered less efficiently from heat stress and persistently activated heat response factors and heat-responsive TEs. Heat stress also induced epimutations at the level of DNA methylation, especially in the CHG sequence context. Regions with modified DNA methylation patterns remained altered for at least 3 wk beyond the stress. These findings urge for disentangling cell type-specific responses under different stress types, especially for stem cells as bridges to the next generation.

A handful of model plants have provided insight into silencing of transposable elements (TEs) through RNA-directed DNA methylation (RdDM). Guided by 24 nt long small-interfering RNAs (siRNAs), this epigenetic regulation installs DNA methylation and histone modifications like H3K9me2, which can be subsequently maintained independently of siRNAs. However, the genome of the clonally propagating duckweed Spirodela polyrhiza (Lemnaceae) has low levels of DNA methylation, very low expression of RdDM components, and near absence of 24 nt siRNAs. Moreover, some genes encoding RdDM factors, DNA methylation maintenance, and RNA silencing mechanisms are missing from the genome. Here, we investigated the distribution of TEs and their epigenetic marks in the Spirodela genome. Although abundant degenerated TEs have largely lost DNA methylation and H3K9me2 is low, they remain marked by the heterochromatin-associated H3K9me1 and H3K27me1 modifications. In contrast, we find high levels of DNA methylation and H3K9me2 in the relatively few intact TEs, which are source of 24 nt siRNAs, like RdDM-controlled TEs in other angiosperms. The data suggest that, potentially as adaptation to vegetative propagation, RdDM extent, silencing components, and targets are different from other angiosperms, preferentially focused on potentially intact TEs. It also provides evidence for heterochromatin maintenance independently of DNA methylation in flowering plants. These discoveries highlight the diversity of silencing mechanisms that exist in plants and the importance of using disparate model species to discover these mechanisms.

Bifacial cambium, which produces xylem and phloem, and monopodial architecture, characterized by apical dominance and lateral branching from axillary buds, are key developmental features of seed plants, consisting of angiosperms and gymnosperms. These allow seed plants to adapt to diverse environments by optimizing resource allocation and structural integrity. In seed plants, SUPPRESSOR OF MAX2-LIKE (SMXL) family members function in phloem development and strigolactone-induced inhibition of axillary bud outgrowth. Although strigolactone signaling regulates most SMXL family members, the only known regulator of SMXL4 and SMXL5 is the RNA-binding protein JULGI. We demonstrate that in angiosperms, by directly regulating SMXL4/5 expression, JULGI uncouples SMXL4/5 activity from strigolactone signaling. JULGI and ancestral SMXLs from seedless vascular plants or SMXL4/5 from seed plants are coexpressed in the phloem tissues of vascular plants, from lycophytes to angiosperms. Core angiosperm SMXL4/5 mRNAs contain a G-rich element in the 5' untranslated region (UTR) that serves as a target sequence for JULGI to negatively regulate SMXL4/5 expression. Heterologous expression of JULGIs from various angiosperms rescued the Arabidopsis jul1 jul2 mutant. Expressing SMXL4/5s from seed plants and ancestral SMXLs rescued Arabidopsis smxl4 smxl5. Angiosperm SMXL4/5s lack an RGKT motif for proteasomal degradation. Indeed, treatment with the synthetic strigolactone analog rac-GR24 induced proteasomal degradation of SMXL from ferns and SMXL5a from gymnosperms, but not SMXL4/5 from angiosperms. These findings suggest that in ancestral angiosperms, the 5' UTR of SMXL4/5 gained G-rich elements, creating a regulatory module with JULGI that allows the phloem development pathway to act independently of strigolactone signaling.

Sexual reproduction results in the development of haploid and diploid cell states during the life cycle. In bryophytes, the dominant multicellular haploid phase produces motile sperm that swim through water to the egg to effect fertilization from which a relatively small diploid phase develops. In angiosperms, the reduced multicellular haploid phase produces non-motile sperm that is delivered to the egg through a pollen tube to effect fertilization from which the dominant diploid phase develops. These different life cycle characteristics are likely to impact the distribution of genetic variation among populations. However, little is known about the distribution of genetic variation among wild populations of bryophytes. To investigate how genetic variation is distributed among populations of a bryophyte and to establish the foundation for population genetics research in bryophytes, we described the genetic diversity of collections of Marchantia polymorpha subsp. ruderalis, a cosmopolitan ruderal liverwort. We identified 78 genetically unique (non-clonal) from a total of 209 sequenced accessions collected from 37 sites in Europe and Japan. There was no detectable population structure among European populations but significant genetic differentiation between Japanese and European populations. By associating genetic variation across the genome with global climate data, we showed that temperature and precipitation influence the frequency of potentially adaptive alleles. This collection establishes the core of an experimental platform that exploits natural genetic variation to answer diverse questions in biology.

The transition from free-living to parasitic lifestyles induces major shifts in evolution, and nowhere is this more acute than in apicomplexans-obligate intracellular parasites of animals that evolved from photosynthetic algae.1 In other cases where photosynthesis has been lost, including most apicomplexans, chlorophyll is also absent, but in coral-infecting apicomplexans (corallicolids), chlorophyll biosynthesis genes are retained in the plastid genome despite their lack of photosystems.2 This suggests that the loss of photosynthesis and chlorophyll were decoupled in this lineage, but because these observations are only based on plastid genomes, two fundamental questions remain unclear. First, how this impacted apicomplexan evolution as a whole is unclear because there are conflicting phylogenetic positions for corallicolids: plastid gene phylogenies place them at the base of the apicomplexans, whereas nuclear rRNA places them with late-branching coccidians (suborder Eimeriorina).2,3 Second, it is unclear if chlorophyll or a metabolic intermediate is synthesized, as most chlorophyll biosynthesis enzymes are encoded in the nucleus. To address these questions, we have sequenced transcriptomes from two corallicolids, infecting Parazoanthus swiftii and Madracis mirabilis hosts. Phylogenomic data strongly support a late-branching relationship closer with coccidians, specifically with the protococcidians and the newly discovered ichthyocolids. We also find evidence for the expression of nucleus-encoded enzymes involved in chlorophyll biosynthesis in corallicolids and protococcidians. Overall, we conclude that chlorophyll synthesis was likely retained through the early evolution of the group and then lost approximately 10 times independently, emphasizing the impact of parallel evolutionary changes in parasitic transitions.

Transport protein particle (TRAPP) tethering complexes are known for their function as Rab GTPase exchange factors. Two versions of the complex are considered functionally separate: TRAPPII, an activator of the Rab11 family (RabA in plants) GTPases that function in post-Golgi sorting, and TRAPPIII, activating Rab1 family (RabD in plants) members that regulate endoplasmic reticulum (ER)-to-Golgi trafficking and autophagy. In Arabidopsis (Arabidopsis thaliana), the TRAPPIII complex has been identified and its subunit composition established, but little is known about its functions. Here, we found that binary subunit interactions of the plant TRAPPIII complex are analogous to those of metazoan TRAPPIII, with the 2 large subunits TRAPPC8 and TRAPPC11 linking the TRAPP core and the small C12 to C13 dimer. To gain insight into the functions of TRAPPIII in plants, we characterized 2 A. thaliana trappc8 mutants. These mutants display abnormalities in plant morphology, particularly in flower and seed development. They also exhibit autophagic defects, a constitutive ER stress response, and elevated levels of the ER lipid dolichol (Dol), which is an indispensable cofactor in protein glycosylation. These results indicate that plant TRAPPC8 is involved in multiple cellular trafficking events and suggest a link between ER stress responses and Dol levels.

Vacuoles are essential for cellular metabolism and growth and the maintenance of internal turgor pressure. They sequester lytic enzymes, ions and secondary metabolites that, if leaked into the cytosol, could lead to cell death. Despite their pivotal roles, quality control pathways that safeguard vacuolar integrity have remained elusive in plants. Here we describe a conserved vacuolar quality control pathway that is activated upon cell wall damage in a turgor-pressure-dependent manner. Cell wall perturbations induce a distinct modification-ATG8ylation-on the vacuolar membrane (tonoplast) that is regulated by the V-ATPase and ATG8 conjugation machinery. Genetic disruption of tonoplast ATG8ylation impairs vacuolar integrity, leading to cell death. Together, our findings reveal a homeostatic pathway that preserves vacuolar integrity upon cell wall damage.

Recent research indicates an involvement of microautophagy in the uptake of seed storage proteins (SSPs) into the plant-specific protein storage vacuole (PSV), particularly in cereal grains. However, because microautophagy plays a vital role in cellular homeostasis by degrading and recycling cellular components, we question whether it is a suitable term for a process involved in long-term storage. Additionally, because fission-type microautophagy shares mechanistic similarities with the intraluminal vesicle (ILV) formation of multivesicular bodies (MVBs), we draw parallels between microautophagy and membrane remodeling facilitated by the endosomal sorting complex required for transport (ESCRT). Finally, we propose that the complex structure of cereal endosperm is an optimal tissue to study microautophagy in a plant- and tissue-specific context to decipher its molecular regulation in anabolism and catabolism.

BACKGROUND: The DNA/H3K9 methylation and Polycomb-group proteins (PcG)-H3K27me3 silencing pathways have long been considered mutually exclusive and specific to transposable elements (TEs) and genes, respectively in mammals, plants, and fungi. However, H3K27me3 can be recruited to many TEs in the absence of DNA/H3K9 methylation machinery and sometimes also co-occur with DNA methylation. RESULTS: In this study, we show that TEs can also be solely targeted and silenced by H3K27me3 in wild-type Arabidopsis plants. These H3K27me3-marked TEs not only comprise degenerate relics but also seemingly intact copies that display the epigenetic features of responsive PcG target genes as well as an active H3K27me3 regulation. We also show that H3K27me3 can be deposited on newly inserted transgenic TE sequences in a TE-specific manner indicating that silencing is determined in cis. Finally, a comparison of Arabidopsis natural accessions reveals the existence of a category of TEs-which we refer to as "bifrons"-that are marked by DNA methylation or H3K27me3 depending on the accession. This variation can be linked to intrinsic TE features and to trans-acting factors and reveals a change in epigenetic status across the TE lifespan. CONCLUSIONS: Our study sheds light on an alternative mode of TE silencing associated with H3K27me3 instead of DNA methylation in flowering plants. It also suggests dynamic switching between the two epigenetic marks at the species level, a new paradigm that might extend to other multicellular eukaryotes.

In organisms ranging from vertebrates to plants, major components of centromeres are rapidly evolving repeat sequences, such as tandem repeats (TRs) and transposable elements (TEs), which harbour centromere-specific histone H3 (CENH3)1,2. Complete centromere structures recently determined in human and Arabidopsis suggest frequent integration and purging of retrotransposons within the TR regions of centromeres3-5. Despite the high impact of 'centrophilic' retrotransposons on the paradox of rapid centromere evolution, the mechanisms involved in centromere targeting remain poorly understood in any organism. Here we show that both Ty3 and Ty1 long terminal repeat retrotransposons rapidly turnover within the centromeric TRs of Arabidopsis species. We demonstrate that the Ty1/Copia element Tal1 (Transposon of Arabidopsis lyrata 1) integrates de novo into regions occupied by CENH3 in Arabidopsis thaliana, and that ectopic expansion of the CENH3 region results in spread of Tal1 integration regions. The integration spectra of chimeric TEs reveal the key structural variations responsible for contrasting chromatin-targeting specificities to centromeres versus gene-rich regions, which have recurrently converted during the evolution of these TEs. Our findings show the impact of centromeric chromatin on TE-mediated rapid centromere evolution, with relevance across eukaryotic genomes.

Dynamic shuttling of proteins between the nucleus and cytoplasm orchestrates vital functions in eukaryotes. Here, we reveal the multifaceted functions of Arabidopsis Sin3-associated protein 18 kDa (SAP18) in the regulation of development and heat-stress tolerance. Proteomic analysis demonstrated that SAP18 is a core component of the nuclear apoptosis- and splicing-associated protein (ASAP) complex in Arabidopsis, contributing to the precise splicing of genes associated with leaf development. Genetic analysis further confirmed the critical role of SAP18 in different developmental processes as part of the ASAP complex, including leaf morphogenesis and flowering time. Interestingly, upon heat shock, SAP18 translocates from the nucleus to cytoplasmic stress granules and processing bodies. The heat-sensitive phenotype of a SAP18 loss-of-function mutant revealed a novel role for SAP18 in plant thermoprotection. These findings significantly expand our understanding of the relevance of SAP18 for plant growth, linking nuclear splicing with cytoplasmic stress responses and providing new perspectives for future exploration of plant thermotolerance mechanisms.

Intracellular recycling via autophagy is governed by post-translational modifications of the autophagy-related (ATG) proteins. One notable example is ATG4-dependent delipidation of ATG8, a process that plays critical but distinct roles in autophagosome formation in yeast and mammals. Here, we aim to elucidate the specific contribution of this process to autophagosome formation in species representative of evolutionarily distant green plant lineages: unicellular green alga Chlamydomonas reinhardtii, with a relatively simple set of ATG genes, and a vascular plant Arabidopsis thaliana, harboring expanded ATG gene families. Remarkably, the more complex autophagy machinery of Arabidopsis renders ATG8 delipidation entirely dispensable for the maturation of autophagosomes, autophagic flux, and related stress tolerance; whereas autophagy in Chlamydomonas strictly depends on the ATG4-mediated delipidation of ATG8. Importantly, we also demonstrate the distinct impact of different Arabidopsis ATG8 orthologs on autophagosome formation, especially prevalent under nitrogen depletion, providing new insight into potential drivers behind the expansion of the ATG8 family in higher plants. Our findings underscore the evolutionary diversification of the molecular mechanism governing the maturation of autophagosomes in eukaryotic lineages and highlight how this conserved pathway is tailored to diverse organisms.

Plants adapted to life on land by developing diverse anatomical features across lineages. The molecular basis of these innovations often involves the emergence of new genes or establishing new connections between conserved elements, though evidence for evolutionary genetic circuit rewiring remains scarce. Here, we show that the thermospermine-dependent pathway regulating vascular cell proliferation in Arabidopsis thaliana operates as two distinct modules with different functions in the bryophyte Marchantia polymorpha. One module controls dichotomous branching at meristems, while the other one modulates gemmae and rhizoid production in the thallus. Heterologous assays and comparative expression analyses reveal that the molecular links between these modules, forming a unified circuit in vascular plants, emerged early in tracheophyte evolution. Our results illustrate how the thermospermine-dependent circuit elements followed two divergent evolutionary trajectories in bryophytes and tracheophytes, eventually influencing distinct developmental processes.

The field of crosslinking mass spectrometry has seen substantial advancements over the past decades, enabling the structural analysis of proteins and protein complexes and serving as a powerful tool in protein-protein interaction studies. However, data analysis of large non-cleavable crosslink studies is still a mostly unsolved problem due to its n-squared complexity. We here introduce an algorithm for the identification of non-cleavable crosslinks implemented in our crosslinking search engine MS Annika that is based on sparse matrix multiplication and allows for proteome-wide searches on commodity hardware. We compare our algorithm to other state-of-the-art crosslinking search engines commonly used in the field and conclude that MS Annika unifies high sensitivity, accurate FDR estimation and computational performance, outperforming competing tools. Application of this algorithm enabled us to employ a proteome-wide search of C. elegans nuclei samples, where we were able to uncover previously unknown protein interactions and conclude a comprehensive structural analysis that provides a detailed view of the Box C/D complex. Moreover, our algorithm will enable researchers to conduct similar studies that were previously unfeasible.

Transposable elements (TEs) are repressed in plants through transcriptional gene silencing (TGS), maintained epigenetic silencing marks such as DNA methylation. However, the mechanisms by which silencing is first installed remain poorly understood in plants. Small interfering (si)RNAs and post-transcriptional gene silencing (PTGS) are believed to mediate the initiation of TGS by guiding the first deposition of DNA methylation. To determine how this silencing installation works, we took advantage of ÉVADÉ (EVD), an endogenous retroelement in Arabidopsis, able to recapitulate true de novo silencing with a sequence of PTGS followed by a TGS. To test whether PTGS is required for TGS, we introduce active EVD into RNA-DEPENDENT-RNA-POLYMERASE-6 (RDR6) mutants, an essential PTGS component. EVD activity and silencing are monitored across several generations. In the absence of PTGS, silencing of EVD is still achieved through installation of RNA-directed DNA methylation (RdDM). Our study shows that PTGS is dispensable for de novo EVD silencing. Although we cannot rule out that PTGS might facilitate TGS, or control TE activity, initiation of epigenetic silencing can take place in its absence.

Meristems are three-dimensional (3D) generative structures that contain stem cells and produce new organs and tissues. Meristems develop in all land plants; however we know little about the spatial and temporal regulation of meristem structure in lineages such as bryophytes. Here, we describe the 3D meristem anatomy during the development of the liverwort Marchantia polymorpha. We show that the apical stem cell of the mature meristem is sub-apical, ventral, and in the outer cell layer. Mature meristem anatomy is therefore asymmetrical in the dorsoventral axis, which is reflected by the domain-specific protein localisation of Class III and Class IV Homeodomain-Leucine-Zippers (MpC3HDZ and MpC4HDZ), and by the promoter activity of MpYUCCA2. The dorsoventral asymmetry that defines the mature meristem is absent in the juvenile meristems of asexual propagules known as gemmae. We discovered that anatomical dorsoventral asymmetry of the meristem forms after 1 to 2 days of gemmaling growth, and is accompanied by expression of the dorsal identity reporter MpC3HDZ. We conclude that the gemma meristem has arrested development and undergoes anatomical rearrangement to develop the 3D meristem structure of the mature plant.

Glyphosate is a broad-spectrum herbicide that kills most vascular plant weeds but is ineffective against many bryophytes. Glyphosate competitively inhibits the enolpyruvyl transferase enzyme 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS). EPSPS catalyzes the production of 5-enolpyruvylshikimate-3-phosphate (EPSP)-an intermediate in the shikimate pathway-from shikimate-3-phosphate (S3P) and phosphoenolpyruvate (PEP) substrates. Here, we show that Marchantia polymorpha mutants with loss-of-function mutations in a closely related enolpyruvyl transferase, MpMurA, were more sensitive to glyphosate than wild-type controls. Overexpression of MpMurA in M. polymorpha increased glyphosate tolerance, and heterologous expression of the M. polymorpha MurA enzyme in Arabidopsis thaliana conferred glyphosate resistance. Furthermore, we demonstrate that MpMurA catalyzes the production of EPSP from S3P and PEP substrates. These data demonstrate that MpMurA contributes to glyphosate tolerance in M. polymorpha. We speculate that the existence of two independent mechanisms for EPSP synthesis-one canonical and EPSPS-dependent, and the other MurA-dependent-may account for glyphosate tolerance in bryophytes. Alternatively, MurA in bryophytes may bind glyphosate, thereby leaving more unbound EPSPS enzymes available, allowing aromatic amino acid biosynthesis to continue.

Genome stability is significantly influenced by the precise coordination of chromatin complexes that facilitate the loading and eviction of histones from chromatin during replication, transcription, and DNA repair processes. In this study, we investigate the role of the Arabidopsis H3 histone chaperones ANTI-SILENCING FUNCTION 1 (ASF1) and HISTONE REGULATOR A (HIRA) in the maintenance of telomeres and 45S rDNA loci, genomic sites that are particularly susceptible to changes in the chromatin structure. We find that both ASF1 and HIRA are essential for telomere length regulation, as telomeres are significantly shorter in asf1a1b and hira mutants. However, these shorter telomeres remain localized around the nucleolus and exhibit a comparable relative H3 occupancy to the wild type. In addition to regulating telomere length, ASF1 and HIRA contribute to silencing 45S rRNA genes and affect their copy number. Besides, ASF1 supports global heterochromatin maintenance. Our findings also indicate that ASF1 transiently binds to the TELOMERE REPEAT BINDING 1 protein and the N terminus of telomerase in vivo, suggesting a physical link between the ASF1 histone chaperone and the telomere maintenance machinery.

Ustilago maydis is a biotrophic pathogen causing smut disease in maize. It secretes a cocktail of effector proteins, which target different host proteins during its biotrophic stages in the host plant. One such class of proteins we identified previously is TOPLESS (TPL) and TOPLESS-RELATED (TPR) transcriptional corepressors. Here, we screened 297 U. maydis effector candidates for their ability to interact with maize TPL protein RAMOSA 1 ENHANCER LOCUS 2 LIKE 2 (RELK2) and their ability to induce auxin signaling and thereby identified three novel TPL-interacting protein effectors (Tip6, Tip7, and Tip8). Structural modeling and mutational analysis allowed the identification of TPL-interaction motifs of Tip6 and Tip7. In planta interaction between Tip6 and Tip7 with RELK2 occurs mainly in nuclear compartments, whereas Tip8 colocalizes with RELK2 in a compartment outside the nucleus. Overexpression of Tip8 in nonhost plants leads to cell death, indicating recognition of the effector or its activity. By performing infection assays with single and multideletion mutants of U. maydis, we demonstrate a positive role of Tip6 and Tip7 in U. maydis virulence. Transcriptional profiling of maize leaves infected with Tip effector mutants in comparison with SG200 strain suggests Tip effector activities are not merely redundant.

The rates of appearance of new mutations play a central role in evolution. However, mutational processes in natural environments and their relationship with growth rates are largely unknown, particular in tropical ecosystems with high biodiversity. Here, we examined the somatic mutation landscapes of two tropical trees, Shorea laevis (slow-growing) and S. leprosula (fast-growing), in central Borneo, Indonesia. Using newly constructed genomes, we identified a greater number of somatic mutations in tropical trees than in temperate trees. In both species, we observed a linear increase in the number of somatic mutations with physical distance between branches. However, we found that the rate of somatic mutation accumulation per meter of growth was 3.7-fold higher in S. laevis than in S. leprosula. This difference in the somatic mutation rate was scaled with the slower growth rate of S. laevis compared to S. leprosula, resulting in a constant somatic mutation rate per year between the two species. We also found that somatic mutations are neutral within an individual, but those mutations transmitted to the next generation are subject to purifying selection. These findings suggest that somatic mutations accumulate with absolute time and older trees have a greater contribution towards generating genetic variation.

Macrophage elaboration of inflammatory responses is dynamically regulated, shifting from acute induction to delayed suppression during the course of infection. Here, we show that such regulation of inflammation is modulated by dynamic shifts in metabolism. In macrophages exposed to the bacterial product lipopolysaccharide (LPS), an initial induction of protein biosynthesis is followed by compensatory induction of the transcription factor nuclear factor erythroid 2-like 1 (NRF1), leading to increased flux through the ubiquitin proteasome system (UPS). A major target of NRF1-mediated UPS flux is the mitochondrial proteome, and in the absence of NRF1, ubiquitinated mitochondrial proteins accumulate to trigger severe mitochondrial stress. Such mitochondrial stress engages the integrated stress response-ATF4 axis, which limits mitochondrial translation to attenuate mitochondrial stress but amplifies inflammatory responses to augment susceptibility to septic shock. Therefore, NRF1 mediates a dynamic regulation of mitochondrial proteostasis in inflammatory macrophages that contributes to curbing inflammatory responses.

Meiosis is a hallmark of sexual reproduction because it represents the transition from one life cycle to the next and, in animals, meiosis produces gametes. Why meiosis evolved has been debated and most studies have focused on recombination of the parental alleles as the main function of meiosis. However, 40 years ago, Robin Holliday proposed that an essential function of meiosis is to oppose the consequence of successive mitoses that cause cellular aging. Cellular aging results from accumulated defective organelles and proteins and modifications of chromatin in the form of DNA methylation and histone modifications referred to collectively as epigenetic marks. Here, recent findings supporting the hypothesis that meiosis opposes cellular aging are reviewed and placed in the context of the diversity of the life cycles of eukaryotes, including animals, yeast, flowering plants and the bryophyte Marchantia.

Much of what we know about eukaryotic transcription stems from animals and yeast; however, plants evolved separately for over a billion years, leaving ample time for divergence in transcriptional regulation. Here we set out to elucidate fundamental properties of cis-regulatory sequences in plants. Using massively parallel reporter assays across four plant species, we demonstrate the central role of sequences downstream of the transcription start site (TSS) in transcriptional regulation. Unlike animal enhancers that are position independent, plant regulatory elements depend on their position, as altering their location relative to the TSS significantly affects transcription. We highlight the importance of the region downstream of the TSS in regulating transcription by identifying a DNA motif that is conserved across vascular plants and is sufficient to enhance gene expression in a dose-dependent manner. The identification of a large number of position-dependent enhancers points to fundamental differences in gene regulation between plants and animals.

Long-term genetic studies of wild populations are very scarce, but are essential for connecting ecological and population genetics models, and for understanding the dynamics of biodiversity. We present a study of a wild wheat population sampled over a 36-year period at high spatial resolution. We genotyped 832 individuals from regular sampling along transects during the course of the experiment. Genotypes were clustered into ecological microhabitats over scales of tens of metres, and this clustering was remarkably stable over the 36 generations of the study. Simulations show that it is difficult to determine whether this spatial and temporal stability reflects extremely limited dispersal or fine-scale local adaptation to ecological parameters. Using a common-garden experiment, we showed that the genotypes found in distinct microhabitats differ phenotypically. Our results provide a rare insight into the population genetics of a natural population over a long monitoring period.

The multicellular haploid stage of land plants develops from a single haploid cell produced by meiosis - the spore. Starting from a non-polar state, these spores develop polarity, divide asymmetrically and establish the first axis of symmetry. Here, we show that the nucleus migrates from the cell centroid to the basal pole during polarisation of the Marchantia polymorpha spore cell. A microtubule organising centre on the leading edge of the nucleus initiates a microtubule array between the nuclear surface and the cortex at the basal pole. Simultaneously, cortical microtubules disappear from the apical hemisphere but persist in the basal hemisphere. This is accompanied by the formation a dense network of fine actin filaments between the nucleus and the basal pole cortex. Experimental depolymerisation of either microtubules or actin filaments disrupts cellular asymmetry. These data demonstrate that the cytoskeleton reorganises during spore polarisation and controls the directed migration of the nucleus to the basal pole. The presence of the nucleus at the basal pole provides the cellular asymmetry for the asymmetric cell division that establishes the apical-basal axis of the plant.

Microtubule organising centres (MTOCs) are sites of localised microtubule nucleation in eukaryotic cells. Regulation of microtubule dynamics often involves KATANIN (KTN): a microtubule severing enzyme that cuts microtubules to generate new negative ends, leading to catastrophic depolymerisation. In Arabidopsis thaliana, KTN is required for the organisation of microtubules in the cell cortex, preprophase band, mitotic spindle and phragmoplast. However, as angiosperms lack MTOCs, the role of KTN in MTOC formation has yet to be studied in plants. Two unique MTOCs - the polar organisers - form on opposing sides of the preprophase nucleus in liverworts. Here, we show that KTN-mediated microtubule depolymerisation regulates the number and organisation of polar organisers formed in Marchantia polymorpha. Mpktn mutants that lacked KTN function had supernumerary disorganised polar organisers compared with wild type. This was in addition to defects in the microtubule organisation in the cell cortex, preprophase band, mitotic spindle and phragmoplast. These data are consistent with the hypothesis that KTN-mediated microtubule dynamics are required for the de novo formation of MTOCs, a previously unreported function in plants.

Rice (Oryza sativa L.) endosperm is essential to provide nutrients for seed germination and determine grain yield. RNA editing, a post-transcriptional modification essential for plant development, unfortunately, is not fully characterized during rice endosperm development. Here, we perform systematic analyses to characterize RNA editome during rice endosperm development. We find that most editing sites are C-to-U CDS-recoding in mitochondria, leading to increased hydrophobic amino acids and changed structures of mitochondrial proteins. Comparative analysis of RNA editome reveals that CDS-recoding sites present higher editing frequencies with lower variabilities and their resultant recoded amino acids tend to exhibit stronger evolutionary conservation across many land plants. Furthermore, we classify mitochondrial genes into three groups, presenting distinct patterns in terms of CDS-recoding events. Besides, we conduct genome-wide screening to detect pentatricopeptide repeat (PPR) proteins and construct PPR-RNA binding profiles, yielding candidate PPR editing factors related to rice endosperm development. Taken together, our findings provide valuable insights for deciphering fundamental mechanisms of rice endosperm development underlying RNA editing machinery.

Land plant bodies develop from stem cells located in meristems. However, we know little about how meristems initiate from non-meristematic cells. The haploid body of bryophytes develops from unicellular spores in isolation from the parental plant, which allows all stages of development to be observed. We discovered that the Marchantia spore undergoes a series of reproducibly oriented cell divisions to generate a flat prothallus on which a meristem later develops de novo. The young sporeling comprises an early cell mass. One cell of the early cell mass elongates and undergoes a formative division that produces the prothalloblast, which initiates prothallus formation. A symmetric division of the prothalloblast followed by two transverse divisions generates a four-celled plate that expands into a flat disc through oblique divisions in three of the four plate-cell-derived quadrants. One quadrant gives rise to a flat flabellum. A notch with a meristem and apical stem cell develops at the margin of the flabellum. The transcription factor Marchantia class III homeodomain-leucine-zipper (MpC3HDZ) is a marker of the first flat prothallus structure and polarizes to the dorsal tissues of flabella and meristems. Mpc3hdz mutants are defective in setting up dorsoventrality and thallus body flatness. We report how a regular set of cell divisions forms the prothallus-the first dorsoventral structure-and how cells on the margin of the prothallus develop a dorsoventralized meristem de novo.

The subgenus Tillandsia (Bromeliaceae) belongs to one of the fastest radiating clades in the plant kingdom and is characterized by the repeated evolution of Crassulacean acid metabolism (CAM). Despite its complex genetic basis, this water-conserving trait has evolved independently across many plant families and is regarded as a key innovation trait and driver of ecological diversification in Bromeliaceae. By producing high-quality genome assemblies of a Tillandsia species pair displaying divergent photosynthetic phenotypes, and combining genome-wide investigations of synteny, transposable element (TE) dynamics, sequence evolution, gene family evolution, and temporal differential expression, we were able to pinpoint the genomic drivers of CAM evolution in Tillandsia. Several large-scale rearrangements associated with karyotype changes between the 2 genomes and a highly dynamic TE landscape shaped the genomes of Tillandsia. However, our analyses show that rewiring of photosynthetic metabolism is mainly obtained through regulatory evolution rather than coding sequence evolution, as CAM-related genes are differentially expressed across a 24-h cycle between the 2 species but are not candidates of positive selection. Gene orthology analyses reveal that CAM-related gene families manifesting differential expression underwent accelerated gene family expansion in the constitutive CAM species, further supporting the view of gene family evolution as a driver of CAM evolution.

DNA methylation is environment-sensitive and can mediate stress responses. In trees, changes in the environment might cumulatively shape the methylome landscape over time. However, because high-resolution methylome studies usually focus on single environmental cues, the stress-specificity and long-term stability of methylation responses remain unclear. Here, we studied the methylome plasticity of a Populus nigra cv. 'Italica' clone widely distributed across Europe. Adult trees from different geographic locations were clonally propagated in a common garden experiment and exposed to cold, heat, drought, herbivory, rust infection, and salicylic acid treatments. Whole-genome bisulfite sequencing revealed stress-induced and naturally occurring DNA methylation variants. In CG/CHG contexts, the same genomic regions were often affected by multiple stresses, suggesting a generic methylome response. Moreover, these variants showed striking overlap with naturally occurring methylation variants between trees from different locations. Drought treatment triggered CHH hypermethylation of transposable elements, affecting entire superfamilies near drought-responsive genes. Thus, we revealed genomic hotspots of methylation change that are not stress-specific and that contribute to natural DNA methylation variation, and identified stress-specific hypermethylation of entire transposon superfamilies with possible functional consequences. Our results underscore the importance of studying multiple stressors in a single experiment for recognizing general versus stress-specific methylome responses.

Cells harbour numerous receptor pathways to respond to diverse stimuli, yet often share common downstream signalling components. Mitogen-activated protein kinase (MPK) cascades are an example of such common hubs in eukaryotes. How such common hubs faithfully transduce distinct signals within the same cell-type is insufficiently understood, yet of fundamental importance for signal integration and processing in plants. We engineered a unique genetic background allowing direct comparisons of a developmental and an immunity pathway in one cell-type, the Arabidopsis root endodermis. We demonstrate that the two pathways maintain distinct functional and transcriptional outputs despite common MPK activity patterns. Nevertheless, activation of different MPK kinases and MPK classes led to distinct functional readouts, matching observed pathway-specific readouts. On the basis of our comprehensive analysis of core MPK signalling elements, we propose that combinatorial activation within the MPK cascade determines the differential regulation of an endodermal master transcription factor, MYB36, that drives pathway-specific gene activation.

Natural populations are subject to selection caused by a range of biotic and abiotic factors in their native habitats. Identifying these agents of selection and quantifying their effects is key to understanding how populations adapt to local conditions. We performed a factorial reciprocal-transplant experiment using locally adapted ecotypes of Arabidopsis thaliana at their native sites to distinguish the contributions of adaptation to soil type and climate. Overall adaptive differentiation was strong at both sites. However, we found only very small differences in the strength of selection on local and non-local soil, and adaptation to soil type at most constituted only a few per cent of overall adaptive differentiation. These results indicate that local climatic conditions rather than soil type are the primary driver of adaptive differentiation between these ecotypes.

Circadian clocks are biochemical timers regulating many physiological and molecular processes according to the day/night cycles. The function of the oscillator relies on negative transcriptional/translational feedback loops operated by the so-called clock genes and the encoded clock proteins. Previously, we identified the small GTPase LIGHT INSENSITIVE PERIOD 1 (LIP1) as a circadian-clock-associated protein that regulates light input to the clock in the model plant Arabidopsis thaliana. We showed that LIP1 is also required for suppressing red and blue light-mediated photomorphogenesis, pavement cell shape determination and tolerance to salt stress. Here, we demonstrate that LIP1 is present in a complex of clock proteins GIGANTEA (GI), ZEITLUPE (ZTL) and TIMING OF CAB 1 (TOC1). LIP1 participates in this complex via GUANINE EX-CHANGE FACTOR 7. Analysis of genetic interactions proved that LIP1 affects the oscillator via modulating the function of GI. We show that LIP1 and GI independently and additively regulate photomorphogenesis and salt stress responses, whereas controlling cell shape and photoperiodic flowering are not shared functions of LIP1 and GI. Collectively, our results suggest that LIP1 affects a specific function of GI, possibly by altering binding of GI to downstream signalling components.

Plants continuously remodel and degrade their organelles due to damage from their metabolic activities and environmental stressors, as well as an integral part of their cell differentiation programs. Whereas certain organelles use local hydrolytic enzymes for limited remodeling, most of the pathways that control the partial or complete dismantling of organelles rely on vacuolar degradation. Specifically, selective autophagic pathways play a crucial role in recognizing and sorting plant organelle cargo for vacuolar clearance, especially under cellular stress conditions induced by factors like heat, drought, and damaging light. In these short reviews, we discuss the mechanisms that control the vacuolar degradation of chloroplasts, mitochondria, endoplasmic reticulum, Golgi, and peroxisomes, with an emphasis on autophagy, recently discovered selective autophagy receptors for plant organelles, and crosstalk with other catabolic pathways.

The variability of proteins at the sequence level creates an enormous potential for proteome complexity. Exploring the depths and limits of this complexity is an ongoing goal in biology. Here, we systematically survey human and plant high-throughput bottom-up native proteomics data for protein truncation variants, where substantial regions of the full-length protein are missing from an observed protein product. In humans, Arabidopsis, and the green alga Chlamydomonas, approximately one percent of observed proteins show a short form, which we can assign by comparison to RNA isoforms as either likely deriving from transcript-directed processes or limited proteolysis. While some detected protein fragments align with known splice forms and protein cleavage events, multiple examples are previously undescribed, such as our observation of fibrocystin proteolysis and nuclear translocation in a green alga. We find that truncations occur almost entirely between structured protein domains, even when short forms are derived from transcript variants. Intriguingly, multiple endogenous protein truncations of phase-separating translational proteins resemble cleaved proteoforms produced by enteroviruses during infection. Some truncated proteins are also observed in both humans and plants, suggesting that they date to the last eukaryotic common ancestor. Finally, we describe novel proteoform-specific protein complexes, where the loss of a domain may accompany complex formation.

The early development of the endosperm is crucial for balancing the allocation of maternal nutrients to offspring. This process is believed to be evolutionarily associated with genomic imprinting, resulting in parentally biased allelic gene expression. Beyond FertilizationIndependentSeed (FIS) genes, the number of imprinted genes involved in early endosperm development and seed size determination remains limited. This study introduces early endosperm-expressed HAIKU (IKU) downstream Candidate F-box 1 (ICF1) and ICF2 as maternally expressed imprinted genes (MEGs) in Arabidopsis thaliana. Although these genes are also demethylated by DEMETER (DME) in the central cell, their activation differs from the direct DME-mediated activation seen in classical MEGs such as the FIS genes. Instead, ICF maternal alleles carry pre-established hypomethylation in their promoters, priming them for activation by the WRKY10 transcription factor in the endosperm. On the contrary, paternal alleles are predominantly suppressed by CG methylation. Furthermore, we find that ICF genes partially contribute to the small seed size observed in iku mutants. Our discovery reveals a two-step regulatory mechanism that highlights the important role of conventional transcription factors in the activation of imprinted genes, which was previously not fully recognized. Therefore, the mechanism provides a new dimension to understand the transcriptional regulation of imprinting in plant reproduction and development.

Root exudates contain specialised metabolites that shape the plant's root microbiome. How host-specific microbes cope with these bioactive compounds, and how this ability affects root microbiomes, remains largely unknown. We investigated how maize root bacteria metabolise benzoxazinoids, the main specialised metabolites of maize. Diverse and abundant bacteria metabolised the major compound in the maize rhizosphere MBOA (6-methoxybenzoxazolin-2(3H)-one) and formed AMPO (2-amino-7-methoxy-phenoxazin-3-one). AMPO forming bacteria were enriched in the rhizosphere of benzoxazinoid-producing maize and could use MBOA as carbon source. We identified a gene cluster associated with AMPO formation in microbacteria. The first gene in this cluster, bxdA encodes a lactonase that converts MBOA to AMPO in vitro. A deletion mutant of the homologous bxdA genes in the genus Sphingobium, did not form AMPO nor was it able to use MBOA as a carbon source. BxdA was identified in different genera of maize root bacteria. Here we show that plant-specialised metabolites select for metabolisation-competent root bacteria. BxdA represents a benzoxazinoid metabolisation gene whose carriers successfully colonize the maize rhizosphere and thereby shape the plant's chemical environmental footprint.

The centromere is an essential chromosome region where the kinetochore is formed to control equal chromosome distribution during cell division. The centromere-specific histone H3 variant CENH3 (also called CENP-A) is a prerequisite for the kinetochore formation. Since CENH3 evolves rapidly, associated factors, including histone chaperones mediating the deposition of CENH3 on the centromere, are thought to act through species-specific amino acid sequences. The functions and interaction networks of CENH3 and histone chaperons have been well-characterized in animals and yeasts. However, molecular mechanisms involved in recognition and deposition of CENH3 are still unclear in plants. Here, we used a swapping strategy between domains of CENH3 of Arabidopsis thaliana and the liverwort Marchantia polymorpha to identify specific regions of CENH3 involved in targeting the centromeres and interacting with the general histone H3 chaperone, nuclear autoantigenic sperm protein (NASP). CENH3's LoopN-α1 region was necessary and sufficient for the centromere targeting in cooperation with the α2 region and was involved in interaction with NASP in cooperation with αN, suggesting a species-specific CENH3 recognition. In addition, by generating an Arabidopsis nasp knock-out mutant in the background of a fully fertile GFP-CENH3/cenh3-1 line, we found that NASP was implicated for de novo CENH3 deposition after fertilization and thus for early embryo development. Our results imply that the NASP mediates the supply of CENH3 in the context of the rapidly evolving centromere identity in land plants.

Nucleosomes are DNA-protein complexes composed of histone proteins that form the basis of eukaryotic chromatin. The nucleosome was a key innovation during eukaryotic evolution, but its origin from histone homologues in Archaea remains unclear. Viral histone repeats, consisting of multiple histone paralogues within a single protein, may reflect an intermediate state. Here we examine the diversity of histones encoded by Nucleocytoviricota viruses. We identified 258 histones from 168 viral metagenomes with variable domain configurations including histone singlets, doublets, triplets and quadruplets, the latter comprising the four core histones arranged in series. Viral histone repeats branch phylogenetically between Archaea and eukaryotes and display intermediate functions in Escherichia coli, self-assembling into eukaryotic-like nucleosomes that stack into archaeal-like oligomers capable of impacting genomic activity and condensing DNA. Histone linkage also facilitates nucleosome formation, promoting eukaryotic histone assembly in E. coli. These data support the hypothesis that viral histone repeats originated in stem-eukaryotes and that nucleosome evolution proceeded through histone repeat intermediates.

The histone H2A variant H2A.W occupies transposons and thus prevents access to them in Arabidopsis thaliana. H2A.W is deposited by the chromatin remodeler DDM1, which also promotes the accessibility of chromatin writers to heterochromatin by an unknown mechanism. To shed light on this question, we solve the cryo-EM structures of nucleosomes containing H2A and H2A.W, and the DDM1-H2A.W nucleosome complex. These structures show that the DNA end flexibility of the H2A nucleosome is higher than that of the H2A.W nucleosome. In the DDM1-H2A.W nucleosome complex, DDM1 binds to the N-terminal tail of H4 and the nucleosomal DNA and increases the DNA end flexibility of H2A.W nucleosomes. Based on these biochemical and structural results, we propose that DDM1 counters the low accessibility caused by nucleosomes containing H2A.W to enable the maintenance of repressive epigenetic marks on transposons and prevent their activity.

Classic genome-wide association studies (GWAS) look for associations between individual SNPs and phenotypes of interest. With the rapid progress of high-throughput genotyping and phenotyping technologies, GWAS have become increasingly powerful for detecting genetic determinants and their molecular mechanisms underpinning natural phenotypic variation. However, GWAS frequently yield results with neither expected nor promising loci, nor any significant associations. This is often because associations between SNPs and a single phenotype are confounded, for example with the environment, other traits, or complex genetic structures. Such confounding can mask true genotype-phenotype associations, or inflate spurious associations. To address these problems, numerous methods have been developed that go beyond the standard model. Such advanced GWAS models are flexible and can offer improved statistical power for understanding the genetics underlying complex traits. Despite this advantage, these models have not been widely adopted and implemented compared to the standard GWAS approach, partly because this literature is diverse and often technical. In this review, our aim is to provide an overview of the application and the benefits of various advanced GWAS models for handling complex traits and genetic structures, targeting plant biologists who wish to carry out GWAS more effectively.

Secondary dormancy is an adaptive trait that increases reproductive success by aligning seed germination with permissive conditions for seedling establishment. Aethionema arabicum is an annual plant and member of the Brassicaceae that grows in environments characterized by hot and dry summers. Aethionema arabicum seeds may germinate in early spring when seedling establishment is permissible. We demonstrate that long-day light regimes induce secondary dormancy in the seeds of Aethionema arabicum (CYP accession), repressing germination in summer when seedling establishment is riskier. Characterization of mutants screened for defective secondary dormancy demonstrated that RGL2 mediates repression of genes involved in gibberellin (GA) signaling. Exposure to high temperature alleviates secondary dormancy, restoring germination potential. These data are consistent with the hypothesis that long-day-induced secondary dormancy and its alleviation by high temperatures may be part of an adaptive response limiting germination to conditions permissive for seedling establishment in spring and autumn.

Plants in habitats with unpredictable conditions often have diversified bet-hedging strategies that ensure fitness over a wider range of variable environmental factors. A striking example is the diaspore (seed and fruit) heteromorphism that evolved to maximize species survival in Aethionema arabicum (Brassicaceae) in which external and endogenous triggers allow the production of two distinct diaspores on the same plant. Using this dimorphic diaspore model, we identified contrasting molecular, biophysical, and ecophysiological mechanisms in the germination responses to different temperatures of the mucilaginous seeds (M+ seed morphs), the dispersed indehiscent fruits (IND fruit morphs), and the bare non-mucilaginous M- seeds obtained by pericarp (fruit coat) removal from IND fruits. Large-scale comparative transcriptome and hormone analyses of M+ seeds, IND fruits, and M- seeds provided comprehensive datasets for their distinct thermal responses. Morph-specific differences in co-expressed gene modules in seeds, as well as in seed and pericarp hormone contents, identified a role of the IND pericarp in imposing coat dormancy by generating hypoxia affecting abscisic acid (ABA) sensitivity. This involved expression of morph-specific transcription factors, hypoxia response, and cell wall remodeling genes, as well as altered ABA metabolism, transport, and signaling. Parental temperature affected ABA contents and ABA-related gene expression and altered IND pericarp biomechanical properties. Elucidating the molecular framework underlying the diaspore heteromorphism can provide insight into developmental responses to globally changing temperatures.

All aerial organs in plants originate from the shoot apical meristem, a specialized tissue at the tip of a plant, enclosing a few stem cells. Understanding developmental dynamics within this tissue in relation to internal and external stimuli is of crucial importance. Imaging the meristem at the cellular level beyond very early stages requires the apex to be detached from the plant body, a procedure that does not allow studies in living, intact plants over longer periods. This protocol describes a new confocal microscopy method with the potential to image the shoot apical meristem of an intact, soil-grown, flowering Arabidopsis plant over several days. The setup opens new avenues to study apical stem cells, their interconnection with the whole plant, and their responses to environmental stimuli. Key features • Novel dissection and imaging method of the shoot apical meristem of Arabidopsis. • Procedure performed with intact, soil-grown, flowering plants. • Possibility of long-term live imaging of the shoot apical meristem. • Protocol can be adapted to different plant species.

The abscission of floral organs and emergence of lateral roots in Arabidopsis is regulated by the peptide ligand inflorescence deficient in abscission (IDA) and the receptor protein kinases HAESA (HAE) and HAESA-like 2 (HSL2). During these cell separation processes, the plant induces defense-associated genes to protect against pathogen invasion. However, the molecular coordination between abscission and immunity has not been thoroughly explored. Here, we show that IDA induces a release of cytosolic calcium ions (Ca2+) and apoplastic production of reactive oxygen species, which are signatures of early defense responses. In addition, we find that IDA promotes late defense responses by the transcriptional upregulation of genes known to be involved in immunity. When comparing the IDA induced early immune responses to known immune responses, such as those elicited by flagellin22 treatment, we observe both similarities and differences. We propose a molecular mechanism by which IDA promotes signatures of an immune response in cells destined for separation to guard them from pathogen attack.

The Diphyllatea (CRuMs) are heterotrophic protists currently divided into three distinct clades (Diphy I-III). Diphy I are biflagellates in the genus Diphylleia, whereas Diphy II and III represent cryptic clades comprising Collodictyon-type quadriflagellates that were recently distinguished based on rRNA gene phylogenies. Here, we isolated Diphyllatea from freshwater crater lakes on two South Pacific islands and generated high-quality transcriptomes from species representing each clade, including the first transcriptomic data from Diphy III. Phylogenomic analyses support the separation of Diphy II and III, while transcriptome completeness highlights the utility of these data for future studies. Lastly, we discuss the biogeography and ecology of Diphyllatea on these remote islands.

In recent years, single-cell proteomics (SCP) has become a valuable addition to other single-cell omics technologies for studying cellular heterogeneity. The amount of protein in a single cell is very limited, and in contrast to sequencing techniques, there are currently no means for protein amplification. Therefore, most single-cell proteomics approaches aim to maximize sample preparation efficiency while minimizing peptide loss. By reducing processing volumes to sub-microliters and avoiding manual transfer steps that could lead to peptide loss, peptide recovery, and the robustness of SCP workflows have been significantly improved. In this chapter, we describe a protocol for label-free SCP sample preparation using the cellenONE® platform and the proteoCHIP LF 48 substrate prior to analysis with high-performance liquid chromatography-mass spectrometry.

Modern life science research is a collaborative effort. Few research groups can single-handedly support the necessary equipment, expertise and personnel needed for the ever-expanding portfolio of technologies that are required across multiple disciplines in today's life science endeavours. Thus, research institutes are increasingly setting up scientific core facilities to provide access and specialised support for cutting-edge technologies. Maintaining the momentum needed to carry out leading research while ensuring high-quality daily operations is an ongoing challenge, regardless of the resources allocated to establish such facilities. Here, we outline and discuss the range of activities required to keep things running once a scientific imaging core facility has been established. These include managing a wide range of equipment and users, handling repairs and service contracts, planning for equipment upgrades, renewals, or decommissioning, and continuously upskilling while balancing innovation and consolidation.

Plants rely on autophagy and membrane trafficking to tolerate stress, combat infections, and maintain cellular homeostasis. However, the molecular interplay between autophagy and membrane trafficking is poorly understood. Using an AI-assisted approach, we identified Rab3GAP-like (Rab3GAPL) as a key membrane trafficking node that controls plant autophagy negatively. Rab3GAPL suppresses autophagy by binding to ATG8, the core autophagy adaptor, and deactivating Rab8a, a small GTPase essential for autophagosome formation and defense-related secretion. Rab3GAPL reduces autophagic flux in three model plant species, suggesting that its negative regulatory role in autophagy is conserved in land plants. Beyond autophagy regulation, Rab3GAPL modulates focal immunity against the oomycete pathogen Phytophthora infestans by preventing defense-related secretion. Altogether, our results suggest that Rab3GAPL acts as a molecular rheostat to coordinate autophagic flux and defense-related secretion by restraining Rab8a-mediated trafficking. This unprecedented interplay between a RabGAP-Rab pair and ATG8 sheds new light on the intricate membrane transport mechanisms underlying plant autophagy and immunity.

Plant morphogenesis relies exclusively on oriented cell expansion and division. Nonetheless, the mechanism(s) determining division plane orientation remain elusive. Here, we studied tissue healing after laser-assisted wounding in roots of Arabidopsis thaliana and uncovered how mechanical forces stabilize and reorient the microtubule cytoskeleton for the orientation of cell division. We identified that root tissue functions as an interconnected cell matrix, with a radial gradient of tissue extendibility causing predictable tissue deformation after wounding. This deformation causes instant redirection of expansion in the surrounding cells and reorientation of microtubule arrays, ultimately predicting cell division orientation. Microtubules are destabilized under low tension, whereas stretching of cells, either through wounding or external aspiration, immediately induces their polymerization. The higher microtubule abundance in the stretched cell parts leads to the reorientation of microtubule arrays and, ultimately, informs cell division planes. This provides a long-sought mechanism for flexible re-arrangement of cell divisions by mechanical forces for tissue reconstruction and plant architecture.

Evolutionarily conserved protein associated with topoisomerase II (PAT1) proteins activate mRNA decay through binding mRNA and recruiting decapping factors to optimize posttranscriptional reprogramming. Here, we generated multiple mutants of pat1, pat1 homolog 1 (path1), and pat1 homolog 2 (path2) and discovered that pat triple mutants exhibit extremely stunted growth and all mutants with pat1 exhibit leaf serration while mutants with pat1 and path1 display short petioles. All three PATs can be found localized to processing bodies and all PATs can target ASYMMETRIC LEAVES 2-LIKE 9 transcripts for decay to finely regulate apical hook and lateral root development. In conclusion, PATs exhibit both specific and redundant functions during different plant growth stages and our observations underpin the selective regulation of the mRNA decay machinery for proper development.

Cell polarity is used to guide asymmetric divisions and create morphologically diverse cells. We find that two oppositely oriented cortical polarity domains present during the asymmetric divisions in the Arabidopsis stomatal lineage are reconfigured into polar domains marking ventral (pore-forming) and outward-facing domains of maturing stomatal guard cells. Proteins that define these opposing polarity domains were used as baits in miniTurboID-based proximity labeling. Among differentially enriched proteins, we find kinases, putative microtubule-interacting proteins, and polar SOSEKIs with their effector ANGUSTIFOLIA. Using AI-facilitated protein structure prediction models, we identify potential protein-protein interaction interfaces among them. Functional and localization analyses of the polarity protein OPL2 and its putative interaction partners suggest a positive interaction with mitotic microtubules and a role in cytokinesis. This combination of proteomics and structural modeling with live-cell imaging provides insights into how polarity is rewired in different cell types and cell-cycle stages.

BACKGROUND: Despite its implications for population dynamics and evolution, the relationship between genetic and phenotypic variation in wild populations remains unclear. Here, we estimated variation and plasticity in life-history traits and fitness of the annual plant Arabidopsis thaliana in two common garden experiments that differed in environmental conditions. We used up to 306 maternal inbred lines from six Iberian populations characterized by low and high genotypic (based on whole-genome sequences) and ecological (vegetation type) diversity. RESULTS: Low and high genotypic and ecological diversity was found in edge and core Iberian environments, respectively. Given that selection is expected to be stronger in edge environments and that ecological diversity may enhance both phenotypic variation and plasticity, we expected genotypic diversity to be positively associated with phenotypic variation and plasticity. However, maternal lines, irrespective of the genotypic and ecological diversity of their population of origin, exhibited a substantial amount of phenotypic variation and plasticity for all traits. Furthermore, all populations harbored maternal lines with canalization (robustness) or sensitivity in response to harsher environmental conditions in one of the two experiments. CONCLUSIONS: Overall, we conclude that the environmental attributes of each population probably determine their genotypic diversity, but all populations maintain substantial phenotypic variation and plasticity for all traits, which represents an asset to endure in changing environments.

Induction of DNA damage triggers rapid phosphorylation of the histone H2A.X (γH2A.X). In animals, mediator of DNA damage checkpoint 1 (MDC1) binds γH2A.X through a tandem BRCA1 carboxyl-terminal (tBRCT) domain and mediates recruitment of downstream effectors of DNA damage response (DDR). However, readers of this modification in plants have remained elusive. We show that from the Arabidopsis BRCT domain proteome, BCP1-4 proteins with tBRCT domains are involved in DDR. Through its tBRCT domain BCP4 binds γH2A.X in vitro and localizes to DNA damage-induced foci in an H2A.X-dependent manner. BCP4 also contains a domain that interacts directly with NBS1 and thus acts as a functional counterpart of MDC1. We also show that BCP1, that contains two tBRCT domains, co-localizes with γH2A.X but it does not bind γH2A.X suggesting functional similarity with human PAXIP1. A phylogenetic analysis supports that PAXIP1 and MDC1 in metazoa and their plant counterparts evolved independently from common ancestors with tBRCT domains. Collectively, our study reveals missing components and provides mechanistic and evolutionary insights into plant DDR.

In sexually propagating organisms, genetic, and epigenetic mutations are evolutionarily relevant only if they occur in the germline and are hence transmitted to the next generation. In contrast to most animals, plants are considered to lack an early segregating germline, implying that somatic cells can contribute genetic information to progeny. Here we demonstrate that 2 ARGONAUTE proteins, AGO5 and AGO9, mark cells associated with sexual reproduction in Arabidopsis (Arabidopsis thaliana) throughout development. Both AGOs are loaded with dynamically changing small RNA populations derived from highly methylated, pericentromeric, long transposons. Sequencing of single stem cell nuclei revealed that many of these transposons are co-expressed within an AGO5/9 expression domain in the shoot apical meristem (SAM). Co-occurrence of transposon expression and specific ARGONAUTE (AGO) expression in the SAM is reminiscent of germline features in animals and supports the existence of an early segregating germline in plants. Our results open the path to investigating transposon biology and epigenome dynamics at cellular resolution in the SAM stem cell niche.

The transition from germinating seeds to emerging seedlings is one of the most vulnerable plant life cycle stages. Heteromorphic diaspores (seed and fruit dispersal units) are an adaptive bet-hedging strategy to cope with spatiotemporally variable environments. While the roles and mechanisms of seedling traits have been studied in monomorphic species, which produce one type of diaspore, very little is known about seedlings in heteromorphic species. Using the dimorphic diaspore model Aethionema arabicum (Brassicaceae), we identified contrasting mechanisms in the germination responses to different temperatures of the mucilaginous seeds (M+ seed morphs), the dispersed indehiscent fruits (IND fruit morphs), and the bare non-mucilaginous M- seeds obtained from IND fruits by pericarp (fruit coat) removal. What follows the completion of germination is the pre-emergence seedling growth phase, which we investigated by comparative growth assays of early seedlings derived from the M+ seeds, bare M- seeds, and IND fruits. The dimorphic seedlings derived from M+ and M- seeds did not differ in their responses to ambient temperature and water potential. The phenotype of seedlings derived from IND fruits differed in that they had bent hypocotyls and their shoot and root growth was slower, but the biomechanical hypocotyl properties of 15-day-old seedlings did not differ between seedlings derived from germinated M+ seeds, M- seeds, or IND fruits. Comparison of the transcriptomes of the natural dimorphic diaspores, M+ seeds and IND fruits, identified 2,682 differentially expressed genes (DEGs) during late germination. During the subsequent 3 days of seedling pre-emergence growth, the number of DEGs was reduced 10-fold to 277 root DEGs and 16-fold to 164 shoot DEGs. Among the DEGs in early seedlings were hormonal regulators, in particular for auxin, ethylene, and gibberellins. Furthermore, DEGs were identified for water and ion transporters, nitrate transporter and assimilation enzymes, and cell wall remodeling protein genes encoding enzymes targeting xyloglucan and pectin. We conclude that the transcriptomes of seedlings derived from the dimorphic diaspores, M+ seeds and IND fruits, undergo transcriptional resetting during the post-germination pre-emergence growth transition phase from germinated diaspores to growing seedlings.

Due to the accelerating climate change, it is crucial to understand how plants adapt to rapid environmental changes. Such adaptation may be mediated by epigenetic mechanisms like DNA methylation, which could heritably alter phenotypes without changing the DNA sequence, especially across clonal generations. However, we are still missing robust evidence of the adaptive potential of DNA methylation in wild clonal populations. Here, we studied genetic, epigenetic and transcriptomic variation of Fragaria vesca, a predominantly clonally reproducing herb. We examined samples from 21 natural populations across three climatically distinct geographic regions, as well as clones of the same individuals grown in a common garden. We found that epigenetic variation was partly associated with climate of origin, particularly in non-CG contexts. Importantly, a large proportion of this variation was heritable across clonal generations. Additionally, a subset of these epigenetic changes affected the expression of genes mainly involved in plant growth and responses to pathogen and abiotic stress. These findings highlight the potential influence of epigenetic changes on phenotypic traits. Our findings indicate that variation in DNA methylation, which can be environmentally inducible and heritable, may enable clonal plant populations to adjust to their environmental conditions even in the absence of genetic adaptation.

We are becoming aware of a growing number of organisms that do not express genetic information equally from both parents as a result of an epigenetic phenomenon called genomic imprinting. Recently, it was shown that the entire paternal genome is repressed during the diploid phase of the life cycle of the liverwort Marchantia polymorpha. The deposition of the repressive epigenetic mark H3K27me3 on the male pronucleus is responsible for the imprinted state, which is reset by the end of meiosis. Here, we put these recent reports in perspective of other forms of imprinting and discuss the potential mechanisms of imprinting in bryophytes and the causes of its evolution.

Comprehensive proteomic analysis is essential to elucidate molecular pathways and protein functions. Despite tremendous progress in proteomics, current studies still suffer from limited proteomic coverage and dynamic range. Here, we utilize micropillar array columns (µPACs) together with wide-window acquisition and the AI-based CHIMERYS search engine to achieve excellent proteomic comprehensiveness for bulk proteomics, affinity purification mass spectrometry and single cell proteomics. Our data show that µPACs identify ≤50% more peptides and ≤24% more proteins, while offering improved throughput, which is critical for large (clinical) proteomics studies. Combining wide precursor isolation widths of m/z 4-12 with the CHIMERYS search engine identified +51-74% and +59-150% more proteins and peptides, respectively, for single cell, co-immunoprecipitation, and multi-species samples over a conventional workflow at well-controlled false discovery rates. The workflow further offers excellent precision, with CVs <7% for low input bulk samples, and accuracy, with deviations <10% from expected fold changes for regular abundance two-proteome mixes. Compared to a conventional workflow, our entire optimized platform discovered 92% more potential interactors in a protein-protein interaction study on the chromatin remodeler Smarca5/Snf2h. These include previously described Smarca5 binding partners and undescribed ones including Arid1a, another chromatin remodeler with key roles in neurodevelopmental and malignant disorders.

Plant viruses account for enormous agricultural losses worldwide, and the most effective way to combat them is to identify genetic material conferring plant resistance to these pathogens. Aiming to identify genetic associations with responses to infection, we screened a large panel of Arabidopsis thaliana natural inbred lines for four disease-related traits caused by infection by A. thaliana-naïve and -adapted isolates of the natural pathogen turnip mosaic virus (TuMV). We detected a strong, replicable association in a 1.5 Mb region on chromosome 2 with a 10-fold increase in relative risk of systemic necrosis. The region contains several plausible causal genes as well as abundant structural variation, including an insertion of a Copia transposon into a Toll/interleukin receptor (TIR-NBS-LRR) coding for a gene involved in defense, that could be either a driver or a consequence of the disease-resistance locus. When inoculated with TuMV, loss-of-function mutant plants of this gene exhibited different symptoms than wild-type plants. The direction and severity of symptom differences depended on the adaptation history of the virus. This increase in symptom severity was specific for infections with the adapted isolate. Necrosis-associated alleles are found worldwide, and their distribution is consistent with a trade-off between resistance during viral outbreaks and a cost of resistance otherwise, leading to negative frequency-dependent selection.

This work demonstrates the utility of high-throughput nanoLC-MS and label-free quantification (LFQ) for sample-limited bottom-up proteomics analysis, including single-cell proteomics (SCP). Conditions were optimized on a 50 μm internal diameter (I.D.) column operated at 100 nL/min in the direct injection workflow to balance method sensitivity and sample throughput from 24 to 72 samples/day. Multiple data acquisition strategies were also evaluated for proteome coverage, including data-dependent acquisition (DDA), wide-window acquisition (WWA), and wide-window data-independent acquisition (WW-DIA). Analyzing 250 pg HeLa digest with a 10-min LC gradient (72 samples/day) provided >900, >1,800, and >3,000 protein group identifications for DDA, WWA, and WW-DIA, respectively. Total method cycle time was further reduced from 20 to 14.4 min (100 samples/day) by employing a trap-and-elute workflow, enabling 70% mass spectrometer utilization. The method was applied to library-free DIA analysis of single-cell samples, yielding >1,700 protein groups identified. In conclusion, this study provides a high-sensitivity, high-throughput nanoLC-MS configuration for sample-limited proteomics.

The algal ancestors of land plants underwent a transition from a unicellular to a multicellular body plan.1 This transition likely took place early in streptophyte evolution, sometime after the divergence of the Chlorokybophyceae/Mesostigmatophyceae lineage, but before the divergence of the Klebsormidiophyceae lineage.2 How this transition was brought about is unknown; however, it was likely facilitated by the evolution of novel mechanisms to spatially regulate morphogenesis. In land plants, RHO of plant (ROP) signaling plays a conserved role in regulating polarized cell growth and cell division orientation to orchestrate morphogenesis.3,4,5,6,7,8 ROP constitutes a plant-specific subfamily of the RHO GTPases, which are more widely conserved throughout eukaryotes.9,10 Although the RHO family originated in early eukaryotes,11,12 how and when the ROP subfamily originated had remained elusive. Here, we demonstrate that ROP signaling was established early in the streptophyte lineage, sometime after the divergence of the Chlorokybophyceae/Mesostigmatophyceae lineage, but before the divergence of the Klebsormidiophyceae lineage. This period corresponds to when the unicellular-to-multicellular transition likely took place in the streptophytes. In addition to being critical for the complex morphogenesis of extant land plants, we speculate that ROP signaling contributed to morphological evolution in early streptophytes.

Exocytosis is a conserved trafficking pathway that transports secretory vesicles to the extracellular space, replenishes the plasma membrane and is essential for establishing cell polarity. Its spatiotemporal regulation is mediated by an evolutionary conserved octameric tethering complex, the exocyst. In plants, certain subunits of this complex have diversified and acquired multiple functions, including a central role in defense against pathogens and pests. Here, I review the latest evidence suggesting the dramatic expansion and functional diversification of the exocyst subunit Exo70 is likely driven by a coevolutionary arms race, in which Exo70 proteins are repeatedly targeted by effectors from multiple pathogens and, in turn, are monitored by plant immune receptors for pathogen perception.

Recent studies have shown that correlations between chromatin modifications and transcription vary among eukaryotes. This is the case for marked differences between the chromatin of the moss Physcomitrium patens and the liverwort Marchantia polymorpha. Mosses and liverworts diverged from hornworts, altogether forming the lineage of bryophytes that shared a common ancestor with land plants. We aimed to describe chromatin in hornworts to establish synapomorphies across bryophytes and approach a definition of the ancestral chromatin organization of land plants. We used genomic methods to define the 3D organization of chromatin and map the chromatin landscape of the model hornwort Anthoceros agrestis. We report that nearly half of the hornwort transposons were associated with facultative heterochromatin and euchromatin and formed the center of topologically associated domains delimited by protein coding genes. Transposons were scattered across autosomes, which contrasted with the dense compartments of constitutive heterochromatin surrounding the centromeres in flowering plants. Most of the features observed in hornworts are also present in liverworts or in mosses but are distinct from flowering plants. Hence, the ancestral genome of bryophytes was likely a patchwork of units of euchromatin interspersed within facultative and constitutive heterochromatin. We propose this genome organization was ancestral to land plants.

Multiplexed and label-free mass spectrometry-based approaches with single-cell resolution have attributed surprising heterogeneity to presumed homogenous cell populations. Even though specialized experimental designs and instrumentation have demonstrated remarkable advances, the efficient sample preparation of single cells still lags. Here, we introduce the proteoCHIP, a universal option for single-cell proteomics sample preparation including multiplexed labeling up to 16-plex with high sensitivity and throughput. The automated processing using a commercial system combining single-cell isolation and picoliter dispensing, the cellenONE, reduces final sample volumes to low nanoliters submerged in a hexadecane layer simultaneously eliminating error-prone manual sample handling and overcoming evaporation. The specialized proteoCHIP design allows direct injection of single cells via a standard autosampler resulting in around 1500 protein groups per TMT10-plex with reduced or eliminated need for a carrier proteome. We evaluated the effect of wider precursor isolation windows at single-cell input levels and found that using 2 Da isolation windows increased overall sensitivity without significantly impacting interference. Using the dedicated mass spectrometry acquisition strategies detailed here, we identified on average close to 2000 proteins per TMT10-plex across 170 multiplexed single cells that readily distinguished human cell types. Overall, our workflow combines highly efficient sample preparation, chromatographic and ion mobility-based filtering, rapid wide-window data-dependent acquisition analysis, and intelligent data analysis for optimal multiplexed single-cell proteomics. This versatile and automated proteoCHIP-based sample preparation approach is sufficiently sensitive to drive biological applications of single-cell proteomics and can be readily adopted by proteomics laboratories.

Immunopeptidomics, as the analysis of antigen peptides being presented to the immune system via major histocompatibility complexes (MHC), is being seen as an imperative tool for identifying epitopes for vaccine development to treat cancer and viral and bacterial infections as well as parasites. The field has made tremendous strides over the last 25 years but currently still faces challenges in sensitivity and throughput for widespread applications in personalized medicine and large vaccine development studies. Cutting-edge technological advancements in sample preparation, liquid chromatography as well as mass spectrometry, and data analysis, however, are currently transforming the field. This perspective showcases how the advent of single-cell proteomics has accelerated this transformation of immunopeptidomics in recent years and will pave the way for even more sensitive and higher-throughput immunopeptidomics analyses.

Despite the wide diversity in male sexual development across land plants, new work reveals the conservation of a heterodimer of transcription factors as master regulators of the male germline.

We propose that several chromatin-mediated regulatory processes are dominated by source-sink relationships in which factors operate as 'sources' to produce or provide a resource and compete with each other to occupy separate 'sinks'. In this model, large portions of genomic DNA operate as 'sinks', which are filled by 'sources', such as available histone variants, covalent modifications to histones, the readers of these modifications and non-coding RNAs. Competing occupation for the sinks by different sources leads to distinct states of genomic equilibrium in differentiated cells. During dynamic developmental events, such as sexual reproduction, we propose that dramatic and rapid reconfiguration of source-sink relationships modifies chromatin states. We envision that re-routing of sources could occur by altering the dimensions of the sink, by reconfiguration of existing sink occupation or by varying the size of the source, providing a central mechanism to explain a plethora of epigenetic phenomena, which contribute to phenotypic variegation, zygotic genome activation and nucleolar dominance.

The mobility of transposable elements (TEs) contributes to evolution of genomes. Their uncontrolled activity causes genomic instability; therefore, expression of TEs is silenced by host genomes. TEs are marked with DNA and H3K9 methylation, which are associated with silencing in flowering plants, animals, and fungi. However, in distantly related groups of eukaryotes, TEs are marked by H3K27me3 deposited by the Polycomb repressive complex 2 (PRC2), an epigenetic mark associated with gene silencing in flowering plants and animals. The direct silencing of TEs by PRC2 has so far only been shown in one species of ciliates. To test if PRC2 silences TEs in a broader range of eukaryotes, we generated mutants with reduced PRC2 activity and analyzed the role of PRC2 in extant species along the lineage of Archaeplastida and in the diatom P. tricornutum. In this diatom and the red alga C. merolae, a greater proportion of TEs than genes were repressed by PRC2, whereas a greater proportion of genes than TEs were repressed by PRC2 in bryophytes. In flowering plants, TEs contained potential cis-elements recognized by transcription factors and associated with neighbor genes as transcriptional units repressed by PRC2. Thus, silencing of TEs by PRC2 is observed not only in Archaeplastida but also in diatoms and ciliates, suggesting that PRC2 deposited H3K27me3 to silence TEs in the last common ancestor of eukaryotes. We hypothesize that during the evolution of Archaeplastida, TE fragments marked with H3K27me3 were selected to shape transcriptional regulation, controlling networks of genes regulated by PRC2.

Stem cells are essential for the development and organ regeneration of multicellular organisms, so their infection by pathogenic viruses must be prevented. Accordingly, mammalian stem cells are highly resistant to viral infection due to dedicated antiviral pathways including RNA interference (RNAi). In plants, a small group of stem cells harbored within the shoot apical meristem generate all postembryonic above-ground tissues, including the germline cells. Many viruses do not proliferate in these cells, yet the molecular bases of this exclusion remain only partially understood. Here, we show that a plant-encoded RNA-dependent RNA polymerase, after activation by the plant hormone salicylic acid, amplifies antiviral RNAi in infected tissues. This provides stem cells with RNA-based virus sequence information, which prevents virus proliferation. Furthermore, we find RNAi to be necessary for stem cell exclusion of several unrelated RNA viruses, despite their ability to efficiently suppress RNAi in the rest of the plant. This work elucidates a molecular pathway of great biological and economic relevance and lays the foundations for our future understanding of the unique systems underlying stem cell immunity.

The chromatin of flowering plants exhibits a wide range of sequence variants of the core and linker histones. Recent studies have demonstrated that specific histone variant enrichment, combined with post-translational modifications (PTMs) of histones, defines distinct chromatin states that impact specific chromatin functions. Chromatin remodelers are emerging as key regulators of histone variant dynamics, contributing to shaping chromatin states and regulating gene transcription in response to environment. Recognizing the histone variants by their specific readers, controlled by histone PTMs, is crucial for maintaining genome and chromatin integrity. In addition, various histone variants have been shown to play essential roles in remodeling chromatin domains to facilitate important programmed transitions throughout the plant life cycle. In this review, we discuss recent findings in this exciting field of research, which holds immense promise for many surprising discoveries related to the evolution of complexity in plant organization through a seemingly simple protein family.

Fertilization in Arabidopsis (Arabidopsis thaliana) is a highly coordinated process that begins with a pollen tube delivering the two sperm cells into the embryo sac. Each sperm cell can then fertilize either the egg or the central cell to initiate embryo or endosperm development, respectively. The success of this double fertilization process requires a tight cell cycle synchrony between the male and female gametes to allow karyogamy (nuclei fusion). However, the cell cycle status of the male and female gametes during fertilization remains elusive as DNA quantification and DNA replication assays have given conflicting results. Here, to reconcile these results, we quantified the DNA replication state by DNA sequencing and performed microscopic analyses of fluorescent markers covering all phases of the cell cycle. We show that male and female Arabidopsis gametes are both arrested prior to DNA replication at maturity and initiate their DNA replication only during fertilization.

Engineering the plant immune system offers genetic solutions to mitigate crop diseases caused by diverse agriculturally significant pathogens and pests. Modification of intracellular plant immune receptors of the nucleotide-binding leucine-rich repeat (NLR) receptor superfamily for expanded recognition of pathogen virulence proteins (effectors) is a promising approach for engineering disease resistance. However, engineering can cause NLR autoactivation, resulting in constitutive defense responses that are deleterious to the plant. This may be due to plant NLRs associating in highly complex signaling networks that coevolve together, and changes through breeding or genetic modification can generate incompatible combinations, resulting in autoimmune phenotypes. The sensor and helper NLRs of the rice (Oryza sativa) NLR pair Pik have coevolved, and mismatching between noncoevolved alleles triggers constitutive activation and cell death. This limits the extent to which protein modifications can be used to engineer pathogen recognition and enhance disease resistance mediated by these NLRs. Here, we dissected incompatibility determinants in the Pik pair in Nicotiana benthamiana and found that heavy metal-associated (HMA) domains integrated in Pik-1 not only evolved to bind pathogen effectors but also likely coevolved with other NLR domains to maintain immune homeostasis. This explains why changes in integrated domains can lead to autoactivation. We then used this knowledge to facilitate engineering of new effector recognition specificities, overcoming initial autoimmune penalties. We show that by mismatching alleles of the rice sensor and helper NLRs Pik-1 and Pik-2, we can enable the integration of synthetic domains with novel and enhanced recognition specificities. Taken together, our results reveal a strategy for engineering NLRs, which has the potential to allow an expanded set of integrations and therefore new disease resistance specificities in plants.

An ancient evolutionary innovation of a novel cell type, the stinging cell (cnidocyte), appeared >600 million years ago in the phylum Cnidaria (sea anemones, corals, hydroids, and jellyfish). A complex bursting nano-injector of venom, the cnidocyst, is embedded in cnidocytes and enables cnidarians to paralyze their prey and predators, contributing to this phylum's evolutionary success. In this work, we show that post-transcriptional regulation by a pan-cnidarian microRNA, miR-2022, is essential for biogenesis of these cells in the sea anemone Nematostella vectensis. By manipulation of miR-2022 levels in a transgenic reporter line of cnidocytes, followed by transcriptomics, single-cell data analysis, prey paralysis assays, and cell sorting of transgenic cnidocytes, we reveal that miR-2022 enables cnidocyte biogenesis in Nematostella, while exhibiting a conserved expression domain with its targets in cnidocytes of other cnidarian species. Thus, here we revealed a functional basis to the conservation of one of nature's most ancient microRNAs.

Cell polarity combined with asymmetric cell divisions (ACDs) generates cellular diversity. In the Arabidopsis stomatal lineage, a single cortical polarity domain marked by BASL orients ACDs and is segregated to the larger daughter to enforce cell fate. We discovered a second, oppositely positioned polarity domain defined by OCTOPUS-LIKE (OPL) proteins, which forms prior to ACD and is segregated to the smaller (meristemoid) daughter. Genetic and misexpression analyses show that OPLs promote meristemoid-amplifying divisions and delay stomatal fate progression. Polarity mediates OPL segregation into meristemoids but is not required for OPL function. OPL localization and activity are largely independent of other stomatal polarity genes and of the brassinosteroid signaling components associated with OPLs in other contexts. While OPLs are unique to seed plants, ectopic expression in the liverwort Marchantia suppressed epidermal fate progression, suggesting that OPLs engage ancient and broadly conserved pathways to regulate cell division and cell fate.

Acclimation and adaptation of metabolism to a changing environment are key processes for plant survival and reproductive success. In the present study, 241 natural accessions of Arabidopsis (Arabidopsis thaliana) were grown under two different temperature regimes, 16 °C and 6 °C, and growth parameters were recorded, together with metabolite profiles, to investigate the natural genome × environment effects on metabolome variation. The plasticity of metabolism, which was captured by metabolic distance measures, varied considerably between accessions. Both relative growth rates and metabolic distances were predictable by the underlying natural genetic variation of accessions. Applying machine learning methods, climatic variables of the original growth habitats were tested for their predictive power of natural metabolic variation among accessions. We found specifically habitat temperature during the first quarter of the year to be the best predictor of the plasticity of primary metabolism, indicating habitat temperature as the causal driver of evolutionary cold adaptation processes. Analyses of epigenome- and genome-wide associations revealed accession-specific differential DNA-methylation levels as potentially linked to the metabolome and identified FUMARASE2 as strongly associated with cold adaptation in Arabidopsis accessions. These findings were supported by calculations of the biochemical Jacobian matrix based on variance and covariance of metabolomics data, which revealed that growth under low temperatures most substantially affects the accession-specific plasticity of fumarate and sugar metabolism. Our findings indicate that the plasticity of metabolic regulation is predictable from the genome and epigenome and driven evolutionarily by Arabidopsis growth habitats.

ATP-binding cassette (ABC) transporters actively transport various substances across membranes, while uridine diphosphate (UDP) glycosyltransferases (UGTs) are proteins that catalyse the chemical modification of various organic compounds. Both of these protein superfamilies have been associated with conferring herbicide resistance in weeds. Little is known about the evolutionary history of these protein families in the Archaeplastida. To infer the evolutionary histories of these protein superfamilies, we compared protein sequences collected from 10 species which represent distinct lineages of the Archaeplastida-the lineage including glaucophyte algae, rhodophyte algae, chlorophyte algae and the streptophytes-and generated phylogenetic trees. We show that ABC transporters were present in the last common ancestor of the Archaeplastida which lived 1.6 billion years ago, and the major clades identified in extant plants were already present then. Conversely, we only identified UGTs in members of the streptophyte lineage, which suggests a loss of these proteins in earlier diverging Archaeplastida lineages or arrival of UGTs into a common ancestor of the streptophyte lineage through horizontal gene transfer from a non-Archaeplastida eukaryote lineage. We found that within the streptophyte lineage, most diversification of the UGT protein family occurred in the vascular lineage, with 17 of the 20 clades identified in extant plants present only in vascular plants. Based on our findings, we conclude that ABC transporters and UGTs are ancient protein families which diversified during Archaeplastida evolution, which may have evolved for developmental functions as plants began to occupy new environmental niches and are now being selected to confer resistance to a diverse range of herbicides in weeds.

Members of the NETWORKED (NET) family are involved in actin-membrane interactions. Here we show that two members of the NET family, NET4A and NET4B, are essential for normal guard cell actin reorganization, which is a process critical for stomatal closure in plant immunity. NET4 proteins interact with F-actin and with members of the Rab7 GTPase RABG3 family through two distinct domains, allowing for simultaneous localization to actin filaments and the tonoplast. NET4 proteins interact with GTP-bound, active RABG3 members, suggesting their function being downstream effectors. We also show that RABG3b is critical for stomatal closure induced by microbial patterns. Taken together, we conclude that the actin cytoskeletal remodelling during stomatal closure involves a molecular link between actin filaments and the tonoplast, which is mediated by the NET4-RABG3b interaction. We propose that stomatal closure to microbial patterns involves the coordinated action of immune-triggered osmotic changes and actin cytoskeletal remodelling likely driving compact vacuolar morphologies.

Individual cells give rise to diverse cell lineages during the development of multicellular organisms. Understanding the contribution of these lineages to mature organisms is a central question of developmental biology. Several techniques to document cell lineages have been used, from marking single cells with mutations that express a visible marker to generating molecular bar codes by CRISPR-induced mutations and subsequent single-cell analysis. Here, we exploit the mutagenic activity of CRISPR to allow lineage tracing within living plants with a single reporter. Cas9-induced mutations are directed to correct a frameshift mutation that restores expression of a nuclear fluorescent protein, labelling the initial cell and all progenitor cells with a strong signal without modifying other phenotypes of the plants. Spatial and temporal control of Cas9 activity can be achieved using tissue-specific and/or inducible promoters. We provide proof of principle for the function of lineage tracing in two model plants. The conserved features of the components and the versatile cloning system, allowing for easy exchange of promoters, are expected to make the system widely applicable.

Long non-coding RNAs (lncRNAs) are understudied and underannotated in plants. In mammals, lncRNA loci are nearly as ubiquitous as protein-coding genes, and their expression is highly variable between individuals of the same species. Using Arabidopsis thaliana as a model, we aimed to elucidate the true scope of lncRNA transcription across plants from different regions and study its natural variation. We used transcriptome deep sequencing datasets spanning hundreds of natural accessions and several developmental stages to create a population-wide annotation of lncRNAs, revealing thousands of previously unannotated lncRNA loci. While lncRNA transcription is ubiquitous in the genome, most loci appear to be actively silenced and their expression is extremely variable between natural accessions. This high expression variability is largely caused by the high variability of repressive chromatin levels at lncRNA loci. High variability was particularly common for intergenic lncRNAs (lincRNAs), where pieces of transposable elements (TEs) present in 50% of these lincRNA loci are associated with increased silencing and variation, and such lncRNAs tend to be targeted by the TE silencing machinery. We created a population-wide lncRNA annotation in Arabidopsis and improve our understanding of plant lncRNA genome biology, raising fundamental questions about what causes transcription and silencing across the genome.

Cross-linking mass spectrometry has become a powerful tool for the identification of protein-protein interactions and for gaining insight into the structures of proteins. We previously published MS Annika, a cross-linking search engine which can accurately identify cross-linked peptides in MS2 spectra from a variety of different MS-cleavable cross-linkers. In this publication, we present MS Annika 2.0, an updated version implementing a new search algorithm that, in addition to MS2 level, only supports the processing of data from MS2-MS3-based approaches for the identification of peptides from MS3 spectra, and introduces a novel scoring function for peptides identified across multiple MS stages. Detected cross-links are validated by estimating the false discovery rate (FDR) using a target-decoy approach. We evaluated the MS3-search-capabilities of MS Annika 2.0 on five different datasets covering a variety of experimental approaches and compared it to XlinkX and MaXLinker, two other cross-linking search engines. We show that MS Annika detects up to 4 times more true unique cross-links while simultaneously yielding less false positive hits and therefore a more accurate FDR estimation than the other two search engines. All mass spectrometry proteomics data along with result files have been deposited to the ProteomeXchange consortium via the PRIDE partner repository with the dataset identifier PXD041955.

Efficient repair of DNA lesions is essential for the faithful transmission of genetic information between somatic cells and for genome integrity across generations. Plants have multiple, partially redundant, and overlapping DNA repair pathways, probably due to the less constricted germline and the inevitable exposure to light including higher energy wavelengths. Many proteins involved in DNA repair and their mode of actions are well described. In contrast, a role for DNA damage-associated RNA components, evident from many other organisms, is less well understood. Here, we have challenged young Arabidopsis thaliana plants with two different types of genotoxic stress and performed de novo assembly and transcriptome analysis. We identified three long noncoding RNAs (lncRNAs) that are lowly or not expressed under regular conditions but up-regulated or induced by DNA damage. We generated CRISPR/Cas deletion mutants and found that the absence of the lncRNAs impairs the recovery capacity of the plants from genotoxic stress. The genetic loci are highly conserved among world-wide distributed Arabidopsis accessions and within related species in the Brassicaceae group. Together, these results suggest that the lncRNAs have a conserved function in connection with DNA damage and provide a basis for mechanistic analysis of their role.

Factors including climate change and increased global exchange are set to escalate the prevalence of plant diseases, posing an unprecedented threat to global food security and making it more challenging to meet the demands of an ever-growing population. As such, new methods of pathogen control are essential to help with the growing danger of crop losses to plant diseases. The intracellular immune system of plants utilizes nucleotide-binding leucine-rich repeat (NLR) receptors to recognize and activate defense responses to pathogen virulence proteins (effectors) delivered to the host. Engineering the recognition properties of plant NLRs toward pathogen effectors is a genetic solution to plant diseases with high specificity, and it is more sustainable than several current methods for pathogen control that frequently rely on agrochemicals. Here, we highlight the pioneering approaches toward enhancing effector recognition in plant NLRs and discuss the barriers and solutions in engineering the plant intracellular immune system.

Cell polarity-broadly defined as the asymmetric distribution of cellular activities and subcellular components within a cell-determines the geometry of cell growth and division during development. RHO GTPase proteins regulate the establishment of cell polarity and are conserved among eukaryotes. RHO of plant (ROP) proteins are a subgroup of RHO GTPases that are required for cellular morphogenesis in plants. However, how ROP proteins modulate the geometry of cell growth and division during the morphogenesis of plant tissues and organs is not well understood. To investigate how ROP proteins function during tissue development and organogenesis, we characterized the function of the single-copy ROP gene of the liverwort Marchantia polymorpha (MpROP). M. polymorpha develops morphologically complex three-dimensional tissues and organs exemplified by air chambers and gemmae, respectively. Mprop loss-of-function mutants form defective air chambers and gemmae, indicating ROP function is required for tissue development and organogenesis. During air chamber and gemma development in wild type, the MpROP protein is enriched to sites of polarized growth at the cell surface and accumulates at the expanding cell plate of dividing cells. Consistent with these observations, polarized cell growth is lost and cell divisions are misoriented in Mprop mutants. We propose that ROP regulates both polarized cell growth and cell division orientation in a coordinated manner to orchestrate tissue development and organogenesis in land plants.

How different intrinsic sequence variations and regulatory modifications of histones combine in nucleosomes remain unclear. To test the importance of histone variants in the organization of chromatin we investigated how histone variants and histone modifications assemble in the Arabidopsis thaliana genome. We showed that a limited number of chromatin states divide euchromatin and heterochromatin into several subdomains. We found that histone variants are as significant as histone modifications in determining the composition of chromatin states. Particularly strong associations were observed between H2A variants and specific combinations of histone modifications. To study the role of H2A variants in organizing chromatin states we determined the role of the chromatin remodeler DECREASED IN DNA METHYLATION (DDM1) in the organization of chromatin states. We showed that the loss of DDM1 prevented the exchange of the histone variant H2A.Z to H2A.W in constitutive heterochromatin, resulting in significant effects on the definition and distribution of chromatin states in and outside of constitutive heterochromatin. We thus propose that dynamic exchanges of histone variants control the organization of histone modifications into chromatin states, acting as molecular landmarks.

1. Core samples from trees are a critical reservoir of ecological information, informing our understanding of past climates, as well as contemporary ecosystem responses to global change. Manual measurements of annual growth rings in trees are slow, labour- intensive and subject to human bias, hindering the generation of big datasets. We present an alternative, neural network- based implementation that automates detection and measurement of tree- ring boundaries from coniferous species.2. We trained our Mask R-CNN extensively on over 8000 manually annotated ring boundaries from microscope-imaged Norway Spruce Picea abies increment cores. We assessed the performance of the trained model after post- processing on real- world data generated from our core processing pipeline.3. The CNN after post- processing performed well, with recognition of over 98% of ring boundaries (recall) with a precision in detection of 96% when tested on real- world data. Additionally, we have implemented automatic measurements based on minimum distance between rings. With minimal editing for missed ring de-tections, these measurements were 98% correlated with human measurements of the same samples. Tests on other three conifer species demonstrate that the CNN generalizes well to other species with similar structure.4. We demonstrate the efficacy of automating the measurement of growth increment in tree core samples. Our CNN- based system provides high predictive performance in terms of both tree- ring detection and growth rate determination. Our applica-tion is readily deployable as a Docker container and requires only basic command line skills. Additionally, an easy re- training option allows users to expand capabili-ties to other wood types. Application outputs include both editable annotations of predictions as well as ring- width measurements in a commonly used .pos format, facilitating the efficient generation of large ring- width measurement datasets from increment core samples, an important source of environmental data.

The ability to map a proteomic fingerprint to transcriptomic data would master the understanding of how gene expression translates into actual phenotype. In contrast to nucleic acid sequencing, in vitro protein amplification is impossible and no single cell proteomic workflow has been established as gold standard yet. Advances in microfluidic sample preparation, multi-dimensional sample separation, sophisticated data acquisition strategies, and intelligent data analysis algorithms have resulted in major improvements to successfully analyze such tiny sample amounts with steadily boosted performance. However, among the broad variation of published approaches, it is commonly accepted that highest possible sensitivity, robustness, and throughput are still the most urgent needs for the field. While many labs have focused on multiplexing to achieve these goals, label-free SCP is a highly promising strategy as well whenever high dynamic range and unbiased accurate quantification are needed. We here focus on recent advances in label-free single-cell mass spectrometry workflows and try to guide our readers to choose the best method or combinations of methods for their specific applications. We further highlight which techniques are most propitious in the future and which applications but also limitations we foresee for the field.

Multicellular organisms perceive and transduce multiple cues to optimize development. Key transcription factors drive developmental changes, but RNA processing also contributes to tissue development. Here, we report that multiple decapping deficient mutants share developmental defects in apical hook, primary and lateral root growth. More specifically, LATERAL ORGAN BOUNDARIES DOMAIN 3 (LBD3)/ASYMMETRIC LEAVES 2-LIKE 9 (ASL9) transcripts accumulate in decapping deficient plants and can be found in complexes with decapping components. Accumulation of ASL9 inhibits apical hook and lateral root formation. Interestingly, exogenous auxin application restores lateral roots formation in both ASL9 over-expressors and mRNA decay-deficient mutants. Likewise, mutations in the cytokinin transcription factors type-B ARABIDOPSIS RESPONSE REGULATORS (B-ARRs) ARR10 and ARR12 restore the developmental defects caused by over-accumulation of capped ASL9 transcript upon ASL9 overexpression. Most importantly, loss-of-function of asl9 partially restores apical hook and lateral root formation in both dcp5-1 and pat triple decapping deficient mutants. Thus, the mRNA decay machinery directly targets ASL9 transcripts for decay, possibly to interfere with cytokinin/auxin responses, during development.

Although the biological utilities of endogenous RNAi (endo-RNAi) have been largely elusive, recent studies reveal its critical role in the non-model fruitfly Drosophila simulans to suppress selfish genes, whose unchecked activities can severely impair spermatogenesis. In particular, hairpin RNA (hpRNA) loci generate endo-siRNAs that suppress evolutionary novel, X-linked, meiotic drive loci. The consequences of deleting even a single hpRNA (Nmy) in males are profound, as such individuals are nearly incapable of siring male progeny. Here, comparative genomic analyses of D. simulans and D. melanogaster mutants of the core RNAi factor dcr-2 reveal a substantially expanded network of recently-emerged hpRNA-target interactions in the former species. The de novo hpRNA regulatory network in D. simulans provides insight into molecular strategies that underlie hpRNA emergence and their potential roles in sex chromosome conflict. In particular, our data support the existence of ongoing rapid evolution of Nmy/Dox-related networks, and recurrent targeting of testis HMG-box loci by hpRNAs. Importantly, the impact of the endo-RNAi network on gene expression flips the convention for regulatory networks, since we observe strong derepression of targets of the youngest hpRNAs, but only mild effects on the targets of the oldest hpRNAs. These data suggest that endo-RNAi are especially critical during incipient stages of intrinsic sex chromosome conflicts, and that continual cycles of distortion and resolution may contribute to speciation.

Centromeres are critical for cell division, loading CENH3 or CENPA histone variant nucleosomes, directing kinetochore formation and allowing chromosome segregation1,2. Despite their conserved function, centromere size and structure are diverse across species. To understand this centromere paradox3,4, it is necessary to know how centromeric diversity is generated and whether it reflects ancient trans-species variation or, instead, rapid post-speciation divergence. To address these questions, we assembled 346 centromeres from 66 Arabidopsis thaliana and 2 Arabidopsis lyrata accessions, which exhibited a remarkable degree of intra- and inter-species diversity. A. thaliana centromere repeat arrays are embedded in linkage blocks, despite ongoing internal satellite turnover, consistent with roles for unidirectional gene conversion or unequal crossover between sister chromatids in sequence diversification. Additionally, centrophilic ATHILA transposons have recently invaded the satellite arrays. To counter ATHILA invasion, chromosome-specific bursts of satellite homogenization generate higher-order repeats and purge transposons, in line with cycles of repeat evolution. Centromeric sequence changes are even more extreme in comparison between A. thaliana and A. lyrata. Together, our findings identify rapid cycles of transposon invasion and purging through satellite homogenization, which drive centromere evolution and ultimately contribute to speciation.

Isothermal titration calorimetry (ITC) is the gold standard for providing quantitative and thermodynamic understanding of the interaction mechanisms between core autophagy machinery, autophagy receptors, and ATG8. Here, we used two model peptides and Arabidopsis thaliana ATG8A to characterize ATG8-peptide interactions. We employed ITC using three different methods (direct ligand titration, displacement, and competition assays) to characterize, directly and indirectly, the interaction of the peptides with ATG8. We then analyzed the ITC data by global and statistical methods and discussed advantages, drawbacks, and negative controls for each approach. We finally provide a thorough description of all the steps, including data analysis and presentation, preparation of recombinant ATG8A from E. coli, and troubleshooting notes for technical problems that can be encountered. Although we used ATG8-peptide interactions here, these assays can be applied to any other one-to-one protein-protein and ligand-protein interactions and competitive binders.

Genome integrity in animals depends on silencing of mobile elements by Piwi-interacting RNAs (piRNAs). A new study in this issue of PLOS Biology reveals recent evolutionary losses of key piRNA biogenesis factors in flies, highlighting adaptability by rapid shift to alternative piRNA biogenesis strategies.

The view on the role of light during seed germination stems mainly from studies with Arabidopsis (Arabidopsis thaliana), where light is required to initiate this process. In contrast, white light is a strong inhibitor of germination in other plants, exemplified by accessions of Aethionema arabicum, another member of Brassicaceae. Their seeds respond to light with gene expression changes of key regulators converse to that of Arabidopsis, resulting in opposite hormone regulation and prevention of germination. However, the photoreceptors involved in this process in A. arabicum remain unknown. Here, we screened a mutant collection of A. arabicum and identified koy-1, a mutant that lost light inhibition of germination due to a deletion in the promoter of HEME OXYGENASE 1, the gene for a key enzyme in the biosynthesis of the phytochrome chromophore. koy-1 seeds were unresponsive to red- and far-red light and hyposensitive under white light. Comparison of hormone and gene expression between wild type and koy-1 revealed that very low light fluence stimulates germination, while high irradiance of red and far-red light is inhibitory, indicating a dual role of phytochromes in light-regulated seed germination. The mutation also affects the ratio between the 2 fruit morphs of A. arabicum, suggesting that light reception via phytochromes can fine-tune several parameters of propagation in adaptation to conditions in the habitat.

UFMylation involves the covalent modification of substrate proteins with UFM1 (Ubiquitin-fold modifier 1) and is important for maintaining ER homeostasis. Stalled translation triggers the UFMylation of ER-bound ribosomes and activates C53-mediated autophagy to clear toxic polypeptides. C53 contains noncanonical shuffled ATG8-interacting motifs (sAIMs) that are essential for ATG8 interaction and autophagy initiation. However, the mechanistic basis of sAIM-mediated ATG8 interaction remains unknown. Here, we show that C53 and sAIMs are conserved across eukaryotes but secondarily lost in fungi and various algal lineages. Biochemical assays showed that the unicellular alga Chlamydomonas reinhardtii has a functional UFMylation pathway, refuting the assumption that UFMylation is linked to multicellularity. Comparative structural analyses revealed that both UFM1 and ATG8 bind sAIMs in C53, but in a distinct way. Conversion of sAIMs into canonical AIMs impaired binding of C53 to UFM1, while strengthening ATG8 binding. Increased ATG8 binding led to the autoactivation of the C53 pathway and sensitization of Arabidopsis thaliana to ER stress. Altogether, our findings reveal an ancestral role of sAIMs in UFMylation-dependent fine-tuning of C53-mediated autophagy activation.

Phenotypic plasticity is the ability of a given genotype to produce multiple phenotypes in response to changing environmental conditions. Understanding the genetic basis of phenotypic plasticity and establishing a predictive model is highly relevant to future agriculture under a changing climate. Here we report findings on the genetic basis of phenotypic plasticity for 23 complex traits using a diverse maize population planted at five sites with distinct environmental conditions. We found that latitude-related environmental factors were the main drivers of across-site variation in flowering time traits but not in plant architecture or yield traits. For the 23 traits, we detected 109 quantitative trait loci (QTLs), 29 for mean values, 66 for plasticity, and 14 for both parameters, and 80% of the QTLs interacted with latitude. The effects of several QTLs changed in magnitude or sign, driving variation in phenotypic plasticity. We experimentally validated one plastic gene, ZmTPS14.1, whose effect was likely mediated by the compensation effect of ZmSPL6 from a downstream pathway. By integrating genetic diversity, environmental variation, and their interaction into a joint model, we could provide site-specific predictions with increased accuracy by as much as 9.9%, 2.2%, and 2.6% for days to tassel, plant height, and ear weight, respectively. This study revealed a complex genetic architecture involving multiple alleles, pleiotropy, and genotype-by-environment interaction that underlies variation in the mean and plasticity of maize complex traits. It provides novel insights into the dynamic genetic architecture of agronomic traits in response to changing environments, paving a practical way toward precision agriculture.

Gene-body methylation (gbM) refers to sparse CG methylation of coding regions, which is especially prominent in evolutionarily conserved house-keeping genes. It is found in both plants and animals, but is directly and stably (epigenetically) inherited over multiple generations in the former. Studies in Arabidopsis thaliana have demonstrated that plants originating from different parts of the world exhibit genome-wide differences in gbM, which could reflect direct selection on gbM, but which could also reflect an epigenetic memory of ancestral genetic and/or environmental factors. Here we look for evidence of such factors in F2 plants resulting from a cross between a southern Swedish line with low gbM and a northern Swedish line with high gbM, grown at two different temperatures. Using bisulfite-sequencing data with nucleotide-level resolution on hundreds of individuals, we confirm that CG sites are either methylated (nearly 100% methylation across sampled cells) or unmethylated (approximately 0% methylation across sampled cells), and show that the higher level of gbM in the northern line is due to more sites being methylated. Furthermore, methylation variants almost always show Mendelian segregation, consistent with their being directly and stably inherited through meiosis. To explore how the differences between the parental lines could have arisen, we focused on somatic deviations from the inherited state, distinguishing between gains (relative to the inherited 0% methylation) and losses (relative to the inherited 100% methylation) at each site in the F2 generation. We demonstrate that deviations predominantly affect sites that differ between the parental lines, consistent with these sites being more mutable. Gains and losses behave very differently in terms of the genomic distribution, and are influenced by the local chromatin state. We find clear evidence for different trans-acting genetic polymorphism affecting gains and losses, with those affecting gains showing strong environmental interactions (G×E). Direct effects of the environment were minimal. In conclusion, we show that genetic and environmental factors can change gbM at a cellular level, and hypothesize that these factors can also lead to transgenerational differences between individuals via the inclusion of such changes in the zygote. If true, this could explain genographic pattern of gbM with selection, and would cast doubt on estimates of epimutation rates from inbred lines in constant environments.

Poly(ADP-ribosyl)ation (PARylation) is a reversible post-translational protein modification that has profound regulatory functions in metabolism, development and immunity, and is conserved throughout the eukaryotic lineage. Contrary to metazoa, many components and mechanistic details of PARylation have remained unidentified in plants. Here we present the transcriptional co-regulator RADICAL-INDUCED CELL DEATH1 (RCD1) as a plant PAR-reader. RCD1 is a multidomain protein with intrinsically disordered regions (IDRs) separating its domains. We have reported earlier that RCD1 regulates plant development and stress-tolerance by interacting with numerous transcription factors (TFs) through its C-terminal RST domain. This study suggests that the N-terminal WWE and PARP-like domains, as well as the connecting IDR play an important regulatory role for RCD1 function. We show that RCD1 binds PAR in vitro via its WWE domain and that PAR-binding determines RCD1 localization to nuclear bodies (NBs) in vivo. Additionally, we found that RCD1 function and stability is controlled by Photoregulatory Protein Kinases (PPKs). PPKs localize with RCD1 in NBs and phosphorylate RCD1 at multiple sites affecting its stability. This work proposes a mechanism for negative transcriptional regulation in plants, in which RCD1 localizes to NBs, binds TFs with its RST domain and is degraded after phosphorylation by PPKs.

Repeated herbicide applications in agricultural fields exert strong selection on weeds such as blackgrass (Alopecurus myosuroides), which is a major threat for temperate climate cereal crops. This inadvertent selection pressure provides an opportunity for investigating the underlying genetic mechanisms and evolutionary processes of rapid adaptation, which can occur both through mutations in the direct targets of herbicides and through changes in other, often metabolic, pathways, known as non-target-site resistance. How much target-site resistance (TSR) relies on de novo mutations vs. standing variation is important for developing strategies to manage herbicide resistance. We first generated a chromosome-level reference genome for A. myosuroides for population genomic studies of herbicide resistance and genome-wide diversity across Europe in this species. Next, through empirical data in the form of highly accurate long-read amplicons of alleles encoding acetyl-CoA carboxylase (ACCase) and acetolactate synthase (ALS) variants, we showed that most populations with resistance due to TSR mutations-23 out of 27 and six out of nine populations for ACCase and ALS, respectively-contained at least two TSR haplotypes, indicating that soft sweeps are the norm. Finally, through forward-in-time simulations, we inferred that TSR is likely to mainly result from standing genetic variation, with only a minor role for de novo mutations.

Archaeological cobs from Paredones and Huaca Prieta (Peru) represent some of the oldest maize known to date, yet they present relevant phenotypic traits corresponding to domesticated maize. This contrasts with the earliest Mexican macro-specimens from Guila Naquitz and San Marcos, which are phenotypically intermediate for these traits, even though they date more recently in time. To gain insights into the origins of ancient Peruvian maize, we sequenced DNA from three Paredones specimens dating ~6700-5000 calibrated years before present (BP), conducting comparative analyses with two teosinte subspecies (Zea mays ssp. mexicana and parviglumis) and extant maize, that include highland and lowland landraces from Mesoamerica and South America. We show that Paredones maize originated from the same domestication event as Mexican maize and was domesticated by ~6700 BP, implying rapid dispersal followed by improvement. Paredones maize shows no relevant gene flow from mexicana, smaller than that observed in teosinte parviglumis. Thus, Paredones samples represent the only maize without confounding mexicana variation found to date. It also harbors significantly fewer alleles previously found to be adaptive to highlands, but not of alleles adaptive to lowlands, supporting a lowland migration route. Our overall results imply that Paredones maize originated in Mesoamerica, arrived in Peru without mexicana introgression through a rapid lowland migration route, and underwent improvements in both Mesoamerica and South America.

A sensing mechanism in mammals perceives xenobiotics and induces the transcription of genes encoding proteins that detoxify these molecules. However, it is unclear if plants sense xenobiotics, and activate an analogous signalling system leading to their detoxification. Using the liverwort Marchantia polymorpha, we tested the hypothesis that there is a sensing system in plants that perceives herbicides resulting in the increased transcription of genes encoding proteins that detoxify these herbicides. Consistent with the hypothesis, we show that chlorsulfuron-treatment induces changes in the M. polymorpha transcriptome. However, these transcriptome changes do not occur in chlorsulfuron (CS)-treated target site resistant mutants, where the gene encoding the target carries a mutation that confers resistance to chlorsulfuron. Instead, we show that inactivation of the chlorsulfuron target, acetolactate synthase (ALS) (also known as acetohydroxyacid synthase (AHAS)), is required for the transcriptome response. These data demonstrate that the transcriptome changes in chlorsulfuron-treated plants are caused by disrupted amino acid synthesis and metabolism resulting from acetolactate synthase inhibition, and indicate that the transcriptome changes are not caused by a herbicide sensing mechanism.

Elongation factors modulate the efficiency of mRNA synthesis by RNA polymerase II (RNAPII) in the context of chromatin, thus contributing to implement proper gene expression programmes. The zinc-finger protein elongation factor 1 (ELF1) is a conserved transcript elongation factor (TEF), whose molecular function so far has not been studied in plants. Using biochemical approaches, we examined the interaction of Arabidopsis ELF1 with DNA and histones in vitro and with the RNAPII elongation complex in vivo. In addition, cytological assays demonstrated the nuclear localisation of the protein, and by means of double-mutant analyses, interplay with genes encoding other elongation factors was explored. The genome-wide distribution of ELF1 was addressed by chromatin immunoprecipitation. ELF1 isolated from Arabidopsis cells robustly copurified with RNAPII and various other elongation factors including SPT4-SPT5, SPT6, IWS1, FACT and PAF1C. Analysis of a CRISPR-Cas9-mediated gene editing mutant of ELF1 revealed distinct genetic interactions with mutants deficient in other elongation factors. Moreover, ELF1 associated with genomic regions actively transcribed by RNAPII. However, ELF1 occupied only c. 33% of the RNAPII transcribed loci with preference for inducible rather than constitutively expressed genes. Collectively, these results establish that Arabidopsis ELF1 shares several characteristic attributes with RNAPII TEFs.

During meiotic prophase I, sister chromatids are arranged in a loop-base array along a proteinaceous structure, called the meiotic chromosome axis. This structure is essential for synapsis and meiotic recombination progression and hence formation of genetically diverse gametes. Proteomic studies in plants aiming to unravel the composition and regulation of meiotic axes are constrained by limited meiotic cells embedded in floral organs. Here we report TurboID (TbID)-based proximity labelling (PL) in meiotic cells of Arabidopsis thaliana. TbID fusion to the two meiotic chromosome axis proteins ASY1 and ASY3 enabled the identification of their proximate 'interactomes' based on affinity purification coupled with mass spectrometry. We identified 39 ASY1 and/or ASY3 proximate candidates covering most known chromosome axis-related proteins. Functional studies of selected candidates demonstrate that not only known meiotic candidates but also new meiotic proteins were uncovered. Hence, TbID-based PL in meiotic cells enables the identification of chromosome axis proximate proteins in A. thaliana.

In contrast to other eukaryotic model organisms, the closely related ubiquitin (Ub)-conjugating enzymes UBC35 and UBC36 are the main sources of K63-linked Ub chains in Arabidopsis.1 Although K63-linked chains have been associated with the regulation of vesicle trafficking, definitive proof for their role in endocytosis was missing. We show that the ubc35 ubc36 mutant has pleiotropic phenotypes related to hormone and immune signaling. Specifically, we reveal that ubc35-1 ubc36-1 plants have altered turnover of integral membrane proteins including FLS2, BRI1, and PIN1 at the plasma membrane. Our data indicates that K63-Ub chains are generally required for endocytic trafficking in plants. In addition, we show that in plants K63-Ub chains are involved in selective autophagy through NBR1, the second major pathway delivering cargoes to the vacuole for degradation. Similar to autophagy-defective mutants, ubc35-1 ubc36-1 plants display an accumulation of autophagy markers. Moreover, autophagy receptor NBR1 interacts with K63-Ub chains, which are required for its delivery to the lytic vacuole.2 Together, we show that K63-Ub chains act as a general signal required for the two main pathways delivering cargo to the vacuole and thus, to maintain proteostasis.

The energy sensor AMP-activated protein kinase (AMPK) can activate autophagy when cellular energy production becomes compromised. However, the degree to which nutrient sensing impinges on the autophagosome closure remains unknown. Here, we provide the mechanism underlying a plant unique protein FREE1, upon autophagy-induced SnRK1α1-mediated phosphorylation, functions as a linkage between ATG conjugation system and ESCRT machinery to regulate the autophagosome closure upon nutrient deprivation. Using high-resolution microscopy, 3D-electron tomography, and protease protection assay, we showed that unclosed autophagosomes accumulated in free1 mutants. Proteomic, cellular and biochemical analysis revealed the mechanistic connection between FREE1 and the ATG conjugation system/ESCRT-III complex in regulating autophagosome closure. Mass spectrometry analysis showed that the evolutionary conserved plant energy sensor SnRK1α1 phosphorylates FREE1 and recruits it to the autophagosomes to promote closure. Mutagenesis of the phosphorylation site on FREE1 caused the autophagosome closure failure. Our findings unveil how cellular energy sensing pathways regulate autophagosome closure to maintain cellular homeostasis.

During development of the human cerebral cortex, multipotent neural progenitors generate excitatory neurons and glial cells. Investigations of the transcriptome and epigenome have revealed important gene regulatory networks underlying this crucial developmental event. However, the posttranscriptional control of gene expression and protein abundance during human corticogenesis remains poorly understood. We addressed this issue by using human telencephalic brain organoids grown using a dual reporter cell line to isolate neural progenitors and neurons and performed cell class and developmental stage-specific transcriptome and proteome analysis. Integrating the two datasets revealed modules of gene expression during human corticogenesis. Investigation of one such module uncovered mTOR-mediated regulation of translation of the 5'TOP element-enriched translation machinery in early progenitor cells. We show that in early progenitors partial inhibition of the translation of ribosomal genes prevents precocious translation of differentiation markers. Overall, our multiomics approach proposes novel posttranscriptional regulatory mechanisms crucial for the fidelity of cortical development.

The Atg8 family of ubiquitin-like proteins play pivotal roles in autophagy and other processes involving vesicle fusion and transport where the lysosome/vacuole is the end station. Nuclear roles of Atg8 proteins are also emerging. Here, we review the structural and functional features of Atg8 family proteins and their protein-protein interaction modes in model organisms such as yeast, Arabidopsis, C. elegans and Drosophila to humans. Although varying in number of homologs, from one in yeast to seven in humans, and more than ten in some plants, there is a strong evolutionary conservation of structural features and interaction modes. The most prominent interaction mode is between the LC3 interacting region (LIR), also called Atg8 interacting motif (AIM), binding to the LIR docking site (LDS) in Atg8 homologs. There are variants of these motifs like "half-LIRs" and helical LIRs. We discuss details of the binding modes and how selectivity is achieved as well as the role of multivalent LIR-LDS interactions in selective autophagy. A number of LIR-LDS interactions are known to be regulated by phosphorylation. New methods to predict LIR motifs in proteins have emerged that will aid in discovery and analyses. There are also other interaction surfaces than the LDS becoming known where we presently lack detailed structural information, like the N-terminal arm region and the UIM-docking site (UDS). More interaction modes are likely to be discovered in future studies.

BACKGROUND: Coccidiosis caused by Eimeria is one of the most severe chicken diseases and poses a great economic threat to the poultry industry. Understanding the evolutionary biology of chicken Eimeria parasites underpins development of new interactions toward the improved prevention and control of this poultry disease. METHODS: We presented an evolutionary blueprint of chicken coccidia by genetically characterizing complete mitogenome assemblies of 33 isolates representing all seven known Eimeria species infecting chickens in China. Further genome- and gene-level phylogenies were also achieved to better understand the evolutionary relationships of these chicken Eimeria at the species level. RESULTS: 33 mitogenomes of chicken eimerian parasites ranged from 6148 bp to 6480 bp in size and encoded typical mitochondrial compositions of apicomplexan parasites including three protein-coding genes (PCGs), seven fragmented small subunit (SSU) and 12/13 fragmented large subunit (LSU) rRNAs. Comparative genomics provided an evolutionary scenario for the genetic diversity of PCGs-cytochrome c oxidase subunits 1 and 3 (cox1 and cox3) and cytochrome b (cytb); all were under purifying selection with cox1 and cox3 being the lowest and highest evolutionary rates, respectively. Genome-wide phylogenies classified the 33 Eimeria isolates into seven subgroups, and furthermore Eimeria tenella and Eimeria necatrix were determined to be more closely related to each other than to the other eight congenic species. Single/concatenated mitochondrial protein gene-based phylogenies supported cox1 as the genetic marker for evolutionary and phylogenetic studies for avain coccidia. CONCLUSIONS: To our knowledge, these are the first population-level mitogenomic data on the genus Eimeria, and its comprehensive molecular characterization provides valuable resources for systematic, population genetic and evolutionary biological studies of apicomplexan parasites in poultry.

Precise organization of growing structures is a fundamental process in developmental biology. In plants, radial growth is mediated by the cambium, a stem cell niche continuously producing wood (xylem) and bast (phloem) in a strictly bidirectional manner. While this process contributes large parts to terrestrial biomass, cambium dynamics eludes direct experimental access due to obstacles in live-cell imaging. Here, we present a cell-based computational model visualizing cambium activity and integrating the function of central cambium regulators. Performing iterative comparisons of plant and model anatomies, we conclude that the receptor-like kinase PXY and its ligand CLE41 are part of a minimal framework sufficient for instructing tissue organization. By integrating tissue-specific cell wall stiffness values, we moreover probe the influence of physical constraints on tissue geometry. Our model highlights the role of intercellular communication within the cambium and shows that a limited number of factors are sufficient to create radial growth by bidirectional tissue production.

BACKGROUND: It is apparent that genomes harbor much structural variation that is largely undetected for technical reasons. Such variation can cause artifacts when short-read sequencing data are mapped to a reference genome. Spurious SNPs may result from mapping of reads to unrecognized duplicated regions. Calling SNP using the raw reads of the 1001 Arabidopsis Genomes Project we identified 3.3 million (44%) heterozygous SNPs. Given that Arabidopsis thaliana (A. thaliana) is highly selfing, and that extensively heterozygous individuals have been removed, we hypothesize that these SNPs reflected cryptic copy number variation. RESULTS: The heterozygosity we observe consists of particular SNPs being heterozygous across individuals in a manner that strongly suggests it reflects shared segregating duplications rather than random tracts of residual heterozygosity due to occasional outcrossing. Focusing on such pseudo-heterozygosity in annotated genes, we use genome-wide association to map the position of the duplicates. We identify 2500 putatively duplicated genes and validate them using de novo genome assemblies from six lines. Specific examples included an annotated gene and nearby transposon that transpose together. We also demonstrate that cryptic structural variation produces highly inaccurate estimates of DNA methylation polymorphism. CONCLUSIONS: Our study confirms that most heterozygous SNP calls in A. thaliana are artifacts and suggest that great caution is needed when analyzing SNP data from short-read sequencing. The finding that 10% of annotated genes exhibit copy-number variation, and the realization that neither gene- nor transposon-annotation necessarily tells us what is actually mobile in the genome suggests that future analyses based on independently assembled genomes will be very informative.

The analysis of ultralow input samples or even individual cells is essential to answering a multitude of biomedical questions, but current proteomic workflows are limited in their sensitivity and reproducibility. Here, we report a comprehensive workflow that includes improved strategies for all steps, from cell lysis to data analysis. Thanks to convenient-to-handle 1 μL sample volume and standardized 384-well plates, the workflow is easy for even novice users to implement. At the same time, it can be performed semi-automatized using CellenONE, which allows for the highest reproducibility. To achieve high throughput, ultrashort gradient lengths down to 5 min were tested using advanced μ-pillar columns. Data-dependent acquisition (DDA), wide-window acquisition (WWA), data-independent acquisition (DIA), and commonly used advanced data analysis algorithms were benchmarked. Using DDA, 1790 proteins covering a dynamic range of four orders of magnitude were identified in a single cell. Using DIA, proteome coverage increased to more than 2200 proteins identified from single-cell level input in a 20 min active gradient. The workflow enabled differentiation of two cell lines, demonstrating its suitability to cellular heterogeneity determination.

Meiosis is a specialized cell division that halves the number of chromosomes in two consecutive rounds of chromosome segregation. In angiosperm plants is meiosis followed by mitotic divisions to form rudimentary haploid gametophytes. In Arabidopsis, termination of meiosis and transition to gametophytic development are governed by TDM1 and SMG7 that mediate inhibition of translation. Mutants deficient in this mechanism do not form tetrads but instead undergo multiple cycles of aberrant nuclear divisions that are likely caused by the failure to downregulate cyclin dependent kinases during meiotic exit. A suppressor screen to identify genes that contribute to meiotic exit uncovered a mutation in cyclin-dependent kinase D;3 (CDKD;3) that alleviates meiotic defects in smg7 deficient plants. The CDKD;3 deficiency prevents aberrant meiotic divisions observed in smg7 mutants or delays their onset after initiation of cytokinesis, which permits formation of functional microspores. Although CDKD;3 acts as an activator of cyclin-dependent kinase A;1 (CDKA;1), the main cyclin dependent kinase that regulates meiosis, cdkd;3 mutation appears to promote meiotic exit independently of CDKA;1. Furthermore, analysis of CDKD;3 interactome revealed enrichment for proteins implicated in cytokinesis, suggesting a more complex function of CDKD;3 in cell cycle regulation.

The shape of modular organisms depends on the branching architecture, which in plants is determined by the fates of generative centers called meristems. The branches of the liverwort Marchantia polymorpha are derived from two adjacent meristems that develop at thallus apices. These meristems may be active and develop branches or may be dormant and do not form branches. The relative number and position of active and dormant meristems define the overall shape and form of the thallus. We show that the clade III SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factor MpSPL1 is required for meristem dormancy. The activity of MpSPL1 is regulated by the liverwort-specific Mpo-MR13 miRNA, which, in turn, is regulated by PIF-mediated signaling. An unrelated PIF-regulated miRNA, MIR156, represses a different SPL gene (belonging to clade IV) that inhibits branching during the shade avoidance response in Arabidopsis thaliana. This suggests that a conserved light signaling mechanism modulates branching architecture in liverworts and angiosperms and therefore is likely operated in the last common ancestor. However, PIF-mediated signaling represses the expression of different miRNA genes with different SPL targets during dichotomous, apical branching in liverworts and during lateral, subapical branching in angiosperms. We speculate that the mechanism that acts downstream of light and regulates meristem dormancy evolved independently in liverworts and angiosperms.

Cytochrome P450 (CYP) monooxygenases and glutathione S-transferases (GST) are enzymes that catalyse chemical modifications of a range of organic compounds. Herbicide resistance has been associated with higher levels of CYP and GST gene expression in some herbicide-resistant weed populations compared to sensitive populations of the same species. By comparing the protein sequences of 9 representative species of the Archaeplastida-the lineage which includes red algae, glaucophyte algae, chlorophyte algae, and streptophytes-and generating phylogenetic trees, we identified the CYP and GST proteins that existed in the common ancestor of the Archaeplastida. All CYP clans and all but one land plant GST classes present in land plants evolved before the divergence of streptophyte algae and land plants from their last common ancestor. We also demonstrate that there are more genes encoding CYP and GST proteins in land plants than in algae. The larger numbers of genes among land plants largely results from gene duplications in CYP clans 71, 72, and 85 and in the GST phi and tau classes [1,2]. Enzymes that either metabolise herbicides or confer herbicide resistance belong to CYP clans 71 and 72 and the GST phi and tau classes. Most CYP proteins that have been shown to confer herbicide resistance are members of the CYP81 family from clan 71. These results demonstrate that the clan and class diversity in extant plant CYP and GST proteins had evolved before the divergence of land plants and streptophyte algae from a last common ancestor estimated to be between 515 and 474 million years ago. Then, early in embryophyte evolution during the Palaeozoic, gene duplication in four of the twelve CYP clans, and in two of the fourteen GST classes, led to the large numbers of CYP and GST proteins found in extant land plants. It is among the genes of CYP clans 71 and 72 and GST classes phi and tau that alleles conferring herbicide resistance evolved in the last fifty years.

After double fertilization, the endosperm in the seeds of many flowering plants undergoes repeated mitotic nuclear divisions without cytokinesis, resulting in a large coenocytic endosperm that then cellularizes. Growth during the coenocytic phase is strongly associated with the final seed size; however, a detailed description of the cellular dynamics controlling the unique coenocytic development in flowering plants has remained elusive. By integrating confocal microscopy live-cell imaging and genetics, we have characterized the entire development of the coenocytic endosperm of Arabidopsis thaliana including nuclear divisions, their timing intervals, nuclear movement and cytoskeleton dynamics. Around each nucleus, microtubules organize into aster-shaped structures that drive actin filament (F-actin) organization. Microtubules promote nuclear movement after division, while F-actin restricts it. F-actin is also involved in controlling the size of both the coenocytic endosperm and the mature seed. The characterization of cytoskeleton dynamics in real time throughout the entire coenocyte endosperm period provides foundational knowledge of plant coenocytic development, insights into the coordination of F-actin and microtubules in nuclear dynamics, and new opportunities to increase seed size and our food security.

The phenotypes of plants can be influenced by the environmental conditions experienced by their parents. However, there is still much uncertainty about how common and how predictable such parental environmental effects really are. We carried out a comprehensive experimental test for parental effects, subjecting plants of multiple Arabidopsis thaliana genotypes to 24 different biotic or abiotic stresses, or combinations thereof, and comparing their offspring phenotypes in a common environment. The majority of environmental stresses caused significant parental effects, with -35% to +38% changes in offspring fitness. The expression of parental effects was strongly genotype-dependent, and multiple environmental stresses often acted nonadditively when combined. The direction and magnitude of parental effects were unrelated to the direct effects on the parents: Some environmental stresses did not affect the parents but caused substantial effects on offspring, while for others, the situation was reversed. Our study demonstrates that parental environmental effects are common and often strong in A. thaliana, but they are genotype-dependent, act nonadditively, and are difficult to predict. We should thus be cautious with generalizing from simple studies with single plant genotypes and/or only few individual environmental stresses. A thorough and general understanding of parental effects requires large multifactorial experiments.

Spectral library search can enable more sensitive peptide identification in tandem mass spectrometry experiments. However, its drawbacks are the limited availability of high-quality libraries and the added difficulty of creating decoy spectra for result validation. We describe MS Ana, a new spectral library search engine that enables high sensitivity peptide identification using either curated or predicted spectral libraries as well as robust false discovery control through its own decoy library generation algorithm. MS Ana identifies on average 36% more spectrum matches and 4% more proteins than database search in a benchmark test on single-shot human cell-line data. Further, we demonstrate the quality of the result validation with tests on synthetic peptide pools and show the importance of library selection through a comparison of library search performance with different configurations of publicly available human spectral libraries.

In plants, the endomembrane system is tightly regulated in response to environmental stresses for maintaining cellular homeostasis. Autophagosomes, the double membrane organelles forming upon nutrient deprivation or stress induction, degrade bulky cytosolic materials for nutrient turnover. Though abiotic stresses have been reported to induce plant autophagy, few receptors or regulators for selective autophagy have been characterized for specific stresses. Here, we have applied immunoprecipitation followed by tandem mass spectrometry using the autophagosome marker protein ATG8 as bait and have identified the E3 ligase of the ufmylation system Ufl1 as a bona fide ATG8 interactor under salt stress. Notably, core components in the ufmylation cascade, Ufl1 and Ufm1, interact with the autophagy kinase complexes proteins ATG1 and ATG6. Cellular and genetic analysis showed that Ufl1 is important for endoplasmic reticulum (ER)-phagy under persisting salt stress. Loss-of-function mutants of Ufl1 display a salt stress hypersensitive phenotype and abnormal ER morphology. Prolonged ER stress responses are detected in ufl1 mutants that phenocopy the autophagy dysfunction atg5 mutants. Consistently, expression of ufmylation cascade components is up-regulated by salt stress. Taken together, our study demonstrates the role of ufmylation in regulating ER homeostasis under salt stress through ER-phagy.

UNLABELLED: The skin is exposed to viral pathogens, but whether they contribute to the oncogenesis of skin cancers has not been systematically explored. Here we investigated 19 skin tumor types by analyzing off-target reads from commonly available next-generation sequencing data for viral pathogens. We identified human papillomavirus 42 (HPV42) in 96% (n = 45/47) of digital papillary adenocarcinoma (DPA), an aggressive cancer occurring on the fingers and toes. We show that HPV42, so far considered a nononcogenic, "low-risk" HPV, recapitulates the molecular hallmarks of oncogenic, "high-risk" HPVs. Using machine learning, we find that HPV-driven transformation elicits a germ cell-like transcriptional program conserved throughout all HPV-driven cancers (DPA, cervical carcinoma, and head and neck cancer). We further show that this germ cell-like transcriptional program, even when reduced to the top two genes (CDKN2A and SYCP2), serves as a fingerprint of oncogenic HPVs with implications for early detection, diagnosis, and therapy of all HPV-driven cancers. SIGNIFICANCE: We identify HPV42 as a uniform driver of DPA and add a new member to the short list of tumorigenic viruses in humans. We discover that all oncogenic HPVs evoke a germ cell-like transcriptional program with important implications for detecting, diagnosing, and treating all HPV-driven cancers. See related commentary by Starrett et al., p. 17. This article is highlighted in the In This Issue feature, p. 1.

Herbicide resistance in weeds is a growing threat to global crop production. Non-target site resistance is problematic because a single resistance allele can confer tolerance to many herbicides (cross resistance), and it is often a polygenic trait so it can be difficult to identify the molecular mechanisms involved. Most characterized molecular mechanisms of non-target site resistance are caused by gain-of-function mutations in genes from a few key gene families-the mechanisms of resistance caused by loss-of-function mutations remain unclear. In this study, we first show that the mechanism of non-target site resistance to the herbicide thaxtomin A conferred by loss-of-function of the gene PAM16 is conserved in Marchantia polymorpha, validating its use as a model species with which to study non-target site resistance. To identify mechanisms of non-target site resistance caused by loss-of-function mutations, we generated 107 UV-B mutagenized M. polymorpha spores and screened for resistance to the herbicide thaxtomin A. We isolated 13 thaxtomin A-resistant mutants and found that 3 mutants carried candidate resistance-conferring SNPs in the MpRTN4IP1L gene. Mprtn4ip1l mutants are defective in coenzyme Q biosynthesis and accumulate higher levels of reactive oxygen species (ROS) than wild-type plants. Mutants are weakly resistant to thaxtomin A and cross resistant to isoxaben, suggesting that loss of MpRTN4IP1L function confers non-target site resistance. Mutants are also defective in thaxtomin A metabolism. We conclude that loss of MpRTN4IP1L function is a novel mechanism of non-target site herbicide resistance and propose that other mutations that increase ROS levels or decrease thaxtomin A metabolism could contribute to thaxtomin A resistance in the field.

Flagellin, the protein subunit of the bacterial flagellum, stimulates the innate immune receptor Toll-like receptor 5 (TLR5) after pattern recognition or evades TLR5 through lack of recognition. This binary response fails to explain the weak agonism of flagellins from commensal bacteria, raising the question of how TLR5 response is tuned. Here, we screened abundant flagellins present in metagenomes from human gut for both TLR5 recognition and activation and uncovered a class of flagellin-TLR5 interaction termed silent recognition. Silent flagellins were weak TLR5 agonists despite pattern recognition. Receptor activity was tuned by a TLR5-flagellin interaction distal to the site of pattern recognition that was present in Salmonella flagellin but absent in silent flagellins. This interaction enabled flagellin binding to preformed TLR5 dimers and increased TLR5 signaling by several orders of magnitude. Silent recognition by TLR5 occurred in human organoids and mice, and silent flagellin proteins were present in human stool. These flagellins were produced primarily by the abundant gut bacteria Lachnospiraceae and were enriched in nonindustrialized populations. Our findings provide a mechanism for the innate immune system to tolerate commensal-derived flagellins while remaining vigilant to the presence of flagellins produced by pathogens.

The acquisition of germination and post-embryonic developmental ability during seed maturation is vital for seed vigor, an important trait for plant propagation and crop production. How seed vigor is established in seeds is still poorly understood. Here, we report the crucial function of Arabidopsis histone variant H3.3 in endowing seeds with post-embryonic developmental potentials. H3.3 is not essential for seed formation, but loss of H3.3 results in severely impaired germination and post-embryonic development. H3.3 exhibits a seed-specific 5' gene end distribution and facilitates chromatin opening at regulatory regions in seeds. During germination, H3.3 is essential for proper gene transcriptional regulation. Moreover, H3.3 is constantly loaded at the 3' gene end, correlating with gene body DNA methylation and the restriction of chromatin accessibility and cryptic transcription at this region. Our results suggest a fundamental role of H3.3 in initiating chromatin accessibility at regulatory regions in seed and licensing the embryonic to post-embryonic transition.

Autophagosomes are double-membraned vesicles that traffic harmful or unwanted cellular macromolecules to the vacuole for recycling. Although autophagosome biogenesis has been extensively studied, autophagosome maturation, i.e., delivery and fusion with the vacuole, remains largely unknown in plants. Here, we have identified an autophagy adaptor, CFS1, that directly interacts with the autophagosome marker ATG8 and localizes on both membranes of the autophagosome. Autophagosomes form normally in Arabidopsis thaliana cfs1 mutants, but their delivery to the vacuole is disrupted. CFS1's function is evolutionarily conserved in plants, as it also localizes to the autophagosomes and plays a role in autophagic flux in the liverwort Marchantia polymorpha. CFS1 regulates autophagic flux by bridging autophagosomes with the multivesicular body-localized ESCRT-I component VPS23A, leading to the formation of amphisomes. Similar to CFS1-ATG8 interaction, disrupting the CFS1-VPS23A interaction blocks autophagic flux and renders plants sensitive to nitrogen starvation. Altogether, our results reveal a conserved vacuolar sorting hub that regulates autophagic flux in plants.

In the field of liquid chromatography-mass spectrometry (LC-MS)-based proteomics, increases in the sampling depth and proteome coverage have mainly been accomplished by rapid advances in mass spectrometer technology. The comprehensiveness and quality of the data that can be generated do, however, also depend on the performance provided by nano-liquid chromatography (nanoLC) separations. Proper selection of reversed-phase separation columns can be important to provide the MS instrument with peptides at the highest possible concentration and separated at the highest possible resolution. In the current contribution, we evaluate the use of the prototype generation 2 μPAC nanoLC columns, which use C18-functionalized superficially porous micropillars as a stationary phase. When compared to traditionally used fully porous silica stationary phases, more precursors could be characterized when performing single shot data-dependent LC-MS/MS analyses of a human cell line tryptic digest. Up to 30% more protein groups and 60% more unique peptides were identified for short gradients (10 min) and limited sample amounts (10-100 ng of cell lysate digest). With LC-MS gradient times of 10, 60, 120, and 180 min, respectively, we identified 2252, 6513, 7382, and 8174 protein groups with 25, 500, 1000, and 2000 ng of the sample loaded on the column. Reduction of sample carryover to the next run (up to 2 to 3%) and decreased levels of methionine oxidation (up to 3-fold) were identified as additional figures of merit. When analyzing a disuccinimidyl dibutyric urea-crosslinked synthetic library, 29 to 59 more unique crosslinked peptides could be identified at an experimentally validated false discovery rate of 1-2%.

Sarcoptes scabiei is the causative agent of a highly contagious skin disease in humans and more than 100 mammals. Here, we reported the first chromosome-level reference genome of S. scabiei isolated from rabbits, with a contig N50 size of 5.92 Mb, a total assembled length of 57.30 Mb, ∼12.65% repetitive sequences and 9333 predicted protein-coding genes. The phylogenetic tree based on 1338 shared high-confidence single-copy orthologous genes estimated that the mammalian ectoparasite S. scabiei and the plant-feeding mite Tetranychus urticae separated ∼340 million years ago. Both neighbour-joining tree and principal component analysis of 20 mite populations isolated from four hosts (humans, pigs, dogs and rabbits) distributed in three countries (China, Australia and the USA) consistently supported genetic subdivisions according to host species rather than geographical location. The demographic history of S. scabiei reconstructed by multiple sequentially Markovian coalescent analysis suggested that S. scabiei isolated from rabbits, humans, dogs and pigs diverged ∼5000 years ago. Investigation of the homeobox (Hox) genes revealed that S. scabiei contains 8 of 10 canonical Hox genes that are present in the arthropod ancestor, and the absence of the Abd-A gene may correlate with the long gap between their front and back legs. Comparative genomics demonstrated that genes specific to scabies mites were mainly enriched in nutrition digestive systems, whereas genes in the families that involved detoxification (cytochrome P450, carboxyl/cholinesterases and the ATP-binding cassette transporter C group) were extremely contracted compared with that of other mites analysed in this study. Selective sweep analysis of mite populations from various pairs of two out of the four host species revealed that the strongest selective sweep signals were mainly enriched in cysteine-type peptidase activity and apoptosis. The results provided clues for the mechanisms of S. scabiei adaptation to a permanent parasitic lifestyle and knowledge that would enable further control of this highly contagious skin disease.

Plant biotrophic pathogens employ secreted molecules, called effectors, to suppress the host immune system and redirect the host's metabolism and development in their favour. Putative effectors of the gall-inducing maize pathogenic fungus Ustilago maydis were analysed for their ability to induce auxin signalling in plants. Using genetic, biochemical, cell-biological, and bioinformatic approaches we functionally elucidate a set of five, genetically linked effectors, called Topless (TPL) interacting protein (Tips) effectors that induce auxin signalling. We show that Tips induce auxin signalling by interfering with central corepressors of the TPL family. CRISPR-Cas9 mutants and deletion strain analysis indicate that the auxin signalling inducing subcluster effectors plays a redundant role in virulence. Although none of the Tips seem to have a conserved interaction motif, four of them bind solely to the N-terminal TPL domain and, for Tip1 and Tip4, we demonstrate direct competition with auxin/indole-3-acetic acid transcriptional repressors for their binding to TPL class of corepressors. Our findings reveal that TPL proteins, key regulators of growth-defence antagonism, are a major target of the U. maydis effectome.

In order to infect a new host species, the pathogen must evolve to enhance infection and transmission in the novel environment. Although we often think of evolution as a process of accumulation, it is also a process of loss. Here, we document an example of regressive evolution of an effector activity in the Irish potato famine pathogen (Phytophthora infestans) lineage, providing evidence that a key sequence motif in the effector PexRD54 has degenerated following a host jump. We began by looking at PexRD54 and PexRD54-like sequences from across Phytophthora species. We found that PexRD54 emerged in the common ancestor of Phytophthora clade 1b and 1c species, and further sequence analysis showed that a key functional motif, the C-terminal ATG8-interacting motif (AIM), was also acquired at this point in the lineage. A closer analysis showed that the P. mirabilis PexRD54 (PmPexRD54) AIM is atypical, the otherwise-conserved central residue mutated from a glutamate to a lysine. We aimed to determine whether this PmPexRD54 AIM polymorphism represented an adaptation to the Mirabilis jalapa host environment. We began by characterizing the M. jalapa ATG8 family, finding that they have a unique evolutionary history compared to previously characterized ATG8s. Then, using co-immunoprecipitation and isothermal titration calorimetry assays, we showed that both full-length PmPexRD54 and the PmPexRD54 AIM peptide bind weakly to the M. jalapa ATG8s. Through a combination of binding assays and structural modelling, we showed that the identity of the residue at the position of the PmPexRD54 AIM polymorphism can underpin high-affinity binding to plant ATG8s. Finally, we conclude that the functionality of the PexRD54 AIM was lost in the P. mirabilis lineage, perhaps owing to as-yet-unknown selection pressure on this effector in the new host environment.

Epigenetic mechanisms form the basis of cellular memory, developmental decisions, and the cellular immune system that defends against transposons and viruses. Organs develop from the shoot apical meristem (SAM) to shape the plant's areal phenotype, and stem cells in the SAM serve as a functional germline. While many details on the regulation of stem cell pool size, organ initiation, and patterning at the meristem periphery are known, we know surprisingly little about the molecular characteristics of SAM cells, including their epigenome and how it changes during development. Here, we summarize information on epigenetic regulation of selected genes necessary for stem cell maintenance. As recent evidence suggests that SAM stem cells might be a hotspot of transposon activation, we discuss this aspect of epigenetic control in the meristem and speculate on mechanisms that maintain the flexibility of SAM stem cells in response to developmental or environmental cues.

Plants show exceptional developmental plasticity and the ability to reprogram cell identities during regeneration. Although regeneration has been used in plant propagation for decades, we only recently gained detailed cellular and molecular insights into this process. Evidently, not all cell types have the same regeneration potential, and only a subset of regeneration-competent cells reach pluripotency. Pluripotent cells exhibit transcriptional similarity to root stem cells. In different plant regeneration systems, transcriptional reprogramming involves transient release of chromatin repression during pluripotency establishment and its restoration during organ or embryo differentiation. Incomplete resetting of the epigenome leads to somaclonal variation in regenerated plants. As single-cell technologies advance, we expect novel, exciting insights into epigenome dynamics during the establishment of pluripotency.

Disentangling the contribution of climatic and edaphic factors to microbiome variation and local adaptation in plants requires an experimental approach to uncouple their effects and test for causality. We used microbial inocula, soil matrices and plant genotypes derived from two natural Arabidopsis thaliana populations in northern and southern Europe in an experiment conducted in climatic chambers mimicking seasonal changes in temperature, day length and light intensity of the home sites of the two genotypes. The southern A. thaliana genotype outperformed the northern genotype in the southern climate chamber, whereas the opposite was true in the northern climate chamber. Recipient soil matrix, but not microbial composition, affected plant fitness, and effects did not differ between genotypes. Differences between chambers significantly affected rhizosphere microbiome assembly, although these effects were small in comparison with the shifts induced by physicochemical differences between soil matrices. The results suggest that differences in seasonal changes in temperature, day length and light intensity between northern and southern Europe have strongly influenced adaptive differentiation between the two A. thaliana populations, whereas effects of differences in soil factors have been weak. By contrast, below-ground differences in soil characteristics were more important than differences in climate for rhizosphere microbiome differentiation.

Natural plant populations often harbour substantial heritable variation in DNA methylation. However, a thorough understanding of the genetic and environmental drivers of this epigenetic variation requires large-scale and high-resolution data, which currently exist only for a few model species. Here, we studied 207 lines of the annual weed Thlaspi arvense (field pennycress), collected across a large latitudinal gradient in Europe and propagated in a common environment. By screening for variation in DNA sequence and DNA methylation using whole-genome (bisulfite) sequencing, we found significant epigenetic population structure across Europe. Average levels of DNA methylation were strongly context-dependent, with highest DNA methylation in CG context, particularly in transposable elements and in intergenic regions. Residual DNA methylation variation within all contexts was associated with genetic variants, which often co-localized with annotated methylation machinery genes but also with new candidates. Variation in DNA methylation was also significantly associated with climate of origin, with methylation levels being lower in colder regions and in more variable climates. Finally, we used variance decomposition to assess genetic versus environmental associations with differentially methylated regions (DMRs). We found that while genetic variation was generally the strongest predictor of DMRs, the strength of environmental associations increased from CG to CHG and CHH, with climate-of-origin as the strongest predictor in about one third of the CHH DMRs. In summary, our data show that natural epigenetic variation in Thlaspi arvense is significantly associated with both DNA sequence and environment of origin, and that the relative importance of the two factors strongly depends on the sequence context of DNA methylation. T. arvense is an emerging biofuel and winter cover crop; our results may hence be relevant for breeding efforts and agricultural practices in the context of rapidly changing environmental conditions.

This study demonstrates how the latest ultrahigh-performance liquid chromatography (UHPLC) technology can be combined with high-resolution accurate-mass (HRAM) mass spectrometry (MS) and long columns packed with fully porous particles to improve bottom-up proteomics analysis with nanoflow liquid chromatography-mass spectrometry (nanoLC-MS) methods. The increased back pressures from the UHPLC system enabled the use of 75 μm I.D. × 75 cm columns packed with 2 μm particles at a typical 300 nL/min flow rate as well as elevated and reduced flow rates. The constant pressure pump operation at 1500 bar reduced sample loading and column washing/equilibration stages and overall overhead time, which maximizes MS utilization time. The versatility of flow rate optimization to balance the sensitivity, throughput with sample loading amount, and capability of using longer gradients contributes to a greater number of peptide and protein identifications for single-shot bottom-up proteomics experiments. The routine proteome profiling and precise quantification of >7000 proteins with single-shot nanoLC-MS analysis open possibilities for large-scale discovery studies with a deep dive into the protein level alterations. Data are available via ProteomeXchange with identifier PXD035665.

The liverwort Marchantia polymorpha has been utilized as a model for biological studies since the 18th century. In the past few decades, there has been a Renaissance in its utilization in genomic and genetic approaches to investigating physiological, developmental, and evolutionary aspects of land plant biology. The reasons for its adoption are similar to those of other genetic models, e.g. simple cultivation, ready access via its worldwide distribution, ease of crossing, facile genetics, and more recently, efficient transformation, genome editing, and genomic resources. The haploid gametophyte dominant life cycle of M. polymorpha is conducive to forward genetic approaches. The lack of ancient whole-genome duplications within liverworts facilitates reverse genetic approaches, and possibly related to this genomic stability, liverworts possess sex chromosomes that evolved in the ancestral liverwort. As a representative of one of the three bryophyte lineages, its phylogenetic position allows comparative approaches to provide insights into ancestral land plants. Given the karyotype and genome stability within liverworts, the resources developed for M. polymorpha have facilitated the development of related species as models for biological processes lacking in M. polymorpha.

Whole-genome bisulfite sequencing (WGBS) is the standard method for profiling DNA methylation at single-nucleotide resolution. Different tools have been developed to extract differentially methylated regions (DMRs), often built upon assumptions from mammalian data. Here, we present MethylScore, a pipeline to analyse WGBS data and to account for the substantially more complex and variable nature of plant DNA methylation. MethylScore uses an unsupervised machine learning approach to segment the genome by classification into states of high and low methylation. It processes data from genomic alignments to DMR output and is designed to be usable by novice and expert users alike. We show how MethylScore can identify DMRs from hundreds of samples and how its data-driven approach can stratify associated samples without prior information. We identify DMRs in the A. thaliana 1,001 Genomes dataset to unveil known and unknown genotype-epigenotype associations .

Genome-wide association studies (GWAS) have revealed that the striking natural variation for DNA CHH-methylation (mCHH; H is A, T, or C) of transposons has oligogenic architecture involving major alleles at a handful of known methylation regulators. Here we use a conditional GWAS approach to show that CHG-methylation (mCHG) has a similar genetic architecture-once mCHH is statistically controlled for. We identify five key trans-regulators that appear to modulate mCHG levels, and show that they interact with a previously identified modifier of mCHH in regulating natural transposon mobilization.

High-throughput sequencing of amplicon libraries is the most widespread and one of the most effective ways to study the taxonomic structure of microbial communities, even despite growing accessibility of whole metagenome sequencing. Due to the targeted amplification, the method provides unparalleled resolution of communities, but at the same time perturbs initial community structure thereby reducing data robustness and compromising downstream analyses. Experimental research of the perturbations is largely limited to comparative studies on different PCR protocols without considering other sources of experimental variation related to characteristics of the initial microbial composition itself. Here we analyse these sources and demonstrate how dramatically they effect the relative abundances of taxa during the PCR cycles. We developed the mathematical model of the PCR amplification assuming the heterogeneity of amplification efficiencies and considering the compositional nature of data. We designed the experiment-five consecutive amplicon cycles (22-26) with 12 replicates for one real human stool microbial sample-and estimated the dynamics of the microbial community in line with the model. We found the high heterogeneity in amplicon efficiencies of taxa that leads to the non-linear and substantial (up to fivefold) changes in relative abundances during PCR. The analysis of possible sources of heterogeneity revealed the significant association between amplicon efficiencies and the energy of secondary structures of the DNA templates. The result of our work highlights non-trivial changes in the dynamics of real-life microbial communities due to their compositional nature. Obtained effects are specific not only for amplicon libraries, but also for any studies of metagenome dynamics.

Complex mechanisms regulate gene dosage throughout eukaryotic life cycles. Mechanisms controlling gene dosage have been extensively studied in animals, however it is unknown how generalizable these mechanisms are to diverse eukaryotes. Here, we use the haploid plant Marchantia polymorpha to assess gene dosage control in its short-lived diploid embryo. We show that throughout embryogenesis, paternal chromosomes are repressed resulting in functional haploidy. The paternal genome is targeted for genomic imprinting by the Polycomb mark H3K27me3 starting at fertilization, rendering the maternal genome in control of embryogenesis. Maintaining haploid gene dosage by this new form of imprinting is essential for embryonic development. Our findings illustrate how haploid-dominant species can regulate gene dosage through paternal chromosome inactivation and initiates the exploration of the link between life cycle history and gene dosage in a broader range of organisms.

BACKGROUND: Plants growing in proximity to other plants are exposed to a variety of metabolites that these neighbors release into the environment. Some species produce allelochemicals to inhibit growth of neighboring plants, which in turn have evolved ways to detoxify these compounds. RESULTS: In order to understand how the allelochemical-receiving target plants respond to chemically diverse compounds, we performed whole-genome transcriptome analysis of Arabidopsis thaliana exposed to either the benzoxazinoid derivative 2-amino- 3H-phenoxazin-3-one (APO) or momilactone B. These two allelochemicals belong to two very different compound classes, benzoxazinoids and diterpenes, respectively, produced by different Poaceae crop species. CONCLUSIONS: Despite their distinct chemical nature, we observed similar molecular responses of A. thaliana to these allelochemicals. In particular, many of the same or closely related genes belonging to the three-phase detoxification pathway were upregulated in both treatments. Further, we observed an overlap between genes upregulated by allelochemicals and those involved in herbicide detoxification. Our findings highlight the overlap in the transcriptional response of a target plant to natural and synthetic phytotoxic compounds and illustrate how herbicide resistance could arise via pathways involved in plant-plant interaction.

Meiosis, at the transition between diploid and haploid life cycle phases, is accompanied by reprograming of cell division machinery and followed by a transition back to mitosis. We show that, in Arabidopsis, this transition is driven by inhibition of translation, achieved by a mechanism that involves processing bodies (P-bodies). During the second meiotic division, the meiosis-specific protein THREE-DIVISION MUTANT 1 (TDM1) is incorporated into P-bodies through interaction with SUPPRESSOR WITH MORPHOGENETIC EFFECTS ON GENITALIA 7 (SMG7). TDM1 attracts eIF4F, the main translation initiation complex, temporarily sequestering it in P-bodies and inhibiting translation. The failure of tdm1 mutants to terminate meiosis can be overcome by chemical inhibition of translation. We propose that TDM1-containing P-bodies down-regulate expression of meiotic transcripts to facilitate transition of cell fates to postmeiotic gametophyte differentiation.

Iron is critical for host-pathogen interactions. While pathogens seek to scavenge iron to spread, the host aims at decreasing iron availability to reduce pathogen virulence. Thus, iron sensing and homeostasis are of particular importance to prevent host infection and part of nutritional immunity. While the link between iron homeostasis and immunity pathways is well established in plants, how iron levels are sensed and integrated with immune response pathways remains unknown. Here we report a receptor kinase SRF3, with a role in coordinating root growth, iron homeostasis and immunity pathways via regulation of callose synthases. These processes are modulated by iron levels and rely on SRF3 extracellular and kinase domains which tune its accumulation and partitioning at the cell surface. Mimicking bacterial elicitation with the flagellin peptide flg22 phenocopies SRF3 regulation upon low iron levels and subsequent SRF3-dependent responses. We propose that SRF3 is part of nutritional immunity responses involved in sensing external iron levels.

We investigated early vegetative growth of natural Arabidopsis thaliana accessions in cold, nonfreezing temperatures, similar to temperatures these plants naturally encounter in fall at northern latitudes. We found that accessions from northern latitudes produced larger seedlings than accessions from southern latitudes, partly as a result of larger seed size. However, their subsequent vegetative growth when exposed to colder temperatures was slower. The difference was too large to be explained by random population differentiation, and is thus suggestive of local adaptation, a notion that is further supported by substantial transcriptome and metabolome changes in northern accessions. We hypothesize that the reduced growth of northern accessions is an adaptive response and a consequence of reallocating resources toward cold acclimation and winter survival.

Although complex interactions between hosts and microbial associates are increasingly well documented, we still know little about how and why hosts shape microbial communities in nature. In addition, host genetic effects on microbial communities vary widely depending on the environment, obscuring conclusions about which microbes are impacted and which plant functions are important. We characterized the leaf microbiota of 200 Arabidopsis thaliana genotypes in eight field experiments and detected consistent host effects on specific, broadly distributed microbial species (operational taxonomic unit [OTUs]). Host genetic effects disproportionately influenced central ecological hubs, with heritability of particular OTUs declining with their distance from the nearest hub within the microbial network. These host effects could reflect either OTUs preferentially associating with specific genotypes or differential microbial success within them. Host genetics associated with microbial hubs explained over 10% of the variation in lifetime seed production among host genotypes across sites and years. We successfully cultured one of these microbial hubs and demonstrated its growth-promoting effects on plants in sterile conditions. Finally, genome-wide association mapping identified many putatively causal genes with small effects on the relative abundance of microbial hubs across sites and years, and these genes were enriched for those involved in the synthesis of specialized metabolites, auxins, and the immune system. Using untargeted metabolomics, we corroborate the consistent association between variation in specialized metabolites and microbial hubs across field sites. Together, our results reveal that host genetic variation impacts the microbial communities in consistent ways across environments and that these effects contribute to fitness variation among host genotypes.

Oligosaccharide fragments of fungal cell wall glycans are important molecular probes for studying both the biology of fungi and fungal infections of humans, animals, and plants. The fungal cell wall contains large amounts of various polysaccharides that are ligands for pattern recognition receptors (PRRs), eliciting an immune response upon recognition. Towards the establishment of a glycan array platform for the identification of new ligands of plant PRRs, tri-, penta-, and heptasaccharide fragments of different cell wall polysaccharides were prepared. Chito- and β-(1→6)-gluco-oligosaccharides were synthesized by automated glycan assembly (AGA), and α-(1→3)- and α-(1→4)-gluco-oligosaccharides were synthesized in solution using a recently reported highly α-selective glycosylation methodology. Incubation of plants with the synthesized oligosaccharides revealed i) length dependence for plant activation by chito-oligosaccharides and ii) β-1,6-glucan oligosaccharides as a new class of glycans capable of triggering plant activation.

Plant traits are an expression of strategic tradeoffs in plant performance that determine variation in allocation of finite resources to alternate physiological functions. Climate factors interact with plant traits to mediate tree survival. This study investigated survival dynamics in Norway spruce (Picea abies) in relation to tree-level morphological traits during a prolonged multi-year outbreak of the bark beetle, Ips typographus, in Central Europe. We acquired datasets describing the trait attributes of individual spruce using remote sensing and field surveys. We used nonlinear regression in a hypothesis-driven framework to quantify survival probability as a function of tree size, crown morphology, intraspecific competition and a growing season water balance. Extant spruce trees that persisted through the outbreak were spatially clustered, suggesting that survival was a non-random process. Larger diameter trees were more susceptible to bark beetles, reflecting either life history tradeoffs or a dynamic interaction between defense capacity and insect aggregation behavior. Competition had a strong negative effect on survival, presumably through resource limitation. Trees with more extensive crowns were buffered against bark beetles, ostensibly by a more robust photosynthetic capability and greater carbon reserves. The outbreak spanned a warming trend and conditions of anomalous aridity. Sustained water limitation during this period amplified the consequences of other factors, rendering even smaller trees vulnerable to colonization by insects. Our results are in agreement with prior research indicating that climate change has the potential to intensify bark beetle activity. However, forest outcomes will depend on complex cross-scale interactions between global climate trends and tree-level trait factors, as well as feedback effects associated with landscape patterns of stand structural diversity.