Rapid selection of P323L in the SARS-CoV-2 polymerase (NSP12) in humans and non-human primate models and confers a large plaque phenotype

The mutational landscape of SARS-CoV-2 varies at both the dominant viral genome sequence and minor genomic variant population. An early change associated with transmissibility was the D614G substitution in the spike protein. This appeared to be accompanied by a P323L substitution in the viral polymerase (NSP12), but this latter change was not under strong selective pressure. Investigation of P323L/D614G changes in the human population showed rapid emergence during the containment phase and early surge phase of wave 1 in the UK. This rapid substitution was from minor genomic variants to become part of the dominant viral genome sequence. A rapid emergence of 323L but not 614G was observed in a non-human primate model of COVID-19 using a starting virus with P323 and D614 in the dominant genome sequence and 323L and 614G in the minor variant population. In cell culture, a recombinant virus with 323L in NSP12 had a larger plaque size than the same recombinant virus with P323. These data suggest that it may be possible to predict the emergence of a new variant based on tracking the distribution and frequency of minor variant genomes at a population level, rather than just focusing on providing information on the dominant viral genome sequence e.g., consensus level reporting. The ability to predict an emerging variant of SARS-CoV-2 in the global landscape may aid in the evaluation of medical countermeasures and non-pharmaceutical interventions.

INSaFLU-TELE-Vir: an open web-based bioinformatics suite for influenza and SARS-CoV-2 genome-based surveillance

A new era of virus surveillance is emerging based on the real-time monitoring of virus evolution at whole-genome scale (World Health Organization 2021). Although national and international health authorities have strongly recommended this technological transition, especially for influenza and SARS-CoV-2 (World Health Organization 2021, Revez et al. 2017), the implementation of genomic surveillance can be particularly challenging due to the lack of bioinformatics infrastructures and/or expertise to process and interpret next-generation sequencing (NGS) data (Oakeson et al. 2017). We developed and implemented INSaFLU-TELE-Vir platform (https://insaflu.insa.pt/) (Borges et al. 2018), which is an influenza- and SARS-CoV-2-oriented bioinformatics free web-based suite that handles primary NGS data (reads) towards the automatic generation of the main “genetic requests'' for effective and timely laboratory surveillance. By handling NGS data collected from any amplicon-based schema (making it applicable for other pathogens), INSaFLU-TELE-Vir enables any laboratory to perform multi-step and intensive bioinformatics analyses in a user-oriented manner without requiring advanced training. INSaFLU-TELE-Vir handles NGS data collected from distinct sequencing technologies (Illumina, Ion Torrent and Oxford Nanopore Technologies), with the possibility of constructing comparative analyses using different technologies. It gives access to user-restricted sample databases and project management, being a transparent and flexible tool specifically designed to automatically update project outputs as more samples are uploaded. Data integration is thus cumulative and scalable, fitting the need for both routine surveillance and outbreak investigation activities. The bioinformatics pipeline consists of six core steps: read quality analysis and improvement, human betacoronaviruses (including SARS-CoV-2 Pango lineages) and influenza type/subtype classification, mutation detection and consensus generation, coverage analysis, alignment/phylogeny, intra-host minor variant detection (and automatic detection of putative mixed infections). read quality analysis and improvement, human betacoronaviruses (including SARS-CoV-2 Pango lineages) and influenza type/subtype classification, mutation detection and consensus generation, coverage analysis, alignment/phylogeny, intra-host minor variant detection (and automatic detection of putative mixed infections). The multiple outputs are provided in nomenclature-stable and standardized formats that can be visualized and explored in situ or through multiple compatible downstream applications for fine-tuned data analysis. Novel features are being implemented into the INSaFLU-TELE-Vir bioinformatics toolkit as part of the OHEJP TELE-Vir (https://onehealthejp.eu/jrp-tele-vir/) project, including rapid detection of selected genotype-phenotype associations, and enhanced geotemporal data visualization. All the code is available in github (https://github.com/INSaFLU) with the possibility of a local docker installation (https://github.com/INSaFLU/docker). A detailed documentation and tutorial is also available (https://insaflu.readthedocs.io/en/latest/). In summary, INSaFLU supplies public health laboratories and researchers with an open and user-friendly framework, potentiating a strengthened and timely multi-country genome-based virus surveillance.

Sequence analysis of SARS-CoV-2 in nasopharyngeal samples from patients with COVID-19 illustrates population variation and diverse phenotypes, placing the in vitro growth properties of B.1.1.7 and B.1.351 lineage viruses in context.

New variants of SARS-CoV-2 are continuing to emerge and dominate the regional and global sequence landscapes. Several variants have been labelled as Variants of Concern (VOCs) because of perceptions or evidence that these may have a transmission advantage, increased risk of morbidly and/or mortality or immune evasion in the context of prior infection or vaccination. Placing the VOCs in context and also the underlying variability of SARS-CoV-2 is essential in understanding virus evolution and selection pressures. Sequences of SARS-CoV-2 in nasopharyngeal swabs from hospitalised patients in the UK were determined and virus isolated. The data indicated the virus existed as a population with a consensus level and non-synonymous changes at a minor variant. For example, viruses containing the nsp12 P323L variation from the Wuhan reference sequence, contained minor variants at the position including P and F and other amino acids. These populations were generally preserved when isolates were amplified in cell culture. In order to place VOCs B.1.1.7 (the UK Kent variant) and B.1.351 (the South African variant) in context their growth was compared to a spread of other clinical isolates. The data indicated that the growth in cell culture of the B.1.1.7 VOC was no different from other variants, suggesting that its apparent transmission advantage was not down to replicating more quickly. Growth of B.1.351 was towards the higher end of the variants. Overall, the study suggested that studying the biology of SARS-CoV-2 is complicated by population dynamics and that these need to be considered with new variants.

Reducing the Number of Mismatches between Hairs and Buccal References When Analysing mtDNA Heteroplasmic Variation by Massively Parallel Sequencing

In forensics, mitochondrial DNA (mtDNA) analysis is foremost applied to rootless hairs often lacking detectable nuclear DNA. Sanger sequencing is the routine mtDNA method in most forensic laboratories, even though interpretation of mixed samples and heteroplasmic sites can be challenging. Individuals may hold cells with low-level heteroplasmy variants below the detection threshold and other cells where this minor variant is the major one. This difference may be interpreted as a mismatch between reference and evidentiary trace samples, such as buccal specimens and rootless hairs. Such mismatches may be solved by Massively Parallel Sequencing (MPS), allowing more sensitive quantitative analysis for mixed positions than Sanger. The mtDNA control region was analysed in buccal reference samples from 26 individuals and 475 corresponding hairs by MPS and compared to Sanger sequencing data generated on the same samples. With MPS, mixed contributions down to 3% were regarded, leading to a substantial increase in the frequency of heteroplasmy. Our results demonstrate that previously reported mismatches between buccal reference and hair shaft samples by Sanger are detected as low-level heteroplasmy by MPS. A detailed overview of buccal and hair heteroplasmy is provided and implications for MPS-based mtDNA analysis in the context of forensic cases are discussed.

Using the Minor Variant Finder software to identify and quantify the allelic burden level of somatic mutations in oncohematologic diseases

Background. There are problems related to both quantitative assessment of an allele burden level of a mutant gene and interpretation of results in DNA samples with the burden level of the mutant allele less than 15–20 %, when using Sanger sequencing for analyzing somatic mutations. Applied Biosystems (USA) has developed new software Minor Variant Finder, which allows determining mutations with the allele burden level from 5 %.The objective: to determine the allele burden level and identification of minor variants of somatic mutations in the ASXL1, JAK2 genes and BCR-ABL oncogene using Minor Variant Finder software in patients with myeloproliferative neoplasms.Materials and methods. The level of mutant allele burden for 15 patients with myeloproliferative neoplasms was determined by the identified mutations using the Minor Variant Finder software, after analysis of point somatic mutations in the ASXL1, JAK2 genes and BCR-ABL oncogene by Sanger sequencing.Results. The allele burden level in all 5 ASXL1-positive samples and BCR-ABL-positive sample was determined as higher than 20 % using the Minor Variant Finder software. The allele burden level in 2 cases was higher than 20 % and in 7 cases lower than 20 %, when we analyzed 9 JAK2-positive samples.Conclusion. Minor Variant Finder software can be used to estimate the allele burden level and to identify minor variants of somatic mutations in the ASXL, JAK2 and BCR-ABL genes.

Topology of the U12-U6atac snRNA complex of the minor spliceosome and binding by NTC-related protein RBM22

ABSTRACTSplicing of precursor messenger RNA is catalyzed by the spliceosome, a dynamic ribonucleoprotein assembly composed of five small nuclear (sn)RNAs and >100 proteins. RNA components catalyze the two transesterification reactions, but proteins perform critical roles in assembly and rearrangement. The catalytic core comprises a paired complex involving U2 and U6 snRNAs for the major form of the spliceosome and U12 and U6atac snRNAs for the minor variant (~0.3% of all spliceosomes in higher eukaryotes); the latter performs identical chemistry, despite limited sequence conservation outside key catalytic elements, and lack of the multi-stem central junction found in the U2-U6 snRNA complex. Here we use solution NMR techniques to show that base pairing patterns of the U12-U6atac snRNA complex of both human and Arabidopsis share key elements with the major spliceosome’s U2-U6 snRNA complex; probing of the single-stranded segment opposing termini of the snRNAs indicates elongation in this region in place of the stacked base pairs at the base of the U6 intramolecular stem loop in the U2-U6 snRNA complex. Binding affinity of RBM22, a protein implicated in remodeling human U2-U6 snRNA prior to catalysis, to U12-U6atac was analyzed by electrophoretic mobility shift assays in which we monitored migration of both protein and RNA components in the same gel. Results indicate that RBM22 binds the U2-U6 and U12-U6atac snRNA complexes specifically and with Kd = 3.5 µM and 8.2 µM, respectively. Similar affinity between RBM22 and each RNA complex suggests that the protein performs the same role in both spliceosomes.

The effect of variant interference on de novo assembly for viral deep sequencing

Viruses have high mutation rates and generally exist as a mixture of variants in biological samples. Next-generation sequencing (NGS) approach has surpassed Sanger for generating long viral sequences, yet how variants affect NGS de novo assembly remains largely unexplored. Our results from >15,000 simulated experiments showed that presence of variants can turn an assembly of one genome into tens to thousands of contigs. This “variant interference” (VI) is highly consistent and reproducible by ten most used de novo assemblers, and occurs independent of genome length, read length, and GC content. The main driver of VI is pairwise identities between viral variants. These findings were further supported by in silico simulations, where selective removal of minor variant reads from clinical datasets allow the “rescue” of full viral genomes from fragmented contigs. These results call for careful interpretation of contigs and contig numbers from de novo assembly in viral deep sequencing.

Current Trends in Applications of Circulatory Microchimerism Detection in Transplantation

Primarily due to recent advances of detection techniques, microchimerism (the proportion of minor variant population is below 1%) has recently gained increasing attention in the field of transplantation. Availability of polymorphic markers, such as deletion insertion or single nucleotide polymorphisms along with a vast array of high sensitivity detection techniques, allow the accurate detection of small quantities of donor- or recipient-related materials. This diagnostic information can improve monitoring of allograft injuries in solid organ transplantations (SOT) as well as facilitate early detection of relapse in allogeneic hematopoietic stem cell transplantation (allo-HSCT). In the present review, genetic marker and detection platform options applicable for microchimerism detection are discussed. Furthermore, current results of relevant clinical studies in the context of microchimerism and SOT or allo-HSCT respectively are also summarized.