The rapid accumulation of mutations in the SARS-CoV-2 Omicron variant that enabled its outbreak raises questions as to whether its proximal origin occurred in humans or another mammalian host. Here, we identified 45 point mutations that Omicron acquired since divergence from the B.1.1 lineage. We found that the Omicron spike protein sequence was subjected to stronger positive selection than that of any reported SARS-CoV-2 variants known to evolve persistently in human hosts, suggesting a possibility of host-jumping. The molecular spectrum of mutations ( i.e. , the relative frequency of the 12 types of base substitutions) acquired by the progenitor of Omicron was significantly different from the spectrum for viruses that evolved in human patients, but resembled the spectra associated with virus evolution in a mouse cellular environment. Furthermore, mutations in the Omicron spike protein significantly overlapped with SARS-CoV-2 mutations known to promote adaptation to mouse hosts, particularly through enhanced spike protein binding affinity for the mouse cell entry receptor. Collectively, our results suggest that the progenitor of Omicron jumped from humans to mice, rapidly accumulated mutations conducive to infecting that host, then jumped back into humans, indicating an inter-species evolutionary trajectory for the Omicron outbreak.

BackgroundMicroRNAs (miRNAs) are endogenous small RNAs having large-scale regulatory effects on plant development and stress responses. Extensive studies of miRNAs have only been performed in a few model plants. Although miRNAs are proved to be involved in plant cold stress responses, little is known for winter-habit monocots. Brachypodium distachyon, with close evolutionary relationship to cool-season cereals, has recently emerged as a novel model plant. There are few reports of Brachypodium miRNAs.ResultsHigh-throughput sequencing and whole-genome-wide data mining led to the identification of 27 conserved miRNAs, as well as 129 predicted miRNAs in Brachypodium. For multiple-member conserved miRNA families, their sizes in Brachypodium were much smaller than those in rice and Populus. The genome organization of miR395 family in Brachypodium was quite different from that in rice. The expression of 3 conserved miRNAs and 25 predicted miRNAs showed significant changes in response to cold stress. Among these miRNAs, some were cold-induced and some were cold-suppressed, but all the conserved miRNAs were up-regulated under cold stress condition.ConclusionOur results suggest that Brachypodium miRNAs are composed of a set of conserved miRNAs and a large proportion of non-conserved miRNAs with low expression levels. Both kinds of miRNAs were involved in cold stress response, but all the conserved miRNAs were up-regulated, implying an important role for cold-induced miRNAs. The different size and genome organization of miRNA families in Brachypodium and rice suggest that the frequency of duplication events or the selection pressure on duplicated miRNAs are different between these two closely related plant species.

Vernalization, sensing of prolonged cold, is important for seasonal flowering in eudicots and monocots. While vernalization silences a repressor (FLC, MADS-box transcription factor) in eudicots, it induces an activator (TaVRN1, an AP1 clade MADS-box transcription factor) in monocots. The mechanism for TaVRN1 induction during vernalization is not well understood. Here we reveal a novel mechanism for controlling TaVRN1 mRNA accumulation in response to prolonged cold sensing in wheat. The carbohydrate-binding protein VER2, a jacalin lectin, promotes TaVRN1 upregulation by physically interacting with the RNA-binding protein TaGRP2. TaGRP2 binds to TaVRN1 pre-mRNA and inhibits TaVRN1 mRNA accumulation. The physical interaction between VER2 and TaGRP2 is controlled by TaGRP2 O-GlcNAc modification, which gradually increases during vernalization. The interaction between VER2 and O-GlcNAc-TaGRP2 reduces TaGRP2 protein accumulation in the nucleus and/or promotes TaGRP2 dissociation from TaVRN1, leading to TaVRN1 mRNA accumulation. Our data reveal a new mechanism for sensing prolonged cold in temperate cereals.

Brassinosteroids (BRs) regulate rice plant architecture, including leaf bending, which affects grain yield. Although BR signaling has been investigated in Arabidopsis thaliana, the components negatively regulating this pathway are less well understood. Here, we demonstrate that Oryza sativa LEAF and TILLER ANGLE INCREASED CONTROLLER (LIC) acts as an antagonistic transcription factor of BRASSINAZOLE-RESISTANT 1 (BZR1) to attenuate the BR signaling pathway. The gain-of-function mutant lic-1 and LIC–overexpressing lines showed erect leaves, similar to BZR1–depleted lines, which indicates the opposite roles of LIC and BZR1 in regulating leaf bending. Quantitative PCR revealed LIC transcription rapidly induced by BR treatment. Image analysis and immunoblotting showed that upon BR treatment LIC proteins translocate from the cytoplasm to the nucleus in a phosphorylation-dependent fashion. Phosphorylation assay in vitro revealed LIC phosphorylated by GSK3–like kinases. For negative feedback, LIC bound to the core element CTCGC in the BZR1 promoter on gel-shift and chromatin immunoprecipitation assay and repressed its transcription on transient transformation assay. LIC directly regulated target genes such as INCREASED LEAF INCLINATION 1 (ILI1) to oppose the action of BZR1. Repression of LIC in ILI1 transcription in protoplasts was partially rescued by BZR1. Phenotypic analysis of the crossed lines depleted in both LIC and BZR1 suggested that BZR1 functionally depends on LIC. Molecular and physiology assays revealed that LIC plays a dominant role at high BR levels, whereas BZR1 is dominant at low levels. Thus, LIC regulates rice leaf bending as an antagonistic transcription factor of BZR1. The phenotypes of lic-1 and LIC–overexpressing lines in erect leaves contribute to ideal plant architecture. Improving this phenotype may be a potential approach to molecular breeding for high yield in rice.

The rapid accumulation of mutations in the SARS-CoV-2 Omicron variant that enabled its outbreak raises questions as to whether its proximal origin occurred in humans or another mammalian host. Here, we identified 45 point mutations that Omicron acquired since divergence from the B.1.1 lineage. We found that the Omicron spike protein sequence was subjected to stronger positive selection than that of any reported SARS-CoV-2 variants known to evolve persistently in human hosts, suggesting the possibility of host-jumping. The molecular spectrum (i.e., the relative frequency of the twelve types of base substitutions) of mutations acquired by the progenitor of Omicron was significantly different from the spectrum for viruses that evolved in human patients, but was highly consistent with spectra associated with evolution in a mouse cellular environment. Furthermore, mutations in the Omicron spike protein significantly overlapped with SARS-CoV-2 mutations known to promote adaptation to mouse hosts, particularly through enhanced spike protein binding affinity for the mouse cell entry receptor. Collectively, our results suggest that the progenitor of Omicron jumped from humans to mice, rapidly accumulated mutations conducive to infecting that host, then jumped back into humans, indicating an inter-species evolutionary trajectory for the Omicron outbreak.

No abstract

Background Polyadenylation plays a key role in producing mature mRNAs in eukaryotes. It is widely believed that the poly(A)-binding proteins (PABs) uniformly bind to poly(A)-tailed mRNAs, regulating their stability and translational efficiency. Results We observe that the homozygous triple mutant of broadly expressed Arabidopsis thaliana PABs, AtPAB2, AtPAB4, and AtPAB8, is embryonic lethal. To understand the molecular basis, we characterize the RNA-binding landscape of these PABs. The AtPAB-binding efficiency varies over one order of magnitude among genes. To identify the sequences accounting for the variation, we perform poly(A)-seq that directly sequences the full-length poly(A) tails. More than 10% of poly(A) tails contain at least one guanosine (G); among them, the G-content varies from 0.8 to 28%. These guanosines frequently divide poly(A) tails into interspersed A-tracts and therefore cause the variation in the AtPAB-binding efficiency among genes. Ribo-seq and genome-wide RNA stability assays show that AtPAB-binding efficiency of a gene is positively correlated with translational efficiency rather than mRNA stability. Consistently, genes with stronger AtPAB binding exhibit a greater reduction in translational efficiency when AtPAB is depleted. Conclusions Our study provides a new mechanism that translational efficiency of a gene can be regulated through the G-content-dependent PAB binding, paving the way for a better understanding of poly(A) tail-associated regulation of gene expression. Electronic supplementary material The online version of this article (10.1186/s13059-019-1799-8) contains supplementary material, which is available to authorized users.

Codon usage bias (CUB) refers to the observation that synonymous codons are not used equally frequently in a genome. CUB is stronger in more highly expressed genes, a phenomenon commonly explained by stronger natural selection on translational accuracy and/or efficiency among these genes. Nevertheless, this phenomenon could also occur if CUB regulates gene expression at the mRNA level, a hypothesis that has not been tested until recently. Here, we attempt to quantify the impact of synonymous mutations on mRNA level in yeast using 3,556 synonymous variants of a heterologous gene encoding green fluorescent protein (GFP) and 523 synonymous variants of an endogenous gene TDH3. We found that mRNA level was positively correlated with CUB among these synonymous variants, demonstrating a direct role of CUB in regulating transcript concentration, likely via regulating mRNA degradation rate, as our additional experiments suggested. More importantly, we quantified the effects of individual synonymous mutations on mRNA level and found them dependent on 1) CUB and 2) mRNA secondary structure, both in proximal sequence contexts. Our study reveals the pleiotropic effects of synonymous codon usage and provides an additional explanation for the well-known correlation between CUB and gene expression level.