Epistatic interaction between ERAP2 and HLA modulates HIV-1 adaptation and disease outcome in an Australian population.

A strong genetic predictor of outcome following untreated HIV-1 infection is the carriage of specific alleles of human leukocyte antigens (HLAs) that present viral epitopes to T cells. Residual variation in outcome measures may be attributed, in part, to viral adaptation to HLA-restricted T cell responses. Variants of the endoplasmic reticulum aminopeptidases (ERAPs) influence the repertoire of T cell epitopes presented by HLA alleles as they trim pathogen-derived peptide precursors to optimal lengths for antigen presentation, along with other functions unrelated to antigen presentation. We investigated whether ERAP variants influence HLA-associated HIV-1 adaptation with demonstrable effects on overall HIV-1 disease outcome. Utilizing host and viral data of 249 West Australian individuals with HIV-1 subtype B infection, we identified a novel association between two linked ERAP2 single nucleotide polymorphisms (SNPs; rs2248374 and rs2549782) with plasma HIV RNA concentration (viral load) (P adjusted = 0.0024 for both SNPs). Greater HLA-associated HIV-1 adaptation in the HIV-1 Gag gene correlated significantly with higher viral load, lower CD4+ T cell count and proportion; P = 0.0103, P = 0.0061, P = 0.0061, respectively). When considered together, there was a significant interaction between the two ERAP2 SNPs and HLA-associated HIV-1 adaptation on viral load (P = 0.0111). In a comprehensive multivariate model, addition of ERAP2 haplotypes and HLA associated adaptation as an interaction term to known HLA and CCR5 determinants and demographic factors, increased the explanatory variance of population viral load from 17.67% to 45.1% in this dataset. These effects were not replicated in publicly available datasets with comparably sized cohorts, suggesting that any true global epistasis may be dependent on specific HLA-ERAP allelic combinations. Our data raises the possibility that ERAP2 variants may shape peptide repertoires presented to HLA class I-restricted T cells to modulate the degree of viral adaptation within individuals, in turn contributing to disease variability at the population level. Analyses of other populations and experimental studies, ideally with locally derived ERAP genotyping and HLA-specific viral adaptations are needed to elucidate this further.

Multiomic single-cell sequencing defines tissue-specific responses in Stevens-Johnson Syndrome and Toxic epidermal necrolysis.

Stevens-Johnson syndrome and toxic epidermal necrolysis (SJS/TEN) is a rare but life-threatening cutaneous drug reaction mediated by human leukocyte antigen (HLA) class I-restricted CD8+ T-cells. To obtain an unbiased assessment of SJS/TEN cellular immunopathogenesis, we performed single-cell (sc) transcriptome, surface proteome, and TCR sequencing on unaffected skin, affected skin, and blister fluid from 17 SJS/TEN patients. From 119,784 total cells, we identified 16 scRNA-defined subsets, confirmed by subset-defining surface protein expression. Keratinocytes upregulated HLA and IFN-response genes in the affected skin. Cytotoxic CD8+ T-cell subpopulations of expanded and unexpanded TCRαß clonotypes were shared in affected skin and blister fluid but absent or unexpanded in SJS/TEN unaffected skin. SJS/TEN blister fluid is a rich reservoir of oligoclonal CD8+ T-cells with an effector phenotype driving SJS/TEN pathogenesis. This multiomic database will act as the basis to define antigen-reactivity, HLA restriction, and signatures of drug-antigen-reactive T-cell clonotypes at a tissue level.

Underrepresentation of activating KIR gene expression in single-cell RNA-seq data is due to KIR gene misassignment.

Standard single-cell RNA-sequencing alignment pipelines exhibit a propensity for misassigning killer immunoglobulin-like receptor (KIR) transcripts, thereby giving rise to inaccuracies in quantifying KIR expression. Alves et al. elucidated that these default workflows frequently misclassify activating KIR transcripts as inhibitory KIR expression, resulting in a skewed representation of the KIR repertoire.

Tracking of activated cTfh cells following sequential influenza vaccinations reveals transcriptional profile of clonotypes driving a vaccine-induced immune response.

Introduction: A vaccine against influenza is available seasonally but is not 100% effective. A predictor of successful seroconversion in adults is an increase in activated circulating T follicular helper (cTfh) cells after vaccination. However, the impact of repeated annual vaccinations on long-term protection and seasonal vaccine efficacy remains unclear. Methods: In this study, we examined the T cell receptor (TCR) repertoire and transcriptional profile of vaccine-induced expanded cTfh cells in individuals who received sequential seasonal influenza vaccines. We measured the magnitude of cTfh and plasmablast cell activation from day 0 (d0) to d7 post-vaccination as an indicator of a vaccine response. To assess TCR diversity and T cell expansion we sorted activated and resting cTfh cells at d0 and d7 post-vaccination and performed TCR sequencing. We also single cell sorted activated and resting cTfh cells for TCR analysis and transcriptome sequencing. Results and discussion: The percent of activated cTfh cells significantly increased from d0 to d7 in each of the 2016-17 (p < 0.0001) and 2017-18 (p = 0.015) vaccine seasons with the magnitude of cTfh activation increase positively correlated with the frequency of circulating plasmablast cells in the 2016-17 (p = 0.0001) and 2017-18 (p = 0.003) seasons. At d7 post-vaccination, higher magnitudes of cTfh activation were associated with increased clonality of cTfh TCR repertoire. The TCRs from vaccine-expanded clonotypes were identified and tracked longitudinally with several TCRs found to be present in both years. The transcriptomic profile of these expanded cTfh cells at the single cell level demonstrated overrepresentation of transcripts of genes involved in the type-I interferon pathway, pathways involved in gene expression, and antigen presentation and recognition. These results identify the expansion and transcriptomic profile of vaccine-induced cTfh cells important for B cell help.

Subset-specific Retention of Donor Myeloid Cells After Major Histocompatibility Complex-matched and Mismatched Liver Transplantation.

BACKGROUND: During solid organ transplantation, donor leukocytes, including myeloid cells, are transferred within the organ to the recipient. Both tolerogenic and alloreactive roles have been attributed to donor myeloid cells; however, their subset-specific retention posttransplantation has not been investigated in detail. METHODS: Major histocompatibility complex (MHC)-matched and mismatched liver transplants were performed in mice, and the fate of donor and recipient myeloid cells was assessed. RESULTS: Following MHC-matched transplantation, a proportion of donor myeloid cells was retained in the graft, whereas others egressed and persisted in the blood, spleen, and bone marrow but not the lymph nodes. In contrast, after MHC-mismatched transplantation, all donor myeloid cells, except Kupffer cells, were depleted. This depletion was caused by recipient T and B cells because all donor myeloid subsets were retained in MHC-mismatched grafts when recipients lacked T and B cells. Recipient myeloid cells rapidly infiltrated MHC-matched and, to a greater extent, MHC-mismatched liver grafts. MHC-mismatched grafts underwent a transient rejection episode on day 7, coinciding with a transition in macrophages to a regulatory phenotype, after which rejection resolved. CONCLUSIONS: Phenotypic and kinetic differences in the myeloid cell responses between MHC-matched and mismatched grafts were identified. A detailed understanding of the dynamics of immune responses to transplantation is critical to improving graft outcomes.

Identification of human progenitors of exhausted CD8+ T cells associated with elevated IFN-γ response in early phase of viral infection.

T cell exhaustion is a hallmark of hepatitis C virus (HCV) infection and limits protective immunity in chronic viral infections and cancer. Limited knowledge exists of the initial viral and immune dynamics that characterise exhaustion in humans. We studied longitudinal blood samples from a unique cohort of individuals with primary infection using single-cell multi-omics to identify the functions and phenotypes of HCV-specific CD8+ T cells. Early elevated IFN-γ response against the transmitted virus is associated with the rate of immune escape, larger clonal expansion, and early onset of exhaustion. Irrespective of disease outcome, we find heterogeneous subsets of progenitors of exhaustion, based on the level of PD-1 expression and loss of AP-1 transcription factors. Intra-clonal analysis shows distinct trajectories with multiple fates and evolutionary plasticity of precursor cells. These findings challenge the current paradigm on the contribution of CD8+ T cells to HCV disease outcome and provide data for future studies on T cell differentiation in human infections.

Adaptation to HLA-associated immune pressure over the course of HIV infection and in circulating HIV-1 strains.

Adaptation to human leukocyte antigen (HLA)-associated immune pressure represents a major driver of human immunodeficiency virus (HIV) evolution at both the individual and population level. To date, there has been limited exploration of the impact of the initial cellular immune response in driving viral adaptation, the dynamics of these changes during infection and their effect on circulating transmitting viruses at the population level. Capturing detailed virological and immunological data from acute and early HIV infection is challenging as this commonly precedes the diagnosis of HIV infection, potentially by many years. In addition, rapid initiation of antiretroviral treatment following a diagnosis is the standard of care, and central to global efforts towards HIV elimination. Yet, acute untreated infection is the critical period in which the diversity of proviral reservoirs is first established within individuals, and associated with greater risk of onward transmissions in a population. Characterizing the viral adaptations evident in the earliest phases of infection, coinciding with the initial cellular immune responses is therefore relevant to understanding which changes are of greatest impact to HIV evolution at the population level. In this study, we utilized three separate cohorts to examine the initial CD8+ T cell immune response to HIV (cross-sectional acute infection cohort), track HIV evolution in response to CD8+ T cell-mediated immunity over time (longitudinal chronic infection cohort) and translate the impact of HLA-driven HIV evolution to the population level (cross-sectional HIV sequence data spanning 30 years). Using next generation viral sequencing and enzyme-linked immunospot interferon-gamma recall responses to peptides representing HLA class I-specific HIV T cell targets, we observed that CD8+ T cell responses can select viral adaptations prior to full antibody seroconversion. Using the longitudinal cohort, we uncover that viral adaptations have the propensity to be retained over time in a non-selective immune environment, which reflects the increasing proportion of pre-adapted HIV strains within the Western Australian population over an approximate 30-year period.

Variants of HCMV UL18 Sequenced Directly from Clinical Specimens Associate with Antibody and T-Cell Responses to HCMV.

Around 80% of adults worldwide carry human cytomegaloviris (HCMV). The HCMV gene UL18 is a homolog of HLA class I genes and encodes a protein with high affinity for the NK and T-cell cytotoxicity inhibitor LIR-1. UL18 was deep sequenced from blood, saliva or urine from Indonesian people with HIV (PWH) (n = 28), Australian renal transplant recipients (RTR) (n = 21), healthy adults (n = 7) and neonates (n = 4). 95% of samples contained more than one variant of HCMV UL18, as defined by carriage of nonsynonymous variations. When aligned with immunological markers of the host's burden of HCMV, the S318N variation associated with high levels of antibody reactive with HCMV lysate in PWH over 12 months on antiretroviral therapy. The A107T variation associated with HCMV antibody levels and inflammatory biomarkers in PWH at early timepoints. Variants D32G, D248N, V250A and E252D aligned with elevated HCMV antibody levels in RTR, while M191K, E196Q and F165L were associated with HCMV-reactive T-cells and proportions of Vδ2- γδ T-cells-populations linked with high burdens of HCMV. We conclude that UL18 is a highly variable gene, where variation may alter the persistent burden of HCMV and/or the host response to that burden.

Altered subgenomic RNA abundance provides unique insight into SARS-CoV-2 B.1.1.7/Alpha variant infections.

B.1.1.7 lineage SARS-CoV-2 is more transmissible, leads to greater clinical severity, and results in modest reductions in antibody neutralization. Subgenomic RNA (sgRNA) is produced by discontinuous transcription of the SARS-CoV-2 genome. Applying our tool (periscope) to ARTIC Network Oxford Nanopore Technologies genomic sequencing data from 4400 SARS-CoV-2 positive clinical samples, we show that normalised sgRNA is significantly increased in B.1.1.7 (alpha) infections (n = 879). This increase is seen over the previous dominant lineage in the UK, B.1.177 (n = 943), which is independent of genomic reads, E cycle threshold and days since symptom onset at sampling. A noncanonical sgRNA which could represent ORF9b is found in 98.4% of B.1.1.7 SARS-CoV-2 infections compared with only 13.8% of other lineages, with a 16-fold increase in median sgRNA abundance. We demonstrate that ORF9b protein levels are increased 6-fold in B.1.1.7 compared to a B lineage virus in vitro. We hypothesise that increased ORF9b in B.1.1.7 is a direct consequence of a triple nucleotide mutation in nucleocapsid (28280:GAT > CAT, D3L) creating a transcription regulatory-like sequence complementary to a region 3' of the genomic leader. These findings provide a unique insight into the biology of B.1.1.7 and support monitoring of sgRNA profiles to evaluate emerging potential variants of concern.

Sequencing of the Viral UL111a Gene Directly from Clinical Specimens Reveals Variants of HCMV-Encoded IL-10 That Are Associated with Altered Immune Responses to HCMV.

Human cytomegalovirus (HCMV) is a beta-herpesvirus carried by ~80% of adults worldwide. Acute infections are often asymptomatic in healthy individuals but generate diverse syndromes in neonates, renal transplant recipients (RTR), and people with HIV (PWH). The HCMV gene UL111a encodes a homolog of human interleukin-10 (IL-10) that interacts with the human IL-10 receptor. Deep sequencing technologies were used to sequence UL111a directly from 59 clinical samples from Indonesian PWH and Australian RTR, healthy adults, and neonates. Overall, 93% of samples contained more than one variant of HCMV, as defined by at least one nonsynonymous variation. Carriage of these variants differed between neonates and adults, Australians and Indonesians, and between saliva and blood leukocytes. The variant alleles of N41D and S71Y occurred together in Australian RTR and were associated with higher T-cell responses to HCMV pp65. The variant P122S was associated with lower levels of antibodies reactive with a lysate of HCMV-infected fibroblasts. L174F was associated with increased levels of antibodies reactive with HCMV lysate, immediate-early 1 (IE-1), and glycoprotein B (gB) in Australian RTR and Indonesians PWH, suggesting a higher viral burden. We conclude that variants of UL111a are common in all populations and may influence systemic responses to HCMV.

On the Evolutionary Trajectory of SARS-CoV-2: Host Immunity as a Driver of Adaptation in RNA Viruses.

Host immunity can exert a complex array of selective pressures on a pathogen, which can drive highly mutable RNA viruses towards viral escape. The plasticity of a virus depends on its rate of mutation, as well as the balance of fitness cost and benefit of mutations, including viral adaptations to the host's immune response. Since its emergence, SARS-CoV-2 has diversified into genetically distinct variants, which are characterised often by clusters of mutations that bolster its capacity to escape human innate and adaptive immunity. Such viral escape is well documented in the context of other pandemic RNA viruses such as the human immunodeficiency virus (HIV) and influenza virus. This review describes the selection pressures the host's antiviral immunity exerts on SARS-CoV-2 and other RNA viruses, resulting in divergence of viral strains into more adapted forms. As RNA viruses obscure themselves from host immunity, they uncover weak points in their own armoury that can inform more comprehensive, long-lasting, and potentially cross-protective vaccine coverage.

Cross-Reactivity to Mutated Viral Immune Targets Can Influence CD8+ T Cell Functionality: An Alternative Viral Adaptation Strategy.

Loss of T cell immunogenicity due to mutations in virally encoded epitopes is a well-described adaptation strategy to limit host anti-viral immunity. Another described, but less understood, adaptation strategy involves the selection of mutations within epitopes that retain immune recognition, suggesting a benefit for the virus despite continued immune pressure (termed non-classical adaptation). To understand this adaptation strategy, we utilized a single cell transcriptomic approach to identify features of the HIV-specific CD8+ T cell responses targeting non-adapted (NAE) and adapted (AE) forms of epitopes containing a non-classical adaptation. T cell receptor (TCR) repertoire and transcriptome were obtained from antigen-specific CD8+ T cells of chronic (n=7) and acute (n=4) HIV-infected subjects identified by either HLA class I tetramers or upregulation of activation markers following peptide stimulation. CD8+ T cells were predominantly dual tetramer+, confirming a large proportion of cross-reactive TCR clonotypes capable of recognizing the NAE and AE form. However, single-reactive CD8+ T cells were identified in acute HIV-infected subjects only, providing the potential for the selection of T cell clones over time. The transcriptomic profile of CD8+ T cells was dependent on the autologous virus: subjects whose virus encoded the NAE form of the epitope (and who transitioned to the AE form at a later timepoint) exhibited an 'effective' immune response, as indicated by expression of transcripts associated with polyfunctionality, cytotoxicity and apoptosis (largely driven by the genes GZMB, IFNÉ£, CCL3, CCL4 and CCL5). These data suggest that viral adaptation at a single amino acid residue can provide an alternative strategy for viral survival by modulating the transcriptome of CD8+ T cells and potentially selecting for less effective T cell clones from the acute to chronic phase.

Sequencing Directly from Clinical Specimens Reveals Genetic Variations in HCMV-Encoded Chemokine Receptor US28 That May Influence Antibody Levels and Interactions with Human Chemokines.

Human cytomegalovirus (HCMV) is a beta-herpesvirus carried by ∼80% of the world's population. Acute infections are asymptomatic in healthy individuals but generate diverse syndromes in neonates, solid organ transplant recipients, and HIV-infected individuals. The HCMV gene US28 encodes a homolog of a human chemokine receptor that is able to bind several chemokines and HIV gp120. Deep sequencing technologies were used to sequence US28 directly from 60 clinical samples from Indonesian HIV patients and Australian renal transplant recipients, healthy adults, and neonates. Molecular modeling approaches were used to predict whether nine nonsynonymous mutations in US28 may alter protein binding to a panel of six chemokines and two variants of HIV gp120. Ninety-two percent of samples contained more than one variant of HCMV, as defined by at least one nonsynonymous mutation. Carriage of these variants differed between neonates and adults, Australian and Indonesian samples, and saliva samples and blood leukocytes. Two nonsynonymous mutations (N170D and R267K) were associated with increased levels of immediate early protein 1 (IE-1) and glycoprotein B (gB) HCMV-reactive antibodies, suggesting a higher viral burden. Seven of the nine mutations were predicted to alter binding of at least one ligand. Overall, HCMV variants are common in all populations and have the potential to affect US28 interactions with human chemokines and/or gp120 and alter responses to the virus. The findings relied on deep sequencing technologies applied directly to clinical samples, so the variants exist in vivo. IMPORTANCE Human cytomegalovirus (HCMV) is a common viral pathogen of solid organ transplant recipients, neonates, and HIV-infected individuals. HCMV encodes homologs of several host genes with the potential to influence viral persistence and/or pathogenesis. Here, we present deep sequencing of an HCMV chemokine receptor homolog, US28, acquired directly from clinical specimens. Carriage of these variants differed between patient groups and was associated with different levels of circulating HCMV-reactive antibodies. These features are consistent with a role for US28 in HCMV persistence and pathogenesis. This was supported by in silico analyses of the variant sequences demonstrating altered ligand-binding profiles. The data delineate a novel approach to understanding the pathogenesis of HCMV and may impact the development of an effective vaccine.

Generation of a Novel SARS-CoV-2 Sub-genomic RNA Due to the R203K/G204R Variant in Nucleocapsid: Homologous Recombination has Potential to Change SARS-CoV-2 at Both Protein and RNA Level.

BACKGROUND: Genetic variations across the SARS-CoV-2 genome may influence transmissibility of the virus and the host's anti-viral immune response, in turn affecting the frequency of variants over time. In this study, we examined the adjacent amino acid polymorphisms in the nucleocapsid (R203K/G204R) of SARS-CoV-2 that arose on the background of the spike D614G change and describe how strains harboring these changes became dominant circulating strains globally. METHODS: Deep-sequencing data of SARS-CoV-2 from public databases and from clinical samples were analyzed to identify and map genetic variants and sub-genomic RNA transcripts across the genome. Results: Sequence analysis suggests that the 3 adjacent nucleotide changes that result in the K203/R204 variant have arisen by homologous recombination from the core sequence of the leader transcription-regulating sequence (TRS) rather than by stepwise mutation. The resulting sequence changes generate a novel sub-genomic RNA transcript for the C-terminal dimerization domain of nucleocapsid. Deep-sequencing data from 981 clinical samples confirmed the presence of the novel TRS-CS-dimerization domain RNA in individuals with the K203/R204 variant. Quantification of sub-genomic RNA indicates that viruses with the K203/R204 variant may also have increased expression of sub-genomic RNA from other open reading frames. CONCLUSIONS: The finding that homologous recombination from the TRS may have occurred since the introduction of SARS-CoV-2 in humans, resulting in both coding changes and novel sub-genomic RNA transcripts, suggests this as a mechanism for diversification and adaptation within its new host.

Manipulating the NKG2D Receptor-Ligand Axis Using CRISPR: Novel Technologies for Improved Host Immunity.

The activating immune receptor natural killer group member D (NKG2D) and its cognate ligands represent a fundamental surveillance system of cellular distress, damage or transformation. Signaling through the NKG2D receptor-ligand axis is critical for early detection of viral infection or oncogenic transformation and the presence of functional NKG2D ligands (NKG2D-L) is associated with tumor rejection and viral clearance. Many viruses and tumors have developed mechanisms to evade NKG2D recognition via transcriptional, post-transcriptional or post-translational interference with NKG2D-L, supporting the concept that circumventing immune evasion of the NKG2D receptor-ligand axis may be an attractive therapeutic avenue for antiviral therapy or cancer immunotherapy. To date, the complexity of the NKG2D receptor-ligand axis and the lack of specificity of current NKG2D-targeting therapies has not allowed for the precise manipulation required to optimally harness NKG2D-mediated immunity. However, with the discovery of clustered regularly interspaced short palindromic repeats (CRISPRs) and CRISPR-associated (Cas) proteins, novel opportunities have arisen in the realm of locus-specific gene editing and regulation. Here, we give a brief overview of the NKG2D receptor-ligand axis in humans and discuss the levels at which NKG2D-L are regulated and dysregulated during viral infection and oncogenesis. Moreover, we explore the potential for CRISPR-based technologies to provide novel therapeutic avenues to improve and maximize NKG2D-mediated immunity.

Reprogramming the anti-tumor immune response via CRISPR genetic and epigenetic editing.

Precise clustered regularly interspaced short palindromic repeats (CRISPR)-mediated genetic and epigenetic manipulation of the immune response has become a promising immunotherapeutic approach toward combating tumorigenesis and tumor progression. CRISPR-based immunologic reprograming in cancer therapy comprises the locus-specific enhancement of host immunity, the improvement of tumor immunogenicity, and the suppression of tumor immunoevasion. To date, the ex vivo re-engineering of immune cells directed to inhibit the expression of immune checkpoints or to express synthetic immune receptors (chimeric antigen receptor therapy) has shown success in some settings, such as in the treatment of melanoma, lymphoma, liver, and lung cancer. However, advancements in nuclease-deactivated CRISPR-associated nuclease-9 (dCas9)-mediated transcriptional activation or repression and Cas13-directed gene suppression present novel avenues for the development of tumor immunotherapies. In this review, the basis for development, mechanism of action, and outcomes from recently published Cas9-based clinical trial (genetic editing) and dCas9/Cas13-based pre-clinical (epigenetic editing) data are discussed. Lastly, we review cancer immunotherapy-specific considerations and barriers surrounding use of these approaches in the clinic.

Dynamic changes to tissue-resident immunity after MHC-matched and MHC-mismatched solid organ transplantation.

The heterogeneous pool of tissue-resident lymphocytes in solid organs mediates infection responses and supports tissue integrity and repair. Their vital functions in normal physiology suggest an important role in solid organ transplantation; however, their detailed examination in this context has not been performed. Here, we report the fate of multiple lymphocyte subsets, including T, B, and innate lymphoid cells, after murine liver and heart transplantation. In major histocompatibility complex (MHC)-matched transplantation, donor lymphocytes are retained in liver grafts and peripheral lymphoid organs of heart and liver transplant recipients. In MHC-mismatched transplantation, increased infiltration of the graft by recipient cells and depletion of donor lymphocytes occur, which can be prevented by removal of recipient T and B cells. Recipient lymphocytes fail to recreate the native organs' phenotypically diverse tissue-resident lymphocyte composition, even in MHC-matched models. These post-transplant changes may leave grafts vulnerable to infection and impair long-term graft function.

Generation of a novel SARS-CoV-2 sub-genomic RNA due to the R203K/G204R variant in nucleocapsid: homologous recombination has potential to change SARS-CoV-2 at both protein and RNA level.

BACKGROUND: Genetic variations across the SARS-CoV-2 genome may influence transmissibility of the virus and the host’s anti-viral immune response, in turn affecting the frequency of variants over-time. In this study, we examined the adjacent amino acid polymorphisms in the nucleocapsid (R203K/G204R) of SARS-CoV-2 that arose on the background of the spike D614G change and describe how strains harboring these changes became dominant circulating strains globally. METHODS: Deep sequencing data of SARS-CoV-2 from public databases and from clinical samples were analyzed to identify and map genetic variants and sub-genomic RNA transcripts across the genome. RESULTS: Sequence analysis suggests that the three adjacent nucleotide changes that result in the K203/R204 variant have arisen by homologous recombination from the core sequence (CS) of the leader transcription-regulating sequence (TRS) rather than by stepwise mutation. The resulting sequence changes generate a novel sub-genomic RNA transcript for the C-terminal dimerization domain of nucleocapsid. Deep sequencing data from 981 clinical samples confirmed the presence of the novel TRS-CS-dimerization domain RNA in individuals with the K203/R204 variant. Quantification of sub-genomic RNA indicates that viruses with the K203/R204 variant may also have increased expression of sub-genomic RNA from other open reading frames. CONCLUSIONS: The finding that homologous recombination from the TRS may have occurred since the introduction of SARS-CoV-2 in humans resulting in both coding changes and novel sub-genomic RNA transcripts suggests this as a mechanism for diversification and adaptation within its new host.

Subgenomic RNA identification in SARS-CoV-2 genomic sequencing data.

We have developed periscope, a tool for the detection and quantification of subgenomic RNA (sgRNA) in SARS-CoV-2 genomic sequence data. The translation of the SARS-CoV-2 RNA genome for most open reading frames (ORFs) occurs via RNA intermediates termed "subgenomic RNAs." sgRNAs are produced through discontinuous transcription, which relies on homology between transcription regulatory sequences (TRS-B) upstream of the ORF start codons and that of the TRS-L, which is located in the 5' UTR. TRS-L is immediately preceded by a leader sequence. This leader sequence is therefore found at the 5' end of all sgRNA. We applied periscope to 1155 SARS-CoV-2 genomes from Sheffield, United Kingdom, and validated our findings using orthogonal data sets and in vitro cell systems. By using a simple local alignment to detect reads that contain the leader sequence, we were able to identify and quantify reads arising from canonical and noncanonical sgRNA. We were able to detect all canonical sgRNAs at the expected abundances, with the exception of ORF10. A number of recurrent noncanonical sgRNAs are detected. We show that the results are reproducible using technical replicates and determine the optimum number of reads for sgRNA analysis. In VeroE6 ACE2+/- cell lines, periscope can detect the changes in the kinetics of sgRNA in orthogonal sequencing data sets. Finally, variants found in genomic RNA are transmitted to sgRNAs with high fidelity in most cases. This tool can be applied to all sequenced COVID-19 samples worldwide to provide comprehensive analysis of SARS-CoV-2 sgRNA.