Cmai: Predicting Antigen-Antibody Interactions from Massive Sequencing Data.

The interaction between antigens and antibodies (B cell receptors, BCRs) is the key step underlying the function of the humoral immune system in various biological contexts. The capability to profile the landscape of antigen-binding affinity of a vast number of BCRs will provide a powerful tool to reveal novel insights at unprecedented levels and will yield powerful tools for translational development. However, current experimental approaches for profiling antibody-antigen interactions are costly and time-consuming, and can only achieve low-to-mid throughput. On the other hand, bioinformatics tools in the field of antibody informatics mostly focus on optimization of antibodies given known binding antigens, which is a very different research question and of limited scope. In this work, we developed an innovative Artificial Intelligence tool, Cmai, to address the prediction of the binding between antibodies and antigens that can be scaled to high-throughput sequencing data. Cmai achieved an AUROC of 0.91 in our validation cohort. We devised a biomarker metric based on the output from Cmai applied to high-throughput BCR sequencing data. We found that, during immune-related adverse events (irAEs) caused by immune-checkpoint inhibitor (ICI) treatment, the humoral immunity is preferentially responsive to intracellular antigens from the organs affected by the irAEs. In contrast, extracellular antigens on malignant tumor cells are inducing B cell infiltrations, and the infiltrating B cells have a greater tendency to co-localize with tumor cells expressing these antigens. We further found that the abundance of tumor antigen-targeting antibodies is predictive of ICI treatment response. Overall, Cmai and our biomarker approach filled in a gap that is not addressed by current antibody optimization works nor works such as AlphaFold3 that predict the structures of complexes of proteins that are known to bind.

Correlation between ABCB1 and OLIG2 polymorphisms and the severity and prognosis of patients with cerebral infarction.

This study investigated the relationship between ATP-binding cassette sub-family B member 1 (ABCB1) and OLIG2 single nucleotide polymorphism (SNP) and neurological injury severity and outcome in cerebral infarction (CI). The neurological injury severity of 298 CI patients was evaluated by the National Institutes of Health Stroke Scale. The prognosis of CI patients at 30 days after admission was evaluated by the modified Rankin Scale. And 322 healthy people were selected as the control group. The SNPs of the ABCB1 gene (rs1045642) and OLIG2 gene (rs1059004 and rs9653711) were detected by TaqMan probe PCR, and the distribution of SNPs genotype was analyzed. SNP rs9653711 was correlated with CI. Recessive models of rs1045642 and rs9653711 were correlated with CI. The genotypes of rs1045642 and rs9653711 and genetic models were associated with CI severity. rs1045642 had no correlation with CI prognosis, while rs9653711 had less correlation. The genotype distribution and recessive model were associated with CI prognosis. SNP rs1059004 was not associated with CI severity and prognosis. Alcohol consumption, hypertension, diabetes, hyperlipidemia, and high levels of homocysteine (HCY) were independent risk factors for CI, while hypertension, hyperlipidemia, and HCY were associated with poor prognosis of CI. ABCB1 rs1045642 and OLOG2 rs9653711 are associated with CI severity.

pan-MHC and cross-Species Prediction of T Cell Receptor-Antigen Binding.

Profiling the binding of T cell receptors (TCRs) of T cells to antigenic peptides presented by MHC proteins is one of the most important unsolved problems in modern immunology. Experimental methods to probe TCR-antigen interactions are slow, labor-intensive, costly, and yield moderate throughput. To address this problem, we developed pMTnet-omni, an Artificial Intelligence (AI) system based on hybrid protein sequence and structure information, to predict the pairing of TCRs of αß T cells with peptide-MHC complexes (pMHCs). pMTnet-omni is capable of handling peptides presented by both class I and II pMHCs, and capable of handling both human and mouse TCR-pMHC pairs, through information sharing enabled this hybrid design. pMTnet-omni achieves a high overall Area Under the Curve of Receiver Operator Characteristics (AUROC) of 0.888, which surpasses competing tools by a large margin. We showed that pMTnet-omni can distinguish binding affinity of TCRs with similar sequences. Across a range of datasets from various biological contexts, pMTnet-omni characterized the longitudinal evolution and spatial heterogeneity of TCR-pMHC interactions and their functional impact. We successfully developed a biomarker based on pMTnet-omni for predicting immune-related adverse events of immune checkpoint inhibitor (ICI) treatment in a cohort of 57 ICI-treated patients. pMTnet-omni represents a major advance towards developing a clinically usable AI system for TCR-pMHC pairing prediction that can aid the design and implementation of TCR-based immunotherapeutics.

A short plus long-amplicon based sequencing approach improves genomic coverage and variant detection in the SARS-CoV-2 genome.

High viral transmission in the COVID-19 pandemic has enabled SARS-CoV-2 to acquire new mutations that may impact genome sequencing methods. The ARTIC.v3 primer pool that amplifies short amplicons in a multiplex-PCR reaction is one of the most widely used methods for sequencing the SARS-CoV-2 genome. We observed that some genomic intervals are poorly captured with ARTIC primers. To improve the genomic coverage and variant detection across these intervals, we designed long amplicon primers and evaluated the performance of a short (ARTIC) plus long amplicon (MRL) sequencing approach. Sequencing assays were optimized on VR-1986D-ATCC RNA followed by sequencing of nasopharyngeal swab specimens from fifteen COVID-19 positive patients. ARTIC data covered 94.47% of the virus genome fraction in the positive control and patient samples. Variant analysis in the ARTIC data detected 217 mutations, including 209 single nucleotide variants (SNVs) and eight insertions & deletions. On the other hand, long-amplicon data detected 156 mutations, of which 80% were concordant with ARTIC data. Combined analysis of ARTIC + MRL data improved the genomic coverage to 97.03% and identified 214 high confidence mutations. The combined final set of 214 mutations included 203 SNVs, 8 deletions and 3 insertions. Analysis showed 26 SARS-CoV-2 lineage defining mutations including 4 known variants of concern K417N, E484K, N501Y, P618H in spike gene. Hybrid analysis identified 7 nonsynonymous and 5 synonymous mutations across the genome that were either ambiguous or not called in ARTIC data. For example, G172V mutation in the ORF3a protein and A2A mutation in Membrane protein were missed by the ARTIC assay. Thus, we show that while the short amplicon (ARTIC) assay provides good genomic coverage with high throughput, complementation of poorly captured intervals with long amplicon data can significantly improve SARS-CoV-2 genomic coverage and variant detection.

IRF1 governs the differential interferon-stimulated gene responses in human monocytes and macrophages by regulating chromatin accessibility.

Myeloid lineage cells use TLRs to recognize and respond to diverse microbial ligands. Although unique transcription factors dictate the outcome of specific TLR signaling, whether lineage-specific differences exist to further modulate the quality of TLR-induced inflammation remains unclear. Comprehensive analysis of global gene transcription in human monocytes, monocyte-derived macrophages, and monocyte-derived dendritic cells stimulated with various TLR ligands identifies multiple lineage-specific, TLR-responsive gene programs. Monocytes are hyperresponsive to TLR7/8 stimulation that correlates with the higher expression of the receptors. While macrophages and monocytes express similar levels of TLR4, macrophages, but not monocytes, upregulate interferon-stimulated genes (ISGs) in response to TLR4 stimulation. We find that TLR4 signaling in macrophages uniquely engages transcription factor IRF1, which facilitates the opening of ISG loci for transcription. This study provides a critical mechanistic basis for lineage-specific TLR responses and uncovers IRF1 as a master regulator for the ISG transcriptional program in human macrophages.

Deep sequencing reveals a DAP1 regulatory haplotype that potentiates autoimmunity in systemic lupus erythematosus.

BACKGROUND: Systemic lupus erythematosus (SLE) is a clinically heterogeneous autoimmune disease characterized by the development of anti-nuclear antibodies. Susceptibility to SLE is multifactorial, with a combination of genetic and environmental risk factors contributing to disease development. Like other polygenic diseases, a significant proportion of estimated SLE heritability is not accounted for by common disease alleles analyzed by SNP array-based GWASs. Death-associated protein 1 (DAP1) was implicated as a candidate gene in a previous familial linkage study of SLE and rheumatoid arthritis, but the association has not been explored further. RESULTS: We perform deep sequencing across the DAP1 genomic segment in 2032 SLE patients, and healthy controls, and discover a low-frequency functional haplotype strongly associated with SLE risk in multiple ethnicities. We find multiple cis-eQTLs embedded in a risk haplotype that progressively downregulates DAP1 transcription in immune cells. Decreased DAP1 transcription results in reduced DAP1 protein in peripheral blood mononuclear cells, monocytes, and lymphoblastoid cell lines, leading to enhanced autophagic flux in immune cells expressing the DAP1 risk haplotype. Patients with DAP1 risk allele exhibit significantly higher autoantibody titers and altered expression of the immune system, autophagy, and apoptosis pathway transcripts, indicating that the DAP1 risk allele mediates enhanced autophagy, leading to the survival of autoreactive lymphocytes and increased autoantibody. CONCLUSIONS: We demonstrate how targeted sequencing captures low-frequency functional risk alleles that are missed by SNP array-based studies. SLE patients with the DAP1 genotype have distinct autoantibody and transcription profiles, supporting the dissection of SLE heterogeneity by genetic analysis.

Co-circulation dynamics and persistence of newly introduced clades of 2012 outbreak associated West Nile Virus in Texas, 2012-2015.

The second largest outbreak of West Nile encephalitis and West Nile fever ever recorded occurred in the United States (U.S) in the summer of 2012. The outbreak was related to the widespread circulation of closely related clades, or groups, of West Nile virus (WNV) into multiple states where they were not previously found. Whether the invading 2012 strains were able to circulate and overwinter in states with their own endemic population of WNV is unknown and the effect of viral genetics on adaptation and persistence in a new ecological niche is unclear. In this study, we sequenced 70 mosquito isolates from multiple counties throughout Texas in 2012-2015. We identified isolates representative of previously described 2012 WNV groups (Groups 8-10) and discovered a novel group which we called Group 11. Although we identified isolates representative of WNV endemic (2/70) to Texas, most isolates (68/70) were related to the invading 2012 strains, and of these Group 10 (45/68) was predominant. We also observed differences among the 2012 WNV groups correlating to their genotype. Group 10 WNV in Texas, which carry two putative positively selected variants, had limited introductions into Texas, wide circulation, and strong evidence of continued persistence perhaps indicative of overwintering. In contrast, Groups 8 and 11, without positively selected variants, had multiple introductions into Texas, limited circulation, and limited persistence. Lastly, we identified a potential transmission source in New York for incoming Group 8 WNV into Texas. Altogether our study suggests that mutations in the WNV genome may influence the range and dynamics of WNV circulation, and the ability of different strains to persist in new ecological niches.

The Flexibility of an Amorphous Cobalt Hydroxide Nanomaterial Promotes the Electrocatalysis of Oxygen Evolution Reaction.

Structural flexibility can be a desirable trait of an operating catalyst because it adapts itself to a given catalytic process for enhanced activity. Here, amorphous cobalt hydroxide nanocages are demonstrated to be a promising electrocatalyst with an overpotential of 0.28 V at 10 mA cm-2 , far outperforming the crystalline counterparts and being in the top rank of the catalysts of their kind, under the condition of electrocatalytic oxygen evolution reaction. From the direct experimental in situ and ex situ results, this enhanced activity is attributed to its high structural flexibility in terms of 1) facile and holistic transformation into catalytic active phase; 2) hosting oxygen vacancies; and 3) structure self-regulation in a real-time process. Significantly, based on plausible catalytic mechanism and computational simulation results, it is disclosed how this structural flexibility facilitates the kinetics of oxygen evolution reaction. This work deepens the understanding of the structure-activity relationship of the Co-based catalysts in electrochemical catalysis, and it inspires more applications that require flexible structures enabled by such amorphous nanomaterials.

Regulatory polymorphisms modulate the expression of HLA class II molecules and promote autoimmunity.

Targeted sequencing of sixteen SLE risk loci among 1349 Caucasian cases and controls produced a comprehensive dataset of the variations causing susceptibility to systemic lupus erythematosus (SLE). Two independent disease association signals in the HLA-D region identified two regulatory regions containing 3562 polymorphisms that modified thirty-seven transcription factor binding sites. These extensive functional variations are a new and potent facet of HLA polymorphism. Variations modifying the consensus binding motifs of IRF4 and CTCF in the XL9 regulatory complex modified the transcription of HLA-DRB1, HLA-DQA1 and HLA-DQB1 in a chromosome-specific manner, resulting in a 2.5-fold increase in the surface expression of HLA-DR and DQ molecules on dendritic cells with SLE risk genotypes, which increases to over 4-fold after stimulation. Similar analyses of fifteen other SLE risk loci identified 1206 functional variants tightly linked with disease-associated SNPs and demonstrated that common disease alleles contain multiple causal variants modulating multiple immune system genes.

DGKE variants cause a glomerular microangiopathy that mimics membranoproliferative GN.

Renal microangiopathies and membranoproliferative GN (MPGN) can manifest similar clinical presentations and histology, suggesting the possibility of a common underlying mechanism in some cases. Here, we performed homozygosity mapping and whole exome sequencing in a Turkish consanguineous family and identified DGKE gene variants as the cause of a membranoproliferative-like glomerular microangiopathy. Furthermore, we identified two additional DGKE variants in a cohort of 142 unrelated patients diagnosed with membranoproliferative GN. This gene encodes the diacylglycerol kinase DGKε, which is an intracellular lipid kinase that phosphorylates diacylglycerol to phosphatidic acid. Immunofluorescence confocal microscopy demonstrated that mouse and rat Dgkε colocalizes with the podocyte marker WT1 but not with the endothelial marker CD31. Patch-clamp experiments in human embryonic kidney (HEK293) cells showed that DGKε variants affect the intracellular concentration of diacylglycerol. Taken together, these results not only identify a genetic cause of a glomerular microangiopathy but also suggest that the phosphatidylinositol cycle, which requires DGKE, is critical to the normal function of podocytes.

Association of two independent functional risk haplotypes in TNIP1 with systemic lupus erythematosus.

OBJECTIVE: Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by autoantibody production and altered type I interferon expression. Genetic surveys and genome-wide association studies have identified >30 SLE susceptibility genes. One of these genes, TNIP1, encodes the ABIN1 protein. ABIN1 functions in the immune system by restricting NF-κB signaling. The present study was undertaken to investigate the genetic factors that influence association with SLE in genes that regulate the NF-κB pathway. METHODS: We analyzed a dense set of genetic markers spanning TNIP1 and TAX1BP1, as well as the TNIP1 homolog TNIP2, in case-control populations of diverse ethnic origins. TNIP1, TNIP2, and TAX1BP1 were fine-mapped in a total of 8,372 SLE cases and 7,492 healthy controls from European-ancestry, African American, Hispanic, East Asian, and African American Gullah populations. Levels of TNIP1 messenger RNA (mRNA) and ABIN1 protein in Epstein-Barr virus-transformed human B cell lines were analyzed by quantitative reverse transcription-polymerase chain reaction and Western blotting, respectively. RESULTS: We found significant associations between SLE and genetic variants within TNIP1, but not in TNIP2 or TAX1BP1. After resequencing and imputation, we identified 2 independent risk haplotypes within TNIP1 in individuals of European ancestry that were also present in African American and Hispanic populations. Levels of TNIP1 mRNA and ABIN1 protein were reduced among subjects with these haplotypes, suggesting that they harbor hypomorphic functional variants that influence susceptibility to SLE by restricting ABIN1 expression. CONCLUSION: Our results confirm the association signals between SLE and TNIP1 variants in multiple populations and provide new insight into the mechanism by which TNIP1 variants may contribute to SLE pathogenesis.

Identification of IRF8, TMEM39A, and IKZF3-ZPBP2 as susceptibility loci for systemic lupus erythematosus in a large-scale multiracial replication study.

Systemic lupus erythematosus (SLE) is a chronic heterogeneous autoimmune disorder characterized by the loss of tolerance to self-antigens and dysregulated interferon responses. The etiology of SLE is complex, involving both heritable and environmental factors. Candidate-gene studies and genome-wide association (GWA) scans have been successful in identifying new loci that contribute to disease susceptibility; however, much of the heritable risk has yet to be identified. In this study, we sought to replicate 1,580 variants showing suggestive association with SLE in a previously published GWA scan of European Americans; we tested a multiethnic population consisting of 7,998 SLE cases and 7,492 controls of European, African American, Asian, Hispanic, Gullah, and Amerindian ancestry to find association with the disease. Several genes relevant to immunological pathways showed association with SLE. Three loci exceeded the genome-wide significance threshold: interferon regulatory factor 8 (IRF8; rs11644034; p(meta-Euro) = 2.08 × 10(-10)), transmembrane protein 39A (TMEM39A; rs1132200; p(meta-all) = 8.62 × 10(-9)), and 17q21 (rs1453560; p(meta-all) = 3.48 × 10(-10)) between IKAROS family of zinc finger 3 (AIOLOS; IKZF3) and zona pellucida binding protein 2 (ZPBP2). Fine mapping, resequencing, imputation, and haplotype analysis of IRF8 indicated that three independent effects tagged by rs8046526, rs450443, and rs4843869, respectively, were required for risk in individuals of European ancestry. Eleven additional replicated effects (5 × 10(-8) < p(meta-Euro) < 9.99 × 10(-5)) were observed with CFHR1, CADM2, LOC730109/IL12A, LPP, LOC63920, SLU7, ADAMTSL1, C10orf64, OR8D4, FAM19A2, and STXBP6. The results of this study increase the number of confirmed SLE risk loci and identify others warranting further investigation.

A Tlr7 translocation accelerates systemic autoimmunity in murine lupus.

The y-linked autoimmune accelerating (yaa) locus is a potent autoimmune disease allele. Transcription profiling of yaa-bearing B cells revealed the overexpression of a cluster of X-linked genes that included Tlr7. FISH analysis demonstrated the translocation of this segment onto the yaa chromosome. The resulting overexpression of Tlr7 increased in vitro responses to Toll-like receptor (TLR) 7 signaling in all yaa-bearing males. B6.yaa mice are not overtly autoimmune, but the addition of Sle1, which contains the autoimmune-predisposing Slam/Cd2 haplotype, causes the development of fatal lupus with numerous immunological aberrations. B6.Sle1yaa CD4 T cells develop the molecular signature for T(FH) cells and also show expression changes in numerous cytokines and chemokines. Disease development and all component autoimmune phenotypes were inhibited by Sles1, a potent suppressor locus. Sles1 had no effect on yaa-enhanced TLR7 signaling in vitro, and these data place Sles1 downstream from the lesion in innate immune responses mediated by TLR7, suggesting that Sles1 modulates the activation of adaptive immunity in response to innate immune signaling.

Sle1ab mediates the aberrant activation of STAT3 and Ras-ERK signaling pathways in B lymphocytes.

The Sle1ab genomic interval on murine chromosome 1 mediates the loss of immune tolerance to chromatin resulting in antinuclear Abs (ANA) production in the lupus-prone NZM2410 mouse. Global gene expression analysis was used to identify the molecular pathways that are dysregulated at the initiation of B lymphocyte autoimmunity in B6.Sle1ab mice. This analysis identified that STAT3 and ras-ERK signaling pathways are aberrantly activated in Sle1ab B lymphocytes, consistent with increased production of IL-6 by splenic B lymphocytes and monocytes in B6.Sle1ab mice. In vitro treatment of splenic mononuclear cells isolated from ANA-positive Sle1ab mice with anti-IL-6 Ab or AG490, an inhibitor of STAT3 signaling pathway, suppressed ANA production in short-term culture, indicating that this pathway was essential to the production of autoantibodies. In vivo treatment of ANA-positive B6.Sle1ab mice with the ras pathway inhibitor, perillyl alcohol, suppressed the increase of ANA. These findings identify IL-6 as a early key cytokine in Sle1ab-mediated disease development and indicate that the STAT3 and ras-ERK signaling pathways are potential therapeutic targets for treating systemic lupus erythematosus.