Genome‐scale metabolic modeling reveals SARS‐CoV‐2‐induced metabolic changes and antiviral targets

Tremendous progress has been made to control the COVID‐19 pandemic caused by the SARS‐CoV‐2 virus. However, effective therapeutic options are still rare. Drug repurposing and combination represent practical strategies to address this urgent unmet medical need. Viruses, including coronaviruses, are known to hijack host metabolism to facilitate viral proliferation, making targeting host metabolism a promising antiviral approach. Here, we describe an integrated analysis of 12 published in vitro and human patient gene expression datasets on SARS‐CoV‐2 infection using genome‐scale metabolic modeling (GEM), revealing complicated host metabolism reprogramming during SARS‐CoV‐2 infection. We next applied the GEM‐based metabolic transformation algorithm to predict anti‐SARS‐CoV‐2 targets that counteract the virus‐induced metabolic changes. We successfully validated these targets using published drug and genetic screen data and by performing an siRNA assay in Caco‐2 cells. Further generating and analyzing RNA‐sequencing data of remdesivir‐treated Vero E6 cell samples, we predicted metabolic targets acting in combination with remdesivir, an approved anti‐SARS‐CoV‐2 drug. Our study provides clinical data‐supported candidate anti‐SARS‐CoV‐2 targets for future evaluation, demonstrating host metabolism targeting as a promising antiviral strategy.     

Yeast metabolic innovations emerged via expanded metabolic network and gene positive selection

Yeasts are known to have versatile metabolic traits, while how these metabolic traits have evolved has not been elucidated systematically. We performed integrative evolution analysis to investigate how genomic evolution determines trait generation by reconstructing genome‐scale metabolic models (GEMs) for 332 yeasts. These GEMs could comprehensively characterize trait diversity and predict enzyme functionality, thereby signifying that sequence‐level evolution has shaped reaction networks towards new metabolic functions. Strikingly, using GEMs, we can mechanistically map different evolutionary events, e.g. horizontal gene transfer and gene duplication, onto relevant subpathways to explain metabolic plasticity. This demonstrates that gene family expansion and enzyme promiscuity are prominent mechanisms for metabolic trait gains, while GEM simulations reveal that additional factors, such as gene loss from distant pathways, contribute to trait losses. Furthermore, our analysis could pinpoint to specific genes and pathways that have been under positive selection and relevant for the formulation of complex metabolic traits, i.e. thermotolerance and the Crabtree effect. Our findings illustrate how multidimensional evolution in both metabolic network structure and individual enzymes drives phenotypic variations.     

Genome‐wide selective detection of Mile red‐bone goat using next‐generation sequencing technology

The ecotype population of goats (Capra hircus) was created by long‐term artificial selection and natural adaptation. Mile red‐bone goat is an indigenous breed with visible red bones, and its special bone structure has received extensive attention. This study aimed to identify genetic variants and candidate genes associated with specific bone phenotypes using next‐generation sequencing technology (NGS). The results revealed that 31,828,206 single nucleotide polymorphisms (SNPs) were obtained from 72 goats (20 Mile red‐bone goats and 52 common goats) by NGS. A total of 100 candidate genes were identified on the basis top 1% window interaction from nucleotide diversity (π), π ratio (π _A/π _B), and pairwise fixation index (F _ST). Exactly 77 known signaling pathways were enriched. Specifically, three coding genes (NMNAT2, LOC102172983, and PNLIP) were annotated in the vitamin metabolism signaling pathways, and NCF2 was annotated to the osteoclast (OC) differentiation pathway. Furthermore, 5862 reliable copy number variations (CNVs) were obtained, and 14 and 24 genes were annotated with the top 1‰ CNV based on F _ST (>0.490) and V _ST (>0.527), respectively. Several pathways related to bone development and metabolism of exogenous substances in vivo, including calcium signaling pathway, OC differentiation, and glycerophospholipid metabolism, were annotated. Specifically, six genes from 19 candidate CNVs, which were obtained by interaction of the top 1‰ CNVs with F _ST and V _ST, were annotated to mucin‐type O‐glycan biosynthesis and metabolic pathways. Briefly, the results implied that pseudopurpurin and specific genetic variants work together to contribute to the red‐bone color and specific bone structure of Mile red‐bone goat. This study is helpful to understanding the genetic basis of the unique bone phenotype of Mile red‐bone goats.     

Genome‐scale gene/reaction essentiality and synthetic lethality analysis

Synthetic lethals are to pairs of non‐essential genes whose simultaneous deletion prohibits growth. One can extend the concept of synthetic lethality by considering gene groups of increasing size where only the simultaneous elimination of all genes is lethal, whereas individual gene deletions are not. We developed optimization‐based procedures for the exhaustive and targeted enumeration of multi‐gene (and by extension multi‐reaction) lethals for genome‐scale metabolic models. Specifically, these approaches are applied to iAF1260, the latest model of Escherichia coli, leading to the complete identification of all double and triple gene and reaction synthetic lethals as well as the targeted identification of quadruples and some higher‐order ones. Graph representations of these synthetic lethals reveal a variety of motifs ranging from hub‐like to highly connected subgraphs providing a birds‐eye view of the avenues available for redirecting metabolism and uncovering complex patterns of gene utilization and interdependence. The procedure also enables the use of falsely predicted synthetic lethals for metabolic model curation. By analyzing the functional classifications of the genes involved in synthetic lethals, we reveal surprising connections within and across clusters of orthologous group functional classifications.     

SARS‐CoV‐2 structural coverage map reveals viral protein assembly,mimicry, and hijacking mechanisms

We modeled 3D structures of all SARS‐CoV‐2 proteins, generating 2,060 models that span 69% of the viral proteome and provide details not available elsewhere. We found that ˜6% of the proteome mimicked human proteins, while ˜7% was implicated in hijacking mechanisms that reverse post‐translational modifications, block host translation, and disable host defenses; a further ˜29% self‐assembled into heteromeric states that provided insight into how the viral replication and translation complex forms. To make these 3D models more accessible, we devised a structural coverage map, a novel visualization method to show what is—and is not—known about the 3D structure of the viral proteome. We integrated the coverage map into an accompanying online resource (https://aquaria.ws/covid) that can be used to find and explore models corresponding to the 79 structural states identified in this work. The resulting Aquaria‐COVID resource helps scientists use emerging structural data to understand the mechanisms underlying coronavirus infection and draws attention to the 31% of the viral proteome that remains structurally unknown or dark.     

FTLD‐TDP assemblies seed neoaggregates with subtype‐specific features via a prion‐like cascade

Morphologically distinct TDP‐43 aggregates occur in clinically different FTLD‐TDP subtypes, yet the mechanism of their emergence and contribution to clinical heterogeneity are poorly understood. Several lines of evidence suggest that pathological TDP‐43 follows a prion‐like cascade, but the molecular determinants of this process remain unknown. We use advanced microscopy techniques to compare the seeding properties of pathological FTLD‐TDP‐A and FTLD‐TDP‐C aggregates. Upon inoculation of patient‐derived aggregates in cells, FTLD‐TDP‐A seeds amplify in a template‐dependent fashion, triggering neoaggregation more efficiently than those extracted from FTLD‐TDP‐C patients, correlating with the respective disease progression rates. Neoaggregates are sequentially phosphorylated with N‐to‐C directionality and with subtype‐specific timelines. The resulting FTLD‐TDP‐A neoaggregates are large and contain densely packed fibrils, reminiscent of the pure compacted fibrils present within cytoplasmic inclusions in postmortem brains. In contrast, FTLD‐TDP‐C dystrophic neurites show less dense fibrils mixed with cellular components, and their respective neoaggregates are small, amorphous protein accumulations. These cellular seeding models replicate aspects of the patient pathological diversity and will be a useful tool in the quest for subtype‐specific therapeutics.     

The atomic portrait of SARS‐CoV‐2 as captured by cryo‐electron microscopy

Transmission electron microscopy has historically been indispensable for virology research, as it offers unique insight into virus function. In the past decade, as cryo‐electron microscopy (cryo‐EM) has matured and become more accessible, we have been able to peer into the structure of viruses at the atomic level and understand how they interact with the host cell, with drugs or with antibodies. Perhaps, there was no time in recent history where cryo‐EM was more needed, as SARS‐CoV‐2 has spread around the globe, causing millions of deaths and almost unquantifiable economic devastation. In this concise review, we aim to mark the most important contributions of cryo‐EM to understanding the structure and function of SARS‐CoV‐2 proteins, from surface spikes to the virus core and from virus‐receptor interactions to antibody binding.     

High‐efficiency delivery of CRISPR‐Cas9 by engineered probiotics enables precise microbiome editing

Antibiotic resistance threatens our ability to treat infectious diseases, spurring interest in alternative antimicrobial technologies. The use of bacterial conjugation to deliver CRISPR‐cas systems programmed to precisely eliminate antibiotic‐resistant bacteria represents a promising approach but requires high in situ DNA transfer rates. We have optimized the transfer efficiency of conjugative plasmid TP114 using accelerated laboratory evolution. We hence generated a potent conjugative delivery vehicle for CRISPR‐cas9 that can eliminate > 99.9% of targeted antibiotic‐resistant Escherichia coli in the mouse gut microbiota using a single dose. We then applied this system to a Citrobacter rodentium infection model, achieving full clearance within four consecutive days of treatment.     

细菌最小基因组研究进展

具有最小基因组的细菌只包含维持自我生命复制所必需的基因,其作为一种潜在的工业生产平台具有诸多优势。由于高通量DNA测序和合成技术的发展,目前已经构建了多种缩减基因组的菌株。本文首先介绍了最小基因组的概念,其次总结了细菌必需基因的相关研究进展,然后梳理了人工缩减与合成微生物基因组的相关工作,最后探讨了在设计和组装基因组的过程中遇到的技术障碍和限制,以期为人工合成基因组的实验与应用提供理论参考。     

Genome‐wide analysis reveals demographic and life‐history patterns associated with habitat modification in landlocked,deep‐spawning sockeye salmon (Oncorhynchus nerka)

Human‐mediated habitat fragmentation in freshwater ecosystems can negatively impact genetic diversity, demography, and life history of native biota, while disrupting the behavior of species that are dependent on spatial connectivity to complete their life cycles. In the Alouette River system (British Columbia, Canada), dam construction in 1928 impacted passage of anadromous sockeye salmon (Oncorhynchus nerka), with the last records of migrants occurring in the 1930s. Since that time, O. nerka persisted as a resident population in Alouette Reservoir until experimental water releases beginning in 2005 created conditions for migration; two years later, returning migrants were observed for the first time in ~70 years, raising important basic and applied questions regarding life‐history variation and population structure in this system. Here, we investigated the genetic distinctiveness and population history of Alouette Reservoir O. nerka using genome‐wide SNP data (n = 7,709 loci) collected for resident and migrant individuals, as well as for neighboring anadromous sockeye salmon and resident kokanee populations within the Fraser River drainage (n = 312 individuals). Bayesian clustering and principal components analyses based on neutral loci revealed five distinct clusters, largely associated with geography, and clearly demonstrated that Alouette Reservoir resident and migrant individuals are genetically distinct from other O. nerka populations in the Fraser River drainage. At a finer level, there was no clear evidence for differentiation between Alouette Reservoir residents and migrants; although we detected eight high‐confidence outlier loci, they all mapped to sex chromosomes suggesting that differences were likely due to uneven sex ratios rather than life history. Taken together, these data suggest that contemporary Alouette Reservoir O. nerka represents a landlocked sockeye salmon population, constituting the first reported instance of deep‐water spawning behavior associated with this life‐history form. This finding punctuates the need for reassessment of conservation status and supports ongoing fisheries management activities in Alouette Reservoir.     

Micro‐CT reconstruction reveals the colony pattern regulations of four dominant reef‐building corals

Colonies are the basic geometric building blocks of coral reefs. However, the forming regulations of both colonies and reefs are still not understood adequately. Therefore, in this study, we reconstructed 25 samples using high‐resolution micro‐computed tomography to investigate coral growth patterns and parameters. Our skeleton and canal reconstructions revealed the characteristics of different coral species, and we further visualized the growth axes and growth rings to understand the coral growth directions. We drew a skeleton grayscale map and calculated the coral skeleton void ratios to ascertain the skeletal diversity, devising a method to quantify coral growth. On the basis of the three‐dimensional (3D) reconstructions and growth parameters, we investigated the growth strategies of different coral species. This research increases the breadth of knowledge on how reef‐building corals grow their colonies, providing information on reef‐forming regulations. The data in this paper contain a large amount of coral growth information, which can be used in further research on reef‐forming patterns under different conditions. The method used in this study can also be applied to animals with porous skeletons.     

Sex‐ and strain‐specific effects of mitochondrial uncoupling on age‐related metabolic diseases in high‐fat diet‐fed mice

Mild uncoupling of oxidative phosphorylation is an intrinsic property of all mitochondria and may have evolved to protect cells against the production of damaging reactive oxygen species. Therefore, compounds that enhance mitochondrial uncoupling are potentially attractive anti‐aging therapies; however, chronic ingestion is associated with a number of unwanted side effects. We have previously developed a controlled‐release mitochondrial protonophore (CRMP) that is functionally liver‐directed and promotes oxidation of hepatic triglycerides by causing a subtle sustained increase in hepatic mitochondrial inefficiency. Here, we sought to leverage the higher therapeutic index of CRMP to test whether mild mitochondrial uncoupling in a liver‐directed fashion could reduce oxidative damage and improve age‐related metabolic disease and lifespan in diet‐induced obese mice. Oral administration of CRMP (20 mg/[kg‐day] × 4 weeks) reduced hepatic lipid content, protein kinase C epsilon activation, and hepatic insulin resistance in aged (74‐week‐old) high‐fat diet (HFD)‐fed C57BL/6J male mice, independently of changes in body weight, whole‐body energy expenditure, food intake, or markers of hepatic mitochondrial biogenesis. CRMP treatment was also associated with a significant reduction in hepatic lipid peroxidation, protein carbonylation, and inflammation. Importantly, long‐term (49 weeks) hepatic mitochondrial uncoupling initiated late in life (94–104 weeks), in conjugation with HFD feeding, protected mice against neoplastic disorders, including hepatocellular carcinoma (HCC), in a strain and sex‐specific manner. Taken together, these studies illustrate the complex variation of aging and provide important proof‐of‐concept data to support further studies investigating the use of liver‐directed mitochondrial uncouplers to promote healthy aging in humans.     

Hydrogen sulphide ameliorates dopamine‐induced astrocytic inflammation and neurodegeneration in minimal hepatic encephalopathy

It has been demonstrated that the action of dopamine (DA) could enhance the production of tumour necrosis factor‐α (TNF‐α) by astrocytes and potentiate neuronal apoptosis in minimal hepatic encephalopathy (MHE). Recently, sodium hydrosulfide (NaHS) has been found to have neuroprotective properties. Our study addressed whether NaHS could rescue DA‐challenged inflammation and apoptosis in neurons to ameliorate memory impairment in MHE rats and in the neuron and astrocyte coculture system. We found that NaHS suppressed DA‐induced p65 acetylation, resulting in reduced TNF‐α production in astrocytes both in vitro and in vivo. Furthermore, decreased apoptosis was observed in neurons exposed to conditioned medium from DA + NaHS‐challenged astrocytes, which was similar to the results obtained in the neurons exposed to TNF‐α + NaHS, suggesting a therapeutic effect of NaHS on the suppression of neuronal apoptosis via the reduction of TNF‐α level. DA triggered the inactivation of p70 S6 ribosomal kinase (S6K1) and dephosphorylation of Bad, resulting in the disaggregation of Bclxl and Bak and the release of cytochrome c (Cyt. c), and this process could be reversed by NaHS administration. Our work demonstrated that NaHS attenuated DA‐induced astrocytic TNF‐α release and ameliorated inflammation‐induced neuronal apoptosis in MHE. Further research into this approach may uncover future potential therapeutic strategies for MHE.     

Diagnostics and correction of batch effects in large‐scale proteomic studies: a tutorial

Advancements in mass spectrometry‐based proteomics have enabled experiments encompassing hundreds of samples. While these large sample sets deliver much‐needed statistical power, handling them introduces technical variability known as batch effects. Here, we present a step‐by‐step protocol for the assessment, normalization, and batch correction of proteomic data. We review established methodologies from related fields and describe solutions specific to proteomic challenges, such as ion intensity drift and missing values in quantitative feature matrices. Finally, we compile a set of techniques that enable control of batch effect adjustment quality. We provide an R package, "proBatch", containing functions required for each step of the protocol. We demonstrate the utility of this methodology on five proteomic datasets each encompassing hundreds of samples and consisting of multiple experimental designs. In conclusion, we provide guidelines and tools to make the extraction of true biological signal from large proteomic studies more robust and transparent, ultimately facilitating reliable and reproducible research in clinical proteomics and systems biology.     

Cryo‐EM reveals mechanisms of angiotensin I‐converting enzyme allostery and dimerization

Hypertension (high blood pressure) is a major risk factor for cardiovascular disease, which is the leading cause of death worldwide. The somatic isoform of angiotensin I‐converting enzyme (sACE) plays a critical role in blood pressure regulation, and ACE inhibitors are thus widely used to treat hypertension and cardiovascular disease. Our current understanding of sACE structure, dynamics, function, and inhibition has been limited because truncated, minimally glycosylated forms of sACE are typically used for X‐ray crystallography and molecular dynamics simulations. Here, we report the first cryo‐EM structures of full‐length, glycosylated, soluble sACE (sACE^S1211). Both monomeric and dimeric forms of the highly flexible apo enzyme were reconstructed from a single dataset. The N‐ and C‐terminal domains of monomeric sACE^S1211 were resolved at 3.7 and 4.1 Å, respectively, while the interacting N‐terminal domains responsible for dimer formation were resolved at 3.8 Å. Mechanisms are proposed for intradomain hinging, cooperativity, and homodimerization. Furthermore, the observation that both domains were in the open conformation has implications for the design of sACE modulators.     

A new family of CRISPR‐type V nucleases with C‐rich PAM recognition

Most CRISPR‐type V nucleases are stimulated to cleave double‐stranded (ds) DNA targets by a T‐rich PAM, which restricts their targeting range. Here, we identify and characterize a new family of type V RNA‐guided nuclease, Cas12l, that exclusively recognizes a C‐rich (5''‐CCY‐3′) PAM. The organization of genes within its CRISPR locus is similar to type II‐B CRISPR‐Cas9 systems, but both sequence analysis and functional studies establish it as a new family of type V effector. Biochemical experiments show that Cas12l nucleases function optimally between 37 and 52°C, depending on the ortholog, and preferentially cut supercoiled DNA. Like other type V nucleases, it exhibits collateral nonspecific ssDNA and ssRNA cleavage activity that is triggered by ssDNA or dsDNA target recognition. Finally, we show that one family member, Asp2Cas12l, functions in a heterologous cellular environment, altogether, suggesting that this new group of CRISPR‐associated nucleases may be harnessed as genome editing reagents.