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.     

CellRegMap: a statistical framework for mapping context‐specific regulatory variants using scRNA‐seq

Single‐cell RNA sequencing (scRNA‐seq) enables characterizing the cellular heterogeneity in human tissues. Recent technological advances have enabled the first population‐scale scRNA‐seq studies in hundreds of individuals, allowing to assay genetic effects with single‐cell resolution. However, existing strategies to analyze these data remain based on principles established for the genetic analysis of bulk RNA‐seq. In particular, current methods depend on a priori definitions of discrete cell types, and hence cannot assess allelic effects across subtle cell types and cell states. To address this, we propose the Cell Regulatory Map (CellRegMap), a statistical framework to test for and quantify genetic effects on gene expression in individual cells. CellRegMap provides a principled approach to identify and characterize genotype–context interactions of known eQTL variants using scRNA‐seq data. This model‐based approach resolves allelic effects across cellular contexts of different granularity, including genetic effects specific to cell subtypes and continuous cell transitions. We validate CellRegMap using simulated data and apply it to previously identified eQTL from two recent studies of differentiating iPSCs, where we uncover hundreds of eQTL displaying heterogeneity of genetic effects across cellular contexts. Finally, we identify fine‐grained genetic regulation in neuronal subtypes for eQTL that are colocalized with human disease variants.     

Reduced genetic diversity associated with the northern expansion of an amphibian species with high habitat specialization,Ascaphus truei,resolved using two types of genetic markers

Reconstruction of historical relationships between geographic regions within a species’ range can indicate dispersal patterns and help predict future responses to shifts in climate. Ascaphus truei (coastal tailed frog) is an indicator species of the health of forests and perennial streams in the Coastal and Cascade Mountains of the Pacific Northwest of North America. We used two genetic techniques—microsatellite and genotype‐by‐sequencing (GBS)—to compare the within‐region genetic diversity of populations near the northern extent of the species’ range (British Columbia, Canada) to two geographic regions in British Columbia and two in Washington, USA, moving toward the core of the range. Allelic richness and heterozygosity declined substantially as latitude increased. The northernmost region had the lowest mean expected heterozygosities for both techniques (microsatellite, M = 0.20, SE = 0.080; GBS, M = 0.025, SE = 0.0010) and the southernmost region had the highest (microsatellite, M = 0.88, SE = 0.054; GBS, M = 0.20, SE = 0.0029). The northernmost regions (NC and MC) clustered together in population structure models for both genetic techniques. Our discovery of reduced diversity may have important conservation and management implications for population connectivity and the response of A. truei to climate change.     

Development of in vivo HDX‐MS with applications to a TonB‐dependent transporter and other proteins

Hydrogen‐deuterium exchange mass spectrometry (HDX‐MS) is a powerful tool that monitors protein dynamics in solution. However, the reversible nature of HDX labels has largely limited the application to in vitro systems. Here, we describe a protocol for measuring HDX‐MS in living Escherichia coli cells applied to BtuB, a TonB‐dependent transporter found in outer membranes (OMs). BtuB is a convenient and biologically interesting system for testing in vivo HDX‐MS due to its controllable HDX behavior and large structural rearrangements that occur during the B12 transport cycle. Our previous HDX‐MS study in native OMs provided evidence for B12 binding and breaking of a salt bridge termed the Ionic Lock, an event that leads to the unfolding of the amino terminus. Although purified OMs provide a more native‐like environment than reconstituted systems, disruption of the cell envelope during lysis perturbs the linkage between BtuB and the TonB complex that drives B12 transport. The in vivo HDX response of BtuB''s plug domain (BtuBp) to B12 binding corroborates our previous in vitro findings that B12 alone is sufficient to break the Ionic Lock. In addition, we still find no evidence of B12 binding‐induced unfolding in other regions of BtuBp that could enable B12 passage. Our protocol was successful in reporting on the HDX of several endogenous E. coli proteins measured in the same measurement. Our success in performing HDX in live cells opens the possibility for future HDX‐MS studies in a native cellular environment.IMPORTANCEWe present a protocol for performing in vivo HDX‐MS, focusing on BtuB, a protein whose native membrane environment is believed to be mechanistically important for B12 transport. The in vivo HDX‐MS data corroborate the conclusions from our previous in vitro HDX‐MS study of the allostery initiated by B12 binding. Our success with BtuB and other proteins opens the possibility for performing additional HDX‐MS studies in a native cellular environment.     

Genetic prediction of ICU hospitalization and mortality in COVID‐19 patients using artificial neural networks

There is an unmet need of models for early prediction of morbidity and mortality of Coronavirus disease‐19 (COVID‐19). We aimed to a) identify complement‐related genetic variants associated with the clinical outcomes of ICU hospitalization and death, b) develop an artificial neural network (ANN) predicting these outcomes and c) validate whether complement‐related variants are associated with an impaired complement phenotype. We prospectively recruited consecutive adult patients of Caucasian origin, hospitalized due to COVID‐19. Through targeted next‐generation sequencing, we identified variants in complement factor H/CFH, CFB, CFH‐related, CFD, CD55, C3, C5, CFI, CD46, thrombomodulin/THBD, and A Disintegrin and Metalloproteinase with Thrombospondin motifs (ADAMTS13). Among 381 variants in 133 patients, we identified 5 critical variants associated with severe COVID‐19: rs2547438 (C3), rs2250656 (C3), rs1042580 (THBD), rs800292 (CFH) and rs414628 (CFHR1). Using age, gender and presence or absence of each variant, we developed an ANN predicting morbidity and mortality in 89.47% of the examined population. Furthermore, THBD and C3a levels were significantly increased in severe COVID‐19 patients and those harbouring relevant variants. Thus, we reveal for the first time an ANN accurately predicting ICU hospitalization and death in COVID‐19 patients, based on genetic variants in complement genes, age and gender. Importantly, we confirm that genetic dysregulation is associated with impaired complement phenotype.     

m6A‐mediated alternative splicing coupled with nonsense‐mediated mRNA decay regulates SAM synthetase homeostasis

Alternative splicing of pre‐mRNAs can regulate gene expression levels by coupling with nonsense‐mediated mRNA decay (NMD). In order to elucidate a repertoire of mRNAs regulated by alternative splicing coupled with NMD (AS‐NMD) in an organism, we performed long‐read RNA sequencing of poly(A)⁺ RNAs from an NMD‐deficient mutant strain of Caenorhabditis elegans, and obtained full‐length sequences for mRNA isoforms from 259 high‐confidence AS‐NMD genes. Among them are the S‐adenosyl‐L‐methionine (SAM) synthetase (sams) genes sams‐3 and sams‐4. SAM synthetase activity autoregulates sams gene expression through AS‐NMD in a negative feedback loop. We furthermore find that METT‐10, the orthologue of human U6 snRNA methyltransferase METTL16, is required for the splicing regulation in␣vivo, and specifically methylates the invariant AG dinucleotide at the distal 3′ splice site (3′SS) in␣vitro. Direct RNA sequencing coupled with machine learning confirms m⁶A modification of endogenous sams mRNAs. Overall, these results indicate that homeostasis of SAM synthetase in C. elegans is maintained by alternative splicing regulation through m⁶A modification at the 3′SS of the sams genes.     

Attenuation of SARS‐CoV‐2 replication and associated inflammation by concomitant targeting of viral and host cap 2'‐O‐ribose methyltransferases

The SARS‐CoV‐2 infection cycle is a multistage process that relies on functional interactions between the host and the pathogen. Here, we repurposed antiviral drugs against both viral and host enzymes to pharmaceutically block methylation of the viral RNA 2''‐O‐ribose cap needed for viral immune escape. We find that the host cap 2''‐O‐ribose methyltransferase MTr1 can compensate for loss of viral NSP16 methyltransferase in facilitating virus replication. Concomitant inhibition of MTr1 and NSP16 efficiently suppresses SARS‐CoV‐2 replication. Using in silico target‐based drug screening, we identify a bispecific MTr1/NSP16 inhibitor with anti‐SARS‐CoV‐2 activity in vitro and in vivo but with unfavorable side effects. We further show antiviral activity of inhibitors that target independent stages of the host SAM cycle providing the methyltransferase co‐substrate. In particular, the adenosylhomocysteinase (AHCY) inhibitor DZNep is antiviral in in vitro, in ex vivo, and in a mouse infection model and synergizes with existing COVID‐19 treatments. Moreover, DZNep exhibits a strong immunomodulatory effect curbing infection‐induced hyperinflammation and reduces lung fibrosis markers ex vivo. Thus, multispecific and metabolic MTase inhibitors constitute yet unexplored treatment options against COVID‐19.     

N‐acetylaspartate release by glutaminolytic ovarian cancer cells sustains protumoral macrophages

Glutaminolysis is known to correlate with ovarian cancer aggressiveness and invasion. However, how this affects the tumor microenvironment is elusive. Here, we show that ovarian cancer cells become addicted to extracellular glutamine when silenced for glutamine synthetase (GS), similar to naturally occurring GS‐low, glutaminolysis‐high ovarian cancer cells. Glutamine addiction elicits a crosstalk mechanism whereby cancer cells release N‐acetylaspartate (NAA) which, through the inhibition of the NMDA receptor, and synergistically with IL‐10, enforces GS expression in macrophages. In turn, GS‐high macrophages acquire M2‐like, tumorigenic features. Supporting this in␣vitro model, in silico data and the analysis of ascitic fluid isolated from ovarian cancer patients prove that an M2‐like macrophage phenotype, IL‐10 release, and NAA levels positively correlate with disease stage. Our study uncovers the unprecedented role of glutamine metabolism in modulating macrophage polarization in highly invasive ovarian cancer and highlights the anti‐inflammatory, protumoral function of NAA.     

Nonproductive exposure of PBMCs to SARS‐CoV‐2 induces cell‐intrinsic innate immune responses

Cell‐intrinsic responses mounted in PBMCs during mild and severe COVID‐19 differ quantitatively and qualitatively. Whether they are triggered by signals emitted by productively infected cells of the respiratory tract or result from physical interaction with virus particles remains unclear. Here, we analyzed susceptibility and expression profiles of PBMCs from healthy donors upon ex vivo exposure to SARS‐CoV and SARS‐CoV‐2. In line with the absence of detectable ACE2 receptor expression, human PBMCs were refractory to productive infection. RT–PCR experiments and single‐cell RNA sequencing revealed JAK/STAT‐dependent induction of interferon‐stimulated genes (ISGs) but not proinflammatory cytokines. This SARS‐CoV‐2‐specific response was most pronounced in monocytes. SARS‐CoV‐2‐RNA‐positive monocytes displayed a lower ISG signature as compared to bystander cells of the identical culture. This suggests a preferential invasion of cells with a low ISG baseline profile or delivery of a SARS‐CoV‐2‐specific sensing antagonist upon efficient particle internalization. Together, nonproductive physical interaction of PBMCs with SARS‐CoV‐2‐ and, to a much lesser extent, SARS‐CoV particles stimulate JAK/STAT‐dependent, monocyte‐accentuated innate immune responses that resemble those detected in vivo in patients with mild COVID‐19.     

A functional assay for serum detection of antibodies against SARS‐CoV‐2 nucleoprotein

The humoral immune response to SARS‐CoV‐2 results in antibodies against spike (S) and nucleoprotein (N). However, whilst there are widely available neutralization assays for S antibodies, there is no assay for N‐antibody activity. Here, we present a simple in vitro method called EDNA (electroporated‐antibody‐dependent neutralization assay) that provides a quantitative measure of N‐antibody activity in unpurified serum from SARS‐CoV‐2 convalescents. We show that N antibodies neutralize SARS‐CoV‐2 intracellularly and cell‐autonomously but require the cytosolic Fc receptor TRIM21. Using EDNA, we show that low N‐antibody titres can be neutralizing, whilst some convalescents possess serum with high titres but weak activity. N‐antibody and N‐specific T‐cell activity correlates within individuals, suggesting N antibodies may protect against SARS‐CoV‐2 by promoting antigen presentation. This work highlights the potential benefits of N‐based vaccines and provides an in vitro assay to allow the antibodies they induce to be tested.     

Loss of enteric neuronal Ndrg4 promotes colorectal cancer via increased release of Nid1 and Fbln2

The N‐Myc Downstream‐Regulated Gene 4 (NDRG4), a prominent biomarker for colorectal cancer (CRC), is specifically expressed by enteric neurons. Considering that nerves are important members of the tumor microenvironment, we here establish different Ndrg4 knockout (Ndrg4 ^−/−) CRC models and an indirect co‐culture of primary enteric nervous system (ENS) cells and intestinal organoids to identify whether the ENS, via NDRG4, affects intestinal tumorigenesis. Linking immunostainings and gastrointestinal motility (GI) assays, we show that the absence of Ndrg4 does not trigger any functional or morphological GI abnormalities. However, combining in vivo, in vitro, and quantitative proteomics data, we uncover that Ndrg4 knockdown is associated with enlarged intestinal adenoma development and that organoid growth is boosted by the Ndrg4 ^−/− ENS cell secretome, which is enriched for Nidogen‐1 (Nid1) and Fibulin‐2 (Fbln2). Moreover, NID1 and FBLN2 are expressed in enteric neurons, enhance migration capacities of CRC cells, and are enriched in human CRC secretomes. Hence, we provide evidence that the ENS, via loss of Ndrg4, is involved in colorectal pathogenesis and that ENS‐derived Nidogen‐1 and Fibulin‐2 enhance colorectal carcinogenesis.     

MicroRNA‐495‐3p diminishes doxorubicin‐induced cardiotoxicity through activating AKT

Doxorubicin (Dox) is a broad‐spectrum antitumour agent; however, its clinical application is impeded due to the cumulative cardiotoxicity. The present study aims to investigate the role and underlying mechanisms of microRNA‐495‐3p (miR‐495‐3p) in Dox‐induced cardiotoxicity. Herein, we found that cardiac miR‐495‐3p expression was significantly decreased in Dox‐treated hearts, and that the miR‐495‐3p agomir could prevent oxidative stress, cell apoptosis, cardiac mass loss, fibrosis and cardiac dysfunction upon Dox stimulation. In contrast, the miR‐495‐3p antagomir dramatically aggravated Dox‐induced cardiotoxicity in mice. Besides, we found that the miR‐495‐3p agomir attenuated, while the miR‐495‐3p antagomir exacerbated Dox‐induced oxidative stress and cellular injury in vitro. Mechanistically, we demonstrated that miR‐495‐3p directly bound to the 3′‐untranslational region of phosphate and tension homology deleted on chromosome ten (PTEN), downregulated PTEN expression and subsequently activated protein kinase B (PKB/AKT) pathway, and that PTEN overexpression or AKT inhibition completely abolished the cardioprotective effects of the miR‐495‐3p agomir. Our study for the first time identify miR‐495‐3p as an endogenous protectant against Dox‐induced cardiotoxicity through activating AKT pathway in vivo and in vitro.     

Cyclic di‐GMP modulates sessile‐motile phenotypes and virulence in Dickeya oryzae via two PilZ domain receptors

Dickeya oryzae is a bacterial pathogen causing the severe rice stem rot disease in China and other rice‐growing countries. We showed recently that the universal bacterial second messenger c‐di‐GMP plays an important role in modulation of bacterial motility and pathogenicity, but the mechanism of regulation remains unknown. In this study, bioinformatics analysis of the D. oryzae EC1 genome led to the identification of two proteins, YcgR and BcsA, both of which contain a conserved c‐di‐GMP receptor domain, known as the PilZ‐domain. By deleting all the genes encoding c‐di‐GMP‐degrading enzymes in D. oryzae EC1, the resultant mutant 7ΔPDE with high c‐di‐GMP levels became nonmotile, formed hyperbiofilm, and lost the ability to colonize and invade rice seeds. These phenotypes were partially reversed by deletion of ycgR in the mutant 7ΔPDE, whereas deletion of bcsA only reversed the hyperbiofilm phenotype of mutant 7ΔPDE. Significantly, double deletion of ycgR and bcsA in mutant 7ΔPDE rescued its motility, biofilm formation, and virulence to levels of wild‐type EC1. In vitro biochemical experiments and in vivo phenotypic assays further validated that YcgR and BcsA proteins are the receptors for c‐di‐GMP, which together play a critical role in regulating the c‐di‐GMP‐associated functionality. The findings from this study fill a gap in our understanding of how c‐di‐GMP modulates bacterial motility and biofilm formation, and provide useful clues for further elucidation of sophisticated virulence regulatory mechanisms in this important plant pathogen.     

Hypophosphorylated pRb knock‐in mice exhibit hallmarks of aging and vitamin C‐preventable diabetes

Despite extensive analysis of pRB phosphorylation in vitro, how this modification influences development and homeostasis in vivo is unclear. Here, we show that homozygous Rb^∆K4 and Rb^∆K7 knock‐in mice, in which either four or all seven phosphorylation sites in the C‐terminal region of pRb, respectively, have been abolished by Ser/Thr‐to‐Ala substitutions, undergo normal embryogenesis and early development, notwithstanding suppressed phosphorylation of additional upstream sites. Whereas Rb^∆K4 mice exhibit telomere attrition but no other abnormalities, Rb^∆K7 mice are smaller and display additional hallmarks of premature aging including infertility, kyphosis, and diabetes, indicating an accumulative effect of blocking pRb phosphorylation. Diabetes in Rb^∆K7 mice is insulin‐sensitive and associated with failure of quiescent pancreatic β‐cells to re‐enter the cell cycle in response to mitogens, resulting in induction of DNA damage response (DDR), senescence‐associated secretory phenotype (SASP), and reduced pancreatic islet mass and circulating insulin level. Pre‐treatment with the epigenetic regulator vitamin C reduces DDR, increases cell cycle re‐entry, improves islet morphology, and attenuates diabetes. These results have direct implications for cell cycle regulation, CDK‐inhibitor therapeutics, diabetes, and longevity.     

DNA metabarcoding quantifies the relative biomass of arthropod taxa in songbird diets: Validation with camera‐recorded diets

Ecological research is often hampered by the inability to quantify animal diets. Diet composition can be tracked through DNA metabarcoding of fecal samples, but whether (complex) diets can be quantitatively determined with metabarcoding is still debated and needs validation using free‐living animals. This study validates that DNA metabarcoding of feces can retrieve actual ingested taxa, and most importantly, that read numbers retrieved from sequencing can also be used to quantify the relative biomass of dietary taxa. Validation was done with the hole‐nesting insectivorous Pied Flycatcher whose diet was quantified using camera footage. Size‐adjusted counts of food items delivered to nestlings were used as a proxy for provided biomass of prey orders and families, and subsequently, nestling feces were assessed through DNA metabarcoding. To explore potential effects of digestion, gizzard and lower intestine samples of freshly collected birds were subjected to DNA metabarcoding. For metabarcoding with Cytochrome Oxidase subunit I (COI), we modified published invertebrate COI primers LCO1490 and HCO1777, which reduced host reads to 0.03%, and amplified Arachnida DNA without significant changing the recovery of other arthropod taxa. DNA metabarcoding retrieved all commonly camera‐recorded taxa. Overall, and in each replicate year (N = 3), the relative scaled biomass of prey taxa and COI read numbers correlated at R = .85 (95CI:0.68–0.94) at order level and at R = .75 (CI:0.67–0.82) at family level. Similarity in arthropod community composition between gizzard and intestines suggested limited digestive bias. This DNA metabarcoding validation demonstrates that quantitative analyses of arthropod diet is possible. We discuss the ecological applications for insectivorous birds.     

A non‐helical region in transmembrane helix 6 of hydrophobic amino acid transporter MhsT mediates substrate recognition

MhsT of Bacillus halodurans is a transporter of hydrophobic amino acids and a homologue of the eukaryotic SLC6 family of Na⁺‐dependent symporters for amino acids, neurotransmitters, osmolytes, or creatine. The broad range of transported amino acids by MhsT prompted the investigation of the substrate recognition mechanism. Here, we report six new substrate‐bound structures of MhsT, which, in conjunction with functional studies, reveal how the flexibility of a Gly‐Met‐Gly (GMG) motif in the unwound region of transmembrane segment 6 (TM6) is central for the recognition of substrates of different size by tailoring the binding site shape and volume. MhsT mutants, harboring substitutions within the unwound GMG loop and substrate binding pocket that mimick the binding sites of eukaryotic SLC6A18/B0AT3 and SLC6A19/B0AT1 transporters of neutral amino acids, exhibited impaired transport of aromatic amino acids that require a large binding site volume. Conservation of a general (G/A/C)ΦG motif among eukaryotic members of SLC6 family suggests a role for this loop in a common mechanism for substrate recognition and translocation by SLC6 transporters of broad substrate specificity.     

Microglial MERTK eliminates phosphatidylserine‐displaying inhibitory post‐synapses

Glia contribute to synapse elimination through phagocytosis in the central nervous system. Despite the important roles of this process in development and neurological disorders, the identity and regulation of the "eat‐me" signal that initiates glia‐mediated phagocytosis of synapses has remained incompletely understood. Here, we generated conditional knockout mice with neuronal‐specific deletion of the flippase chaperone Cdc50a, to induce stable exposure of phosphatidylserine, a well‐known "eat‐me" signal for apoptotic cells, on the neuronal outer membrane. Surprisingly, acute Cdc50a deletion in mature neurons causes preferential phosphatidylserine exposure in neuronal somas and specific loss of inhibitory post‐synapses without effects on other synapses, resulting in abnormal excitability and seizures. Ablation of microglia or the deletion of microglial phagocytic receptor Mertk prevents the loss of inhibitory post‐synapses and the seizure phenotype, indicating that microglial phagocytosis is responsible for inhibitory post‐synapse elimination. Moreover, we found that phosphatidylserine is used for microglia‐mediated pruning of inhibitory post‐synapses in normal brains, suggesting that phosphatidylserine serves as a general "eat‐me" signal for inhibitory post‐synapse elimination.