The butyrophilin (BTN) gene family: from milk fat to the regulation of the immune response

Phosphoantigen recognition by Vγ9Vδ2 T cells

Vγ9Vδ2 T cells comprise 1-10% of human peripheral blood T cells. As multifunctional T cells with a strong antimicrobial and antitumor potential, they are of strong interest for immunotherapeutic development. Their hallmark is the eponymous Vγ9Vδ2 T-cell antigen receptor (TCR), which mediates activation by so-called "phosphoantigens" (PAg). PAg are small pyrophosphorylated intermediates of isoprenoid synthesis of microbial or host origin, with the latter elevated in some tumors and after administration of aminobisphosphonates. This review summarizes the progress in understanding PAg-recognition, with emphasis on the interaction between butyrophilins (BTN) and PAg and insights gained by phylogenetic studies on BTNs and Vγ9Vδ2 T cells, especially the comparison of human and alpaca. It proposes a composite ligand model in which BTN3A1-A2/A3-heteromers and BTN2A1 homodimers form a Vγ9Vδ2 TCR activating complex. An initiating step is the binding of PAg to the intracellular BTN3A1-B30.2 domain and formation of a complex with the B30.2 domains of BTN2A1. On the extracellular surface this results in BTN2A1-IgV binding to Vγ9-TCR framework determinants and BTN3A-IgV to additional complementarity determining regions of both TCR chains. Unresolved questions of this model are discussed, as well as questions on the structural basis and the physiological consequences of PAg-recognition.

Cancer immunotherapy by γδ T cells

The premise of cancer immunotherapy is that cancers are specifically visible to an immune system tolerized to healthy self. The promise of cancer immunotherapy is that immune effector mechanisms and immunological memory can jointly eradicate cancers and inoperable metastases and de facto vaccinate against recurrence. For some patients with hitherto incurable diseases, including metastatic melanoma, this promise is being realized by game-changing immunotherapies based on αβ T cells. Today's challenges are to bring benefit to greater numbers of patients of diverse ethnicities, target more cancer types, and achieve a cure while incurring fewer adverse events. In meeting those challenges, specific benefits may be offered by γδ T cells, which compose a second T cell lineage with distinct recognition capabilities and functional traits that bridge innate and adaptive immunity. γδ T cell-based clinical trials, including off-the-shelf adoptive cell therapy and agonist antibodies, are yielding promising results, although identifiable problems remain. In addressing those problems, we advocate that immunotherapies be guided by the distinctive biology of γδ T cells, as elucidated by ongoing research.

Primate-specific BTN3A2 protects against SARS-CoV-2 infection by interacting with and reducing ACE2

BACKGROUND

Coronavirus disease 2019 (COVID-19) is an immune-related disorder caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The complete pathogenesis of the virus remains to be determined. Unraveling the molecular mechanisms governing SARS-CoV-2 interactions with host cells is crucial for the formulation of effective prophylactic measures and the advancement of COVID-19 therapeutics.

METHODS

We analyzed human lung single-cell RNA sequencing dataset to discern the association of butyrophilin subfamily 3 member A2 (BTN3A2) expression with COVID-19. The BTN3A2 gene edited cell lines and transgenic mice were infected by live SARS-CoV-2 in a biosafety level 3 (BSL-3) laboratory. Immunoprecipitation, flow cytometry, biolayer interferometry and competition ELISA assays were performed in BTN3A2 gene edited cells. We performed quantitative real-time PCR, histological and/or immunohistochemical analyses for tissue samples from mice with or without SARS-CoV-2 infection.

FINDINGS

The BTN3A2 mRNA level was correlated with COVID-19 severity. BTN3A2 expression was predominantly identified in epithelial cells, elevated in pathological epithelial cells from COVID-19 patients and co-occurred with ACE2 expression in the same lung cell subtypes. BTN3A2 targeted the early stage of the viral life cycle by inhibiting SARS-CoV-2 attachment through interactions with the receptor-binding domain (RBD) of the Spike protein and ACE2. BTN3A2 inhibited ACE2-mediated SARS-CoV-2 infection by reducing ACE2 in vitro and in vivo.

INTERPRETATION

These results reveal a key role of BTN3A2 in the fight against COVID-19. Identifying potential monoclonal antibodies which mimic BTN3A2 may facilitate disruption of SARS-CoV-2 infection, providing a therapeutic avenue for COVID-19.

FUNDING

This study was supported by the National Natural Science Foundation of China (32070569, U1902215, and 32371017), the CAS "Light of West China" Program, and Yunnan Province (202305AH340006).

Inferring Characteristics of the Tumor Immune Microenvironment of Patients with HNSCC from Single-Cell Transcriptomics of Peripheral Blood

In this study, we explore the possibility of inferring characteristics of the tumor immune microenvironment from the blood. Specifically, we investigate two datasets of patients with head and neck squamous cell carcinoma with matched single-cell RNA sequencing (scRNA-seq) from peripheral blood mononuclear cells (PBMCs) and tumor tissues. Our analysis shows that the immune cell fractions and gene expression profiles of various immune cells within the tumor microenvironment can be inferred from the matched PBMC scRNA-seq data. We find that the established exhausted T-cell signature can be predicted from the blood and serve as a valuable prognostic blood biomarker of immunotherapy response. Additionally, our study reveals that the inferred ratio between tumor memory B- and regulatory T-cell fractions is predictive of immunotherapy response and is superior to the well-established cytolytic and exhausted T-cell signatures. These results highlight the promising potential of PBMC scRNA-seq in cancer immunotherapy and warrant, and will hopefully facilitate, further investigations on a larger scale. The code for predicting tumor immune microenvironment from PBMC scRNA-seq, TIMEP, is provided, offering other researchers the opportunity to investigate its prospective applications in various other indications.

SIGNIFICANCE

Our work offers a new and promising paradigm in liquid biopsies to unlock the power of blood single-cell transcriptomics in cancer immunotherapy.

Prognostic model incorporating immune checkpoint genes to predict the immunotherapy efficacy for lung adenocarcinoma: a cohort study integrating machine learning algorithms

This study aimed to develop and validate a nomogram based on immune checkpoint genes (ICGs) for predicting prognosis and immune checkpoint blockade (ICB) efficacy in lung adenocarcinoma (LUAD) patients. A total of 385 LUAD patients from the TCGA database and 269 LUAD patients in the combined dataset (GSE41272 + GSE50081) were divided into training and validation cohorts, respectively. Three different machine learning algorithms including random forest (RF), least absolute shrinkage and selection operator (LASSO) logistic regression analysis, and support vector machine (SVM) were employed to select the predictive markers from 82 ICGs to construct the prognostic nomogram. The X-tile software was used to stratify patients into high- and low-risk subgroups based on the nomogram-derived risk scores. Differences in functional enrichment and immune infiltration between the two subgroups were assessed using gene set variation analysis (GSVA) and various algorithms. Additionally, three lung cancer cohorts receiving ICB therapy were utilized to evaluate the ability of the model to predict ICB efficacy in the real world. Five ICGs were identified as predictive markers across all three machine learning algorithms, leading to the construction of a nomogram with strong potential for prognosis prediction in both the training and validation cohorts (all AUC values close to 0.800). The patients were divided into high- (risk score ≥ 185.0) and low-risk subgroups (risk score < 185.0). Compared to the high-risk subgroup, the low-risk subgroup exhibited enrichment in immune activation pathways and increased infiltration of activated immune cells, such as CD8 + T cells and M1 macrophages (P < 0.05). Furthermore, the low-risk subgroup had a greater likelihood of benefiting from ICB therapy and longer progression-free survival (PFS) than did the high-risk subgroup (P < 0.05) in the two cohorts receiving ICB therapy. A nomogram based on ICGs was constructed and validated to aid in predicting prognosis and ICB treatment efficacy in LUAD patients.

A comprehensive meta-analysis of tissue resident memory T cells and their roles in shaping immune microenvironment and patient prognosis in non-small cell lung cancer

Tissue-resident memory T cells (TRM) are a specialized subset of long-lived memory T cells that reside in peripheral tissues. However, the impact of TRM-related immunosurveillance on the tumor-immune microenvironment (TIME) and tumor progression across various non-small-cell lung cancer (NSCLC) patient populations is yet to be elucidated. Our comprehensive analysis of multiple independent single-cell and bulk RNA-seq datasets of patient NSCLC samples generated reliable, unique TRM signatures, through which we inferred the abundance of TRM in NSCLC. We discovered that TRM abundance is consistently positively correlated with CD4+ T helper 1 cells, M1 macrophages, and resting dendritic cells in the TIME. In addition, TRM signatures are strongly associated with immune checkpoint and stimulatory genes and the prognosis of NSCLC patients. A TRM-based machine learning model to predict patient survival was validated and an 18-gene risk score was further developed to effectively stratify patients into low-risk and high-risk categories, wherein patients with high-risk scores had significantly lower overall survival than patients with low-risk. The prognostic value of the risk score was independently validated by the Cancer Genome Atlas Program (TCGA) dataset and multiple independent NSCLC patient datasets. Notably, low-risk NSCLC patients with higher TRM infiltration exhibited enhanced T-cell immunity, nature killer cell activation, and other TIME immune responses related pathways, indicating a more active immune profile benefitting from immunotherapy. However, the TRM signature revealed low TRM abundance and a lack of prognostic association among lung squamous cell carcinoma patients in contrast to adenocarcinoma, indicating that the two NSCLC subtypes are driven by distinct TIMEs. Altogether, this study provides valuable insights into the complex interactions between TRM and TIME and their impact on NSCLC patient prognosis. The development of a simplified 18-gene risk score provides a practical prognostic marker for risk stratification.

A Ménage à trois: NLRC5, immunity, and metabolism

The nucleotide-binding and oligomerization domain-like receptors (NLRs) NLR family CARD domain-containing protein 5 (NLRC5) and Class II Major Histocompatibility Complex Transactivator (CIITA) are transcriptional regulators of major histocompatibility complex (MHC) class I and class II genes, respectively. MHC molecules are central players in our immune system, allowing the detection of hazardous 'non-self' antigens and, thus, the recognition and elimination of infected or transformed cells from the organism. Recently, CIITA and NLRC5 have emerged as regulators of selected genes of the butyrophilin (BTN) family that interestingly are located in the extended MHC locus. BTNs are transmembrane proteins exhibiting structural similarities to B7 family co-modulatory molecules. The family member BTN2A2, which indeed contributes to the control of T cell activation, was found to be transcriptionally regulated by CIITA. NLRC5 emerged instead as an important regulator of the BTN3A1, BTN3A2, and BTN3A3 genes. Together with BTN2A1, BTN3As regulate non-conventional Vγ9Vδ2 T cell responses triggered by selected metabolites of microbial origin or accumulating in hematologic cancer cells. Even if endogenous metabolites conform to the canonical definition of 'self', metabolically abnormal cells can represent a danger for the organism and should be recognized and controlled by immune system cells. Collectively, new data on the role of NLRC5 in the expression of BTN3As link the mechanisms regulating canonical 'non-self' presentation and those marking cells with abnormal metabolic configurations for immune recognition, an evolutionary parallel that we discuss in this perspective review.

From immunomodulation to therapeutic prospects: Unveiling the biology of butyrophilins in cancer

Butyrophilin (BTN) proteins are a type of membrane protein that belongs to the Ig superfamily. They exhibit a high degree of structural similarity to molecules in the B7 family. They fulfill a complex function in regulating immune responses, including immunomodulatory roles, as they influence γδ T cells. The biology of BTN molecules indicates that they are capable of inhibiting the immune system's ability to detect antigens within tumors. A dynamic association between BTN molecules and cellular surfaces is also recognized in specific contexts, influencing their biology. Notably, the dynamism of BTN3A1 is associated with the immunosuppression of T cells or the activation of Vγ9Vδ2 T cells. Cancer immunotherapy relies heavily on T cells to modulate immune function within the intricate interaction of the tumor microenvironment (TME). A significant interaction between the TME and antitumor immunity involves the presence of BTN, which should be taken into account when developing immunotherapy. This review explores potential therapeutic applications of BTN molecules, based on the current understanding of their biology.

Interferon-γ and its response are determinants of antibody-mediated rejection and clinical outcomes in patients after renal transplantation

Interferon-γ (IFN-γ) is an important cytokine in tissue homeostasis and immune response, while studies about it in antibody-mediated rejection (ABMR) are very limited. This study aims to comprehensively elucidate the role of IFN-γ in ABMR after renal transplantation. In six renal transplantation cohorts, the IFN-γ responses (IFNGR) biological process was consistently top up-regulated in ABMR compared to stable renal function or even T cell-mediated rejection in both allografts and peripheral blood. According to single-cell analysis, IFNGR levels were found to be broadly elevated in most cell types in allografts and peripheral blood with ABMR. In allografts with ABMR, M1 macrophages had the highest IFNGR levels and were heavily infiltrated, while kidney resident M2 macrophages were nearly absent. In peripheral blood, CD14+ monocytes had the top IFNGR level and were significantly increased in ABMR. Immunofluorescence assay showed that levels of IFN-γ and M1 macrophages were sharply elevated in allografts with ABMR than non-rejection. Importantly, the IFNGR level in allografts was identified as a strong risk factor for long-term renal graft survival. Together, this study systematically analyzed multi-omics from thirteen independent cohorts and identified IFN-γ and IFNGR as determinants of ABMR and clinical outcomes in patients after renal transplantation.

Milk Fat Globules: 2024 Updates

Milk fat globules (MFGs) are a remarkable example of nature's ingenuity. Human milk (HM) carries contains 3-5% fat, 0.8-0.9% protein, 6.9-7.2% carbohydrate calculated as lactose, and 0.2% mineral constituents. Most of these nutrients are carried in these MFGs, which are composed of an energy-rich triacylglycerol (TAG) core surrounded by a triple membrane structure. The membrane contains polar lipids, specialized proteins, glycoproteins, and cholesterol. Each of these bioactive components serves important nutritional, immunological, neurological, and digestive functions. These MFGs are designed to release energy rapidly in the upper gastrointestinal tract and then persist for some time in the gut lumen so that the protective bioactive molecules are conveyed to the colon. These properties may shape the microbial colonization and innate immune properties of the developing gastrointestinal tract. Milk fat globules in milk from humans and ruminants may resemble in structure but there are considerable differences in size, profile, composition, and specific constituents. There are possibilities to not only enhance the nutritional composition in a goal-oriented fashion to correct specific deficiencies in the infant but also to use these fat globules as a nutraceutical in infants who require specific treatments. To mention a few, there might be possibilities in enhancing neurodevelopment, in defense against gastrointestinal and respiratory tract infections, improving insulin sensitivity, treating chronic inflammation, and altering plasma lipids. This review provides an overview of the composition, structure, and biological activities of the various components of the MFGs. We have assimilated research findings from our own laboratory with an extensive review of the literature utilizing key terms in multiple databases including PubMed, EMBASE, and Science Direct. To avoid bias in the identification of studies, keywords were short-listed a priori from anecdotal experience and PubMed's Medical Subject Heading (MeSH) thesaurus.

Significance of an altered lncRNA landscape in schizophrenia and cognition: clues from a case-control association study

Genetic etiology of schizophrenia is poorly understood despite large genome-wide association data. Long non-coding RNAs (lncRNAs) with a probable regulatory role are emerging as important players in neuro-psychiatric disorders including schizophrenia. Prioritising important lncRNAs and analyses of their holistic interaction with their target genes may provide insights into disease biology/etiology. Of the 3843 lncRNA SNPs reported in schizophrenia GWASs extracted using lincSNP 2.0, we prioritised n = 247 based on association strength, minor allele frequency and regulatory potential and mapped them to lncRNAs. lncRNAs were then prioritised based on their expression in brain using lncRBase, epigenetic role using 3D SNP and functional relevance to schizophrenia etiology. 18 SNPs were finally tested for association with schizophrenia (n = 930) and its endophenotypes-tardive dyskinesia (n = 176) and cognition (n = 565) using a case-control approach. Associated SNPs were characterised by ChIP seq, eQTL, and transcription factor binding site (TFBS) data using FeatSNP. Of the eight SNPs significantly associated, rs2072806 in lncRNA hsaLB_IO39983 with regulatory effect on BTN3A2 was associated with schizophrenia (p = 0.006); rs2710323 in hsaLB_IO_2331 with role in dysregulation of ITIH1 with tardive dyskinesia (p < 0.05); and four SNPs with significant cognition score reduction (p < 0.05) in cases. Two of these with two additional variants in eQTL were observed among controls (p < 0.05), acting likely as enhancer SNPs and/or altering TFBS of eQTL mapped downstream genes. This study highlights important lncRNAs in schizophrenia and provides a proof of concept of novel interactions of lncRNAs with protein-coding genes to elicit alterations in immune/inflammatory pathways of schizophrenia.

A distinct topology of BTN3A IgV and B30.2 domains controlled by juxtamembrane regions favors optimal human γδ T cell phosphoantigen sensing

Butyrophilin (BTN)-3A and BTN2A1 molecules control the activation of human Vγ9Vδ2 T cells during T cell receptor (TCR)-mediated sensing of phosphoantigens (PAg) derived from microbes and tumors. However, the molecular rules governing PAg sensing remain largely unknown. Here, we establish three mechanistic principles of PAg-mediated γδ T cell activation. First, in humans, following PAg binding to the intracellular BTN3A1-B30.2 domain, Vγ9Vδ2 TCR triggering involves the extracellular V-domain of BTN3A2/BTN3A3. Moreover, the localization of both protein domains on different chains of the BTN3A homo-or heteromers is essential for efficient PAg-mediated activation. Second, the formation of BTN3A homo-or heteromers, which differ in intracellular trafficking and conformation, is controlled by molecular interactions between the juxtamembrane regions of the BTN3A chains. Finally, the ability of PAg not simply to bind BTN3A-B30.2, but to promote its subsequent interaction with the BTN2A1-B30.2 domain, is essential for T-cell activation. Defining these determinants of cooperation and the division of labor in BTN proteins improves our understanding of PAg sensing and elucidates a mode of action that may apply to other BTN family members.

Role of the viral polymerase during adaptation of influenza A viruses to new hosts

As a group, influenza-A viruses (IAV) infect a wide range of animal hosts, however, they are constrained to infecting selected host species by species-specific interactions between the host and virus, that are required for efficient replication of the viral RNA genome. When IAV cross the species barrier, they acquire mutations in the viral genome to enable interactions with the new host factors, or to compensate for their loss. The viral polymerase genes polymerase basic 1, polymerase basic 2, and polymerase-acidic are important sites of host adaptation. In this review, we discuss why the viral polymerase is so vital to the process of host adaptation, look at some of the known viral mutations, and host factors involved in adaptation, particularly of avian IAV to mammalian hosts.

BTN3A3 evasion promotes the zoonotic potential of influenza A viruses

Spillover events of avian influenza A viruses (IAVs) to humans could represent the first step in a future pandemic1. Several factors that limit the transmission and replication of avian IAVs in mammals have been identified. There are several gaps in our understanding to predict which virus lineages are more likely to cross the species barrier and cause disease in humans1. Here, we identified human BTN3A3 (butyrophilin subfamily 3 member A3)2 as a potent inhibitor of avian IAVs but not human IAVs. We determined that BTN3A3 is expressed in human airways and its antiviral activity evolved in primates. We show that BTN3A3 restriction acts primarily at the early stages of the virus life cycle by inhibiting avian IAV RNA replication. We identified residue 313 in the viral nucleoprotein (NP) as the genetic determinant of BTN3A3 sensitivity (313F or, rarely, 313L in avian viruses) or evasion (313Y or 313V in human viruses). However, avian IAV serotypes, such as H7 and H9, that spilled over into humans also evade BTN3A3 restriction. In these cases, BTN3A3 evasion is due to substitutions (N, H or Q) in NP residue 52 that is adjacent to residue 313 in the NP structure3. Thus, sensitivity or resistance to BTN3A3 is another factor to consider in the risk assessment of the zoonotic potential of avian influenza viruses.

The intestinal γδ T cells: functions in the gut and in the distant organs

Located in the frontline against the largest population of microbiota, the intestinal mucosa of mammals has evolved to become an effective immune system. γδ T cells, a unique T cell subpopulation, are rare in circulation blood and lymphoid tissues, but rich in the intestinal mucosa, particularly in the epithelium. Via rapid production of cytokines and growth factors, intestinal γδ T cells are key contributors to epithelial homeostasis and immune surveillance of infection. Intriguingly, recent studies have revealed that the intestinal γδ T cells may play novel exciting functions ranging from epithelial plasticity and remodeling in response to carbohydrate diets to the recovery of ischemic stroke. In this review article, we update regulatory molecules newly defined in lymphopoiesis of the intestinal γδ T cells and their novel functions locally in the intestinal mucosa, such as epithelial remodeling, and distantly in pathological setting, e.g., ischemic brain injury repair, psychosocial stress responses, and fracture repair. The challenges and potential revenues in intestinal γδ T cell studies are discussed.

Integrative analysis of risk factors for immune-related adverse events of checkpoint blockade therapy in cancer

Immune-related adverse events (irAEs) induced by checkpoint inhibitors involve a multitude of different risk factors. Here, to interrogate the multifaceted underlying mechanisms, we compiled germline exomes and blood transcriptomes with clinical data, before and after checkpoint inhibitor treatment, from 672 patients with cancer. Overall, irAE samples showed a substantially lower contribution of neutrophils in terms of baseline and on-therapy cell counts and gene expression markers related to neutrophil function. Allelic variation of HLA-B correlated with overall irAE risk. Analysis of germline coding variants identified a nonsense mutation in an immunoglobulin superfamily protein, TMEM162. In our cohort and the Cancer Genome Atlas (TCGA) data, TMEM162 alteration was associated with higher peripheral and tumor-infiltrating B cell counts and suppression of regulatory T cells in response to therapy. We developed machine learning models for irAE prediction, validated using additional data from 169 patients. Our results provide valuable insights into risk factors of irAE and their clinical utility.

Genome-wide identification and skin expression of immunoglobulin superfamily in discus fish (Symphysodon aequifasciatus) reveal common genes associated with vertebrate lactation

Discus Symphysodon spp. employs an unusual parental care behavior where fry feed on parental skin mucus after hatching. Studies on discus immunoglobulin superfamily (IgSF) especially during parental care are scarce. Here, a total of 518 IgSF members were identified based on discus genome and clustered into 12 groups, unevenly distributing on 30 linkage groups. A total of 92 pairs of tandem duplication and 40 pairs of segmental duplication that underwent purifying selection were identified. IgSF genes expressed differentially in discus skin during different care stages and between male and female parents. Specifically, the transcription of btn1a1, similar with mammalian lactation, increased after spawning, reached a peak when fry started biting on parents' skin mucus, and then decreased. The expression of btn2a1 and other immune members, e.g., nect4, fcl5 and cd22, were up-regulated when fry stopped biting on mucus. These results suggest the expression differentiation of IgSF genes in skin of discus fish during parental care.

Division of labor and cooperation between different butyrophilin proteins controls phosphoantigen-mediated activation of human γδ T cells

Butyrophilin (BTN)-3A and BTN2A1 molecules control TCR-mediated activation of human Vγ9Vδ2 T-cells triggered by phosphoantigens (PAg) from microbes and tumors, but the molecular rules governing antigen sensing are unknown. Here we establish three mechanistic principles of PAg-action. Firstly, in humans, following PAg binding to the BTN3A1-B30.2 domain, Vγ9Vδ2 TCR triggering involves the V-domain of BTN3A2/BTN3A3. Moreover, PAg/B30.2 interaction, and the critical γδ-T-cell-activating V-domain, localize to different molecules. Secondly, this distinct topology as well as intracellular trafficking and conformation of BTN3A heteromers or ancestral-like BTN3A homomers are controlled by molecular interactions of the BTN3 juxtamembrane region. Finally, the ability of PAg not simply to bind BTN3A-B30.2, but to promote its subsequent interaction with the BTN2A1-B30.2 domain, is essential for T-cell activation. Defining these determinants of cooperation and division of labor in BTN proteins deepens understanding of PAg sensing and elucidates a mode of action potentially applicable to other BTN/BTNL family members.

Evolution of mitogen-activated protein kinase family and their immune function in Apostichopus japonicus

The mitogen-activated protein kinase family plays an important role in cell differentiation, growth, proliferation, and survival. However, the current research on the mitogen-activated protein kinase (MAPK) family in invertebrates is limited to the individual gene, and the analysis has not been conducted at the family level. In the present study, echinoderm MAPK family was identified by genomic screening, and five members, including three ERK subfamily members, one c-Jun N-terminal kinase (JNK) subfamily, and one p38-MAPK member were detected. Phylogenetic analysis showed that three MAPK subfamilies were separated into three separated clusters, and ERK subfamily appeared earlier than the other two subfamilies. Synteny analysis revealed that the p38 subfamily might be derived from the continuous gene duplication events of MAPK14 subfamily in invertebrates, which displayed genome expansion via gene duplication in vertebrates. The role of MAPK family in echinoderm immune defense was determined by investigating the expression profiles of MAPKs in Vibrio splendidus-challenged Apostichopus japonicus and LPS-exposed coelomocytes. The result showed that five MAPK members displayed induced expression profiles both in vitro and in vivo, and the peak expression was detected at different time points. Our study provides new insights into the evolutionary history of the MAPK family and show the similar immune function among MAPK members.

Predicting tumor immune microenvironment and checkpoint therapy response of head & neck cancer patients from blood immune single-cell transcriptomics

The immune state of tumor microenvironment is crucial for determining immunotherapy response but is not readily accessible. Here we investigate if we can infer the tumor immune state from the blood and further predict immunotherapy response. First, we analyze a dataset of head and neck squamous cell carcinoma (HNSCC) patients with matched scRNA-Seq of peripheral blood mononuclear cells (PBMCs) and tumor tissues. We find that the tumor immune cell fractions of different immune cell types and many of the genes they express can be inferred from the matched PBMC scRNA-Seq. Second, analyzing another HNSCC dataset with PBMC scRNA-Seq and immunotherapy response, we find that the inferred ratio between tumor memory B and regulatory T cell fractions is predictive of immunotherapy response and is superior to the well-established cytolytic and exhausted T-cell signatures. Overall, these results showcase the potential of scRNA-Seq liquid biopsies in cancer immunotherapy, calling for their larger-scale testing.

SIGNIFICANCE

This head and neck cancer study demonstrates the potential of using blood single-cell transcriptomics to (1) infer the tumor immune status and (2) predict immunotherapy response from the tumor immune status inferred from blood. These results showcase the potential of single-cell transcriptomics liquid biopsies for further advancing personalized cancer immunotherapy.

A stop-gain variant in BTNL9 is associated with atherogenic lipid profiles

Current understanding of lipid genetics has come mainly from studies in European-ancestry populations; limited effort has focused on Polynesian populations, whose unique population history and high prevalence of dyslipidemia may provide insight into the biological foundations of variation in lipid levels. Here, we performed an association study to fine map a suggestive association on 5q35 with high-density lipoprotein cholesterol (HDL-C) seen in Micronesian and Polynesian populations. Fine-mapping analyses in a cohort of 2,851 Samoan adults highlighted an association between a stop-gain variant (rs200884524; c.652C>T, p.R218∗; posterior probability = 0.9987) in BTNL9 and both lower HDL-C and greater triglycerides (TGs). Meta-analysis across this and several other cohorts of Polynesian ancestry from Samoa, American Samoa, and Aotearoa New Zealand confirmed the presence of this association (βHDL-C = -1.60 mg/dL, p HDL-C = 7.63 × 10-10; βTG = 12.00 mg/dL, p TG = 3.82 × 10-7). While this variant appears to be Polynesian specific, there is also evidence of association from other multiancestry analyses in this region. This work provides evidence of a previously unexplored contributor to the genetic architecture of lipid levels and underscores the importance of genetic analyses in understudied populations.

Genomic diversity and relationship analyses of endangered German Black Pied cattle (DSN) to 68 other taurine breeds based on whole-genome sequencing

German Black Pied cattle (Deutsches Schwarzbuntes Niederungsrind, DSN) are an endangered dual-purpose cattle breed originating from the North Sea region. The population comprises about 2,500 cattle and is considered one of the ancestral populations of the modern Holstein breed. The current study aimed at defining the breeds closest related to DSN cattle, characterizing their genomic diversity and inbreeding. In addition, the detection of selection signatures between DSN and Holstein was a goal. Relationship analyses using fixation index (F_ST), phylogenetic, and admixture analyses were performed between DSN and 68 other breeds from the 1000 Bull Genomes Project. Nucleotide diversity, observed heterozygosity, and expected heterozygosity were calculated as metrics for genomic diversity. Inbreeding was measured as excess of homozygosity (F_Hom) and genomic inbreeding (F_RoH) through runs of homozygosity (RoHs). Region-wide F_ST and cross-population-extended haplotype homozygosity (XP-EHH) between DSN and Holstein were used to detect selection signatures between the two breeds, and RoH islands were used to detect selection signatures within DSN and Holstein. DSN showed a close genetic relationship with breeds from the Netherlands, Belgium, Northern Germany, and Scandinavia, such as Dutch Friesian Red, Dutch Improved Red, Belgian Red White Campine, Red White Dual Purpose, Modern Angler, Modern Danish Red, and Holstein. The nucleotide diversity in DSN (0.151%) was higher than in Holstein (0.147%) and other breeds, e.g., Norwegian Red (0.149%), Red White Dual Purpose (0.149%), Swedish Red (0.149%), Hereford (0.145%), Angus (0.143%), and Jersey (0.136%). The F_Hom and F_RoH values in DSN were among the lowest. Regions with high F_ST between DSN and Holstein, significant XP-EHH regions, and RoH islands detected in both breeds harbor candidate genes that were previously reported for milk, meat, fertility, production, and health traits, including one QTL detected in DSN for endoparasite infection resistance. The selection signatures between DSN and Holstein provide evidence of regions responsible for the dual-purpose properties of DSN and the milk type of Holstein. Despite the small population size, DSN has a high level of diversity and low inbreeding. F_ST supports its relatedness to breeds from the same geographic origin and provides information on potential gene pools that could be used to maintain diversity in DSN.

Identification of tumor-specific neoantigens and immune clusters of hepatocellular carcinoma for mRNA vaccine development

BACKGROUND

To screen efficacious neoantigens for the development of LIHC mRNA vaccines, construct LIHC immune clusters, and therefore select patients who might benefit from vaccination.

METHODS

RNA-seq data and clinical information of 371 TCGA-LIHC and 231 ICGC-LIHC cohorts were downloaded. Differentially expressed genes and their associations with prognosis were analyzed by GEPIA, genetic alterations were examined in the cBioPortal portal, and the association between genes and immune infiltrating cells was explored by TIMER. The immune clusters were constructed by consistency clustering, and the immune landscape was described using CIBERSORT.

RESULTS

POLR3C and KPNA2 were identified as LIHC tumor neoantigens related to inferior prognosis and antigen-presenting cell infiltration. In addition, three immune clusters (IC1, IC2 and IC3) with significant differences in molecular, immune cytological, and clinical features were identified in both the TCGA and ICGC LIHC cohorts. Immune "hot" phenotype IC3 displayed a better survival than IC2, and immune "cold" phenotype IC1 exhibited a high tumor mutation burden.

CONCLUSION

In conclusion, for the development of anti-LIHC mRNA vaccines, we identified efficacious neoantigens POLR3C and KPNA2, profiled the tumor microenvironment of LIHC, and identified IC1 patients as the subgroup who might not most benefit from vaccination.

Genetics of varicose veins reveals polygenic architecture and genetic overlap with arterial and venous disease

Varicose veins represent a common cause of cardiovascular morbidity, with limited available medical therapies. Although varicose veins are heritable and epidemiologic studies have identified several candidate varicose vein risk factors, the molecular and genetic basis remains uncertain. Here we analyzed the contribution of common genetic variants to varicose veins using data from the Veterans Affairs Million Veteran Program and four other large biobanks. Among 49,765 individuals with varicose veins and 1,334,301 disease-free controls, we identified 139 risk loci. We identified genetic overlap between varicose veins, other vascular diseases and dozens of anthropometric factors. Using Mendelian randomization, we prioritized therapeutic targets via integration of proteomic and transcriptomic data. Finally, topological enrichment analyses confirmed the biologic roles of endothelial shear flow disruption, inflammation, vascular remodeling and angiogenesis. These findings may facilitate future efforts to develop nonsurgical therapies for varicose veins.

Construction of a metabolism-related gene prognostic model to predict survival of pancreatic cancer patients

Pancreatic cancer (PC) is one of the most fatal malignant tumors, and is commonly diagnosed at an advanced stage with no effective therapy. Metabolism-related genes (MRGs) and immune-related genes (IRGs) play considerable roles in the tumor microenvironment. Therefore, an effective prediction model based on MRGs and IRGs could aid in the prognosis of PC. In this study, differential expression analysis was performed to gain 25 intersectional genes from 857 differentially expressed MRGs (DEMRGs), and 1353 differentially expressed IRGs, from The Cancer Genome Atlas database of PC. Cox and Lasso regression were applied and a five-DEMRGs prognostic model constructed. Survival analysis, ROC values, risk curve and validation analysis showed that the model could independently predict PC prognosis. In addition, the correlation analysis suggested that the five-DEMRGs prognostic model could reflect the status of the immune microenvironment, including Tregs, M1 macrophages and Mast cell resting. Therefore, our study provides new underlying predictive biomarkers and associated immunotherapy targets.

Structure and evolution of the squamate major histocompatibility complex as revealed by two Anolis lizard genomes

The major histocompatibility complex (MHC) is an important genomic region for adaptive immunity and has long been studied in ecological and evolutionary contexts, such as disease resistance and mate and kin selection. The MHC has been investigated extensively in mammals and birds but far less so in squamate reptiles, the third major radiation of amniotes. We localized the core MHC genomic region in two squamate species, the green anole (Anolis carolinensis) and brown anole (A. sagrei), and provide the first detailed characterization of the squamate MHC, including the presence and ordering of known MHC genes in these species and comparative assessments of genomic structure and composition in MHC regions. We find that the Anolis MHC, located on chromosome 2 in both species, contains homologs of many previously-identified mammalian MHC genes in a single core MHC region. The repetitive element composition in anole MHC regions was similar to those observed in mammals but had important distinctions, such as higher proportions of DNA transposons. Moreover, longer introns and intergenic regions result in a much larger squamate MHC region (11.7 Mb and 24.6 Mb in the green and brown anole, respectively). Evolutionary analyses of MHC homologs of anoles and other representative amniotes uncovered generally monophyletic relationships between species-specific homologs and a loss of the peptide-binding domain exon 2 in one of two mhc2β gene homologs of each anole species. Signals of diversifying selection in each anole species was evident across codons of mhc1, many of which appear functionally relevant given known structures of this protein from the green anole, chicken, and human. Altogether, our investigation fills a major gap in understanding of amniote MHC diversity and evolution and provides an important foundation for future squamate-specific or vertebrate-wide investigations of the MHC.

BTN3A: A Promising Immune Checkpoint for Cancer Prognosis and Treatment

Butyrophilin-3A (BTN3A) subfamily members are a group of immunoglobulins present on the surface of different cell types, including innate and cancer cells. Due to their high similarity with the B7 family members, different studies have been conducted and revealed the involvement of BTN3A molecules in modulating T cell activity within the tumor microenvironment (TME). However, a great part of this research focused on γδ T cells and how BTN3A contributes to their functions. In this review, we will depict the roles and various aspects of BTN3A molecules in distinct tumor microenvironments and review how BTN3A receptors modulate diverse immune effector functions including those of CD4+ (Th1), cytotoxic CD8+ T cells, and NK cells. We will also highlight the potential of BTN3A molecules as therapeutic targets for effective immunotherapy and successful cancer control, which could represent a bright future for patient treatment.

A novel chromatin regulator-related immune checkpoint related gene prognostic signature and potential candidate drugs for endometrial cancer patients

Background

Endometrial cancer (EC) is the most common gynecologic malignancy in developed countries and its prevalence is increasing. As an emerging therapy with a promising efficacy, immunotherapy has been extensively applied in the treatment of solid tumors. In addition, chromatin regulators (CRs), as essential upstream regulators of epigenetics, play a significant role in tumorigenesis and cancer development.

Methods

CRs and immune checkpoint-related genes (ICRGs) were obtained from the previous top research. The Genome Cancer Atlas (TCGA) was utilized to acquire the mRNA expression and clinical information of patients with EC. Correlation analysis was utilized for screen CRs-related ICRGs (CRRICRGs). By Cox regression and least absolute shrinkage and selection operator (LASSO) analysis, prognosis related CRRICRGs were screened out and risk model was constructed. The Kaplan–Meier curve was used to estimate the prognosis between high- and low-risk group. By comparing the IC50 value, the drugs sensitivity difference was explored. We obtained small molecule drugs for the treatment of UCEC patients based on CAMP dataset.

Results

We successfully constructed a 9 CRRICRs-based prognostic signature for patients with UCEC and found the riskscore was an independent prognostic factor. The results of functional analysis suggested that CRRICRGs may be involved in immune processes associated with cancer. Immune characteristics analysis provided further evidence that the CRRICRGs-based model was correlated with immune cells infiltration and immune checkpoint. Eight small molecule drugs that may be effective for the treatment of UCEC patients were screened. Effective drugs identified by drug sensitivity profiling in high- and low-risk groups.

Conclusion

In summary, our study provided novel insights into the function of CRRICRGs in UCEC. We also developed a reliable prognostic panel for the survival of patients with UCEC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41065-022-00253-w.

BTN3A3 inhibits the proliferation, migration and invasion of ovarian cancer cells by regulating ERK1/2 phosphorylation

Butyrophilin Subfamily 3 Member A3 (BTN3A3) is a type I transmembrane protein belonging to the immunoglobulin (Ig) superfamily, which is expressed in many cancers. Clinical data show that ovarian cancer patients with high expression of BTN3A3 have a longer survival time, but the mechanism of BTN3A3 in the occurrence and progression of ovarian cancer is still unclear. Here, we found that BTN3A3 knockdown can promote the proliferation, migration and invasion of ovarian cancer cells, while overexpression of BTN3A3 can inhibit the proliferation, migration and invasion of ovarian cancer cells. We analyzed the immunoprecipitated BTN3A3 complex by mass spectrometry and found that BTN3A3 binds to FGF2, and the overexpression of BTN3A3 leads to a decrease in the protein level of FGF2, which in turn leads to a decrease in the level of phosphorylation of ERK1/2. By increasing the protein level of FGF2, it was found that the level of ERK1/2 phosphorylation also increased. Finally, the cancer promotion phenomenon caused by BTN3A3 knockdown can be improved by using ERK1/2 inhibitor SCH772984. To sum up, BTN3A3 interacts with FGF2, which inhibits FGF2/ERK1/2 axis and ultimately inhibits the proliferation, migration and invasion of ovarian cancer cells. Our results suggest that BTN3A3 may be a prognostic marker and a potential therapeutic target for ovarian cancer.

Missense Variant rs28362680 in BTNL2 Reduces Risk of Coronary Heart Disease

Background

The pathological basis of coronary heart disease (CHD) is atherosclerosis. BTNL2 can inhibit the activation of T cells. We aimed to explore the association between BTNL2 genetic variants and CHD risk in the southern Chinese Han population.

Methods

We recruited 1419 participants to perform an association analysis between missense variants in BTNL2 and CHD risk through SNPStats online software. Genotyping of all candidate SNPs were completed by the Agena MassARRAY. In addition, we used false-positive report probability analysis to detect whether the positive findings were noteworthy observations. We also used Haploview 4.2 software and SNPStats online software to conduct the haplotype analysis and analysis of linkage disequilibrium (LD). Finally, the interaction of SNP-SNP in CHD risk was evaluated by multi-factor dimensionality reduction (MDR).

Results

The results showed that BTNL2-rs35624343, -rs117896888, -rs41441651, -rs41417449, -rs28362680 and -rs2076523 were significantly associated with the CHD susceptibility. Especially for BTNL2-rs28362680, the allele A (OR = 0.68, p < 0.0001), genotype AA (OR = 0.40, p = 0.001) or GA (OR = 0.68, p < 0.0001) were associated with the reducing CHD risk. And -rs28362680 significantly reduced the CHD risk under all genetic models (dominant: OR = 0.64, p < 0.0001; recessive: OR = 0.47, p = 0.003; overdominant: OR = 0.73, p = 0.004; log-additive: OR = 0.66, p < 0.0001). And -rs28362680 was also closely associated with CHD risk reduction in all stratified analyses (age, gender, smoking, drinking, hypertension and diabetes). In addition, haplotype analysis showed that the “C_rs117896888C_rs41441651T_rs41417449A_rs28362680” (OR = 0.65, p < 0.0001) and “G_rs117896888T_rs41441651C_rs41417449A_rs28362680” (OR = 0.68, p = 0.013) may reduce CHD risk.

Conclusion

Missense variants (rs35624343, rs117896888, rs41441651, rs41417449, rs28362680, rs2076523) may be protective factors for the CHD risk. In particular, there were sufficient evidences that BTNL2-rs28362680 can protective CHD risk.

A novel immune checkpoint-related gene signature for hepatocellular carcinoma to predict clinical outcomes and therapeutic response

Immune checkpoint genes (ICGs) have recently been proven to perform instrumental functions in the maintenance of immune homeostasis and represent a promising therapeutic strategy; however, their expression patterns and prognostic values are not fully elucidated in hepatocellular carcinoma (HCC). In this investigation, we focused on establishing and validating a prognostic gene signature to facilitate decision-making in clinical practice. Clinical information, as well as transcriptome data, was obtained from the Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) database. Univariate Cox regression and least absolute shrinkage and selection operator (LASSO) Cox method were employed to build a multi-gene signature in the TCGA database, while the ICGC database was used for validation. Subsequently, utilizing the six-gene signature, we were able to categorize patients into high- and low-risk groups. In two cohorts, survival analysis findings revealed a dismal outlook for the high-risk group. The receiver operating characteristic curves were utilized to estimate the gene signature's prediction ability. Moreover, correlation analysis showed high-risk group was linked to advanced pathological stage, infiltration of immune cells and therapeutic response. In summary, this unique gene profile might serve not only as a useful prognostic indicator but also as a marker of therapy responsiveness in HCC.

Identifying a Novel Defined Pyroptosis-Associated Long Noncoding RNA Signature Contributes to Predicting Prognosis and Tumor Microenvironment of Bladder Cancer

Background

Bladder cancer (BLCA) is a common malignant tumor of the urinary tract, which is the sixth most common cancer among men. Numerous studies suggested that pyroptosis and long noncoding RNAs (lncRNAs) played an essential role in the development of cancers. However, the role of pyroptosis-related lncRNAs in BLCA and their prognostic value are still unclear.

Methods

In this study, we constructed a signature model through least absolute shrinkage and selection operator (LASSO) Cox regression analysis and Cox univariate analysis based on The Cancer Genome Atlas (TCGA) database. The expression of 12 pyroptosis-related lncRNAs was also confirmed by qRT-PCR in BLCA cell lines. TIMER, XCELL, QUANTISEQ, MCPCOUNTER, EPIC, and CIBERSORT R script were applied to quantify the relative proportions of infiltrating immune cells. Correlation coefficients were computed by Spearman analyses. The Kaplan–Meier method, Cox regression model, and log-rank tests were used to evaluate the prognostic value. The R package of pRRophetic was used to predict IC50 of common chemotherapeutic agents.

Results

A total of 12 pyroptosis-related lncRNAs with great prognosis value were identified. The expression was investigated by qRT-PCR in four BLCA cell lines. Then, 126 cases were identified as high-risk group, and 277 cases were identified as low-risk group based on the cutoff point. Patients in the low-risk group showed a significant survival advantage. Furthermore, we found that clinical features were significantly related to the risk score. As well, based on the C-index values, a nomogram was constructed. The gene set enrichment analysis (GSEA) results showed that mitogen-activated protein kinase (MAPK) signaling, transforming growth factor (TGF)-β signaling, and WNT signaling were with important significance in the high-risk group. Moreover, we found that riskscore was positively correlated with M0 macrophages and M2 macrophages.

Conclusions

In conclusion, our study indicated that pyroptosis is closely connected to BLCA. The riskscore generated from the expression of 12 pyroptosis-related lncRNAs was evaluated by various clinical features including survival status, tumor microenvironment, clinicopathological characteristic, and chemotherapy. It may offer a significant basis for future studies.

Recurrent expansions of B30.2-associated immune receptor families in fish

B30.2 domains, also known as PRY/SPRY, are key components of specific subsets of two large families of proteins involved in innate immunity: the tripartite motif proteins (TRIMs) and the Nod-like receptors (NLRs). TRIM proteins are important, often inducible factors of antiviral innate immunity, targeting multiple steps of viral cycles through a variety of mechanisms. NLRs prime and regulate systemic innate defenses, especially against bacteria, and control inflammation. Large TRIM and NLR subsets characterized by the presence of a B30.2 domain have been reported from a few fish species including zebrafish and seem to be strongly prone to gene duplication/expansion. Here, we performed a large-scale survey of these receptors across about 150 fish genomes, focusing on ray-finned fishes. We assessed the number and genomic distribution of domains and domain combinations associated with TRIMs, NLRs, and other genes containing B30.2 domains and looked for gene expansion patterns across fish groups. We then used a model to test the impact of taxonomy, genome size, and environmental variables on the copy numbers of these genes. Our findings reveal novel domain structures, clade-specific gains and losses. They also assist with the timing of the gene expansions, reveal patterns associated with the MHC, and lay the groundwork for further studies delving deeper into the forces that drive the copy number variation of immune genes on a species level.

Development of ICT01, a first-in-class, anti-BTN3A antibody for activating Vγ9Vδ2 T cell-mediated antitumor immune response

[Figure: see text].

BTN2A2 protein negatively regulates T cells to ameliorate collagen-induced arthritis in mice

Rheumatoid arthritis (RA) is an autoimmune disorder characterized by persistent inflammatory responses in target tissues and organs, resulting in the destruction of joints. Collagen type II (CII)-induced arthritis (CIA) is the most used animal model for human RA. Although BTN2A2 protein has been previously shown to inhibit T cell functions in vitro, its effect on autoimmune arthritis has not been reported. In this study, we investigate the ability of a recombinant BTN2A2-IgG2a Fc (BTN2A2-Ig) fusion protein to treat CIA. We show here that administration of BTN2A2-Ig attenuates established CIA, as compared with control Ig protein treatment. This is associated with reduced activation, proliferation and Th1/Th17 cytokine production of T cells in BTN2A2-Ig-treated CIA mice. BTN2A2-Ig also inhibits CII-specific T cell proliferation and Th1/Th17 cytokine production. Although the percentage of effector T cells is decreased in BTN2A2-Ig-treated CIA mice, the proportions of naive T cells and regulatory T cells is increased. Furthermore, BTN2A2-Ig reduces the percentage of proinflammatory M1 macrophages but increases the percentage of anti-inflammatory M2 macrophages in the CIA mice. Our results suggest that BTN2A2-Ig protein has the potential to be used in the treatment of collagen-induced arthritis models.

The extracellular region of bovine milk butyrophilin exhibits closer structural similarity to human myelin oligodendrocyte glycoprotein than to immunological BTN family receptors

Bovine butyrophilin (BTN1A1) is an abundant type I transmembrane glycoprotein exposed on the surface of milk fat globules. We have solved the crystal structure of its extracellular region via multiple wavelength anomalous dispersion after incorporation of selenomethionine into the bacterially produced protein. The butyrophilin ectodomain exhibits two subdomains with immunoglobulin fold, each comprising a β-sandwich with a central disulfide bridge as well as one N-linked glycosylation. The fifth Cys residue at position 193 is unpaired and prone to forming disulfide crosslinks. The apparent lack of a ligand-binding site or receptor activity suggests a function predominantly as hydrophilic coat protein to prevent coagulation of the milk fat droplets. While there is less structural resemblance to members of the human butyrophilin family such as BTN3A, which play a role as immune receptors, the N-terminal bovine butyrophilin subdomain shows surprising similarity to the human myelin oligodendrocyte glycoprotein, a protein exposed on the surface of myelin sheaths. Thus, our study lends structural support to earlier hypotheses of a correlation between the consumption of cow milk and prevalence of neurological autoimmune diseases and may offer guidance for the breeding of cattle strains that express modified butyrophilin showing less immunological cross-reactivity.

Multi-omics analysis of the expression and prognostic value of the butyrophilins in breast cancer

Butyrophilins (BTNs) belong to the immunoglobulin superfamily of transmembrane proteins and play a role in the regulation of lymphocyte activation, several autoimmune diseases, and the progression of human cancers. However, the associated clinicopathologic characteristics and prognostic value of BTNs in breast cancer remain unknown. This study aimed to discover potential key related BTN genes and signaling pathways in breast cancer, which could provide new insights for immune-based strategies. In the present study, the mRNA expression level and prognostic value of BTN2A1, BTN3A1, BTN3A2, BTN3A3, BTNL2, BTNL9, ERMAP, and MOG were measured. Up-regulation of these genes was significantly correlated with improved overall and relapse-free survival. We then analyzed the prognostic outcomes of breast cancer subtypes, genetic alterations, interaction networks, and the functional enrichment of eight BTN family genes. Our results showed that these eight genes played essential roles in tumor progression. Furthermore, an immune infiltration analysis indicated that most candidate BTN family members were associated with intratumoral immune cell infiltration, especially that of γδ T cells. Finally, gene set enrichment analysis for a single hub gene revealed that each BTN gene played a vital role in tumor progression through immune signaling pathways. These findings provided new insights into breast cancer pathogenesis and identified eight potential biomarkers for breast cancer.

RANKL-Induced Btn2a2 - A T Cell Immunomodulatory Molecule - During Osteoclast Differentiation Fine-Tunes Bone Resorption

Butyrophilins, which are members of the extended B7 family of immunoregulators structurally related to the B7 family, have diverse functions on immune cells as co-stimulatory and co-inhibitory molecules. Despite recent advances in the understanding on butyrophilins' role on adaptive immune cells during infectious or autoimmune diseases, nothing is known about their role in bone homeostasis. Here, we analyzed the role of one specific butyrophilin, namely Btn2a2, as we have recently shown that Btn2a2 is expressed on the monocyte/macrophage lineage that also gives rise to bone degrading osteoclasts. We found that expression of Btn2a2 on monocytes and pre-osteoclasts is upregulated by the receptor activator of nuclear factor κ-B ligand (RANKL), an essential protein required for osteoclast formation. Interestingly, in Btn2a2-deficient osteoclasts, typical osteoclast marker genes (Nfatc1, cathepsin K, TRAP, and RANK) were downregulated following RANKL stimulation. In vitro osteoclast assays resulted in decreased TRAP positive osteoclast numbers in Btn2a2-deficient cells. However, Btn2a2-deficient osteoclasts revealed abnormal fusion processes shown by their increased size. In vivo steady state µCT and histological analysis of bone architecture in complete Btn2a2-deficient mice showed differences in bone parameters further highlighting the fine-tuning effect of BTN2a2. Moreover, in rheumatoid arthritis patients and experimental arthritis, we detected significantly decreased serum levels of the secreted soluble Btn2a2 protein. Taken together, we identified the involvement of the immunomodulatory molecule Btn2a2 in osteoclast differentiation with potential future implications in basic and translational osteoimmunology.

ERMAP is a B7 family-related molecule that negatively regulates T cell and macrophage responses

T cell activation and tolerance are tightly regulated by costimulatory and coinhibitory molecules. B7 family members play a crucial role in regulating immune responses. In this study, we identified erythroid membrane-associated protein (ERMAP) as a novel T cell inhibitory molecule. ERMAP shares significant sequence and structural homology with existing B7 family members in its extracellular domain. The ERMAP protein is expressed on the cell surface of resting and activated antigen-presenting cells (APCs) and in some tumor tissues. The putative ERMAP receptor is expressed on activated CD4 and CD8 T cells and macrophages. Both mouse and human ERMAP-IgG2a Fc (ERMAP-Ig) fusion proteins inhibit T cell functions in vitro. Administration of ERMAP-Ig protein ameliorates autoimmune diseases, including experimental autoimmune encephalomyelitis and type 1 diabetes, in mice. Anti-ERMAP antibody enhances macrophage phagocytosis of cancer cells in vitro. Furthermore, administration of an anti-ERMAP antibody inhibits tumor growth in mice likely by blocking the inhibitory effects of ERMAP on T cells and macrophages. Our results suggest that therapeutic interaction with the ERMAP inhibitory pathway may represent a novel strategy for treating patients with autoimmune disease or cancer.

Nutritional and ecological perspectives of the interrelationships between diet and the gut microbiome in multiple sclerosis: Insights from marmosets

Studies in experimental autoimmune encephalomyelitis (EAE), the animal model of multiple sclerosis, have shown potential links between diet components, microbiome composition, and modulation of immune responses. In this review, we reanalyze and discuss findings in an outbred marmoset EAE model in which a yogurt-based dietary supplement decreased disease frequency and severity. We show that although diet has detectable effects on the fecal microbiome, microbiome changes are more strongly associated with the EAE development. Using an ecological framework, we further show that the dominant factors influencing the gut microbiota were marmoset sibling pair and experimental time point. These findings emphasize challenges in assigning cause-and-effect relationships in studies of diet-microbiome-host interactions and differentiating the diet effects from other environmental, stochastic, and host-related factors. We advocate for animal experiments to be designed to allow causal inferences of the microbiota's role in pathology while considering the complex ecological processes that shape microbial communities.

Prognostic and Therapeutic Significance of BTN3A Proteins in Tumors

The Butyrophilin 3A (BTN3A) family is a type I transmembrane protein belonging to the immunoglobulin (Ig) superfamily. The family contains three members: BTN3A1, BTN3A2 and BTN3A3, which share 95% homology in the extracellular domain. The expression of BTN3A family members is different in different types of tumors, which plays an important role in tumor prognosis. Among them, there are many studies on tumor immunity of BTN3A1, which shows that it is essential for the activation of Vγ9Vδ2 T cells, while BTN3A3 is expected to become a potential therapeutic target for breast cancer. Recent studies have shown that the BTN3A family is closely related to the occurrence and development of tumors. Now the BTN3A family has become one of the research hotspots and is expected to become new tumor prediction and treatment targets.

Low expression of BTN3A3 indicates poor prognosis and promotes cell proliferation, migration and invasion in non-small cell lung cancer

Background

The butyrophilin (BTN) family has many members with diverse functions related to immunomodulation, initiation and progression of tumors. BTN3A3 belongs to the BTN family, and exploring its expression and correlation with the prognosis of non-small cell lung cancer (NSCLC) patients has great clinical significance.

Methods

Clinical specimens were used to detect BTN3A3 expression. Small interfering RNA (siRNA) was used to knock down BTN3A3 and analyze the proliferative, migratory and invading ability of the transfected NSCLC cells. Multiplex immunofluorescence staining was used to detect the expression of BTN3A3 protein in the tumor microenvironment (TME). We analyzed the relationship between the expression of BTN3A3 and the clinicopathological features and prognosis of NSCLC patients.

Results

The expression of BTN3A3 in NSCLC tissues was significantly lower than in adjacent tissues, and patients with low expression of BTN3A3 had late clinical stages and lower degree of tumor differentiation. Knocking down BTN3A3 promoted the proliferation, migration and invasion of NSCLC. In the TME, the density of BTN3A3+ tumor cells positively correlated with the density of CD8+ T cells, and the patients with low expression of BTN3A3 had poor overall survival (OS).

Conclusions

Changes in the BTN3A3 expression level may play a potential key role in the carcinogenesis and development of NSCLC. Patients with low expression of BTN3A3 showed a more aggressive and invasive phenotype and a lower level of CD8+ T-cell infiltration, which may be an important factor affecting the OS of NSCLC patients.

Gut γδ T cells as guardians, disruptors, and instigators of cancer

Colorectal cancer is the third most common cancer worldwide with nearly 2 million cases per year. Immune cells and inflammation are a critical component of colorectal cancer progression, and they are used as reliable prognostic indicators of patient outcome. With the growing appreciation for immunology in colorectal cancer, interest is growing on the role γδ T cells have to play, as they represent one of the most prominent immune cell populations in gut tissue. This group of cells consists of both resident populations-γδ intraepithelial lymphocytes (γδ IELs)-and transient populations that each has unique functions. The homeostatic role of these γδ T cell subsets is to maintain barrier integrity and prevent microorganisms from breaching the mucosal layer, which is accomplished through crosstalk with enterocytes and other immune cells. Recent years have seen a surge in discoveries regarding the regulation of γδ IELs in the intestine and the colon with particular new insights into the butyrophilin family. In this review, we discuss the development, specialities, and functions of γδ T cell subsets during cancer progression. We discuss how these cells may be used to predict patient outcome, as well as how to exploit their behavior for cancer immunotherapy.

γδ T cells in tissue physiology and surveillance

γδ T cells are a unique T cell subpopulation that are rare in secondary lymphoid organs but enriched in many peripheral tissues, such as the skin, intestines and lungs. By rapidly producing large amounts of cytokines, γδ T cells make key contributions to immune responses in these tissues. In addition to their immune surveillance activities, recent reports have unravelled exciting new roles for γδ T cells in steady-state tissue physiology, with functions ranging from the regulation of thermogenesis in adipose tissue to the control of neuronal synaptic plasticity in the central nervous system. Here, we review the roles of γδ T cells in tissue homeostasis and in surveillance of infection, aiming to illustrate their major impact on tissue integrity, tissue repair and immune protection.

A glance over the fence: Using phylogeny and species comparison for a better understanding of antigen recognition by human γδ T-cells

Both, jawless and jawed vertebrates possess three lymphocyte lineages defined by highly diverse antigen receptors: Two T-cell- and one B-cell-like lineage. In both phylogenetic groups, the theoretically possible number of individual antigen receptor specificities can even outnumber that of lymphocytes of a whole organism. Despite fundamental differences in structure and genetics of these antigen receptors, convergent evolution led to functional similarities between the lineages. Jawed vertebrates possess αβ and γδ T-cells defined by eponymous αβ and γδ T-cell antigen receptors (TCRs). "Conventional" αβ T-cells recognize complexes of Major Histocompatibility Complex (MHC) class I and II molecules and peptides. Non-conventional T-cells, which can be αβ or γδ T-cells, recognize a large variety of ligands and differ strongly in phenotype and function between species and within an organism. This review describes similarities and differences of non-conventional T-cells of various species and discusses ligands and functions of their TCRs. A special focus is laid on Vγ9Vδ2 T-cells whose TCRs act as sensors for phosphorylated isoprenoid metabolites, so-called phosphoantigens (PAg), associated with microbial infections or altered host metabolism in cancer or after drug treatment. We discuss the role of butyrophilin (BTN)3A and BTN2A1 in PAg-sensing and how species comparison can help in a better understanding of this human Vγ9Vδ2 T-cell subset.

Harnessing Unconventional T Cells for Immunotherapy of Tuberculosis

Even if the incidence of tuberculosis (TB) has been decreasing over the last years, the number of patients with TB is increasing worldwide. The emergence of multidrug-resistant and extensively drug-resistant TB is making control of TB more difficult. Mycobacterium bovis bacillus Calmette–Guérin vaccine fails to prevent pulmonary TB in adults, and there is an urgent need for a vaccine that is also effective in patients with human immunodeficiency virus (HIV) coinfection. Therefore, TB control may benefit on novel therapeutic options beyond antimicrobial treatment. Host-directed immunotherapies could offer therapeutic strategies for patients with drug-resistant TB or with HIV and TB coinfection. In the last years, the use of donor lymphocytes after hematopoietic stem cell transplantation has emerged as a new strategy in the cure of hematologic malignancies in order to induce graft-versus leukemia and graft-versus-infection effects. Moreover, adoptive therapy has proven to be effective in controlling cytomegalovirus and Epstein-Barr virus reactivation in immunocompromised patients with ex vivo expanded viral antigen-specific T cells. Unconventional T cells are a heterogeneous group of T lymphocytes with limited diversity. One of their characteristics is that antigen recognition is not restricted by the classical major histocompatibility complex (MHC). They include CD1 (cluster of differentiation 1)–restricted T cells, MHC-related protein-1–restricted mucosal-associated invariant T (MAIT) cells, MHC class Ib–reactive T cells, and γδ T cells. Because these T cells are genotype-independent, they are also termed “donor unrestricted” T cells. The combined features of low donor diversity and the lack of genetic restriction make these cells suitable candidates for T cell–based immunotherapy of TB.

Expression profile of immune checkpoint genes and their roles in predicting immunotherapy response

Immune checkpoint genes (ICGs) play critical roles in circumventing self-reactivity and represent a novel target to develop treatments for cancers. However, a comprehensive analysis for the expression profile of ICGs at a pan-cancer level and their correlation with patient response to immune checkpoint blockade (ICB) based therapy is still lacking. In this study, we defined three expression patterns of ICGs using a comprehensive survey of RNA-seq data of tumor and immune cells from the functional annotation of the mammalian genome (FANTOM5) project. The correlation between the expression patterns of ICGs and patients survival and response to ICB therapy was investigated. The expression patterns of ICGs were robust across cancers, and upregulation of ICGs was positively correlated with high lymphocyte infiltration and good prognosis. Furthermore, we built a model (ICGe) to predict the response of patients to ICB therapy using five features of ICG expression. A validation scenario of six independent datasets containing data of 261 patients with CTLA-4 and PD-1 blockade immunotherapies demonstrated that ICGe achieved area under the curves of 0.64-0.82 and showed a robust performance and outperformed other mRNA-based predictors. In conclusion, this work revealed expression patterns of ICGs and underlying correlations between ICGs and response to ICB, which helps to understand the mechanisms of ICGs in ICB signal pathways and other anticancer treatments.

An Update on the Molecular Basis of Phosphoantigen Recognition by Vγ9Vδ2 T Cells

About 1-5% of human blood T cells are Vγ9Vδ2 T cells. Their hallmark is the expression of T cell antigen receptors (TCR) whose γ-chains contain a rearrangement of Vγ9 with JP (TRGV9JP or Vγ2Jγ1.2) and are paired with Vδ2 (TRDV2)-containing δ-chains. These TCRs respond to phosphoantigens (PAg) such as (E)-4-hydroxy-3-methyl-but-2-enyl pyrophosphate (HMBPP), which is found in many pathogens, and isopentenyl pyrophosphate (IPP), which accumulates in certain tumors or cells treated with aminobisphosphonates such as zoledronate. Until recently, these cells were believed to be restricted to primates, while no such cells are found in rodents. The identification of three genes pivotal for PAg recognition encoding for Vγ9, Vδ2, and butyrophilin (BTN) 3 in various non-primate species identified candidate species possessing PAg-reactive Vγ9Vδ2 T cells. Here, we review the current knowledge of the molecular basis of PAg recognition. This not only includes human Vγ9Vδ2 T cells and the recent discovery of BTN2A1 as Vγ9-binding protein mandatory for the PAg response but also insights gained from the identification of functional PAg-reactive Vγ9Vδ2 T cells and BTN3 in the alpaca and phylogenetic comparisons. Finally, we discuss models of the molecular basis of PAg recognition and implications for the development of transgenic mouse models for PAg-reactive Vγ9Vδ2 T cells.

Analysis of colon-infiltrating γδ T cells in chronic inflammatory bowel disease and in colitis-associated cancer

Inflammatory bowel disease (IBD) remains a global health problem with a significant percentage of patients progressing to chronic inflammation and colitis-associated cancer (CAC). Whether or not γδ T cells contribute to initiation and maintenance of inflammation in IBD and in the development of CAC is not known. We have evaluated the frequency, phenotype, and functions of γδ T cells among tissue-infiltrating lymphocytes in healthy donors and IBD and CAC patients. Results show that Vδ1 T cells are the dominant γδ T-cell population in healthy tissue, whereas Vδ2 T significantly abound in chronic IBD. Vδ2 T cells produce more IFN-γ, TNF-α, and IL-17 than Vδ1 T cells in chronic inflamed IBD. In CAC patients no significant cytokine production was detected in tissue-resident Vδ1 T cells, but Vδ2 T cells produced remarkable amounts of IFN-γ and TNF-α; these data were confirmed by the analysis of an independent cohort of IBD transcriptomes. Moreover, transcriptomes of IBD patients revealed a clear-cut clusterization of genes related with the maintenance of the inflammatory status. In conclusion, our results demonstrating that Vδ2 T cells have a proinflammatory profile in chronic IBD are suggestive of their participation in IBD and CAC pathogenesis.