Quantitative Prediction of the Landscape of T Cell Epitope Immunogenicity in Sequence Space

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

Genome-wide detection of human intronic AG-gain variants located between splicing branchpoints and canonical splice acceptor sites

Human genetic variants that introduce an AG into the intronic region between the branchpoint (BP) and the canonical splice acceptor site (ACC) of protein-coding genes can disrupt pre-mRNA splicing. Using our genome-wide BP database, we delineated the BP-ACC segments of all human introns and found extreme depletion of AG/YAG in the [BP+8, ACC-4] high-risk region. We developed AGAIN as a genome-wide computational approach to systematically and precisely pinpoint intronic AG-gain variants within the BP-ACC regions. AGAIN identified 350 AG-gain variants from the Human Gene Mutation Database, all of which alter splicing and cause disease. Among them, 74% created new acceptor sites, whereas 31% resulted in complete exon skipping. AGAIN also predicts the protein-level products resulting from these two consequences. We performed AGAIN on our exome/genomes database of patients with severe infectious diseases but without known genetic etiology and identified a private homozygous intronic AG-gain variant in the antimycobacterial gene SPPL2A in a patient with mycobacterial disease. AGAIN also predicts a retention of six intronic nucleotides that encode an in-frame stop codon, turning AG-gain into stop-gain. This allele was then confirmed experimentally to lead to loss of function by disrupting splicing. We further showed that AG-gain variants inside the high-risk region led to misspliced products, while those outside the region did not, by two case studies in genes STAT1 and IRF7. We finally evaluated AGAIN on our 14 paired exome-RNAseq samples and found that 82% of AG-gain variants in high-risk regions showed evidence of missplicing

Human MCTS1-dependent translation of JAK2 is essential for IFN-γ immunity to mycobacteria.

Human inherited disorders of interferon-gamma (IFN-γ) immunity underlie severe mycobacterial diseases. We report X-linked recessive MCTS1 deficiency in men with mycobacterial disease from kindreds of different ancestries (from China, Finland, Iran, and Saudi Arabia). Complete deficiency of this translation re-initiation factor impairs the translation of a subset of proteins, including the kinase JAK2 in all cell types tested, including T lymphocytes and phagocytes. JAK2 expression is sufficiently low to impair cellular responses to interleukin-23 (IL-23) and partially IL-12, but not other JAK2-dependent cytokines. Defective responses to IL-23 preferentially impair the production of IFN-γ by innate-like adaptive mucosal-associated invariant T cells (MAIT) and γδ T lymphocytes upon mycobacterial challenge. Surprisingly, the lack of MCTS1-dependent translation re-initiation and ribosome recycling seems to be otherwise physiologically redundant in these patients. These findings suggest that X-linked recessive human MCTS1 deficiency underlies isolated mycobacterial disease by impairing JAK2 translation in innate-like adaptive T lymphocytes, thereby impairing the IL-23-dependent induction of IFN-γ

Human OTULIN haploinsufficiency impairs cell-intrinsic immunity to staphylococcal alpha-toxin

The molecular basis of interindividual clinical variability upon infection with Staphylococcus aureus is unclear. We describe patients with haploinsufficiency for the linear deubiquitinase OTULIN, encoded by a gene on chromosome 5p. Patients suffer from episodes of life-threatening necrosis, typically triggered by S. aureus infection. The disorder is phenocopied in patients with the 5p- (Cri-du-Chat) chromosomal deletion syndrome. OTULIN haploinsufficiency causes an accumulation of linear ubiquitin in dermal fibroblasts, but tumor necrosis factor receptor-mediated nuclear factor kappa B signaling remains intact. Blood leukocyte subsets are unaffected. The OTULIN-dependent accumulation of caveolin-1 in dermal fibroblasts, but not leukocytes, facilitates the cytotoxic damage inflicted by the staphylococcal virulence factor alpha-toxin. Naturally elicited antibodies against alpha-toxin contribute to incomplete clinical penetrance. Human OTULIN haploinsufficiency underlies life-threatening staphylococcal disease by disrupting cell-intrinsic immunity to alpha-toxin in nonleukocytic cells.Peer reviewe

Impaired IL-23-dependent induction of IFN-gamma underlies mycobacterial disease in patients with inherited TYK2 deficiency

Human cells homozygous for rare loss-of-expression (LOE) TYK2 alleles have impaired, but not abolished, cellular responses to IFN-alpha/beta (underlying viral diseases in the patients) and to IL-12 and IL-23 (underlying mycobacterial diseases). Cells homozygous for the common P1104A TYK2 allele have selectively impaired responses to IL-23 (underlying isolated mycobacterial disease). We report three new forms of TYK2 deficiency in six patients from five families homozygous for rare TYK2 alleles (R864C, G996R, G634E, or G1010D) or compound heterozygous for P1104A and a rare allele (A928V). All these missense alleles encode detectable proteins. The R864C and G1010D alleles are hypomorphic and loss-of-function (LOF), respectively, across signaling pathways. By contrast, hypomorphic G996R, G634E, and A928V mutations selectively impair responses to IL-23, like P1104A. Impairment of the IL-23-dependent induction of IFN-gamma is the only mechanism of mycobacterial disease common to patients with complete TYK2 deficiency with or without TYK2 expression, partial TYK2 deficiency across signaling pathways, or rare or common partial TYK2 deficiency specific for IL-23 signaling.ANRS Nord-Sud ; CIBSS ; CODI ; Comité para el Desarrollo de la Investigación ; Fulbright Future Scholarshi

Inherited Human ITK Deficiency Impairs IFN-γ Immunity and Underlies Tuberculosis

Inborn errors of IFN-γ immunity can underlie tuberculosis (TB). We report three patients from two kindreds without EBV viremia or disease but with severe TB and inherited complete ITK deficiency, a condition associated with severe EBV disease that renders immunological studies challenging. They have CD4+ αβ T lymphocytopenia with a concomitant expansion of CD4-CD8- double-negative (DN) αβ and Vδ2- γδ T lymphocytes, both displaying a unique CD38+CD45RA+T-bet+EOMES- phenotype. Itk-deficient mice recapitulated an expansion of the γδ T and DN αβ T lymphocyte populations in the thymus and spleen, respectively. Moreover, the patients\u27 T lymphocytes secrete small amounts of IFN-γ in response to TCR crosslinking, mitogens, or forced synapse formation with autologous B lymphocytes. Finally, the patients\u27 total lymphocytes secrete small amounts of IFN-γ, and CD4+, CD8+, DN αβ T, Vδ2+ γδ T, and MAIT cells display impaired IFN-γ production in response to BCG. Inherited ITK deficiency undermines the development and function of various IFN-γ-producing T cell subsets, thereby underlying TB

Prediction of HIV-associated neurocognitive disorder (HAND) from three genetic features of envelope gp120 glycoprotein

Abstract Background 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. Results 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. Conclusions 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

Differential proteogenomics applied to HCV-induced cryoglobulinemia

Graphical summaries of the concept and the findings of our study. Illustrations were modified from the resources distributed in the Togo picture gallery (<a href="http://g86.dbcls.jp/~togoriv/">http://g86.dbcls.jp/~togoriv/</a>), licensed under CC-BY 4.0 ©Togo picture gallery by the Database Center for Life Science (DBCLS), Japan

MOESM1 of Prediction of HIV-associated neurocognitive disorder (HAND) from three genetic features of envelope gp120 glycoprotein

Additional file 1

Homozygosity for TYK2 P1104A underlies tuberculosis in about 1% of patients in a cohort of European ancestry

International audienceThe 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