Template switching and duplications in SARS-CoV-2 genomes give rise to insertion variants that merit monitoring

The appearance of multiple new SARS-CoV-2 variants during the COVID-19 pandemic is a matter of grave concern. Some of these variants, such as B.1.617.2, B.1.1.7, and B.1.351, manifest higher infectivity and virulence than the earlier SARS-CoV-2 variants, with potential dramatic effects on the course of the pandemic. So far, analysis of new SARS-CoV-2 variants focused primarily on nucleotide substitutions and short deletions that are readily identifiable by comparison to consensus genome sequences. In contrast, insertions have largely escaped the attention of researchers although the furin site insert in the Spike (S) protein is thought to be a determinant of SARS-CoV-2 virulence. Here, we identify 346 unique inserts of different lengths in SARS-CoV-2 genomes and present evidence that these inserts reflect actual virus variance rather than sequencing artifacts. Two principal mechanisms appear to account for the inserts in the SARS-CoV-2 genomes, polymerase slippage and template switch that might be associated with the synthesis of subgenomic RNAs. At least three inserts in the N-terminal domain of the S protein are predicted to lead to escape from neutralizing antibodies, whereas other inserts might result in escape from T-cell immunity. Thus, inserts in the S protein can affect its antigenic properties and merit monitoring.

Insertions in SARS-CoV-2 genome caused by template switch and duplications give rise to new variants of potential concern

The appearance of multiple new SARS-CoV-2 variants during the winter of 2020-2021 is a matter of grave concern. Some of these new variants, such as B.1.351 and B.1.1.17, manifest higher infectivity and virulence than the earlier SARS-CoV-2 variants, with potential dramatic effects on the course of the COVID-19 pandemic. So far, analysis of new SARS-CoV-2 variants focused primarily on point nucleotide substitutions and short deletions that are readily identifiable by comparison to consensus genome sequences. In contrast, insertions have largely escaped the attention of researchers although the furin site insert in the spike protein is thought to be a determinant of SARS-CoV-2 virulence and other inserts might have contributed to coronavirus pathogenicity as well. Here, we investigate insertions in SARS-CoV-2 genomes and identify 141 unique inserts of different lengths. We present evidence that these inserts reflect actual virus variance rather than sequencing errors. Two principal mechanisms appear to account for the inserts in the SARS-CoV-2 genomes, polymerase slippage and template switch that might be associated with the synthesis of subgenomic RNAs. We show that inserts in the Spike glycoprotein can affect its antigenic properties and thus have to be monitored. At least, two inserts in the N-terminal domain of the Spike (ins246DSWG and ins15ATLRI) that were first detected in January 2021 are predicted to lead to escape from neutralizing antibodies whereas other inserts might result in escape from T-cell immunity.

DNA polymerase stalling at structured DNA constrains the expansion of Short Tandem Repeats

BackgroundShort tandem repeats (STRs) contribute significantly to de novo mutagenesis, driving phenotypic diversity and genetic disease. Although highly diverse, their repetitive sequences induce DNA polymerase slippage and stalling, leading to length and sequence variation. However, current studies of DNA synthesis through STRs are restricted to a handful of selected sequences, limiting our broader understanding of their evolutionary behaviour and hampering the characterisation of the determinants of their abundance and stability in eukaryotic genomes.ResultsWe perform a comprehensive analysis of DNA synthesis at all STR permutations and interrogate the impact of STR sequence and secondary structure on their genomic representation and mutability. To do so, we developed a high-throughput primer extension assay that allows monitoring of the kinetics and fidelity of DNA synthesis through 20,000 sequences comprising all STR permutations in different lengths. By combining these measurements with population-scale genomic data, we show that the response of a model replicative DNA polymerase to variously structured DNA is sufficient to predict the complex genomic behaviour of STRs, including abundance and mutational constraints. We demonstrate that DNA polymerase stalling at DNA structures induces error-prone DNA synthesis, which constrains STR expansion.ConclusionsOur data support a model in which STR length in eukaryotic genomes results from a balance between expansion due to polymerase slippage at repeated DNA sequences and point mutations caused by error-prone DNA synthesis at DNA structures.

Landscape of multi-nucleotide variants in 125,748 human exomes and 15,708 genomes

Multi-nucleotide variants (MNVs), defined as two or more nearby variants existing on the same haplotype in an individual, are a clinically and biologically important class of genetic variation. However, existing tools for variant interpretation typically do not accurately classify MNVs, and understanding of their mutational origins remains limited. Here, we systematically survey MNVs in 125,748 whole exomes and 15,708 whole genomes from the Genome Aggregation Database (gnomAD). We identify 1,996,125 MNVs across the genome with constituent variants falling within 2 bp distance of one another, of which 31,510 exist within the same codon, including 405 predicted to result in gain of a nonsense mutation, 1,818 predicted to rescue a nonsense mutation event that would otherwise be caused by one of the constituent variants, and 16,481 additional variants predicted to alter protein sequences. We show that the distribution of MNVs is highly non-uniform across the genome, and that this non-uniformity can be largely explained by a variety of known mutational mechanisms, such as CpG deamination, replication error by polymerase zeta, or polymerase slippage at repeat junctions. We also provide an estimate of the dinucleotide mutation rate caused by polymerase zeta. Finally, we show that differential CpG methylation drives MNV differences across functional categories. Our results demonstrate the importance of incorporating haplotype-aware annotation for accurate functional interpretation of genetic variation, and refine our understanding of genome-wide mutational mechanisms of MNVs.

The P1N-PISPOtrans-Frame Gene of Sweet Potato Feathery Mottle Potyvirus Is Produced during Virus Infection and Functions as an RNA Silencing Suppressor

ABSTRACTThe positive-sense RNA genome ofSweet potato feathery mottle virus(SPFMV) (genusPotyvirus, familyPotyviridae) contains a large open reading frame (ORF) of 3,494 codons translatable as a polyprotein and two embedded shorter ORFs in the −1 frame: PISPO, of 230 codons, and PIPO, of 66 codons, located in the P1 and P3 regions, respectively. PISPO is specific to some sweet potato-infecting potyviruses, while PIPO is present in all potyvirids. In SPFMV these two extra ORFs are preceded by conserved G2A6motifs. We have shown recently that a polymerase slippage mechanism at these sites could produce transcripts bringing these ORFs in frame with the upstream polyprotein, thus leading to P1N-PISPO and P3N-PIPO products (B. Rodamilans, A. Valli, A. Mingot, D. San Leon, D. B. Baulcombe, J. J. Lopez-Moya, and J.A. Garcia, J Virol 89:6965–6967, 2015, doi:10.1128/JVI.00337-15). Here, we demonstrate by liquid chromatography coupled to mass spectrometry that both P1 and P1N-PISPO are produced during viral infection and coexist in SPFMV-infectedIpomoea batatasplants. Interestingly, transient expression of SPFMV gene products coagroinfiltrated with a reporter gene inNicotiana benthamianarevealed that P1N-PISPO acts as an RNA silencing suppressor, a role normally associated with HCPro in other potyviruses. Moreover, mutation of WG/GW motifs present in P1N-PISPO abolished its silencing suppression activity, suggesting that the function might require interaction with Argonaute components of the silencing machinery, as was shown for other viral suppressors. Altogether, our results reveal a further layer of complexity of the RNA silencing suppression activity within thePotyviridaefamily.IMPORTANCEGene products of potyviruses include P1, HCPro, P3, 6K1, CI, 6K2, VPg/NIaPro, NIb, and CP, all derived from the proteolytic processing of a large polyprotein, and an additional P3N-PIPO product, with the PIPO segment encoded in a different frame within the P3 cistron. In sweet potato feathery mottle virus (SPFMV), another out-of-frame element (PISPO) was predicted within the P1 region. We have shown recently that a polymerase slippage mechanism can generate the transcript variants with extra nucleotides that could be translated into P1N-PISPO and P3N-PIPO. Now, we demonstrate by mass spectrometry analysis that P1N-PISPO is indeed produced in SPFMV-infected plants, in addition to P1. Interestingly, while in other potyviruses the suppressor of RNA silencing is HCPro, we show here that P1N-PISPO exhibited this activity in SPFMV, revealing how the complexity of the gene content could contribute to supply this essential function in members of thePotyviridaefamily.

Length, orientation, and plant host influence the mutation frequency in microsatellites

Microsatellites are simple, tandem DNA repeats that represent unstable regions of the genome. They undergo frequent changes in tract length by base additions or deletions due to DNA polymerase slippage during replication. To characterize factors affecting the frequency of spontaneous mutations occurring in microsatellites in plants, a reporter system was used in Arabidopsis thaliana and tomato ( Lycopersicon esculentum ). The β-glucuronidase (GUS) reporter system was used to measure the mutation frequency in various microsatellites (G7, G10, G13, G16, and C16) in somatic tissues. Our results indicate that this frequency increases with the number of repeats: a G16 tract was almost 80-fold more mutable than a G7 tract. Furthermore, the frequency of mutations depends on repeat orientation, as G16 was 3-fold more mutable than C16. The mutation rate was also found to differ markedly in Arabidopsis and tomato for an identical microsatellite. Indeed, Arabidopsis showed a 5-fold higher mutation frequency than tomato with the same G7 reporter construct. Finally, mutation in a G16 tract was frequent enough that mutations transmitted germinally to the next generation could be detected at a relatively high frequency.