A giant virus genome is densely packaged by stable nucleosomes within virions.

The two doublet histones of Marseillevirus are distantly related to the four eukaryotic core histones and wrap 121 base pairs of DNA to form remarkably similar nucleosomes. By permeabilizing Marseillevirus virions and performing genome-wide nuclease digestion, chemical cleavage, and mass spectrometry assays, we find that the higher-order organization of Marseillevirus chromatin fundamentally differs from that of eukaryotes. Marseillevirus nucleosomes fully protect DNA within virions as closely abutted 121-bp DNA-wrapped cores without linker DNA or phasing along genes. Likewise, we observed that nucleosomes reconstituted onto multi-copy tandem repeats of a nucleosome-positioning sequence are tightly packed. Dense promiscuous packing of fully wrapped nucleosomes rather than "beads on a string" with genic punctuation represents a distinct mode of DNA packaging by histones. We suggest that doublet histones have evolved for viral genome protection and may resemble an early stage of histone differentiation leading to the eukaryotic octameric nucleosome.

Vγ9Vδ2 T-cells Are Potent Inhibitors of SARS-CoV-2 Replication and Represent Effector Phenotypes in Patients With COVID-19.

Vγ9Vδ2 T cells play a key role in the innate immune response to viral infections through butyrophilin 3A (BTN3A). Here, we report blood Vγ9Vδ2 T cells decreased in clinically mild COVID-19 compared to healthy volunteers, and this was maintained up to 28 days and in the recovery period. Terminally differentiated Vγ9Vδ2 T cells tended to be enriched on the day of diagnosis, 28 days after, and during the recovery period. These cells showed cytotoxic and inflammatory activities following anti-BTN3A activation. BTN3A upregulation and Vγ9Vδ2 T-cell infiltration were observed in a lung biopsy from a fatal SARS-CoV-2 infection. In vitro, SARS-CoV-2 infection increased BTN3A expression in macrophages and lung cells that enhanced the anti-SARS-CoV-2 Vγ9Vδ2 T-cell cytotoxicity and interferon-γ and tumor necrosis factor-α. Increasing concentrations of anti-BTN3A lead to viral replication inhibition. Altogether, we report Vγ9Vδ2 T cells are important in the immune response against SARS-CoV-2 infection and activation by anti-BTN3A antibody may enhance their response. Clinical Trials Registration. NCT04816760.

Large-Scale Outbreak of Mycoplasma pneumoniae Infection, Marseille, France, 2023-2024.

We report a large-scale outbreak of Mycoplasma pneumoniae respiratory infections encompassing 218 cases (0.8% of 26,449 patients tested) during 2023-2024 in Marseille, France. The bacterium is currently circulating and primarily affects children <15 years of age. High prevalence of co-infections warrants the use of a syndromic diagnostic strategy.

Resistance-associated mutations to the anti-SARS-CoV-2 agent nirmatrelvir: Selection not induction.

Mutations associated with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) resistance to antiprotease nirmatrelvir were reported. We aimed to detect them in SARS-CoV-2 genomes and quasispecies retrieved in our institute before drug availability in January 2022 and to analyze the impact of mutations on protease (3CLpro) structure. We sought for 38 3CLpro nirmatrelvir resistance mutations in a set of 62 673 SARS-CoV-2 genomes obtained in our institute from respiratory samples collected between 2020 and 2023 and for these mutations in SARS-CoV-2 quasispecies for 90 samples collected in 2020, using Python. SARS-CoV-2 protease with major mutation E166V was generated with Swiss Pdb Viewer and Molegro Molecular Viewer. We detected 22 (58%) of the resistance-associated mutations in 417 (0.67%) of the genomes analyzed; 325 (78%) of these genomes had been obtained from samples collected in 2020-2021. APOBEC signatures were found for 12/22 mutations. We also detected among viral quasispecies from 90 samples some minority reads harboring any of 15 nirmatrelvir resistance mutations, including E166V. Also, we predicted that E166V has a very limited effect on 3CLpro structure but may prevent drug attachment. Thus, we evidenced that mutations associated with nirmatrelvir resistance pre-existed in SARS-CoV-2 before drug availability. These findings further warrant SARS-CoV-2 genomic surveillance and SARS-CoV-2 quasispecies characterization.

Pan-microscopic examination of monkeypox virus in trophoblasts cells reveals new insights into virions release through filopodia-like projections.

Vertical transmission has been described following monkeypox virus (MPXV) infection in pregnant women. The presence of MPXV has been reported in the placenta from infected women, but whether pathogens colonize placenta remains unexplored. We identify trophoblasts as a target cell for MPXV replication. In a pan-microscopy approach, we decipher the specific infectious cycle of MPXV and inner cellular structures in trophoblasts. We identified the formation of a specialized region for viral morphogenesis and replication in placental cells. We also reported infection-induced cellular remodeling. We found that MPXV stimulates cytoskeleton reorganization with intercellular extensions for MPXV cell spreading specifically to trophoblastic cells. Altogether, the specific infectious cycle of MPXV in trophoblast cells and these protrusions that were structurally and morphologically similar to filopodia reveal new insights into the infection of MPXV.

Sequencing and characterization of human bocavirus genomes from patients diagnosed in Southern France between 2017 and 2022.

The diversity and evolution of the genomes of human bocavirus (HBoV), which causes respiratory diseases, have been scarcely studied. Here, we aimed to obtain and characterize HBoV genomes from patients's nasopharyngeal samples collected between 2017 and 2022 period (5 years and 7 months). Next-generation sequencing (NGS) used Illumina technology after having implemented using GEMI an in-house multiplex PCR amplification strategy. Genomes were assembled and analyzed with CLC Genomics, Mafft, BioEdit, MeV, Nextclade, MEGA, and iTol. A total of 213 genomes were obtained. Phylogeny classified them all as of Bocavirus 1 (HBoV1) species. Five HBoV1 genotypic clusters determined by hierarchical clustering analysis of 27 variable genome positions were scattered over the study period although with differences in yearly prevalence. A total of 167 amino acid substitutions were detected. Besides, coinfection was observed for 52% of the samples, rhinoviruses then adenoviruses (HAdVs) being the most common viruses. Principal component analysis showed that HBoV1 genotypic cluster α tended to be correlated with HAdV co-infection. Subsequent HAdV typing for HBoV1-positive samples and negative controls demonstrated that HAdVC species predominated but HAdVB was that significantly HBoV1-associated. Overall, we described here the first HBoV1 genomes sequenced for France. HBoV1 and HAdVB association deserves further investigation.

The return of the "Mistigri" (virus adaptative gain by gene loss) through the SARS-CoV-2 XBB.1.5 chimera that predominated in 2023.

Severe acute respiratory syndrome coronavirus 2 XBB.1.5 is the first recombinant lineage to predominate at the country and global scales. Very interestingly, like the Marseille-4B subvariant (or B.1.160) and the pandemic variant B.1.1.7 (or Alpha) previously, it has its ORF8 gene inactivated by a stop codon. We aimed here to study the distribution of stop codons in ORF8 of XBB.1.5 and non-XBB.1.5 genomes. We identified that a stop codon was present at 89 (74%) ORF8 codons in ≥1 of 15 222 404 genomes available in GISAID. The mean proportion of genomes with a stop codon per codon was 0.11% (range, 0%-7.8%). In addition, a stop codon was detected at 15 (12%) codons in at least 1000 genomes. These 15 codons are notably located on seven stem-loop hairpin regions and in the signal peptide region for the case of the XBB.1.5 lineage (codon 8). Thus, it is very likely that stop codons in ORF8 gene contributed on at least three occasions and independently during the pandemic to the evolutionary success of a lineage that became transiently predominant. Such association of gene loss with evolutionary success, which suits the recently described Mistigri rule, is an important biological phenomenon very unknown in virology while largely described in cellular organisms.

Emergence of a second SARS-CoV-2 variant with a tremendous genetic leap from its ancestors.

The on-going emergence of the Omicron BA.2.86 variant is one of the major events in SARS-CoV-2 genetic evolution that remain enigmatic regarding the overall virus mutation rate, together with the emergence of the initial Omicron variant, BA.1. Indeed, the genomes of the Omicron BA.2.86 lineage, an offspring of the second major Omicron subvariant, BA.2, harbor 39 additional mutations in the spike compared to this ancestor. Here we comment on the phylogeny of BA.2.86, on the positions, and frequencies in other SARS-CoV-2, of mutations in its spike, and on the structural model of this protein that concentrates most of BA.2.86 additional mutations.

Sequencing of monkeypox virus from infected patients reveals viral genomes with APOBEC3-like editing, gene inactivation, and bacterial agents of skin superinfection.

A large outbreak of Monkeypox virus (MPXV) infections has arisen in May 2022 in nonendemic countries. Here, we performed DNA metagenomics using next-generation sequencing with Illumina or Nanopore technologies for clinical samples from MPXV-infected patients diagnosed between June and July 2022. Classification of the MPXV genomes and determination of their mutational patterns were performed using Nextclade. Twenty-five samples from 25 patients were studied. A MPXV genome was obtained for 18 patients, essentially from skin lesions and rectal swabbing. All 18 genomes were classified in clade IIb, lineage B.1, and we identified four B.1 sublineages (B.1.1, B.1.10, B.1.12, B.1.14). We detected a high number of mutations (range, 64-73) relatively to a 2018 Nigerian genome (genome GenBank Accession no. NC_063383.1), which were harbored by a large part of a set of 3184 MPXV genomes of lineage B.1 recovered from GenBank and Nextstrain; and we detected 35 mutations relatively to genome ON563414.3 (a B.1 lineage reference genome). Nonsynonymous mutations occurred in genes encoding central proteins, among which transcription factors and core and envelope proteins, and included two mutations that would truncate a RNA polymerase subunit and a phospholipase d-like protein, suggesting an alternative start codon and gene inactivation, respectively. A large majority (94%) of nucleotide substitutions were G > A or C > U, suggesting the action of human APOBEC3 enzymes. Finally, >1000 reads were identified as from Staphylococcus aureus and Streptococcus pyogenes for 3 and 6 samples, respectively. These findings warrant a close genomic monitoring of MPXV to get a better picture of the genetic micro-evolution and mutational patterns of this virus, and a close clinical monitoring of skin bacterial superinfection in monkeypox patients.

How to distinguish mpox from its mimickers: An observational retrospective cohort study.

During the current global outbreak of mpox (formerly monkeypox), atypical features were frequently described outside endemic areas, raising concerns around differential diagnosis. In this study, we included 372 adult patients who had clinical signs consistent with mpox and who were screened using non-variola orthopoxvirus specific quantitative polymerase chain reaction (PCR) between 15 May and 15 November 2022 at the University Hospital Institute Méditerranée Infection, Marseille, France. At least one clinical sample was positive for 143 (38.4%) of these patients and 229 (61.6%) were negative. Clinically, patients who had mpox presented more frequently with systemic signs (69.9% vs. 31.0%, p < 10-6 ) including fever (51.0% vs. 30.1%, p < 10-3 ), myalgia (33.5% vs. 17.9%, p = 0.002), and lymphadenopathy (38.5% vs. 13.1%, p < 10-6 ). Among the patients who were negative for the non-variola orthopoxvirus, an alternative diagnosis was identified in 58 of them (25.3%), including chickenpox (n = 30, 13.1%), syphilis (n = 9, 4%), bacterial skin infection (n = 8, 3.5%), gonococcus (n = 5, 2.2%), HSV infection (n = 5, 2.2%), and histoplasmosis (n = 1, 0.4%). Overall, in the current outbreak, we show that mpox has a poorly specific clinical presentation. This reinforces the importance of microbiological confirmation. In symptomatic patients who are negative for the monkeypox virus by PCR, a broad differential diagnosis should be maintained.

The emergence, spread and vanishing of a French SARS-CoV-2 variant exemplifies the fate of RNA virus epidemics and obeys the Mistigri rule.

The nature and dynamics of mutations associated with the emergence, spread, and vanishing of SARS-CoV-2 variants causing successive waves are complex. We determined the kinetics of the most common French variant ("Marseille-4") for 10 months since its onset in July 2020. Here, we analyzed and classified into subvariants and lineages 7453 genomes obtained by next-generation sequencing. We identified two subvariants, Marseille-4A, which contains 22 different lineages of at least 50 genomes, and Marseille-4B. Their average lifetime was 4.1 ± 1.4 months, during which 4.1 ± 2.6 mutations accumulated. Growth rate was 0.079 ± 0.045, varying from 0.010 to 0.173. Most of the lineages exhibited a bell-shaped distribution. Several beneficial mutations at unpredicted sites initiated a new outbreak, while the accumulation of other mutations resulted in more viral heterogenicity, increased diversity and vanishing of the lineages. Marseille-4B emerged when the other Marseille-4 lineages vanished. Its ORF8 gene was knocked out by a stop codon, as reported in SARS-CoV-2 of mink and in the Alpha variant. This subvariant was associated with increased hospitalization and death rates, suggesting that ORF8 is a nonvirulence gene. We speculate that the observed heterogenicity of a lineage may predict the end of the outbreak.

New putative periodontopathogens and periodontal health-associated species: A systematic review and meta-analysis.

To investigate the existence of any association between new putative periodontal pathogens and periodontitis. Two independent reviewers conducted electronic literature searches in the MEDLINE (PubMed), EMBASE, DOSS and Google Scholar databases as well as a manual search to identify eligible clinical studies prior to November 2022. Studies comparing the prevalence of microorganisms other than the already-known periodontal pathogens in subgingival plaque and/or saliva samples between subjects with periodontitis and subject with periodontal health were included. Meta-analyses were performed on data provided by the included studies. Fifty studies including a total of 2739 periodontitis subjects and 1747 subjects with periodontal health were included. The Archaea domain and 25 bacterial species (Anaeroglobus geminatus, Bacteroidales [G-2] bacterium HMT 274, Desulfobulbus sp. HMT 041, Dialister invisus, Dialister pneumosintes, Eubacterium brachy, Enterococcus faecalis, Eubacterium nodatum, Eubacterium saphenum, Filifactor alocis, Fretibacterium sp. HMT 360, Fretibacterium sp. HMT 362, Mogibacterium timidum, Peptoniphilaceae sp. HMT 113, Peptostreptococcus stomatis, Porphyromonas endodontalis, Slackia exigua, Streptococcus gordonii, Selenomonas sputigena, Treponema amylovorum, Treponema lecithinolyticum, Treponema maltophilum, Treponema medium, Treponema parvum and Treponema socranskii) were found to be statistically significantly associated with periodontitis. Network studies should be conducted to investigate the role of these newly identified periodontitis-associated microorganisms through interspecies interaction and host-microbe crosstalk analyses.

Yaravirus: A novel 80-nm virus infecting Acanthamoeba castellanii.

Here we report the discovery of Yaravirus, a lineage of amoebal virus with a puzzling origin and evolution. Yaravirus presents 80-nm-sized particles and a 44,924-bp dsDNA genome encoding for 74 predicted proteins. Yaravirus genome annotation showed that none of its genes matched with sequences of known organisms at the nucleotide level; at the amino acid level, six predicted proteins had distant matches in the nr database. Complimentary prediction of three-dimensional structures indicated possible function of 17 proteins in total. Furthermore, we were not able to retrieve viral genomes closely related to Yaravirus in 8,535 publicly available metagenomes spanning diverse habitats around the globe. The Yaravirus genome also contained six types of tRNAs that did not match commonly used codons. Proteomics revealed that Yaravirus particles contain 26 viral proteins, one of which potentially representing a divergent major capsid protein (MCP) with a predicted double jelly-roll domain. Structure-guided phylogeny of MCP suggests that Yaravirus groups together with the MCPs of Pleurochrysis endemic viruses. Yaravirus expands our knowledge of the diversity of DNA viruses. The phylogenetic distance between Yaravirus and all other viruses highlights our still preliminary assessment of the genomic diversity of eukaryotic viruses, reinforcing the need for the isolation of new viruses of protists.

Correlation of Redondovirus and Entamoeba gingivalis Detections in the Human Oral Cavity Suggests That This Amoeba Is Possibly the Redondovirus Host.

The virome of the human oral cavity and the relationships between viruses and diseases such as periodontitis are scarcely deciphered. Redondoviruses were reported in the human oral cavity in 2019, including in periodontitis patients. Here, we aimed at detecting redondoviruses and at searching for a potential viral host in human saliva. Non-stimulated saliva was collected between December 2020 and June 2021. These samples were tested using real-time PCR regarding the presence of redondovirus and Entamoeba gingivalis DNA. Similarity searches were performed using BLAST against eukaryotic and prokaryotic sequences from GenBank. The redondovirus DNA was detected in 46% of the 28 human saliva samples. In addition, short fragments of redondovirus genomes were detected in silico within Entamoeba sequences. Finally, Entamoeba gingivalis DNA was detected in 46% of the 28 saliva samples, with a strong correlation between redondovirus DNA and E. gingivalis DNA detections, in 93% of the cases. Regarded together, these findings and previous ones strongly support the presence of redondoviruses in the human oral cavity and their association to E. gingivalis as their likely host.

Observational Cohort Study of Evolving Epidemiologic, Clinical, and Virologic Features of Monkeypox in Southern France.

We enrolled 136 patients with laboratory-confirmed monkeypox during June 4-August 31, 2022, at the University Hospital Institute Méditerranée Infection in Marseille, France. The median patient age was 36 years (interquartile range 31-42 years). Of 136 patients, 125 (92%) were men who have sex with men, 15 (11%) reported previous smallpox vaccinations, and 21 (15.5%) were HIV-positive. The most frequent lesion locations were the genitals (68 patients, 53%), perianal region (65 patients, 49%), and oral/perioral area (22 patients, 17%). Lesion locations largely corresponded with the route of contamination. Most (68%) patients had isolated anal, genital, or oral lesions when they were first seen, including 56 (61%) who had >1 positive site without a visible lesion. Concurrent sexually transmitted infections were diagnosed in 19 (15%) patients, and 7 patients (5%) were asymptomatic. We recommend vaccination campaigns, intensified testing for sexually transmitted infections, and increased contact tracing to control the ongoing monkeypox outbreak.

The puzzling mutational landscape of the SARS-2-variant Omicron.

The recently emerging SARS-CoV-2 variant omicron displays an unusual association of 30 mutations, 3 deletions, and 1 insertion. To analyze the impact of this atypic mutational landscape, we constructed a complete structure of the omicron spike protein. Compared with the delta variant, the receptor-binding domain (RBD) of omicron has an increased electrostatic surface potential, but a decreased affinity for the ACE-2 receptor. The N-terminal domain (NTD) has both a decreased surface potential and a lower affinity for lipid rafts. The omicron variant is predicted to be less fusogenic and thus less pathogenic than delta, due to a geometric reorganization of the S1-S2 cleavage site. Overall, these virological parameters suggest that omicron does not have a significant infectivity advantage over the delta variant. However, in omicron, neutralizing epitopes are greatly affected, suggesting that current vaccines will probably confer little protection against this variant. In conclusion, the puzzling mutational pattern of the omicron variant combines contradictory properties which may either decrease (virological properties) or increase (immunological escape/facilitation) the transmission of this variant in the human population. This Janus-like phenotype may explain some conflicting reports on the initial assessment of omicron and provide new insights about the molecular mechanisms controlling its dissemination and pathogenesis worldwide.

Culture and identification of a "Deltamicron" SARS-CoV-2 in a three cases cluster in southern France.

Multiple SARS-CoV-2 variants have successively, or concomitantly spread worldwide since the summer of 2020. A few co-infections with different variants were reported and genetic recombinations, common among coronaviruses, were reported or suspected based on co-detection of signature mutations of different variants in a given genome. Here we report three infections in southern France with a Delta 21J_AY.4-Omicron 21K/BA.1 "Deltamicron" recombinant. The hybrid genome harbors signature mutations of the two lineages, supported by a mean sequencing depth of 1163-1421 reads and a mean nucleotide diversity of 0.1%-0.6%. It is composed of the near full-length spike gene (from codons 156-179) of an Omicron 21K/BA.1 variant in a Delta 21J/AY.4 lineage backbone. Importantly, we cultured an isolate of this recombinant and sequenced its genome. It was observed by scanning electron microscopy. As it is misidentified with current variant screening quantitative polymerase chain reaction (qPCR), we designed and implemented for routine diagnosis a specific duplex qPCR. Finally, structural analysis of the recombinant spike suggested its hybrid content could optimize viral binding to the host cell membrane. These findings prompt further studies of the virological, epidemiological, and clinical features of this recombinant.

First cases of infection with the 21L/BA.2 Omicron variant in Marseille, France.

The SARS-CoV-2 21K/BA.1, 21L/BA.2, and BA.3 Omicron variants have recently emerged worldwide. To date, the 21L/BA.2 Omicron variant has remained very minority globally but became predominant in Denmark instead of the 21K/BA.1 variant. Here, we describe the first cases diagnosed with this variant in south-eastern France. We identified 13 cases using variant-specific qPCR and next-generation sequencing between 28/11/2021 and 31/01/2022, the first two cases being diagnosed in travelers returning from Tanzania. Overall, viral genomes displayed a mean (±standard deviation) number of 65.9 ± 2.5 (range, 61-69) nucleotide substitutions and 31.0 ± 8.3 (27-50) nucleotide deletions, resulting in 49.6 ± 2.2 (45-52) amino acid substitutions (including 28 in the spike protein) and 12.4 ± 1.1 (12-15) amino acid deletions. Phylogeny showed the distribution in three different clusters of these genomes, which were most closely related to genomes from England and South Africa, from Singapore and Nepal, or from France and Denmark. Structural predictions highlighted a significant enlargement and flattening of the surface of the 21L/BA.2 N-terminal domain of the spike protein compared to that of the 21K/BA.1 Omicron variant, which may facilitate initial viral interactions with lipid rafts. Close surveillance is needed at global, country, and center scales to monitor the incidence and clinical outcome of the 21L/BA.2 Omicron variant.

Emergence in southern France of a new SARS-CoV-2 variant harbouring both N501Y and E484K substitutions in the spike protein.

SARS-CoV-2 variants have become a major virological, epidemiological, and clinical concern, particularly with regard to the risk of escape from vaccine-induced immunity. Here, we describe the emergence of a new variant, with the index case returning from travel in Cameroon. For 13 SARS-CoV-2-positive patients living in the same geographical area of southeastern France, a qPCR test for screening variant-associated mutations showed an atypical combination. The genome sequences were obtained by next-generation sequencing with Oxford Nanopore Technologies on GridION instruments within about 8 h. Analysis revealed 46 nucleotide substitutions and 37 deletions, resulting in 30 amino acid substitutions and 12 deletions. Fourteen of the amino acid substitutions, including N501Y and E484K, and nine deletions are located in the spike protein. This genotype pattern led to the establishment of a new Pangolin lineage, named B.1.640.2, that is a phylogenetic sister group to the old B.1.640 lineage, which has now been renamed B.1.640.1. The lineages differ by 25 nucleotide substitutions and 33 deletions. The combination of mutations in these isolates and their phylogenetic position indicate, based on our previous definition, that they represent a new variant, which we have named "IHU". These data are a further example of the unpredictability of the emergence of SARS-CoV-2 variants, and of their possible introduction into a given geographical area from abroad.

Limited spread of a rare spike E484K-harboring SARS-CoV-2 in Marseille, France.

We detected SARS-CoV-2 of PANGO lineage R.1 with the spike substitution E484K in three patients. Eleven other sequences in France and 8,831 worldwide were available from GISAID, 92% originating from Japan. The three genome sequences from our institute were phylogenetically closest to another from Guinea-Conakry, where one of the patients had travelled. These viruses did not exhibit any unusual features in cell culture. Spike structural predictions indicated a 1.3-time higher transmissibility index than for the globally spread B.1.1.7 variant but also an affinity loss for gangliosides that might have slowed dissemination. The spread of new SARS-CoV-2 mutants/variants is still not well understood and therefore difficult to predict, and this hinders implementation of effective preventive measures, including adapted vaccines.