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.

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.

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.

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.

A noncanonical vacuolar sugar transferase required for biosynthesis of antimicrobial defense compounds in oat.

Plants produce an array of natural products with important ecological functions. These compounds are often decorated with oligosaccharide groups that influence bioactivity, but the biosynthesis of such sugar chains is not well understood. Triterpene glycosides (saponins) are a large family of plant natural products that determine important agronomic traits, as exemplified by avenacins, antimicrobial defense compounds produced by oats. Avenacins have a branched trisaccharide moiety consisting of l-arabinose linked to 2 d-glucose molecules that is critical for antifungal activity. Plant natural product glycosylation is usually performed by uridine diphosphate-dependent glycosyltransferases (UGTs). We previously characterized the arabinosyltransferase that initiates the avenacin sugar chain; however, the enzymes that add the 2 remaining d-glucose molecules have remained elusive. Here we characterize the enzymes that catalyze these last 2 glucosylation steps. AsUGT91G16 is a classical cytosolic UGT that adds a 1,2-linked d-glucose molecule to l-arabinose. Unexpectedly, the enzyme that adds the final 1,4-linked d-glucose (AsTG1) is not a UGT, but rather a sugar transferase belonging to Glycosyl Hydrolase family 1 (GH1). Unlike classical UGTs, AsTG1 is vacuolar. Analysis of oat mutants reveals that AsTG1 corresponds to Sad3, a previously uncharacterized locus shown by mutation to be required for avenacin biosynthesis. AsTG1 and AsUGT91G16 form part of the avenacin biosynthetic gene cluster. Our demonstration that a vacuolar transglucosidase family member plays a critical role in triterpene biosynthesis highlights the importance of considering other classes of carbohydrate-active enzymes in addition to UGTs as candidates when elucidating pathways for the biosynthesis of glycosylated natural products in plants.

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.

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.

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.

FoldAtlas: a repository for genome-wide RNA structure probing data.

Most RNA molecules form internal base pairs, leading to a folded secondary structure. Some of these structures have been demonstrated to be functionally significant. High-throughput RNA structure chemical probing methods generate millions of sequencing reads to provide structural constraints for RNA secondary structure prediction. At present, processed data from these experiments are difficult to access without computational expertise. Here we present FoldAtlas, a web interface for accessing raw and processed structural data across thousands of transcripts. FoldAtlas allows a researcher to easily locate, view, and retrieve probing data for a given RNA molecule. We also provide in silico and in vivo secondary structure predictions for comparison, visualized in the browser as circle plots and topology diagrams. Data currently integrated into FoldAtlas are from a new high-depth Structure-seq data analysis in Arabidopsis thaliana, released with this work. AVAILABILITY AND IMPLEMENTATION: The FoldAtlas website can be accessed at www.foldatlas.com Source code is freely available at github.com/mnori/foldatlas under the MIT license. Raw reads data are available under the NCBI SRA accession SRP066985. CONTACT: yiliang.ding@jic.ac.uk or matthew.norris@jic.ac.ukSupplementary information: Supplementary data are available at Bioinformatics online.

Algal ancestor of land plants was preadapted for symbiosis.

Colonization of land by plants was a major transition on Earth, but the developmental and genetic innovations required for this transition remain unknown. Physiological studies and the fossil record strongly suggest that the ability of the first land plants to form symbiotic associations with beneficial fungi was one of these critical innovations. In angiosperms, genes required for the perception and transduction of diffusible fungal signals for root colonization and for nutrient exchange have been characterized. However, the origin of these genes and their potential correlation with land colonization remain elusive. A comprehensive phylogenetic analysis of 259 transcriptomes and 10 green algal and basal land plant genomes, coupled with the characterization of the evolutionary path leading to the appearance of a key regulator, a calcium- and calmodulin-dependent protein kinase, showed that the symbiotic signaling pathway predated the first land plants. In contrast, downstream genes required for root colonization and their specific expression pattern probably appeared subsequent to the colonization of land. We conclude that the most recent common ancestor of extant land plants and green algae was preadapted for symbiotic associations. Subsequent improvement of this precursor stage in early land plants through rounds of gene duplication led to the acquisition of additional pathways and the ability to form a fully functional arbuscular mycorrhizal symbiosis.

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.

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.

THREaD Mapper Studio: a novel, visual web server for the estimation of genetic linkage maps.

The estimation of genetic linkage maps is a key component in plant and animal research, providing both an indication of the genetic structure of an organism and a mechanism for identifying candidate genes associated with traits of interest. Because of this importance, several computational solutions to genetic map estimation exist, mostly implemented as stand-alone software packages. However, the estimation process is often largely hidden from the user. Consequently, problems such as a program crashing may occur that leave a user baffled. THREaD Mapper Studio (http://cbr.jic.ac.uk/threadmapper) is a new web site that implements a novel, visual and interactive method for the estimation of genetic linkage maps from DNA markers. The rationale behind the web site is to make the estimation process as transparent and robust as possible, while also allowing users to use their expert knowledge during analysis. Indeed, the 3D visual nature of the tool allows users to spot features in a data set, such as outlying markers and potential structural rearrangements that could cause problems with the estimation procedure and to account for them in their analysis. Furthermore, THREaD Mapper Studio facilitates the visual comparison of genetic map solutions from third party software, aiding users in developing robust solutions for their data sets.

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.

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.

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.

Computational approaches and software tools for genetic linkage map estimation in plants.

Genetic maps are an important component within the plant biologist's toolkit, underpinning crop plant improvement programs. The estimation of plant genetic maps is a conceptually simple yet computationally complex problem, growing ever more so with the development of inexpensive, high-throughput DNA markers. The challenge for bioinformaticians is to develop analytical methods and accompanying software tools that can cope with datasets of differing sizes, from tens to thousands of markers, that can incorporate the expert knowledge that plant biologists typically use when developing their maps, and that facilitate user-friendly approaches to achieving these goals. Here, we aim to give a flavour of computational approaches for genetic map estimation, discussing briefly many of the key concepts involved, and describing a selection of software tools that employ them. This review is intended both for plant geneticists as an introduction to software tools with which to estimate genetic maps, and for bioinformaticians as an introduction to the underlying computational approaches.

MAPS: an interactive web server for membrane annotation of transmembrane protein structures.

The exact positioning of the membrane in transmembrane (TM) proteins plays important functional roles. Yet, the structures of TM proteins in protein data bank (pdb) have no information about the explicit position of the membrane. Using a simple hydrophobic lipid-protein mismatch energy function and a flexible lipid/water boundary, the position of lipid bilayer for representative TM proteins in pdb have been annotated. A web server called MAPS (Membrane Annotation of Protein Structures; available at: http://www.boseinst.ernet.in/gautam/maps) has been set up that allows the user to interactively analyze membrane-protein orientations of any uploaded pdb structure with user-defined membrane flexibility parameters.

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.

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.