Profiling immunoglobulin repertoires across multiple human tissues using RNA sequencing

Pathobiological signatures of dysbiotic lung injury in pediatric patients undergoing stem cell transplantation

Hematopoietic cell transplantation (HCT) uses cytotoxic chemotherapy and/or radiation followed by intravenous infusion of stem cells to cure malignancies, bone marrow failure and inborn errors of immunity, hemoglobin and metabolism. Lung injury is a known complication of the process, due in part to disruption in the pulmonary microenvironment by insults such as infection, alloreactive inflammation and cellular toxicity. How microorganisms, immunity and the respiratory epithelium interact to contribute to lung injury is uncertain, limiting the development of prevention and treatment strategies. Here we used 278 bronchoalveolar lavage (BAL) fluid samples to study the lung microenvironment in 229 pediatric patients who have undergone HCT treated at 32 children's hospitals between 2014 and 2022. By leveraging paired microbiome and human gene expression data, we identified high-risk BAL compositions associated with in-hospital mortality (P = 0.007). Disadvantageous profiles included bacterial overgrowth with neutrophilic inflammation, microbiome contraction with epithelial fibroproliferation and profound commensal depletion with viral and staphylococcal enrichment, lymphocytic activation and cellular injury, and were replicated in an independent cohort from the Netherlands (P = 0.022). In addition, a broad array of previously occult pathogens was identified, as well as a strong link between antibiotic exposure, commensal bacterial depletion and enrichment of viruses and fungi. Together these lung-immune system-microorganism interactions clarify the important drivers of fatal lung injury in pediatric patients who have undergone HCT. Further investigation is needed to determine how personalized interpretation of heterogeneous pulmonary microenvironments may be used to improve pediatric HCT outcomes.

Elucidating immunological characteristics of the adenoma-carcinoma sequence in colorectal cancer patients in South Korea using a bioinformatics approach

Colorectal cancer (CRC) is one of the top five most common and life-threatening malignancies worldwide. Most CRC develops from advanced colorectal adenoma (ACA), a precancerous stage, through the adenoma-carcinoma sequence. However, its underlying mechanisms, including how the tumor microenvironment changes, remain elusive. Therefore, we conducted an integrative analysis comparing RNA-seq data collected from 40 ACA patients who visited Dongguk University Ilsan Hospital with normal adjacent colons and tumor samples from 18 CRC patients collected from a public database. Differential expression analysis identified 21 and 79 sequentially up- or down-regulated genes across the continuum, respectively. The functional centrality of the continuum genes was assessed through network analysis, identifying 11 up- and 13 down-regulated hub-genes. Subsequently, we validated the prognostic effects of hub-genes using the Kaplan-Meier survival analysis. To estimate the immunological transition of the adenoma-carcinoma sequence, single-cell deconvolution and immune repertoire analyses were conducted. Significant composition changes for innate immunity cells and decreased plasma B-cells with immunoglobulin diversity were observed, along with distinctive immunoglobulin recombination patterns. Taken together, we believe our findings suggest underlying transcriptional and immunological changes during the adenoma-carcinoma sequence, contributing to the further development of pre-diagnostic markers for CRC.

Patterns of Unwanted Biological and Technical Expression Variation Among 49 Human Tissues

Tissue gene expression studies are impacted by biological and technical sources of variation, which can be broadly classified into wanted and unwanted variation. The latter, if not addressed, results in misleading biological conclusions. Methods have been proposed to reduce unwanted variation, such as normalization and batch correction. A more accurate understanding of all causes of variation could significantly improve the ability of these methods to remove unwanted variation while retaining variation corresponding to the biological question of interest. We used 17,282 samples from 49 human tissues in the Genotype-Tissue Expression data set (v8) to investigate patterns and causes of expression variation. Transcript expression was transformed to z-scores, and only the most variable 2% of transcripts were evaluated and clustered based on coexpression patterns. Clustered gene sets were assigned to different biological or technical causes based on histologic appearances and metadata elements. We identified 522 variable transcript clusters (median: 11 per tissue) among the samples. Of these, 63% were confidently explained, 16% were likely explained, 7% were low confidence explanations, and 14% had no clear cause. Histologic analysis annotated 46 clusters. Other common causes of variability included sex, sequencing contamination, immunoglobulin diversity, and compositional tissue differences. Less common biological causes included death interval (Hardy score), disease status, and age. Technical causes included blood draw timing and harvesting differences. Many of the causes of variation in bulk tissue expression were identifiable in the Tabula Sapiens data set of single-cell expression. This is among the largest explorations of the underlying sources of tissue expression variation. It uncovered expected and unexpected causes of variable gene expression and demonstrated the utility of matched histologic specimens. It further demonstrated the value of acquiring meaningful tissue harvesting metadata elements to use for improved normalization, batch correction, and analysis of both bulk and single-cell RNA-seq data.

Toolkit for mapping the clonal landscape of tumor-infiltrating B cells

Our current understanding of whether B cell involvement in the tumor microenvironment benefits the patient or the tumor - in distinct cancers, subcohorts and individual patients - is quite limited. Both statements are probably true in most cases: certain clonal B cell populations contribute to the antitumor response, while others steer the immune response away from the desired mechanics. To step up to a new level of understanding and managing B cell behaviors in the tumor microenvironment, we need to rationally discern these roles, which are cumulatively defined by B cell clonal functional programs, specificities of their B cell receptors, specificities and isotypes of the antibodies they produce, and their spatial interactions within the tumor environment. Comprehensive analysis of these characteristics of clonal B cell populations is now becoming feasible with the development of a whole arsenal of advanced technical approaches, which include (1) methods of single-cell and spatial transcriptomics, genomics, and proteomics; (2) methods of massive identification of B cell specificities; (3) methods of deep error-free profiling of B cell receptor repertoires. Here we overview existing techniques, summarize their current application for B cells studies and propose promising future directions in advancing B cells exploration.

Large-scale integrative analysis of juvenile idiopathic arthritis for new insight into its pathogenesis

BACKGROUND

Juvenile idiopathic arthritis (JIA) is one of the most prevalent rheumatic disorders in children and is classified as an autoimmune disease (AID). While a robust genetic contribution to JIA etiology has been established, the exact pathogenesis remains unclear.

METHODS

To prioritize biologically interpretable susceptibility genes and proteins for JIA, we conducted transcriptome-wide and proteome-wide association studies (TWAS/PWAS). Then, to understand the genetic architecture of JIA, we systematically analyzed single-nucleotide polymorphism (SNP)-based heritability, a signature of natural selection, and polygenicity. Next, we conducted HLA typing using multi-ethnicity RNA sequencing data. Additionally, we examined the T cell receptor (TCR) repertoire at a single-cell level to explore the potential links between immunity and JIA risk.

RESULTS

We have identified 19 TWAS genes and two PWAS proteins associated with JIA risks. Furthermore, we observe that the heritability and cell type enrichment analysis of JIA are enriched in T lymphocytes and HLA regions and that JIA shows higher polygenicity compared to other AIDs. In multi-ancestry HLA typing, B*45:01 is more prevalent in African JIA patients than in European JIA patients, whereas DQA1*01:01, DQA1*03:01, and DRB1*04:01 exhibit a higher frequency in European JIA patients. Using single-cell immune repertoire analysis, we identify clonally expanded T cell subpopulations in JIA patients, including CXCL13+BHLHE40+ TH cells which are significantly associated with JIA risks.

CONCLUSION

Our findings shed new light on the pathogenesis of JIA and provide a strong foundation for future mechanistic studies aimed at uncovering the molecular drivers of JIA.

Multiomics Data Reveal the Important Role of ANXA2R in T Cell-mediated Rejection After Renal Transplantation

BACKGROUND

T cell-mediated rejection (TCMR) is a severe issue after renal transplantation, but research on its T cell-receptor (TCR) repertoire is lacking. This study intended to elucidate the TCR repertoire landscape in TCMR and hence identify novel potential targets.

METHODS

A total of 12 multiomics data sets were collected. The TRUST4 algorithm was used to construct and analyze the TCR repertoire in renal allografts with TCMR and stable renal function. Then, novel TCR-related key genes were identified through various criteria and literature research. In bulk transcriptome, cell line, single-cell transcriptome data sets, multiple immune cell infiltration algorithms, and gene set enrichment analysis were used to analyze potential mechanisms of the identified key gene. Twenty-three pathological sections were collected for immunofluorescence staining in the clinical cohort. Finally, the diagnostic and prognostic values of ANXA2R were evaluated in multiple renal transplant data sets.

RESULTS

Allografts with TCMR showed significantly increased clonotype and specific clonal expansion. ANXA2R was found to be a novel key gene for TCMR and showed strong positive connections with the TCR complex and lymphocyte cells, especially CD8 + T cells. Immunofluorescence staining confirmed the existence of ANXA2R + CD8 + T cells, with their percentage significantly elevated in TCMR compared with stable renal function. Finally, both mRNA and protein levels of ANXA2R showed promising diagnostic and prognostic value for renal transplant recipients.

CONCLUSIONS

ANXA2R , identified as a novel TCR-related gene, had critical roles in clinicopathology, diagnosis, and prognosis in renal transplantation, which offered promising potential therapeutic targets.

Comprehensive characterization and database construction of immune repertoire in the largest Chinese glioma cohort

Immune receptor repertoire is valuable for developing immunotherapeutic interventions, but remains poorly understood across glioma subtypes including IDH wild type, IDH mutation without 1p/19q codeletion (IDHmut-noncodel) and IDH mutation with 1p/19q codeletion (IDHmut-codel). We assembled over 320,000 TCR/BCR clonotypes from the largest glioma cohort of 913 RNA sequencing samples in the Chinese population, finding that immune repertoire diversity was more prominent in the IDH wild type (the most aggressive glioma). Fewer clonotypes were shared within each glioma subtype, indicating high heterogeneity of the immune repertoire. The TRA-CDR3 was longer in private than in public clonotypes in IDH wild type. CDR3 variable motifs had higher proportions of hydrophobic residues in private than in public clonotypes, suggesting private CDR3 sequences have greater potential for tumor antigen recognition. Finally, we developed GTABdb, a web-based database designed for hosting, exploring, visualizing, and analyzing glioma immune repertoire. Our study will facilitate developing glioma immunotherapy.

Computational immunogenomic approaches to predict response to cancer immunotherapies

Cancer immunogenomics is an emerging field that bridges genomics and immunology. The establishment of large-scale genomic collaborative efforts along with the development of new single-cell transcriptomic techniques and multi-omics approaches have enabled characterization of the mutational and transcriptional profiles of many cancer types and helped to identify clinically actionable alterations as well as predictive and prognostic biomarkers. Researchers have developed computational approaches and machine learning algorithms to accurately obtain clinically useful information from genomic and transcriptomic sequencing data from bulk tissue or single cells and explore tumours and their microenvironment. The rapid growth in sequencing and computational approaches has resulted in the unmet need to understand their true potential and limitations in enabling improvements in the management of patients with cancer who are receiving immunotherapies. In this Review, we describe the computational approaches currently available to analyse bulk tissue and single-cell sequencing data from cancer, stromal and immune cells, as well as how best to select the most appropriate tool to address various clinical questions and, ultimately, improve patient outcomes.

Whole blood transcriptomics identifies subclasses of pediatric septic shock

BACKGROUND

Sepsis is a highly heterogeneous syndrome, which has hindered the development of effective therapies. This has prompted investigators to develop a precision medicine approach aimed at identifying biologically homogenous subgroups of patients with septic shock and critical illnesses. Transcriptomic analysis can identify subclasses derived from differences in underlying pathophysiological processes that may provide the basis for new targeted therapies. The goal of this study was to elucidate pathophysiological pathways and identify pediatric septic shock subclasses based on whole blood RNA expression profiles.

METHODS

The subjects were critically ill children with cardiopulmonary failure who were a part of a prospective randomized insulin titration trial to treat hyperglycemia. Genome-wide expression profiling was conducted using RNA sequencing from whole blood samples obtained from 46 children with septic shock and 52 mechanically ventilated noninfected controls without shock. Patients with septic shock were allocated to subclasses based on hierarchical clustering of gene expression profiles, and we then compared clinical characteristics, plasma inflammatory markers, cell compositions using GEDIT, and immune repertoires using Imrep between the two subclasses.

RESULTS

Patients with septic shock depicted alterations in innate and adaptive immune pathways. Among patients with septic shock, we identified two subtypes based on gene expression patterns. Compared with Subclass 2, Subclass 1 was characterized by upregulation of innate immunity pathways and downregulation of adaptive immunity pathways. Subclass 1 had significantly worse clinical outcomes despite the two classes having similar illness severity on initial clinical presentation. Subclass 1 had elevated levels of plasma inflammatory cytokines and endothelial injury biomarkers and demonstrated decreased percentages of CD4 T cells and B cells and less diverse T cell receptor repertoires.

CONCLUSIONS

Two subclasses of pediatric septic shock patients were discovered through genome-wide expression profiling based on whole blood RNA sequencing with major biological and clinical differences. Trial Registration This is a secondary analysis of data generated as part of the observational CAF-PINT ancillary of the HALF-PINT study (NCT01565941). Registered March 29, 2012.

Human IgM-expressing memory B cells

A hallmark of T cell dependent (TD) humoral immune responses is the generation of long-lived memory B cells. The generation of these cells occurs primarily in the germinal center (GC) reaction, where antigen-activated B cells undergo affinity maturation as a major consequence of the combined processes of proliferation, somatic hypermutation of their immunoglobulin V (IgV) region genes, and selection for improved affinity of their B-cell antigen receptors. As many B cells also undergo class-switching to IgG or IgA in these TD responses, there was traditionally a focus on class-switched memory B cells in both murine and human studies on memory B cells. However, it has become clear that there is also a large subset of IgM-expressing memory B cells, which have important phenotypic and functional similarities but also differences to class-switched memory B cells. There is an ongoing discussion about the origin of distinct subsets of human IgM+ B cells with somatically mutated IgV genes. We argue here that the vast majority of human IgM-expressing B cells with somatically mutated IgV genes in adults is indeed derived from GC reactions, even though a generation of some mostly lowly mutated IgM+ B cells from other differentiation pathways, mainly in early life, may exist.

Pulmonary microbiome and transcriptome signatures reveal distinct pathobiologic states associated with mortality in two cohorts of pediatric stem cell transplant patients

Lung injury is a major determinant of survival after pediatric hematopoietic cell transplantation (HCT). A deeper understanding of the relationship between pulmonary microbes, immunity, and the lung epithelium is needed to improve outcomes. In this multicenter study, we collected 278 bronchoalveolar lavage (BAL) samples from 229 patients treated at 32 children's hospitals between 2014-2022. Using paired metatranscriptomes and human gene expression data, we identified 4 patient clusters with varying BAL composition. Among those requiring respiratory support prior to sampling, in-hospital mortality varied from 22-60% depending on the cluster (p=0.007). The most common patient subtype, Cluster 1, showed a moderate quantity and high diversity of commensal microbes with robust metabolic activity, low rates of infection, gene expression indicating alveolar macrophage predominance, and low mortality. The second most common cluster showed a very high burden of airway microbes, gene expression enriched for neutrophil signaling, frequent bacterial infections, and moderate mortality. Cluster 3 showed significant depletion of commensal microbes, a loss of biodiversity, gene expression indicative of fibroproliferative pathways, increased viral and fungal pathogens, and high mortality. Finally, Cluster 4 showed profound microbiome depletion with enrichment of Staphylococci and viruses, gene expression driven by lymphocyte activation and cellular injury, and the highest mortality. BAL clusters were modeled with a random forest classifier and reproduced in a geographically distinct validation cohort of 57 patients from The Netherlands, recapitulating similar cluster-based mortality differences (p=0.022). Degree of antibiotic exposure was strongly associated with depletion of BAL microbes and enrichment of fungi. Potential pathogens were parsed from all detected microbes by analyzing each BAL microbe relative to the overall microbiome composition, which yielded increased sensitivity for numerous previously occult pathogens. These findings support personalized interpretation of the pulmonary microenvironment in pediatric HCT, which may facilitate biology-targeted interventions to improve outcomes.

Comment on 'rigorous benchmarking of T cell receptor repertoire profiling methods for cancer RNA sequencing'

Transcriptome sequencing has become common in cancer research, resulting in the generation of a substantial volume of RNA sequencing (RNA-Seq) data. The ability to extract immune repertoires from these data is crucial for obtaining information on infiltrating T- and B-lymphocyte clones when dedicated amplicon T-cell/B-cell receptors sequencing (TCR-Seq/BCR-Seq) methods are unavailable. In response to this demand, several dedicated computational methods have been developed, including MiXCR, TRUST and ImRep. In the recent publication in Briefings in Bioinformatics, Peng et al. have conducted an intensive, systematic comparison of the three previously mentioned tools. Although their effort is commendable, we do have a few constructive critiques regarding technical elements of their analysis.

Whole Blood Transcriptomics Identifies Subclasses of Pediatric Septic Shock

Background

Sepsis is a highly heterogeneous syndrome, that has hindered the development of effective therapies. This has prompted investigators to develop a precision medicine approach aimed at identifying biologically homogenous subgroups of patients with septic shock and critical illnesses. Transcriptomic analysis can identify subclasses derived from differences in underlying pathophysiological processes that may provide the basis for new targeted therapies. The goal of this study was to elucidate pathophysiological pathways and identify pediatric septic shock subclasses based on whole blood RNA expression profiles.

Methods

The subjects were critically ill children with cardiopulmonary failure who were a part of a prospective randomized insulin titration trial to treat hyperglycemia. Genome-wide expression profiling was conducted using RNA-sequencing from whole blood samples obtained from 46 children with septic shock and 52 mechanically ventilated noninfected controls without shock. Patients with septic shock were allocated to subclasses based on hierarchical clustering of gene expression profiles, and we then compared clinical characteristics, plasma inflammatory markers, cell compositions using GEDIT, and immune repertoires using Imrep between the two subclasses.

Results

Patients with septic shock depicted alterations in innate and adaptive immune pathways. Among patients with septic shock, we identified two subtypes based on gene expression patterns. Compared with Subclass 2, Subclass 1 was characterized by upregulation of innate immunity pathways and downregulation of adaptive immunity pathways. Subclass 1 had significantly worse clinical outcomes despite the two classes having similar illness severity on initial clinical presentation. Subclass 1 had elevated levels of plasma inflammatory cytokines and endothelial injury biomarkers and demonstrated decreased percentages of CD4 T cells and B cells, and less diverse T-Cell receptor repertoires.

Conclusions

Two subclasses of pediatric septic shock patients were discovered through genome-wide expression profiling based on whole blood RNA sequencing with major biological and clinical differences.

Trial Registration

This is a secondary analysis of data generated as part of the observational CAF PINT ancillary of the HALF PINT study (NCT01565941). Registered 29 March 2012.

Overexpression of CD73 is associated with recurrence and poor prognosis of gingivobuccal oral cancer as revealed by transcriptome and deep immune profiling of paired tumor and margin tissues

BACKGROUND

For various cancers, differences in response to treatment and subsequent survival period have been reported to be associated with variation in immune contextures.

AIM

We sought to identify whether such association exists in respect of gingivobuccal oral cancer.

MATERIALS AND METHODS

We performed deep immune profiling of tumor and margin tissues collected from 46 treatment naïve, Human Papillomavirus (HPV) negative, patients. Each patient was followed for 24 months and prognosis (recurrence/death) noted. Key findings were validated by comparing with TCGA-HNSC cohort data.

RESULTS

About 28% of patients showed poor post-treatment prognosis. These patients exhibited a high probability of recurrence even within 1 year and death within 2 years. There was restricted immune cell infiltration in tumor, but not in margin, among these patients. Reduced expression of eight immune-related genes (IRGs) (NT5E, THRA, RBP1, TLR4, ITGA6, BMPR1B, ITGAV, SSTR1) in tumor strongly predicted better quality of prognosis, both in our patient cohort and in TCGA-HNSC cohort. Tumors of patients with better prognosis were associated with (a) lower CD73+ cells with concomitant lower expression level of NT5E/CD73, (b) higher proportions of CD4+ and CD8+ T cells, B cells, NK cells, M1 macrophages, (c) higher %Granzyme+ cells, (d) higher TCR and BCR repertoire diversities. CD73 expression in tumor was associated with low CD8+ and CD4+ T cells, low immune repertoire diversity, and advanced cancer stage.

DISCUSSION AND CONCLUSION

High infiltration of anti-tumor immune cells in both tumors and margins results in good prognosis, while in patients with minimal infiltration in tumors in spite of high infiltration in margins results in poor prognosis. Targeted CD73 immune-checkpoint inhibition may improve clinical outcome.

Rigorous benchmarking of T-cell receptor repertoire profiling methods for cancer RNA sequencing

The ability to identify and track T-cell receptor (TCR) sequences from patient samples is becoming central to the field of cancer research and immunotherapy. Tracking genetically engineered T cells expressing TCRs that target specific tumor antigens is important to determine the persistence of these cells and quantify tumor responses. The available high-throughput method to profile TCR repertoires is generally referred to as TCR sequencing (TCR-Seq). However, the available TCR-Seq data are limited compared with RNA sequencing (RNA-Seq). In this paper, we have benchmarked the ability of RNA-Seq-based methods to profile TCR repertoires by examining 19 bulk RNA-Seq samples across 4 cancer cohorts including both T-cell-rich and T-cell-poor tissue types. We have performed a comprehensive evaluation of the existing RNA-Seq-based repertoire profiling methods using targeted TCR-Seq as the gold standard. We also highlighted scenarios under which the RNA-Seq approach is suitable and can provide comparable accuracy to the TCR-Seq approach. Our results show that RNA-Seq-based methods are able to effectively capture the clonotypes and estimate the diversity of TCR repertoires, as well as provide relative frequencies of clonotypes in T-cell-rich tissues and low-diversity repertoires. However, RNA-Seq-based TCR profiling methods have limited power in T-cell-poor tissues, especially in highly diverse repertoires of T-cell-poor tissues. The results of our benchmarking provide an additional appealing argument to incorporate RNA-Seq into the immune repertoire screening of cancer patients as it offers broader knowledge into the transcriptomic changes that exceed the limited information provided by TCR-Seq.

Large-scale antibody immune response mapping of splenic B cells and bone marrow plasma cells in a transgenic mouse model

Molecular characterization of antibody immunity and human antibody discovery is mainly carried out using peripheral memory B cells, and occasionally plasmablasts, that express B cell receptors (BCRs) on their cell surface. Despite the importance of plasma cells (PCs) as the dominant source of circulating antibodies in serum, PCs are rarely utilized because they do not express surface BCRs and cannot be analyzed using antigen-based fluorescence-activated cell sorting. Here, we studied the antibodies encoded by the entire mature B cell populations, including PCs, and compared the antibody repertoires of bone marrow and spleen compartments elicited by immunization in a human immunoglobulin transgenic mouse strain. To circumvent prior technical limitations for analysis of plasma cells, we applied single-cell antibody heavy and light chain gene capture from the entire mature B cell repertoires followed by yeast display functional analysis using a cytokine as a model immunogen. We performed affinity-based sorting of antibody yeast display libraries and large-scale next-generation sequencing analyses to follow antibody lineage performance, with experimental validation of 76 monoclonal antibodies against the cytokine antigen that identified three antibodies with exquisite double-digit picomolar binding affinity. We observed that spleen B cell populations generated higher affinity antibodies compared to bone marrow PCs and that antigen-specific splenic B cells had higher average levels of somatic hypermutation. A degree of clonal overlap was also observed between bone marrow and spleen antibody repertoires, indicating common origins of certain clones across lymphoid compartments. These data demonstrate a new capacity to functionally analyze antigen-specific B cell populations of different lymphoid organs, including PCs, for high-affinity antibody discovery and detailed fundamental studies of antibody immunity.

Antigen-specificity measurements are the key to understanding T cell responses

Antigen-specific T cells play a central role in the adaptive immune response and come in a wide range of phenotypes. T cell receptors (TCRs) mediate the antigen-specificities found in T cells. Importantly, high-throughput TCR sequencing provides a fingerprint which allows tracking of specific T cells and their clonal expansion in response to particular antigens. As a result, many studies have leveraged TCR sequencing in an attempt to elucidate the role of antigen-specific T cells in various contexts. Here, we discuss the published approaches to studying antigen-specific T cells and their specific TCR repertoire. Further, we discuss how these methods have been applied to study the TCR repertoire in various diseases in order to characterize the antigen-specific T cells involved in the immune control of disease.

Age-related changes of the human splenic marginal zone B cell compartment

In this review, we will summarize the growing body of knowledge on the age-related changes of human splenic B cell composition and molecular evidence of immune maturation and discuss the contribution of these changes on splenic protective function. From birth on, the splenic marginal zone (sMZ) contains a specialized B cell subpopulation, which recruits and archives memory B cells from immune responses throughout the organism. The quality of sMZ B cell responses is augmented by germinal center (GC)-dependent maturation of memory B cells during childhood, however, in old age, these mechanisms likely contribute to waning of splenic protective function.

Łabaj P, Mangul S. RNA-seq data science: From raw data to effective interpretation

RNA sequencing (RNA-seq) has become an exemplary technology in modern biology and clinical science. Its immense popularity is due in large part to the continuous efforts of the bioinformatics community to develop accurate and scalable computational tools to analyze the enormous amounts of transcriptomic data that it produces. RNA-seq analysis enables genes and their corresponding transcripts to be probed for a variety of purposes, such as detecting novel exons or whole transcripts, assessing expression of genes and alternative transcripts, and studying alternative splicing structure. It can be a challenge, however, to obtain meaningful biological signals from raw RNA-seq data because of the enormous scale of the data as well as the inherent limitations of different sequencing technologies, such as amplification bias or biases of library preparation. The need to overcome these technical challenges has pushed the rapid development of novel computational tools, which have evolved and diversified in accordance with technological advancements, leading to the current myriad of RNA-seq tools. These tools, combined with the diverse computational skill sets of biomedical researchers, help to unlock the full potential of RNA-seq. The purpose of this review is to explain basic concepts in the computational analysis of RNA-seq data and define discipline-specific jargon.

Patterns of unwanted biological and technical expression variation across 49 human tissues

All tissue-based gene expression studies are impacted by biological and technical sources of variation. Numerous methods are used to normalize and batch correct these datasets. A more accurate understanding of all causes of variation could further optimize these approaches. We used 17,282 samples from 49 tissues in the Genotype Tissue Expression (GTEx) dataset (v8) to investigate patterns and causes of expression variation. Transcript expression was normalized to Z-scores and only the most variable 2% of transcripts were evaluated and clustered based on co-expression patterns. Clustered gene sets were solved to different biological or technical causes related to metadata elements and histologic images. We identified 522 variable transcript clusters (median 11 per tissue) across the samples. Of these, 64% were confidently explained, 15% were likely explained, 7% were low confidence explanations and 14% had no clear cause. Common causes included sex, sequencing contamination, immunoglobulin diversity, and compositional tissue differences. Less common biological causes included death interval (Hardy score), muscle atrophy, diabetes status, and menopause. Technical causes included brain pH and harvesting differences. Many of the causes of variation in bulk tissue expression were identifiable in the Tabula Sapiens dataset of single cell expression. This is the largest exploration of the underlying sources of tissue expression variation. It uncovered expected and unexpected causes of variable gene expression. These identified sources of variation will inform which metadata to acquire with tissue harvesting and can be used to improve normalization, batch correction, and analysis of both bulk and single cell RNA-seq data.

Tumour infiltrating B cells discriminate checkpoint blockade-induced responses

BACKGROUND

Immune cell-driven anti-cancer activity is paramount for effective responses to checkpoint inhibitors (ICB). However, the contribution of the different immune cell subsets in the circulation and within the tumour is poorly understood.

MATERIALS AND METHODS

To elucidate the role of the different cell subsets in anti-tumour responses elicited by ICB, we performed single-cell analysis of the transcriptome and surface proteome of paired pre- and early on-treatment metastatic melanoma tumour biopsies and matched peripheral blood mononuclear cell samples. We next compared the survival of metastatic melanoma patients treated with ICB according to the abundance of pre-treatment tumour-infiltrating B cell clonotypes.

RESULTS

We identified cell clusters associated with disease control or progression, defined differential expression of biological pathways likely involved in the immune awakening against the tumour and examined how cell-cell communication patterns between the tumour cell subsets change during treatment. Furthermore, we discovered that B cells (immunoglobulin expression and abundance of B cell clonotypes) discriminate the clinical response after ICB and propose that B cells likely contribute to anti-tumour immunity by antigen presentation through major histocompatibility complex molecules. Finally, we demonstrated that the abundance of tumour-infiltrating B cell clonotypes at baseline identifies two distinct risk groups, a finding that we confirmed in an independent cohort.

CONCLUSIONS

Our exploratory translational study provides new insights on the mechanistic role of B cells in anti-melanoma immunity during treatment with ICB. Additionally, we support pre-treatment B cell tumour infiltration as a promising prognostic biomarker to be further validated as a tool for clinical risk stratification.

Vaccination History, Body Mass Index, Age, and Baseline Gene Expression Predict Influenza Vaccination Outcomes

Seasonal influenza is a primary public health burden in the USA and globally. Annual vaccination programs are designed on the basis of circulating influenza viral strains. However, the effectiveness of the seasonal influenza vaccine is highly variable between seasons and among individuals. A number of factors are known to influence vaccination effectiveness including age, sex, and comorbidities. Here, we sought to determine whether whole blood gene expression profiling prior to vaccination is informative about pre-existing immunological status and the immunological response to vaccine. We performed whole transcriptome analysis using RNA sequencing (RNAseq) of whole blood samples obtained prior to vaccination from 275 participants enrolled in an annual influenza vaccine trial. Serological status prior to vaccination and 28 days following vaccination was assessed using the hemagglutination inhibition assay (HAI) to define baseline immune status and the response to vaccination. We find evidence that genes with immunological functions are increased in expression in individuals with higher pre-existing immunity and in those individuals who mount a greater response to vaccination. Using a random forest model, we find that this set of genes can be used to predict vaccine response with a performance similar to a model that incorporates physiological and prior vaccination status alone. A model using both gene expression and physiological factors has the greatest predictive power demonstrating the potential utility of molecular profiling for enhancing prediction of vaccine response. Moreover, expression of genes that are associated with enhanced vaccination response may point to additional biological pathways that contribute to mounting a robust immunological response to the seasonal influenza vaccine.

Changes in the Immune Cell Repertoire for the Treatment of Malignant Melanoma

Immune checkpoint inhibitors (ICIs) have been used for the treatment of various types of cancers, including malignant melanoma. Mechanistic exploration of tumor immune responses is essential to improve the therapeutic efficacy of ICIs. Since tumor immune responses are based on antigen-specific immune responses, investigators have focused on T cell receptors (TCRs) and have analyzed changes in the TCR repertoire. The proliferation of T cell clones against tumor antigens is detected in patients who respond to treatment with ICIs. The proliferation of these T cell clones is observed within tumors as well as in the peripheral blood. Clonal proliferation has been detected not only in CD8-positive T cells but also in CD4-positive T cells, resident memory T cells, and B cells. Moreover, changes in the repertoire at an early stage of treatment seem to be useful for predicting the therapeutic efficacy of ICIs. Further analyses of the repertoire of immune cells are desirable to improve and predict the therapeutic efficacy of ICIs.

Memory persistence and differentiation into antibody-secreting cells accompanied by positive selection in longitudinal BCR repertoires

The stability and plasticity of B cell-mediated immune memory ensures the ability to respond to the repeated challenges. We have analyzed the longitudinal dynamics of immunoglobulin heavy chain repertoires from memory B cells, plasmablasts, and plasma cells from the peripheral blood of generally healthy volunteers. We reveal a high degree of clonal persistence in individual memory B cell subsets, with inter-individual convergence in memory and antibody-secreting cells (ASCs). ASC clonotypes demonstrate clonal relatedness to memory B cells, and are transient in peripheral blood. We identify two clusters of expanded clonal lineages with differing prevalence of memory B cells, isotypes, and persistence. Phylogenetic analysis revealed signs of reactivation of persisting memory B cell-enriched clonal lineages, accompanied by new rounds of affinity maturation during proliferation and differentiation into ASCs. Negative selection contributes to both persisting and reactivated lineages, preserving the functionality and specificity of B cell receptors (BCRs) to protect against current and future pathogens.

Benchmarking computational methods for B-cell receptor reconstruction from single-cell RNA-seq data

Multiple methods have recently been developed to reconstruct full-length B-cell receptors (BCRs) from single-cell RNA sequencing (scRNA-seq) data. This need emerged from the expansion of scRNA-seq techniques, the increasing interest in antibody-based drug development and the importance of BCR repertoire changes in cancer and autoimmune disease progression. However, a comprehensive assessment of performance-influencing factors such as the sequencing depth, read length or number of somatic hypermutations (SHMs) as well as guidance regarding the choice of methodology is still lacking. In this work, we evaluated the ability of six available methods to reconstruct full-length BCRs using one simulated and three experimental SMART-seq datasets. In addition, we validated that the BCRs assembled in silico recognize their intended targets when expressed as monoclonal antibodies. We observed that methods such as BALDR, BASIC and BRACER showed the best overall performance across the tested datasets and conditions, whereas only BASIC demonstrated acceptable results on very short read libraries. Furthermore, the de novo assembly-based methods BRACER and BALDR were the most accurate in reconstructing BCRs harboring different degrees of SHMs in the variable domain, while TRUST4, MiXCR and BASIC were the fastest. Finally, we propose guidelines to select the best method based on the given data characteristics.

Comprehensive Characterizations of Immune Receptor Repertoire in Tumors and Cancer Immunotherapy Studies

We applied our computational algorithm TRUST4 to assemble immune receptor (T-cell receptor/B-cell receptor) repertoires from approximately 12,000 RNA sequencing samples from The Cancer Genome Atlas and seven immunotherapy studies. From over 35 million assembled complete complementary-determining region 3 sequences, we observed that the expression of CCL5 and MZB1 is the most positively correlated genes with T-cell clonal expansion and B-cell clonal expansion, respectively. We analyzed amino acid evolution during B-cell receptor somatic hypermutation and identified tyrosine as the preferred residue. We found that IgG1+IgG3 antibodies together with FcRn were associated with complement-dependent cytotoxicity and antibody-dependent cellular cytotoxicity or phagocytosis. In addition to B-cell infiltration, we discovered that B-cell clonal expansion and IgG1+IgG3 antibodies are also correlated with better patient outcomes. Finally, we created a website, VisualizIRR, for users to interactively explore and visualize the immune repertoires in this study. See related Spotlight by Liu and Han, p. 786.

Antibody repertoire sequencing analysis

High-throughput sequencing for B cell receptor (BCR) repertoire provides useful insights for the adaptive immune system. With the continuous development of the BCR-seq technology, many efforts have been made to develop methods for analyzing the ever-increasing BCR repertoire data. In this review, we comprehensively outline different BCR repertoire library preparation protocols and summarize three major steps of BCR-seq data analysis, i. e., V(D)J sequence annotation, clonal phylogenetic inference, and BCR repertoire profiling and mining. Different from other reviews in this field, we emphasize background intuition and the statistical principle of each method to help biologists better understand it. Finally, we discuss data mining problems for BCR-seq data and with a highlight on recently emerging multiple-sample analysis.

The contribution of uncharted RNA sequences to tumor identity in lung adenocarcinoma

The identity of cancer cells is defined by the interplay between genetic, epigenetic transcriptional and post-transcriptional variation. A lot of this variation is present in RNA-seq data and can be captured at once using reference-free, k-mer analysis. An important issue with k-mer analysis, however, is the difficulty of distinguishing signal from noise. Here, we use two independent lung adenocarcinoma datasets to identify all reproducible events at the k-mer level, in a tumor versus normal setting. We find reproducible events in many different locations (introns, intergenic, repeats) and forms (spliced, polyadenylated, chimeric etc.). We systematically analyze events that are ignored in conventional transcriptomics and assess their value as biomarkers and for tumor classification, survival prediction, neoantigen prediction and correlation with the immune microenvironment. We find that unannotated lincRNAs, novel splice variants, endogenous HERV, Line1 and Alu repeats and bacterial RNAs each contribute to different, important aspects of tumor identity. We argue that differential RNA-seq analysis of tumor/normal sample collections would benefit from this type k-mer analysis to cast a wider net on important cancer-related events. The code is available at https://github.com/Transipedia/dekupl-lung-cancer-inter-cohort.

A Pan-Cancer Analysis of Tumor-Infiltrating B Cell Repertoires

Tumor-infiltrating B cells can play an important role in anti-tumor responses but their presence is not well understood. In this study, we extracted the B cell receptor repertoires from 9522 tumor and adjacent non-tumor samples across 28 tumor types in the Cancer Genome Atlas project and performed diversity and network analysis. We identified differences in diversity and network statistics across tumor types and subtypes and observed a trend towards increased clonality in primary tumors compared to adjacent non-tumor tissues. We also found significant associations between the repertoire features and mutation load, tumor stage, and age. Our V-gene usage analysis identified similar V-gene usage patterns in colorectal and endometrial cancers. Lastly, we evaluated the prognostic value of the repertoire features and identified significant associations with survival in seven tumor types. This study warrants further research into better understanding the role of tumor-infiltrating B cells across a wide range of tumor types.

Macrophage-Derived Immunoglobulin M Inhibits Inflammatory Responses via Modulating Endoplasmic Reticulum Stress

Immunoglobulin (Ig), a characteristic marker of B cells, is a multifunctional evolutionary conserved antibody critical for maintaining tissue homeostasis and developing fully protective humoral responses to pathogens. Increasing evidence revealed that Ig is widely expressed in non-immune cells; moreover, Ig produced by different lineages cells plays different biological roles. Recently, it has been reported that monocytes or macrophages also express Ig. However, its function remains unclear. In this study, we further identified that Ig, especially Ig mu heavy chain (IgM), was mainly expressed in mice macrophages. We also analyzed the IgM repertoire characteristic in macrophages and found that the V_HDJ_H rearrangements of macrophage-derived IgM showed a restricted and conservative V_HDJ_H pattern, which differed from the diverse V_HDJ_H rearrangement pattern of the B cell-expressed IgM in an individual. Functional investigation showed that IgM knockdown significantly promoted macrophage migration and FAK/Src-Akt axis activation. Furthermore, some inflammatory cytokines such as MCP1 and IL-6 increased after IgM knockdown under LPS stimulation. A mechanism study revealed that the IgM interacted with binding immunoglobulin protein (Bip) and inhibited inflammatory response and unfolded protein response (UPR) activation in macrophages. Our data elucidate a previously unknown function of IgM in macrophages that explains its ability to act as a novel regulator of Bip to participate in endoplasmic reticulum stress and further regulate the inflammatory response.

Different B cell subpopulations show distinct patterns in their IgH repertoire metrics

Several human B cell subpopulations are recognised in the peripheral blood, which play distinct roles in the humoral immune response. These cells undergo developmental and maturational changes involving VDJ recombination, somatic hypermutation and class switch recombination, altogether shaping their immunoglobulin heavy chain (IgH) repertoire. Here, we sequenced the IgH repertoire of naïve, marginal zone, switched and plasma cells from 10 healthy adults along with matched unsorted and in silico separated CD19⁺ bulk B cells. Using advanced bioinformatic analysis and machine learning, we show that sorted B cell subpopulations are characterised by distinct repertoire characteristics on both the individual sequence and the repertoire level. Sorted subpopulations shared similar repertoire characteristics with their corresponding in silico separated subsets. Furthermore, certain IgH repertoire characteristics correlated with the position of the constant region on the IgH locus. Overall, this study provides unprecedented insight over mechanisms of B cell repertoire control in peripherally circulating B cell subpopulations.

Systematic memory B cell archiving and random display shape the human splenic marginal zone throughout life

Human memory B cells (MBCs) are generated and diversified in secondary lymphoid tissues throughout the organism. A paired immunoglobulin (Ig)-gene repertoire analysis of peripheral blood (PB) and splenic MBCs from infant, adult, and elderly humans revealed that throughout life, circulating MBCs are comprehensively archived in the spleen. Archive MBC clones are systematically preserved and uncoupled from class-switching. Clonality in the spleen increases steadily, but boosts at midlife, thereby outcompeting small clones. The splenic marginal zone (sMZ) represents a primed MBC compartment, generated from a stochastic exchange within the archive memory pool. This is supported by functional assays, showing that PB and splenic CD21⁺ MBCs acquire transient CD21^high expression upon NOTCH2-stimulation. Our study provides insight that the human MBC system in PB and spleen is composed of three interwoven compartments: the dynamic relationship of circulating, archive, and its subset of primed (sMZ) memory changes with age, thereby contributing to immune aging.

The T cell receptor repertoire of tumor infiltrating T cells is predictive and prognostic for cancer survival

Tumor infiltration by T cells is paramount for effective anti-cancer immune responses. We hypothesized that the T cell receptor (TCR) repertoire of tumor infiltrating T lymphocytes could therefore be indicative of the functional state of these cells and determine disease course at different stages in cancer progression. Here we show that the diversity of the TCR of tumor infiltrating T cell at baseline is prognostic in various cancers, whereas the TCR clonality of T cell infiltrating metastatic melanoma pre-treatment is predictive for activity and efficacy of PD1 blockade immunotherapy.

Kmerator Suite: design of specific k-mer signatures and automatic metadata discovery in large RNA-seq datasets

The huge body of publicly available RNA-sequencing (RNA-seq) libraries is a treasure of functional information allowing to quantify the expression of known or novel transcripts in tissues. However, transcript quantification commonly relies on alignment methods requiring a lot of computational resources and processing time, which does not scale easily to large datasets. K-mer decomposition constitutes a new way to process RNA-seq data for the identification of transcriptional signatures, as k-mers can be used to quantify accurately gene expression in a less resource-consuming way. We present the Kmerator Suite, a set of three tools designed to extract specific k-mer signatures, quantify these k-mers into RNA-seq datasets and quickly visualize large dataset characteristics. The core tool, Kmerator, produces specific k-mers for 97% of human genes, enabling the measure of gene expression with high accuracy in simulated datasets. KmerExploR, a direct application of Kmerator, uses a set of predictor gene-specific k-mers to infer metadata including library protocol, sample features or contaminations from RNA-seq datasets. KmerExploR results are visualized through a user-friendly interface. Moreover, we demonstrate that the Kmerator Suite can be used for advanced queries targeting known or new biomarkers such as mutations, gene fusions or long non-coding RNAs for human health applications.

TRUST4: immune repertoire reconstruction from bulk and single-cell RNA-seq data

We introduce the TRUST4 open-source algorithm for reconstruction of immune receptor repertoires in αβ/γδ T cells and B cells from RNA-sequencing (RNA-seq) data. Compared with competing methods, TRUST4 supports both FASTQ and BAM format and is faster and more sensitive in assembling longer-even full-length-receptor repertoires. TRUST4 can also call repertoire sequences from single-cell RNA-seq (scRNA-seq) data without V(D)J enrichment, and is compatible with both SMART-seq and 5' 10x Genomics platforms.