Search

Key messageSpecific combinations of LFY and WAPO1 natural alleles maximize spikelet number per spike in wheat. Spikelet number per spike (SNS) is an important yield component in wheat that determines the maximum number of grains that can be formed in a wheat spike. In wheat, loss-of-function mutations in LEAFY (LFY) or its interacting protein WHEAT ORTHOLOG OF APO1 (WAPO1) significantly reduce SNS by reducing the rate of formation of spikelet meristems. In previous studies, we identified a natural amino acid change in WAPO1 (C47F) that significantly increases SNS in hexaploid wheat. In this study, we searched for natural variants in LFY that were associated with differences in SNS and detected significant effects in the LFY-B region in a nested association mapping population. We generated a large mapping population and confirmed that the LFY-B polymorphism R80S is linked with the differences in SNS, suggesting that LFY-B is the likely causal gene. A haplotype analysis revealed two amino acid changes P34L and R80S, which were both enriched during wheat domestication and breeding suggesting positive selection. We also explored the interactions between the LFY and WAPO1 natural variants for SNS using biparental populations and identified significant interaction, in which the positive effect of the 80S and 34L alleles from LFY-B was only detected in the WAPO-A1 47F background but not in the 47C background. Based on these results, we propose that the allele combination WAPO-A1-47F/LFY-B 34L 80S can be used in wheat breeding programs to maximize SNS and increase grain yield potential in wheat.

Faithful genome partition during cell division relies on proper congression of chromosomes to the spindle equator before sister chromatid segregation. Here we uncover a kinesin-7 motor, kinetochore-associated kinesin 1 (KAK1), that is required for mitotic chromosome congression in Arabidopsis. KAK1 associates dynamically with kinetochores throughout mitosis. Loss of KAK1 results in severe defects in chromosome congression at metaphase, yet segregation errors at anaphase are rarely observed. KAK1 specifically interacts with the spindle assembly checkpoint protein BUB3.3 and both proteins show interdependent kinetochore localization. Chromosome misalignment in BUB3.3-depleted plants can be rescued by artificial tethering of KAK1 to kinetochores but not vice versa, demonstrating that KAK1 acts downstream of BUB3.3 to orchestrate microtubule-based chromosome transport at kinetochores. Moreover, we show that KAK1's motor activity is essential for driving chromosome congression to the metaphase plate. Thus, our findings reveal that plants have repurposed BUB3.3 to interface with a specialized kinesin adapted to integrate proper chromosome congression and checkpoint control through a distinct kinetochore design.

Improving photosynthesis, the fundamental process by which plants convert light energy into chemical energy, is a key area of research with great potential for enhancing sustainable agricultural productivity and addressing global food security challenges. This perspective delves into the latest advancements and approaches aimed at optimizing photosynthetic efficiency. Our discussion encompasses the entire process, beginning with light harvesting and its regulation and progressing through the bottleneck of electron transfer. We then delve into the carbon reactions of photosynthesis, focusing on strategies targeting the enzymes of the Calvin-Benson-Bassham (CBB) cycle. Additionally, we explore methods to increase carbon dioxide (CO2) concentration near the Rubisco, the enzyme responsible for the first step of CBB cycle, drawing inspiration from various photosynthetic organisms, and conclude this section by examining ways to enhance CO2 delivery into leaves. Moving beyond individual processes, we discuss two approaches to identifying key targets for photosynthesis improvement: systems modeling and the study of natural variation. Finally, we revisit some of the strategies mentioned above to provide a holistic view of the improvements, analyzing their impact on nitrogen use efficiency and on canopy photosynthesis.

This article focuses on screening the major secreted proteins by the ischemia-challenged cardiac stromal fibroblasts (CF), the assessment of their expression status and functional role in the post-ischemic left ventricle (LV) and in the ischemia-challenged CF culture and to phenotype CF at single cell resolution based on the positivity of the identified mediators. The expression level of CRSP2, HSP27, IL-8, Cofilin-1, and HSP90 in the LV tissues following coronary artery bypass graft (CABG) and myocardial infarction (MI) and CF cells followed the screening profile derived from the MS/MS findings. The histology data unveiled ECM disorganization, inflammation and fibrosis reflecting the ischemic pathology. CRSP2, HSP27, and HSP90 were significantly upregulated in the LV-CABG tissues with a concomitant reduction ion LV-MI whereas Cofilin-1, IL8, Nrf2, and Troponin I were downregulated in LV-CABG and increased in LV-MI. Similar trends were exhibited by ischemic CF. Single cell transcriptomics revealed multiple sub-phenotypes of CF based on their respective upregulation of CRSP2, HSP27, IL-8, Cofilin-1, HSP90, Troponin I and Nrf2 unveiling pathological and pro-healing phenotypes. Further investigations regarding the underlying signaling mechanisms and validation of sub-populations would offer novel translational avenues for the management of cardiac diseases.

Low iron (Fe) bioavailability can limit the biosynthesis of Fe-containing proteins, which are especially abundant in photosynthetic organisms, thus negatively affecting global primary productivity. Understanding cellular coping mechanisms under Fe limitation is therefore of great interest. We surveyed the temporal responses of Chlamydomonas (Chlamydomonas reinhardtii) cells transitioning from an Fe-rich to an Fe-free medium to document their short and long-term adjustments. While slower growth, chlorosis and lower photosynthetic parameters are evident only after one or more days in Fe-free medium, the abundance of some transcripts, such as those for genes encoding transporters and enzymes involved in Fe assimilation, change within minutes, before changes in intracellular Fe content are noticeable, suggestive of a sensitive mechanism for sensing Fe. Promoter reporter constructs indicate a transcriptional component to this immediate primary response. With acetate provided as a source of reduced carbon, transcripts encoding respiratory components are maintained relative to transcripts encoding components of photosynthesis and tetrapyrrole biosynthesis, indicating metabolic prioritization of respiration over photosynthesis. In contrast to the loss of chlorophyll, carotenoid content is maintained under Fe limitation despite a decrease in the transcripts for carotenoid biosynthesis genes, indicating carotenoid stability. These changes occur more slowly, only after the intracellular Fe quota responds, indicating a phased response in Chlamydomonas, involving both primary and secondary responses during acclimation to poor Fe nutrition.

Gene expression can be highly plastic in response to environmental variation. However, we know little about how expression plasticity is shaped by natural selection and evolves in wild and domesticated species. We used genotypic selection analysis to characterize selection on drought-induced plasticity of over 7,500 leaf transcripts of 118 rice accessions (genotypes) from different environmental conditions grown in a field experiment. Gene expression plasticity was neutral for most gradually plastic transcripts, but transcripts with discrete patterns of expression showed stronger selection on expression plasticity. Whether plasticity was adaptive and co-gradient or maladaptive and counter-gradient varied among varietal groups. No transcripts that experienced selection for plasticity across environments showed selection against plasticity within environments, indicating a lack of evidence for costs of adaptive plasticity that may constrain its evolution. Selection on expression plasticity was influenced by degree of plasticity, transcript length and gene body methylation. We observed positive selection on plasticity of co-expression modules containing transcripts involved in photosynthesis, translation and responsiveness to abiotic stress. Taken together, these results indicate that patterns of selection on expression plasticity were context-dependent and likely associated with environmental conditions of varietal groups, but that the evolution of adaptive plasticity would likely not be constrained by opposing patterns of selection on plasticity within compared to across environments. These results offer a genome-wide view of patterns of selection and ecological constraints on gene expression plasticity and provide insights into the interplay between plastic and evolutionary responses to drought at the molecular level.

Modification of lignin in feedstocks via genetic engineering aims to reduce biomass recalcitrance to facilitate efficient conversion processes. These improvements can be achieved by expressing exogenous enzymes that interfere with native biosynthetic pathways responsible for the production of the lignin precursors. In planta expression of a bacterial 3-dehydroshikimate dehydratase in poplar trees reduced lignin content and altered the monomer composition, which enabled higher yields of sugars after cell wall polysaccharide hydrolysis. Understanding how plants respond to such genetic modifications at the transcriptional and metabolic levels is needed to facilitate further improvement and field deployment. In this work, we acquired fundamental knowledge on lignin-modified poplar expressing 3-dehydroshikimate dehydratase using RNA-seq and metabolomics. The data clearly demonstrate that changes in gene expression and metabolite abundance can occur in a strict spatiotemporal fashion, revealing tissue-specific responses in the xylem, phloem, or periderm. In the poplar line that exhibited the strongest reduction in lignin, we found that 3% of the transcripts had altered expression levels and ~19% of the detected metabolites had differential abundance in the xylem from older stems. The changes affected predominantly the shikimate and phenylpropanoid pathways as well as secondary cell wall metabolism, and resulted in significant accumulation of hydroxybenzoates derived from protocatechuate and salicylate.

Dissecting plant responses to the environment is key to understanding whether and how plants adapt to anthropogenic climate change. Stomata, plants' pores for gas exchange, are expected to decrease in density following increased CO2 concentrations, a trend already observed in multiple plant species. However, it is unclear whether such responses are based on genetic changes and evolutionary adaptation. Here we make use of extensive knowledge of 43 genes in the stomatal development pathway and newly generated genome information of 191 Arabidopsis thaliana historical herbarium specimens collected over 193 years to directly link genetic variation with climate change. While we find that the essential transcription factors SPCH, MUTE and FAMA, central to stomatal development, are under strong evolutionary constraints, several regulators of stomatal development show signs of local adaptation in contemporary samples from different geographic regions. We then develop a functional score based on known effects of gene knock-out on stomatal development that recovers a classic pattern of stomatal density decrease over the past centuries, suggesting a genetic component contributing to this change. This approach combining historical genomics with functional experimental knowledge could allow further investigations of how different, even in historical samples unmeasurable, cellular plant phenotypes may have already responded to climate change through adaptive evolution.

Abstract: RNA methylation plays a central regulatory role in plant biology and is a relatively new target for plant improvement efforts. In nearly all cases, perturbation of the RNA methylation machinery results in deleterious phenotypes. However, a recent landmark paper reported that transcriptome‐wide use of the human RNA demethylase FTO substantially increased the yield of rice and potatoes. Here, we have performed the first independent replication of those results and demonstrated broader transferability of the trait, finding increased flower and fruit count in the model species Arabidopsis thaliana. We also performed RNA‐seq of our FTO‐transgenic plants, which we analyzed in conjunction with previously published datasets to detect several previously unrecognized patterns in the functional and structural classification of the upregulated and downregulated genes. From these, we present mechanistic hypotheses to explain these surprising results with the goal of spurring more widespread interest in this promising new approach to plant engineering.

Schizophyllum commune is a mushroom-forming fungus notable for its distinctive fruiting bodies with split gills. It is used as a model organism to study mushroom development, lignocellulose degradation and mating type loci. It is a hypervariable species with considerable genetic and phenotypic diversity between the strains. In this study, we systematically phenotyped 16 dikaryotic strains for aspects of mushroom development and 18 monokaryotic strains for lignocellulose degradation. There was considerable heterogeneity among the strains regarding these phenotypes. The majority of the strains developed mushrooms with varying morphologies, although some strains only grew vegetatively under the tested conditions. Growth on various carbon sources showed strain-specific profiles. The genomes of seven monokaryotic strains were sequenced and analyzed together with six previously published genome sequences. Moreover, the related species Schizophyllum fasciatum was sequenced. Although there was considerable genetic variation between the genome assemblies, the genes related to mushroom formation and lignocellulose degradation were well conserved. These sequenced genomes, in combination with the high phenotypic diversity, will provide a solid basis for functional genomics analyses of the strains of S. commune.

The red and far-red light photoreceptor phytochrome B (phyB) transmits light signals following cytosol-to-nuclear translocation to regulate transcriptional networks therein. This necessitates changes in protein-protein interactions of phyB in the cytosol, about which little is presently known. Via introduction of a nucleus-excluding G767R mutation into the dominant, constitutively active phyBY276H (YHB) allele, we explore the functional consequences of expressing a cytosol-localized YHBG767R variant in transgenic Arabidopsis seedlings. We show that YHBG767R elicits selective constitutive photomorphogenic phenotypes in dark-grown phyABCDE null mutants, wild type and other phy-deficient genotypes. These responses include light-independent apical hook opening, cotyledon unfolding, seed germination and agravitropic hypocotyl growth with minimal suppression of hypocotyl elongation. Such phenotypes correlate with reduced PIF3 levels, which implicates cytosolic targeting of PIF3 turnover or PIF3 translational inhibition by YHBG767R. However, as expected for a cytoplasm-tethered phyB, YHBG767R elicits reduced light-mediated signaling activity compared with similarly expressed wild-type phyB in phyABCDE mutant backgrounds. YHBG767R also interferes with wild-type phyB light signaling, presumably by formation of cytosol-retained and/or otherwise inactivated heterodimers. Our results suggest that cytosolic interactions with PIFs play an important role in phyB signaling even under physiological conditions.

Tobacco mild green mosaic virus (TMGMV)-like nanocarriers were designed for gene delivery to plant cells. High aspect ratio TMGMVs were coated with a polycationic biopolymer, poly(allylamine) hydrochloride (PAH), to generate highly charged nanomaterials (TMGMV-PAH; 56.20 ± 4.7 mV) that efficiently load (1:6 TMGMV:DNA mass ratio) and deliver single-stranded and plasmid DNA to plant cells. The TMGMV-PAH were taken up through energy-independent mechanisms in Arabidopsis protoplasts. TMGMV-PAH delivered a plasmid DNA encoding a green fluorescent protein (GFP) to the protoplast nucleus (70% viability), as evidenced by GFP expression using confocal microscopy and Western blot analysis. TMGMV-PAH were inactivated (iTMGMV-PAH) using UV cross-linking to prevent systemic infection in intact plants. Inactivated iTMGMV-PAH-mediated pDNA delivery and gene expression of GFP in vivo was determined using confocal microscopy and RT-qPCR. Virus-like nanocarrier-mediated gene delivery can act as a facile and biocompatible tool for advancing genetic engineering in plants.

Plant-microbe symbioses require intense interaction and genetic coordination to successfully establish in specific cell types of the host and symbiont. Traditional RNA-seq methodologies lack the cellular resolution to fully capture these complexities, but single-cell and spatial transcriptomics (ST) are now allowing scientists to probe symbiotic interactions at an unprecedented level of detail. Here, we discuss the advantages that novel spatial and single-cell transcriptomic technologies provide in studying plant-microbe endosymbioses and highlight key recent studies. Finally, we consider the remaining limitations of applying these approaches to symbiosis research, which are mainly related to the simultaneous capture of both plant and microbial transcripts within the same cells.

Arsenic contamination of soils threatens the health of millions globally through accumulation in crops. While plants detoxify arsenic via phytochelatin (PC) complexation and efflux of arsenite from roots, arsenite efflux mechanisms are not fully understood. Here, white lupin (Lupinus albus) was grown in semi-hydroponics and exudation of glutathione (GSH) derivatives and PCs in response to arsenic was scrutinised using LC-MS/MS. Inhibiting synthesis of PC precursor GSH with L-buthionine sulfoximine (BSO) or ABC transporters with vanadate drastically reduced (>22%) GSH-derivative and PC2 exudation, but not PC3 exudation. This was accompanied by arsenic hypersensitivity in plants treated with BSO and moderate sensitivity with vanadate treatment. Investigating arsenic-phytochelatin (As-PC) complexation revealed two distinct As-PC complexes, As bound to GSH and PC2 (GS-As-PC2) and As bound to PC3 (As-PC3), in exudates of As-treated lupin. Vanadate inhibited As-PC exudation, while BSO inhibited both the synthesis and exudation of As-PC complexes. These results demonstrate a role of GSH-derivatives and PC exudation in lupin arsenic tolerance and reveal As-PC exudation as a new potential mechanism contributing to active arsenic efflux in plants. Overall, this study uncovers insight into rhizosphere arsenic detoxification with potential to help mitigate pollution and reduce arsenic accumulation in crops.

Zygnematophyceae are the algal sisters of land plants. Here we sequenced four genomes of filamentous Zygnematophyceae, including chromosome-scale assemblies for three strains of Zygnema circumcarinatum. We inferred traits in the ancestor of Zygnematophyceae and land plants that might have ushered in the conquest of land by plants: expanded genes for signaling cascades, environmental response, and multicellular growth. Zygnematophyceae and land plants share all the major enzymes for cell wall synthesis and remodifications, and gene gains shaped this toolkit. Co-expression network analyses uncover gene cohorts that unite environmental signaling with multicellular developmental programs. Our data shed light on a molecular chassis that balances environmental response and growth modulation across more than 600 million years of streptophyte evolution.

Soil fungi belonging to different functional guilds, such as saprotrophs, pathogens, and mycorrhizal symbionts, play key roles in forest ecosystems. To date, no study has compared the actual gene expression of these guilds in different forest soils. We used metatranscriptomics to study the competition for organic resources by these fungal groups in boreal, temperate, and Mediterranean forest soils. Using a dedicated mRNA annotation pipeline combined with the JGI MycoCosm database, we compared the transcripts of these three fungal guilds, targeting enzymes involved in C- and N mobilization from plant and microbial cell walls. Genes encoding enzymes involved in the degradation of plant cell walls were expressed at a higher level in saprotrophic fungi than in ectomycorrhizal and pathogenic fungi. However, ectomycorrhizal and saprotrophic fungi showed similarly high expression levels of genes encoding enzymes involved in fungal cell wall degradation. Transcripts for N-related transporters were more highly expressed in ectomycorrhizal fungi than in other groups. We showed that ectomycorrhizal and saprotrophic fungi compete for N in soil organic matter, suggesting that their interactions could decelerate C cycling. Metatranscriptomics provides a unique tool to test controversial ecological hypotheses and to better understand the underlying ecological processes involved in soil functioning and carbon stabilization.

Abstract: Tomato fruit ripening is triggered by the demethylation of key genes, which alters their transcriptional levels thereby initiating and propagating a cascade of physiological events. What is unknown, is how these processes are altered when fruit are ripened using postharvest practices to extend shelf-life, as these practices often reduce fruit quality. To address this, postharvest handling-induced changes in the fruit DNA methylome and transcriptome, and how they correlate with ripening speed, and ripening indicators such as ethylene, ABA, and carotenoids, were assessed. This study comprehensively connected changes in physiological events with dynamic molecular changes. Ripening fruit that reached ‘Turning’ (T) after dark-storage at 20°C, 12.5°C, or 5°C chilling (followed by 20°C rewarming), were compared to fresh-harvest fruit ‘FHT’. Fruit stored at 12.5°C, had the biggest epigenetic marks and alterations in gene expression, exceeding changes induced by postharvest chilling. Fruit physiological and chronological age were uncoupled at 12.5°C, as the time-to-ripening was the longest. Fruit ripening to Turning at 12.5°C was not climacteric; there was no respiratory or ethylene burst, rather, fruit were high in ABA. Clear differentiation between postharvest-ripened and ‘FHT’ was evident in the methylome and transcriptome. Higher expression of photosynthetic genes and chlorophyll levels in ‘FHT’ fruit, pointed to light as influencing the molecular changes in fruit ripening. Finally, correlative analyses of the -omics data putatively identified genes regulated by DNA methylation. Collectively, these data improve our interpretation of how tomato fruit ripening patterns are altered by postharvest practices, and long-term are expected to help improve fruit quality.

The symbiotic interaction of plants with arbuscular mycorrhizal (AM) fungi is ancient and widespread. Plants provide AM fungi with carbon in exchange for nutrients and water, making this interaction a prime target for crop improvement. However, plant-fungal interactions are restricted to a small subset of root cells, precluding the application of most conventional functional genomic techniques to study the molecular bases of these interactions. Here we used single-nucleus and spatial RNA sequencing to explore both Medicago truncatula and Rhizophagus irregularis transcriptomes in AM symbiosis at cellular and spatial resolution. Integrated, spatially registered single-cell maps revealed infected and uninfected plant root cell types. We observed that cortex cells exhibit distinct transcriptome profiles during different stages of colonization by AM fungi, indicating dynamic interplay between both organisms during establishment of the cellular interface enabling successful symbiosis. Our study provides insight into a symbiotic relationship of major agricultural and environmental importance and demonstrates a paradigm combining single-cell and spatial transcriptomics for the analysis of complex organismal interactions.

The Pseudomonas syringae species complex is a heterogeneous group of plant pathogenic bacteria associated with a wide distribution of plant species. Advances in genomics are revealing the complex evolutionary history of this species complex and the wide array of genetic adaptations underpinning their diverse lifestyles. Here, we genomically characterize two P. syringae isolates collected from diseased Callery pears (Pyrus calleryana) in Berkeley, California in 2019 and 2022. We also isolated a lytic bacteriophage, which we characterized and evaluated for biocontrol efficiency. Using a multilocus sequence analysis and core genome alignment, we classified the P. syringae isolates as members of phylogroup 2, related to other strains previously isolated from Pyrus and Prunus. An analysis of effector proteins demonstrated an evolutionary conservation of effectoromes across isolates classified in PG2 and yet uncovered unique effector profiles for each, including the two newly identified isolates. Whole-genome sequencing of the associated phage uncovered a novel phage genus related to Pseudomonas syringae pv. actinidiae phage PHB09 and the Flaumdravirus genus. Finally, using in planta infection assays, we demonstrate that the phage was equally useful in symptom mitigation of immature pear fruit regardless of the Pss strain tested. Overall, this study demonstrates the diversity of P. syringae and their viruses associated with ornamental pear trees, posing spill-over risks to commercial pear trees and the possibility of using phages as biocontrol agents to reduce the impact of disease.IMPORTANCEGlobal change exacerbates the spread and impact of pathogens, especially in agricultural settings. There is a clear need to better monitor the spread and diversity of plant pathogens, including in potential spillover hosts, and for the development of novel and sustainable control strategies. In this study, we characterize the first described strains of Pseudomonas syringae pv. syringae isolated from Callery pear in Berkeley, California from diseased tissues in an urban environment. We show that these strains have divergent virulence profiles from previously described strains and that they can cause disease in commercial pears. Additionally, we describe a novel bacteriophage that is associated with these strains and explore its potential to act as a biocontrol agent. Together, the data presented here demonstrate that ornamental pear trees harbor novel P. syringae pv. syringae isolates that potentially pose a risk to local fruit production, or vice versa-but also provide us with novel associated phages, effective in disease mitigation.

The entomopathogenic fungus Conoideocrella luteorostrata has recently been implicated in natural epizootics among exotic elongate hemlock scale (EHS) insects in Christmas tree farms in the eastern United States. Since 1913, C. luteorostrata has been reported from various plant feeding Hemiptera in the southeastern United States, but comprehensive morphological and phylogenetic studies of U.S. populations are lacking. The recovery of multiple strains of C. luteorostrata from mycosed EHS in North Carolina provided an opportunity to conduct pathogenicity assays and morphological and phylogenetic studies to investigate genus- and species-level boundaries among the Clavicipitaceae. Pathogenicity assays confirmed that C. luteorostrata causes mortality of EHS crawlers, an essential first step in developing this fungus as a biocontrol. Morphological studies revealed that conidia aligned with previous measurements of the Paecilomyces-like asexual state of C. luteorostrata, with conidiophore morphology consistent with historical observations. Additionally, a Hirsutella-like synanamorph was observed in select C. luteorostrata strains. In both a four-locus, 54-taxon Clavicipitaceae-wide phylogenetic analysis including D1-D2 domains of the nuclear 28S rRNA region (28S), elongation factor 1 alpha (EF1-α), DNA-directed RNA polymerase II subunit 1 (RPB1), and DNA-directed RNA polymerase II subunit 2 (RPB2) and a two-locus, 38-taxon (28S and EF1-α) phylogenetic analysis, all three Conoideocrella species were resolved as strongly supported monophyletic lineages across all loci and both methods (maximum likelihood and Bayesian inference) of phylogenetic inference except for 28S for C. tenuis. Despite the strong support for individual Conoideocrella species, none of the analyses supported the monophyly of Conoideocrella with the inclusion of Dussiella. Due to the paucity of RPB1 and RPB2 sequence data, EF1-α provided superior delimitation of intraspecies groupings for Conoideocrella and should be used in future studies. Further development of C. luteorostrata as a biocontrol against EHS will require additional surveys across diverse Hemiptera and expanded pathogenicity testing to clarify host range and efficacy of this fungus.

Long terminal repeat retrotransposons (LTR-RTs) are powerful mutagens regarded as a major source of genetic novelty and important drivers of evolution. Yet, the uncontrolled and potentially selfish proliferation of LTR-RTs can lead to deleterious mutations and genome instability, with large fitness costs for their host. While population genomics data suggest that an ongoing LTR-RT mobility is common in many species, the understanding of their dual role in evolution is limited. Here, we harness the genetic diversity of 320 sequenced natural accessions of the Mediterranean grass Brachypodium distachyon to characterize how genetic and environmental factors influence plant LTR-RT dynamics in the wild. When combining a coverage-based approach to estimate global LTR-RT copy number variations with mobilome-sequencing of nine accessions exposed to eight different stresses, we find little evidence for a major role of environmental factors in LTR-RT accumulations in B. distachyon natural accessions. Instead, we show that loss of RNA polymerase IV (Pol IV), which mediates RNA-directed DNA methylation in plants, results in high transcriptional and transpositional activities of RLC_BdisC024 (HOPPLA) LTR-RT family elements, and that these effects are not stress-specific. This work supports findings indicating an ongoing mobility in B. distachyon and reveals that host RNA-directed DNA methylation rather than environmental factors controls their mobility in this wild grass model.

Mutant populations are crucial for functional genomics and discovering novel traits for crop breeding. Sorghum, a drought and heat-tolerant C4 species, requires a vast, large-scale, annotated, and sequenced mutant resource to enhance crop improvement through functional genomics research. Here, we report a sorghum large-scale sequenced mutant population with 9.5 million ethyl methane sulfonate (EMS)-induced mutations that covered 98% of sorghum's annotated genes using inbred line BTx623. Remarkably, a total of 610 320 mutations within the promoter and enhancer regions of 18 000 and 11 790 genes, respectively, can be leveraged for novel research of cis-regulatory elements. A comparison of the distribution of mutations in the large-scale mutant library and sorghum association panel (SAP) provides insights into the influence of selection. EMS-induced mutations appeared to be random across different regions of the genome without significant enrichment in different sections of a gene, including the 5' UTR, gene body, and 3'-UTR. In contrast, there were low variation density in the coding and UTR regions in the SAP. Based on the Ka /Ks value, the mutant library (~1) experienced little selection, unlike the SAP (0.40), which has been strongly selected through breeding. All mutation data are publicly searchable through SorbMutDB (https://www.depts.ttu.edu/igcast/sorbmutdb.php) and SorghumBase (https://sorghumbase.org/). This current large-scale sequence-indexed sorghum mutant population is a crucial resource that enriched the sorghum gene pool with novel diversity and a highly valuable tool for the Poaceae family, that will advance plant biology research and crop breeding.

Several closely related Myb-like activator proteins are known to have partially redundant functions within the plant circadian clock, but their specific roles are not well understood. To clarify the function of the REVEILLE 4, REVEILLE 6, and REVEILLE 8 transcriptional activators, we characterized the growth and clock phenotypes of CRISPR-Cas9-generated single, double, and triple rve mutants. We found that these genes act synergistically to regulate flowering time, redundantly to regulate leaf growth, and antagonistically to regulate hypocotyl elongation. We previously reported that increasing intensities of monochromatic blue and red light have opposite effects on the period of triple rve468 mutants. Here, we further examined light quality-specific phenotypes of rve mutants and report that rve468 mutants lack the blue light-specific increase in expression of some circadian clock genes observed in wild type. To investigate the basis of these blue light-specific circadian phenotypes, we examined RVE protein abundances and degradation rates in blue and red light and found no significant differences between these conditions. We next examined genetic interactions between RVE genes and ZEITLUPE and ELONGATED HYPOCOTYL5, two factors with blue light-specific functions in the clock. We found that the RVEs interact additively with both ZEITLUPE and ELONGATED HYPOCOTYL5 to regulate circadian period, which suggests that neither of these factors are required for the blue light-specific differences that we observed. Overall, our results suggest that the RVEs have separable functions in plant growth and circadian regulation and that they are involved in blue light-specific circadian signaling via a novel mechanism.

Cowpea, Vigna unguiculata L. Walp., is a diploid warm-season legume of critical importance as both food and fodder in sub-Saharan Africa. This species is also grown in Northern Africa, Europe, Latin America, North America, and East to Southeast Asia. To capture the genomic diversity of domesticates of this important legume, de novo genome assemblies were produced for representatives of six subpopulations of cultivated cowpea identified previously from genotyping of several hundred diverse accessions. In the most complete assembly (IT97K-499-35), 26,026 core and 4963 noncore genes were identified, with 35,436 pan genes when considering all seven accessions. GO terms associated with response to stress and defense response were highly enriched among the noncore genes, while core genes were enriched in terms related to transcription factor activity, and transport and metabolic processes. Over 5 million single nucleotide polymorphisms (SNPs) relative to each assembly and over 40 structural variants >1 Mb in size were identified by comparing genomes. Vu10 was the chromosome with the highest frequency of SNPs, and Vu04 had the most structural variants. Noncore genes harbor a larger proportion of potentially disruptive variants than core genes, including missense, stop gain, and frameshift mutations; this suggests that noncore genes substantially contribute to diversity within domesticated cowpea.

We present chromosome-level genome assemblies from representative species of three independently evolved seagrass lineages: Posidonia oceanica, Cymodocea nodosa, Thalassia testudinum and Zostera marina. We also include a draft genome of Potamogeton acutifolius, belonging to a freshwater sister lineage to Zosteraceae. All seagrass species share an ancient whole-genome triplication, while additional whole-genome duplications were uncovered for C. nodosa, Z. marina and P. acutifolius. Comparative analysis of selected gene families suggests that the transition from submerged-freshwater to submerged-marine environments mainly involved fine-tuning of multiple processes (such as osmoregulation, salinity, light capture, carbon acquisition and temperature) that all had to happen in parallel, probably explaining why adaptation to a marine lifestyle has been exceedingly rare. Major gene losses related to stomata, volatiles, defence and lignification are probably a consequence of the return to the sea rather than the cause of it. These new genomes will accelerate functional studies and solutions, as continuing losses of the 'savannahs of the sea' are of major concern in times of climate change and loss of biodiversity.

Bud dormancy is a crucial strategy for perennial plants to withstand adverse winter conditions. However, the regulatory mechanism of bud dormancy in tree peony (Paeonia suffruticosa) remains largely unknown. Here, we observed dramatically reduced and increased accumulation of abscisic acid (ABA) and bioactive gibberellins (GAs) GA1 and GA3, respectively, during bud endo-dormancy release of tree peony under prolonged chilling treatment. An Illumina RNA sequencing study was performed to identify potential genes involved in the bud endo-dormancy regulation in tree peony. Correlation matrix, principal component, and interaction network analyses identified a down-regulated MYB transcription factor gene, PsMYB306, the expression of which positively correlated with 9-CIS-EPOXYCAROTENOID DIOXYGENASE 3 (PsNCED3) expression. Protein modeling analysis revealed four residues within the R2R3 domain of PsMYB306 to possess DNA binding capability. Transcription of PsMYB306 was increased by ABA treatment. Overexpression of PsMYB306 in petunia (Petunia hybrida) inhibited seed germination and plant growth, concomitant with elevated ABA and decreased GA contents. Silencing of PsMYB306 accelerated cold-triggered tree peony bud burst, and influenced the production of ABA and GAs and the expression of their biosynthetic genes. ABA application reduced bud dormancy release and transcription of ENT-KAURENOIC ACID OXIDASE 1 (PsKAO1), GA20-OXIDASE 1 (PsGA20ox1), and GA3-OXIDASE 1 (PsGA3ox1) associated with GA biosynthesis in PsMYB306-silenced buds. In vivo and in vitro binding assays confirmed that PsMYB306 specifically transactivated the promoter of PsNCED3. Silencing of PsNCED3 also promoted bud break and growth. Altogether, our findings suggest that PsMYB306 negatively modulates cold-induced bud endo-dormancy release by regulating ABA production.

The kinetochore scaffold 1 (KNL1) protein recruits spindle assembly checkpoint (SAC) proteins to ensure accurate chromosome segregation during mitosis. Despite such a conserved function among eukaryotic organisms, its molecular architectures have rapidly evolved so that the functional mode of plant KNL1 is largely unknown. To understand how SAC signaling is regulated at kinetochores, we characterized the function of the KNL1 gene in Arabidopsis thaliana. The KNL1 protein was detected at kinetochores throughout the mitotic cell cycle, and null knl1 mutants were viable and fertile but exhibited severe vegetative and reproductive defects. The mutant cells showed serious impairments of chromosome congression and segregation, that resulted in the formation of micronuclei. In the absence of KNL1, core SAC proteins were no longer detected at the kinetochores, and the SAC was not activated by unattached or misaligned chromosomes. Arabidopsis KNL1 interacted with SAC essential proteins BUB3.3 and BMF3 through specific regions that were not found in known KNL1 proteins of other species, and recruited them independently to kinetochores. Furthermore, we demonstrated that upon ectopic expression, the KNL1 homolog from the dicot tomato was able to functionally substitute KNL1 in A. thaliana, while others from the monocot rice or moss associated with kinetochores but were not functional, as reflected by sequence variations of the kinetochore proteins in different plant lineages. Our results brought insights into understanding the rapid evolution and lineage-specific connection between KNL1 and the SAC signaling molecules.

Photoprotection mechanisms are ubiquitous among photosynthetic organisms. The photoprotection capacity of the green alga Chlamydomonas reinhardtii is correlated with protein levels of stress-related light-harvesting complex (LHCSR) proteins, which are strongly induced by high light (HL). However, the dynamic response of overall thylakoid structure during acclimation to growth in HL has not been fully understood. Here, we combined live-cell super-resolution microscopy and analytical membrane subfractionation to investigate macroscale structural changes of thylakoid membranes during HL acclimation in Chlamydomonas. Subdiffraction-resolution live-cell imaging revealed that the overall thylakoid structures became thinned and shrunken during HL acclimation. The stromal space around the pyrenoid also became enlarged. Analytical density-dependent membrane fractionation indicated that the structural changes were partly a consequence of membrane unstacking. The analysis of both an LHCSR loss-of-function mutant, npq4 lhcsr1, and a regulatory mutant that over-expresses LHCSR, spa1-1, showed that structural changes occurred independently of LHCSR protein levels, demonstrating that LHCSR was neither necessary nor sufficient to induce the thylakoid structural changes associated with HL acclimation. In contrast, stt7-9, a mutant lacking a kinase of major light-harvesting antenna proteins, had a slower thylakoid structural response to HL relative to all other lines tested but still showed membrane unstacking. These results indicate that neither LHCSR- nor antenna-phosphorylation-dependent HL acclimation are required for the observed macroscale structural changes of thylakoid membranes in HL conditions.

BackgroundColletotrichum fungi infect a wide diversity of monocot and dicot hosts, causing diseases on almost all economically important plants worldwide. Colletotrichum is also a suitable model for studying gene family evolution on a fine scale to uncover events in the genome associated with biological changes.ResultsHere we present the genome sequences of 30 Colletotrichum species covering the diversity within the genus. Evolutionary analyses revealed that the Colletotrichum ancestor diverged in the late Cretaceous in parallel with the diversification of flowering plants. We provide evidence of independent host jumps from dicots to monocots during the evolution of Colletotrichum, coinciding with a progressive shrinking of the plant cell wall degradative arsenal and expansions in lineage-specific gene families. Comparative transcriptomics of 4 species adapted to different hosts revealed similarity in gene content but high diversity in the modulation of their transcription profiles on different plant substrates. Combining genomics and transcriptomics, we identified a set of core genes such as specific transcription factors, putatively involved in plant cell wall degradation.ConclusionsThese results indicate that the ancestral Colletotrichum were associated with dicot plants and certain branches progressively adapted to different monocot hosts, reshaping the gene content and its regulation.

Photosystem I (PSI) serves as a model system for studying fundamental processes such as electron transfer (ET) and energy conversion, which are not only central to photosynthesis but also have broader implications for bioenergy production and biomimetic device design. In this study, we employed electron paramagnetic resonance (EPR) spectroscopy to investigate key light-induced charge separation steps in PSI isolated from several green algal and cyanobacterial species. Following photoexcitation, rapid sequential ET occurs through either of two quasi-symmetric branches of donor/acceptor cofactors embedded within the protein core, termed the A and B branches. Using high-frequency (130 GHz) time-resolved EPR (TR-EPR) and deuteration techniques to enhance spectral resolution, we observed that at low temperatures prokaryotic PSI exhibits reversible ET in the A branch and irreversible ET in the B branch, while PSI from eukaryotic counterparts displays either reversible ET in both branches or exclusively in the B branch. Furthermore, we observed a notable correlation between low-temperature charge separation to the terminal [4Fe-4S] clusters of PSI, termed FA and FB, as reflected in the measured FA/FB ratio. These findings enhance our understanding of the mechanistic diversity of PSI's ET across different species and underscore the importance of experimental design in resolving these differences. Though further research is necessary to elucidate the underlying mechanisms and the evolutionary significance of these variations in PSI charge separation, this study sets the stage for future investigations into the complex interplay between protein structure, ET pathways, and the environmental adaptations of photosynthetic organisms.

Phytophthora cinnamomi is a hemibiotrophic oomycete causing Phytophthora root rot in over 5,000 plant species, threatening natural ecosystems, forestry, and agriculture. Genomic studies of P. cinnamomi are limited compared to other Phytophthora spp. despite the importance of this destructive and highly invasive pathogen. The genome of two genetically and phenotypically distinct P. cinnamomi isolates collected from avocado orchards in California were sequenced using PacBio and Illumina sequencing. Genome sizes were estimated by flow cytometry and assembled de novo to 140-141 Mb genomes with 21,111-21,402 gene models. Genome analyses revealed that both isolates exhibited complex heterozygous genomes fitting the two-speed genome model. The more virulent isolate encodes a larger secretome and more RXLR effectors when compared to the less virulent isolate. Transcriptome analysis after P. cinnamomi infection in Arabidopsis thaliana, Nicotiana benthamiana, and Persea americana de Mill (avocado) showed that this pathogen deploys common gene repertoires in all hosts and host-specific subsets, especially among effectors. Overall, our results suggested that clonal P. cinnamomi isolates employ similar strategies as other Phytophthora spp. to increase phenotypic diversity (e.g., polyploidization, gene duplications, and a bipartite genome architecture) to cope with environmental changes. Our study also provides insights into common and host-specific P. cinnamomi infection strategies and may serve as a method for narrowing and selecting key candidate effectors for functional studies to determine their contributions to plant resistance or susceptibility.

Recently, the siRNAs pathways, and especially reproductive phasiRNAs, have attracted attention in eudicots since their biological roles are still unknown and their biogenesis took different evolutionary pathways compared to monocots. In this work, we used Coffea arabica L., a recent allotetraploid formed from the hybridization of Coffea canephora and C. eugenioides unreduced gametes, to explore microsporogenesis and small RNAs-related pathways in a eudicot crop. First, we identified the microsporogenesis stages during anther development revealing that pre-meiosis occurs in anthers of 1.5 mm inside floral buds (FBs), whereas meiosis between 1.5 and 4.2 mm FBs, and post-meiosis in FBs larger than 4.2 mm. These stages coincide with the Brazilian winter, a period of FBs reduced growth which suggests temperature sensitivity. Next, we identified and quantified the expression of reproductive 21- and 24-nt phasiRNAs during coffee anther development together with their canonical and novel miRNA triggers, and characterized the DCL and AGO families. Our results showed that the pattern of reproductive phasiRNA abundance in C. arabica is unique among described eudicots and the canonical trigger car-miR2275 is involved in the processing of both 21- and 24-nt phasiRNAs. Fourteen DCL genes were identified, but DCL5, related to phasiRNA biosynthesis in monocots, was not, according to its specificity for monocots. Thus, our work explored the knowledge gap about microsporogenesis and related siRNAs pathways in coffee, contributing to the control of reproductive development and the improvement of fertility in eudicots.

Giant viruses (Nucleocytoviricota) are significant lethality agents of various eukaryotic hosts. Although metagenomics indicates their ubiquitous distribution, available giant virus isolates are restricted to a very small number of protist and algal hosts. Here we report on the first viral isolate that replicates in the amoeboflagellate Naegleria. This genus comprises the notorious human pathogen Naegleria fowleri, the causative agent of the rare but fatal primary amoebic meningoencephalitis. We have elucidated the structure and infection cycle of this giant virus, Catovirus naegleriensis (a.k.a. Naegleriavirus, NiV), and show its unique adaptations to its Naegleria host using fluorescence in situ hybridization, electron microscopy, genomics, and proteomics. Naegleriavirus is only the fourth isolate of the highly diverse subfamily Klosneuvirinae, and like its relatives the NiV genome contains a large number of translation genes, but lacks transfer RNAs (tRNAs). NiV has acquired genes from its Naegleria host, which code for heat shock proteins and apoptosis inhibiting factors, presumably for host interactions. Notably, NiV infection was lethal to all Naegleria species tested, including the human pathogen N. fowleri. This study expands our experimental framework for investigating giant viruses and may help to better understand the basic biology of the human pathogen N. fowleri.

Bridging molecular information to ecosystem-level processes would provide the capacity to understand system vulnerability and, potentially, a means for assessing ecosystem health. Here, we present an integrated dataset containing environmental and metagenomic information from plant-associated microbial communities, plant transcriptomics, plant and soil metabolomics, and soil chemistry and activity characterization measurements derived from the model tree species Populus trichocarpa. Soil, rhizosphere, root endosphere, and leaf samples were collected from 27 different P. trichocarpa genotypes grown in two different environments leading to an integrated dataset of 318 metagenomes, 98 plant transcriptomes, and 314 metabolomic profiles that are supported by diverse soil measurements. This expansive dataset will provide insights into causal linkages that relate genomic features and molecular level events to system-level properties and their environmental influences.

Gene cloning in repeat-rich polyploid genomes remains challenging. Here, we describe a strategy for overcoming major bottlenecks in cloning of the powdery mildew resistance gene (R-gene) Pm69 derived from tetraploid wild emmer wheat. A conventional positional cloning approach was not effective owing to suppressed recombination. Chromosome sorting was compromised by insufficient purity. A Pm69 physical map, constructed by assembling Oxford Nanopore Technology (ONT) long-read genome sequences, revealed a rapidly evolving nucleotide-binding leucine-rich repeat (NLR) R-gene cluster with structural variations. A single candidate NLR was identified by anchoring RNA sequencing reads from susceptible mutants to ONT contigs and was validated by virus-induced gene silencing. Pm69 is likely a newly evolved NLR and was discovered in only one location across the wild emmer wheat distribution range in Israel. Pm69 was successfully introgressed into cultivated wheat, and a diagnostic molecular marker was used to accelerate its deployment and pyramiding with other R-genes.

Organic bilayers or amorphous silica films of a few nanometer thickness featuring embedded molecular wires offer opportunities for chemically separating while at the same time electronically connecting photo- or electrocatalytic components. Such ultrathin membranes enable the integration of components for which direct coupling is not sufficiently efficient or stable. Photoelectrocatalytic systems for the generation or utilization of renewable energy are among the most prominent ones for which ultrathin separation layers open up new approaches for component integration for improving efficiency. Recent advances in the assembly and spectroscopic, microscopic, and photoelectrochemical characterization have enabled the systematic optimization of the structure, energetics, and density of embedded molecular wires for maximum charge transfer efficiency. The progress enables interfacial designs for the nanoscale integration of the incompatible oxidation and reduction catalysis environments of artificial photosystems and of microbial (or biomolecular)-abiotic systems for renewable energy.

Intimate associations between fungi and intracellular bacterial endosymbionts are becoming increasingly well understood. Phylogenetic analyses demonstrate that bacterial endosymbionts of Mucoromycota fungi are related either to free-living Burkholderia or Mollicutes species. The so-called Burkholderia-related endosymbionts or BRE comprise Mycoavidus, Mycetohabitans and Candidatus Glomeribacter gigasporarum. These endosymbionts are marked by genome contraction thought to be associated with intracellular selection. However, the conclusions drawn thus far are based on a very small subset of endosymbiont genomes, and the mechanisms leading to genome streamlining are not well understood. The purpose of this study was to better understand how intracellular existence shapes Mycoavidus and BRE functionally at the genome level. To this end we generated and analyzed 14 novel draft genomes for Mycoavidus living within the hyphae of Mortierellomycotina fungi. We found that our novel Mycoavidus genomes were significantly reduced compared to free-living Burkholderiales relatives. Using a genome-scale phylogenetic approach including the novel and available existing genomes of Mycoavidus, we show that the genus is an assemblage composed of two independently derived lineages including three well supported clades of Mycoavidus. Using a comparative genomic approach, we shed light on the functional implications of genome reduction, documenting shared and unique gene loss patterns between the three Mycoavidus clades. We found that many endosymbiont isolates demonstrate patterns of vertical transmission and host-specificity, but others are present in phylogenetically disparate hosts. We discuss how reductive evolution and host specificity reflect convergent adaptation to the intrahyphal selective landscape, and commonalities of eukaryotic endosymbiont genome evolution.

Marine algae are central to global carbon fixation, and their productivity is dictated largely by resource availability. Reduced nutrient availability is predicted for vast oceanic regions as an outcome of climate change; however, there is much to learn regarding response mechanisms of the tiny picoplankton that thrive in these environments, especially eukaryotic phytoplankton. Here, we investigate responses of the picoeukaryote Micromonas commoda, a green alga found throughout subtropical and tropical oceans. Under shifting phosphate availability scenarios, transcriptomic analyses revealed altered expression of transfer RNA modification enzymes and biased codon usage of transcripts more abundant during phosphate-limiting versus phosphate-replete conditions, consistent with the role of transfer RNA modifications in regulating codon recognition. To associate the observed shift in the expression of the transfer RNA modification enzyme complement with the transfer RNAs encoded by M. commoda, we also determined the transfer RNA repertoire of this alga revealing potential targets of the modification enzymes. Codon usage bias was particularly pronounced in transcripts encoding proteins with direct roles in managing phosphate limitation and photosystem-associated proteins that have ill-characterized putative functions in "light stress." The observed codon usage bias corresponds to a proposed stress response mechanism in which the interplay between stress-induced changes in transfer RNA modifications and skewed codon usage in certain essential response genes drives preferential translation of the encoded proteins. Collectively, we expose a potential underlying mechanism for achieving growth under enhanced nutrient limitation that extends beyond the catalog of up- or downregulated protein-encoding genes to the cell biological controls that underpin acclimation to changing environmental conditions.

Fungal effectors play critical roles in manipulating plant immune responses and promoting colonization. Sphaerulina musiva is a heterothallic ascomycete fungus that causes Septoria leaf spot and stem canker disease in poplar (Populus spp.) plantations. This disease can result in premature defoliation, branch and stem breakage, increased mortality, and plantation failure. However, little is known about the interaction between S. musiva and poplar. Previous work predicted 142 candidate secreted effector proteins in S. musiva (SmCSEPs), 19 of which were selected for further functional characterization in this study. SmCSEP3 induced plant cell death in Nicotiana benthamiana, while 8 out of 19 tested SmCSEPs suppressed cell death. The signal peptides of these eight SmCSEPs exhibited secretory activity in a yeast signal sequence trap assay. Confocal microscopy revealed that four of these eight SmCSEPs target both the cytoplasm and the nucleus, whereas four predominantly localize to discrete punctate structures. Pathogen challenge assays in N. benthamiana demonstrated that the transient expression of six SmCSEPs promoted Fusarium proliferatum infection. The expression of these six SmCSEP genes were induced during infection. SmCSEP2, SmCSEP13, and SmCSEP25 suppressed chitin-triggered reactive oxygen species burst and callose deposition in N. benthamiana. The candidate secreted effector proteins of S. musiva target multiple compartments in the plant cell and modulate different pattern-triggered immunity pathways. [Formula: see text] The author(s) have dedicated the work to the public domain under the Creative Commons CC0 "No Rights Reserved" license by waiving all of his or her rights to the work worldwide under copyright law, including all related and neighboring rights, to the extent allowed by law, 2023.

Abnormal expression of the transcriptional regulator and hedgehog (Hh) signaling pathway effector Gli3 is known to trigger congenital disease, most frequently affecting the central nervous system (CNS) and the limbs. Accurate delineation of the genomic cis-regulatory landscape controlling Gli3 transcription during embryonic development is critical for the interpretation of noncoding variants associated with congenital defects. Here, we employed a comparative genomic analysis on fish species with a slow rate of molecular evolution to identify seven previously unknown conserved noncoding elements (CNEs) in Gli3 intronic intervals (CNE15-21). Transgenic assays in zebrafish revealed that most of these elements drive activities in Gli3 expressing tissues, predominantly the fins, CNS, and the heart. Intersection of these CNEs with human disease associated SNPs identified CNE15 as a putative mammalian craniofacial enhancer, with conserved activity in vertebrates and potentially affected by mutation associated with human craniofacial morphology. Finally, comparative functional dissection of an appendage-specific CNE conserved in slowly evolving fish (elephant shark), but not in teleost (CNE14/hs1586) indicates co-option of limb specificity from other tissues prior to the divergence of amniotes and lobe-finned fish. These results uncover a novel subset of intronic Gli3 enhancers that arose in the common ancestor of gnathostomes and whose sequence components were likely gradually modified in other species during the process of evolutionary diversification.

One of the reasons for the high efficiency and selectivity of biological catalysts arise from their ability to control the pathways of substrates and products using protein channels, and by modulating the transport in the channels using the interaction with the protein residues and the water/hydrogen-bonding network. This process is clearly demonstrated in Photosystem II (PS II), where its light-driven water oxidation reaction catalyzed by the Mn4CaO5 cluster occurs deep inside the protein complex and thus requires the transport of two water molecules to and four protons from the metal center to the bulk water. Based on the recent advances in structural studies of PS II from X-ray crystallography and cryo-electron microscopy, in this review we compare the channels that have been proposed to facilitate this mass transport in cyanobacteria, red and green algae, diatoms, and higher plants. The three major channels (O1, O4, and Cl1 channels) are present in all species investigated; however, some differences exist in the reported structures that arise from the different composition and arrangement of membrane extrinsic subunits between the species. Among the three channels, the Cl1 channel, including the proton gate, is the most conserved among all photosynthetic species. We also found at least one branch for the O1 channel in all organisms, extending all the way from Ca/O1 via the 'water wheel' to the lumen. However, the extending path after the water wheel varies between most species. The O4 channel is, like the Cl1 channel, highly conserved among all species while having different orientations at the end of the path near the bulk. The comparison suggests that the previously proposed functionality of the channels in T. vestitus (Ibrahim et al., Proc Natl Acad Sci USA 117:12624-12635, 2020; Hussein et al., Nat Commun 12:6531, 2021) is conserved through the species, i.e. the O1-like channel is used for substrate water intake, and the tighter Cl1 and O4 channels for proton release. The comparison does not eliminate the potential role of O4 channel as a water intake channel. However, the highly ordered hydrogen-bonded water wire connected to the Mn4CaO5 cluster via the O4 may strongly suggest that it functions in proton release, especially during the S0 → S1 transition (Saito et al., Nat Commun 6:8488, 2015; Kern et al., Nature 563:421-425, 2018; Ibrahim et al., Proc Natl Acad Sci USA 117:12624-12635, 2020; Sakashita et al., Phys Chem Chem Phys 22:15831-15841, 2020; Hussein et al., Nat Commun 12:6531, 2021).

Previous studies have shown that abscisic acid (ABA) and ethylene are involved in pulp maturation and peel coloration in the nonclimacteric citrus fruits. There are also signs indicating that other plant hormones may play some roles in citrus fruit ripening. In this study, we compared profiles of genome-wide gene expression and changes in hormones and peel pigments between fruits of Shatangju mandarin (Citrus reticulata Blanco, designated WT) and its natural mutant, Yuenongwanju (designated MT). The MT fruit matures ~2 months later than the WT fruit. Significant differences in fruit diameter, total soluble solids, titratable acid content, chlorophylls and carotenoids were detected between the fruits of the two genotypes at the sampled time points. Genome-wide transcriptome profiling showed that many genes involved in auxin and ABA metabolism and/or signaling pathways were differentially expressed between the MT and the WT fruits. Importantly, the expression of CrYUCCA8 was significantly lower and the expression of CrNCED5 was significantly higher in WT than in MT fruits at 230 and 250 DPA, respectively. In addition, the indole-3-acetic acid (IAA) level in the MT fruit was significantly higher than that in the WT counterpart, whereas a significantly lower level of ABA was detected in the mutant. Treatment of the WT fruit with exogenous IAA significantly delayed fruit maturation. Our results provide experimental evidence supporting the notion that auxin is a negative regulator of fruit maturation in citrus.

The functions of closely related Myb-like repressor and Myb-like activator proteins within the plant circadian oscillator have been well-studied as separate groups, but the genetic interactions between them are less clear. We hypothesized that these repressors and activators would interact additively to regulate both circadian and growth phenotypes. We used CRISPR-Cas9 to generate new mutant alleles and performed physiological and molecular characterization of plant mutants for five of these core Myb-like clock factors compared with a repressor mutant and an activator mutant. We first examined circadian clock function in plants likely null for both the repressor proteins, CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) and LATE ELONGATED HYPOCOTYL (LHY), and the activator proteins, REVEILLE 4 (RVE4), REVEILLE (RVE6), and REVEILLE (RVE8). The rve468 triple mutant has a long period and flowers late, while cca1 lhy rve468 quintuple mutants, similarly to cca1 lhy mutants, have poor circadian rhythms and flower early. This suggests that CCA1 and LHY are epistatic to RVE4, RVE6, and RVE8 for circadian clock and flowering time function. We next examined hypocotyl elongation and rosette leaf size in these mutants. The cca1 lhy rve468 mutants have growth phenotypes intermediate between cca1 lhy and rve468 mutants, suggesting that CCA1, LHY, RVE4, RVE6, and RVE8 interact additively to regulate growth. Together, our data suggest that these five Myb-like factors interact differently in regulation of the circadian clock versus growth. More generally, the near-norm al seedling phenotypes observed in the largely arrhythmic quintuple mutant demonstrate that circadian-regulated output processes, like control of hypocotyl elongation, do not always depend upon rhythmic oscillator function.

Molecular clocks are the basis for dating the divergence between lineages over macroevolutionary timescales (~105 to 108 years). However, classical DNA-based clocks tick too slowly to inform us about the recent past. Here, we demonstrate that stochastic DNA methylation changes at a subset of cytosines in plant genomes display a clocklike behavior. This "epimutation clock" is orders of magnitude faster than DNA-based clocks and enables phylogenetic explorations on a scale of years to centuries. We show experimentally that epimutation clocks recapitulate known topologies and branching times of intraspecies phylogenetic trees in the self-fertilizing plant Arabidopsis thaliana and the clonal seagrass Zostera marina, which represent two major modes of plant reproduction. This discovery will open new possibilities for high-resolution temporal studies of plant biodiversity.

Species of the phylum Blastocladiomycota, early-diverging zoosporic (flagellated) lineages of fungi, are vastly understudied. This phylum includes the genus Coelomomyces, which consists of more than 80 fungal species that are obligate parasites of arthropods. Known Coelomomyces species lack a complete asexual life cycle, instead surviving through an obligate heterecious alternation of generations life cycle. Despite their global distribution and interesting life cycle, little is known about the genomics of any Coelomomyces species. To address this, we generated three draft-level genomes and annotations for C. lativittatus representing its haploid meiospore, orange gamete, and amber gamete life stages. These draft genome assemblies ranged in size from 5002 to 5799 contigs, with a total length of 19.8-22.8 Mb and a mean of 7416 protein-coding genes. We then demonstrated the utility of these genomes by combining the draft annotations as a reference for analysis of C. lativittatus transcriptomes. We analyzed transcriptomes from across host-associated life stages, including infected larvae and excised mature sporangia from the mosquito Anopheles quadrimaculatus. We identified differentially expressed genes and enriched GO terms both across and within life stages and used these to make hypotheses about C. lativittatus biology. Generally, we found the C. lativittatus transcriptome to be a complex and dynamic expression landscape; GO terms related to metabolism and transport processes were enriched during infection and terms related to dispersal were enriched during sporulation. We further identified five high mobility group (HMG)-box genes in C. lativittatus, three belonging to clades with mating type (MAT) loci from other fungi, as well as four ortholog expansions in C. lativittatus compared with other fungi. The C. lativittatus genomes and transcriptomes reported here are a valuable resource and may be leveraged toward furthering understanding of the biology of these and other early-diverging fungal lineages.

Spray-induced gene silencing (SIGS) is an emerging tool for crop pest protection. It utilizes exogenously applied double-stranded RNA to specifically reduce pest target gene expression using endogenous RNA interference machinery. In this study, SIGS methods were developed and optimized for powdery mildew fungi, which are widespread obligate biotrophic fungi that infect agricultural crops, using the known azole-fungicide target cytochrome P450 51 (CYP51) in the Golovinomyces orontii-Arabidopsis thaliana pathosystem. Additional screening resulted in the identification of conserved gene targets and processes important to powdery mildew proliferation: apoptosis-antagonizing transcription factor in essential cellular metabolism and stress response; lipid catabolism genes lipase a, lipase 1, and acetyl-CoA oxidase in energy production; and genes involved in manipulation of the plant host via abscisic acid metabolism (9-cis-epoxycarotenoid dioxygenase, xanthoxin dehydrogenase, and a putative abscisic acid G-protein coupled receptor) and secretion of the effector protein, effector candidate 2. Powdery mildew is the dominant disease impacting grapes and extensive powdery mildew resistance to applied fungicides has been reported. We therefore developed SIGS for the Erysiphe necator-Vitis vinifera system and tested six successful targets identified using the G. orontii-A. thaliana system. For all targets tested, a similar reduction in powdery mildew disease was observed between systems. This indicates screening of broadly conserved targets in the G. orontii-A. thaliana pathosystem identifies targets and processes for the successful control of other powdery mildew fungi. The efficacy of SIGS on powdery mildew fungi makes SIGS an exciting prospect for commercial powdery mildew control.

The DELLA genes, also known as 'Green Revolution' genes, encode conserved master growth regulators that control plant development in response to internal and environmental cues. Functioning as nuclear-localized transcription regulators, DELLAs modulate expression of target genes via direct protein-protein interaction of their carboxy-terminal GRAS domain with hundreds of transcription factors (TFs) and epigenetic regulators. However, the molecular mechanism of DELLA-mediated transcription reprogramming remains unclear. Here by characterizing new missense alleles of an Arabidopsis DELLA, repressor of ga1-3 (RGA), and co-immunoprecipitation assays, we show that RGA binds histone H2A via the PFYRE subdomain within its GRAS domain to form a TF-RGA-H2A complex at the target chromatin. Chromatin immunoprecipitation followed by sequencing analysis further shows that this activity is essential for RGA association with its target chromatin globally. Our results indicate that, although DELLAs are recruited to target promoters by binding to TFs via the LHR1 subdomain, DELLA-H2A interaction via the PFYRE subdomain is necessary to stabilize the TF-DELLA-H2A complex at the target chromatin. This study provides insights into the two distinct key modular functions in DELLA for its genome-wide transcription regulation in plants.

Despite their global presence and ubiquity, members of the class Geoglossomycetes (Pezizomycotina, Ascomycota) are understudied systematically and ecologically. These fungi have long been presumed saprobic due to their occurrence in or near leaf litter and soils. Additionally, they lack an apparent association with other organisms, reinforcing this perception. However, observations of sporocarps near ericaceous shrubs have given rise to an alternative hypothesis that members of Geoglossomycetes may form ericoid mycorrhizae or ectomycorrhizae. This claim, however, has yet to be confirmed via microscopy or amplicon-based studies examining root communities. As a result, our current understanding of their ecology is based on cursory observations. This study presents a comparative analysis of genomic signatures related to ecological niche to investigate the hypothesis of an ericoid mycorrhizal or ectomycorrhizal ecology in the class. We compared the carbohydrate-active enzyme (CAZyme) and secondary metabolite contents of six newly sequenced Geoglossomycetes genomes with those of fungi representing specific ecologies across Pezizomycotina. Our analysis reveals CAZyme and secondary metabolite content patterns consistent with ectomycorrhizal (EcM) members of Pezizomycotina. Specifically, we found a reduction in CAZyme-encoding genes and secondary metabolite clusters that suggests a mutualistic ecology. Our work includes the broadest taxon sampling for a phylogenomic study of Pezizomycotina to date. It represents the first functional genomic and genome-scale phylogenetic study of the class Geoglossomycetes and improves the foundational knowledge of the ecology and evolution of these understudied fungi.

The morphological diversity of the inflorescence determines flower and seed production, which is critical for plant adaptation. Hall's panicgrass (Panicum hallii, P. hallii) is a wild perennial grass that has been developed as a model to study perennial grass biology and adaptive evolution. Highly divergent inflorescences have evolved between the 2 major ecotypes in P. hallii, the upland ecotype (P. hallii var hallii, HAL2 genotype) with compact inflorescence and large seed and the lowland ecotype (P. hallii var filipes, FIL2 genotype) with an open inflorescence and small seed. Here we conducted a comparative analysis of the transcriptome and DNA methylome, an epigenetic mark that influences gene expression regulation, across different stages of inflorescence development using genomic references for each ecotype. Global transcriptome analysis of differentially expressed genes (DEGs) and co-expression modules underlying the inflorescence divergence revealed the potential role of cytokinin signaling in heterochronic changes. Comparing DNA methylome profiles revealed a remarkable level of differential DNA methylation associated with the evolution of P. hallii inflorescence. We found that a large proportion of differentially methylated regions (DMRs) were located in the flanking regulatory regions of genes. Intriguingly, we observed a substantial bias of CHH hypermethylation in the promoters of FIL2 genes. The integration of DEGs, DMRs, and Ka/Ks ratio results characterized the evolutionary features of DMR-associated DEGs that contribute to the divergence of the P. hallii inflorescence. This study provides insights into the transcriptome and epigenetic landscape of inflorescence divergence in P. hallii and a genomic resource for perennial grass biology.