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.     

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.     

Degrees of compositional shift in tree communities vary along a gradient of temperature change rates over one decade: Application of an individual‐based temporal beta‐diversity concept

Temporal patterns in communities have gained widespread attention recently, to the extent that temporal changes in community composition are now termed “temporal beta‐diversity.” Previous studies of beta‐diversity have made use of two classes of dissimilarity indices: incidence‐based (e.g., Sørensen and Jaccard dissimilarity) and abundance‐based (e.g., Bray–Curtis and Ružička dissimilarity). However, in the context of temporal beta‐diversity, the persistence of identical individuals and turnover among other individuals within the same species over time have not been considered, despite the fact that both will affect compositional changes in communities. To address this issue, I propose new index concepts for beta‐diversity and the relative speed of compositional shifts in relation to individual turnover based on individual identity information. Individual‐based beta‐diversity indices are novel dissimilarity indices that consider individual identity information to quantitatively evaluate temporal change in individual turnover and community composition. I applied these new indices to individually tracked tree monitoring data in deciduous and evergreen broad‐leaved forests across the Japanese archipelago with the objective of quantifying the effect of climate change trends (i.e., rates of change in both annual mean temperature and annual precipitation) on individual turnover and compositional shifts at each site. A new index explored the relative contributions of mortality and recruitment processes to temporal changes in community composition. Clear patterns emerged showing that an increase in the temperature change rate facilitated the relative contribution of mortality components. The relative speed of compositional shift increased with increasing temperature change rates in deciduous forests but decreased with increasing warming rates in evergreen forests. These new concepts provide a way to identify novel and high‐resolution temporal patterns in communities.     

circRNA 001306 enhances hepatocellular carcinoma growth by up‐regulating CDK16 expression via sponging miR‐584‐5p

Circular RNAs (circRNAs) have been demonstrated to play important roles in cancer progress. However, the roles in hepatocellular carcinoma (HCC) are still unclear. Here, we found has_circRNA_001306 (circ_1306) was up‐regulated in HCC tissues and cell lines. Knockdown the expression circ_1306 significantly suppressed HCC cell proliferation and induced the cell apoptosis in vitro and in vivo. Furthermore, we identified circ_1306 could up‐regulate the expression of CDK16 by sponging miR‐584‐5p. The expression of miR‐584‐5p was decreased, and the expression of CDK16 was increased in HCC tissues and cell lines. Meanwhile, either knockdown of miR‐584‐5p or overexpression of CDK16 could suppress the HCC cell proliferation. In vivo, overexpression of miR‐584‐5p or knockdown of circ_1306 could inhibit the expression of CDK16, and suppress tumour growth. Altogether, our findings suggested that circ_1306 could promoter HCC progress by miR‐584‐5p/CDK16 axis, which provided a novel marker for HCC diagnosis and treatment.     

Establishment of a five‐enzalutamide‐resistance‐related‐gene‐based classifier for recurrence‐free survival predicting of prostate cancer

To identify prostate cancer (PCa) patients with a high risk of recurrence is critical before delivering adjuvant treatment. We developed a classifier based on the Enzalutamide treatment resistance‐related genes to assist the currently available staging system in predicting the recurrence‐free survival (RFS) prognosis of PCa patients. We overlapped the DEGs from two datasets to obtain a more convincing Enzalutamide‐resistance‐related‐gene (ERRG) cluster. The five‐ERRG‐based classifier obtained good predictive values in both the training and validation cohorts. The classifier precisely predicted RFS of patients in four cohorts, independent of patient age, pathological tumour stage, Gleason score and PSA levels. The classifier and the clinicopathological factors were combined to construct a nomogram, which had an increased predictive accuracy than that of each variable alone. Besides, we also compared the differences between high‐ and low‐risk subgroups and found their differences were enriched in cancer progression‐related pathways. The five‐ERRG‐based classifier is a practical and reliable predictor, which adds value to the existing staging system for predicting the RFS prognosis of PCa after radical prostatectomy, enabling physicians to make more informed treatment decisions concerning adjuvant therapy.     

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.     

High diversity in Delta variant across countries revealed by genome‐wide analysis of SARS‐CoV‐2 beyond the Spike protein

The highly contagious Delta variant of SARS‐CoV‐2 has become a prevalent strain globally and poses a public health challenge around the world. While there has been extensive focus on understanding the amino acid mutations in the Delta variant’s Spike protein, the mutational landscape of the rest of the SARS‐CoV‐2 proteome (25 proteins) remains poorly understood. To this end, we performed a systematic analysis of mutations in all the SARS‐CoV‐2 proteins from nearly 2 million SARS‐CoV‐2 genomes from 176 countries/territories. Six highly prevalent missense mutations in the viral life cycle‐associated Membrane (I82T), Nucleocapsid (R203M, D377Y), NS3 (S26L), and NS7a (V82A, T120I) proteins are almost exclusive to the Delta variant compared to other variants of concern (mean prevalence across genomes: Delta = 99.74%, Alpha = 0.06%, Beta = 0.09%, and Gamma = 0.22%). Furthermore, we find that the Delta variant harbors a more diverse repertoire of mutations across countries compared to the previously dominant Alpha variant. Overall, our study underscores the high diversity of the Delta variant between countries and identifies a list of amino acid mutations in the Delta variant’s proteome for probing the mechanistic basis of pathogenic features such as high viral loads, high transmissibility, and reduced susceptibility against neutralization by vaccines.     

Causes and consequences of variation in early‐life telomere length in a bird metapopulation

Environmental conditions during early‐life development can have lasting effects shaping individual heterogeneity in fitness and fitness‐related traits. The length of telomeres, the DNA sequences protecting chromosome ends, may be affected by early‐life conditions, and telomere length (TL) has been associated with individual performance within some wild animal populations. Thus, knowledge of the mechanisms that generate variation in TL, and the relationship between TL and fitness, is important in understanding the role of telomeres in ecology and life‐history evolution. Here, we investigate how environmental conditions and morphological traits are associated with early‐life blood TL and if TL predicts natal dispersal probability or components of fitness in 2746 wild house sparrow (Passer domesticus) nestlings from two populations sampled across 20 years (1994–2013). We retrieved weather data and we monitored population fluctuations, individual survival, and reproductive output using field observations and genetic pedigrees. We found a negative effect of population density on TL, but only in one of the populations. There was a curvilinear association between TL and the maximum daily North Atlantic Oscillation index during incubation, suggesting that there are optimal weather conditions that result in the longest TL. Dispersers tended to have shorter telomeres than non‐dispersers. TL did not predict survival, but we found a tendency for individuals with short telomeres to have higher annual reproductive success. Our study showed how early‐life TL is shaped by effects of growth, weather conditions, and population density, supporting that environmental stressors negatively affect TL in wild populations. In addition, shorter telomeres may be associated with a faster pace‐of‐life, as individuals with higher dispersal rates and annual reproduction tended to have shorter early‐life TL.     

Aging‐related cell type‐specific pathophysiologic immune responses that exacerbate disease severity in aged COVID‐19 patients

Coronavirus disease 2019 (COVID‐19) is especially severe in aged patients, defined as 65 years or older, for reasons that are currently unknown. To investigate the underlying basis for this vulnerability, we performed multimodal data analyses on immunity, inflammation, and COVID‐19 incidence and severity as a function of age. Our analysis leveraged age‐specific COVID‐19 mortality and laboratory testing from a large COVID‐19 registry, along with epidemiological data of ~3.4 million individuals, large‐scale deep immune cell profiling data, and single‐cell RNA‐sequencing data from aged COVID‐19 patients across diverse populations. We found that decreased lymphocyte count and elevated inflammatory markers (C‐reactive protein, D‐dimer, and neutrophil–lymphocyte ratio) are significantly associated with age‐specific COVID‐19 severities. We identified the reduced abundance of naïve CD8 T cells with decreased expression of antiviral defense genes (i.e., IFITM3 and TRIM22) in aged severe COVID‐19 patients. Older individuals with severe COVID‐19 displayed type I and II interferon deficiencies, which is correlated with SARS‐CoV‐2 viral load. Elevated expression of SARS‐CoV‐2 entry factors and reduced expression of antiviral defense genes (LY6E and IFNAR1) in the secretory cells are associated with critical COVID‐19 in aged individuals. Mechanistically, we identified strong TGF‐beta‐mediated immune–epithelial cell interactions (i.e., secretory‐non‐resident macrophages) in aged individuals with critical COVID‐19. Taken together, our findings point to immuno‐inflammatory factors that could be targeted therapeutically to reduce morbidity and mortality in aged COVID‐19 patients.     

Cancer‐associated mutations in VAV1 trigger variegated signaling outputs and T‐cell lymphomagenesis

Mutations in VAV1, a gene that encodes a multifunctional protein important for lymphocytes, are found at different frequencies in peripheral T‐cell lymphoma (PTCL), non‐small cell lung cancer, and other tumors. However, their pathobiological significance remains unsettled. After cataloguing 51 cancer‐associated VAV1 mutations, we show here that they can be classified in five subtypes according to functional impact on the three main VAV1 signaling branches, GEF‐dependent activation of RAC1, GEF‐independent adaptor‐like, and tumor suppressor functions. These mutations target new and previously established regulatory layers of the protein, leading to quantitative and qualitative changes in VAV1 signaling output. We also demonstrate that the most frequent VAV1 mutant subtype drives PTCL formation in mice. This process requires the concurrent engagement of two downstream signaling branches that promote the chronic activation and transformation of follicular helper T cells. Collectively, these data reveal the genetic constraints associated with the lymphomagenic potential of VAV1 mutant subsets, similarities with other PTCL driver genes, and potential therapeutic vulnerabilities.     

Single‐cell transcriptomics of suprachiasmatic nuclei reveal a Prokineticin‐driven circadian network

Circadian rhythms in mammals are governed by the hypothalamic suprachiasmatic nucleus (SCN), in which 20,000 clock cells are connected together into a powerful time‐keeping network. In the absence of network‐level cellular interactions, the SCN fails as a clock. The topology and specific roles of its distinct cell populations (nodes) that direct network functions are, however, not understood. To characterise its component cells and network structure, we conducted single‐cell sequencing of SCN organotypic slices and identified eleven distinct neuronal sub‐populations across circadian day and night. We defined neuropeptidergic signalling axes between these nodes, and built neuropeptide‐specific network topologies. This revealed their temporal plasticity, being up‐regulated in circadian day. Through intersectional genetics and real‐time imaging, we interrogated the contribution of the Prok2‐ProkR2 neuropeptidergic axis to network‐wide time‐keeping. We showed that Prok2‐ProkR2 signalling acts as a key regulator of SCN period and rhythmicity and contributes to defining the network‐level properties that underpin robust circadian co‐ordination. These results highlight the diverse and distinct contributions of neuropeptide‐modulated communication of temporal information across the SCN.     

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.     

Male‐biased dispersal in a fungus‐gardening ant symbiosis

For nearly all organisms, dispersal is a fundamental life‐history trait that can shape their ecology and evolution. Variation in dispersal capabilities within a species exists and can influence population genetic structure and ecological interactions. In fungus‐gardening (attine) ants, co‐dispersal of ants and mutualistic fungi is crucial to the success of this obligate symbiosis. Female‐biased dispersal (and gene flow) may be favored in attines because virgin queens carry the responsibility of dispersing the fungi, but a paucity of research has made this conclusion difficult. Here, we investigate dispersal of the fungus‐gardening ant Trachymyrmex septentrionalis using a combination of maternally (mitochondrial DNA) and biparentally inherited (microsatellites) markers. We found three distinct, spatially isolated mitochondrial DNA haplotypes; two were found in the Florida panhandle and the other in the Florida peninsula. In contrast, biparental markers illustrated significant gene flow across this region and minimal spatial structure. The differential patterns uncovered from mitochondrial DNA and microsatellite markers suggest that most long‐distance ant dispersal is male‐biased and that females (and concomitantly the fungus) have more limited dispersal capabilities. Consequently, the limited female dispersal is likely an important bottleneck for the fungal symbiont. This bottleneck could slow fungal genetic diversification, which has significant implications for both ant hosts and fungal symbionts regarding population genetics, species distributions, adaptive responses to environmental change, and coevolutionary patterns.     

Opposing activities of IFITM proteins in SARS‐CoV‐2 infection

Interferon‐induced transmembrane proteins (IFITMs) restrict infections by many viruses, but a subset of IFITMs enhance infections by specific coronaviruses through currently unknown mechanisms. We show that SARS‐CoV‐2 Spike‐pseudotyped virus and genuine SARS‐CoV‐2 infections are generally restricted by human and mouse IFITM1, IFITM2, and IFITM3, using gain‐ and loss‐of‐function approaches. Mechanistically, SARS‐CoV‐2 restriction occurred independently of IFITM3 S‐palmitoylation, indicating a restrictive capacity distinct from reported inhibition of other viruses. In contrast, the IFITM3 amphipathic helix and its amphipathic properties were required for virus restriction. Mutation of residues within the IFITM3 endocytosis‐promoting YxxФ motif converted human IFITM3 into an enhancer of SARS‐CoV‐2 infection, and cell‐to‐cell fusion assays confirmed the ability of endocytic mutants to enhance Spike‐mediated fusion with the plasma membrane. Overexpression of TMPRSS2, which increases plasma membrane fusion versus endosome fusion of SARS‐CoV‐2, attenuated IFITM3 restriction and converted amphipathic helix mutants into infection enhancers. In sum, we uncover new pro‐ and anti‐viral mechanisms of IFITM3, with clear distinctions drawn between enhancement of viral infection at the plasma membrane and amphipathicity‐based mechanisms used for endosomal SARS‐CoV‐2 restriction.     

Engineering defensin α‐helix to produce high‐affinity SARS‐CoV‐2 spike protein binding ligands

The binding of severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) spike protein to the angiotensin‐converting enzyme 2 (ACE2) receptor expressed on the host cells is a critical initial step for viral infection. This interaction is blocked through competitive inhibition by soluble ACE2 protein. Therefore, developing high‐affinity and cost‐effective ACE2 mimetic ligands that disrupt this protein–protein interaction is a promising strategy for viral diagnostics and therapy. We employed human and plant defensins, a class of small (2–5 kDa) and highly stable proteins containing solvent‐exposed alpha‐helix, conformationally constrained by two disulfide bonds. Therefore, we engineered the amino acid residues on the constrained alpha‐helix of defensins to mimic the critical residues on the ACE2 helix 1 that interact with the SARS‐CoV‐2 spike protein. The engineered proteins (h‐deface2, p‐deface2, and p‐deface2‐MUT) were soluble and purified to homogeneity with a high yield from a bacterial expression system. The proteins demonstrated exceptional thermostability (Tm 70.7°C), high‐affinity binding to the spike protein with apparent K _d values of 54.4 ± 11.3, 33.5 ± 8.2, and 14.4 ± 3.5 nM for h‐deface2, p‐deface2, and p‐deface2‐MUT, respectively, and were used in a diagnostic assay that detected SARS‐CoV‐2 neutralizing antibodies. This work addresses the challenge of developing helical ACE2 mimetics by demonstrating that defensins provide promising scaffolds to engineer alpha‐helices in a constrained form for designing of high‐affinity ligands.     

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

Pangolins have been suggested as potential reservoir of zoonotic viruses, including SARS‐CoV‐2 causing the global COVID‐19 outbreak. Here, we study the binding of two SARS‐CoV‐2‐like viruses isolated from pangolins, GX/P2V/2017 and GD/1/2019, to human angiotensin‐converting enzyme 2 (hACE2), the receptor of SARS‐CoV‐2. We find that the spike protein receptor‐binding domain (RBD) of pangolin CoVs binds to hACE2 as efficiently as the SARS‐CoV‐2 RBD in vitro. Furthermore, incorporation of pangolin CoV RBDs allows entry of pseudotyped VSV particles into hACE2‐expressing cells. A screen for binding of pangolin CoV RBDs to ACE2 orthologs from various species suggests a broader host range than that of SARS‐CoV‐2. Additionally, cryo‐EM structures of GX/P2V/2017 and GD/1/2019 RBDs in complex with hACE2 show their molecular binding in modes similar to SARS‐CoV‐2 RBD. Introducing the Q498H substitution found in pangolin CoVs into the SARS‐CoV‐2 RBD expands its binding capacity to ACE2 homologs of mouse, rat, and European hedgehog. These findings suggest that these two pangolin CoVs may infect humans, highlighting the necessity of further surveillance of pangolin CoVs.