New Insights on Metabolic Features of Bacillus subtilis Based on Multistrain Genome-Scale Metabolic Modeling.

Bacillus subtilis is an effective workhorse for the production of many industrial products. The high interest aroused by B. subtilis has guided a large metabolic modeling effort of this species. Genome-scale metabolic models (GEMs) are powerful tools for predicting the metabolic capabilities of a given organism. However, high-quality GEMs are required in order to provide accurate predictions. In this work, we construct a high-quality, mostly manually curated genome-scale model for B. subtilis (iBB1018). The model was validated by means of growth performance and carbon flux distribution and provided significantly more accurate predictions than previous models. iBB1018 was able to predict carbon source utilization with great accuracy while identifying up to 28 metabolites as potential novel carbon sources. The constructed model was further used as a tool for the construction of the panphenome of B. subtilis as a species, by means of multistrain genome-scale reconstruction. The panphenome space was defined in the context of 183 GEMs representative of 183 B. subtilis strains and the array of carbon sources sustaining growth. Our analysis highlights the large metabolic versatility of the species and the important role of the accessory metabolism as a driver of the panphenome, at a species level.

sRNA Analysis Evidenced the Involvement of Different Plant Viruses in the Activation of RNA Silencing-Related Genes and the Defensive Response Against Plum pox virus of 'GF305' Peach Grafted with 'Garrigues' Almond.

Plum pox virus (PPV) causes sharka disease in Prunus trees. Peach (P. persica) trees are severely affected by PPV, and no definitive source of genetic resistance has been identified. However, previous results showed that PPV-resistant 'Garrigues' almond (P. dulcis) was able to transfer its resistance to 'GF305' peach through grafting, reducing symptoms and viral load in PPV-infected plants. A recent study tried to identify genes responsible for this effect by studying messenger RNA expression through RNA sequencing in peach and almond plants, before and after grafting and before and after PPV infection. In this work, we used the same peach and almond samples but focused the high-throughput analyses on small RNA (sRNA) expression. We studied massive sequencing data and found an interesting pattern of sRNA overexpression linked to antiviral defense genes that suggested activation of these genes followed by downregulation to basal levels. We also discovered that 'Garrigues' almond plants were infected by different plant viruses that were transferred to peach plants. The large amounts of viral sRNA found in grafted peaches indicated a strong RNA silencing antiviral response and led us to postulate that these plant viruses could be collaborating in the observed "Garrigues effect."

Breaking-Cas-interactive design of guide RNAs for CRISPR-Cas experiments for ENSEMBL genomes.

The CRISPR/Cas technology is enabling targeted genome editing in multiple organisms with unprecedented accuracy and specificity by using RNA-guided nucleases. A critical point when planning a CRISPR/Cas experiment is the design of the guide RNA (gRNA), which directs the nuclease and associated machinery to the desired genomic location. This gRNA has to fulfil the requirements of the nuclease and lack homology with other genome sites that could lead to off-target effects. Here we introduce the Breaking-Cas system for the design of gRNAs for CRISPR/Cas experiments, including those based in the Cas9 nuclease as well as others recently introduced. The server has unique features not available in other tools, including the possibility of using all eukaryotic genomes available in ENSEMBL (currently around 700), placing variable PAM sequences at 5' or 3' and setting the guide RNA length and the scores per nucleotides. It can be freely accessed at: http://bioinfogp.cnb.csic.es/tools/breakingcas, and the code is available upon request.

The P1N-PISPO trans-Frame Gene of Sweet Potato Feathery Mottle Potyvirus Is Produced during Virus Infection and Functions as an RNA Silencing Suppressor.

UNLABELLED: The positive-sense RNA genome of Sweet potato feathery mottle virus (SPFMV) (genus Potyvirus, family Potyviridae) contains a large open reading frame (ORF) of 3,494 codons translatable as a polyprotein and two embedded shorter ORFs in the -1 frame: PISPO, of 230 codons, and PIPO, of 66 codons, located in the P1 and P3 regions, respectively. PISPO is specific to some sweet potato-infecting potyviruses, while PIPO is present in all potyvirids. In SPFMV these two extra ORFs are preceded by conserved G2A6 motifs. We have shown recently that a polymerase slippage mechanism at these sites could produce transcripts bringing these ORFs in frame with the upstream polyprotein, thus leading to P1N-PISPO and P3N-PIPO products (B. Rodamilans, A. Valli, A. Mingot, D. San Leon, D. B. Baulcombe, J. J. Lopez-Moya, and J.A. Garcia, J Virol 89:6965-6967, 2015, doi:10.1128/JVI.00337-15). Here, we demonstrate by liquid chromatography coupled to mass spectrometry that both P1 and P1N-PISPO are produced during viral infection and coexist in SPFMV-infected Ipomoea batatas plants. Interestingly, transient expression of SPFMV gene products coagroinfiltrated with a reporter gene in Nicotiana benthamiana revealed that P1N-PISPO acts as an RNA silencing suppressor, a role normally associated with HCPro in other potyviruses. Moreover, mutation of WG/GW motifs present in P1N-PISPO abolished its silencing suppression activity, suggesting that the function might require interaction with Argonaute components of the silencing machinery, as was shown for other viral suppressors. Altogether, our results reveal a further layer of complexity of the RNA silencing suppression activity within the Potyviridae family. IMPORTANCE: Gene products of potyviruses include P1, HCPro, P3, 6K1, CI, 6K2, VPg/NIaPro, NIb, and CP, all derived from the proteolytic processing of a large polyprotein, and an additional P3N-PIPO product, with the PIPO segment encoded in a different frame within the P3 cistron. In sweet potato feathery mottle virus (SPFMV), another out-of-frame element (PISPO) was predicted within the P1 region. We have shown recently that a polymerase slippage mechanism can generate the transcript variants with extra nucleotides that could be translated into P1N-PISPO and P3N-PIPO. Now, we demonstrate by mass spectrometry analysis that P1N-PISPO is indeed produced in SPFMV-infected plants, in addition to P1. Interestingly, while in other potyviruses the suppressor of RNA silencing is HCPro, we show here that P1N-PISPO exhibited this activity in SPFMV, revealing how the complexity of the gene content could contribute to supply this essential function in members of the Potyviridae family.

Virus-induced gene silencing in transgenic plants: transgene silencing and reactivation associate with two patterns of transgene body methylation.

We used bisulfite sequencing to study the methylation of a viral transgene whose expression was silenced upon plum pox virus infection of the transgenic plant and its subsequent recovery as a consequence of so-called virus-induced gene silencing (VIGS). VIGS was associated with a general increase in the accumulation of small RNAs corresponding to the coding region of the viral transgene. After VIGS, the transgene promoter was not methylated and the coding region showed uneven methylation, with the 5' end being mostly unmethylated in the recovered tissue or mainly methylated at CG sites in regenerated silenced plants. The methylation increased towards the 3' end, which showed dense methylation in all three contexts (CG, CHG and CHH). This methylation pattern and the corresponding silenced status were maintained after plant regeneration from recovered silenced tissue and did not spread into the promoter region, but were not inherited in the sexual offspring. Instead, a new pattern of methylation was observed in the progeny plants consisting of disappearance of the CHH methylation, similar CHG methylation at the 3' end, and an overall increase in CG methylation in the 5' end. The latter epigenetic state was inherited over several generations and did not correlate with transgene silencing and hence virus resistance. These results suggest that the widespread CG methylation pattern found in body gene bodies located in euchromatic regions of plant genomes may reflect an older silencing event, and most likely these genes are no longer silenced.

Methyltransferase-like (METTL) homologues participate in Nicotiana benthamiana antiviral responses.

Methyltransferase (MTase) enzymes catalyze the addition of a methyl group to a variety of biological substrates. MTase-like (METTL) proteins are Class I MTases whose enzymatic activities contribute to the epigenetic and epitranscriptomic regulation of multiple cellular processes. N6-adenosine methylation (m6A) is a common chemical modification of eukaryotic and viral RNA whose abundance is jointly regulated by MTases and METTLs, demethylases, and m6A binding proteins. m6A affects various cellular processes including RNA degradation, post-transcriptional processing, and antiviral immunity. Here, we used Nicotiana benthamiana and plum pox virus (PPV), an RNA virus of the Potyviridae family, to investigated the roles of MTases in plant-virus interaction. RNA sequencing analysis identified MTase transcripts that are differentially expressed during PPV infection; among these, accumulation of a METTL gene was significantly downregulated. Two N. benthamiana METTL transcripts (NbMETTL1 and NbMETTL2) were cloned and further characterized. Sequence and structural analyses of the two encoded proteins identified a conserved S-adenosyl methionine (SAM) binding domain, showing they are SAM-dependent MTases phylogenetically related to human METTL16 and Arabidopsis thaliana FIONA1. Overexpression of NbMETTL1 and NbMETTL2 caused a decrease of PPV accumulation. In sum, our results indicate that METTL homologues participate in plant antiviral responses.

New strategy for bioplastic and exopolysaccharides production: Enrichment of field microbiomes with cyanobacteria.

Seven photosynthethic microbiomes were collected from field environmental samples to test their potential in polyhydroxybutyrate (PHB) and exopolysaccharides (EPS) production, both alternatives to chemical-based polymers. Microscope observations together with microbial sequence analysis revealed the microbiome enrichment in cyanobacteria after culture growth under phosphorus limitation. PHB and EPS production were studied under three culture factors (phototrophy, mixotrophy and heterotrophy) by evaluating and optimizing the effect of three parameters (organic and inorganic carbon and days under light:dark cycles) by Box-Behnken design. Results showed that optimal conditions for both biopolymers synthesis were microbiome-dependent; however, the addition of organic carbon boosted PHB production in all the tested microbiomes, producing up to 14 %dcw PHB with the addition of 1.2 g acetate·L-1 and seven days under light:dark photoperiods. The highest EPS production was 59 mg·L-1 with the addition of 1.2 g acetate·L-1 and four days under light:dark photoperiods. The methodology used is suitable for enriching microbiomes in cyanobacteria, and for testing the best conditions for bioproduct synthesis for further scale up.

Toward merging bottom-up and top-down model-based designing of synthetic microbial communities.

The increasing interest of microbial communities as promising biocatalyst is leading an intense effort into the development of computational frameworks assisting the analysis and rational engineering of such complex ecosystems. Here, we critically review the recent computational and model-guided advances in the system-level engineering of microbiome, including both the rational bottom-up and the evolutionary top-down approaches. Furthermore, we highlight modeling and computational methods supporting both engineering paradigms. Finally, we discuss the advantages of combining both strategies into a hybrid top-down/bottom-up (middle-out) strategy to engineer synthetic microbial communities with improved performance and scope.

Repositioning microbial biotechnology against COVID-19: the case of microbial production of flavonoids.

Coronavirus-related disease 2019 (COVID-19) became a pandemic in February 2020, and worldwide researchers try to tackle the disease with approved drugs of all kinds, or to develop novel compounds inhibiting viral spreading. Flavonoids, already investigated as antivirals in general, also might bear activities specific for the viral agent causing COVID-19, SARS-CoV-2. Microbial biotechnology and especially synthetic biology may help to produce flavonoids, which are exclusive plant secondary metabolites, at a larger scale or indeed to find novel pharmaceutically active flavonoids. Here, we review the state of the art in (i) antiviral activity of flavonoids specific for coronaviruses and (ii) results derived from computational studies, mostly docking studies mainly inhibiting specific coronaviral proteins such as the 3CL (main) protease, the spike protein or the RNA-dependent RNA polymerase. In the end, we strive towards a synthetic biology pipeline making the fast and tailored production of valuable antiviral flavonoids possible by applying the last concepts of division of labour through co-cultivation/microbial community approaches to the DBTL (Design, Build, Test, Learn) principle.

Plant Virus Genome Is Shaped by Specific Dinucleotide Restrictions That Influence Viral Infection.

The presence of CpG and UpA dinucleotides is restricted in the genomes of animal RNA viruses to avoid specific host defenses. We wondered whether a similar phenomenon exists in nonanimal RNA viruses. Here, we show that these two dinucleotides, especially UpA, are underrepresented in the family Potyviridae, the most important group of plant RNA viruses. Using plum pox virus (PPV; Potyviridae family) as a model, we show that an increase in UpA frequency strongly diminishes virus accumulation. Remarkably, unlike previous observations in animal viruses, PPV variants harboring CpG-rich fragments display just faint (or no) attenuation. The anticorrelation between UpA frequency and viral fitness additionally demonstrates the relevance of this particular dinucleotide: UpA-high mutants are attenuated in a dose-dependent manner, whereas a UpA-low variant displays better fitness than its parental control. Using high-throughput sequencing, we also show that UpA-rich PPV variants are genetically stable, without apparent changes in sequence that revert and/or compensate for the dinucleotide modification despite its attenuation. In addition, we also demonstrate here that the PPV restriction of UpA-rich variants works independently of the classical RNA silencing pathway. Finally, we show that the anticorrelation between UpA frequency and RNA accumulation applies to mRNA-like fragments produced by the host RNA polymerase II. Together, our results inform us about a dinucleotide-based system in plant cells that controls diverse RNAs, including RNA viruses.IMPORTANCE Dinucleotides (combinations of two consecutive nucleotides) are not randomly present in RNA viruses; in fact, the presence of CpG and UpA is significantly repressed in their genomes. Although the meaning of this phenomenon remains obscure, recent studies with animal-infecting viruses have revealed that their low CpG/UpA frequency prevents virus restriction via a host antiviral system that recognizes, and promotes the degradation of, CpG/UpA-rich RNAs. Whether similar systems act in organisms from other life kingdoms has been unknown. To fill this gap in our knowledge, we built several synthetic variants of a plant RNA virus with deoptimized dinucleotide frequencies and analyzed their viral fitness and genome adaptation. In brief, our results inform us for the first time about an effective dinucleotide-based system that acts in plants against viruses. Remarkably, this viral restriction in plants is reminiscent of, but not identical to, the equivalent antiviral response in animals.

Abscisic Acid Connects Phytohormone Signaling with RNA Metabolic Pathways and Promotes an Antiviral Response that Is Evaded by a Self-Controlled RNA Virus.

A complex network of cellular receptors, RNA targeting pathways, and small-molecule signaling provides robust plant immunity and tolerance to viruses. To maximize their fitness, viruses must evolve control mechanisms to balance host immune evasion and plant-damaging effects. The genus Potyvirus comprises plant viruses characterized by RNA genomes that encode large polyproteins led by the P1 protease. A P1 autoinhibitory domain controls polyprotein processing, the release of a downstream functional RNA-silencing suppressor, and viral replication. Here, we show that P1Pro, a plum pox virus clone that lacks the P1 autoinhibitory domain, triggers complex reprogramming of the host transcriptome and high levels of abscisic acid (ABA) accumulation. A meta-analysis highlighted ABA connections with host pathways known to control RNA stability, turnover, maturation, and translation. Transcriptomic changes triggered by P1Pro infection or ABA showed similarities in host RNA abundance and diversity. Genetic and hormone treatment assays showed that ABA promotes plant resistance to potyviral infection. Finally, quantitative mathematical modeling of viral replication in the presence of defense pathways supported self-control of polyprotein processing kinetics as a viral mechanism that attenuates the magnitude of the host antiviral response. Overall, our findings indicate that ABA is an active player in plant antiviral immunity, which is nonetheless evaded by a self-controlled RNA virus.

An atypical RNA silencing suppression strategy provides a snapshot of the evolution of sweet potato-infecting potyviruses.

Plant viruses usually encode proteins with RNA silencing suppression (RSS) activity to counteract plant defenses. In Potyvirus, the largest genus in the family Potyviridae, this role is taken over by the multifunctional HCPro, also involved in aphid transmission, polyprotein processing and virion formation. Recently, the large P1 of Sweet potato feathery mottle virus (SPFMV) was characterized finding an extra ORF produced after polymerase slippage, which originates the product P1N-PISPO. Transient expression assays showed that SPFMV P1 and P1N-PISPO presented RSS activity, while HCPro did not. In this work, we analyze possible differences between HCPro of SPFMV and other potyviruses, testing HCPro RSS activity in a transient expression assay, and using a Plum pox virus-based system to test the ability of SPFMV P1N-PISPO and HCPro to serve as RNA silencing suppressors in the context of a viral infection. Our results indicate that not only P1 and P1N-PISPO, but also HCPro display RSS activity when expressed in a suitable context, stressing the importance of the selected experimental system for testing anti-silencing capacity of proteins. The presence of multiple viral silencing suppressors in SPFMV adds complexity to an already intricate RSS system, and provides insight into the hypothetical evolution of sweet potato-infecting potyvirids.

Assorted Processing of Synthetic Trans-Acting siRNAs and Its Activity in Antiviral Resistance.

The use of syn-tasiRNAs has been proposed as an RNA interference technique alternative to those previously described: hairpin based, virus induced gene silencing or artificial miRNAs. In this study we engineered the TAS1c locus to impair Plum pox virus (PPV) infection by replacing the five native siRNAs with two 210-bp fragments from the CP and the 3´NCR regions of the PPV genome. Deep sequencing analysis of the small RNA species produced by both constructs in planta has shown that phased processing of the syn-tasiRNAs is construct-specific. While in syn-tasiR-CP construct the processing was as predicted 21-nt phased in register with miR173-guided cleavage, the processing of syn-tasiR-3NCR is far from what was expected. A 22-nt species from the miR173-guided cleavage was a guide of two series of phased small RNAs, one of them in an exact 21-nt register, and the other one in a mixed of 21-/22-nt frame. In addition, both constructs produced abundant PPV-derived small RNAs in the absence of miR173 as a consequence of a strong sense post-transcriptional gene silencing induction. The antiviral effect of both constructs was also evaluated in the presence or absence of miR173 and showed that the impairment of PPV infection was not significantly higher when miR173 was present. The results show that syn-tasiRNAs processing depends on construct-specific factors that should be further studied before the so-called MIGS (miRNA-induced gene silencing) technology can be used reliably.

Transcriptomic analysis of Prunus domestica undergoing hypersensitive response to plum pox virus infection.

Plum pox virus (PPV) infects Prunus trees around the globe, posing serious fruit production problems and causing severe economic losses. One variety of Prunus domestica, named 'Jojo', develops a hypersensitive response to viral infection. Here we compared infected and non-infected samples using next-generation RNA sequencing to characterize the genetic complexity of the viral population in infected samples and to identify genes involved in development of the resistance response. Analysis of viral reads from the infected samples allowed reconstruction of a PPV-D consensus sequence. De novo reconstruction showed a second viral isolate of the PPV-Rec strain. RNA-seq analysis of PPV-infected 'Jojo' trees identified 2,234 and 786 unigenes that were significantly up- or downregulated, respectively (false discovery rate; FDR≤0.01). Expression of genes associated with defense was generally enhanced, while expression of those related to photosynthesis was repressed. Of the total of 3,020 differentially expressed unigenes, 154 were characterized as potential resistance genes, 10 of which were included in the NBS-LRR type. Given their possible role in plant defense, we selected 75 additional unigenes as candidates for further study. The combination of next-generation sequencing and a Prunus variety that develops a hypersensitive response to PPV infection provided an opportunity to study the factors involved in this plant defense mechanism. Transcriptomic analysis presented an overview of the changes that occur during PPV infection as a whole, and identified candidates suitable for further functional characterization.