Akaluc bioluminescence offers superior sensitivity to track in vivo dynamics of SARS-CoV-2 infection

Monitoring in vivo viral dynamics can improve our understanding of pathogenicity and tissue tropism. Because the gene size of RNA viruses is typically small, NanoLuc is the primary choice for accommodation within viral genome. However, NanoLuc/Furimazine and also the conventional firefly luciferase/D-luciferin are known to exhibit relatively low tissue permeability and thus less sensitivity for visualization of deep tissue including lungs. Here, we demonstrated in vivo sufficient visualization of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection using the pair of a codon-optimized Akaluc and AkaLumine. We engineered the codon-optimized Akaluc gene possessing the similar GC ratio of SARS-CoV-2. Using the SARS-CoV-2 recombinants carrying the codon-optimized Akaluc, we visualized in vivo infection of respiratory organs, including the tissue-specific differences associated with particular variants. Additionally, we could evaluate the efficacy of antivirals by monitoring changes in Akaluc signals. Overall, we offer an effective technology for monitoring viral dynamics in live animals.

Longitudinal analysis of genomic mutations in SARS-CoV-2 isolates from persistent COVID-19 patient

A primary reason for the ongoing spread of coronavirus disease 2019 (COVID-19) is the continuous acquisition of mutations by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, the mechanism of acquiring mutations is not fully understood. In this study, we isolated SARS-CoV-2 from an immunocompromized patient persistently infected with Omicron strain BF.5 for approximately 4 months to analyze its genome and evaluate drug resistance. Although the patient was administered the antiviral drug remdesivir (RDV), there were no acquired mutations in RDV binding site, and all isolates exhibited susceptibility to RDV. Notably, upon analyzing the S protein sequence of the day 119 isolate, we identified mutations acquired by mutant strains emerging from the BF.5 variant, suggesting that viral genome analysis in persistent COVID-19 patients may be useful in predicting viral evolution. These results suggest mutations in SARS-CoV-2 are acquired during long-term viral replication rather than in response to antiviral drugs.

Multiple mutations of SARS-CoV-2 Omicron BA.2 variant orchestrate its virological characteristics

IMPORTANCE

Most studies investigating the characteristics of emerging SARS-CoV-2 variants have been focusing on mutations in the spike proteins that affect viral infectivity, fusogenicity, and pathogenicity. However, few studies have addressed how naturally occurring mutations in the non-spike regions of the SARS-CoV-2 genome impact virological properties. In this study, we proved that multiple SARS-CoV-2 Omicron BA.2 mutations, one in the spike protein and another downstream of the spike gene, orchestrally characterize this variant, shedding light on the importance of Omicron BA.2 mutations out of the spike protein.

HTLV-1 bZIP factor impairs DNA mismatch repair system

Adult T-cell leukemia (ATL) is a peripheral T-cell malignancy caused by human T-cell leukemia virus type 1 (HTLV-1). Microsatellite instability (MSI) has been observed in ATL cells. Although MSI results from impaired mismatch repair (MMR) pathway, no null mutations in the genes encoding MMR factors are detectable in ATL cells. Thus, it is unclear whether or not impairment of MMR causes the MSI in ATL cells. HTLV-1 bZIP factor (HBZ) protein interacts with numerous host transcription factors and significantly contributes to disease pathogenesis and progression. Here we investigated the effect of HBZ on MMR in normal cells. The ectopic expression of HBZ in MMR-proficient cells induced MSI, and also suppressed the expression of several MMR factors. We then hypothesized that the HBZ compromises MMR by interfering with a transcription factor, nuclear respiratory factor 1 (NRF-1), and identified the consensus NRF-1 binding site at the promoter of the gene encoding MutS homologue 2 (MSH2), an essential MMR factor. The luciferase reporter assay revealed that NRF-1 overexpression enhanced MSH2 promoter activity, while co-expression of HBZ reversed this enhancement. These results supported the idea that HBZ suppresses the transcription of MSH2 by inhibiting NRF-1. Our data demonstrate that HBZ causes impaired MMR, and may imply a novel oncogenesis driven by HTLV-1.

Pharmacologic characterization of TBP1901, a prodrug form of aglycone curcumin, and CRISPR-Cas9 screen for therapeutic targets of aglycone curcumin

Curcumin (aglycone curcumin) has antitumor properties in a variety of malignancies via the alteration of multiple cancer-related biological pathways; however, its clinical application has been hampered due to its poor bioavailability. To overcome this limitation, we have developed a synthesized curcumin β-D-glucuronide sodium salt (TBP1901), a prodrug form of aglycone curcumin. In this study, we aimed to clarify the pharmacologic characteristics of TBP1901. In β-glucuronidase (GUSB)-proficient mice, both curcumin β-D-glucuronide and its active metabolite, aglycone curcumin, were detected in the blood after TBP1901 injection, whereas only curcumin β-D-glucuronide was detected in GUSB-impaired mice, suggesting that GUSB plays a pivotal role in the conversion of TBP1901 into aglycone curcumin in vivo. TBP1901 itself had minimal antitumor effects in vitro, whereas it demonstrated significant antitumor effects in vivo. Genome-wide clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 screen disclosed the genes associated with NF-κB signaling pathway and mitochondria were among the highest hit. In vitro, aglycone curcumin inhibited NF-kappa B signaling pathways whereas it caused production of reactive oxygen species (ROS). ROS scavenger, N-acetyl-L-cysteine, partially reversed antitumor effects of aglycone curcumin. In summary, TBP1901 can exert antitumor effects as a prodrug of aglycone curcumin through GUSB-dependent activation.

Virological characteristics of the SARS-CoV-2 Omicron BA.2.75 variant

The SARS-CoV-2 Omicron BA.2.75 variant emerged in May 2022. BA.2.75 is a BA.2 descendant but is phylogenetically distinct from BA.5, the currently predominant BA.2 descendant. Here, we show that BA.2.75 has a greater effective reproduction number and different immunogenicity profile than BA.5. We determined the sensitivity of BA.2.75 to vaccinee and convalescent sera as well as a panel of clinically available antiviral drugs and antibodies. Antiviral drugs largely retained potency, but antibody sensitivity varied depending on several key BA.2.75-specific substitutions. The BA.2.75 spike exhibited a profoundly higher affinity for its human receptor, ACE2. Additionally, the fusogenicity, growth efficiency in human alveolar epithelial cells, and intrinsic pathogenicity in hamsters of BA.2.75 were greater than those of BA.2. Our multilevel investigations suggest that BA.2.75 acquired virological properties independent of BA.5, and the potential risk of BA.2.75 to global health is greater than that of BA.5.

Virological characteristics of the SARS-CoV-2 Omicron BA.2 subvariants, including BA.4 and BA.5

After the global spread of the SARS-CoV-2 Omicron BA.2, some BA.2 subvariants, including BA.2.9.1, BA.2.11, BA.2.12.1, BA.4, and BA.5, emerged in multiple countries. Our statistical analysis showed that the effective reproduction numbers of these BA.2 subvariants are greater than that of the original BA.2. Neutralization experiments revealed that the immunity induced by BA.1/2 infections is less effective against BA.4/5. Cell culture experiments showed that BA.2.12.1 and BA.4/5 replicate more efficiently in human alveolar epithelial cells than BA.2, and particularly, BA.4/5 is more fusogenic than BA.2. We further provided the structure of the BA.4/5 spike receptor-binding domain that binds to human ACE2 and considered how the substitutions in the BA.4/5 spike play roles in ACE2 binding and immune evasion. Moreover, experiments using hamsters suggested that BA.4/5 is more pathogenic than BA.2. Our multiscale investigations suggest that the risk of BA.2 subvariants, particularly BA.4/5, to global health is greater than that of original BA.2.

A panel of nanobodies recognizing conserved hidden clefts of all SARS-CoV-2 spike variants including Omicron

We are amid the historic coronavirus infectious disease 2019 (COVID-19) pandemic. Imbalances in the accessibility of vaccines, medicines, and diagnostics among countries, regions, and populations, and those in war crises, have been problematic. Nanobodies are small, stable, customizable, and inexpensive to produce. Herein, we present a panel of nanobodies that can detect the spike proteins of five SARS-CoV-2 variants of concern (VOCs) including Omicron. Here we show via ELISA, lateral flow, kinetic, flow cytometric, microscopy, and Western blotting assays that our nanobodies can quantify the spike variants. This panel of nanobodies broadly neutralizes viral infection caused by pseudotyped and authentic SARS-CoV-2 VOCs. Structural analyses show that the P86 clone targets epitopes that are conserved yet unclassified on the receptor-binding domain (RBD) and contacts the N-terminal domain (NTD). Human antibodies rarely access both regions; consequently, the clone buries hidden crevasses of SARS-CoV-2 spike proteins that go undetected by conventional antibodies.

A panel of nanobodies are presented that can detect the spike proteins of five SARS-CoV-2 variants of concern and structural analyses show that one clone targets conserved epitopes on the receptor-binding domain and contacts the N-terminal domain.

Virological characteristics of the SARS-CoV-2 Omicron BA.2 spike

Soon after the emergence and global spread of the SARS-CoV-2 Omicron lineage BA.1, another Omicron lineage, BA.2, began outcompeting BA.1. The results of statistical analysis showed that the effective reproduction number of BA.2 is 1.4-fold higher than that of BA.1. Neutralization experiments revealed that immunity induced by COVID vaccines widely administered to human populations is not effective against BA.2, similar to BA.1, and that the antigenicity of BA.2 is notably different from that of BA.1. Cell culture experiments showed that the BA.2 spike confers higher replication efficacy in human nasal epithelial cells and is more efficient in mediating syncytia formation than the BA.1 spike. Furthermore, infection experiments using hamsters indicated that the BA.2 spike-bearing virus is more pathogenic than the BA.1 spike-bearing virus. Altogether, the results of our multiscale investigations suggest that the risk of BA.2 to global health is potentially higher than that of BA.1.

Virological characteristics of the SARS-CoV-2 Omicron BA.2 spike

Soon after the emergence and global spread of the SARS-CoV-2 Omicron lineage BA.1, another Omicron lineage, BA.2, began outcompeting BA.1. The results of statistical analysis showed that the effective reproduction number of BA.2 is 1.4-fold higher than that of BA.1. Neutralization experiments revealed that immunity induced by COVID vaccines widely administered to human populations is not effective against BA.2, similar to BA.1, and that the antigenicity of BA.2 is notably different from that of BA.1. Cell culture experiments showed that the BA.2 spike confers higher replication efficacy in human nasal epithelial cells and is more efficient in mediating syncytia formation than the BA.1 spike. Furthermore, infection experiments using hamsters indicated that the BA.2 spike-bearing virus is more pathogenic than the BA.1 spike-bearing virus. Altogether, the results of our multiscale investigations suggest that the risk of BA.2 to global health is potentially higher than that of BA.1.

Altered TMPRSS2 usage by SARS-CoV-2 Omicron impacts infectivity and fusogenicity

The SARS-CoV-2 Omicron BA.1 variant emerged in 20211 and has multiple mutations in its spike protein2. Here we show that the spike protein of Omicron has a higher affinity for ACE2 compared with Delta, and a marked change in its antigenicity increases Omicron's evasion of therapeutic monoclonal and vaccine-elicited polyclonal neutralizing antibodies after two doses. mRNA vaccination as a third vaccine dose rescues and broadens neutralization. Importantly, the antiviral drugs remdesivir and molnupiravir retain efficacy against Omicron BA.1. Replication was similar for Omicron and Delta virus isolates in human nasal epithelial cultures. However, in lung cells and gut cells, Omicron demonstrated lower replication. Omicron spike protein was less efficiently cleaved compared with Delta. The differences in replication were mapped to the entry efficiency of the virus on the basis of spike-pseudotyped virus assays. The defect in entry of Omicron pseudotyped virus to specific cell types effectively correlated with higher cellular RNA expression of TMPRSS2, and deletion of TMPRSS2 affected Delta entry to a greater extent than Omicron. Furthermore, drug inhibitors targeting specific entry pathways3 demonstrated that the Omicron spike inefficiently uses the cellular protease TMPRSS2, which promotes cell entry through plasma membrane fusion, with greater dependency on cell entry through the endocytic pathway. Consistent with suboptimal S1/S2 cleavage and inability to use TMPRSS2, syncytium formation by the Omicron spike was substantially impaired compared with the Delta spike. The less efficient spike cleavage of Omicron at S1/S2 is associated with a shift in cellular tropism away from TMPRSS2-expressing cells, with implications for altered pathogenesis.

Characterization of the immune resistance of SARS-CoV-2 Mu variant and the robust immunity induced by Mu infection

We have recently revealed that the SARS-CoV-2 Mu variant shows a pronounced resistance to antibodies elicited by natural SARS-CoV-2 infection and vaccination. However, it remains unclear which mutations determine the resistance of SARS-CoV-2 Mu to antiviral sera. Also, it is unclear how SARS-CoV-2 Mu infection induces antiviral immunity. Here we reveal that the two mutations in the SARS-CoV-2 Mu spike protein, YY144-145TSN and E484K, are responsible for the resistance to COVID-19 convalescent sera during early 2020 and vaccine sera. Notably, the convalescent sera of SARS-CoV-2 Mu-infected individuals are broadly antiviral against Mu as well as other SARS-CoV-2 variants of concern/interest.

The SARS-CoV-2 Lambda variant exhibits enhanced infectivity and immune resistance

SARS-CoV-2 Lambda, a variant of interest, has spread in some South American countries; however, its virological features and evolutionary traits remain unclear. In this study, we use pseudoviruses and reveal that the spike protein of the Lambda variant is more infectious than that of other variants due to the T76I and L452Q mutations. The RSYLTPGD246-253N mutation, a unique 7-amino acid deletion in the N-terminal domain of the Lambda spike protein, is responsible for evasion from neutralizing antibodies and further augments antibody-mediated enhancement of infection. Although this mutation generates a nascent N-linked glycosylation site, the additional N-linked glycan is dispensable for the virological property conferred by this mutation. Since the Lambda variant has dominantly spread according to the increasing frequency of the isolates harboring the RSYLTPGD246-253N mutation, our data suggest that the RSYLTPGD246-253N mutation is closely associated with the substantial spread of the Lambda variant in South America.

Enhanced fusogenicity and pathogenicity of SARS-CoV-2 Delta P681R mutation

During the current coronavirus disease 2019 (COVID-19) pandemic, a variety of mutations have accumulated in the viral genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and, at the time of writing, four variants of concern are considered to be potentially hazardous to human society¹. The recently emerged B.1.617.2/Delta variant of concern is closely associated with the COVID-19 surge that occurred in India in the spring of 2021 (ref. ²). However, the virological properties of B.1.617.2/Delta remain unclear. Here we show that the B.1.617.2/Delta variant is highly fusogenic and notably more pathogenic than prototypic SARS-CoV-2 in infected hamsters. The P681R mutation in the spike protein, which is highly conserved in this lineage, facilitates cleavage of the spike protein and enhances viral fusogenicity. Moreover, we demonstrate that the P681R-bearing virus exhibits higher pathogenicity compared with its parental virus. Our data suggest that the P681R mutation is a hallmark of the virological phenotype of the B.1.617.2/Delta variant and is associated with enhanced pathogenicity.

Neutralization of the SARS-CoV-2 Mu Variant by Convalescent and Vaccine Serum

On August 30, 2021, the WHO classified the SARS-CoV-2 Mu variant (B.1.621 lineage) as a new variant of interest. The WHO defines “comparative assessment of virus characteristics and public health risks” as primary action in response to the emergence of new SARS-CoV-2 variants. Here, we demonstrate that the Mu variant is highly resistant to sera from COVID-19 convalescents and BNT162b2-vaccinated individuals. Direct comparison of different SARS-CoV-2 spike proteins revealed that Mu spike is more resistant to serum-mediated neutralization than all other currently recognized variants of interest (VOI) and concern (VOC). This includes the Beta variant (B.1.351) that has been suggested to represent the most resistant variant to convalescent and vaccinated sera to date (e.g., Collier et al, Nature, 2021; Wang et al, Nature, 2021). Since breakthrough infection by newly emerging variants is a major concern during the current COVID-19 pandemic (Bergwerk et al., NEJM, 2021), we believe that our findings are of significant public health interest. Our results will help to better assess the risk posed by the Mu variant for vaccinated, previously infected and naïve populations.

SARS-CoV-2 B.1.617 Mutations L452R and E484Q Are Not Synergistic for Antibody Evasion

The SARS-CoV-2 B.1.617 variant emerged in the Indian state of Maharashtra in late 2020. There have been fears that 2 key mutations seen in the receptor-binding domain, L452R and E484Q, would have additive effects on evasion of neutralizing antibodies. We report that spike bearing L452R and E484Q confers modestly reduced sensitivity to BNT162b2 mRNA vaccine-elicited antibodies following either first or second dose. The effect is similar in magnitude to the loss of sensitivity conferred by L452R or E484Q alone. These data demonstrate reduced sensitivity to vaccine-elicited neutralizing antibodies by L452R and E484Q but lack of synergistic loss of sensitivity.

FRET-Based Detection and Quantification of HIV-1 Virion Maturation

HIV-1 infectivity is achieved through virion maturation. Virus particles undergo structural changes via cleavage of the Gag polyprotein mediated by the viral protease, causing the transition from an uninfectious to an infectious status. The majority of proviruses in people living with HIV-1 treated with combination antiretroviral therapy are defective with large internal deletions. Defective proviral DNA frequently preserves intact sequences capable of expressing viral structural proteins to form virus-like particles whose maturation status is an important factor for chronic antigen-mediated immune stimulation and inflammation. Thus, novel methods to study the maturation capability of defective virus particles are needed to characterize their immunogenicity. To build a quantitative tool to study virion maturation in vitro, we developed a novel single virion visualization technique based on fluorescence resonance energy transfer (FRET). We inserted an optimized intramolecular CFP-YPF FRET donor-acceptor pair bridged with an HIV-1 protease cleavage sequence between the Gag MA-CA domains. This system allowed us to microscopically distinguish mature and immature virions via their FRET signal when the FRET donor and acceptor proteins were separated by the viral protease during maturation. We found that approximately 80% of the FRET labeled virus particles were mature with equivalent infectivity to wild type. The proportion of immature virions was increased by treatment of virus producer cells with a protease inhibitor in a dose-dependent manner, which corresponded to a relative decrease in infectivity. Potential areas of application for this tool are assessing maturation efficiency in different cell type settings of intact or deficient proviral DNA integrated cells. We believe that this FRET-based single-virion imaging platform will facilitate estimating the impact on the immune system of both extracellular intact and defective viruses by quantifying the Gag maturation status.

APOBEC3B is preferentially expressed at the G2/M phase of cell cycle

APOBEC3B (A3B) is a cytosine deaminase that converts cytosine to uracil in single-stranded DNA. Cytosine-to-thymine and cytosine-to-guanine base substitution mutations in trinucleotide motifs (APOBEC mutational signatures) were found in various cancers including lymphoid hematological malignancies such as multiple myeloma and A3B has been shown to be an enzymatic source of mutations in those cancers. Although the importance of A3B is being increasingly recognized, it is unclear how A3B expression is regulated in cancer cells as well as normal cells. To answer these fundamental questions, we analyzed 1276 primary myeloma cells using single-cell RNA-sequencing (scRNA-seq) and found that A3B was preferentially expressed at the G2/M phase, in sharp contrast to the expression patterns of other APOBEC3 genes. Consistently, we demonstrated that A3B protein was preferentially expressed at the G2/M phase in myeloma cells by cell sorting. We also demonstrated that normal blood cells expressing A3B were also enriched in G2/M-phase cells by analyzing scRNA-seq data from 86,493 normal bone marrow mononuclear cells. Furthermore, we revealed that A3B was expressed mainly in plasma cells, CD10⁺ B cells and erythroid cells, but not in granulocyte-macrophage progenitors. A3B expression profiling in normal blood cells may contribute to understanding the defense mechanism of A3B against viruses, and partially explain the bias of APOBEC mutational signatures in lymphoid but not myeloid malignancies. This study identified the cells and cellular phase in which A3B is highly expressed, which may help reveal the mechanisms behind carcinogenesis and cancer heterogeneity, as well as the biological functions of A3B in normal blood cells.

APOBEC3B reporter myeloma cell lines identify DNA damage response pathways leading to APOBEC3B expression

Apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like (APOBEC) DNA cytosine deaminase 3B (A3B) is a DNA editing enzyme which induces genomic DNA mutations in multiple myeloma and in various other cancers. APOBEC family proteins are highly homologous so it is especially difficult to investigate the biology of specifically A3B in cancer cells. To easily and comprehensively investigate A3B function in myeloma cells, we used CRISPR/Cas9 to generate A3B reporter cells that contain 3×FLAG tag and IRES-EGFP sequences integrated at the end of the A3B gene. These reporter cells stably express 3xFLAG tagged A3B and the reporter EGFP and this expression is enhanced by known stimuli, such as PMA. Conversely, shRNA knockdown of A3B decreased EGFP fluorescence and 3xFLAG tagged A3B protein levels. We screened a series of anticancer treatments using these cell lines and identified that most conventional therapies, such as antimetabolites or radiation, exacerbated endogenous A3B expression, but recent molecular targeted therapeutics, including bortezomib, lenalidomide and elotuzumab, did not. Furthermore, chemical inhibition of ATM, ATR and DNA-PK suppressed EGFP expression upon treatment with antimetabolites. These results suggest that DNA damage triggers A3B expression through ATM, ATR and DNA-PK signaling.

Coherence between oscillations in the cardiorespiratory system and tissue oxygen index in muscle recovering from intensive exercise in humans

It has been shown that the tissue oxygen index (TOI) measured by near-infrared spectroscopy oscillates at very low frequencies during recovery after exercise and that this oscillation is derived from interactions among biochemical substances involved in oxidative metabolism in skeletal muscle. As a further step, we examined whether TOI in muscle interacts through oscillation with factors related to oxygen in the cardiorespiratory system. For this examination, coherence and phase difference between the TOI in the vastus lateralis and heart rate (HR) and between TOI and arterial oxygen saturation (SpO₂) were sequentially determined during recovery (2-60 min) after severe cycle exercise with a workload of 7.5% of body weight for 20 s. Significant coherence between TOI and HR was obtained in the very low-frequency band (approximate range: 0.002-0.03 Hz) and in the low-frequency band (approximate range: 0.06-0.12 Hz). The phase difference was negative in the low-frequency band and positive in the very low-frequency band. The coherence between TOI and SpO₂ was significant in the very low-frequency band. The phase difference was negative. There were no sequential changes in these coherences and phase differences. The results suggest that TOI in skeletal muscle interrelates with factors related to the heart and lungs.

Structural Determinants of the APOBEC3G N-Terminal Domain for HIV-1 RNA Association

APOBEC3G (A3G) is a cellular protein that inhibits HIV-1 infection through virion incorporation. The interaction of the A3G N-terminal domain (NTD) with RNA is essential for A3G incorporation in the HIV-1 virion. The interaction between A3G-NTD and RNA is not completely understood. The A3G-NTD is also recognized by HIV-1 Viral infectivity factor (Vif) and A3G-Vif binding leads to A3G degradation. Therefore, the A3G-Vif interaction is a target for the development of antiviral therapies that block HIV-1 replication. However, targeting the A3G-Vif interactions could disrupt the A3G-RNA interactions that are required for A3G's antiviral activity. To better understand A3G-RNA binding, we generated in silico docking models to simulate the RNA-binding propensity of A3G-NTD. We simulated the A3G-NTD residues with high RNA-binding propensity, experimentally validated our prediction by testing A3G-NTD mutations, and identified structural determinants of A3G-RNA binding. In addition, we found a novel amino acid residue, I26 responsible for RNA interaction. The new structural insights provided here will facilitate the design of pharmaceuticals that inhibit A3G-Vif interactions without negatively impacting A3G-RNA interactions.