Genome‐scale metabolic modeling reveals key features of a minimal gene set

Mesoplasma florum, a fast‐growing near‐minimal organism, is a compelling model to explore rational genome designs. Using sequence and structural homology, the set of metabolic functions its genome encodes was identified, allowing the reconstruction of a metabolic network representing ˜ 30% of its protein‐coding genes. Growth medium simplification enabled substrate uptake and product secretion rate quantification which, along with experimental biomass composition, were integrated as species‐specific constraints to produce the functional iJL208 genome‐scale model (GEM) of metabolism. Genome‐wide expression and essentiality datasets as well as growth data on various carbohydrates were used to validate and refine iJL208. Discrepancies between model predictions and observations were mechanistically explained using protein structures and network analysis. iJL208 was also used to propose an in silico reduced genome. Comparing this prediction to the minimal cell JCVI‐syn3.0 and its parent JCVI‐syn1.0 revealed key features of a minimal gene set. iJL208 is a stepping‐stone toward model‐driven whole‐genome engineering.     

Synergy: A Web Resource for Exploring Gene Regulation in Synechocystis sp. PCC6803

Despite being a highly studied model organism, most genes of the cyanobacterium Synechocystis sp. PCC 6803 encode proteins with completely unknown function. To facilitate studies of gene regulation in Synechocystis, we have developed Synergy (http://synergy.plantgenie.org), a web application integrating co-expression networks and regulatory motif analysis. Co-expression networks were inferred from publicly available microarray experiments, while regulatory motifs were identified using a phylogenetic footprinting approach. Automatically discovered motifs were shown to be enriched in the network neighborhoods of regulatory proteins much more often than in the neighborhoods of non-regulatory genes, showing that the data provide a sound starting point for studying gene regulation in Synechocystis. Concordantly, we provide several case studies demonstrating that Synergy can be used to find biologically relevant regulatory mechanisms in Synechocystis. Synergy can be used to interactively perform analyses such as gene/motif search, network visualization and motif/function enrichment. Considering the importance of Synechocystis for photosynthesis and biofuel research, we believe that Synergy will become a valuable resource to the research community.     

Protein structure,amino acid composition and sequence determine proteome vulnerability to oxidation‐induced damage

Oxidative stress alters cell viability, from microorganism irradiation sensitivity to human aging and neurodegeneration. Deleterious effects of protein carbonylation by reactive oxygen species (ROS) make understanding molecular properties determining ROS susceptibility essential. The radiation‐resistant bacterium Deinococcus radiodurans accumulates less carbonylation than sensitive organisms, making it a key model for deciphering properties governing oxidative stress resistance. We integrated shotgun redox proteomics, structural systems biology, and machine learning to resolve properties determining protein damage by γ‐irradiation in Escherichia coli and D. radiodurans at multiple scales. Local accessibility, charge, and lysine enrichment accurately predict ROS susceptibility. Lysine, methionine, and cysteine usage also contribute to ROS resistance of the D. radiodurans proteome. Our model predicts proteome maintenance machinery, and proteins protecting against ROS are more resistant in D. radiodurans. Our findings substantiate that protein‐intrinsic protection impacts oxidative stress resistance, identifying causal molecular properties.     

TrOn: An Anatomical Ontology for the Beetle Tribolium castaneum

In a morphological ontology the expert’s knowledge is represented in terms, which describe morphological structures and how these structures relate to each other. With the assistance of ontologies this expert knowledge is made processable by machines, through a formal and standardized representation of terms and their relations to each other. The red flour beetle Tribolium castaneum, a representative of the most species rich animal taxon on earth (the Coleoptera), is an emerging model organism for development, evolution, physiology, and pest control. In order to foster Tribolium research, we have initiated the Tribolium Ontology (TrOn), which describes the morphology of the red flour beetle. The content of this ontology comprises so far most external morphological structures as well as some internal ones. All modeled structures are consistently annotated for the developmental stages larva, pupa and adult. In TrOn all terms are grouped into three categories: Generic terms represent morphological structures, which are independent of a developmental stage. In contrast, downstream of such terms are concrete terms which stand for a dissectible structure of a beetle at a specific life stage. Finally, there are mixed terms describing structures that are only found at one developmental stage. These terms combine the characteristics of generic and concrete terms with features of both. These annotation principles take into account the changing morphology of the beetle during development and provide generic terms to be used in applications or for cross linking with other ontologies and data resources. We use the ontology for implementing an intuitive search function at the electronic iBeetle-Base, which stores morphological defects found in a genome wide RNA interference (RNAi) screen. The ontology is available for download at http://ibeetle-base.uni-goettingen.de.     

Membranome 3.0: Database of single‐pass membrane proteins with AlphaFold models

The Membranome database provides comprehensive structural information on single‐pass (i.e., bitopic) membrane proteins from six evolutionarily distant organisms, including protein–protein interactions, complexes, mutations, experimental structures, and models of transmembrane α‐helical dimers. We present a new version of this database, Membranome 3.0, which was significantly updated by revising the set of 5,758 bitopic proteins and incorporating models generated by AlphaFold 2 in the database. The AlphaFold models were parsed into structural domains located at the different membrane sides, modified to exclude low‐confidence unstructured terminal regions and signal sequences, validated through comparison with available experimental structures, and positioned with respect to membrane boundaries. Membranome 3.0 was re‐developed to facilitate visualization and comparative analysis of multiple 3D structures of proteins that belong to a specified family, complex, biological pathway, or membrane type. New tools for advanced search and analysis of proteins, their interactions, complexes, and mutations were included. The database is freely accessible at https://membranome.org.     

Evolutionary patterns of nucleotide substitution rates in plastid genomes of Quercus

Molecular evolution, including nucleotide substitutions, plays an important role in understanding the dynamics and mechanisms of species evolution. Here, we sequenced whole plastid genomes (plastomes) of Quercus fabri, Quercus semecarpifolia, Quercus engleriana, and Quercus phellos and compared them with 14 other Quercus plastomes to explore their evolutionary relationships using 67 shared protein‐coding sequences. While many previously identified evolutionary relationships were found, our findings do not support previous research which retrieve Quercus subg. Cerris sect. Ilex as a monophyletic group, with sect. Ilex found to be polyphyletic and composed of three strongly supported lineages inserted between sections Cerris and Cyclobalanposis. Compared with gymnosperms, Quercus plastomes showed higher evolutionary rates (Dn/Ds = 0.3793). Most protein‐coding genes experienced relaxed purifying selection, and the high Dn value (0.1927) indicated that gene functions adjusted to environmental changes effectively. Our findings suggest that gene interval regions play an important role in Quercus evolution. We detected greater variation in the intergenic regions (trnH‐psbA, trnK_UUU‐rps16, trnfM_CAU‐rps14, trnS_GCU‐trnG_GCC, and atpF‐atpH), intron losses (petB and petD), and pseudogene loss and degradation (ycf15). Additionally, the loss of some genes suggested the existence of gene exchanges between plastid and nuclear genomes, which affects the evolutionary rate of the former. However, the connective mechanism between these two genomes is still unclear.     

Integrative analysis provides multi‐omics evidence for the pathogenesis of placenta percreta

Pernicious placenta previa with placenta percreta (PP) is a catastrophic condition during pregnancy. However, the underlying pathogenesis remains unclear. In the present study, the placental tissues of normal cases and PP tissues of pernicious placenta previa cases were collected to determine the expression profile of protein‐coding genes, miRNAs, and lncRNAs through sequencing. Weighted gene co‐expression network analysis (WGCNA), accompanied by miRNA target prediction and correlation analysis, were employed to select potential hub protein‐coding genes and lncRNAs. The expression levels of selected protein‐coding genes, Wnt5A and MAPK13, were determined by quantitative PCR and immunohistochemical staining, and lncRNA PTCHD1‐AS and PAPPA‐AS1 expression levels were determined by quantitative PCR and fluorescence in situ hybridization. The results indicated that 790 protein‐coding genes, 382 miRNAs, and 541 lncRNAs were dysregulated in PP tissues, compared with normal tissues. WGCNA identified coding genes in the module (ME) black and ME turquoise modules that may be involved in the pathogenesis of PP. The selected potential hub protein‐coding genes, Wnt5A and MAPK13, were down‐regulated in PP tissues, and their expression levels were positively correlated with the expression levels of PTCHD1‐AS and PAPPA‐AS1. Further analysis demonstrated that PTCHD1‐AS and PAPPA‐AS1 regulated Wnt5A and MAPK13 expression by interacting with specific miRNAs. Collectively, our results provided multi‐omics data to better understand the pathogenesis of PP and help identify predictive biomarkers and therapeutic targets for PP.     

Engineering a genome‐reduced bacterium to eliminate Staphylococcus aureus biofilms in vivo

Bacteria present a promising delivery system for treating human diseases. Here, we engineered the genome‐reduced human lung pathogen Mycoplasma pneumoniae as a live biotherapeutic to treat biofilm‐associated bacterial infections. This strain has a unique genetic code, which hinders gene transfer to most other bacterial genera, and it lacks a cell wall, which allows it to express proteins that target peptidoglycans of pathogenic bacteria. We first determined that removal of the pathogenic factors fully attenuated the chassis strain in vivo. We then designed synthetic promoters and identified an endogenous peptide signal sequence that, when fused to heterologous proteins, promotes efficient secretion. Based on this, we equipped the chassis strain with a genetic platform designed to secrete antibiofilm and bactericidal enzymes, resulting in a strain capable of dissolving Staphylococcus aureus biofilms preformed on catheters in vitro, ex vivo, and in vivo. To our knowledge, this is the first engineered genome‐reduced bacterium that can fight against clinically relevant biofilm‐associated bacterial infections.     

NMR hawk‐eyed view of AlphaFold2 structures

The prediction of the three‐dimensional (3D) structure of proteins from the amino acid sequence made a stunning breakthrough reaching atomic accuracy. Using the neural network‐based method AlphaFold2, 3D structures of almost the entire human proteome have been predicted and made available (https://www.alphafold.ebi.ac.uk). To gain insight into how well AlphaFold2 structures represent the conformation of proteins in solution, I here compare the AlphaFold2 structures of selected small proteins with their 3D structures that were determined by nuclear magnetic resonance (NMR) spectroscopy. Proteins were selected for which the 3D solution structures were determined on the basis of a very large number of distance restraints and residual dipolar couplings and are thus some of the best‐resolved solution structures of proteins to date. The quality of the backbone conformation of the AlphaFold2 structures is assessed by fitting a large set of experimental residual dipolar couplings (RDCs). The analysis shows that experimental RDCs fit extremely well to the AlphaFold2 structures predicted for GB3, DinI, and ubiquitin. In the case of GB3, the accuracy of the AlphaFold2 structure even surpasses that of a 1.1 Å crystal structure. Fitting of experimental RDCs furthermore allows identification of AlphaFold2 structures that are best representative of the protein''s conformation in solution as seen for the EF hands of the N‐terminal domain of Ca²⁺‐ligated calmodulin. Taken together, the analysis shows that structures predicted by AlphaFold2 can be highly representative of the solution conformation of proteins. The combination of AlphaFold2 structures with RDCs promises to be a powerful approach to study structural changes in proteins.     

Clinical translation of mesenchymal stromal cell extracellular vesicles: Considerations on scientific rationale and production requisites

The present paper is a commentary to ‘Identification and characterization of hADSC‐derived exosome proteins from different isolation methods’ (Huang et al. 2021; 10.1111/jcmm.16775). Given the enthusiasm for the potential of mesenchymal stromal cell‐derived extracellular vesicles (MSC‐EVs), some considerations deserve attention as they move through successive stages of research and application into humans. We herein remark the prerequisite of generating that evidence ensuring a high consistency in safety, composition and biological activity of the intended MSC‐EV preparations, and the suitability of disparate isolation techniques to produce efficacious EV preparations and fulfil requirements for standardized clinical‐grade biomanufacturing.     

Characterization of an aminotransferase from Acanthamoeba polyphaga Mimivirus

Acanthamoeba polyphaga Mimivirus, a complex virus that infects amoeba, was first reported in 2003. It is now known that its DNA genome encodes for nearly 1,000 proteins including enzymes that are required for the biosynthesis of the unusual sugar 4‐amino‐4,6‐dideoxy‐d‐glucose, also known as d‐viosamine. As observed in some bacteria, the pathway for the production of this sugar initiates with a nucleotide‐linked sugar, which in the Mimivirus is thought to be UDP‐d‐glucose. The enzyme required for the installment of the amino group at the C‐4′ position of the pyranosyl moiety is encoded in the Mimivirus by the L136 gene. Here, we describe a structural and functional analysis of this pyridoxal 5′‐phosphate‐dependent enzyme, referred to as L136. For this analysis, three high‐resolution X‐ray structures were determined: the wildtype enzyme/pyridoxamine 5′‐phosphate/dTDP complex and the site‐directed mutant variant K185A in the presence of either UDP‐4‐amino‐4,6‐dideoxy‐d‐glucose or dTDP‐4‐amino‐4,6‐dideoxy‐d‐glucose. Additionally, the kinetic parameters of the enzyme utilizing either UDP‐d‐glucose or dTDP‐d‐glucose were measured and demonstrated that L136 is efficient with both substrates. This is in sharp contrast to the structurally related DesI from Streptomyces venezuelae, whose three‐dimensional architecture was previously reported by this laboratory. As determined in this investigation,DesI shows a profound preference in its catalytic efficiency for the dTDP‐linked sugar substrate. This difference can be explained in part by a hydrophobic patch in DesI that is missing in L136. Notably, the structure of L136 reported here represents the first three‐dimensional model for a virally encoded PLP‐dependent enzyme and thus provides new information on sugar aminotransferases in general.     

Identification and Analysis of the Polyhydroxyalkanoate-Specific β-Ketothiolase and Acetoacetyl Coenzyme A Reductase Genes in the Cyanobacterium Synechocystis sp. Strain PCC6803

Synechocystis sp. strain PCC6803 possesses a polyhydroxyalkanoate (PHA)-specific β-ketothiolase encoded by phaA_Syn and an acetoacetyl-coenzyme A (CoA) reductase encoded by phaB_Syn. A similarity search of the entire Synechocystis genome sequence identified a cluster of two putative open reading frames (ORFs) for these genes, slr1993 and slr1994. Sequence analysis showed that the ORFs encode proteins having 409 and 240 amino acids, respectively. The two ORFs are colinear and most probably coexpressed, as revealed by sequence analysis of the promoter regions. Heterologous transformation of Escherichia coli with the two genes and the PHA synthase of Synechocystis resulted in accumulation of PHAs that accounted for up to 12.3% of the cell dry weight under high-glucose growth conditions. Targeted disruption of the above gene cluster in Synechocystis eliminated the accumulation of PHAs. ORFs slr1993 and slr1994 thus encode the PHA-specific β-ketothiolase and acetoacetyl-CoA reductase of Synechocystis and, together with the recently characterized PHA synthase genes in this organism (S. Hein, H. Tran, and A. Steinbüchel, Arch. Microbiol. 170:162–170, 1998), form the first complete PHA biosynthesis pathway known in cyanobacteria. Sequence alignment of all known short-chain-length PHA-specific acetoacetyl-CoA reductases also suggests an extended signature sequence, VTGXXXGIG, for this group of proteins. Phylogenetic analysis further places the origin of phaA_Syn and phaB_Syn in the γ subdivision of the division Proteobacteria.     

Xanthomonas hortorum – beyond gardens: Current taxonomy,genomics, and virulence repertoires

TaxonomyBacteria; Phylum Proteobacteria; Class Gammaproteobacteria; Order Lysobacterales (earlier synonym of Xanthomonadales); Family Lysobacteraceae (earlier synonym of Xanthomonadaceae); Genus Xanthomonas; Species X. hortorum; Pathovars: pv. carotae, pv. vitians, pv. hederae, pv. pelargonii, pv. taraxaci, pv. cynarae, and pv. gardneri.Host range Xanthomonas hortorum affects agricultural crops, and horticultural and wild plants. Tomato, carrot, artichoke, lettuce, pelargonium, ivy, and dandelion were originally described as the main natural hosts of the seven separate pathovars. Artificial inoculation experiments also revealed other hosts. The natural and experimental host ranges are expected to be broader than initially assumed. Additionally, several strains, yet to be assigned to a pathovar within X. hortorum, cause diseases on several other plant species such as peony, sweet wormwood, lavender, and oak‐leaf hydrangea.Epidemiology and control X. hortorum pathovars are mainly disseminated by infected seeds (e.g., X. hortorum pvs carotae and vitians) or cuttings (e.g., X. hortorum pv. pelargonii) and can be further dispersed by wind and rain, or mechanically transferred during planting and cultivation. Global trade of plants, seeds, and other propagating material constitutes a major pathway for their introduction and spread into new geographical areas. The propagules of some pathovars (e.g., X. horturum pv. pelargonii) are spread by insect vectors, while those of others can survive in crop residues and soils, and overwinter until the following growing season (e.g., X. hortorum pvs vitians and carotae). Control measures against X. hortorum pathovars are varied and include exclusion strategies (i.e., by using certification programmes and quarantine regulations) to multiple agricultural practices such as the application of phytosanitary products. Copper‐based compounds against X. hortorum are used, but the emergence of copper‐tolerant strains represents a major threat for their effective management. With the current lack of efficient chemical or biological disease management strategies, host resistance appears promising, but is not without challenges. The intrastrain genetic variability within the same pathovar poses a challenge for breeding cultivars with durable resistance.Useful websites https://gd.eppo.int/taxon/XANTGA, https://gd.eppo.int/taxon/XANTCR, https://gd.eppo.int/taxon/XANTPE, https://www.euroxanth.eu, http://www.xanthomonas.org, http://www.xanthomonas.org/dokuwiki     

PDBsum extras: SARS‐CoV‐2 and AlphaFold models

The PDBsum web server provides structural analyses of the entries in the Protein Data Bank (PDB). Two recent additions are described here. The first is the detailed analysis of the SARS‐CoV‐2 virus protein structures in the PDB. These include the variants of concern, which are shown both on the sequences and 3D structures of the proteins. The second addition is the inclusion of the available AlphaFold models for human proteins. The pages allow a search of the protein against existing structures in the PDB via the Sequence Annotated by Structure (SAS) server, so one can easily compare the predicted model against experimentally determined structures. The server is freely accessible to all at http://www.ebi.ac.uk/pdbsum.     

α‐synuclein suppresses microglial autophagy and promotes neurodegeneration in a mouse model of Parkinson’s disease

The cell‐to‐cell transfer of α‐synuclein (α‐Syn) greatly contributes to Parkinson''s disease (PD) pathogenesis and underlies the spread of α‐Syn pathology. During this process, extracellular α‐Syn can activate microglia and neuroinflammation, which plays an important role in PD. However, the effect of extracellular α‐Syn on microglia autophagy is poorly understood. In the present study, we reported that extracellular α‐Syn inhibited the autophagy initiation, as indicated by LC3‐II reduction and p62 protein elevation in BV2 and cultured primary microglia. The in vitro findings were verified in microglia‐enriched population isolated from α‐Syn‐overexpressing mice induced by adeno‐associated virus (AAV2/9)‐encoded wildtype human α‐Syn injection into the substantia nigra (SN). Mechanistically, α‐Syn led to microglial autophagic impairment through activating toll‐like receptor 4 (Tlr4) and its downstream p38 and Akt‐mTOR signaling because Tlr4 knockout and inhibition of p38, Akt as well as mTOR prevented α‐Syn‐induced autophagy inhibition. Moreover, inhibition of Akt reversed the mTOR activation but failed to affect p38 phosphorylation triggered by α‐Syn. Functionally, the in vivo evidence showed that lysozyme 2 Cre (Lyz2 ^cre)‐mediated depletion of autophagy‐related gene 5 (Atg5) in microglia aggravated the neuroinflammation and dopaminergic neuron losses in the SN and exacerbated the locomotor deficit in α‐Syn‐overexpressing mice. Taken together, the results suggest that extracellular α‐Syn, via Tlr4‐dependent p38 and Akt‐mTOR signaling cascades, disrupts microglial autophagy activity which synergistically contributes to neuroinflammation and PD development.     

Molecular basis of cross‐species ACE2 interactions with SARS‐CoV‐2‐like viruses of pangolin origin

Correction to: The EMBO Journal (2021) 40: e107786. DOI 10.15252/embj.2021107786 | Published online 8 June 2021The authors would like to add three references to the paper: Starr et al and Zahradník et al also reported that the Q498H or Q498R mutation has enhanced binding affinity to ACE2; and Liu et al reported on the binding of bat coronavirus to ACE2.Starr et al and Zahradník et al have now been cited in the Discussion section, and the following sentence has been corrected from:“According to our data, the SARS‐CoV‐2 RBD with Q498H increases the binding strength to hACE2 by 5‐fold, suggesting the Q498H mutant is more ready to interact with human receptor than the wildtype and highlighting the necessity for more strict control of virus and virus‐infected animals”.to“Here, according to our data and two recently published papers, the SARS‐CoV‐2 RBD with Q498H or Q498R increases the binding strength to hACE2 (Starr et al, 2020; Zahradník et al, 2021), suggesting the mutant with Q498H or Q498R is more ready to interact with human receptor than the wild type and highlighting the necessity for more strict control of virus and virus‐infected animals”.The Liu et al citation has been added to the following sentence:“In another paper published by our group recently, RaTG13 RBD was found to bind to hACE2 with much lower binding affinity than SARS‐CoV‐2 though RaTG13 displays the highest whole‐genome sequence identity (96.2%) with the SARS‐CoV‐2 (Liu et al, 2021)”.Additionally, the authors have added the GISAID accession IDs to the sequence names of the SARS‐CoV‐2 in two human samples (Discussion section). To make identification unambiguous, the sequence names have been updated from “SA‐lsf‐27 and SA‐lsf‐37” to “GISAID accession ID: EPI_ISL_672581 and EPI_ISL_672589”.Lastly, the authors declare in the Materials and Methods section that all experiments employed SARS‐CoV‐2 pseudovirus in cultured cells. These experiments were performed in a BSL‐2‐level laboratory and approved by Science and Technology Conditions Platform Office, Institute of Microbiology, Chinese Academy of Sciences.These changes are herewith incorporated into the paper.