Insights into SARS-CoV-2 in Angola during the COVID-19 peak: Molecular epidemiology and genome surveillance.

Background: In Angola, COVID-19 cases have been reported in all provinces, resulting in >105,000 cases and >1900 deaths. However, no detailed genomic surveillance into the introduction and spread of the SARS-CoV-2 virus has been conducted in Angola. We aimed to investigate the emergence and epidemic progression during the peak of the COVID-19 pandemic in Angola. Methods: We generated 1210 whole-genome SARS-CoV-2 sequences, contributing West African data to the global context, that were phylogenetically compared against global strains. Virus movement events were inferred using ancestral state reconstruction. Results: The epidemic in Angola was marked by four distinct waves of infection, dominated by 12 virus lineages, including VOCs, VOIs, and the VUM C.16, which was unique to South-Western Africa and circulated for an extended period within the region. Virus exchanges occurred between Angola and its neighboring countries, and strong links with Brazil and Portugal reflected the historical and cultural ties shared between these countries. The first case likely originated from southern Africa. Conclusion: A lack of a robust genome surveillance network and strong dependence on out-of-country sequencing limit real-time data generation to achieve timely disease outbreak responses, which remains of the utmost importance to mitigate future disease outbreaks in Angola.

SARS-CoV-2 Genetic Diversity and Lineage Dynamics in Egypt during the First 18 Months of the Pandemic.

COVID-19 was first diagnosed in Egypt on 14 February 2020. By the end of November 2021, over 333,840 cases and 18,832 deaths had been reported. As part of the national genomic surveillance, 1027 SARS-CoV-2 near whole-genomes were generated and published by the end of July 2021. Here we describe the genomic epidemiology of SARS-CoV-2 in Egypt over this period using a subset of 976 high-quality Egyptian genomes analyzed together with a representative set of global sequences within a phylogenetic framework. A single lineage, C.36, introduced early in the pandemic was responsible for most of the cases in Egypt. Furthermore, to remain dominant in the face of mounting immunity from previous infections and vaccinations, this lineage acquired several mutations known to confer an adaptive advantage. These results highlight the value of continuous genomic surveillance in regions where VOCs are not predominant and the need for enforcement of public health measures to prevent expansion of the existing lineages.

Targeted Sanger sequencing to recover key mutations in SARS-CoV-2 variant genome assemblies produced by next-generation sequencing.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is adaptively evolving to ensure its persistence within human hosts. It is therefore necessary to continuously monitor the emergence and prevalence of novel variants that arise. Importantly, some mutations have been associated with both molecular diagnostic failures and reduced or abrogated next-generation sequencing (NGS) read coverage in some genomic regions. Such impacts are particularly problematic when they occur in genomic regions such as those that encode the spike (S) protein, which are crucial for identifying and tracking the prevalence and dissemination dynamics of concerning viral variants. Targeted Sanger sequencing presents a fast and cost-effective means to accurately extend the coverage of whole-genome sequences. We designed a custom set of primers to amplify a 401 bp segment of the receptor-binding domain (RBD) (between positions 22698 and 23098 relative to the Wuhan-Hu-1 reference). We then designed a Sanger sequencing wet-laboratory protocol. We applied the primer set and wet-laboratory protocol to sequence 222 samples that were missing positions with key mutations K417N, E484K, and N501Y due to poor coverage after NGS sequencing. Finally, we developed SeqPatcher, a Python-based computational tool to analyse the trace files yielded by Sanger sequencing to generate consensus sequences, or take preanalysed consensus sequences in fasta format, and merge them with their corresponding whole-genome assemblies. We successfully sequenced 153 samples of 222 (69 %) using Sanger sequencing and confirmed the occurrence of key beta variant mutations (K417N, E484K, N501Y) in the S genes of 142 of 153 (93 %) samples. Additionally, one sample had the Y508F mutation and four samples the S477N. Samples with RT-PCR Ct scores ranging from 13.85 to 37.47 (mean=25.70) could be Sanger sequenced efficiently. These results show that our method and pipeline can be used to improve the quality of whole-genome assemblies produced using NGS and can be used with any pairs of the most used NGS and Sanger sequencing platforms.

Viral genome sequence datasets display pervasive evidence of strand-specific substitution biases that are best described using non-reversible nucleotide substitution models.

Background: The vast majority of phylogenetic trees are inferred from molecular sequence data (nucleotides or amino acids) using time-reversible evolutionary models which assume that, for any pair of nucleotide or amino acid characters, the relative rate of X to Y substitution is the same as the relative rate of Y to X substitution. However, this reversibility assumption is unlikely to accurately reflect the actual underlying biochemical and/or evolutionary processes that lead to the fixation of substitutions. Here, we use empirical viral genome sequence data to reveal that evolutionary non-reversibility is pervasive among most groups of viruses. Specifically, we consider two non-reversible nucleotide substitution models: (1) a 6-rate non-reversible model (NREV6) in which Watson-Crick complementary substitutions occur at identical relative rates and which might therefor be most applicable to analyzing the evolution of genomes where both complementary strands are subject to the same mutational processes (such as might be expected for double-stranded (ds) RNA or dsDNA genomes); and (2) a 12-rate non-reversible model (NREV12) in which all relative substitution types are free to occur at different rates and which might therefore be applicable to analyzing the evolution of genomes where the complementary genome strands are subject to different mutational processes (such as might be expected for viruses with single-stranded (ss) RNA or ssDNA genomes). Results: Using likelihood ratio and Akaike Information Criterion-based model tests, we show that, surprisingly, NREV12 provided a significantly better fit to 21/31 dsRNA and 20/30 dsDNA datasets than did the general time reversible (GTR) and NREV6 models with NREV6 providing a better fit than NREV12 and GTR in only 5/30 dsDNA and 2/31 dsRNA datasets. As expected, NREV12 provided a significantly better fit to 24/33 ssDNA and 40/47 ssRNA datasets. Next, we used simulations to show that increasing degrees of strand-specific substitution bias decrease the accuracy of phylogenetic inference irrespective of whether GTR or NREV12 is used to describe mutational processes. However, in cases where strand-specific substitution biases are extreme (such as in SARS-CoV-2 and Torque teno sus virus datasets) NREV12 tends to yield more accurate phylogenetic trees than those obtained using GTR. Conclusion: We show that NREV12 should, be seriously considered during the model selection phase of phylogenetic analyses involving viral genomic sequences.

Integrated molecular characterisation of the MAPK pathways in human cancers reveals pharmacologically vulnerable mutations and gene dependencies.

The mitogen-activated protein kinase (MAPK) pathways are crucial regulators of the cellular processes that fuel the malignant transformation of normal cells. The molecular aberrations which lead to cancer involve mutations in, and transcription variations of, various MAPK pathway genes. Here, we examine the genome sequences of 40,848 patient-derived tumours representing 101 distinct human cancers to identify cancer-associated mutations in MAPK signalling pathway genes. We show that patients with tumours that have mutations within genes of the ERK-1/2 pathway, the p38 pathways, or multiple MAPK pathway modules, tend to have worse disease outcomes than patients with tumours that have no mutations within the MAPK pathways genes. Furthermore, by integrating information extracted from various large-scale molecular datasets, we expose the relationship between the fitness of cancer cells after CRISPR mediated gene knockout of MAPK pathway genes, and their dose-responses to MAPK pathway inhibitors. Besides providing new insights into MAPK pathways, we unearth vulnerabilities in specific pathway genes that are reflected in the re sponses of cancer cells to MAPK targeting drugs: a revelation with great potential for guiding the development of innovative therapies.

Integrative landscape of dysregulated signaling pathways of clinically distinct pancreatic cancer subtypes.

Despite modern therapeutic advances, the survival prospects of pancreatic cancer patients have remained poor. Besides being highly metastatic, pancreatic cancer is challenging to treat because it is caused by a heterogeneous array of somatic mutations that impact a variety of signaling pathways and cellular regulatory systems. Here we use publicly available transcriptomic, copy number alteration and mutation profiling datasets from pancreatic cancer patients together with data on disease outcomes to show that the three major pancreatic cancer subtypes each display distinctive aberrations in cell signaling and metabolic pathways. Importantly, patients afflicted with these different pancreatic cancer subtypes also exhibit distinctive survival profiles. Within these patients, we find that dysregulation of the phosphoinositide 3-kinase and mitogen-activated protein kinase pathways, and p53 mediated disruptions of cell cycle processes are apparently drivers of disease. Further, we identify for the first time the molecular perturbations of signalling networks that are likely the primary causes of oncogenesis in each of the three pancreatic cancer subtypes.

Zika virus complete genome from Salvador, Bahia, Brazil.

In May 2015 the first autochthonous Zika virus infection was reported in Brazil. Rapid and urgent measures are needed to contain the ongoing outbreak. Here we report the full-length ZIKV coding sequence from Bahia. Genetic analysis of outbreak sequences will be essential for characterizing the diversity of circulating strains, identifying hotspots of virus transmission and guiding public health control. Rapid and urgent measures are needed to contain the ongoing outbreak.

Phylogenetic Exploration of Nosocomial Transmission Chains of 2009 Influenza A/H1N1 among Children Admitted at Red Cross War Memorial Children's Hospital, Cape Town, South Africa in 2011.

Traditional modes of investigating influenza nosocomial transmission have entailed a combination of confirmatory molecular diagnostic testing and epidemiological investigation. Common hospital-acquired infections like influenza require a discerning ability to distinguish between viral isolates to accurately identify patient transmission chains. We assessed whether influenza hemagglutinin sequence phylogenies can be used to enrich epidemiological data when investigating the extent of nosocomial transmission over a four-month period within a paediatric Hospital in Cape Town South Africa. Possible transmission chains/channels were initially determined through basic patient admission data combined with Maximum likelihood and time-scaled Bayesian phylogenetic analyses. These analyses suggested that most instances of potential hospital-acquired infections resulted from multiple introductions of Influenza A into the hospital, which included instances where virus hemagglutinin sequences were identical between different patients. Furthermore, a general inability to establish epidemiological transmission linkage of patients/viral isolates implied that identified isolates could have originated from asymptomatic hospital patients, visitors or hospital staff. In contrast, a traditional epidemiological investigation that used no viral phylogenetic analyses, based on patient co-admission into specific wards during a particular time-frame, suggested that multiple hospital acquired infection instances may have stemmed from a limited number of identifiable index viral isolates/patients. This traditional epidemiological analysis by itself could incorrectly suggest linkage between unrelated cases, underestimate the number of unique infections and may overlook the possible diffuse nature of hospital transmission, which was suggested by sequencing data to be caused by multiple unique introductions of influenza A isolates into individual hospital wards. We have demonstrated a functional role for viral sequence data in nosocomial transmission investigation through its ability to enrich traditional, non-molecular observational epidemiological investigation by teasing out possible transmission pathways and working toward more accurately enumerating the number of possible transmission events.

SDT: a virus classification tool based on pairwise sequence alignment and identity calculation.

The perpetually increasing rate at which viral full-genome sequences are being determined is creating a pressing demand for computational tools that will aid the objective classification of these genome sequences. Taxonomic classification approaches that are based on pairwise genetic identity measures are potentially highly automatable and are progressively gaining favour with the International Committee on Taxonomy of Viruses (ICTV). There are, however, various issues with the calculation of such measures that could potentially undermine the accuracy and consistency with which they can be applied to virus classification. Firstly, pairwise sequence identities computed based on multiple sequence alignments rather than on multiple independent pairwise alignments can lead to the deflation of identity scores with increasing dataset sizes. Also, when gap-characters need to be introduced during sequence alignments to account for insertions and deletions, methodological variations in the way that these characters are introduced and handled during pairwise genetic identity calculations can cause high degrees of inconsistency in the way that different methods classify the same sets of sequences. Here we present Sequence Demarcation Tool (SDT), a free user-friendly computer program that aims to provide a robust and highly reproducible means of objectively using pairwise genetic identity calculations to classify any set of nucleotide or amino acid sequences. SDT can produce publication quality pairwise identity plots and colour-coded distance matrices to further aid the classification of sequences according to ICTV approved taxonomic demarcation criteria. Besides a graphical interface version of the program for Windows computers, command-line versions of the program are available for a variety of different operating systems (including a parallel version for cluster computing platforms).

Evidence of pervasive biologically functional secondary structures within the genomes of eukaryotic single-stranded DNA viruses.

Single-stranded DNA (ssDNA) viruses have genomes that are potentially capable of forming complex secondary structures through Watson-Crick base pairing between their constituent nucleotides. A few of the structural elements formed by such base pairings are, in fact, known to have important functions during the replication of many ssDNA viruses. Unknown, however, are (i) whether numerous additional ssDNA virus genomic structural elements predicted to exist by computational DNA folding methods actually exist and (ii) whether those structures that do exist have any biological relevance. We therefore computationally inferred lists of the most evolutionarily conserved structures within a diverse selection of animal- and plant-infecting ssDNA viruses drawn from the families Circoviridae, Anelloviridae, Parvoviridae, Nanoviridae, and Geminiviridae and analyzed these for evidence of natural selection favoring the maintenance of these structures. While we find evidence that is consistent with purifying selection being stronger at nucleotide sites that are predicted to be base paired than at sites predicted to be unpaired, we also find strong associations between sites that are predicted to pair with one another and site pairs that are apparently coevolving in a complementary fashion. Collectively, these results indicate that natural selection actively preserves much of the pervasive secondary structure that is evident within eukaryote-infecting ssDNA virus genomes and, therefore, that much of this structure is biologically functional. Lastly, we provide examples of various highly conserved but completely uncharacterized structural elements that likely have important functions within some of the ssDNA virus genomes analyzed here.

Detection of recombination breakpoints in the genomes of human enterovirus 71 strains isolated in the Netherlands in epidemic and non-epidemic years, 1963-2010.

Evolutionary history of human enterovirus 71 (EV71) in the Netherlands shows displacement of virus subgenogroups, that only partly can be explained by antigenic changes. Additionally, occasional epidemics have occurred that remain to be explained. Previous studies have shown subgenogroup specific recombination events in the genome of Asian EV71 strains. To find clues on the role of genome recombination in evolution of the EV71 subgenogroups found in Europe and in the evolution of strains capable of causing outbreaks, we analyzed the genomes of 19 strains representing the genetic diversity of EV71 in the Netherlands between 1963 and 2010. We selected viruses from EV71 endemic and epidemic years (1986 and 2007). Subgenogroup specific genome recombination events were detected for subgenogroup B0, B1 and B2 viruses, in line with observed genome recombination events in Asian subgenogroup B3 and B4 viruses. Considering recombination events distinguishing strains from epidemic years from those of non-epidemic years, breakpoints for recombination were detected in the 5'UTR of B2 viruses from the outbreak in 1986, with highest similarity of the 5'UTR to B4 and B3 strains isolated during outbreaks in the Asian Pacific region. No indications for recombination were found in genogroup C isolates. Except for the '86 B2 isolates' Dutch isolates phylogenetically interspersed with international reference strains of the same subgenogroup, indicating a global dissemination of (recombinant) EV71 viruses. The difference observed in the 5'UTR of EV71 strains isolated in endemic versus epidemic years suggests that changes in the 5'UTR cause evolution of strains capable of causing outbreaks.

Genome analysis of a severe and a mild isolate of Papaya ringspot virus-type W found in Brazil.

Papaya ringspot virus-type W (PRSV-W) is one of the most economically threatening viruses of cucurbits in Brazil. Premunization is one of the most effective PRSV control measures currently applied in squash and zucchini crops. PRSV-W-1, a mild and premunizing strain of PRSV has been successfully used to protect cucurbits against both the severe PRSV-W-C strain and other Brazilian PRSVs. To aid in understanding the mechanism by which PRSV-W-1 premunization operates, the complete genome sequences of PRSV-W-1 and PRSV-W-C were determined. PRSV-W-1 had a genome size of 10,332 nucleotides, whereas indels within the coat protein encoding gene meant that the genome size of PRSV-W-C was six nucleotides shorter than that of the mild strain. The genomes of the two strains shared 94.63% nucleotide sequence identity, with the 5' UTR and P1 being the most variable regions, and the coat protein and 3' UTR being the most conserved. Rigorous recombination analysis revealed that neither PRSV-W-1 nor PRSV-W-C was obviously recombinant, there was significant evidence that many other fully sequenced PRSV genomes were recombinant.