Immunodominant T cell epitopes preferentially targeted in multiple individuals are the critical element of successful vaccines and targeted immunotherapies. However, the underlying principles of this ..."convergence" of adaptive immunity among different individuals remain poorly understood. To quantitatively describe epitope immunogenicity, here we propose a supervised machine learning framework generating probabilistic estimates of immunogenicity, termed "immunogenicity scores," based on the numerical features computed through sequence-based simulation approximating the molecular scanning process of peptides presented onto major histocompatibility complex (MHC) by the human T cell receptor (TCR) repertoire. Notably, overlapping sets of intermolecular interaction parameters were commonly utilized in MHC-I and MHC-II prediction. Moreover, a similar simulation of individual TCR-peptide interaction using the same set of interaction parameters yielded correlates of TCR affinity. Pathogen-derived epitopes and tumor-associated epitopes with positive T cell reactivity generally had higher immunogenicity scores than non-immunogenic counterparts, whereas thymically expressed self-epitopes were assigned relatively low scores regardless of their immunogenicity annotation. Immunogenicity score dynamics among single amino acid mutants delineated the landscape of position- and residue-specific mutational impacts. Simulation of position-specific immunogenicity score dynamics detected residues with high escape potential in multiple epitopes, consistent with known escape mutations in the literature. This study indicates that targeting of epitopes by human adaptive immunity is to some extent directed by defined thermodynamic principles. The proposed framework also has a practical implication in that it may enable to more efficiently prioritize epitope candidates highly prone to T cell recognition in multiple individuals, warranting prospective validation across different cohorts.
HIV-associated neurocognitive disorder (HAND) remains an important and yet potentially underdiagnosed manifestation despite the fact that the modern combination antiretroviral therapy (cART) has ...achieved effective viral suppression and greatly reduced the incidence of life-threatening events. Although HIV neurotoxicity is thought to play a central role, the potential of viral genetic signature as diagnostic and/or prognostic biomarker has yet to be fully explored.
Using a manually curated sequence metadataset (80 specimens, 2349 sequences), we demonstrated that only three genetic features are sufficient to predict HAND status regardless of sampling tissues; the accuracy reached 100 and 94% in the hold-out testing subdataset and the entire dataset, respectively. The three genetic features stratified HAND into four distinct clusters. Extrapolating the classification to the 1619 specimens registered in the Los Alamos HIV Sequence Database, the global HAND prevalence was estimated to be 46%, with significant regional variations (30-71%). The R package HANDPrediction was implemented to ensure public availability of key codes.
Our analysis revealed three amino acid positions in gp120 glycoprotein, providing the basis of the development of novel cART regimens specifically optimized for HAND-associated quasispecies. Moreover, the classifier can readily be translated into a diagnostic biomarker, warranting prospective validation.
The pathophysiology of adverse events following programmed cell death protein 1 (PD-1) blockade, including tuberculosis (TB) and autoimmunity, remains poorly characterized. We studied a patient with ...inherited PD-1 deficiency and TB who died of pulmonary autoimmunity. The patient's leukocytes did not express PD-1 or respond to PD-1-mediated suppression. The patient's lymphocytes produced only small amounts of interferon (IFN)-γ upon mycobacterial stimuli, similarly to patients with inborn errors of IFN-γ production who are vulnerable to TB. This phenotype resulted from a combined depletion of Vδ2
γδ T, mucosal-associated invariant T and CD56
natural killer lymphocytes and dysfunction of other T lymphocyte subsets. Moreover, the patient displayed hepatosplenomegaly and an expansion of total, activated and RORγT
CD4
CD8
double-negative αβ T cells, similar to patients with STAT3 gain-of-function mutations who display lymphoproliferative autoimmunity. This phenotype resulted from excessive amounts of STAT3-activating cytokines interleukin (IL)-6 and IL-23 produced by activated T lymphocytes and monocytes, and the STAT3-dependent expression of RORγT by activated T lymphocytes. Our work highlights the indispensable role of human PD-1 in governing both antimycobacterial immunity and self-tolerance, while identifying potentially actionable molecular targets for the diagnostic and therapeutic management of TB and autoimmunity in patients on PD-1 blockade.
High-dimensional cytometry is a powerful technique for deciphering the immunopathological factors common to multiple individuals. However, rational comparisons of multiple batches of experiments ...performed on different occasions or at different sites are challenging because of batch effects. In this study, we describe the integration of multibatch cytometry datasets (iMUBAC), a flexible, scalable, and robust computational framework for unsupervised cell-type identification across multiple batches of high-dimensional cytometry datasets, even without technical replicates. After overlaying cells from multiple healthy controls across batches, iMUBAC learns batch-specific cell-type classification boundaries and identifies aberrant immunophenotypes in patient samples from multiple batches in a unified manner. We illustrate unbiased and streamlined immunophenotyping using both public and in-house mass cytometry and spectral flow cytometry datasets. The method is available as the R package iMUBAC (https://github.com/casanova-lab/iMUBAC).
Antibodies cross-reactive to pathogens and autoantigens are considered pivotal in both infection control and accompanying autoimmunity. However, the pathogenic roles of autoantibodies largely remain ...elusive without a priori knowledge of disease-specific autoantigens. Here, through a novel quantitative proteogenomics approach, we demonstrated a successful identification of immunoglobulin variable heavy chain (VH) sequences highly enriched in pathological immune complex from clinical specimens obtained from a patient with hepatitis C virus-induced cryoglobulinemia (HCV-CG). Reconstructed single-domain antibodies were reactive to both HCV antigens and potentially liver-derived human proteins. Moreover, over the course of antiviral therapy, a substantial "de-evolution" of a distinct sub-repertoire was discovered, to which proteomically identified cryoprecipitation-prone autoantibodies belonged. This sub-repertoire was characterized by IGHJ6*03-derived, long, hydrophobic complementarity determining region (CDR-H3). This study provides a proof-of-concept of de novo mining of autoantibodies and corresponding autoantigen candidates in a disease-specific context in human, thus facilitating future reverse-translational research for the discovery of novel biomarkers and the development of antigen-specific immunotherapy against various autoantibody-related disorders.
Since 2017, hepatitis A virus (HAV) infection has been an epidemic among men who have sex with men (MSM) in Japan. We have come across 11 MSM patients with hepatitis A who were also infected with ...HIV. In 1999–2000, we came across 5 HIV-infected patients with hepatitis A. Since the conditions of current HIV-infected patients have changed owing to the recent progress in anti-HIV therapies, we compared clinical features of hepatitis A between patients in 2017–2018 and those in 1999–2000. By comparing the background characteristics of the patients, we found that the CD4/CD8 ratio was significantly higher in the 2017–2018 group. After the onset of hepatitis, peak levels of hepatic transaminases were found to be higher in the 2017–2018 group, suggesting severe hepatocellular damage. In contrast, neither the peak level of total bilirubin nor the nadir of prothrombin time was significantly different among the 2 groups. We also analyzed the HAV genome derived from some of the recently infected patients, and found that the HAV strains were almost the same among these patients; slight differences were observed from the previously identified strain. Thus, we concluded that the recovery of immunity by recent anti-HIV therapies may result in more severe hepatocellular damages and differences in clinical features.
Natural immunoglobulin M (IgM) is reactive to autoantigens and is believed to be important for autoimmunity. Blood pentameric IgM loaded with antigens forms a large immune complex (IC) that contains ...various elements, including apoptosis inhibitor of macrophage (AIM). Here we demonstrate that this IgM-AIM association contributes to autoantibody production under obese conditions. In mice fed a high-fat diet, natural IgM increased through B cell TLR4 stimulation. AIM associated with IgM and protected AIM from renal excretion, increasing blood AIM levels along with the obesity-induced IgM augmentation. Meanwhile, the AIM association inhibited IgM binding to the Fcα/μ receptor on splenic follicular dendritic cells, thereby protecting the IgM IC from Fcα/μ receptor-mediated internalization. This supported IgM-dependent autoantigen presentation to B cells, stimulating IgG autoantibody production. Accordingly, in obese AIM-deficient (AIM−/−) mice, the increase of multiple IgG autoantibodies observed in obese wild-type mice was abrogated. Thus, the AIM-IgM association plays a critical role in the obesity-associated autoimmune process.
Display omitted
► AIM is protected from renal excretion via binding to IgM pentamers in blood ► Association with AIM interferes with Fcα/μ receptor-mediated internalization of IgM ► IgM-dependent autoantigen presentation on FDCs is preserved by AIM ► Absence of AIM tempers obesity-associated multiple IgG autoantibody production
In addition to inducing metabolic and cardiovascular diseases, obesity also increases blood immunoglobulin M (IgM) levels, followed by multiple autoantibody expansion leading to autoimmune disease. Here, Arai, Miyazaki, and colleagues show that apoptosis inhibitor of macrophage (AIM) associates with IgM and abrogates IgM internalization through the Fcα/μ receptor. This response supports IgM-dependent autoantigen presentation to B cells on splenic follicular dendritic cells, thereby enhancing the development of high-affinity IgG autoantibody-producing plasma cells. Thus, AIM suppression may be a promising candidate for the next generation of therapy for preventing obesity-associated autoimmune diseases.
Pre-existing low-frequency resistance-associated variants (RAVs) may jeopardize successful sustained virological responses (SVR) to HCV treatment with direct-acting antivirals (DAAs). However, the ...potential impact of low-frequency (∼0.1%) mutations, concatenated mutations (haplotypes), and their association with genotypes (Gts) on the treatment outcome has not yet been elucidated, most probably owing to the difficulty in detecting pre-existing minor haplotypes with sufficient length and accuracy. Herein, we characterize a methodological framework based on Illumina MiSeq next-generation sequencing (NGS) coupled with bioinformatics of quasispecies reconstruction (QSR) to realize highly accurate variant calling and genotype-haplotype detection. The core-to-NS3 protease coding sequences in 10 HCV monoinfected patients, 5 of whom had a history of blood transfusion, and 11 HCV/HIV coinfected patients with hemophilia, were studied. Simulation experiments showed that, for minor variants constituting more than 1%, our framework achieved a positive predictive value (PPV) of 100% and sensitivities of 91.7-100% for genotyping and 80.6% for RAV screening. Genotyping analysis indicated the prevalence of dominant Gt1a infection in coinfected patients (6/11 vs 0/10, p = 0.01). For clinical samples, minor genotype overlapping infection was prevalent in HCV/HIV coinfected hemophiliacs (10/11) and patients who experienced whole-blood transfusion (4/5) but none in patients without exposure to blood (0/5). As for RAV screening, the Q80K/R and S122K/R variants were particularly prevalent among minor RAVs observed, detected in 12/21 and 6/21 cases, respectively. Q80K was detected only in coinfected patients, whereas Q80R was predominantly detected in monoinfected patients (1/11 vs 7/10, p < 0.01). Multivariate interdependence analysis revealed the previously unrecognized prevalence of Gt1b-Q80K, in HCV/HIV coinfected hemophiliacs Odds ratio = 13.4 (3.48-51.9), p < 0.01. Our study revealed the distinct characteristics of viral quasispecies between the subgroups specified above and the feasibility of NGS and QSR-based genetic deconvolution of pre-existing minor Gts, RAVs, and their interrelationships.
The human genetic basis of tuberculosis (TB) has long remained elusive. We recently reported a high level of enrichment in homozygosity for the common TYK2 P1104A variant in a heterogeneous cohort of ...patients with TB from non-European countries in which TB is endemic. This variant is homozygous in ∼1/600 Europeans and ∼1/5,000 people from other countries outside East Asia and sub-Saharan Africa. We report a study of this variant in the UK Biobank cohort. The frequency of P1104A homozygotes was much higher in patients with TB (6/620, 1%) than in controls (228/114,473, 0.2%), with an odds ratio (OR) adjusted for ancestry of 5.0 95% confidence interval (CI): 1.96–10.31, P = 2 × 10−3. Conversely,we did not observe enrichment for P1104A heterozygosity, or for TYK2 I684S or V362F homozygosity or heterozygosity. Moreover, it is unlikely that more than 10% of controls were infected with Mycobacterium tuberculosis, as 97% were of European genetic ancestry, born between 1939 and 1970, and resided in the United Kingdom. Had all of them been infected, the OR for developing TB upon infection would be higher. These findings suggest that homozygosity for TYK2 P1104A may account for ∼1% of TB cases in Europeans.