Early Introduction and Community Transmission of SARS-CoV-2 Omicron Variant, New York, New York, USA.

The Omicron variant of SARS-CoV-2 has become dominant in most countries and has raised significant global health concerns. As a global commerce center, New York, New York, USA, constantly faces the risk for multiple variant introductions of SARS-CoV-2. To elucidate the introduction and transmission of the Omicron variant in the city of New York, we created a comprehensive genomic and epidemiologic analysis of 392 Omicron virus specimens collected during November 25-December 11, 2021. We found evidence of 4 independent introductions of Omicron subclades, including the Omicron subclade BA.1.1 with defining substitution of R346K in the spike protein. The continuous genetic divergence within each Omicron subclade revealed their local community transmission and co-circulation in New York, including both household and workplace transmissions supported by epidemiologic evidence. Our study highlights the urgent need for enhanced genomic surveillance and effective response planning for better prevention and management of emerging SARS-CoV-2 variants.

sitePath: a visual tool to identify polymorphism clades and help find fixed and parallel mutations.

BACKGROUND: Identifying polymorphism clades on phylogenetic trees could help detect punctual mutations that are associated with viral functions. With visualization tools coloring the tree, it is easy to visually find clades where most sequences have the same polymorphism state. However, with the fast accumulation of viral sequences, a computational tool to automate this process is urgently needed. RESULTS: Here, by implementing a branch-and-bound-like search method, we developed an R package named sitePath to identify polymorphism clades automatically. Based on the identified polymorphism clades, fixed and parallel mutations could be inferred. Furthermore, sitePath also integrated visualization tools to generate figures of the calculated results. In an example with the influenza A virus H3N2 dataset, the detected fixed mutations coincide with antigenic shift mutations. The highly specificity and sensitivity of sitePath in finding fixed mutations were achieved for a range of parameters and different phylogenetic tree inference software. CONCLUSIONS: The result suggests that sitePath can identify polymorphism clades per site. The clustering of sequences on a phylogenetic tree can be used to infer fixed and parallel mutations. High-quality figures of the calculated results could also be generated by sitePath.

Potential intervariant and intravariant recombination of Delta and Omicron variants.

Among numerous severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concerns, Omicron is more infectious and immune-escaping, while Delta is more pathogenic. Here, we provide evidence for both intervariant and intravariant recombination of the rapidly evolving new SARS-CoV-2 genomes, including XD/XE/XF and BA.3, raising concerns of potential more infectious, immune-escaping, and disease-causing Omicron and Delta-Omicron variants.

SARS-CoV-2 Continuous Genetic Divergence and Changes in Multiplex RT-PCR Detection Pattern on Positive Retesting Median 150 Days after Initial Infection.

Being in the epicenter of the COVID-19 pandemic, our lab tested 193,054 specimens for SARS-CoV-2 RNA by diagnostic multiplex reverse transcription polymerase chain reaction (mRT-PCR) starting in March 2020, of which 17,196 specimens resulted positive. To investigate the dynamics of virus molecular evolution and epidemiology, whole genome amplification (WGA) and Next Generation Sequencing (NGS) were performed on 9516 isolates. 7586 isolates with a high quality were further analyzed for the mutation frequency and spectrum. Lastly, we evaluated the utility of the mRT-PCR detection pattern among 26 reinfected patients with repeat positive testing three months after testing negative from the initial infection. Our results show a continuation of the genetic divergence in viral genomes. Furthermore, our results indicate that independent mutations in the primer and probe regions of the nucleocapsid gene amplicon and envelope gene amplicon accumulate over time. Some of these mutations correlate with the changes of detection pattern of viral targets of mRT-PCR. Our data highlight the significance of a continuous genetic divergence on a gene amplification-based assay, the value of the mRT-PCR detection pattern for complementing the clinical diagnosis of reinfection, and the potential for WGA and NGS to identify mutation hotspots throughout the entire viral genome to optimize the design of the PCR-based gene amplification assay.

Web Resources for SARS-CoV-2 Genomic Database, Annotation, Analysis and Variant Tracking.

The SARS-CoV-2 genomic data continue to grow, providing valuable information for researchers and public health officials. Genomic analysis of these data sheds light on the transmission and evolution of the virus. To aid in SARS-CoV-2 genomic analysis, many web resources have been developed to store, collate, analyze, and visualize the genomic data. This review summarizes web resources used for the SARS-CoV-2 genomic epidemiology, covering data management and sharing, genomic annotation, analysis, and variant tracking. The challenges and further expectations for these web resources are also discussed. Finally, we highlight the importance and need for continued development and improvement of related web resources to effectively track the spread and understand the evolution of the virus.

Genetic differentiation and diversity of SARS-CoV-2 Omicron variant in its early outbreak.

The recently emerged Omicron variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has quickly spread around the world. Although many consensus mutations of the Omicron variant have been recognized, little is known about its genetic variation during its transmission in the population. Here, we comprehensively analyzed the genetic differentiation and diversity of the Omicron variant during its early outbreak. We found that Omicron achieved more structural variations, especially deletions, on the SARS-CoV-2 genome than the other four variants of concern (Alpha, Beta, Gamma, and Delta) in the same timescale. In addition, the Omicron variant acquired, except for 50 consensus mutations, seven great new non-synonymous nucleotide substitutions during its spread. Three of them are on the S protein, including S_A701V, S_L1081V, and S_R346K, which belong to the receptor-binding domain (RBD). The Omicron BA.1 branch could be divided into five divergent groups spreading across different countries and regions based on these seven novel mutations. Furthermore, we found that the Omicron variant possesses more mutations related to a faster transmission rate than the other SARS-CoV-2 variants by assessing the relationship between the genetic diversity and transmission rate. The findings indicated that more attention should be paid to the significant genetic differentiation and diversity of the Omicron variant for better disease prevention and control.

Detecting Potentially Adaptive Mutations from the Parallel and Fixed Patterns in SARS-CoV-2 Evolution.

Early identification of adaptive mutations could provide timely help for the control and prevention of the COVID-19 pandemic. The fast accumulation of SARS-CoV-2 sequencing data provides important support, while also raising a great challenge for the recognition of adaptive mutations. Here, we proposed a computational strategy to detect potentially adaptive mutations from their fixed and parallel patterns in the phylogenetic trajectory. We found that the biological meanings of fixed substitution and parallel mutation are highly complementary, and can reasonably be integrated as a fixed and parallel (paraFix) mutation, to identify potentially adaptive mutations. Tracking the dynamic evolution of SARS-CoV-2, 37 sites in spike protein were identified as having experienced paraFix mutations. Interestingly, 70% (26/37) of them have already been experimentally confirmed as adaptive mutations. Moreover, most of the mutations could be inferred as paraFix mutations one month earlier than when they became regionally dominant. Overall, we believe that the concept of paraFix mutations will help researchers to identify potentially adaptive mutations quickly and accurately, which will provide invaluable clues for disease control and prevention.

covSampler: A subsampling method with balanced genetic diversity for large-scale SARS-CoV-2 genome data sets.

Phylogenetic analysis has been widely used to describe, display, and infer the evolutionary patterns of viruses. The unprecedented accumulation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes has provided valuable materials for the real-time study of SARS-CoV-2 evolution. However, the large number of SARS-CoV-2 genome sequences also poses great challenges for data analysis. Several methods for subsampling these large data sets have been introduced. However, current methods mainly focus on the spatiotemporal distribution of genomes without considering their genetic diversity, which might lead to post-subsampling bias. In this study, a subsampling method named covSampler was developed for the subsampling of SARS-CoV-2 genomes with consideration of both their spatiotemporal distribution and their genetic diversity. First, covSampler clusters all genomes according to their spatiotemporal distribution and genetic variation into groups that we call divergent pathways. Then, based on these divergent pathways, two kinds of subsampling strategies, representative subsampling and comprehensive subsampling, were provided with adjustable parameters to meet different users' requirements. Our performance and validation tests indicate that covSampler is efficient and stable, with an abundance of options for user customization. Overall, our work has developed an easy-to-use tool and a webserver (https://www.covsampler.net) for the subsampling of SARS-CoV-2 genome sequences.

Genomic evidence for divergent co-infections of co-circulating SARS-CoV-2 lineages.

Co-infection of RNA viruses may contribute to their recombination and cause severe clinical symptoms. However, the tracking and identification of SARS-CoV-2 co-infection persist as challenges. Due to the lack of methods for detecting co-infected samples in a large amount of deep sequencing data, the lineage composition, spatial-temporal distribution, and frequency of SARS-CoV-2 co-infection events in the population remains unclear. Here, we propose a hypergeometric distribution-based method named Cov2Coinfect with the ability to decode the lineage composition from 50,809 deep sequencing data. By resolving the mutational patterns in each sample, Cov2Coinfect can precisely determine the co-infected SARS-CoV-2 variants from deep sequencing data. Results from two independent and parallel projects in the United States achieved a similar co-infection rate of 0.3-0.5 % in SARS-CoV-2 positive samples. Notably, all co-infected variants were highly consistent with the co-circulating SARS-CoV-2 lineages in the regional epidemiology, demonstrating that the co-circulation of different variants is an essential prerequisite for co-infection. Overall, our study not only provides a robust method to identify the co-infected SARS-CoV-2 variants from sequencing samples, but also highlights the urgent need to pay more attention to co-infected patients for better disease prevention and control.

MicroRNA-148a Acts as a Tumor Suppressor in Osteosarcoma via Targeting Rho-Associated Coiled-Coil Kinase.

MicroRNAs (miRs) have been demonstrated to be involved in the development and progression of osteosarcoma (OS), but the molecular mechanism still remains to be fully investigated. The present study investigated the function of miR-148a in OS, as well as its underlying mechanism. Our data showed that miR-148a was significantly downregulated in OS tissues compared to their matched adjacent normal tissues, and also in OS cell lines compared to normal human osteoblast cells. Low expression of miR-148a was significantly associated with tumor progression and a poor prognosis for OS patients. Rho-associated coiled-coil kinase 1 (ROCK1) was then identified as a target of miR-148a in Saos-2 and U2OS cells, and the expression of ROCK1 was significantly increased in OS tissues and cell lines. Moreover, the protein expression of ROCK1 was markedly reduced in miR-148a-overexpressing Saos-2 and U2OS cells, but significantly increased in miR-148a-downregulated Saos-2 and U2OS cells. Further investigation indicated that miR-148a had a suppressive effect on the proliferative, migratory, and invasive capacities of Saos-2 and U2OS cells. Moreover, overexpression of ROCK1 attenuated the inhibitory effects of miR-148a upregulation on the malignant phenotypes of Saos-2 and U2OS cells. In addition, overexpression of miR-148a significantly inhibited the tumor growth of U2OS cells in nude mice. Taken together, these data demonstrate that miR-148a acts as a tumor suppressor in OS, at least partly, via targeting ROCK1. Therefore, the miR-148a/ROCK1 axis may become a potential therapeutic target for OS.