- 1Division of Rheumatology, Mayo Clinic, Rochester, MN, United States
- 2Division of Internal Medicine, Dupuytren University Hospital, Limoges, France
- 3INSERM U1308, Faculty of Medicine, University of Limoges, Limoges, France
- 4Stem Cell Genomics Core, Stem Cell Program, University of California, San Diego, La Jolla, CA, United States
- 5Next Generation Sequencing Core, Salk Institute for Biological Studies, La Jolla, CA, United States
- 6Division of Internal Medicine, Nîmes University Hospital, University of Montpellier, Nîmes, France
- 7Division of Gastroenterology, Department of Medicine, University of California, San Diego, San Diego, CA, United States
- 8Division of Laboratory Medicine and Pathology, Mayo Clinic, Rochester, MN, United States
- 9Cell Biology, Dupuytren University Hospital, Limoges, France
Objective: To identify the key coding genes underlying the biomarkers and pathways associated with giant cell arteritis (GCA), we performed an in situ spatial profiling of molecules involved in the temporal arteries of GCA patients and controls. Furthermore, we performed pharmacogenomic network analysis to identify potential treatment targets.
Methods: Using human formalin-fixed paraffin-embedded temporal artery biopsy samples (GCA, n = 9; controls, n = 7), we performed a whole transcriptome analysis using the NanoString GeoMx Digital Spatial Profiler. In total, 59 regions of interest were selected in the intima, media, adventitia, and perivascular adipose tissue (PVAT). Differentially expressed genes (DEGs) (fold-change > 2 or < −2, p-adjusted < 0.01) were compared across each layer to build a spatial and pharmacogenomic network and to explore the pathophysiological mechanisms of GCA.
Results: Most of the transcriptome (12,076 genes) was upregulated in GCA arteries, compared to control arteries. Among the screened genes, 282, 227, 40, and 5 DEGs were identified in the intima, media, adventitia, and PVAT, respectively. Genes involved in the immune process and vascular remodeling were upregulated within GCA temporal arteries but differed across the arterial layers. The immune-related functions and vascular remodeling were limited to the intima and media.
Conclusion: This study is the first to perform an in situ spatial profiling characterization of the molecules involved in GCA. The pharmacogenomic network analysis identified potential target genes for approved and novel immunotherapies.
Introduction
Giant cell arteritis (GCA), also called temporal arteritis or Horton’s disease, is the most common systemic large-vessel vasculitis; it primarily affects the aorta and its major branches in individuals aged over the age of 50s (1). This inflammatory disease is related to immune system aging (2) and a high risk of arterial occlusion, which can lead to irreversible blindness, strokes, and death (3, 4). Despite advances in understanding the underlying disease mechanisms of GCA over the past three decades, the treatment options are still limited. Furthermore, a high proportion of patients relapse during treatment with corticosteroids, methotrexate, or tocilizumab, which are the current recommended treatments for this vasculitis (5–8). In almost half of the patients, arterial inflammation persists even 1 year after corticosteroid treatment (9).
The management of GCA is challenging owing to its histopathological complexity, which occurs because it involves three arterial layers, from outer to inner: adventitia, media, and intima. Despite the potential significance of perivascular adipose tissue (PVAT) in the development of GCA, its role has not been thoroughly investigated (10). The disease process begins in the outer layer of the vessel, the adventitia, where resident dendritic cells are stimulated by an unknown trigger. This initiates an inflammatory response that involves infiltration by lymphocytes, macrophages, and giant cells, which damage the elastic tissue through release of metalloproteases and reactive oxygen species. As a result, vascular smooth muscle cells differentiate into myofibroblasts, migrate, and eventually cause intimal hyperplasia, leading to vessel obstruction (11).
Although the characterization of GCA through immunoprofiling of arterial tissue sections is essential for understanding immune cell interactions and evaluating the effects of immunosuppressive and anti-inflammatory drugs, previous studies had some limitations such as tissue heterogeneity or small sample sizes (12–15). Multiple studies have analyzed the transcriptomic profile of the entire artery without considering the anatomical characteristics of GCA lesions. As a result, transcriptional events specific to a particular arterial layer may have been missed when the entire artery was analyzed. Therefore, we used NanoString GeoMx® Digital Spatial Profiler (DSP) technology to explore the molecular events involved in inflammatory and vascular remodeling in each sub-layer of GCA (16). Thanks to this technology that allows complex RNA profiling and quantification of the transcripts in specific regions within a tissue, we were able to study spatially the transcriptome of GCA and non-GCA areas of human temporal arteries using formalin-fixed, paraffin-embedded tissues. Whole tissue spatial transcriptomic analyses can better elucidate functional mechanisms and gene pathways in pathological tissues with complex organization, such as human temporal artery, which is the main artery involved in patients with GCA.
The present study is the first to provide in situ spatial transcriptomic analyses for the identification of layer-specific genes and cellular functions related to GCA, potentially leading to the identification of novel therapeutic targets. Here, we proposed a model in which spatiotemporal transcriptomic analysis can be used to identify previously unknown pathways and new drug targets for GCA. Therefore, our results can be useful for drug development strategies and may improve patient outcomes.
Patients and methods
Study design and participants
In this retrospective observational study, we recruited non-consecutive patients referred to the internal medicine department of the Limoges University Hospital Dupuytren, France, between January 1, 2014, and December 31, 2018. The participants underwent temporal artery biopsy (TAB) for suspected GCA, which was diagnosed based on the 2022 ACR/EULAR criteria (17) and histologically confirmed based on the presence of panarteritis with an inflammatory infiltrate (e.g., macrophages and lymphocytes) and fragmentation of the internal elastic lamella, corresponding to the most common histological pattern described by Hernández-Rodríguez et al. (18, 19). This study included patients who were treated with glucocorticoids for 1–10 days before TAB as this treatment does not affect the histopathological findings of TAB (20, 21). Participants with a negative TAB and a final diagnosis other than GCA were considered controls. We recorded the age, sex, symptoms, treatments, and outcomes of each participant at the time of diagnosis.
Tissue collection and preparation
Unilateral TABs were performed under local anesthesia by an ophthalmologic surgeon using standard techniques (22). A 1–4-cm-long specimen of the temporal artery was obtained for histological analysis. The temporal artery specimens were analyzed using standard protocols. The specimens were subjected to gross evaluation and were cross sectioned serially. Tissues were fixed in 10% neutral buffered formalin for 24 h, followed by transfer to 70% ethanol before paraffin embedding. For each patient and control, the most representative segments were grouped into formalin-fixed paraffin-embedded (FFPE) temporal artery blocks using the tissue microarray method. Tissue microarrays were constructed from 2-mm-diameter cores punched from each FFPE tissue block.
NanoString GeoMx Digital Spatial Profiler
The NanoString GeoMx DSP platform quantifies RNA abundance by counting unique indexing oligonucleotides assigned to each target of interest. Those oligonucleotides are covalently attached to mRNA hybridization probes with an ultraviolet-photocleavable linker (16). First, FFPE temporal artery blocks from GCA and control patients were sectioned to a thickness of 5 μm. The sections were mounted on an adhesive or positively charged slides (cat# S21.2113.A; Bond Plus slides; Leica Biosystems, Richmond, VA, USA) with a maximum of 24 individual tissues per slide at the San Diego Pathology Department (San Diego, CA, USA). Some slides were stored at 4°C, while others were deparaffinized and used for hematoxylin & eosin (H&E) staining (San Diego Pathology Department) for subsequent microscopy imaging (Keyence BZ-X800). A cardiovascular pathologist (M.C.B.) and two GCA specialists (S.P. and C.M.W.) reviewed all H&E-stained tissues. Overall, one GCA patient was classified as healed and excluded from this study because the artery exhibited no inflammatory infiltrate, and one control sample was also excluded because it showed tertiary lymphoid structure in the adventitia.
The freshly cut sections were stained with fluorescently labeled imaging reagents and mRNA markers as described previously (16). Briefly, sections were labeled with the cell markers, α-smooth muscle actin, CD3, and CD68, and the nuclear stain SYTO 13. The NanoString commercial Human Whole Transcriptome Atlas (GeoMx WTA) panel (18,000 protein-coding genes) was selected. Four slides (two with GCA patients and two with controls) were scanned on a GeoMx DSP instrument (NanoString Technologies Inc., Seattle, WA, USA) and individual regions of interest (ROIs) with a maximum diameter of 400 μm were created within each of the four arterial layers in GCA and controls.
In this study, we considered a sample as a control where no inflammatory infiltrate was observed, and the 4 artery layers (intima, media, adventitia, and PVAT) were distinguished on the basis of H&E staining as well as internal elastic lamina between intima and media; concentric smooth muscle cells in the media unlike the adventitia, and, finally, the presence of adipocytes in the PVAT unlike the adventitia. A representative portion of each layer was selected by two GCA specialists (S.P. and R.G.) and confirmed by a pathologist (M.C.B.) (Supplemental Figure 1). Once each ROI was compartmentalized, ultraviolet-cleaved indexing oligonucleotides were collected into a 96-well plate. Libraries were prepared according to the manufacturer’s instructions (protocol 01/2019). The library quality and quantity were assessed using the Bioanalyzer DNA High Sensitivity Analysis. Samples were sequenced on the NovaSeq 6000 platform (Illumina, San Diego, CA, USA) and reads were digitally quantified and normalized using GeoMx DSP Data Analysis software (Merritt).
Data processing and analysis
Data obtained from the NanoString GeoMx DSP platform were analyzed using a dedicated software (GeoMx Analysis Suite 2.3). The quality control analysis showed that all ROIs had raw read counts above 1 million and good alignment rate and sequencing saturation. The local outliers were negative control probes and were removed from the ROIs. ROIs with deduplicated probe counts Q3 (upper quartile) < 2 were excluded. Furthermore, ROIs were removed if < 1% of their genes had read counts higher than the limit of quantification (LOQ). LOQ was defined as geomean (NegProbe) × geoSD(NegProbe)² for each ROI. Target filtering was applied to retain gene targets with read counts above LOQ in at least 10% of ROIs. Q3 normalization was applied on the filtered ROIs and gene targets.
The raw data were deposited at the NCBI Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE237441) under the accession number GSE237441.
Differential gene expression across groups was analyzed using linear mixed effect models; differentially expressed genes (DEGs) were defined as fold-change > 2 or < −2, and Benjamini-Hochberg-adjusted p < 0.01. Analyses were conducted using R software (version 3.2.2; R Foundation for Statistical Computing, Vienna, Austria).
The upregulated and downregulated DEGs were analyzed using DAVID Bioinformatics Resources software (version 2021; https://david.ncifcrf.gov/) for functional annotations and gene set enrichment analysis. DAVID was used to obtain statistically enriched Gene Ontology (GO) terms for biological processes (GOTERM_BP_FAT). GO terms with EASE score < 0.05 were considered statistically enriched. Based on the known mechanisms in GCA, pathways of interest, signaling cascades, and molecules were selected from among the dysregulated GO terms. These selections are presented as dot plots using ggplot2 package of R software (version 3.2.2).
Dysregulated genes per layer with a fold-change > 2 or < −2 were organized according to the STRING network (https://string-db.org/). The one or two main networks were represented for each layer. Genes disconnected from the main network were removed. An interaction score with the highest confidence (0.900) was used. The interaction sources excluded text mining. Using networks generated by the STRING network, a pharmacogenomic network was constructed by targeting existing druggable genes (https://www.dgidb.org/). Drugs with Food and Drug Administration (FDA) approval were selected.
Fluorescence immunohistochemistry
From block, 4-μm-thick slices were cut and dewaxed by successive immersion in pure xylene (Sigma). Rehydration and unmasking of antigens were carried out by sequentially immersion into 100%, 90%, 80% and 70% ethanol baths (Sigma), phosphate-buffered saline (PBS; Gibco), and citrate buffer (pH 7). Saturation of nonspecific antigenic sites was achieved by incubation in PBS containing 3% bovine serum albumin (BSA). The primary antibodies were diluted in 3% BSA and incubated overnight at 4°C. The antibodies used were as follows: anti-CD74 (1:50, mouse, Abcam), anti-CD90 (1:50, rabbit, Abcam), anti-CD3 (dilution recommended by the supplier, rabbit, Ventana). The sections were subsequently incubated with a secondary antibody coupled to Alexa Fluor 594 or 488 (Molecular Probes) for 60 min. After washing with PBS, nuclear counterstaining was carried out with 4′,6-diamidino- 2-phenylindole (Sigma) and the sections were mounted in an aqueous medium (Vectashield, Sigma). Images were acquired with Leica microscope (DMi8, Leica).
Role of the funding source
The funders of the study had no role in study design, data collection, analysis, or interpretation, or in writing the manuscript.
Ethics
The study protocol (87RI21_0065) was approved by the Ethics Committee of the University Hospital of Limoges, France (23). Temporal artery tissues were obtained from the biological collection DC-2010-1079. All participants provided informed consent.
Results
Participant characteristics
The characteristics of GCA patients (n = 9) and controls (n = 7) are presented in the Supplemental Table 1. Among the controls, the final diagnoses were non-arteritic ischemic visual impairment (n = 4), peripheral seronegative polyarthritis (no argument for polymyalgia rheumatica) (n = 2), and infection-related encephalopathy (n = 1). Three controls had an inflammatory syndrome related to an acute infection. Among the entire cohort, one patient had diabetes (control), one had a previous cardiovascular complication (control), and two patients were active smokers (one control and one GCA). Corticoids were used before TAB for 8/9 GCA patients (mean 4.3 days before TAB, mean dose of 0.85 mg/kg/d) and 3/7 controls (mean 18.3 days before TAB, mean dose of 0.7 mg/kg/day).
All GCA patients had panarteritis involving the three arterial layers with destruction of the internal lamella on histological examination. The histological sections of the artery are presented in Supplemental Figure 1. Then temporal arteries from the 16 participants were subjected to transcriptional analysis for each arterial layer (intima, media, and adventitia) and the PVAT following technical pipeline indicated in Figure 1. Because the control arteries did not present CD3 and CD68 staining, layers were selected under the control of a senior pathologist depending on their macro characteristics. Four ROIs involving the media, adventitia, and PVAT were excluded from analysis due to poor quality, respectively. Altogether, we analyzed 59 ROIs, checked for quality and correspondence to each layer, from 16 participants to conduct our spatial transcriptomic study.
Figure 1 Schematic summary of the protocol study. Tissue microarrays (TMA) are prepared from a cohort of human temporal artery biopsies. Freshly cut sections of TMA are run through NanoString GeoMx DSP instrument. Synthetic DNA oligonucleotide barcodes are attached to in situ probes for mRNA detection via a UV photocleavable linker. A single tissue section is labeled with a cocktail of these probes and fluorescently labeled morphology markers. Regions of interest (ROI) for molecular profiling are selected based on the fluorescent images, and UV light illuminates each region independently, which releases the barcodes from that region for subsequent counting using Next-Generation sequencing.
Temporal artery whole transcriptome analysis exhibits different profiles depending on layers in giant cell arteritis
As attempted, correlation analysis demonstrated that disease status (GCA and control) only differentiated samples (Figure 2A). Interestingly, Pearson correlation coefficient showed that GCA transcriptome remains more homogeneous than in controls as illustrated by increasing intensity of blue reflecting correlation index tending toward 1.
Figure 2 Temporal Artery Whole Transcriptome Analysis across layers. (A) Hierarchical clustering according to the correlation between the 59 ROIs studied with identification of status (Giant Cell Arteritis or Control), gender (Male or Female) and layers (Intima, Media, Adventitia, Perivascular adipose tissue). (B) A total of 337 differentially expressed genes (DEGs) were identified for the intima, media, and/or adventitia of GCA compared with controls (fold change < -2 or > 2, adjusted p-value < 0.01). Some genes are only dysregulated in one layer and others are shared by two or three arterial layers. (C) Hierarchical clustering of patients by DEGs and principal component analysis on gene expression profiles (12,076 genes in total) from temporal artery biopsy samples in both patient groups (Giant Cell Arteritis and control temporal arteries) for each different layer. GCA, giant cell arteritis; F, female; M, male; PVAT, perivascular adipose tissue; DEGs, differential expressed genes; PC, principal component.
A total of 12,076 genes were detected and 337 different genes (2.8%) were dysregulated between GCA and control participants (fold-change < −2 or > 2, p-adjusted < 0.01). Venn diagram presents unique and common DEGs for the intima, media, and adventitia (Figure 2B). The number of DEGs increased from external layer, adventitia with 40 genes, to media and intima with respectively 227 and 282 genes. In total, 32 common DEGs were significantly identified among the three layers (intima, media and adventitia), corresponding to the major histocompatibility complex (CD74, B2M, HLA-DRB1, HLA-A, HLA-B, HLA-C, and HLA-DRA), macrophages (CD68), immunoglobulins (IGKC and IGHG1-4), apolipoprotein (APOC1), cathepsins (CTSB and CTSZ), and tissue remodeling (ADAM15, COL1A2, COL3A1, TMSB10, and ARHGDIB). The other common genes were NPIPB6, PABPC1, and those associated with ribosomes (RPL23, RPL28, RPL31, RPL37, RPS12, RPS28, RPS29, RPS3, and RPS4X). Patterns for the four layers in each participant are shown in Figure 2C. Principal component analysis (PCA) highlighted that only intima present a distribution profile more homogeneous and reproducible than other GCA layers and control samples.
Altogether, these results highlight that the intima is the layer where the transcriptome is the most dysregulated depicting a gradient of decreased DEG from the internal to the periphery of the artery.
Top differentially expressed genes across layers in giant cell arteritis compared to control temporal arteries
Overall, most DEGs were upregulated in ROIs from GCA compared to controls. The top 20 upregulated and top 5 downregulated genes per arterial layer in GCA patients compared to the controls are shown in Figure 3A. The highest numbers of DEGs were observed in the media and intima.
Figure 3 Top differentially expressed genes (DEGs) across layers in Giant Cell Arteritis compared to control temporal arteries. (A) Top 20 upregulated genes and top 5 downregulated genes for each layer (intima, media, adventitia, and perivascular adipose tissue) of Giant Cell Arteritis temporal arteries. (B) Hierarchical clustering among each ROI participant for the top 30 dysregulated genes per layer. Up, up-regulated genes; ns, non statistically significative dysregulated genes; Down, down-regulated genes; GCA, giant cell arteritis; PVAT, perivascular adipose tissue.
Figure 3B presents the top 30 dysregulated genes across arterial layers between both studied groups. In the media and intima, CD74 was the most upregulated gene between GCA and control specimens. Healthy temporal arteries do not exhibit CD74 staining on immunofluorescence analysis (Supplemental Figure 2). In contrary, GCA arteries has CD74+ cells, particularly within the intima, close to the media border, as well as in the media and the adventitia. Also, cells expressing CD74 are also CD90+; however, they do not express CD3, marker for lymphocytes. Finally, migration inhibitory factor (MIF), a receptor for CD74, was significantly overexpressed in the intima of GCA compared to controls (p-adj = 0.0031).
Altogether, these results unveiled that most dysregulated genes in the three arterial layers belong to the antigen presentation gene family and vascular remodeling functions.
Immune related functions and vascular remodeling pathways are predominant in media and intima
We classified the top 30 DEGs by arterial layers according to their main functions (Figure 4A). Most top DEGs per layer were associated with immune or remodeling functions. Immune function-related genes were associated with macrophages (CD68 and cathepsins), major histocompatibility complex, complement system, and immunoglobulins. Remodeling-related genes were mainly associated with metalloproteases, reactive oxygen species, cytoskeleton organization, and collagen.
Figure 4 Functions and pathways related top DEG in Giant Cell Arteritis temporal arteries. (A) Top 30 DEG per layer (intima, media, adventitia) categorized by immune and remodeling functions (B) Dot plot showing significant remodeling related pathways according to layer (intima, media, adventitia). MHC, major histocompatibility complex; MMP, metalloproteinase; ROS, reactive oxygen species.
Using the DAVID tool, remodeling-related pathways of interest with gene dysregulation were selected (Figure 4B). PVAT was excluded due to the low number of dysregulated genes between GCA and control specimens. Remodeling pathways associated with GCA primarily involve the intima and media and include the reactive oxygen species and collagen metabolic process, fibroblast proliferation, muscle cell differentiation, ossification, and angiogenesis. The dysregulated genes for each signaling pathway are detailed in Supplemental Figure 3.
Altogether these results suggest that activated immune cells are actively and mainly present in the intima and the media layers. Likewise, vascular remodeling process following a gradient from the outside to inside part of the vessel.
Meta-analysis studies similar expression profiles among genes of interest across arterial layers
The DEGs among the four layers were compared (Figure 5A). The meta-analysis identified 412 upregulated genes and 101 downregulated genes. Genes whose expression profiles were similar to CD74, the top DEG, are shown in Figure 5B. These genes were mainly involved in antigen presentation, including the major histocompatibility complex (B2M, HLA-B, HLA-DRA, and HLA-DRB1), cathepsin (CTSB), ferritin (FTH1), immunoglobulin (IGHG2), and lysozyme (LYZ).
Figure 5 Meta-analysis across layers of Giant Cell Arteritis temporal arteries. (A) Hierarchical clustering of 513 dysregulated genes in giant cell arteritis (GCA) across layers. (B–F) Genes presenting the same expression profiles across layer for: (B) CD74, the most dysregulated gene in GCA compared to control; (C) CD68, macrophage marker; the two main genes dysregulated and implicated in vascular remodeling, (D) MMP2, a matrix metalloproteinase type 2 marker, (E) COL3A1, collagen type III alpha-1 marker. PVAT, perivascular adipose tissue; DEG, differential expressed gene; FC, fold-change.
As macrophages are the key immune cells in GCA, we presented the genes with a similar expression profile within the four layers like CD68, a marker for macrophages (Figure 5C). The DEGs included those related to cytokines (TNF), chemokines (CCL5, CXCL9, and CXCR4), YKL40 (CHI3L1), metalloproteases (MMP2 and MMP9) secreted by macrophages, membrane receptors (CD63 and CD206), and iron metabolism-related molecules (FTH1 and SOD2). Similarly, we evaluated vascular remodeling during GCA based on the two mains associated DEGs: MMP2 (Figure 5D) and collagen COL3A1 (Figure 5E). We identified the top genes with the same expression profile across the layers.
Potential novel drug-targeted pathways for GCA treatment
Transcriptomic profiling is an efficient and essential tool for drug target discovery. Several studies have combined gene expression profiling with systematic approaches for drug re-repurposing or biomarker discovery (24, 25). Among strategies for novel drug target identification, the network-based framework uses curated biological network topology to investigate the association across genes, diseases, and drugs (26). We constructed a pharmacogenomic network by integrating DEGs from each artery layer, a high-confidence protein-protein interaction network, and drug-gene interactions. Then, the top dysregulated genes are presented as a STRING analysis to highlight their relationships. The significantly connected genes are shown in Figure 6. Networks show the strong relationships between several genes. Two mains separate networks exist for each layer corresponding to global immune function (yellow) and vascular remodeling (blue). Key genes are present at nodes of those interactions and connected to other genes.
Figure 6 Pharmacogenomic network analysis uncovers drug-gene interaction landscape in Giant Cell Arteritis. The two largest connected components in the DEG-based protein-protein interaction network are composed of 43 up-regulated genes (nodes) in intima (A) and 42 up-regulated genes in media (B). Using pharmacological information, genes were identified as druggable with FDA-approved, pharmaceutical drugs. Pharmaceutical drugs are divided into 2 types: specific target of one layer and intima-media common drugs. GCA, giant cell arteritis.
Approved therapies targeting these genes are presented, of which have not been evaluated in relation to GCA previously. The spatial approach highlights therapies targeting top DEGs common to the intima and the media (red). Therapies targeting a gene only present in one network layer are also represented (black).
Discussion
In the present study, we performed a spatial transcriptome-wide expression analysis to identify key genes and biomolecular processes involved in temporal arteritis. By comparing gene expression patterns within the different layers of the temporal arteries among GCA and non-GCA patients, we saw a significant upregulation - rather than down-regulation - of genes in all arterial layers from GCA in comparison to controls. This could suggest some insufficient regulatory mechanisms occurring during GCA, and overriding of regulatory functions. Intima and media layers were particularly affected by gene expression dysregulation in GCA in comparison to the other layers, as well as compared to controls. Differences between gene expressions were less obvious in the adventitia, which may be due to the low number of inflammatory cells in the adventitia from the two TABs of GCA patients (GCA#4 and GCA#5). Furthermore, the initial mechanism of GCA involves adventitial damage and stimulation of resident dendritic cells via toll-like receptors by an unknown stimulus (27). As TABs are often performed at an advanced stage of GCA when vascular remodeling has already occurred, the adventitia may be less affected by dysregulation than the media and intima at this stage. The PVAT did not show major transcriptomic differences between the groups which might be due to the low number of cells present within the ROIs compared to the other layers. Although the perivascular layer has been poorly studied in GCA, it may be involved in the inflammation of the adjacent artery like it has been shown in other vascular diseases (28, 29). A previous study has suggested the involvement of the periadventitial tissue in a subgroup of GCA (30), however, in our study, all the GCA patients did not exhibit the histological pattern they described.
CD74 was the most dysregulated gene across layers in GCA. This gene encodes a cell transmembrane receptor of the macrophage migration inhibitory factor (31), which is expressed on antigen presenting cells (dendritic cells, macrophages, and B lymphocytes) and non-immune cells (epithelial and endothelial cells). CD74 is known to be involved in apoptosis, immune response, and cell migration. Notably, CD74 plays a role in tissue repair and wound healing. For instance, CD74 deficiency protects against glomerulonephritis in lupus (32). Interestingly, CD74 is a ligand for the macrophage migration inhibitory factor (MIF), which is over-expressed in the serum of patients with GCA (33).
In the present study, most genes with the same pattern as CD74 expression across different layers were involved in the immune response and belonged to the major histocompatibility complex classes I (B2M, HLA-A, and HLA-B) and II (HLA-DRA and HLA-DRB1) (34).
Next, we focused on macrophages, as they are the main immune cells involved in GCA. Macrophages can form granulomas and fuse together to produce giant cells. Lymphocytes recruit and activate macrophages in arteries through secretion of cytokines and growth factors (IL-6, IL-17, TNF-α, IFN-γ, and GM-CSF). Two main populations of CD68+ macrophages have been identified in GCA: adventitial macrophages, which are involved in TGF-β1, IL-6, and IL-1β secretion, and intimal and medial macrophages, which are involved in reactive oxygen species, nitric oxide, and metalloproteases secretion (13). They can also produce PDGF and VEGF. Under the combined action of MMP, PDGF, and VEGF, the internal elastic lamella of the media is destroyed, and vascular remodeling occurs. In the present study, the main DEGs with the same expression pattern as CD68 were related to antigen presentation. The FTH1 and SOD2 genes, related to M1 macrophage polarization, exhibit the same expression pattern as CD68 among the arterial layers, as compared to CD163, CD206, and CHI3L1, which are linked to M2 macrophages (35). Previous studies have demonstrated heterogeneity among CD68 macrophages in GCA temporal arteries (36, 37). Interestingly, the expression profiles of CTSB and TMSB10 are similar to the one of CD68, which is related to the promotion of M2 macrophages (38, 39). On the other hand, FTH1 encodes the major intracellular iron storage protein (ferritin heavy chain 1) and is expressed abundantly by M1 macrophages (35). We note that lymphocytic activation was less marked than that of macrophages. As most of the GCA patients in our study had glucocorticoids before TAB, the Th17 population was probably reduced, which could explain a small lymphocyte activation as demonstrated by Deng et al. (40).
Vascular remodeling plays a key role in GCA. Intimal hyperplasia is the main cause of arterial occlusion and leads to often-irreversible ischemic complications, which depend on vascular muscle cell differentiation and fibroblast proliferation, as described previously (24). No significant differences were observed between GCA patients and controls in terms of the expressions of IL-1 and IL-6 genes, although these cytokines are produced in the temporal arteries (19). These findings suggest that vascular remodeling has already started when TAB was performed. In this study, dysregulated remodeling is mainly correlated with the intima and media layers. The metalloproteinase genes MMP2 and MMP9 were overexpressed, which was also shown in previous studies (41). The release of reactive oxygen species by macrophages leads to lipid peroxidation of the membrane phospholipids of vascular smooth muscle cells, resulting in their combination with MMP and apoptosis (42). Furthermore, PDGF secreted by macrophages and giant cells participate in intimal hyperplasia (43). Neoangiogenesis comes from the action of VEGF secreted by macrophages. Our study demonstrated that angiogenesis occurs in the intima and media as previously published (43). Interestingly, ossification and calcification have never been thoroughly studied in GCA. Here, we demonstrated that these processes are present within the intima and the media, and are likely to play a role in the pathophysiology of GCA (44). Moreover, ossification cluster data are consistent with the overexpression of COL1A1, MMP9 and COL2A1, which have been previously found in GCA (45–47).
Finally, we built a pharmacogenomic network based on DEGs and Food and Drug Administration (FDA) approved drugs targeting these genes. By providing a system-wide view of mechanistic gene interactions, our spatial pharmacogenomic network can facilitate the design of novel pharmacologic intervention strategies. Our results suggest that macrophages are a promising target in GCA. In particular, granulocyte-macrophage colony-stimulating factor can be blocked by mavrilimumab, a growth factor for macrophage maturation (48, 49). Our pharmacogenomic network shows the potential benefit of blocking several macrophage-related genes using MMP2/MMP9- and CXCR4-targeting drugs. CD74 being, in this study, the top DEG common to the three arterial layers in GCA, this gene could be another potential drug target to consider. Interestingly, the anti-CD74 drug milatuzumab has been evaluated in hematological diseases with promising results so far (50).
The present study has several strengths. First, we used an innovative technique for human artery spatial analysis. Second, this is the first whole-transcriptome analysis of TABs in GCA using the NanoString GeoMx DSP. Third, we identified a large number of DEGs, even after multiple hypothesis correction, indicating that, despite the limited number of patients studied here, the gene expression differences between GCA and control temporal arteries are robust and reliable (fold-change < −2 or > 2 with p-adjusted value < 0.01). There were a few limitations to this study. First, the study had a relatively small sample size (9 GCA patients and 7 controls) compared to what it is commonly conducted in GCA research. Also, all GCA patients exhibited panarteritis but no other histological patterns were observed. Because of the sample size, clinical characteristics of patients (large-vessel vasculitis, cranial disease, relapse, or remission) could not be significantly analyzed. Some controls had an inflammatory background (2 peripheral seronegative arthritis and 1 infection related encephalopathy), but these 3 controls had a negative TAB with no inflammatory infiltrate. Finally, our study used FFPE samples from TAB specimens but did not analyze cells within circulation.
In conclusion, we uncovered the transcriptomic signatures of GCA arteries on a regional scale spatially according to the arterial layers for the first time using the GeoMx DSP. This novel technique allows to better understand GCA and identify new therapeutic targets. For that, we used a pharmacogenomic network-based approach to highlight the targeted genes and therapeutics for a potential multimodal treatment of GCA. We hope our findings will be useful for further whole-tissue and single-cell multiomics profiling studies investigating the biomolecular pathways, networks, and biomarkers in GCA.
Data availability statement
The data presented in the study are deposited in the NCBI GEO repository, accession number GSE237441.
Ethics statement
The studies involving humans were approved by Ethics Committee of the University Hospital of Limoges. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
SP and EM designed this study. EL collected data on clinical and patient characteristics. SP, MB, and CW reviewed the temporal artery sections. SP and RG selected the ROIs. FT and HA performed immunostaining. SP, EM, RG, and TN wrote the manuscript. EL, A-LF, M-OJ, CW, and K-HL reviewed the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This project was funded by NanoString GeoMx Whole Transcriptome Atlas Regional Grant (G-035 UCSD Molina). The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article, or the decision to submit it for publication.
Acknowledgments
The authors thank the study participants, and the UC San Diego Stem Cell Genomics Core, the Stem Cell Program, La Jolla, CA, and Cristian Quintero for his assistance with the Keyence microscope. We also thank San Diego Pathology and David Tobias Vargas for his technical assistance, as well as NanoString for the grant opportunity and the GeoMx DSP specialists in Seattle, WA. SP received a fellowship grant from the Servier Institute France.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1237986/full#supplementary-material
Supplementary Figure 1 | Sections of arteries with selected region of interest (intima, media, adventitia, perivascular adipose tissue) from each study participant: left hematoxylin eosin staining; right CD68 (yellow) CD4 (red) alpha-smooth muscle actin (green) SYTO 13 (blue) staining. GCA, giant cell arteritis; A, adventitia; M, media; I, intima; PVAT, perivascular adipose tissue.
Supplementary Figure 2 | Immunofluorescence analysis by microscopy of healthy and GCA temporal arteries. (A): In healthy artery, cells do not express CD74. (B): In GCA arteries, CD90+ cells express CD74, particularly in the outer part of the neointima near the media. CD3+ cells do not express CD74. GCA, giant cell arteritis, A, adventitia, M, media, NI, neointima.
Supplementary Figure 3 | Dysregulated genes associated with remodeling pathways across layers in Giant Cell Arteritis compared to control temporal arteries. GCA, giant cell arteritis. The English in this document has been checked by at least two professional editors, both native speakers of English. For a certificate, please see: http://www.textcheck.com/certificate/xvD6Vy.
References
1. Weyand CM, Goronzy JJ. Clinical practice. Giant-cell arteritis and polymyalgia rheumatica. N Engl J Med (2014) 371(1):50–7. doi: 10.1056/NEJMcp1214825
2. Gloor AD, Berry GJ, Goronzy JJ, Weyand CM. Age as a risk factor in vasculitis. Semin Immunopathol (2022) 44(3):281–301. doi: 10.1007/s00281-022-00911-1
3. Liozon E, Dalmay F, Lalloue F, Gondran G, Bezanahary H, Fauchais A-L, et al. Risk factors for permanent visual loss in biopsy-proven giant cell arteritis: A study of 339 patients. J Rheumatol (2016) 43(7):1393–9. doi: 10.3899/jrheum.151135
4. Parreau S, Dumonteil S, Montoro FM, Gondran G, Bezanahary H, Palat S, et al. Giant cell arteritis-related stroke in a large inception cohort: A comparative study. Semin Arthritis Rheum (2022) 55:152020. doi: 10.1016/j.semarthrit.2022.152020
5. Maz M, Chung SA, Abril A, Langford CA, Gorelik M, Guyatt G, et al. 2021 American college of rheumatology/Vasculitis foundation guideline for the management of giant cell arteritis and takayasu arteritis. Arthritis Rheumatol (2021) 73(8):1349–65. doi: 10.1002/art.41774
6. Hellmich B, Agueda A, Monti S, Buttgereit F, de Boysson H, Brouwer E, et al. 2018 Update of the EULAR recommendations for the management of large vessel vasculitis. Ann Rheum Dis (2020) 79(1):19–30. doi: 10.1136/annrheumdis-2019-215672
7. Stone JH, Tuckwell K, DiMonaco S, Klearman M, Aringer M, Blockmans D, et al. Trial of tocilizumab in giant-cell arteritis. N Engl J Med (2017) 377(4):317–28. doi: 10.1056/NEJMoa1613849
8. Jover JA, Hernández-García C, Morado IC, Vargas E, Bafiares A, Fernandez Gutierrez B. Combined treatment of giant-cell arteritis with methotrexate and prednisone. a randomized, double-blind, placebo-controlled trial. Ann Intern Med (2001) 134(2):106–14. doi: 10.7326/0003-4819-134-2-200101160-00010
9. Maleszewski JJ, Younge BR, Fritzlen JT, Hunder GG, Goronzy JJ, Warrington KJ, et al. Clinical and pathological evolution of giant cell arteritis: a prospective study of follow-up temporal artery biopsies in 40 treated patients. Mod Pathol (2017) 30(6):788–96. doi: 10.1038/modpathol.2017.10
10. Bley TA, Wieben O, Uhl M, Thiel J, Schmidt D, Langer M. High-resolution MRI in giant cell arteritis: imaging of the wall of the superficial temporal artery. AJR Am J Roentgenol (2005) 184(1):283–7. doi: 10.2214/ajr.184.1.01840283
11. Weyand CM, Goronzy JJ. Immunology of giant cell arteritis. Circ Res (2023) 132(2):238–50. doi: 10.1161/CIRCRESAHA.122.322128
12. Bolha L, Hočevar A, Suljič A, Jurčič V. Inflammatory cell composition and immune-related microRNA signature of temporal artery biopsies from patients with giant cell arteritis. Front Immunol (2021) 12:791099. doi: 10.3389/fimmu.2021.791099
13. Reitsema RD, van der Geest KSM, Sandovici M, Jiemy WF, Graver JC, Abdulahad WH, et al. Phenotypic, transcriptomic and functional profiling reveal reduced activation thresholds of CD8+ T cells in giant cell arteritis. Rheumatol (Oxford) (2022) 62(1):417–27. doi: 10.1093/rheumatology/keac250
14. Bubak AN, Mescher T, Mariani M, Frietze SE, Hassell JE Jr, Niemeyer CS, et al. Targeted RNA sequencing of forMalin-fixed, paraffin-embedded temporal arteries from giant cell arteritis cases reveals viral signatures. Neurol Neuroimmunol Neuroinflamm (2021) 8(6):e0178. doi: 10.1212/NXI.0000000000001078
15. Kuret T, Lakota K, Čučnik S, Jurčič V, Distler O, Rotar Ž, et al. Dysregulated expression of arterial microRNAs and their target gene networks in temporal arteries of treatment-naïve patients with giant cell arteritis. Int J Mol Sci (2021) 22(12):6520. doi: 10.3390/ijms22126520
16. Merritt CR, Ong GT, Church SE, Barker K, Danaher P, Geiss G, et al. Multiplex digital spatial profiling of proteins and RNA in fixed tissue. Nat Biotechnol (2020) 38(5):586–99. doi: 10.1038/s41587-020-0472-9
17. Ponte C, Grayson PC, Robson JC, Suppiah R, Gribbons KB, Judge A, et al. 2022 American College of Rheumatology/EULAR classification criteria for giant cell arteritis. Ann Rheum Dis (2022) 81(12):1647–53. doi: 10.1136/ard-2022-223480
18. Lie JT. Illustrated histopathologic classification criteria for selected vasculitis syndromes. American College of Rheumatology Subcommittee on Classification of Vasculitis. Arthritis Rheum (1990) 33(8):1074–87. doi: 10.1002/art.1780330804
19. Hernández-Rodríguez J, Murgia G, Villar I, Campo E, Mackie SL, Chakrabarty A, et al. Description and validation of histological patterns and proposal of a dynamic model of inflammatory infiltration in giant-cell arteritis. Med (Baltimore) (2016) 95(8):e2368. doi: 10.1097/MD.0000000000002368
20. Jakobsson K, Jacobsson L, Mohammad AJ, Nilsson JA, Warrington K, Matteson EL, et al. The effect of clinical features and glucocorticoids on biopsy findings in giant cell arteritis. BMC Musculoskelet Disord (2016) 17(1):363. doi: 10.1186/s12891-016-1225-2
21. Narváez J, Bernad B, Roig-Vilaseca D, Garcia-Gomez C, Gomez-Vaquero C, Juanola X, et al. Influence of previous corticosteroid therapy on temporal artery biopsy yield in giant cell arteritis. Semin Arthritis Rheum (2007) 37(1):13–9. doi: 10.1016/j.semarthrit.2006.12.005
22. Parreau S, Liozon E, Chen JJ, Curumthaullee MF, Fauchais A-L, Warrington KJ, et al. Temporal artery biopsy: A technical guide and review of its importance and indications. Surv Ophthalmol (2023) 68(1):104–12. doi: 10.1016/j.survophthal.2022.08.008
23. Parreau S, Vedrenne N, Regent A, Richard L, Sindou P, Mouthon L, et al. An immunohistochemical analysis of fibroblasts in giant cell arteritis. Ann Diagn Pathol (2021) 52:151728. doi: 10.1016/j.anndiagpath.2021.151728
24. Yang X, Kui L, Tang M, Li D, Wei K, Chen W, et al. High-throughput transcriptome profiling in drug and biomarker discovery. Front Genet (2020) 11:19. doi: 10.3389/fgene.2020.00019
25. Feng Y, Wang Q, Wang T. Drug target protein-protein interaction networks: A systematic perspective. BioMed Res Int (2017) 2017:1289259. doi: 10.1155/2017/1289259
26. Barabási AL, Gulbahce N, Loscalzo J. Network medicine: a network-based approach to human disease. Nat Rev Genet (2011) 12(1):56–68. doi: 10.1038/nrg2918
27. Pryshchep O, Ma-Krupa W, Younge BR, Goronzy JJ, Weyand CM. Vessel-specific Toll-like receptor profiles in human medium and large arteries. Circulation (2008) 118(12):1276–84. doi: 10.1161/CIRCULATIONAHA.108.789172
28. Adachi Y, Ueda K, Nomura S, Ito K, Katoh M, Katagiri M, et al. Beiging of perivascular adipose tissue regulates its inflammation and vascular remodeling. Nat Commun (2022) 13(1):5117. doi: 10.1038/s41467-022-32658-6
29. Horimatsu T, Kim HW, Weintraub NL. The role of perivascular adipose tissue in non-atherosclerotic vascular disease. Front Physiol (2017) 8:969. doi: 10.3389/fphys.2017.00969
30. McDonald HM, Farmer JP, Blanco PL. Periadventitial tissue examination in temporal artery biopsies for suspected giant cell arteritis: a case series and literature review. Can J Ophthalmol (2019) 54(5):615–20. doi: 10.1016/j.jcjo.2018.12.011
31. Farr L, Ghosh S, Moonah S. Role of MIF cytokine/CD74 receptor pathway in protecting against injury and promoting repair. Front Immunol (2020) 11:1273. doi: 10.3389/fimmu.2020.01273
32. Zhou Y, Chen H, Liu L, Yu X, Sukhova GK, Yang M, et al. CD74 deficiency mitigates systemic lupus erythematous-like autoimmunity and pathological findings in mice. J Immunol (2017) 198(7):2568–77. doi: 10.4049/jimmunol.1600028
33. Kanemitsu H, Matsunawa M, Wakabayashi K, Sato M, Takahashi R, Odai T, et al. Increased serum levels of macrophage migration inhibitory factor (MIF) in patients with microscopic polyangiitis. Open Access Rheumatol (2009) 13:1–8. doi: 10.2147/oarrr.s4906
34. Abualrous ET, Sticht J, Freund C. Major histocompatibility complex (MHC) class I and class II proteins: impact of polymorphism on antigen presentation. Curr Opin Immunol (2021) 70:95–104. doi: 10.1016/j.coi.2021.04.009
35. Corna G, Campana L, Pignatti E, Castiglioni A, Tagliafico E, Bosurgi L, et al. Polarization dictates iron handling by inflammatory and alternatively activated macrophages. Haematologica (2010) 95(11):1814–22. doi: 10.3324/haematol.2010.02387
36. Weyand CM, Wagner AD, Björnsson J, Goronzy JJ. Correlation of the topographical arrangement and the functional pattern of tissue-infiltrating macrophages in giant cell arteritis. J Clin Invest (1996) 98(7):1642–9. doi: 10.1172/JCI118959
37. Esen I, Jiemy WF, van Sleen Y, van der Geest KSM, Sandovici M, Heeringa P, et al. Functionally heterogenous macrophage subsets in the pathogenesis of giant cell arteritis: Novel targets for disease monitoring and treatment. J Clin Med (2021) 10(21):4958. doi: 10.3390/jcm10214958
38. Oelschlaegel D, Weiss Sadan T, Salpeter S, Krug S, Blum G, Schmitz W, et al. Cathepsin inhibition modulates metabolism and polarization of tumor-associated macrophages. Cancers (Basel) (2020) 12(9):279. doi: 10.3390/cancers12092579
39. Zeng H, Yang X, Yang L, Li W, Zheng Y. Thymosin β10 promotes tumor-associated macrophages M2 conversion and proliferation via the PI3K/Akt pathway in lung adenocarcinoma. Respir Res (2020) 21(1):328. doi: 10.1186/s12931-020-01589-7
40. Deng J, Younge BR, Olshen RA, Goronzy JJ, Weyand CM. Th17 and Th1 T-cell responses in giant cell arteritis. Circulation (2010) 23(7):906–15. doi: 10.1161/CIRCULATIONAHA.109.872903
41. Rodríguez-Pla A, Bosch-Gil JA, Rosselló-Urgell J, Huguet-Redecilla P, Stone JH, Vilardell-Tarres M. Metalloproteinase-2 and -9 in giant cell arteritis: involvement in vascular remodeling. Circulation (2005) 112(2):264–9. doi: 10.1161/CIRCULATIONAHA
42. Rittner HL, Kaiser M, Brack A, Szweda LI, Goronzy JJ, Weyand CM. Tissue-destructive macrophages in giant cell arteritis. Circ Res (1999) 84(9):1050–8. doi: 10.1161/01.res.84.9.1050
43. Kaiser M, Weyand CM, Björnsson J, Goronzy JJ. Platelet-derived growth factor, intimal hyperplasia, and ischemic complications in giant cell arteritis. Arthritis Rheum (1998) 41(4):623–33. doi: 10.1002/1529-0131(199804)41:4<623::AID-ART9>3.0.CO;2-6
44. Banerjee S, Bagheri M, Sandfort V, Ahlman MA, Malayeri AA, Bluemke DA, et al. Vascular calcification in patients with large-vessel vasculitis compared to patients with hyperlipidemia. Semin Arthritis Rheum (2019) 48(6):1068–73. doi: 10.1016/j.semarthrit.2018.09.001
45. Segarra M, García-Martínez A, Sánchez M, Hernandez-Rodriguez J, Lozano E, Grau JM, et al. Gelatinase expression and proteolytic activity in giant-cell arteritis. Ann Rheum Dis (2007) 66(11):1429–35. doi: 10.1136/ard.2006.068148
46. Samson M, Genet C, Corbera-Bellalta M, Greigert H, Espigol-Frigole G, Gerard C, et al. Human monocyte-derived suppressive cells (HuMoSC) for cell therapy in giant cell arteritis. Front Immunol (2023) 14:1137794. doi: 10.3389/fimmu.2023.1137794
47. Corbera-Bellalta M, García-Martínez A, Lozano E, Planas-Rigol E, Tavera-Bahillo I, Alba MA, et al. Changes in biomarkers after therapeutic intervention in temporal arteries cultured in Matrigel: a new model for preclinical studies in giant-cell arteritis. Ann Rheum Dis (2014) 73(3):616–23. doi: 10.1136/annrheumdis-2012-202883
48. Corbera-Bellalta M, Alba-Rovira R, Muralidharan S, Espigol-Frigole G, Rios Garces R, Marco-Hernandez J, et al. Blocking GM-CSF receptor α with mavrilimumab reduces infiltrating cells, pro-inflammatory markers and neoangiogenesis in ex vivo cultured arteries from patients with giant cell arteritis. Ann Rheum Dis (2022) 81(4):524–36. doi: 10.1136/annrheumdis-2021-220873
49. Cid MC, Unizony SH, Blockmans D, Brouwer E, Dagna L, Dasgupta B, et al. Efficacy and safety of mavrilimumab in giant cell arteritis: a phase 2, randomised, double-blind, placebo-controlled trial. Ann Rheum Dis (2022) 81(5):653–61. doi: 10.1136/annrheumdis-2021-221865
50. Christian BA, Poi M, Jones JA, Porcu P, Maddocks K, Flynn JM, et al. The combination of milatuzumab, a humanized anti-CD74 antibody, and veltuzumab, a humanized anti-CD20 antibody, demonstrates activity in patients with relapsed and refractory B-cell non-Hodgkin lymphoma. Br J Haematol (2015) 169(5):701–10. doi: 10.1111/bjh.13354
Keywords: giant cell arteritis, vasculitis, transcriptomic analysis, temporal artery biopsy, gene expression
Citation: Parreau S, Molina E, Dumonteil S, Goulabchand R, Naves T, Bois MC, Akil H, Terro F, Fauchais A-L, Liozon E, Jauberteau M-O, Weyand CM and Ly K-H (2023) Use of high-plex data provides novel insights into the temporal artery processes of giant cell arteritis. Front. Immunol. 14:1237986. doi: 10.3389/fimmu.2023.1237986
Received: 10 June 2023; Accepted: 21 August 2023;
Published: 06 September 2023.
Edited by:
Gilles Kaplanski, Assistance Publique Hôpitaux de Marseille, FranceReviewed by:
Maxime Samson, Université de Bourgogne, FranceAlessandro Mangogna, Institute for Maternal and Child Health Burlo Garofolo (IRCCS), Italy
Copyright © 2023 Parreau, Molina, Dumonteil, Goulabchand, Naves, Bois, Akil, Terro, Fauchais, Liozon, Jauberteau, Weyand and Ly. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Simon Parreau, c2ltb24ucGFycmVhdUBjaHUtbGltb2dlcy5mcg==
†These authors share junior authorship
‡These authors share senior authorship