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.

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.

In chronic infection, HIV gag-specific CD4+ T cell receptor diversity is higher than CD8+ T cell receptor diversity and is associated with less HIV quasispecies diversity.

Cellular immune responses to Gag correlate with improved HIV viral control. The full extent of cellular immune responses comprise both the number of epitopes recognized by CD4+ and CD8+ T cells, as well as the diversity of the T cell receptor (TCR) repertoire directed against each epitope. The optimal diversity of the responsive TCR repertoire is unclear. Therefore, we evaluated the TCR diversity of CD4+ and CD8+ T cells responding to HIV-1 Gag to determine if TCR diversity correlates with clinical or virologic metrics. Previous studies of TCR repertoires have been limited primarily to CD8+ T cell responses directed against a small number of well-characterized T cell epitopes restricted by specific human leucocyte antigens. We stimulated peripheral blood mononuclear cells from 21chronic HIV-infected individuals overnight with a pool of HIV-1 Gag peptides, followed by sorting of activated CD4+ and CD8+ T cells and TCR deep sequencing. We found Gag-reactive CD8+ T cells to be more oligoclonal, with a few dominant TCRs comprising the bulk of the repertoire, compared to the highly diverse TCR repertoires of Gag-reactive CD4+ T cells. HIV viral sequencing of the same donors revealed that high CD4+ T cell TCR diversity was strongly associated with lower HIV Gag genetic diversity. We conclude that the TCR repertoire of Gag-reactive CD4+ T helper cells display substantial diversity without a clearly dominant circulating TCR clonotype, in contrast to a hierarchy of dominant TCR clonotypes in the Gag-reactive CD8+ T cells, and may serve to limit HIV diversity during chronic infection.IMPORTANCE Human T cells recognize portions of viral proteins bound to host molecules (human leucocyte antigens) on the surface of infected cells. T cells recognize these foreign proteins through their T cell receptors (TCRs), which are formed by the assortment of several available V, D and J genes to create millions of combinations of unique TCRs. We measured the diversity of T cells responding to the HIV Gag protein. We found the CD8+ T cell response is primarily made up of a few dominant unique TCRs whereas the CD4+ T cell subset has a much more diverse repertoire of TCRs. We also found there was less change in the virus sequences in subjects with more diverse TCR repertoires. HIV has a high mutation rate, which allows it to evade the immune response. Our findings describe the characteristics of a virus-specific T cell response that may allow it to limit viral evolution.

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.

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.

Repair of an Attenuated Low-Passage Murine Cytomegalovirus Bacterial Artificial Chromosome Identifies a Novel Spliced Gene Essential for Salivary Gland Tropism.

The cloning of herpesviruses as bacterial artificial chromosomes (BACs) has revolutionized the study of herpesvirus biology, allowing rapid and precise manipulation of viral genomes. Several clinical strains of human cytomegalovirus (HCMV) have been cloned as BACs; however, no low-passage strains of murine CMV (MCMV), which provide a model mimicking these isolates, have been cloned. Here, the low-passage G4 strain of was BAC cloned. G4 carries an m157 gene that does not ligate the natural killer (NK) cell-activating receptor, Ly49H, meaning that unlike laboratory strains of MCMV, this virus replicates well in C57BL/6 mice. This BAC clone exhibited normal replication during acute infection in the spleen and liver but was attenuated for salivary gland tropism. Next-generation sequencing revealed a C-to-A mutation at nucleotide position 188422, located in the 3' untranslated region of sgg1, a spliced gene critical for salivary gland tropism. Repair of this mutation restored tropism for the salivary glands. Transcriptional analysis revealed a novel spliced gene within the sgg1 locus. This small open reading frame (ORF), sgg1.1, starts at the 3' end of the first exon of sgg1 and extends exon 2 of sgg1. This shorter spliced gene is prematurely terminated by the nonsense mutation at nt 188422. Sequence analysis of tissue culture-passaged virus demonstrated that sgg1.1 was stable, although other mutational hot spots were identified. The G4 BAC will allow in vivo studies in a broader range of mice, avoiding the strong NK cell responses seen in B6 mice with other MCMV BAC-derived MCMVs.IMPORTANCE Murine cytomegalovirus (MCMV) is widely used as a model of human CMV (HCMV) infection. However, this model relies on strains of MCMV that have been serially passaged in the laboratory for over four decades. These laboratory strains have been cloned as bacterial artificial chromosomes (BACs), which permits rapid and precise manipulation. Low-passage strains of MCMV add to the utility of the mouse model of HCMV infection but do not exist as cloned BACs. This study describes the first such low-passage MCMV BAC. This BAC-derived G4 was initially attenuated in vivo, with subsequent full genomic sequencing revealing a novel spliced transcript required for salivary gland tropism. These data suggest that MCMV, like HCMV, undergoes tissue culture adaptation that can limit in vivo growth and supports the use of BAC clones as a way of standardizing viral strains and minimizing interlaboratory strain variation.

Deep sequence analysis of HIV adaptation following vertical transmission reveals the impact of immune pressure on the evolution of HIV.

Human immunodeficiency virus (HIV) can adapt to an individual's T cell immune response via genomic mutations that affect antigen recognition and impact disease outcome. These viral adaptations are specific to the host's human leucocyte antigen (HLA) alleles, as these molecules determine which peptides are presented to T cells. As HLA molecules are highly polymorphic at the population level, horizontal transmission events are most commonly between HLA-mismatched donor/recipient pairs, representing new immune selection environments for the transmitted virus. In this study, we utilised a deep sequencing approach to determine the HIV quasispecies in 26 mother-to-child transmission pairs where the potential for founder viruses to be pre-adapted is high due to the pairs being haplo-identical at HLA loci. This scenario allowed the assessment of specific HIV adaptations following transmission in either a non-selective immune environment, due to recipient HLA mismatched to original selecting HLA, or a selective immune environment, mediated by matched donor/recipient HLA. We show that the pattern of reversion or fixation of HIV adaptations following transmission provides insight into the replicative cost, and likely compensatory networks, associated with specific adaptations in vivo. Furthermore, although transmitted viruses were commonly heavily pre-adapted to the child's HLA genotype, we found evidence of de novo post-transmission adaptation, representing new epitopes targeted by the child's T cell response. High-resolution analysis of HIV adaptation is relevant when considering vaccine and cure strategies for individuals exposed to adapted viruses via transmission or reactivated from reservoirs.

Pan-genomic analyses identify key Helicobacter pylori pathogenic loci modified by carcinogenic host microenvironments.

OBJECTIVE: Helicobacter pylori is the strongest risk factor for gastric cancer; however, the majority of infected individuals do not develop disease. Pathological outcomes are mediated by complex interactions among bacterial, host and environmental constituents, and two dietary factors linked with gastric cancer risk are iron deficiency and high salt. We hypothesised that prolonged adaptation of H. pylori to in vivo carcinogenic microenvironments results in genetic modification important for disease. DESIGN: Whole genome sequencing of genetically related H. pylori strains that differ in virulence and targeted H. pylori sequencing following prolonged exposure of bacteria to in vitro carcinogenic conditions were performed. RESULTS: A total of 180 unique single nucleotide polymorphisms (SNPs) were identified among the collective genomes when compared with a reference H. pylori genome. Importantly, common SNPs were identified in isolates harvested from iron-depleted and high salt carcinogenic microenvironments, including an SNP within fur (FurR88H). To investigate the direct role of low iron and/or high salt, H. pylori was continuously cultured in vitro under low iron or high salt conditions to assess fur genetic variation. Exposure to low iron or high salt selected for the FurR88H variant after only 5 days. To extend these results, fur was sequenced in 339 clinical H. pylori strains. Among the isolates examined, 17% (40/232) of strains isolated from patients with premalignant lesions harboured the FurR88H variant, compared with only 6% (6/107) of strains from patients with non-atrophic gastritis alone (p=0.0034). CONCLUSION: These results indicate that specific genetic variation arises within H. pylori strains during in vivo adaptation to conditions conducive for gastric carcinogenesis.

HIV integration sites in latently infected cell lines: evidence of ongoing replication.

BACKGROUND: Assessing the location and frequency of HIV integration sites in latently infected cells can potentially inform our understanding of how HIV persists during combination antiretroviral therapy. We developed a novel high throughput sequencing method to evaluate HIV integration sites in latently infected cell lines to determine whether there was virus replication or clonal expansion in these cell lines observed as multiple integration events at the same position. RESULTS: We modified a previously reported method using random DNA shearing and PCR to allow for high throughput robotic processing to identify the site and frequency of HIV integration in latently infected cell lines. Latently infected cell lines infected with intact virus demonstrated multiple distinct HIV integration sites (28 different sites in U1, 110 in ACH-2 and 117 in J1.1 per 150,000 cells). In contrast, cell lines infected with replication-incompetent viruses (J-Lat cells) demonstrated single integration sites. Following in vitro passaging of the ACH-2 cell line, we observed a significant increase in the frequency of unique HIV integration sites and there were multiple mutations and large deletions in the proviral DNA. When the ACH-2 cell line was cultured with the integrase inhibitor raltegravir, there was a significant decrease in the number of unique HIV integration sites and a transient increase in the frequency of 2-LTR circles consistent with virus replication in these cells. CONCLUSION: Cell lines latently infected with intact HIV demonstrated multiple unique HIV integration sites indicating that these cell lines are not clonal and in the ACH-2 cell line there was evidence of low level virus replication. These findings have implications for the use of latently infected cell lines as models of HIV latency and for the use of these cells as standards.

Do variations in the HLA-E ligand encoded by UL40 distinguish individuals susceptible to HCMV disease?

Human cytomegalovirus (HCMV) is carried lifelong by ∼80 % of adults worldwide, generating distinct disease syndromes in transplant recipients, people with HIV (PWH) and neonates. Amino acids 15-23 encoded by the HCMV gene UL40 match positions 3-11 of HLA-A and HLA-C, and constitute a "signal peptide" able to stabilise cell surface HLA-E as a restriction element and a ligand of NKG2A and NKG2C. We present next generation sequencing of UL40 amplified from 15 Australian renal transplant recipients (RTR), six healthy adults and four neonates, and 21 Indonesian PWH. We found no groupwise associations between the presence of multiple sequences and HCMV burden (highest in PWH) or HCMV-associated symptoms in neonates. Homology between UL40 and corresponding HLA-C and HLA-A peptides in 11 RTR revealed perfect matches with HLA-C in three individuals, all carrying HCMV encoding only VMAPRTLIL - a peptide previously associated with viremia. However indices of the burden of HCMV did not segregate in our cohort.

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.

Abacavir inhibits but does not cause self-reactivity to HLA-B*57:01-restricted EBV specific T cell receptors.

Pre-existing pathogen-specific memory T cell responses can contribute to multiple adverse outcomes including autoimmunity and drug hypersensitivity. How the specificity of the T cell receptor (TCR) is subverted or seconded in many of these diseases remains unclear. Here, we apply abacavir hypersensitivity (AHS) as a model to address this question because the disease is linked to memory T cell responses and the HLA risk allele, HLA-B*57:01, and the initiating insult, abacavir, are known. To investigate the role of pathogen-specific TCR specificity in mediating AHS we performed a genome-wide screen for HLA-B*57:01 restricted T cell responses to Epstein-Barr virus (EBV), one of the most prevalent human pathogens. T cell epitope mapping revealed HLA-B*57:01 restricted responses to 17 EBV open reading frames and identified an epitope encoded by EBNA3C. Using these data, we cloned the dominant TCR for EBNA3C and a previously defined epitope within EBNA3B. TCR specificity to each epitope was confirmed, however, cloned TCRs did not cross-react with abacavir plus self-peptide. Nevertheless, abacavir inhibited TCR interactions with their cognate ligands, demonstrating that TCR specificity may be subverted by a drug molecule. These results provide an experimental road map for future studies addressing the heterologous immune responses of TCRs including T cell mediated adverse drug reactions.

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.

Visual Genomics Analysis Studio as a Tool to Analyze Multiomic Data.

Type B adverse drug reactions (ADRs) are iatrogenic immune-mediated syndromes with mechanistic etiologies that remain incompletely understood. Some of the most severe ADRs, including delayed drug hypersensitivity reactions, are T-cell mediated, restricted by specific human leukocyte antigen risk alleles and sometimes by public or oligoclonal T-cell receptors (TCRs), central to the immunopathogenesis of tissue-damaging response. However, the specific cellular signatures of effector, regulatory, and accessory immune populations that mediate disease, define reaction phenotype, and determine severity have not been defined. Recent development of single-cell platforms bringing together advances in genomics and immunology provides the tools to simultaneously examine the full transcriptome, TCRs, and surface protein markers of highly heterogeneous immune cell populations at the site of the pathological response at a single-cell level. However, the requirement for advanced bioinformatics expertise and computational hardware and software has often limited the ability of investigators with the understanding of diseases and biological models to exploit these new approaches. Here we describe the features and use of a state-of-the-art, fully integrated application for analysis and visualization of multiomic single-cell data called Visual Genomics Analysis Studio (VGAS). This unique user-friendly, Windows-based graphical user interface is specifically designed to enable investigators to interrogate their own data. While VGAS also includes tools for sequence alignment and identification of associations with host or organism genetic polymorphisms, in this review we focus on its application for analysis of single-cell TCR-RNA-Cellular Indexing of Transcriptomes and Epitopes by Sequencing (CITE)-seq, enabling holistic cellular characterization by unbiased transcriptome and select surface proteome. Critically, VGAS does not require user-directed coding or access to high-performance computers, instead incorporating performance-optimized hidden code to provide application-based fast and intuitive tools for data analyses and production of high-resolution publication-ready graphics on standard specification laptops. Specifically, it allows analyses of comprehensive single-cell TCR sequencing (scTCR-seq) data, detailing (i) functional pairings of α-ß heterodimer TCRs, (ii) one-click histograms to display entropy and gene rearrangements, and (iii) Circos and Sankey plots to visualize clonality and dominance. For unbiased single-cell RNA sequencing (scRNA-seq) analyses, users extract cell transcriptome signatures according to global structure via principal component analysis, t-distributed stochastic neighborhood embedding, or uniform manifold approximation and projection plots, with overlay of scTCR-seq enabling identification and selection of the immunodominant TCR-expressing populations. Further integration with similar sequence-based detection of surface protein markers using oligo-labeled antibodies (CITE-seq) provides comparative understanding of surface protein expression, with differential gene or protein analyses visualized using volcano plot or heatmap functions. These data can be compared to reference cell atlases or suitable controls to reveal discrete disease-specific subsets, from epithelial to tissue-resident memory T-cells, and activation status, from senescence through exhaustion, with more finite transcript expression displayed as violin and box plots. Importantly, guided tutorial videos are available, as are regular application updates based on the latest advances in bioinformatics and user feedback.

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 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.

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.

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.