Genomic characterization of Influenza A (H1N1)pdm09 and SARS-CoV-2 from Influenza Like Illness (ILI) and Severe Acute Respiratory Illness (SARI) cases reported between July-December, 2022.

Influenza Like Illness (ILI) and Severe Acute Respiratory Infection (SARI) cases are more prone to Influenza and SARS-CoV-2 infection. Accordingly, we genetically characterized Influenza and SARS-CoV-2 in 633 ILI and SARI cases by rRT-PCR and WGS. ILI and SARI cases showed H1N1pdm09 prevalence of 20.9% and 23.2% respectively. 135 (21.3%) H1N1pdm09 and 23 (3.6%) H3N2 and 5 coinfection (0.78%) of H1N1pdm09 and SARS-CoV-2 were detected. Phylogenetic analysis revealed H1N1pdm09 resemblance to clade 6B.1A.5a.2 and their genetic relatedness to InfA/Perth/34/2020, InfA/Victoria/88/2020 and InfA/Victoria/2570/2019. Pan 24 HA and 26 NA nonsynonymous mutations and novel HA (G6D, Y7F, Y78H, P212L, G339R, T508K and S523T) and NA (S229A) mutations were observed. S74R, N129D, N156K, S162N, K163Q and S164T alter HA Cb and Sa antibody recognizing site. Similarly, M19T, V13T substitution and multiple mutations in transmembrane and NA head domain drive antigenic drift. SARS-CoV-2 strains genetically characterized to Omicron BA.2.75 lineage containing thirty nonsynonymous spike mutations exhibited enhanced virulence and transmission rates. Coinfection although detected very minimal, the mutational changes in H1N1pdm09 and SARS-CoV-2 virus infected individuals could alter antibody receptor binding sites, allowing the viruses to escape immune response resulting in better adaptability and transmission. Thus continuous genomic surveillance is required to tackle any future outbreak.

Exploring the diversity and evolutionary strategies of prophages in Hyphomicrobiales, comparing animal-associated with non-animal-associated bacteria.

The Hyphomicrobiales bacterial order (previously Rhizobiales) exhibits a wide range of lifestyle characteristics, including free-living, plant-association, nitrogen-fixing, and association with animals (Bartonella and Brucella). This study explores the diversity and evolutionary strategies of bacteriophages within the Hyphomicrobiales order, comparing animal-associated (AAB) with non-animal-associated bacteria (NAAB). We curated 560 high-quality complete genomes of 58 genera from this order and used the PHASTER server for prophage annotation and classification. For 19 genera with representative genomes, we curated 96 genomes and used the Defense-Finder server to summarize the type of anti-phage systems (APS) found in this order. We analyzed the genetic repertoire and length distributions of prophages, estimating evolutionary rates and comparing intact, questionable, and incomplete prophages in both groups. Analyses of best-fit parameters and bootstrap sensitivity were used to understand the evolutionary processes driving prophage gene content. A total of 1860 prophages distributed in Hyphomicrobiales were found, 695 in AAB and 1165 in the NAAB genera. The results revealed a similar number of prophages per genome in AAB and NAAB and a similar length distribution, suggesting shared mechanisms of genetic acquisition of prophage genes. Changes in the frequency of specific gene classes were observed between incomplete and intact prophages, indicating preferential loss or enrichment in both groups. The analysis of best-fit parameters and bootstrap sensitivity tests indicated a higher selection coefficient, induction rate, and turnover in NAAB genomes. We found 68 types of APS in Hyphomicrobiales; restriction modification (RM) and abortive infection (Abi) were the most frequent APS found for all Hyphomicrobiales, and within the AAB group. This classification of APS showed that NAAB genomes have a greater diversity of defense systems compared to AAB, which could be related to the higher rates of prophage induction and turnover in the latter group. Our study provides insights into the distributions of both prophages and APS in Hyphomicrobiales genomes, demonstrating that NAAB carry more defense systems against phages, while AAB show increased prophage stability and an increased number of incomplete prophages. These results suggest a greater role for domesticated prophages within animal-associated bacteria in Hyphomicrobiales.

Whole genome sequence analysis of shallot virus X from India reveals it to be a natural recombinant with positive selection pressure.

BACKGROUND: Shallots are infected by various viruses like Onion yellow dwarf virus (OYDV), Leek yellow stripe virus (LYSV), Shallot latent virus (SLV) and Shallot virus X (ShVX). In India, they have been found to be persistently infected by ShVX. ShVX also infects onion and garlic in combination with other carlaviruses and potyviruses. ShVX is a member of genus Allexivirus of family Alphaflexiviridae. ShVX has a monopartite genome, which is represented by positive sense single-stranded RNA. Globally, only six complete and 3 nearly complete genome sequences of ShV X are reported to date. This number is insufficient to measure a taxon's true molecular diversity. Moreover, the complete genome sequence of ShVX from Asia has not been reported as yet. Therefore, this study was undertaken to generate a complete genome sequence of ShVX from India. RESULTS: Shallot virus X (ShVX) is one of the significant threats to Allium crop production. In this study, we report the first complete genome sequence of the ShVX from India through Next-generation sequencing (NGS). The complete genome of the ShVX (Accession No. OK104171), from this study comprised 8911 nucleotides. In-silico analysis of the sequence revealed variability between this isolate and isolates from other countries. The dissimilarities are spread all over the genome specifically some non-coding intergenic regions. Statistical analysis of individual genes for site-specific selection indicates a positive selection in NABP region. The presence of a recombination event was detected in coat protein region. The sequence similarity percentage and phylogenetic analysis indicate ShVX Indian isolate is a distinctly different isolate. Recombination and site-specific selection may have a function in the evolution of this isolate. This is the first detailed study of the ShVX complete genome sequence from Southeast Asia. CONCLUSION: This study presents the first report of the entire genome sequence of an Indian isolate of ShVX along with an in-depth exploration of its evolutionary traits. The findings highlight the Indian variant as a naturally occurring recombinant, emphasizing the substantial role of recombination in the evolution of this viral species. This insight into the molecular diversity of strains within a specific geographical region holds immense significance for comprehending and forecasting potential epidemics. Consequently, the insights garnered from this research hold practical value for shaping ShVX management strategies and providing a foundation for forthcoming studies delving into its evolutionary trajectory.

Long-read powered viral metagenomics in the oligotrophic Sargasso Sea.

Dominant microorganisms of the Sargasso Sea are key drivers of the global carbon cycle. However, associated viruses that shape microbial community structure and function are not well characterised. Here, we combined short and long read sequencing to survey Sargasso Sea phage communities in virus- and cellular fractions at viral maximum (80 m) and mesopelagic (200 m) depths. We identified 2,301 Sargasso Sea phage populations from 186 genera. Over half of the phage populations identified here lacked representation in global ocean viral metagenomes, whilst 177 of the 186 identified genera lacked representation in genomic databases of phage isolates. Viral fraction and cell-associated viral communities were decoupled, indicating viral turnover occurred across periods longer than the sampling period of three days. Inclusion of long-read data was critical for capturing the breadth of viral diversity. Phage isolates that infect the dominant bacterial taxa Prochlorococcus and Pelagibacter, usually regarded as cosmopolitan and abundant, were poorly represented.

A novel birnavirus identified as the causative agent of summer atrophy of pearl oyster (Pinctada fucata (Gould)).

The Akoya pearl oyster (Pinctada fucata (Gould)) is the most important species for pearl cultivation in Japan. Mass mortality of 0-year-old juvenile oysters and anomalies in adults, known as summer atrophy, have been observed in major pearl farming areas during the season when seawater temperatures exceed about 20 °C since 2019. In this study, we identified a novel birnavirus as the pathogen of summer atrophy and named it Pinctada birnavirus (PiBV). PiBV was first presumed to be the causative agent when it was detected specifically and frequently in the infected oysters in a comparative metatranscriptomics of experimentally infected and healthy pearl oysters. Subsequently, the symptoms of summer atrophy were reproduced by infection tests using purified PiBV. Infection of juvenile oysters with PiBV resulted in an increase in the PiBV genome followed by the atrophy of soft body and subsequent mortality. Immunostaining with a mouse antiserum against a recombinant PiBV protein showed that the virus antigen was localized mainly in the epithelial cells on the outer surface of the mantle. Although the phylogenetic analysis using maximum likelihood method placed PiBV at the root of the genus Entomobirnavirus, the identity of the bi-segmented, genomic RNA to that of known birnaviruses at the full-length amino acid level was low, suggesting that PiBV forms a new genus. The discovery of PiBV will be the basis for research to control this emerging disease.

Molecular characterization of emerging Echovirus 11 (E11) shed light on the recombinant origin of a variant associated with severe hepatitis in neonates.

Echovirus 11 (E11) has gained attention owing to its association with severe neonatal infections. Due to the limited data available, the World Health Organization (WHO) considers public health risk to the general population to be low. The present study investigated the genetic variation and molecular evolution of E11 genomes collected from May to December 2023. Whole genome sequencing (WGS) was performed for 16 E11 strains. Phylogenetic analysis on WG showed how all Italian strains belonged to genogroup D5, similarly to other E11 strains recently reported in France and Germany all together aggregated into separate clusters. A cluster-specific recombination pattern was also identified using phylogenetic analysis of different genome regions. Echovirus 6 was identified as the major recombinant virus in 3Cpro and 3Dpol regions. The molecular clock analysis revealed that the recombination event probably occurred in June 2018 (95% HPD interval: Jan 2016-Jan 2020). Shannon entropy analyses, within P1 region, showed how 11 amino acids exhibited relatively high entropy. Five of them were exposed on the canyon region which is responsible for receptor binding with the neonatal Fc receptor. The present study showed the recombinant origin of a new lineage of E11 associated with severe neonatal infections.

Prevalence and diversity of imported severe acute respiratory syndrome coronavirus 2 variants in China, 2021-2022.

The causative agent of coronavirus disease 2019 (COVID-19), known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has spread accumulatively to 240 countries and continues to evolve. To gain a comprehensive understanding of the epidemiological characteristics of imported variants in China and their correlation with global circulating variants, genomic surveillance data from 11 139 imported COVID-19 cases submitted by Chinese provincial CDC laboratories between 2021 and 2022 were analyzed. Consensus sequences underwent rigorous quality checks, followed by amino acid mutations analysis using Nextclade. Sequences with satisfactory quality control status were classified according to the Pango nomenclature. The results showed that the dominant variants in imported cases reflected the global epidemic trend. An increase in the number of imported SARS-CoV-2 lineages monitored in China in the second half of 2022, and the circulating Omicron subvariants changed from the ancestral lineages of BA.5 and BA.2 into the lineages containing key amino acid mutations of spike protein. There was significant variation in the detection of Omicron subvariants among continents (χ2 = 321.968, p < 0.001) in the second half of 2022, with four lineages (BA.2.3.7, BA.2.2, BA.5.2.7, and XBB.1.2) identified through imported surveillance mainly prevalent respectively in Taiwan, China, Hong Kong SAR, China, Russian Federation, and Singapore. These findings revealed the alterations in circulating imported variants from 2021 to 2022 in China, reflecting the higher diversity of lineages in the second half of 2022, and revealed the predominant lineages of countries or regions that are in close contacts to China, providing new insights into the global prevalence of SARS-CoV-2.

Molecular detection and genetic characteristics of porcine circovirus 3 and porcine circovirus 4 in central China.

Porcine circoviruses (PCVs) are a significant cause of concern for swine health, with four genotypes currently recognized. Two of these, PCV3 and PCV4, have been detected in pigs across all age groups, in both healthy and diseased animals. These viruses have been associated with various clinical manifestations, including porcine dermatitis and nephropathy syndrome (PDNS) and respiratory and enteric signs. In this study, we detected PCV3 and PCV4 in central China between January 2022 and February 2023. We tested fecal swabs and tissue samples from growing-finishing and suckling pigs with or without respiratory and systemic manifestations and found the prevalence of PCV3 to be 15.15% (15/99) and that of PCV3/PCV4 coinfection to be 4.04% (4/99). This relatively low prevalence might be attributed to the fact that most of the clinical samples were collected from pigs exhibiting respiratory signs, with only a few samples having been obtained from pigs with diarrhea. In some cases, PCV2 was also detected, and the coinfection rates of PCV2/3, PCV2/4, and PCV2/3/4 were 6.06% (6/99), 5.05% (5/99), and 3.03% (3/99), respectively. The complete genomic sequences of four PCV3 and two PCV4 isolates were determined. All four of the PCV3 isolates were of subtype PCV3b, and the two PCV4 isolates were of subtype PCV4b. Two mutations (A24V and R27K) were found in antibody recognition domains of PCV3, suggesting that they might be associated with immune escape. This study provides valuable insights into the molecular epidemiology and evolution of PCV3 and PCV4 that will be useful in future investigations of genotyping, immunogenicity, and immune evasion strategies.

Nanopore adaptive sampling of a metagenomic sample derived from a human monkeypox case.

In 2022, a series of human monkeypox cases in multiple countries led to the largest and most widespread outbreak outside the known endemic areas. Setup of proper genomic surveillance is of utmost importance to control such outbreaks. To this end, we performed Nanopore (PromethION P24) and Illumina (NextSeq. 2000) Whole Genome Sequencing (WGS) of a monkeypox sample. Adaptive sampling was applied for in silico depletion of the human host genome, allowing for the enrichment of low abundance viral DNA without a priori knowledge of sample composition. Nanopore sequencing allowed for high viral genome coverage, tracking of sample composition during sequencing, strain determination, and preliminary assessment of mutational pattern. In addition to that, only Nanopore data allowed us to resolve the entire monkeypox virus genome, with respect to two structural variants belonging to the genes OPG015 and OPG208. These SVs in important host range genes seem stable throughout the outbreak and are frequently misassembled and/or misannotated due to the prevalence of short read sequencing or short read first assembly. Ideally, standalone standard Illumina sequencing should not be used for Monkeypox WGS and de novo assembly, since it will obfuscate the structure of the genome, which has an impact on the quality and completeness of the genomes deposited in public databases and thus possibly on the ability to evaluate the complete genetic reason for the host range change of monkeypox in the current pandemic.

Genomic characterization of human respiratory syncytial virus circulating in Islamabad, Pakistan, during an outbreak in 2022-2023.

In this study, conducted at the National Institute of Health, Islamabad, during an outbreak of human respiratory syncytial virus (hRSV) from December 2022 to January 2023, the first whole-genome sequences of hRSV isolates from Islamabad, Pakistan, were determined. Out of 10 positive samples, five were sequenced, revealing the presence of two genotypes: RSV-A (GA2.3.5, ON1 strain) and RSV-B (GB5.0.5.a, BA-10 strain). A rare non-synonymous substitution (E232G) in G the protein and N276S in the F protein were found in RSV-A. In RSV-B, the unique mutations K191R, Q209R, and I206M were found in the F protein. These mutations could potentially influence vaccine efficacy and viral pathogenicity. This research underscores the importance of genomic surveillance for understanding RSV diversity and guiding public health responses in Pakistan.

Substantial viral diversity in bats and rodents from East Africa: insights into evolution, recombination, and cocirculation.

BACKGROUND: Zoonotic viruses cause substantial public health and socioeconomic problems worldwide. Understanding how viruses evolve and spread within and among wildlife species is a critical step when aiming for proactive identification of viral threats to prevent future pandemics. Despite the many proposed factors influencing viral diversity, the genomic diversity and structure of viral communities in East Africa are largely unknown. RESULTS: Using 38.3 Tb of metatranscriptomic data obtained via ultradeep sequencing, we screened vertebrate-associated viromes from 844 bats and 250 rodents from Kenya and Uganda collected from the wild. The 251 vertebrate-associated viral genomes of bats (212) and rodents (39) revealed the vast diversity, host-related variability, and high geographic specificity of viruses in East Africa. Among the surveyed viral families, Coronaviridae and Circoviridae showed low host specificity, high conservation of replication-associated proteins, high divergence among viral entry proteins, and frequent recombination. Despite major dispersal limitations, recurrent mutations, cocirculation, and occasional gene flow contribute to the high local diversity of viral genomes. CONCLUSIONS: The present study not only shows the landscape of bat and rodent viromes in this zoonotic hotspot but also reveals genomic signatures driven by the evolution and dispersal of the viral community, laying solid groundwork for future proactive surveillance of emerging zoonotic pathogens in wildlife. Video Abstract.

Genomic analyses of Symbiomonas scintillans show no evidence for endosymbiotic bacteria but does reveal the presence of giant viruses.

Symbiomonas scintillans Guillou et Chrétiennot-Dinet, 1999 is a tiny (1.4 µm) heterotrophic microbial eukaryote. The genus was named based on the presence of endosymbiotic bacteria in its endoplasmic reticulum, however, like most such endosymbionts neither the identity nor functional association with its host were known. We generated both amplification-free shotgun metagenomics and whole genome amplification sequencing data from S. scintillans strains RCC257 and RCC24, but were unable to detect any sequences from known lineages of endosymbiotic bacteria. The absence of endobacteria was further verified with FISH analyses. Instead, numerous contigs in assemblies from both RCC24 and RCC257 were closely related to prasinoviruses infecting the green algae Ostreococcus lucimarinus, Bathycoccus prasinos, and Micromonas pusilla (OlV, BpV, and MpV, respectively). Using the BpV genome as a reference, we assembled a near-complete 190 kbp draft genome encoding all hallmark prasinovirus genes, as well as two additional incomplete assemblies of closely related but distinct viruses from RCC257, and three similar draft viral genomes from RCC24, which we collectively call SsVs. A multi-gene tree showed the three SsV genome types branched within highly supported clades with each of BpV2, OlVs, and MpVs, respectively. Interestingly, transmission electron microscopy also revealed a 190 nm virus-like particle similar the morphology and size of the endosymbiont originally reported in S. scintillans. Overall, we conclude that S. scintillans currently does not harbour an endosymbiotic bacterium, but is associated with giant viruses.

Early detection of emerging viral variants through analysis of community structure of coordinated substitution networks.

The emergence of viral variants with altered phenotypes is a public health challenge underscoring the need for advanced evolutionary forecasting methods. Given extensive epistatic interactions within viral genomes and known viral evolutionary history, efficient genomic surveillance necessitates early detection of emerging viral haplotypes rather than commonly targeted single mutations. Haplotype inference, however, is a significantly more challenging problem precluding the use of traditional approaches. Here, using SARS-CoV-2 evolutionary dynamics as a case study, we show that emerging haplotypes with altered transmissibility can be linked to dense communities in coordinated substitution networks, which become discernible significantly earlier than the haplotypes become prevalent. From these insights, we develop a computational framework for inference of viral variants and validate it by successful early detection of known SARS-CoV-2 strains. Our methodology offers greater scalability than phylogenetic lineage tracing and can be applied to any rapidly evolving pathogen with adequate genomic surveillance data.

Molecular characterization and host reaction to tomato mottle mosaic virus isolated from sweet pepper seeds in Japan.

Many countries have identified tomato mottle mosaic virus (ToMMV) as a serious threat to tomato production. Here, we constructed and characterized infectious clones of ToMMV isolated from Japanese sweet pepper seeds. The genome of the Japanese isolate is 6399 nucleotides in length and exhibits the highest identity with previously characterized isolates. For example, it is 99.7% identical to that of the Mauritius isolate, which occurs worldwide. Phylogenetic analysis based on complete genome sequences revealed that the Japanese isolates clustered in the same clade as those from other countries. When homozygous tomato cultivars with tobamovirus resistance genes were inoculated with an infectious cDNA clone of ToMMV, the virus systemically infected tomato plants with symptoms typical of Tm-1-carrying tomato cultivars. In contrast, tomato cultivars carrying Tm-2 or Tm-22 showed symptoms only on the inoculated leaves. Furthermore, when commercial cultivars of Tm-22 heterozygous tomato were inoculated with ToMMV, systemic infections were observed in all cultivars, with infection frequencies ranging from 25 to 100%. Inoculation of heterozygous sweet pepper cultivars with tobamovirus resistance genes (L1, L3, and L4) with ToMMV resulted in an infection frequency of about 70%, but most of the infected L1, L3, and L4 cultivars were symptomless, and 10-20% showed symptoms of necrosis and yellowing. Tomato mosaic virus strain L11A, an attenuated virus, did not provide cross-protection against ToMMV and led to systemic infection with typical symptoms. These results suggest that ToMMV might cause extensive damage to existing tomato and sweet pepper cultivars commonly grown in Japan.

Locations and structures of influenza A virus packaging-associated signals and other functional elements via an in silico pipeline for predicting constrained features in RNA viruses.

Influenza A virus contains regions of its segmented genome associated with ability to package the segments into virions, but many such regions are poorly characterised. We provide detailed predictions of the key locations within these packaging-associated regions, and their structures, by applying a recently-improved pipeline for delineating constrained regions in RNA viruses and applying structural prediction algorithms. We find and characterise other known constrained regions within influenza A genomes, including the region associated with the PA-X frameshift, regions associated with alternative splicing, and constraint around the initiation motif for a truncated PB1 protein, PB1-N92, associated with avian viruses. We further predict the presence of constrained regions that have not previously been described. The extra characterisation our work provides allows investigation of these key regions for drug target potential, and points towards determinants of packaging compatibility between segments.

RNA genome packaging and capsid assembly of bluetongue virus visualized in host cells.

Unlike those of double-stranded DNA (dsDNA), single-stranded DNA (ssDNA), and ssRNA viruses, the mechanism of genome packaging of dsRNA viruses is poorly understood. Here, we combined the techniques of high-resolution cryoelectron microscopy (cryo-EM), cellular cryoelectron tomography (cryo-ET), and structure-guided mutagenesis to investigate genome packaging and capsid assembly of bluetongue virus (BTV), a member of the Reoviridae family of dsRNA viruses. A total of eleven assembly states of BTV capsid were captured, with resolutions up to 2.8 Å, with most visualized in the host cytoplasm. ATPase VP6 was found underneath the vertices of capsid shell protein VP3 as an RNA-harboring pentamer, facilitating RNA packaging. RNA packaging expands the VP3 shell, which then engages middle- and outer-layer proteins to generate infectious virions. These revealed "duality" characteristics of the BTV assembly mechanism reconcile previous contradictory co-assembly and core-filling models and provide insights into the mysterious RNA packaging and capsid assembly of Reoviridae members and beyond.

A giant virus infecting the amoeboflagellate Naegleria.

Giant viruses (Nucleocytoviricota) are significant lethality agents of various eukaryotic hosts. Although metagenomics indicates their ubiquitous distribution, available giant virus isolates are restricted to a very small number of protist and algal hosts. Here we report on the first viral isolate that replicates in the amoeboflagellate Naegleria. This genus comprises the notorious human pathogen Naegleria fowleri, the causative agent of the rare but fatal primary amoebic meningoencephalitis. We have elucidated the structure and infection cycle of this giant virus, Catovirus naegleriensis (a.k.a. Naegleriavirus, NiV), and show its unique adaptations to its Naegleria host using fluorescence in situ hybridization, electron microscopy, genomics, and proteomics. Naegleriavirus is only the fourth isolate of the highly diverse subfamily Klosneuvirinae, and like its relatives the NiV genome contains a large number of translation genes, but lacks transfer RNAs (tRNAs). NiV has acquired genes from its Naegleria host, which code for heat shock proteins and apoptosis inhibiting factors, presumably for host interactions. Notably, NiV infection was lethal to all Naegleria species tested, including the human pathogen N. fowleri. This study expands our experimental framework for investigating giant viruses and may help to better understand the basic biology of the human pathogen N. fowleri.

Unveiling hidden structural patterns in the SARS-CoV-2 genome: Computational insights and comparative analysis.

SARS-CoV-2, the causative agent of COVID-19, is known to exhibit secondary structures in its 5' and 3' untranslated regions, along with the frameshifting stimulatory element situated between ORF1a and 1b. To identify additional regions containing conserved structures, we utilized a multiple sequence alignment with related coronaviruses as a starting point. We applied a computational pipeline developed for identifying non-coding RNA elements. Our pipeline employed three different RNA structural prediction approaches. We identified forty genomic regions likely to harbor structures, with ten of them showing three-way consensus substructure predictions among our predictive utilities. We conducted intracomparisons of the predictive utilities within the pipeline and intercomparisons with four previously published SARS-CoV-2 structural datasets. While there was limited agreement on the precise structure, different approaches seemed to converge on regions likely to contain structures in the viral genome. By comparing and combining various computational approaches, we can predict regions most likely to form structures, as well as a probable structure or ensemble of structures. These predictions can be used to guide surveillance, prophylactic measures, or therapeutic efforts. Data and scripts employed in this study may be found at https://doi.org/10.5281/zenodo.8298680.

Phylogenetic analysis of the promoter element 2 of paramyxo- and filoviruses.

Paramyxo- and filovirus genomes are equipped with bipartite promoters at their 3' ends to initiate RNA synthesis. The two elements, the primary promoter element 1 (PE1) and the secondary promoter element 2 (PE2), are separated by a spacer region that must be precisely a multiple of 6 nucleotides (nts), indicating these viruses adhere to the "rule of six." However, our knowledge of PE2 has been limited to a narrow spectrum of virus species. In this study, a comparative analysis of 1,647 paramyxoviral genomes from a public database revealed that the paramyxovirus PE2 can be clearly categorized into two distinct subcategories: one marked by C repeats at every six bases (exclusive to the subfamily Orthoparamyxovirinae) and another characterized by CG repeats every 6 nts (observed in the subfamilies Avulavirinae and Rubulavirinae). This unique pattern collectively mirrors the evolutionary lineage of these subfamilies. Furthermore, we showed that PE2 of the Rubulavirinae, with the exception of mumps virus, serves as part of the gene-coding region. This may be due to the fact that the Rubulavirinae are the only paramyxoviruses that cannot propagate without RNA editing. Filoviruses have three to eight consecutive uracil repeats every six bases (UN5) in PE2, which is located in the 3' end region of the genome. We obtained PE2 sequences from 2,195 filoviruses in a public database and analyzed the sequence conservation among virus species. Our results indicate that the continuity of UN5 hexamers is consistently maintained with a high degree of conservation across virus species. IMPORTANCE: The genomic intricacies of paramyxo- and filoviruses are highlighted by the bipartite promoters-promoter element 1 (PE1) and promoter element 2 (PE2)-at their 3' termini. The spacer region between these elements follows the "rule of six," crucial for genome replication. By a comprehensive analysis of paramyxoviral genome sequences, we identified distinct subcategories of PE2 based on C and CG repeats that were specific to Orthoparamyxovirinae and Avulavirinae/Rubulavirinae, respectively, mirroring their evolutionary lineages. Notably, the PE2 of Rubulavirinae is integrated into the gene-coding region, a unique trait potentially linked to its strict dependence on RNA editing for virus growth. This study also focused on the PE2 sequences in filovirus genomes. The strict conservation of the continuity of UN5 among virus species emphasizes its crucial role in viral genome replication.

Optimisation and evaluation of viral genomic sequencing of SARS-CoV-2 rapid diagnostic tests: a laboratory and cohort-based study.

BACKGROUND: Sequencing of SARS-CoV-2 from rapid diagnostic tests (RDTs) can bolster viral genomic surveillance efforts; however, approaches to maximise and standardise pathogen genome recovery from RDTs remain underdeveloped. We aimed to systematically optimise the elution of genetic material from RDT components and to evaluate the efficacy of RDT sequencing for outbreak investigation. METHODS: In this laboratory and cohort-based study we seeded RDTs with inactivated SARS-CoV-2 to optimise the elution of genomic material from RDT lateral flow strips. We measured the effect of changes in buffer type, time in buffer, and rotation on PCR cycle threshold (Ct) value. We recruited individuals older than 18 years residing in the greater Boston area, MA, USA, from July 18 to Nov 5, 2022, via email advertising to students and staff at Harvard University, MA, USA, and via broad social media advertising. All individuals recruited were within 5 days of a positive diagnostic test for SARS-CoV-2; no other relevant exclusion criteria were applied. Each individual completed two RDTs and one PCR swab. On Dec 29, 2022, we also collected RDTs from a convenience sample of individuals who were positive for SARS-CoV-2 and associated with an outbreak at a senior housing facility in MA, USA. We extracted all returned PCR swabs and RDT components (ie, swab, strip, or buffer); samples with a Ct of less than 40 were subject to amplicon sequencing. We compared the efficacy of elution and sequencing across RDT brands and components and used RDT-derived sequences to infer transmission links within the outbreak at the senior housing facility. We conducted metagenomic sequencing of negative RDTs from symptomatic individuals living in the senior housing facility. FINDINGS: Neither elution duration of greater than 10 min nor rotation during elution impacted viral titres. Elution in Buffer AVL (Ct=31·4) and Tris-EDTA Buffer (Ct=30·8) were equivalent (p=0·34); AVL outperformed elution in lysis buffer and 50% lysis buffer (Ct=40·0, p=0·0029 for both) as well as Universal Viral Transport Medium (Ct=36·7, p=0·079). Performance of RDT strips was poorer than that of matched PCR swabs (mean Ct difference 10·2 [SD 4·3], p<0·0001); however, RDT swabs performed similarly to PCR swabs (mean Ct difference 4·1 [5·2], p=0·055). No RDT brand significantly outperformed another. Across sample types, viral load predicted the viral genome assembly length. We assembled greater than 80% complete genomes from 12 of 17 RDT-derived swabs, three of 18 strips, and four of 11 residual buffers. We generated outbreak-associated SARS-CoV-2 genomes using both amplicon and metagenomic sequencing and identified multiple introductions of the virus that resulted in downstream transmission. INTERPRETATION: RDT-derived swabs are a reasonable alternative to PCR swabs for viral genomic surveillance and outbreak investigation. RDT-derived lateral flow strips yield accurate, but significantly fewer, viral reads than matched PCR swabs. Metagenomic sequencing of negative RDTs can identify viruses that might underlie patient symptoms. FUNDING: The National Science Foundation, the Hertz Foundation, the National Institute of General Medical Sciences, Harvard Medical School, the Howard Hughes Medical Institute, the US Centers for Disease Control and Prevention, the Broad Institute and the National Institute of Allergy and Infectious Diseases.