The genomic and phylogenetic analysis of Marseillevirus cajuinensis raises questions about the evolution of Marseilleviridae lineages and their taxonomical organization.

Marseilleviruses (MsV) are a group of viruses that compose the Marseilleviridae family within the Nucleocytoviricota phylum. They have been found in different samples, mainly in freshwater. MsV are classically organized into five phylogenetic lineages (A/B/C/D/E), but the current taxonomy does not fully represent all the diversity of the MsV lineages. Here, we describe a novel strain isolated from a Brazilian saltwater sample named Marseillevirus cajuinensis. Based on genomics and phylogenetic analyses, M. cajuinensis exhibits a 380,653-bp genome that encodes 515 open reading frames. Additionally, M. cajuinensis encodes a transfer RNA, a feature that is rarely described for Marseilleviridae. Phylogeny suggests that M. cajuinensis forms a divergent branch within the MsV lineage A. Furthermore, our analysis suggests that the common ancestor for the five classical lineages of MsV diversified into three major groups. The organization of MsV into three main groups is reinforced by a comprehensive analysis of clusters of orthologous groups, sequence identities, and evolutionary distances considering several MsV isolates. Taken together, our results highlight the importance of discovering new viruses to expand the knowledge about known viruses that belong to the same lineages or families. This work proposes a new perspective on the Marseilleviridae lineages organization that could be helpful to a future update in the taxonomy of the Marseilleviridae family. IMPORTANCE: Marseilleviridae is a family of viruses whose members were mostly isolated from freshwater samples. In this work, we describe the first Marseillevirus isolated from saltwater samples, which we called Marseillevirus cajuinensis. Most of M. cajuinensis genomic features are comparable to other Marseilleviridae members, such as its high number of unknown proteins. On the other hand, M. cajuinensis encodes a transfer RNA, which is a gene category involved in protein translation that is rarely described in this viral family. Additionally, our phylogenetic analyses suggested the existence of, at least, three major Marseilleviridae groups. These observations provide a new perspective on Marseilleviridae lineages organization, which will be valuable in future updates to the taxonomy of the family since the current official classification does not capture all the Marseilleviridae known diversity.

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.

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 question of screening organ donors for hepatitis e virus: a case report of transmission by kidney transplantation in France and a review of the literature.

BACKGROUND: Hepatitis E is a potentially serious infection in organ recipients, with an estimated two-thirds of cases becoming chronic, and with a subsequent risk of cirrhosis and death. In Europe, transmission occurs most often through the consumption of raw or undercooked pork, more rarely through blood transfusion, but also after solid organ transplantation. Here we describe a case of Hepatitis E virus (HEV) infection transmitted following kidney transplantation and review the literature describing cases of HEV infection transmitted by solid organ transplantation. CASE PRESENTATION: Three weeks after kidney transplantation, the patient presented with an isolated minimal increase in GGT and hepatic cytolysis 6 months later, leading to the diagnosis of genotype 3c hepatitis E, with a plasma viral load of 6.5 log10IU/mL. In retrospect, HEV RNA was detected in the patient's serum from the onset of hepatitis, and in the donor's serum on the day of donation, with 100% identity between the viral sequences, confirming donor-derived HEV infection. Hepatitis E had a chronic course, was treated by ribavirin, and relapsed 10 months after the end of treatment. DISCUSSION: Seven cases of transmission of HEV by solid organ transplantation have been described since 2012 without systematic screening for donors, all diagnosed at the chronic infection stage; two patients died. HEV organ donor transmission may be underestimated and there is insufficient focus on immunocompromised patients in whom mild liver function test impairment is potentially related to hepatitis E. However, since HEV infection is potentially severe in these patients, and as evidence accumulates, we believe that systematic screening of organ donors should be implemented for deceased and living donors regardless of liver function abnormalities, as is already the case in the UK and Spain. In January 2024, the French regulatory agency of transplantation has implemented mandatory screening of organ donors for HEV RNA.

Management of immunosuppression in lung transplant recipients and COVID-19 outcomes: an observational retrospective cohort-study.

BACKGROUND: The aim of this study was to assess the impact of immunosuppression management on coronavirus disease 2019 (COVID-19) outcomes. METHODS: We performed a single-center retrospective study in a cohort of 358 lung transplant recipients (LTx) over the period from March 2020 to April 2022. All included symptomatic patients had at least one positive SARS-CoV-2 rt-PCR. We used a composite primary outcome for COVID-19 including increased need for oxygen since the hospital admission, ICU transfer, and in-hospital mortality. We assessed by univariate and multivariate analyses the risk factors for poor outcomes. RESULTS: Overall, we included 91 LTx who contracted COVID-19. The COVID-19 in-hospital mortality rate reached 4.4%. By hierarchical clustering, we found a strong and independent association between the composite poor outcome and the discontinuation of at least one immunosuppressive molecule among tacrolimus, cyclosporine, mycophenolate mofetil, and everolimus. Obesity (OR = 16, 95%CI (1.96; 167), p = 0.01) and chronic renal failure (OR = 4.6, 95%CI (1.4; 18), p = 0.01) were also independently associated with the composite poor outcome. Conversely, full vaccination was protective (OR = 0.23, 95%CI (0.046; 0.89), p = 0.047). CONCLUSION: The administration of immunosuppressive drugs such as tacrolimus, cyclocporine or everolimus can have a protective effect in LTx with COVID-19, probably related to their intrinsic antiviral capacity.

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.

From viral democratic genomes to viral wild bunch of quasispecies.

The tremendous majority of RNA genomes from pathogenic viruses analyzed and deposited in databases are consensus or "democratic" genomes. They represent the genomes most frequently found in the clinical samples of patients but do not account for the huge genetic diversity of coexisting genomes, which is better described as quasispecies. A viral quasispecies is defined as the dynamic distribution of nonidentical but closely related mutants, variants, recombinant, or reassortant viral genomes. Viral quasispecies have collective behavior and dynamics and are the subject of internal interactions that comprise interference, complementation, or cooperation. In the setting of SARS-CoV-2 infection, intrahost SARS-CoV-2 genetic diversity was recently notably reported for immunocompromised, chronically infected patients, for patients treated with monoclonal antibodies targeting the viral spike protein, and for different body compartments of a single patient. A question that deserves attention is whether such diversity is generated postinfection from a clonal genome in response to selection pressure or is already present at the time of infection as a quasispecies. In the present review, we summarize the data supporting that hosts are infected by a "wild bunch" of viruses rather than by multiple virions sharing the same genome. Each virion in the "wild bunch" may have different virulence and tissue tropisms. As the number of viruses replicated during host infections is huge, a viral quasispecies at any time of infection is wide and is also influenced by host-specific selection pressure after infection, which accounts for the difficulty in deciphering and predicting the appearance of more fit variants and the evolution of epidemics of novel RNA viruses.

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.

Scurvy complicated with Capnocytophaga sputigena sepsis as a possible cause of death of king Saint-Louis of France (1270 AD).

The cause of death of Saint-Louis is not known, but recent findings indicated that he presented scurvy and inflammatory jaw disease, which has been associated with infection by oral commensals. Here, we have the exceptional opportunity to analyze the relics of the viscera of King Saint-Louis. A 4.3 g sample from the viscera relics of King Saint-Louis conserved in Versailles' cathedral was subjected to radiocarbon dating, electronic and optic microscopy, and elementary, palynological, molecular, proteomics and microbiological analyses including specific PCR and v3v4 16 S rRNA gene amplification prior to large-scale sequencing using an Illumina MiSeq instrument. The measured radiocarbon age was Cal 1290 CE-1400, which was compatible with that of the viscera of St Louis viscera, considering the addition of lime, incense and vegetables within the human organs. Elemental and palynological analyses confirmed a medieval embalming process. Proteomics analysis identified mainly human muscle and blood proteins. Specific PCR for plague, amoebiasis, shigellosis and typhoid fever was negative. C. sputigena was identified as the main pathogenic species representing 10.8 % of all microbial sequences. In contrast, C. sputigena was found in only 0.001 % of samples sequenced in our center, and the 23 positive human samples showed a dramatically lower abundance (0.02-2.6 %). In the literature, human infections with C. sputigena included odontitis, dental abscess, sinusitis, thoracic infections and bacteremia, particularly in immunocompromised patients with oral and dental diseases consistent with recent analysis of King Saint-Louis' jaw. C. sputigena, a commensal of the mouth that is potentially pathogenic and responsible for fatal bacteremia, may have been the cause of the king's death.

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.

Reinfections with Different SARS-CoV-2 Omicron Subvariants, France.

We describe 188 patients in France who were successively infected with different SARS-CoV-2 Omicron subvariants, including BA.1, BA.2, and BA.5. Time between 2 infections was <90 days for 50 (26.6%) patients and <60 days for 28 (14.9%) patients. This finding suggests that definitions for SARS-CoV-2 reinfection require revision.

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.

Characteristics of the first 1119 SARS-CoV-2 Omicron variant cases, in Marseille, France, November-December 2021.

One thousand one hundred and nineteen cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant cases have been diagnosed at the Institut Hospitalo-Universitaire Méditerranée Infection, Marseille, France, between November 28, 2021, and December 31, 2021. Among the 825 patients with known vaccination status, 383 (46.4%) were vaccinated, of whom 91.9% had received at least two doses of the vaccine. Interestingly, 26.3% of cases developed SARS-CoV-2 infection within 21 days following the last dose of vaccine suggesting possible early production of anti-SARS-CoV-2 facilitating antibodies. Twenty-one patients have been hospitalized, one patient required intensive care, and another patient who received a vaccine booster dose died.

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.