Methods for sequence and structural analysis of B and T cell receptor repertoires

Pan-cancer secreted proteome and skeletal muscle regulation: insight from a proteogenomic data-driven knowledge base

Large-scale, pan-cancer analysis is enabled by data driven knowledge bases that link tumor molecular profiles with phenotypes. A debilitating cancer-related phenotype is skeletal muscle loss, or cachexia, which occurs partly from tumor products secreted into circulation. Using the LinkedOmicsKB knowledge base assembled from the Clinical Proteomics Tumor Analysis Consortium proteogenomic analysis, along with catalogs of human secretome proteins, ligand-receptor pairs and molecular signatures, we sought to identify candidate pan-cancer proteins secreted to blood that could regulate skeletal muscle phenotypes in multiple solid cancers. Tumor proteins having significant pan-cancer associations with muscle were referenced against secretome proteins secreted to blood from the Human Protein Atlas, then verified as increased in paired tumor vs. normal tissues in pan-cancer manner. This workflow revealed seven secreted proteins from cancers afflicting kidneys, head and neck, lungs and pancreas that classified as protein-binding activity modulator, extracellular matrix protein or intercellular signaling molecule. Concordance of these biomarkers with validated molecular signatures of cachexia and senescence supported relevance to muscle and cachexia disease biology, and high tumor expression of the biomarker set associated with lower overall survival. In this article, we discuss avenues by which skeletal muscle and cachexia may be regulated by these candidate pan-cancer proteins secreted to blood, and conceptualize a strategy that considers them collectively as a biomarker signature with potential for refinement by data analytics and radiogenomics for predictive testing of future risk in a non-invasive, blood-based panel amenable to broad uptake and early management.

Rational computational design and development of an immunogenic multiepitope vaccine incorporating transmembrane proteins of Staphylococcus lugdunensis

Staphylococcus lugdunensis has emerged as a significant human pathogen, responsible for a range of infections from skin and soft tissue infections to endocarditis and bacteremia. Notably, abscess formation is a common manifestation, reflecting its potential shift from a benign skin commensal to a serious pathogen, akin to infective endocarditis. With the rising prevalence of antibiotic resistance, there is a pressing need for novel therapeutic strategies. This study addresses this need by exploring the development of an effective S. lugdunensis vaccine. Multiepitope vaccines, which incorporate various antigenic fragments from S. lugdunensis proteins, offer a promising approach to elicit a robust immune response. Computational tools are instrumental in selecting epitopes based on their predicted immunogenicity and non-toxicity. Molecular docking and molecular dynamics (MD) simulations further elucidate the interactions between vaccine constructs and immune system molecules, such as B-cell and T-cell receptors, providing detailed insights into binding affinity, specificity, and stability. This study highlights the potential of integrating multiepitope vaccine design with advanced computational methods to expedite and enhance vaccine development, addressing a critical gap amid escalating antibiotic resistance.

Similarity Network Analysis of the Adaptive Immune Response in the Proximal Airway

OBJECTIVES

Recent immunologic study of the adaptive immune repertoire in the subglottic airway demonstrated high-frequency T cell clones that do not overlap between individuals. However, the anatomic distribution and antigenic target of the T cell repertoire in the proximal airway mucosa remain unresolved.

METHODS

Single-cell RNA sequencing of matched scar and unaffected mucosa from idiopathic subglottic stenosis patients (iSGS, n = 32) was performed and compared with airway mucosa from healthy controls (n = 10). T cell receptor (TCR) sequences were interrogated via similarity network analysis to explore antigenic targets using the published algorithm: Grouping of Lymphocyte Interactions by Paratope Hotspots (GLIPH2).

RESULTS

The mucosal T cell repertoire in healthy control airways consisted of highly expressed T cell clones conserved across anatomic subsites (trachea, bronchi, bronchioles, and lung). In iSGS, high-frequency clones were equally represented in both scar and adjacent non-scar tissue. Significant differences in repertoire structure between iSGS scar and unaffected mucosa was observed, driven by unique low-frequency clones. GLIPH2 results suggest low-frequency clones share targets between multiple iSGS patients.

CONCLUSION

Healthy airway mucosa has a highly conserved T cell repertoire across multiple anatomic subsites. Similarly, iSGS patients have highly expressed T cell clones present in both scar and unaffected mucosa. iSGS airway scar possesses an abundance of less highly expanded clones with predicted antigen targets shared between patients. Interrogation of these shared motifs suggests abundant adaptive immunity to viral targets in iSGS airway scar. These results provide insight into disease pathogenesis and illuminate new treatment strategies in iSGS.

LEVEL OF EVIDENCE

NA Laryngoscope, 134:3245-3252, 2024.

Identification of PLAC8 as a Potential Biomarker for the Diagnosis of Interstitial Cystitis

BACKGROUND

Interstitial cystitis is a diagnosis of exclusion due to the complexity of its etiology and pathology, which is a chronic disease with an unknown etiology. To our knowledge, few studies were performed to identify predictive biomarkers for interstitial cystitis.

OBJECTIVE

This study aimed to identify and validate potential biomarkers for Interstitial Cystitis (IC).

METHODS

The interstitial cystitis datasets were retrieved from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were identified by using the R package and were subjected to functional and pathway enrichment analysis. Key biomarkers of interstitial cystitis were identified by using Lasso regression analysis and the SVM-RFE algorithm. The diagnostic value of key biomarkers was validated in internal and external datasets, and pathways that relate to biomarkers of interstitial cystitis were screened. The ssGSEA was employed to identify the immune cells closely related to biomarkers. The expression of PLAC8 in patients with interstitial cystitis was detected by Immune-Histochemistry (IHC).

RESULTS

Sixteen differentially expressed genes associated with interstitial cystitis were identified, which were primarily linked to the biological process of the chemokine signaling pathway. PLAC8, identified as a biomarker for interstitial cystitis, was validated to express a significantly different between IC and normal bladder tissues. PLAC8-related pathways were analyzed, with a focus on NF-κB, TNF, Toll-like receptor, chemokine, IL-17, and JAK-STAT signaling pathways. PLAC8 was proved to be closely related to immune activations, which is similar to the pathogenesis of IC, which is a chronic dysregulated immune disease. Meanwhile, we also observed a higher level of PLAC8 in IC tissues.

CONCLUSION

PLAC8 has promising application prospects as a biomarker for interstitial cystitis diagnosis. These findings could aid in the diagnosis and treatment of interstitial cystitis.

SPLASH: A statistical, reference-free genomic algorithm unifies biological discovery

Today's genomics workflows typically require alignment to a reference sequence, which limits discovery. We introduce a unifying paradigm, SPLASH (Statistically Primary aLignment Agnostic Sequence Homing), which directly analyzes raw sequencing data, using a statistical test to detect a signature of regulation: sample-specific sequence variation. SPLASH detects many types of variation and can be efficiently run at scale. We show that SPLASH identifies complex mutation patterns in SARS-CoV-2, discovers regulated RNA isoforms at the single-cell level, detects the vast sequence diversity of adaptive immune receptors, and uncovers biology in non-model organisms undocumented in their reference genomes: geographic and seasonal variation and diatom association in eelgrass, an oceanic plant impacted by climate change, and tissue-specific transcripts in octopus. SPLASH is a unifying approach to genomic analysis that enables expansive discovery without metadata or references.

Improved prediction of MHC-peptide binding using protein language models

Major histocompatibility complex Class I (MHC-I) molecules bind to peptides derived from intracellular antigens and present them on the surface of cells, allowing the immune system (T cells) to detect them. Elucidating the process of this presentation is essential for regulation and potential manipulation of the cellular immune system. Predicting whether a given peptide binds to an MHC molecule is an important step in the above process and has motivated the introduction of many computational approaches to address this problem. NetMHCPan, a pan-specific model for predicting binding of peptides to any MHC molecule, is one of the most widely used methods which focuses on solving this binary classification problem using shallow neural networks. The recent successful results of Deep Learning (DL) methods, especially Natural Language Processing (NLP-based) pretrained models in various applications, including protein structure determination, motivated us to explore their use in this problem. Specifically, we consider the application of deep learning models pretrained on large datasets of protein sequences to predict MHC Class I-peptide binding. Using the standard performance metrics in this area, and the same training and test sets, we show that our models outperform NetMHCpan4.1, currently considered as the-state-of-the-art.

DeepAIR: A deep learning framework for effective integration of sequence and 3D structure to enable adaptive immune receptor analysis

Structural docking between the adaptive immune receptors (AIRs), including T cell receptors (TCRs) and B cell receptors (BCRs), and their cognate antigens are one of the most fundamental processes in adaptive immunity. However, current methods for predicting AIR-antigen binding largely rely on sequence-derived features of AIRs, omitting the structure features that are essential for binding affinity. In this study, we present a deep learning framework, termed DeepAIR, for the accurate prediction of AIR-antigen binding by integrating both sequence and structure features of AIRs. DeepAIR achieves a Pearson's correlation of 0.813 in predicting the binding affinity of TCR, and a median area under the receiver-operating characteristic curve (AUC) of 0.904 and 0.942 in predicting the binding reactivity of TCR and BCR, respectively. Meanwhile, using TCR and BCR repertoire, DeepAIR correctly identifies every patient with nasopharyngeal carcinoma and inflammatory bowel disease in test data. Thus, DeepAIR improves the AIR-antigen binding prediction that facilitates the study of adaptive immunity.

Tutorial: integrative computational analysis of bulk RNA-sequencing data to characterize tumor immunity using RIMA

RNA-sequencing (RNA-seq) has become an increasingly cost-effective technique for molecular profiling and immune characterization of tumors. In the past decade, many computational tools have been developed to characterize tumor immunity from gene expression data. However, the analysis of large-scale RNA-seq data requires bioinformatics proficiency, large computational resources and cancer genomics and immunology knowledge. In this tutorial, we provide an overview of computational analysis of bulk RNA-seq data for immune characterization of tumors and introduce commonly used computational tools with relevance to cancer immunology and immunotherapy. These tools have diverse functions such as evaluation of expression signatures, estimation of immune infiltration, inference of the immune repertoire, prediction of immunotherapy response, neoantigen detection and microbiome quantification. We describe the RNA-seq IMmune Analysis (RIMA) pipeline integrating many of these tools to streamline RNA-seq analysis. We also developed a comprehensive and user-friendly guide in the form of a GitBook with text and video demos to assist users in analyzing bulk RNA-seq data for immune characterization at both individual sample and cohort levels by using RIMA.

[WITHDRAWN] SPLASH: a statistical, reference-free genomic algorithm unifies biological discovery

The authors have withdrawn this manuscript due to a duplicate posting of manuscript number BIORXIV/2022/497555. Therefore, the authors do not wish this work to be cited as reference for the project. If you have any questions, please contact the corresponding author. The correct preprint can be found at doi: https://doi.org/10.1101/2022.06.24.497555.

TPBTE: A model based on convolutional Transformer for predicting the binding of TCR to epitope

T cell receptors (TCRs) selectively bind to antigens to fight pathogens with specific immunity. Current tools focus on the nature of amino acids within sequences and take less into account the nature of amino acids far apart and the relationship between sequences, leading to significant differences in the results from different datasets. We propose TPBTE, a model based on convolutional Transformer for Predicting the Binding of TCR to Epitope. It takes epitope sequences and the complementary decision region 3 (CDR3) sequences of TCRβ chain as inputs. And it uses a convolutional attention mechanism to learn amino acid representations between different positions of the sequences based on learning local features of the sequences. At the same time, it uses cross attention to learn the interaction information between TCR sequences and epitope sequences. A comprehensive evaluation of the TCR-epitope data shows that the average area under the curve of TPBTE outperforms the baseline model, and demonstrate an intentional performance. In addition, TPBTE can give the probability of binding TCR to epitopes, which can be used as the first step of epitope screening, narrowing the scope of epitope search and reducing the time of epitope search.

Progress in kidney transplantation: The role for systems immunology

The development of systems biology represents an immense breakthrough in our ability to perform translational research and deliver personalized and precision medicine. A multidisciplinary approach in combination with use of novel techniques allows for the extraction and analysis of vast quantities of data even from the volume and source limited samples that can be obtained from human subjects. Continued advances in microfluidics, scalability and affordability of sequencing technologies, and development of data analysis tools have made the application of a multi-omics, or systems, approach more accessible for use outside of specialized centers. The study of alloimmune and protective immune responses after solid organ transplant offers innumerable opportunities for a multi-omics approach, however, transplant immunology labs are only just beginning to adopt the systems methodology. In this review, we focus on advances in biological techniques and how they are improving our understanding of the immune system and its interactions, highlighting potential applications in transplant immunology. First, we describe the techniques that are available, with emphasis on major advances that allow for increased scalability. Then, we review initial applications in the field of transplantation with a focus on topics that are nearing clinical integration. Finally, we examine major barriers to adapting these methods and discuss potential future developments.

Peripheral B cells from patients with hepatitis C virus-associated lymphoma exhibit clonal expansion and an anergic-like transcriptional profile

Chronic HCV infection remains a global health concern due to its involvement in hepatic and extrahepatic diseases, including B cell non-Hodgkin lymphoma (BNHL). Clinical and epidemiological evidence support a causal role for HCV in BNHL development, although mechanistic insight is lacking. We performed RNA-sequencing on peripheral B cells from patients with HCV alone, BNHL alone, and HCV-associated BNHL to identify unique and shared transcriptional profiles associated with transformation. In patients with HCV-associated BNHL, we observed the enrichment of an anergic-like gene signature and evidence of clonal expansion that was correlated with the expression of epigenetic regulatory genes. Our data support a role for viral-mediated clonal expansion of anergic-like B cells in HCV-associated BNHL development and suggest epigenetic dysregulation as a potential mechanism driving expansion. We propose epigenetic mechanisms may be involved in both HCV-associated lymphoma and regulation of B cell anergy, representing an attractive target for clinical interventions.

Machine Learning Approaches to TCR Repertoire Analysis

Sparked by the development of genome sequencing technology, the quantity and quality of data handled in immunological research have been changing dramatically. Various data and database platforms are now driving the rapid progress of machine learning for immunological data analysis. Of various topics in immunology, T cell receptor repertoire analysis is one of the most important targets of machine learning for assessing the state and abnormalities of immune systems. In this paper, we review recent repertoire analysis methods based on machine learning and deep learning and discuss their prospects.

A SARS-CoV-2 neutralizing antibody discovery by single cell sequencing and molecular modeling

Although vaccines have been developed, mutations of SARS-CoV-2, especially the dominant B.1.617.2 (delta) and B.1.529 (omicron) strains with more than 30 mutations on their spike protein, have caused a significant decline in prophylaxis, calling for the need for drug improvement. Antibodies are drugs preferentially used in infectious diseases and are easy to get from immunized organisms. The current study combined molecular modeling and single memory B cell sequencing to assess candidate sequences before experiments, providing a strategy for the fabrication of SARS-CoV-2 neutralizing antibodies. A total of 128 sequences were obtained after sequencing 196 memory B cells, and 42 sequences were left after merging extremely similar ones and discarding incomplete ones, followed by homology modeling of the antibody variable region. Thirteen candidate sequences were expressed, of which three were tested positive for receptor binding domain recognition but only one was confirmed as having broad neutralization against several SARS-CoV-2 variants. The current study successfully obtained a SARS-CoV-2 antibody with broad neutralizing abilities and provided a strategy for antibody development in emerging infectious diseases using single memory B cell BCR sequencing and computer assistance in antibody fabrication.

Advances in Lymphoma Molecular Diagnostics

Lymphomas encompass a diverse group of malignant lymphoid neoplasms. Over recent years much scientific effort has been undertaken to identify and understand molecular changes in lymphomas, resulting in a wide range of genetic alterations that have been reported across all types of lymphomas. As many of these changes are now incorporated into the World Health Organization’s defined criteria for the diagnostic evaluation of patients with lymphoid neoplasms, their accurate identification is crucial. Even if many alterations are not routinely evaluated in daily clinical practice, they may still have implications in risk stratification, treatment, prognosis or disease monitoring. Moreover, some alterations can be used for targeted treatment. Therefore, these advances in lymphoma molecular diagnostics in some cases have led to changes in treatment algorithms. Here, we give an overview of and discuss advances in molecular techniques in current clinical practice, as well as highlight some of them in a clinical context.

Brain Immunoinformatics: A Symmetrical Link between Informatics, Wet Lab and the Clinic

Breakthrough advances in informatics over the last decade have thoroughly influenced the field of immunology. The intermingling of machine learning with wet lab applications and clinical results has hatched the newly defined immunoinformatics society. Immunoinformatics of the central neural system, referred to as neuroimmunoinformatics (NII), investigates symmetrical and asymmetrical interactions of the brain-immune interface. This interdisciplinary overview on NII is addressed to bioscientists and computer scientists. We delineate the dominating trajectories and field-shaping achievements and elaborate on future directions using bridging language and terminology. Computation, varying from linear modeling to complex deep learning approaches, fuels neuroimmunology through three core directions. Firstly, by providing big-data analysis software for high-throughput methods such as next-generation sequencing and genome-wide association studies. Secondly, by designing models for the prediction of protein morphology, functions, and symmetrical and asymmetrical protein–protein interactions. Finally, NII boosts the output of quantitative pathology by enabling the automatization of tedious processes such as cell counting, tracing, and arbor analysis. The new classification of microglia, the brain’s innate immune cells, was an NII achievement. Deep sequencing classifies microglia in “sensotypes” to accurately describe the versatility of immune responses to physiological and pathological challenges, as well as to experimental conditions such as xenografting and organoids. NII approaches complex tasks in the brain-immune interface, recognizes patterns and allows for hypothesis-free predictions with ultimate targeted individualized treatment strategies, and personalizes disease prognosis and treatment response.

Progress and challenges in mass spectrometry-based analysis of antibody repertoires

Humoral immunity is divided into the cellular B cell and protein-level antibody responses. High-throughput sequencing has advanced our understanding of both these fundamental aspects of B cell immunology as well as aspects pertaining to vaccine and therapeutics biotechnology. Although the protein-level serum and mucosal antibody repertoire make major contributions to humoral protection, the sequence composition and dynamics of antibody repertoires remain underexplored. This limits insight into important immunological and biotechnological parameters such as the number of antigen-specific antibodies, which are for example, relevant for pathogen neutralization, microbiota regulation, severity of autoimmunity, and therapeutic efficacy. High-resolution mass spectrometry (MS) has allowed initial insights into the antibody repertoire. We outline current challenges in MS-based sequence analysis of antibody repertoires and propose strategies for their resolution.

Predicting T Cell Receptor Antigen Specificity From Structural Features Derived From Homology Models of Receptor-Peptide-Major Histocompatibility Complexes

The physical interaction between the T cell receptor (TCR) and its cognate antigen causes T cells to activate and participate in the immune response. Understanding this physical interaction is important in predicting TCR binding to a target epitope, as well as potential cross-reactivity. Here, we propose a way of collecting informative features of the binding interface from homology models of T cell receptor-peptide-major histocompatibility complex (TCR-pMHC) complexes. The information collected from these structures is sufficient to discriminate binding from non-binding TCR-pMHC pairs in multiple independent datasets. The classifier is limited by the number of crystal structures available for the homology modelling and by the size of the training set. However, the classifier shows comparable performance to sequence-based classifiers requiring much larger training sets.

Receptor repertoires of murine follicular T helper cells reveal a high clonal overlap in separate lymph nodes in autoimmunity

Follicular T helper cells (Tfh) are a specialized subset of CD4 effector T cells that are crucial for germinal center (GC) reactions and for selecting B cells to undergo affinity maturation. Despite this central role for humoral immunity, only few data exist about their clonal distribution when multiple lymphoid organs are exposed to the same antigen (Ag) as it is the case in autoimmunity. Here, we used an autoantibody-mediated disease model of the skin and injected one auto-Ag into the two footpads of the same mouse and analyzed the T cell receptor (TCR)β sequences of Tfh located in GCs of both contralateral draining lymph nodes. We found that over 90% of the dominant GC-Tfh clonotypes were shared in both lymph nodes but only transiently. The initially dominant Tfh clonotypes especially declined after establishment of chronic disease while GC reaction and autoimmune disease continued. Our data demonstrates a dynamic behavior of Tfh clonotypes under autoimmune conditions and emphasizes the importance of the time point for distinguishing auto-Ag-specific Tfh clonotypes from potential bystander activated ones.

Predictive Markers of Immunogenicity and Efficacy for Human Vaccines

Vaccines represent one of the major advances of modern medicine. Despite the many successes of vaccination, continuous efforts to design new vaccines are needed to fight "old" pandemics, such as tuberculosis and malaria, as well as emerging pathogens, such as Zika virus and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Vaccination aims at reaching sterilizing immunity, however assessing vaccine efficacy is still challenging and underscores the need for a better understanding of immune protective responses. Identifying reliable predictive markers of immunogenicity can help to select and develop promising vaccine candidates during early preclinical studies and can lead to improved, personalized, vaccination strategies. A systems biology approach is increasingly being adopted to address these major challenges using multiple high-dimensional technologies combined with in silico models. Although the goal is to develop predictive models of vaccine efficacy in humans, applying this approach to animal models empowers basic and translational vaccine research. In this review, we provide an overview of vaccine immune signatures in preclinical models, as well as in target human populations. We also discuss high-throughput technologies used to probe vaccine-induced responses, along with data analysis and computational methodologies applied to the predictive modeling of vaccine efficacy.

Characterization of amino acid residues of T-cell receptors interacting with HLA-A*02-restricted antigen peptides

Background

The present study aimed to explore residues’ properties interacting with HLA-A*02-restricted peptides on T-cell receptors (TCRs) and their effects on bond types of interaction and binding free energy.

Methods

We searched the crystal structures of HLA-A*02-restricted peptide-TCR complexes from the Protein Data Bank (PDB) database and subsequently collected relevant parameters. We then employed Schrodinger to analyze the bond types of interaction and Gromacs 2019 to evaluate the TCR-antigen peptide complex’s molecular dynamics simulation. Finally, we compared the changes of bond types of interaction and binding free energy before and after residue substitution to ensure consistency of the conditions before and after residue substitution.

Results

The main sites on the antigen peptides that formed the intermolecular interaction [hydrogen bond (HB) and pi stack] with TCRs were P4, P8, P2, and P6. The hydrophobicity of the amino acids inside or outside the disulfide bond of TCRs may be related to the intermolecular interaction and binding free energy between TCRs and peptides. Residues located outside the disulfide bond of TCR α or β chains and forming pi stack force played favorable roles in the complex intermolecular interaction and binding free energy. The residues of the TCR α or β chains that interacted with peptides were replaced by alanine (Ala) or glycine (Gly), and their intermolecular binding free energy of the complex had been improved. However, it had nothing to do with the formation of HB.

Conclusions

The findings of this study suggest that the hydrophobic nature of the amino acids inside or outside the disulfide bonds on the TCR may be associated with the intermolecular interaction and binding between the TCR and polypeptide. The residues located outside the TCR α or β single-chain disulfide bond and forming the pi-stack force showed a beneficial effect on the intermolecular interaction and binding of the complex. In addition, the part of the residues on the TCR α or β single chain that produced bond types of interaction with the polypeptide after being replaced by Ala or Gly, the intermolecular binding free energy of the complex was increased, regardless of whether HB was formed.

Progress toward improved understanding of antibody maturation

Upon encountering an antigen, antibodies mature through various rounds of somatic mutations, resulting in higher affinities and specificities to the particular antigen. We review recent progress in four areas of antibody maturation studies. (1) Next-generation and single-cell sequencing have revolutionized the analysis of antibody repertoires by dramatically increasing the sequences available to study the state and evolution of the immune system. Computational methods, including machine learning tools, have been developed for reconstituting antibody clonal lineages and for general repertoire analysis. (2) The availability of X-ray structures, thermodynamic and kinetic data, and molecular dynamics simulations provide information on the biophysical mechanisms responsible for improved affinity. (3) In addition to improved binding to a specific antigen, providing affinity-independent diversity and self/nonself discrimination are fundamental functions of the immune system. Recent studies, including X-ray structures, yield improved understanding of both mechanisms. (4) Results from in vivo maturation help to develop methods of in vitro maturation to improve antibody properties for therapeutic applications, frequently combining computational and experimental approaches.

Approaching Shared Pathophysiology in Immune-Mediated Diseases through Functional Genomics

Immune-mediated diseases (IMDs) are complex pathologies that are strongly influenced by environmental and genetic factors. Associations between genetic loci and susceptibility to these diseases have been widely studied, and hundreds of risk variants have emerged during the last two decades, with researchers observing a shared genetic pattern among them. Nevertheless, the pathological mechanism behind these associations remains a challenge that has just started to be understood thanks to functional genomic approaches. Transcriptomics, regulatory elements, chromatin interactome, as well as the experimental characterization of genomic findings, constitute key elements in the emerging understandings of how genetics affects the etiopathogenesis of IMDs. In this review, we will focus on the latest advances in the field of functional genomics, centering our attention on systemic rheumatic IMDs.