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
+/- 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.
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.
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.
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.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
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.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Abacavir hypersensitivity syndrome (ABC HSS) is strongly associated with carriage of human leukocyte antigen (HLA)-B*57:01, which has a 100% negative predictive value for the development of ABC HSS. ...However, 45% of individuals who carry HLA-B*57:01 can tolerate ABC. We investigated immune and non-immune related genes in ABC HSS (n = 95) and ABC tolerant (n = 43) HLA-B*57:01 + patients to determine other factors required for the development of ABC HSS. Assignment of phenotype showed that ABC HSS subjects were significantly less likely than tolerants to carry only ERAP1 hypoactive trimming allotypes (p = 0.02). An altered self-peptide repertoire model by which abacavir activates T cells is in keeping with observation that endoplasmic reticulum aminopeptidase 1 (ERAP1) allotypes that favour efficient peptide trimming are more common in ABC HSS patients compared to patients who tolerate ABC. Independently, non-specific immune activation via soluble cluster of differentiation antigen 14 (sCD14) may also influence susceptibility to ABC HSS.
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.
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.
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.
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.
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.
DNA sequence-based typing at the HLA-A, -B, -C, -DPB1, -DQA1, -DQB1, and -DRB1 loci was performed on 496 healthy adult donors from San Diego, California, to characterize allele frequencies in support ...of studies of T cell responses to common allergens. Deviations from Hardy Weinberg proportions were detected at each locus except A and C. Several alleles were found in more than 15% of individuals, including the class II alleles DPB1∗02:01, DPB1∗04:01, DQA1∗01:02, DQA1∗05:01, DQB1∗03:01, and the class I allele A∗02:01. Genotype data will be available in the Allele Frequencies Net Database (AFND 3562).