A fine-scale Arabidopsis chromatin landscape reveals chromatin conformation-associated transcriptional dynamics.

Plants, as sessile organisms, deploy transcriptional dynamics for adapting to extreme growth conditions such as cold stress. Emerging evidence suggests that chromatin architecture contributes to transcriptional regulation. However, the relationship between chromatin architectural dynamics and transcriptional reprogramming in response to cold stress remains unclear. Here, we apply a chemical-crosslinking assisted proximity capture (CAP-C) method to elucidate the fine-scale chromatin landscape, revealing chromatin interactions within gene bodies closely associated with RNA polymerase II (Pol II) densities across initiation, pausing, and termination sites. We observe dynamic changes in chromatin interactions alongside Pol II activity alterations during cold stress, suggesting local chromatin dynamics may regulate Pol II activity. Notably, cold stress does not affect large-scale chromatin conformations. We further identify a comprehensive promoter-promoter interaction (PPI) network across the genome, potentially facilitating co-regulation of gene expression in response to cold stress. Our study deepens the understanding of chromatin conformation-associated gene regulation in plant response to cold.

Distinct accessory roles of Arabidopsis VEL proteins in Polycomb silencing.

Polycomb repressive complex 2 (PRC2) mediates epigenetic silencing of target genes in animals and plants. In Arabidopsis, PRC2 is required for the cold-induced epigenetic silencing of the FLC floral repressor locus to align flowering with spring. During this process, PRC2 relies on VEL accessory factors, including the constitutively expressed VRN5 and the cold-induced VIN3. The VEL proteins are physically associated with PRC2, but their individual functions remain unclear. Here, we show an intimate association between recombinant VRN5 and multiple components within a reconstituted PRC2, dependent on a compact conformation of VRN5 central domains. Key residues mediating this compact conformation are conserved among VRN5 orthologs across the plant kingdom. In contrast, VIN3 interacts with VAL1, a transcriptional repressor that binds directly to FLC These associations differentially affect their role in H3K27me deposition: Both proteins are required for H3K27me3, but only VRN5 is necessary for H3K27me2. Although originally defined as vernalization regulators, VIN3 and VRN5 coassociate with many targets in the Arabidopsis genome that are modified with H3K27me3. Our work therefore reveals the distinct accessory roles for VEL proteins in conferring cold-induced silencing on FLC, with broad relevance for PRC2 targets generally.

A Possible Aquatic Origin of the Thiaminase TenA of the Human Gut Symbiont Bacteroides thetaiotaomicron.

TenA thiamin-degrading enzymes are commonly found in prokaryotes, plants, fungi and algae and are involved in the thiamin salvage pathway. The gut symbiont Bacteroides thetaiotaomicron (Bt) produces a TenA protein (BtTenA) which is packaged into its extracellular vesicles. An alignment of BtTenA protein sequence with proteins from different databases using the basic local alignment search tool (BLAST) and the generation of a phylogenetic tree revealed that BtTenA is related to TenA-like proteins not only found in a small number of intestinal bacterial species but also in some aquatic bacteria, aquatic invertebrates, and freshwater fish. This is, to our knowledge, the first report describing the presence of TenA-encoding genes in the genome of members of the animal kingdom. By searching metagenomic databases of diverse host-associated microbial communities, we found that BtTenA homologues were mostly represented in biofilms present on the surface of macroalgae found in Australian coral reefs. We also confirmed the ability of a recombinant BtTenA to degrade thiamin. Our study shows that BttenA-like genes which encode a novel sub-class of TenA proteins are sparingly distributed across two kingdoms of life, a feature of accessory genes known for their ability to spread between species through horizontal gene transfer.

Genomics and biochemical analyses reveal a metabolon key to ß-L-ODAP biosynthesis in Lathyrus sativus.

Grass pea (Lathyrus sativus L.) is a rich source of protein cultivated as an insurance crop in Ethiopia, Eritrea, India, Bangladesh, and Nepal. Its resilience to both drought and flooding makes it a promising crop for ensuring food security in a changing climate. The lack of genetic resources and the crop's association with the disease neurolathyrism have limited the cultivation of grass pea. Here, we present an annotated, long read-based assembly of the 6.5 Gbp L. sativus genome. Using this genome sequence, we have elucidated the biosynthetic pathway leading to the formation of the neurotoxin, ß-L-oxalyl-2,3-diaminopropionic acid (ß-L-ODAP). The final reaction of the pathway depends on an interaction between L. sativus acyl-activating enzyme 3 (LsAAE3) and a BAHD-acyltransferase (LsBOS) that form a metabolon activated by CoA to produce ß-L-ODAP. This provides valuable insight into the best approaches for developing varieties which produce substantially less toxin.

RNA G-quadruplex structure contributes to cold adaptation in plants.

Nucleotide composition is suggested to infer gene functionality and ecological adaptation of species to distinct environments. However, the underlying biological function of nucleotide composition dictating environmental adaptations is largely unknown. Here, we systematically analyze the nucleotide composition of transcriptomes across 1000 plants (1KP) and their corresponding habitats. Intriguingly, we find that plants growing in cold climates have guanine (G)-enriched transcriptomes, which are prone to forming RNA G-quadruplex structures. Both immunofluorescence detection and in vivo structure profiling reveal that RNA G-quadruplex formation in plants is globally enhanced in response to cold. Cold-responsive RNA G-quadruplexes strongly enhanced mRNA stability, rather than affecting translation. Disruption of individual RNA G-quadruplex promotes mRNA decay in the cold, leading to impaired plant cold response. Therefore, we propose that plants adopted RNA G-quadruplex structure as a molecular signature to facilitate their adaptation to the cold during evolution.

In vivo single-molecule analysis reveals COOLAIR RNA structural diversity.

Cellular RNAs are heterogeneous with respect to their alternative processing and secondary structures, but the functional importance of this complexity is still poorly understood. A set of alternatively processed antisense non-coding transcripts, which are collectively called COOLAIR, are generated at the Arabidopsis floral-repressor locus FLOWERING LOCUS C (FLC)1. Different isoforms of COOLAIR influence FLC transcriptional output in warm and cold conditions2-7. Here, to further investigate the function of COOLAIR, we developed an RNA structure-profiling method to determine the in vivo structure of single RNA molecules rather than the RNA population average. This revealed that individual isoforms of the COOLAIR transcript adopt multiple structures with different conformational dynamics. The major distally polyadenylated COOLAIR isoform in warm conditions adopts three predominant structural conformations, the proportions and conformations of which change after cold exposure. An alternatively spliced, strongly cold-upregulated distal COOLAIR isoform6 shows high structural diversity, in contrast to proximally polyadenylated COOLAIR. A hyper-variable COOLAIR structural element was identified that was complementary to the FLC transcription start site. Mutations altering the structure of this region changed FLC expression and flowering time, consistent with an important regulatory role of the COOLAIR structure in FLC transcription. Our work demonstrates that isoforms of non-coding RNA transcripts adopt multiple distinct and functionally relevant structural conformations, which change in abundance and shape in response to external conditions.

An Integrated Linkage Map of Three Recombinant Inbred Populations of Pea (Pisum sativum L.).

Biparental recombinant inbred line (RIL) populations are sets of genetically stable lines and have a simple population structure that facilitates the dissection of the genetics of interesting traits. On the other hand, populations derived from multiparent intercrosses combine both greater diversity and higher numbers of recombination events than RILs. Here, we describe a simple population structure: a three-way recombinant inbred population combination. This structure was easy to produce and was a compromise between biparental and multiparent populations. We show that this structure had advantages when analyzing cultivar crosses, and could achieve a mapping resolution of a few genes.

Population genomic analysis of Aegilops tauschii identifies targets for bread wheat improvement.

Aegilops tauschii, the diploid wild progenitor of the D subgenome of bread wheat, is a reservoir of genetic diversity for improving bread wheat performance and environmental resilience. Here we sequenced 242 Ae. tauschii accessions and compared them to the wheat D subgenome to characterize genomic diversity. We found that a rare lineage of Ae. tauschii geographically restricted to present-day Georgia contributed to the wheat D subgenome in the independent hybridizations that gave rise to modern bread wheat. Through k-mer-based association mapping, we identified discrete genomic regions with candidate genes for disease and pest resistance and demonstrated their functional transfer into wheat by transgenesis and wide crossing, including the generation of a library of hexaploids incorporating diverse Ae. tauschii genomes. Exploiting the genomic diversity of the Ae. tauschii ancestral diploid genome permits rapid trait discovery and functional genetic validation in a hexaploid background amenable to breeding.

Wheat in vivo RNA structure landscape reveals a prevalent role of RNA structure in modulating translational subgenome expression asymmetry.

BACKGROUND: Polyploidy, especially allopolyploidy, which entails merging divergent genomes via hybridization and whole-genome duplication (WGD), is a major route to speciation in plants. The duplication among the parental genomes (subgenomes) often leads to one subgenome becoming dominant over the other(s), resulting in subgenome asymmetry in gene content and expression. Polyploid wheats are allopolyploids with most genes present in two (tetraploid) or three (hexaploid) functional copies, which commonly show subgenome expression asymmetry. It is unknown whether a similar subgenome asymmetry exists during translation. We aim to address this key biological question and explore the major contributing factors to subgenome translation asymmetry. RESULTS: Here, we obtain the first tetraploid wheat translatome and reveal that subgenome expression asymmetry exists at the translational level. We further perform in vivo RNA structure profiling to obtain the wheat RNA structure landscape and find that mRNA structure has a strong impact on translation, independent of GC content. We discover a previously uncharacterized contribution of RNA structure in subgenome translation asymmetry. We identify 3564 single-nucleotide variations (SNVs) across the transcriptomes between the two tetraploid wheat subgenomes, which induce large RNA structure disparities. These SNVs are highly conserved within durum wheat cultivars but are divergent in both domesticated and wild emmer wheat. CONCLUSIONS: We successfully determine both the translatome and in vivo RNA structurome in tetraploid wheat. We reveal that RNA structure serves as an important modulator of translational subgenome expression asymmetry in polyploids. Our work provides a new perspective for molecular breeding of major polyploid crops.

Pan-genome analysis identifies intersecting roles for Pseudomonas specialized metabolites in potato pathogen inhibition.

Agricultural soil harbors a diverse microbiome that can form beneficial relationships with plants, including the inhibition of plant pathogens. Pseudomonas spp. are one of the most abundant bacterial genera in the soil and rhizosphere and play important roles in promoting plant health. However, the genetic determinants of this beneficial activity are only partially understood. Here, we genetically and phenotypically characterize the Pseudomonas fluorescens population in a commercial potato field, where we identify strong correlations between specialized metabolite biosynthesis and antagonism of the potato pathogens Streptomyces scabies and Phytophthora infestans. Genetic and chemical analyses identified hydrogen cyanide and cyclic lipopeptides as key specialized metabolites associated with S. scabies inhibition, which was supported by in planta biocontrol experiments. We show that a single potato field contains a hugely diverse and dynamic population of Pseudomonas bacteria, whose capacity to produce specialized metabolites is shaped both by plant colonization and defined environmental inputs.

Potato scab and blight are two major diseases which can cause heavy crop losses. They are caused, respectively, by the bacterium Streptomyces scabies and an oomycete (a fungus-like organism) known as Phytophthora infestans. Fighting these disease-causing microorganisms can involve crop management techniques ­ for example, ensuring that a field is well irrigated helps to keep S. scabies at bay. Harnessing biological control agents can also offer ways to control disease while respecting the environment. Biocontrol bacteria, such as Pseudomonas, can produce compounds that keep S. scabies and P. infestans in check. However, the identity of these molecules and how irrigation can influence Pseudomonas population remains unknown. To examine these questions, Pacheco-Moreno et al. sampled and isolated hundreds of Pseudomonas strains from a commercial potato field, closely examining the genomes of 69 of these. Comparing the genetic information of strains based on whether they could control the growth of S. scabies revealed that compounds known as cyclic lipopeptides are key to controlling the growth of S. scabies and P. infestans. Whether the field was irrigated also had a large impact on the strains forming the Pseudomonas population. Working out how Pseudomonas bacteria block disease could speed up the search for biological control agents. The approach developed by Pacheco-Moreno et al. could help to predict which strains might be most effective based on their genetic features. Similar experiments could also work for other combinations of plants and diseases.

Appressorium-mediated plant infection by Magnaporthe oryzae is regulated by a Pmk1-dependent hierarchical transcriptional network.

Rice blast is a devastating disease caused by the fungal pathogen Magnaporthe oryzae that threatens rice production around the world. The fungus produces a specialized infection cell, called the appressorium, that enables penetration through the plant cell wall in response to surface signals from the rice leaf. The underlying biology of plant infection, including the regulation of appressorium formation, is not completely understood. Here we report the identification of a network of temporally coregulated transcription factors that act downstream of the Pmk1 mitogen-activated protein kinase pathway to regulate gene expression during appressorium-mediated plant infection. We show that this tiered regulatory mechanism involves Pmk1-dependent phosphorylation of the Hox7 homeobox transcription factor, which regulates genes associated with induction of major physiological changes required for appressorium development-including cell-cycle control, autophagic cell death, turgor generation and melanin biosynthesis-as well as controlling a additional set of virulence-associated transcription factor-encoding genes. Pmk1-dependent phosphorylation of Mst12 then regulates gene functions involved in septin-dependent cytoskeletal re-organization, polarized exocytosis and effector gene expression, which are necessary for plant tissue invasion. Identification of this regulatory cascade provides new potential targets for disease intervention.

Draft Genome Sequence Resource for the Orange Rust Pathogen of Sugarcane Puccinia kuehnii.

Puccinia kuehnii is an obligate biotrophic fungal pathogen that causes orange rust of sugarcane, which is prevalent in many countries around the globe. In the United States, orange rust was first detected in sugarcane in Florida in 2007 and poses a persistent and economically damaging threat to the sugarcane industry in this region. Here, we generated the first genome assemblies for two isolates of P. kuehnii (1040 and 2143) collected in Florida in 2017 from two sugarcane cultivars, CL85-1040 and CP89-2143, respectively. These two rust genome resources will be of immense value for future genomic studies, particularly further exploration of the predicted secretomes that may help define key pathogenicity determinants for this economically important pathogen.

In vivo nuclear RNA structurome reveals RNA-structure regulation of mRNA processing in plants.

BACKGROUND: mRNA processing is critical for gene expression. A challenge in regulating mRNA processing is how to recognize the actual mRNA processing sites, such as splice and polyadenylation sites, when the sequence content is insufficient for this purpose. Previous studies suggested that RNA structure affects mRNA processing. However, the regulatory role of RNA structure in mRNA processing remains unclear. RESULTS: Here, we perform in vivo selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemical profiling on Arabidopsis and generate the in vivo nuclear RNA structure landscape. We find that nuclear mRNAs fold differently from cytosolic mRNAs across translation start and stop sites. Notably, we discover a two-nucleotide single-stranded RNA structure feature upstream of 5' splice sites that is strongly associated with splicing and the selection of alternative 5' splice sites. The regulatory role of this RNA structure feature is further confirmed by experimental validation. Moreover, we find the single-strandedness of branch sites is also associated with 3' splice site recognition. We also identify an RNA structure feature comprising two close-by single-stranded regions that is specifically associated with both polyadenylation and alternative polyadenylation events. CONCLUSIONS: We successfully identify pre-mRNA structure features associated with splicing and polyadenylation at whole-genome scale and validate an RNA structure feature which can regulate splicing. Our study unveils a new RNA structure regulatory mechanism for mRNA processing.

Formation and diversification of a paradigm biosynthetic gene cluster in plants.

Numerous examples of biosynthetic gene clusters (BGCs), including for compounds of agricultural and medicinal importance, have now been discovered in plant genomes. However, little is known about how these complex traits are assembled and diversified. Here, we examine a large number of variants within and between species for a paradigm BGC (the thalianol cluster), which has evolved recently in a common ancestor of the Arabidopsis genus. Comparisons at the species level reveal differences in BGC organization and involvement of auxiliary genes, resulting in production of species-specific triterpenes. Within species, the thalianol cluster is primarily fixed, showing a low frequency of deleterious haplotypes. We further identify chromosomal inversion as a molecular mechanism that may shuffle more distant genes into the cluster, so enabling cluster compaction. Antagonistic natural selection pressures are likely involved in shaping the occurrence and maintenance of this BGC. Our work sheds light on the birth, life and death of complex genetic and metabolic traits in plants.

RNA G-quadruplex structures exist and function in vivo in plants.

BACKGROUND: Guanine-rich sequences are able to form complex RNA structures termed RNA G-quadruplexes in vitro. Because of their high stability, RNA G-quadruplexes are proposed to exist in vivo and are suggested to be associated with important biological relevance. However, there is a lack of direct evidence for RNA G-quadruplex formation in living eukaryotic cells. Therefore, it is unclear whether any purported functions are associated with the specific sequence content or the formation of an RNA G-quadruplex structure. RESULTS: Using rG4-seq, we profile the landscape of those guanine-rich regions with the in vitro folding potential in the Arabidopsis transcriptome. We find a global enrichment of RNA G-quadruplexes with two G-quartets whereby the folding potential is strongly influenced by RNA secondary structures. Using in vitro and in vivo RNA chemical structure profiling, we determine that hundreds of RNA G-quadruplex structures are strongly folded in both Arabidopsis and rice, providing direct evidence of RNA G-quadruplex formation in living eukaryotic cells. Subsequent genetic and biochemical analyses show that RNA G-quadruplex folding is able to regulate translation and modulate plant growth. CONCLUSIONS: Our study reveals the existence of RNA G-quadruplex in vivo and indicates that RNA G-quadruplex structures act as important regulators of plant development and growth.

Intact RNA structurome reveals mRNA structure-mediated regulation of miRNA cleavage in vivo.

MicroRNA (miRNA)-mediated cleavage is involved in numerous essential cellular pathways. miRNAs recognize target RNAs via sequence complementarity. In addition to complementarity, in vitro and in silico studies have suggested that RNA structure may influence the accessibility of mRNAs to miRNA-induced silencing complexes (miRISCs), thereby affecting RNA silencing. However, the regulatory mechanism of mRNA structure in miRNA cleavage remains elusive. We investigated the role of in vivo RNA secondary structure in miRNA cleavage by developing the new CAP-STRUCTURE-seq method to capture the intact mRNA structurome in Arabidopsis thaliana. This approach revealed that miRNA target sites were not structurally accessible for miRISC binding prior to cleavage in vivo. Instead, we found that the unfolding of the target site structure plays a key role in miRISC activity in vivo. We found that the single-strandedness of the two nucleotides immediately downstream of the target site, named Target Adjacent nucleotide Motif, can promote miRNA cleavage but not miRNA binding, thus decoupling target site binding from cleavage. Our findings demonstrate that mRNA structure in vivo can modulate miRNA cleavage, providing evidence of mRNA structure-dependent regulation of biological processes.

Fulvic acid increases forage legume growth inducing preferential up-regulation of nodulation and signalling-related genes.

The use of potential biostimulants is of broad interest in plant science for improving yields. The application of a humic derivative called fulvic acid (FA) may improve forage crop production. FA is an uncharacterized mixture of chemicals and, although it has been reported to increase growth parameters in many species including legumes, its mode of action remains unclear. Previous studies of the action of FA have lacked appropriate controls, and few have included field trials. Here we report yield increases due to FA application in three European Medicago sativa cultivars, in studies which include the appropriate nutritional controls which hitherto have not been used. No significant growth stimulation was seen after FA treatment in grass species in this study at the treatment rate tested. Direct application to bacteria increased Rhizobium growth and, in M. sativa trials, root nodulation was stimulated. RNA transcriptional analysis of FA-treated plants revealed up-regulation of many important early nodulation signalling genes after only 3 d. Experiments in plate, glasshouse, and field environments showed yield increases, providing substantial evidence for the use of FA to benefit M. sativa forage production.

Sparse Epistatic Patterns in the Evolution of Terpene Synthases.

We explore sequence determinants of enzyme activity and specificity in a major enzyme family of terpene synthases. Most enzymes in this family catalyze reactions that produce cyclic terpenes-complex hydrocarbons widely used by plants and insects in diverse biological processes such as defense, communication, and symbiosis. To analyze the molecular mechanisms of emergence of terpene cyclization, we have carried out in-depth examination of mutational space around (E)-ß-farnesene synthase, an Artemisia annua enzyme which catalyzes production of a linear hydrocarbon chain. Each mutant enzyme in our synthetic libraries was characterized biochemically, and the resulting reaction rate data were used as input to the Michaelis-Menten model of enzyme kinetics, in which free energies were represented as sums of one-amino-acid contributions and two-amino-acid couplings. Our model predicts measured reaction rates with high accuracy and yields free energy landscapes characterized by relatively few coupling terms. As a result, the Michaelis-Menten free energy landscapes have simple, interpretable structure and exhibit little epistasis. We have also developed biophysical fitness models based on the assumption that highly fit enzymes have evolved to maximize the output of correct products, such as cyclic products or a specific product of interest, while minimizing the output of byproducts. This approach results in nonlinear fitness landscapes that are considerably more epistatic. Overall, our experimental and computational framework provides focused characterization of evolutionary emergence of novel enzymatic functions in the context of microevolutionary exploration of sequence space around naturally occurring enzymes.

Agricultural Selection of Wheat Has Been Shaped by Plant-Microbe Interactions.

The influence of wheat (modern wheat, both bread and pasta, their wild ancestors and synthetic hybrids) on the microbiota of their roots and surrounding soil is characterized. We isolated lines of bread wheat by hybridizing diploid (Aegilops tauschii) with tetraploid Triticum durum and crossed it with a modern cultivar of Triticum aestivum. The newly created, synthetic hybrid wheat, which recapitulate the breeding history of wheat through artificial selection, is found to support a microbiome enriched in beneficial Glomeromycetes fungi, but also in, potentially detrimental, Nematoda. We hypothesize that during wheat domestication this plant-microbe interaction diminished, suggesting an evolutionary tradeoff; sacrificing advantageous nutrient acquisition through fungal interactions to minimize interaction with pathogenic fungi. Increased plant selection for Glomeromycetes and Nematoda is correlated with the D genome derived from A. tauschii. Despite differences in their soil microbiota communities, overall wheat plants consistently show a low ratio of eukaryotes to prokaryotes. We propose that this is a mechanism for protection against soil-borne fungal disease and appears to be deeply rooted in the wheat genome. We suggest that the influence of plants on the composition of their associated microbiota is an integral factor, hitherto overlooked, but intrinsic to selection during wheat domestication.