HLA-C: An Accomplice in Rheumatic Diseases

HLA variations in patients with diffuse large B-cell lymphoma and association with disease risk and prognosis: a case-control study

Introduction

Human leukocyte antigen (HLA) polymorphisms have been associated with the development of various autoimmune diseases, as well as malignant neoplasms. Non-Hodgkin lymphomas (NHLs) are a heterogenous group of lymphoid malignancies in which a genetic substrate has been established and is deemed to play a crucial role in disease pathogenesis. This study aimed to identify whether variations in the HLA gene region were associated with diffuse large B-cell lymphoma (DLBCL) risk and prognosis.

Methods

We defined HLA class I (HLA-A, HLA-B, HLA-C) and class II (HLA-DRB1, HLA-DQB1) alleles in 60 patients with DLBCL and compared the results to those found by 236 healthy adult donors from the bone marrow bank of Northern Greece. HLA typing was performed by two molecular methods, Sequence - Specific Oligonucleotide HLA typing (SSO) and Sequence - Specific Primer HLA typing (SSP), from white blood cells recovered from peripheral blood. The phenotypic frequencies of HLA-A, HLA-B, HLA-C, HLA-DRB1 and HLA-DQB1 between patients and controls were compared with the 2-sided Fisher's exact test. Results with p-value <0.05 were considered statistically significant. Odds Ratios with 95% Confidence Intervals were calculated to further strengthen the results. The 2-sided Fisher's exact test was also applied to alleles found only in one of the two groups, while the odds ratios together with the confidence intervals were corrected with Haldane-Anscombe method.

Results

Among the studied HLA polymorphisms, the frequency HLA-C*12 allele was significantly lower in patients with DLBCL compared with control subjects (6.7% vs. 34.7%, OR = 0.16, 95% CI: 0.04-0.44). Frequency of HLA-B*39 was significantly lower in patients with DLBCL compared with controls, but due to the low frequency of this polymorphism in the studied population and small sample size, determinations regarding the significance of this findings were limited. Survival analysis revealed that the presence of HLA-C*12 was not associated with improved or worsened overall and progression-free survival. No statistically significant associations were observed in the phenotypic frequencies of HLA-A, HLA-DQB1, HLA-DRB1 and the rest of HLA-B alleles between the control and DLBCL groups.

Discussion

Collectively, our results provide valuable insight regarding the role of HLA variations on DLBCL risk. Further studies are required to consolidate our findings and ascertain the clinical implications of these genetic variations on DLBCL management and prognosis.

High-resolution HLA genotyping identifies risk alleles in both class I and II for primary autoimmune neutropenia in early childhood in a Danish cohort

HLA studies in patients with autoimmune neutropenia (AIN) have shown very consistent results for the association with HLA class II alleles at low resolution. This study aimed to examine the association of both HLA class I and class II at high resolution to clarify the contribution of risk alleles to the disease. A total of 107 AIN patients were genotyped for six loci of HLA class I (HLA-A, -B and -C) and class II (HLA-DRB1, -DQB1, and -DPB1) genes by a high-resolution (3-field, 6-digit) analysis and compared with HLA typing of 1000 healthy controls. Compared with the controls, the allele frequencies were significantly higher in AIN patients for A*02:17:01G, C*01:02:01G, DRB1*10:01:01G, DRB1*14:01:01G, DRB1*16:01:01G, DQB1*05:02:01G, and DQB1*05:03:01G but lower significant for C*03:04:01G, DRB1*04:01:01G, DRB1*13:02:01G, DQB1*03:02:01G, and DQB1*06:04:01G. Frequently associated two-locus haplotypes were found to be DRB1*10:01:01G-DQB1*05:01:01G and DRB1*16:01:01G-DQB1*05:02:01G, while the S2 (Q- or D-KRAA) shared epitope (SE) was associated with lower risk. A unique association with HLA alleles was observed between patients with specific anti-HNA-1a antibodies and broad-reacting anti-FcγRIIIb. Anti-HNA-1a antibody-positive patients were associated with C*01:02:01G, DRB1*01:01:01G, DRB1*16:01:01G, DQB1*05:01:01G, DQB1*05:02:01G, DQB1*06:04:01G, and DPB1*10:01:01G; the two-locus haplotypes DRB1*01:01:01G-DQB1*05:01:01G and DRB1*16:01:01G-DQB1*05:02:01G; and the S3P (Q- or R-RRAA) SE. Anti-FcγRIIIb antibody-positive patients were associated with the alleles A*02:17:01G, DRB1*10:01:01G, and DQB1*05:02:01G; the haplotypes DRB1*10:01:01G-DQB1*05:01:01G and DRB1*11:01:02G-DQB1*05:02:01G; and the S3D (DRRAA) SE. The different associations regarding FcγRIIIb antibody specificities could indicate disease heterogeneity.

Neuroimmune Mechanisms Underlying Post-acute Sequelae of SARS-CoV-2 (PASC) Pain, Predictions from a Ligand-Receptor Interactome

PURPOSE OF REVIEW

Individuals with post-acute sequelae of SARS-CoV-2 (PASC) complain of persistent musculoskeletal pain. Determining how COVID-19 infection produces persistent pain would be valuable for the development of therapeutics aimed at alleviating these symptoms.

RECENT FINDINGS

To generate hypotheses regarding neuroimmune interactions in PASC, we used a ligand-receptor interactome to make predictions about how ligands from PBMCs in individuals with COVID-19 communicate with dorsal root ganglia (DRG) neurons to induce persistent pain. In a structured literature review of -omics COVID-19 studies, we identified ligands capable of binding to receptors on DRG neurons, which stimulate signaling pathways including immune cell activation and chemotaxis, the complement system, and type I interferon signaling. The most consistent finding across immune cell types was an upregulation of genes encoding the alarmins S100A8/9 and MHC-I. This ligand-receptor interactome, from our hypothesis-generating literature review, can be used to guide future research surrounding mechanisms of PASC-induced pain.

Clinical trait-specific genetic analysis in Behçet's disease identifies novel loci associated with ocular and neurological involvement

Behçet's disease is a complex inflammatory vasculitis with a broad spectrum of clinical manifestations. The purpose of this study was to investigate the genetics underlying specific clinical features of Behçet's disease. A total of 436 patients with Behçet's disease from Turkey were studied. Genotyping was performed using the Infinium ImmunoArray-24 BeadChip. After imputation and quality control measures, logistic regressions adjusting for sex and the first five principal components were performed for each clinical trait using a case-case genetic analysis approach. A weighted genetic risk score was calculated for each clinical feature. Genetic association analyses of previously identified susceptibility loci in Behçet's disease revealed a genetic association between ocular lesions and HLA-B/MICA (rs116799036: OR = 1.85 [95% CI = 1.35-2.52], p-value = 1.1 × 10-4). The genetic risk score was significantly higher in Behçet's disease patients with ocular lesions compared to those without ocular involvement, which is explained by the genetic variation in the HLA region. New genetic loci predisposing to specific clinical features in Behçet's disease were suggested when genome-wide variants were evaluated. The most significant associations were observed in ocular involvement with SLCO4A1 (rs6062789: OR = 0.41 [95% CI = 0.30-0.58], p-value = 1.92 × 10-7), and neurological involvement with DDX60L (rs62334264: OR = 4.12 [95% CI 2.34 to 7.24], p-value = 8.85 × 10-7). Our results emphasize the role of genetic factors in predisposing to specific clinical manifestations in Behçet's disease, and might shed additional light into disease heterogeneity, pathogenesis, and variability of Behçet's disease presentation across populations.

The novel HLA-C allele, HLA-C*07:02:141 was identified in a Chinese individual

HLA-C*07:02:141 shows one nucleotide substitution at position 4 when compared to HLA-C*07:02:01:01.

"HLA-C: evolution, epigenetics, and pathological implications in the major histocompatibility complex"

HLA-C, a gene located within the major histocompatibility complex, has emerged as a prominent target in biomedical research due to its involvement in various diseases, including cancer and autoimmune disorders; even though its recent addition to the MHC, the interaction between HLA-C and KIR is crucial for immune responses, particularly in viral infections. This review provides an overview of the structure, origin, function, and pathological implications of HLA-C in the major histocompatibility complex. In the last decade, we systematically reviewed original publications from Pubmed, ScienceDirect, Scopus, and Google Scholar. Our findings reveal that genetic variations in HLA-C can determine susceptibility or resistance to certain diseases. However, the first four exons of HLA-C are particularly susceptible to epigenetic modifications, which can lead to gene silencing and alterations in immune function. These alterations can manifest in diseases such as alopecia areata and psoriasis and can also impact susceptibility to cancer and the effectiveness of cancer treatments. By comprehending the intricate interplay between genetic and epigenetic factors that regulate HLA-C expression, researchers may develop novel strategies for preventing and treating diseases associated with HLA-C dysregulation.

Role of HLA-I Structural Variants and the Polyreactive Antibodies They Generate in Immune Homeostasis

Cell-surface HLA-I molecules consisting of β2-microglobulin (β2m) associated heavy chains (HCs), referred to as Face-1, primarily present peptides to CD8+ T-cells. HCs consist of three α-domains, with selected amino acid sequences shared by all alleles of all six isoforms. The cell-surface HLA undergoes changes upon activation by pathological conditions with the expression of β2m-free HCs (Face-2) resulting in exposure of β2m-masked sequences shared by almost all alleles and the generation of HLA-polyreactive antibodies (Abs) against them. Face-2 may homodimerize or heterodimerize with the same (Face-3) or different alleles (Face-4) preventing exposure of shared epitopes. Non-allo immunized males naturally carry HLA-polyreactive Abs. The therapeutic intravenous immunoglobulin (IVIg) purified from plasma of thousands of donors contains HLA-polyreactive Abs, admixed with non-HLA Abs. Purified HLA-polyreactive monoclonal Abs (TFL-006/007) generated in mice after immunizing with Face-2 are documented to be immunoregulatory by suppressing or activating different human lymphocytes, much better than IVIg. Our objectives are (a) to elucidate the complexity of the HLA-I structural variants, and their Abs that bind to both shared and uncommon epitopes on different variants, and (b) to examine the roles of those Abs against HLA-variants in maintaining immune homeostasis. These may enable the development of personalized therapeutic strategies for various pathological conditions.

Association of anti-HSC70 autoantibodies with cutaneous ulceration and severe disease in juvenile dermatomyositis

OBJECTIVES

JDM is an inflammatory myopathy characterized by prominent vasculopathy. AECAs are frequently detected in inflammatory and autoimmune diseases. We sought to determine whether AECAs correlate with clinical features of JDM, and thus serve as biomarkers to guide therapy or predict outcome.

METHODS

Plasma samples from 63 patients with JDM, 49 patients with polyarticular JIA and 40 juvenile healthy controls were used to detect anti-heat shock cognate 71 kDa protein (HSC70) autoantibodies, a newly identified AECA, in ELISA assays. Clinical features were compared between JDM patients with and without anti-HSC70 autoantibodies.

RESULTS

Anti-HSC70 autoantibodies were detected in 35% of patients with JDM, in 0% of patients with JIA (P < 0.0001) and in 0% of healthy donors (P < 0.0001). Both the presence of cutaneous ulcers (59% vs 17%, P < 0.002) and the use of wheelchairs and/or assistive devices (64% vs 27%, P < 0.007) were strongly associated with anti-HSC70 autoantibodies in JDM. High scores on the severity of myositis damage measures at the time of measurement of anti-HSC70 autoantibodies and an increased number of hospitalizations were also associated with anti-HSC70 autoantibodies. Intravenous immunoglobulin therapy was used more often in anti-HSC70 autoantibody-positive patients.

CONCLUSION

Anti-HCS70 autoantibodies are detected frequently in children with JDM and are novel myositis-associated autoantibodies correlating with disease severity.

A Summary on the Genetics of Systemic Lupus Erythematosus, Rheumatoid Arthritis, Systemic Sclerosis, and Sjögren's Syndrome

Systemic lupus erythematosus, systemic sclerosis, rheumatoid arthritis, and Sjögren's syndrome are four major autoimmune rheumatic diseases characterized by the presence of autoantibodies, caused by a dysregulation of the immune system that leads to a wide variety of clinical manifestations. These conditions present complex etiologies strongly influenced by multiple environmental and genetic factors. The human leukocyte antigen (HLA) region was the first locus identified to be associated and still represents the strongest susceptibility factor for each of these conditions, particularly the HLA class II genes, including DQA1, DQB1, and DRB1, but class I genes have also been associated. Over the last two decades, the genetic component of these disorders has been extensively investigated and hundreds of non-HLA risk genetic variants have been uncovered. Furthermore, it is widely accepted that autoimmune rheumatic diseases share molecular disease pathways, such as the interferon (IFN) type I pathways, which are reflected in a common genetic background. Some examples of well-known pleiotropic loci for autoimmune rheumatic diseases are the HLA region, DNASEL13, TNIP1, and IRF5, among others. The identification of the causal molecular mechanisms behind the genetic associations is still a challenge. However, recent advances have been achieved through mouse models and functional studies of the loci. Here, we provide an updated overview of the genetic architecture underlying these four autoimmune rheumatic diseases, with a special focus on the HLA region.

Comprehensive Analysis of Colorectal Cancer Immunity and Identification of Immune-Related Prognostic Targets

Colorectal cancer (COAD) is ranked as the third most common cancer and second in terms of cancer-related deaths worldwide. Due to its poor overall survival and prognosis, the incidents of COAD are significantly increasing. Although treatment methods have greatly been improved in the last decade, it is still not good enough to have satisfactory treatment outcomes. In recent years, immunotherapy has been successful to some extent in the treatment of many cancers but still, many patients do not respond to immunotherapy. Therefore, it is essential to have a deeper understanding of the immune characteristics of the tumor microenvironment and identify meaningful immune targets. In terms of immune targets, COAD has been poorly explored; thus, in the current study, based on the immune cell infiltration score and differentially expressed genes, COAD tumors were classified into hot and cold tumors. The Least Absolute Shrinkage and Selection Operator (LASSO) regression analysis was used to identify hub genes, construct a prognostic model, and screen potential immune targets. In total, 12 genes (CLK3, CYSLTR2, GJA10, CYP4Z1, FAM185A, LINC00324, EEF1A1P34, EEF1B2P8, PTCSC3, MIR6780A, LINC01666, and RNU6.661P) differentially expressed between hot and cold tumors were screened out. Among them, CYSLTR2 was considered as a potential candidate gene, because it showed a significant positive correlation with immune cell infiltration and immune checkpoints (PDCD1, CD274, and CTLA4). Finally, we constructed and validated a new prognostic model for COAD showing 0.854 AUC for the ROC curve, and these results provide sufficient potential to choose CYSLTR2 as an important immune target for the prognosis of COAD.

Functional Analysis of Autoantibody Signatures in Rheumatoid Arthritis

For the identification of antigenic protein biomarkers for rheumatoid arthritis (RA), we conducted IgG profiling on high density protein microarrays. Plasma IgG of 96 human samples (healthy controls, osteoarthritis, seropositive and seronegative RA, n = 24 each) and time-series plasma of a pristane-induced arthritis (PIA) rat model (n = 24 total) were probed on AIT’s 16k protein microarray. To investigate the analogy of underlying disease pathways, differential reactivity analysis was conducted. A total of n = 602 differentially reactive antigens (DIRAGs) at a significance cutoff of p < 0.05 were identified between seropositive and seronegative RA for the human samples. Correlation with the clinical disease activity index revealed an inverse correlation of antibodies against self-proteins found in pathways relevant for antigen presentation and immune regulation. The PIA model showed n = 1291 significant DIRAGs within acute disease. Significant DIRAGs for (I) seropositive, (II) seronegative and (III) PIA were subjected to the Reactome pathway browser which also revealed pathways relevant for antigen presentation and immune regulation; of these, seven overlapping pathways had high significance. We therefore conclude that the PIA model reflects the biological similarities of the disease pathogenesis. Our data show that protein array analysis can elucidate biological differences and pathways relevant in disease as well be a useful additional layer of omics information.

Monocytes and Macrophages in Spondyloarthritis: Functional Roles and Effects of Current Therapies

Spondyloarthritis (SpA) is a family of chronic inflammatory diseases, being the most prevalent ankylosing spondylitis (AS) and psoriatic arthritis (PsA). These diseases share genetic, clinical and immunological features, such as the implication of human leukocyte antigen (HLA) class I molecule 27 (HLA-B27), the inflammation of peripheral, spine and sacroiliac joints and the presence of extra-articular manifestations (psoriasis, anterior uveitis, enthesitis and inflammatory bowel disease). Monocytes and macrophages are essential cells of the innate immune system and are the first line of defence against external agents. In rheumatic diseases including SpA, the frequency and phenotypic and functional characteristics of both cell types are deregulated and are involved in the pathogenesis of these diseases. In fact, monocytes and macrophages play key roles in the inflammatory processes characteristics of SpA. The aim of this review is analysing the characteristics and functional roles of monocytes and macrophages in these diseases, as well as the impact of different current therapies on these cell types.

Primary Thyroid NUT Carcinoma With High PD-L1 Expression and Novel Massive IGKV Gene Fusions: A Case Report With Treatment Implications and Literature Review

Background

Nuclear protein in testis (NUT) carcinoma (NC) is a rare and aggressive undifferentiated carcinoma that typically arises from midline supradiaphragmatic structures. It is uniquely driven by a NUT gene rearrangement on chromosome 15q14. Few thyroid NCs have been reported and there are no established treatment guidelines for NUT carcinoma.

Method

Ultrasound-guided fine needle aspiration smear was performed for the preoperative diagnosis of thyroid lesions. Cytopathology, histology, and immunochemical staining all indicated NC. Fluorescence in situ hybridization (FISH), qRT-PCR, and next-generation sequencing (NGS) were used to analyze the genetic characteristics of NC.

Results

We describe a rare case of thyrogenic NC in a 38-year-old male with cytological, histological, immunohistochemical, and genetic features. Cytological smears and histopathological specimens showed typical features of NC. Immunohistochemistry confirmed strong immunoreactivity with NUT, EMA, P63, TTF-1, and c-myc. CK19 was positive exclusively in sudden keratosis. No immunoreactivity was found for neuroendocrine markers. FISH was applied to isolate the NUT gene on chromosome 15q14. The NGS results revealed a BRD4-NUT gene fusion, which was further confirmed by RT-qPCR. Structural variation (SV) of NUTM1 occurred in the exon region, and the mutation site was 15q14. Moreover, BRD4 single-nucleotide variation (SNV) occurs in the 3′ UTR at mutation site 19p13.12. The PD-L1 combined predictive score was over 30%. The patient received chemotherapy, followed by programmed cell death 1 (PD-1) inhibition with camrelizumab, and died 10 months after surgery.

Conclusion

Thyroid NC is an extremely rare and fatal malignant tumor. It is necessary to consider NC when squamous differentiation is observed cytologically or histologically. NGS is an effective tool for obtaining the final diagnosis and obtaining a better understanding of tumor pathogenesis. A large number of IGKV gene fusions in addition to the BRD4-NUT fusion may play a role in the pathogenesis and immunotherapy response of NC. Immunotherapy for NC remains to be explored due to the rarity of this aggressive malignancy.

Autoimmunity in psoriatic arthritis: pathophysiological and clinical aspects

Psoriatic arthritis (PsA) is an underdiagnosed entity with a broad impact on the quality of life. Although the pathogenesis is largely unknown, autoimmune footprints of the inflammation in PsA have increasingly been recognized. Most of the genetic variation predisposing to PsA is mapped to the class I major histocompatibility complex (MHC) region and shared by a variety of autoimmune diseases. Polymorphisms in the genes IL12B, IL23R, IL13, TNIP1, TRAF3IP2, TYK2, and many others explain the non- HLA genetic risk with little known functional consequences. Entheseal and synovial cellular infiltrate with oligoclonal CD8+ T cells and occasional germinal centers, loss of regulatory T cell function, and specific autoantibodies such as anti-PsA peptide, anti-LL-37, and anti-ADAMTSL5 are the immunopathological findings suggestive of autoimmunity. These were supported by clinical observations of autoimmune multimorbidity and treatment response to calcineurin/mTOR and co-stimulation inhibition.

Kinetics of Abacavir-Induced Remodelling of the Major Histocompatibility Complex Class I Peptide Repertoire

Abacavir hypersensitivity syndrome can occur in individuals expressing the HLA-B*57:01 major histocompatibility complex class I allotype when utilising the drug abacavir as a part of their anti-retroviral regimen. The drug is known to bind within the HLA-B*57:01 antigen binding cleft, leading to the selection of novel self-peptide ligands, thus provoking life-threatening immune responses. However, the sub-cellular location of abacavir binding and the mechanics of altered peptide selection are not well understood. Here, we probed the impact of abacavir on the assembly of HLA-B*57:01 peptide complexes. We show that whilst abacavir had minimal impact on the maturation or average stability of HLA-B*57:01 molecules, abacavir was able to differentially enhance the formation, selectively decrease the dissociation, and alter tapasin loading dependency of certain HLA-B*57:01-peptide complexes. Our data reveals a spectrum of abacavir mediated effects on the immunopeptidome which reconciles the heterogeneous functional T cell data reported in the literature.

HLA-Cw1 and Psoriasis

Psoriasis is a chronic inflammatory skin condition with regional and ethnic differences in its prevalence and clinical manifestations. Human leukocyte antigen (HLA)-Cw6 is the disease allele conferring the greatest risk to psoriasis, but its prevalence is lower in Asian individuals. Recent studies have found associations between HLA-Cw1 and some Asian populations with psoriasis, especially Southern Chinese. HLA-Cw6 was associated with type I early-onset psoriasis, guttate psoriasis, Koebner phenomenon, and better response to methotrexate, interleukin (IL)-12/23, IL-17, and IL-23 targeting drugs. In contrast, HLA-Cw1 positivity has been associated with erythrodermic psoriasis, pustular psoriasis, and the axial type of psoriatic arthritis. Furthermore, HLA-Cw1 was more frequently associated with high-need patients who did not respond to conventional therapies. No known trigger factor nor autoantigen has been identified for HLA-Cw1 positivity. However, HLA-Cw1 has been linked to some viral agents. For example, cytotoxic T lymphocytes recognize multiple cytomegalovirus pp65-derived epitopes presented by HLA alleles, including HLA-C*01:02. In addition, cytomegalovirus can lead to severe exacerbation of psoriatic skin disease. The proposed interaction between viral infection, HLA-Cw1, and psoriasis is through the killer cell immunoglobulin-like receptors of natural killer cells. Given the diverse nature of psoriasis pathogenesis and the difference in HLA-Cw prevalence in different racial groups, more studies are needed to confirm the role of HLA-Cw1 in psoriasis.

Identification of the novel HLA-C*03:04:79 allele in a Chinese bone marrow donor

HLA-C*03:04:79 differs from HLA-C*03:04:01:01 by a single nucleotide substitution at position 834 G>A.

Is gene expression among women with rheumatoid arthritis dysregulated during a postpartum flare?

Background

To evaluate our hypotheses that, when rheumatoid arthritis (RA) flares postpartum, gene expression patterns are altered compared to (a) healthy women, (b) RA women whose disease activity is low or in remission postpartum, and (c) pre-pregnancy expression profiles.

Methods

Twelve women with RA and five healthy women were included in this pilot study. RA disease activity and postpartum flare were assessed using the Clinical Disease Activity Index (CDAI). Total RNA from frozen whole blood was used for RNA sequencing. Differential gene expression within the same women (within-group) over time, i.e., postpartum vs. third trimester (T3) or pre-pregnancy (T0), were examined, using a significance threshold of q < 0.05 and fold-change ≥ 2.

Results

Nine of the women with RA experienced a flare postpartum (RA_Flare), while three had low disease activity or were in remission (RA_NoFlare) during that time frame. Numerous immune-related genes were differentially expressed postpartum (vs. T3) during a flare. Fold-changes in expression from T3 to postpartum were mostly comparable between the RA_Flare and healthy groups. At 3 months postpartum, compared to healthy women, several genes were significantly differentially expressed only among the RA_Flare women, and not among the RA_NoFlare women. Some of these genes were among those whose “normal” expression was significantly modulated postpartum, and the postpartum expression patterns were significantly altered during the RA flare. There were also some genes that were significantly differentially expressed in RA_Flare compared to both healthy and RA_NoFlare women, even though their expression was not significantly modulated postpartum. Furthermore, while postpartum expression profiles were similar to those at pre-pregnancy among healthy women, significant differences were found between those time points among the RA_Flare women.

Conclusions

The large majority of gene expression changes between T3 and 3 months postpartum among RA women who flared postpartum reflected normal postpartum changes also seen among healthy women. Nonetheless, during a postpartum flare, a set of immune-related genes showed dysregulated expression compared to healthy women and women with RA whose disease activity was low or in remission during the same time frame, while other genes demonstrated significant differences in expression compared to RA pre-pregnancy levels.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13075-021-02418-w.