Mutational Analysis of Circulating Omicron SARS-CoV-2 Lineages in the Al-Baha Region of Saudi Arabia.

Purpose: Omicron (B.1.1.529) is one of the highly mutated variants of concern of SARS-CoV-2. Lineages of Omicron bear a remarkable degree of mutations leading to enhanced pathogenicity and upward transmission trajectory. Mutating Omicron lineages may trigger a fresh COVID-19 wave at any time in any region. We aimed at the whole-genome sequencing of SARS-CoV-2 to determine variants/subvariants and significant mutations which can foster virus evolution, monitoring of disease spread, and outbreak management. Methods: We used Illumina-NovaSeq 6000 for SARS-CoV-2 genome sequencing, MEGA 10.2 and nextstrain tools for phylogeny; CD-HIT program (version 4.8.1) and MUSCLE program for clustering and alignment. At the same time, UCSF Chimera was employed for protein visualization. Results: Predominant Omicron pango lineages in Al-Baha were BA.5.2/B22 (n=4, 57%), and other lineages were BA.2.12/21L (n=1, 14.28%), BV.1/22B (n=1, 14.28%) and BA.5.2.18/22B (n=1, 14.28%). 22B nextstrain clade was predominant, while only one lineage showed 21L. BA.5.2/22B, BA.5.2/22B harbored a maximum of n=24 mutations in the spike region. Twelve crucial RBD mutations: D405N, R408S, K417N, N440K, L452R, S477N, T478K, E484A, F486V, Q498R, N501Y, and Y505H were identified except the lineage BA.5.2/22B in which F486V mutation was not observed. Critical deletions S106 in membrane protein NSP6, E31in nucleocapsid, and L24 in spike region were observed in all the lineages. Furthermore, we identified common mutations of Omicron variants of SARS-CoV-2 in therapeutic hot spot spike region: T19I, D405N, R408S, K417N, N440K, L452R, S477N, T478K, E484A, F486V, Q498R, N501Y, Y505H, D614G, A653V, H655Y, N679K, P681H, N764K, D796Y, Q954H, N969K, D1146D, L452R, F486V, N679K and D796Y. The effect of RBD-targeted mutations on neutralizing (NAbs) binding was considerable. Conclusion: The outcome of this first report on SARS-CoV-2 variants identification and mutation in the Al-Baha region could be used to lay down the policies to manage and impede the regional outbreak of COVID-19 effectively.

Analysis of Thaumatotibia leucotreta (Lepidoptera: Tortricidae: Olethreutinae) mitochondrial genomes in the context of a recent host range expansion.

BACKGROUND: The false codling moth (FCM), Thaumatotibia leucotreta (Meyrick, 1913), is a significant pest of various important economic crops and is a EU quarantine pest. In the last decade the pest has been reported on Rosa spp. In this study we determined whether this shift occurred within specific FCM populations across seven eastern sub-Saharan countries or whether the species opportunistically switches to this novel host as it presents itself. To achieve this, we assessed the genetic diversity of complete mitogenomes of T. leucotreta specimens intercepted at import and analysed potential linkages with the geographical origin and host species. RESULTS: Genomic, geographical and host information were integrated into a T. leucotreta Nextstrain build which contains 95 complete mitogenomes generated from material intercepted at import between January 2013 and December 2018. Samples represented seven sub-Saharan countries and mitogenomic sequences grouped in six main clades. DISCUSSION: If host strains of FCM would exist, specialization from a single haplotype towards the novel host is expected. Instead, we find specimens intercepted on Rosa spp. in all six clades. The absence of linkage between genotype and host suggests opportunistic expansion to the new host plant. This underlines risks of introducing new plant species to an area as the effect of pests already present on the new plant might be unpredictable with current knowledge.

Evaluation of EPISEQ SARS-CoV-2 and a Fully Integrated Application to Identify SARS-CoV-2 Variants from Several Next-Generation Sequencing Approaches.

Whole-genome sequencing has become an essential tool for real-time genomic surveillance of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) worldwide. The handling of raw next-generation sequencing (NGS) data is a major challenge for sequencing laboratories. We developed an easy-to-use web-based application (EPISEQ SARS-CoV-2) to analyse SARS-CoV-2 NGS data generated on common sequencing platforms using a variety of commercially available reagents. This application performs in one click a quality check, a reference-based genome assembly, and the analysis of the generated consensus sequence as to coverage of the reference genome, mutation screening and variant identification according to the up-to-date Nextstrain clade and Pango lineage. In this study, we validated the EPISEQ SARS-CoV-2 pipeline against a reference pipeline and compared the performance of NGS data generated by different sequencing protocols using EPISEQ SARS-CoV-2. We showed a strong agreement in SARS-CoV-2 clade and lineage identification (>99%) and in spike mutation detection (>99%) between EPISEQ SARS-CoV-2 and the reference pipeline. The comparison of several sequencing approaches using EPISEQ SARS-CoV-2 revealed 100% concordance in clade and lineage classification. It also uncovered reagent-related sequencing issues with a potential impact on SARS-CoV-2 mutation reporting. Altogether, EPISEQ SARS-CoV-2 allows an easy, rapid and reliable analysis of raw NGS data to support the sequencing efforts of laboratories with limited bioinformatics capacity and those willing to accelerate genomic surveillance of SARS-CoV-2.

Genetic and Antigenic Characterization of an Expanding H3 Influenza A Virus Clade in U.S. Swine Visualized by Nextstrain.

Defining factors that influence spatial and temporal patterns of influenza A virus (IAV) is essential to inform vaccine strain selection and strategies to reduce the spread of potentially zoonotic swine-origin IAV. The relative frequency of detection of the H3 phylogenetic clade 1990.4.a (colloquially known as C-IVA) in U.S. swine declined to 7% in 2017 but increased to 32% in 2019. We conducted phylogenetic and phenotypic analyses to determine putative mechanisms associated with increased detection. We created an implementation of Nextstrain to visualize the emergence, spatial spread, and genetic evolution of H3 IAV in swine, identifying two C-IVA clades that emerged in 2017 and cocirculated in multiple U.S. states. Phylodynamic analysis of the hemagglutinin (HA) gene documented low relative genetic diversity from 2017 to 2019, suggesting clonal expansion. The major H3 C-IVA clade contained an N156H amino acid substitution, but hemagglutination inhibition (HI) assays demonstrated no significant antigenic drift. The minor HA clade was paired with the neuraminidase (NA) clade N2-2002B prior to 2016 but acquired and maintained an N2-2002A in 2016, resulting in a loss of antigenic cross-reactivity between N2-2002B- and -2002A-containing H3N2 strains. The major C-IVA clade viruses acquired a nucleoprotein (NP) of the H1N1pdm09 lineage through reassortment in the replacement of the North American swine-lineage NP. Instead of genetic or antigenic diversity within the C-IVA HA, our data suggest that population immunity to H3 2010.1 along with the antigenic diversity of the NA and the acquisition of the H1N1pdm09 NP gene likely explain the reemergence and transmission of C-IVA H3N2 in swine. IMPORTANCE Genetically distinct clades of influenza A virus (IAV) in swine undermine efforts to control the disease. Swine producers commonly use vaccines, and vaccine strains are selected by identifying the most common hemagglutinin (HA) gene from viruses detected in a farm or a region. In 2019, we identified an increase in the detection frequency of an H3 phylogenetic clade, C-IVA, which was previously circulating at much lower levels in U.S. swine. Our study identified genetic and antigenic factors contributing to its resurgence by linking comprehensive phylodynamic analyses with empirical wet-lab experiments and visualized these evolutionary analyses in a Nextstrain implementation. The contemporary C-IVA HA genes did not demonstrate an increase in genetic diversity or significant antigenic changes. N2 genes did demonstrate antigenic diversity, and the expanding C-IVA clade acquired a nucleoprotein (NP) gene segment via reassortment. Virus phenotype and vaccination targeting prior dominant HA clades likely contributed to the clade's success.

Whole-Genome Sequence Approach and Phylogenomic Stratification Improve the Association Analysis of Mutations With Patient Data in Influenza Surveillance.

Each year, seasonal influenza results in high mortality and morbidity. The current classification of circulating influenza viruses is mainly focused on the hemagglutinin gene. Whole-genome sequencing (WGS) enables tracking mutations across all influenza segments allowing a better understanding of the epidemiological effects of intra- and inter-seasonal evolutionary dynamics, and exploring potential associations between mutations across the viral genome and patient's clinical data. In this study, mutations were identified in 253 Influenza A (H3N2) clinical isolates from the 2016-2017 influenza season in Belgium. As a proof of concept, available patient data were integrated with this genomic data, resulting in statistically significant associations that could be relevant to improve the vaccine and clinical management of infected patients. Several mutations were significantly associated with the sampling period. A new approach was proposed for exploring mutational effects in highly diverse Influenza A (H3N2) strains through considering the viral genetic background by using phylogenetic classification to stratify the samples. This resulted in several mutations that were significantly associated with patients suffering from renal insufficiency. This study demonstrates the usefulness of using WGS data for tracking mutations across the complete genome and linking these to patient data, and illustrates the importance of accounting for the viral genetic background in association studies. A limitation of this association study, especially when analyzing stratified groups, relates to the number of samples, especially in the context of national surveillance of small countries. Therefore, we investigated if international databases like GISAID may help to verify whether observed associations in the Belgium A (H3N2) samples, could be extrapolated to a global level. This work highlights the need to construct international databases with both information of viral genome sequences and patient data.

Whole-genome-based phylogenomic analysis of the Belgian 2016-2017 influenza A(H3N2) outbreak season allows improved surveillance.

Seasonal influenza epidemics are associated with high mortality and morbidity in the human population. Influenza surveillance is critical for providing information to national influenza programmes and for making vaccine composition predictions. Vaccination prevents viral infections, but rapid influenza evolution results in emerging mutants that differ antigenically from vaccine strains. Current influenza surveillance relies on Sanger sequencing of the haemagglutinin (HA) gene. Its classification according to World Health Organization (WHO) and European Centre for Disease Prevention and Control (ECDC) guidelines is based on combining certain genotypic amino acid mutations and phylogenetic analysis. Next-generation sequencing technologies enable a shift to whole-genome sequencing (WGS) for influenza surveillance, but this requires laboratory workflow adaptations and advanced bioinformatics workflows. In this study, 253 influenza A(H3N2) positive clinical specimens from the 2016-2017 Belgian season underwent WGS using the Illumina MiSeq system. HA-based classification according to WHO/ECDC guidelines did not allow classification of all samples. A new approach, considering the whole genome, was investigated based on using powerful phylogenomic tools including beast and Nextstrain, which substantially improved phylogenetic classification. Moreover, Bayesian inference via beast facilitated reassortment detection by both manual inspection and computational methods, detecting intra-subtype reassortants at an estimated rate of 15 %. Real-time analysis (i.e. as an outbreak is ongoing) via Nextstrain allowed positioning of the Belgian isolates into the globally circulating context. Finally, integration of patient data with phylogenetic groups and reassortment status allowed detection of several associations that would have been missed when solely considering HA, such as hospitalized patients being more likely to be infected with A(H3N2) reassortants, and the possibility to link several phylogenetic groups to disease severity indicators could be relevant for epidemiological monitoring. Our study demonstrates that WGS offers multiple advantages for influenza monitoring in (inter)national influenza surveillance, and proposes an improved methodology. This allows leveraging all information contained in influenza genomes, and allows for more accurate genetic characterization and reassortment detection.

Nextstrain: una herramienta que analiza la epidemiología molecular del SARS-CoV-2 / Nextstrain: a tool to analyze the molecular epidemiology of SARS-CoV-2

La preocupación mundial por el nuevo coronavirus (2019-nCoV), como una amenaza global para la salud pública, fue el motor para que los análisis filogenéticos sufrieran un crecimiento exponencial. El objetivo de esta revisión fue describir el modo de funcionamiento y las bondades de la herramienta Nextstrain, así como el secuenciamiento del virus SARS-CoV-2 en el mundo. Se uso la interfaz de la página de Nextstrain para mostrar sus funcionalidades y los modos de visualización de datos, y se descargaron estos de la web GISAID para mostrar la cantidad de secuenciamientos del SARS-CoV-2 hasta la fecha. Nextstrain es un proyecto de código abierto creado por biólogos bioinformáticos, para aprovechar el potencial científico y de salud pública de los datos de genomas de patógenos. Nextstrain consiste en un conjunto de herramientas que toman secuencias sin procesar (en formato FASTA). Nextstrain realiza una alineación de secuencia de los datos de entrada en alineación de secuencia múltiple basada en la transformación rápida de Fourier. Se basa en el uso de dos softwares: Augur y Auspice. Nextstrain es una herramienta eficaz para mostrar datos epidemiológicos de manera simple para un público no especializado. Puede ser usado en la salud pública, ya que muestra datos en tiempo real de las epidemias y su distribución geográfica. Se puede usar para dar seguimiento a los brotes como es el caso del COVID-19(AU)

Worldwide concern about the novel coronavirus (2019-nCoV) as a global threat to public health is the reason for the exponential growth of phylogenetic analyses. The purpose of this review was to describe the mode of operation and advantages of the tool Nextstrain, as well as the sequencing of the SARS-CoV-2 virus worldwide. The interface of the Nextstrain page was used to show its functions and data visualization modes. These were downloaded from the website GISAID to show the number of SARS-CoV-2 sequencing processes performed so far. Nextstrain is an open code project created by bioinformatics biologists to make good use of the scientific and public health potential of data about genomes of pathogens. Nextstrain consists in a set of tools operating with unprocessed sequences (in FASTA format). Nextstrain performs a sequence alignment of the input data into a multiple sequence alignment based on fast Fourier transform. Its use is based on two software applications: Augur and Auspice. Nextstrain is an efficient tool by which lay people may obtain epidemiological data in a simple manner. It may be used in the public health sector, since it shows real time data about epidemics and their geographic distribution. It may also be used to follow-up outbreaks, as is the case with COVID-19(AU)

Molecular Epidemiology and Evolutionary Trajectory of Emerging Echovirus 30, Europe.

In 2018, an upsurge in echovirus 30 (E30) infections was reported in Europe. We conducted a large-scale epidemiologic and evolutionary study of 1,329 E30 strains collected in 22 countries in Europe during 2016-2018. Most E30 cases affected persons 0-4 years of age (29%) and 25-34 years of age (27%). Sequences were divided into 6 genetic clades (G1-G6). Most (53%) sequences belonged to G1, followed by G6 (23%), G2 (17%), G4 (4%), G3 (0.3%), and G5 (0.2%). Each clade encompassed unique individual recombinant forms; G1 and G4 displayed >2 unique recombinant forms. Rapid turnover of new clades and recombinant forms occurred over time. Clades G1 and G6 dominated in 2018, suggesting the E30 upsurge was caused by emergence of 2 distinct clades circulating in Europe. Investigation into the mechanisms behind the rapid turnover of E30 is crucial for clarifying the epidemiology and evolution of these enterovirus infections.

Genomic Surveillance of Circulating SARS-CoV-2 in South East Italy: A One-Year Retrospective Genetic Study.

In order to provide insights into the evolutionary and epidemiological viral dynamics during the current COVID-19 pandemic in South Eastern Italy, a total of 298 genomes of SARS-CoV-2 strains collected in the Apulia and Basilicata regions, between March 2020 and January 2021, were sequenced. The genomic analysis performed on the draft genomes allowed us to assign the genetic clades and lineages of belonging to each sample and provide an overview of the main circulating viral variants. Our data showed the spread in Apulia and Basilicata of SARS-CoV-2 variants which have emerged during the second wave of infections and are being currently monitored worldwide for their increased transmission rate and their possible impact on vaccines and therapies. These results emphasize the importance of genome sequencing for the epidemiological surveillance of the new SARS-CoV-2 variants' spread.

Geographical Landscape and Transmission Dynamics of SARS-CoV-2 Variants Across India: A Longitudinal Perspective.

Globally, SARS-CoV-2 has moved from one tide to another with ebbs in between. Genomic surveillance has greatly aided the detection and tracking of the virus and the identification of the variants of concern (VOC). The knowledge and understanding from genomic surveillance is important for a populous country like India for public health and healthcare officials for advance planning. An integrative analysis of the publicly available datasets in GISAID from India reveals the differential distribution of clades, lineages, gender, and age over a year (Apr 2020-Mar 2021). The significant insights include the early evidence towards B.1.617 and B.1.1.7 lineages in the specific states of India. Pan-India longitudinal data highlighted that B.1.36* was the predominant clade in India until January-February 2021 after which it has gradually been replaced by the B.1.617.1 lineage, from December 2020 onward. Regional analysis of the spread of SARS-CoV-2 indicated that B.1.617.3 was first seen in India in the month of October in the state of Maharashtra, while the now most prevalent strain B.1.617.2 was first seen in Bihar and subsequently spread to the states of Maharashtra, Gujarat, and West Bengal. To enable a real time understanding of the transmission and evolution of the SARS-CoV-2 genomes, we built a transmission map available on https://covid19-indiana.soic.iupui.edu/India/EmergingLineages/April2020/to/March2021. Based on our analysis, the rate estimate for divergence in our dataset was 9.48 e-4 substitutions per site/year for SARS-CoV-2. This would enable pandemic preparedness with the addition of future sequencing data from India available in the public repositories for tracking and monitoring the VOCs and variants of interest (VOI). This would help aid decision making from the public health perspective.

Infectious Bronchitis Virus Surveillance in Broilers in California (2012-20).

Infectious bronchitis virus (IBV) causes severe economic losses among chicken flocks worldwide. Although IBV molecular surveillance has been conducted in California broilers, seasonal and spatial-temporal trends in IBV prevalence are poorly defined. The goals of this study were to evaluate seasonal and spatial-temporal trends in IBV prevalence and to determine the predominant IBV genotypes obtained over the last 8 yr from a broiler company located in the California Central Valley. In total, 3439 broilers with a suspicion of IBV infection were submitted to the California Animal Health and Food Safety laboratories between January 2012 and February 2020. Swabs from tracheas, kidneys, and cecal tonsils from each submission were independently pooled and screened for IBV using reverse transcriptase quantitative PCR (RT-qPCR). Positive samples were submitted for virus isolation. Viral isolates were subject to a conventional RT-PCR targeting the S1 gene hypervariable region. Positive samples from this RT-PCR were sequenced, and phylogenetic analyses were performed. In total, 1243 pooled swab samples were positive for IBV. Positive results were more frequently detected in fall and winter months compared to spring. Spatial analyses revealed an IBV hot spot in the vicinity of Livingston, and two areas with a low prevalence (i.e., cold spots) around Riverdale. The IBV spatial-temporal distribution identified three significant clusters: one hot spot around Turlock from 2015 to 2016, a second hot spot around Merced from 2012 to 2016, and a cold spot around Fresno from 2017 to 2020. Predominant genotypes changed over time from IBV Cal 99, which was predominant between 2012 and 2014, to IBV 3099 in 2019. Vaccination efforts were initiated in 2018, and as a result, we detected an emerging variant with 92% similarity to CA 3099 in 2020. This work highlights the importance of ongoing surveillance in IBV prevention programs. Surveillance strategies are necessary to monitor trends in diseases such as infectious bronchitis, and the tools used for surveillance need to be sensitive enough to detect new variants and identify spatial-temporal trends.

Vigilancia del virus de la bronquitis infecciosa en pollos de engorde en California (2012­2020). El virus de la bronquitis infecciosa (IBV) causa graves pérdidas económicas entre las parvadas de pollos en todo el mundo. Aunque la vigilancia molecular del virus de la bronquitis infecciosa se ha realizado en pollos de engorde en el estado de California, las tendencias estacionales y espacio-temporales sobre la prevalencia de este virus están mal definidas. Los objetivos de este estudio fueron evaluar las tendencias estacionales y espacio-temporales sobre la prevalencia del virus de la bronquitis infecciosa y determinar los genotipos predominantes de este virus obtenidos durante los últimos ocho años de una empresa de pollos de engorde ubicada en el Valle Central de California. En total, 3439 pollos de engorde con sospecha de infección por el virus de la bronquitis infecciosa se enviaron a los laboratorios de Salud Animal y Seguridad Alimentaria del estado de California entre enero del 2012 y febrero del 2020. Los hisopos de tráqueas, riñones y tonsilas cecales de cada caso se combinaron de forma independiente y se examinaron para detectar al virus de la bronquitis utilizando transcripción reversa y un método cuantitativo de PCR (RT-qPCR). Se enviaron muestras positivas para aislamiento del virus. Los aislados virales se sometieron a un método convencional de RT-PCR dirigido a la región hipervariable del gene S1. Se secuenciaron muestras positivas mediante la prueba RT-PCR y se realizaron análisis filogenéticos. Un total de 1243 muestras combinadas de hisopos dieron positivo para el virus de la bronquitis infecciosa. Los resultados positivos se detectaron con mayor frecuencia en los meses de otoño e invierno en comparación con la primavera. Los análisis espaciales revelaron un punto activo para el virus de la bronquitis infecciosa en las cercanías de Livingston y dos áreas con una baja prevalencia (es decir, puntos fríos) alrededor de Riverdale. La distribución espacio-temporal del virus de la bronquitis identificó tres grupos importantes: un punto activo alrededor de Turlock entre los años 2015 a 2016, un segundo punto activo alrededor de Merced entre los años 2012 a 2016 y un punto frío alrededor de Fresno entre los años 2017 a 2020. Los genotipos predominantes cambiaron con el tiempo, consideraron el subtipo IBV Cal 99, que fue predominante entre 2012 y 2014, a el tipo IBV 3099 en 2019. Los esfuerzos de vacunación se iniciaron en el 2018, y como resultado, detectamos una variante emergente con un 92% de similitud con el virus CA 3099 en 2020. Este trabajo destaca la importancia de la vigilancia continua en los programas de prevención para la bronquitis infecciosa. Las estrategias de vigilancia son necesarias para monitorear las tendencias en enfermedades como la bronquitis infecciosa, y las herramientas utilizadas para la vigilancia deben ser lo suficientemente sensibles como para detectar nuevas variantes e identificar tendencias espacio-temporales.