Microbiomes of Inflammatory Thoracic Aortic Aneurysms Due to Giant Cell Arteritis and Clinically Isolated Aortitis Differ From Those of Non-Inflammatory Aneurysms

Authors

Ted M. Getz¹*, Gary S. Hoffman^2,*, Roshan Padmanabhan¹, Alexandra Villa-Forte², Eric E. Roselli^3,, Eugene Blackstone^3,, Douglas Johnston ^3,, Gosta Pettersson^3,, Edward Soltesz^3,, Lars G. Svensson^3,, Leonard H. Calabrese², Alison H. Clifford^2,5, Charis Eng^1,4,6,7

Affiliated Institutions

¹ Genomic Medicine Institute; Lerner Research Institute; Cleveland Clinic; Cleveland, Ohio
² Center for Vasculitis Care and Research; Department of Rheumatic and Immunologic Diseases; Cleveland Clinic; Cleveland, Ohio
³ Center for Aortic Diseases; Heart Vascular Institute; Cleveland Clinic; Cleveland, Ohio
⁴ Taussig Cancer Institute; Cleveland Clinic; Cleveland, Ohio
⁵ Division of Rheumatology; Department of Medicine; University of Alberta, Edmonton, Alberta T6G 2R7, Canada
⁶ Department of Genetics and Genome Sciences; Case Western Reserve University School of Medicine; Cleveland, Ohio
⁷ Germline High Risk Focus Group; CASE Comprehensive Cancer Center; Case Western Reserve University School of Medicine; Cleveland, Ohio

*Joint first authors

Corresponding Authors

Gary S. Hoffman
Cleveland Clinic Department of Rheumatic and Immunologic Diseases
9500 Euclid Avenue
A50, Cleveland, OH 44195
(216) 445 6996
[email protected]

Charis Eng
Cleveland Clinic Genomic Medicine Institute
9500 Euclid Avenue
NE50
Cleveland, OH 44195
(216) 444 3440
[email protected]

DOI

Suggested Citation

Abstract

Objective: We sought to characterize microbiomes of thoracic aortas from patients with non-infectious aortitis due to giant cell arteritis (GCA) and clinically isolated aortitis (CIA) and to compare them to non-inflammatory aorta aneurysm controls. We also compared microbiomes from concurrently processed and separately reported temporal arteries (TA) and aortas.

Methods: From 220 prospectively enrolled patients undergoing surgery for thoracic aorta aneurysm, 49 were selected. Inflammatory and non-inflammatory cases were selected based on ability to match for age (+/-10 years), gender, and race. Biopsies were collected under aseptic conditions and snap-frozen. Taxonomic classification of bacterial sequences was performed to the genus level and relative abundances were calculated. Microbiome differential abundances were analyzed by principal coordinates analysis.

Results : Forty-nine patients with thoracic aortic aneurysms (12 CIA, 14 GCA, 23 non-inflammatory aneurysms) were enrolled. Alpha (P = 0.018) and beta (P = 0.024) diversity differed between specimens from aortitis cases and controls. There were no significant differences between CIA and GCA (P > 0.7). The largest differential abundances between non-infectious aortitis and non-inflammatory control samples includedEnterobacteriaceae, Phascolarctobacterium, Acinetobactor, Klebsiella, and Prevotella. Functional metagenomic predictions with PICRUSt revealed enrichment of oxidative phosphorylation and porphyrin metabolism pathways and downregulation of transcription factor pathways in aortitis compared to controls. Microbiomes of aortic samples differed significantly from temporal artery samples from a companion study, in both control and GCA groups (P = 0.0002).

Conclusion: Thoracic aorta aneurysms, far from being sterile, contain unique microbiomes that differ from those found in temporal arteries. The aorta microbiomes are most similar between aneurysms that were associated with inflammation, GCA, and CIA, but differed from those associated with non-inflammatory etiologies. These findings are promising in that they indicate that microbes may play a role in the pathogenesis of aortitis-associated aneurysms or non-inflammatory aneurysms by promoting or protecting against inflammation. However, we cannot rule out that these changes are related to alterations in tissue substrate that favor secondary changes in microbial communities.

Keywords: aorta; aneurysms; aortitis; microbiome; clinically isolated aortitis; giant cell arteritis

Introduction

Aortitis is a clinical-histopathological diagnosis. Pre-operatively, it is often unrecognized; the diagnosis is usually made only after aortic aneurysm surgery and microscopic review [1]. Clinical and imaging evaluation is then required to determine diagnostic subsets of aortitis. Most forms of non-infectious aortitis have been presumed to be autoimmune and the result of a primary systemic large vessel vasculitis (LVV), such as giant cell arteritis (GCA) or Takayasu’s arteritis (TAK), or other systemic autoimmune diseases (eg, rheumatoid arthritis, the spondyloarthritides, Behcet’s disease, Cogan syndrome, relapsing polychondritis, systemic lupus erythematosus, Sjogren syndrome, granulomatosis with polyangiitis, sarcoidosis, or IgG4-related disease) [1, 2]. Aortitis may also appear in isolation, as a topographically limited lesion. This latter group has been the focus of recent consensus guidelines and coining the term “clinically isolated aortitis” (CIA) to indicate a focal-regional non-infectious form of vasculitis usually restricted to the proximal thoracic aorta [3, 4]. In any given case of CIA, it is uncertain whether vasculitis is pathologically truly isolated and potentially cured by surgical resection, or if it represents a more extensive vasculitis that is subclinical, undiscovered by imaging beyond the surgical site, and destined to be expressed as a systemic LVV in the future [3]. While some previous studies have found that patients with CIA are not at increased risk of future vascular events, most have found that many patients will develop new distal aortic aneurysms and/or other large vessel lesions [5, 6]. In a recent study of 196 patients with non-infectious aortitis, 66% had diagnoses of CIA, 21% had GCA, 7% had TAK, and 6% had other systemic inflammatory diseases, at the time of surgery. During a mean follow-up of over 4.5 years, 19% of patients with CIA developed new symptoms, 45% developed new radiographic vascular lesions, 40% underwent additional vascular surgery, and 12% died. Fifteen percent of patients initially classified with CIA developed features of a systemic disease, most often GCA [7]. Among numerous questions raised by such studies are the following: First, is CIA a limited presentation of diseases within a spectrum of LVV, including GCA? Second, do common conditions exist within CIA-affected and GCA-affected aortas and GCA-affected temporal arteries that favor developing inflammatory aneurysms or temporal artery inflammation and stenoses? And third, if this were the case, would those conditions differ from that present in non-inflammatory thoracic aortic aneurysms? Because of our interests in infectious triggers and vascular properties that favor the genesis of vasculitis, we examined the vascular microbiome in subsets of aortitis and non-inflammatory aorta aneurysm controls [8-11]. Unlike prior microbiome studies of large vessels, we collected surgically aseptic specimens that were snap-frozen and not fixed in formalin or paraffin-embedded.

Materials and Methods

Sample Accrual and Collection

We prospectively enrolled 220 consecutive consenting patients undergoing thoracic aorta aneurysm surgery. Aortic samples were collected under strictly aseptic conditions. All aorta and temporal artery specimens (from a linked study [12]) were processed at the same time and analyzed in a blinded fashion.

The study protocol was conducted in compliance with the Helsinki Declaration and approved by the Institutional Review Board for Human Subjects’ Protection at the Cleveland Clinic. Aorta samples were split into thirds, with one-third sent for routine histopathological review, one-third snap-frozen under sterile technique for microbiome analysis, and one-third utilized for a related biochemistry study [13]. Patients were classified according to clinical phenotype and histopathology as having either biopsy-positive and clinically compatible GCA (N = 14) or clinically isolated aortitis (CIA, N = 12), or non-inflammatory aneurysms (controls). Control aortas (N = 190) were obtained from patients with aneurysms associated with diseases such as uncontrolled hypertension, bicuspid aortic valves, Marfan syndrome, and cystic medial degeneration of uncertain etiology. Controls (N = 23) were selected to match aortitis cases (N=26) based on same gender, race and age (+/-10 years) (Table 1). Following > 3 years of collection, all aorta and temporal artery samples, which were blinded to diagnosis, were processed at the same time.

Table 1. Demographic Characteristics, Comorbid Conditions, Substance Use History, and Relevant Laboratory Values of Patients with CIA, GCA, and Non-Inflammatory Aortitis Aneurysms.

*Statistical analysis by 2-way ANOVA with P < 0.05 considered significant. All quantitative values are shown as mean +/- standard deviation, and all qualitative values are shown as number of patients with characteristic (percentage of group in column with characteristic). **One each with rheumatoid arthritis and Sjogren syndrome.

16S Ribosomal RNA Gene Sequencing of Aorta Biopsies

Total deoxyribonucleic acid (DNA) was isolated using the RNeasy PowerMag KitAllPrep PowerViral DNA/RNA Isolation Kit according to the manufacturer’s protocol (Qiagen, Valencia, CA) with minor modifications [14]. Briefly, all beads, tubes, and non-enzymatic reagents were treated with UV light for 30 minutes prior to use; samples were digested with 20 μL of 20 ng/μL Proteinase K (Roche Diagnostics Corp., Indianapolis, IN) at 65^oC degrees for 1 hour, then transferred to 0.1 mm glass bead-containing tubes, after which the samples were homogenized using the TissueLyser LT (Qiagen). The quality and purity of the isolated total DNA were confirmed spectrophotometrically using a NanoDrop 2000 device (Fisher Scientific SAS, Illkirch, France). DNA concentration was quantified using the Qubit 2.0 instrument applying the Qubit dsDNA HS Assay (Life Technologies, USA). Extracted DNA samples were stored at -20°C.

Bacterial 16S rRNA gene amplification and library construction were performed according to the 16S Metagenomic Sequencing Library Preparation guide from Illumina (Forest City, CA). In brief, 2 μL total DNA was amplified using primers targeting the 16S V3 and V4 region (Illumina) at 95^oC for 5 minutes, followed by 35 cycles of 95^oC for 30 seconds, 72^oC for 30 seconds, and a final extension of 72^oC for 10 minutes. Then 16S rDNA amplicons were run out on a 1% agarose gel, size-selected at 450-500 bp, and gel-purified using QIAquick Gel Purification kit (Qiagen). A second round of PCR was performed to add Nextera XT indices (Illumina) to purified amplicons. Indexed PCR products were cleaned with Ampure XP beads (Beckman Coulter, Inc., Brea, CA) and resulting libraries quantified with the QuantiFluor dsDNA system according to the manufacturer’s protocol (Promega, Madison, WI). Samples were then normalized to 10nM and pooled into sequencing libraries. Pooled V3-V4 amplicon libraries were sequenced using the Illumina MiSeq platform with V3 reagent kit. The 300-bp paired-end reads for each sample were demultiplexed and quality checked using FastQC 0.11.3.

Because the specimens were stored in aluminum foil and sterile containers, 16S rRNA sequencing was performed from both and revealed identical results. The containers had sterilized water aliquoted into them which was then utilized as a template for 16S rRNA gene sequencing, revealing no contamination.

Microbiome Analysis

All microbiome analysis was performed in a blinded fashion, in the first instance, after which a co-author not involved in the microbiome analysis unblinded meta-data.

A hybrid post-sequencing analysis methodology using QIIME and MICCA was adopted, in which preprocessing was performed in QIIME and open-reference operational taxonomic unit (OTU) picking was performed with MICCA and Phyloseq [15]. After the biom files were created, downstream analysis was performed with QIIME. Then 250-bp Illumina Paired-end reads were merged with FLASH [16] and low quality reads (Phred < 20) were filtered out using the split_libraries.py command in QIIME (version 1.9) [17]. MICCA vsearch (version 1.9.5)[18] was utilized for clustering the sequences with a threshold of 97% similarity, and representative sequences were classified using RDP classifier (version 2.11) [19]. Multiple sequence alignments were done on the comprehensive sequence using MUSCLE (version 3.8.31) [20, 21] against the Greengenes database (version 13.8) [22], filtered at 97% similarity, and FastTree (version 2.1.8) was used for phylogenetic tree construction [23].

Data clean-up was performed by removing the singletons and discarding taxa represented in fewer than 5% of total samples. The rarefaction value was set to 1107 reads per sample to reduce sampling heterogeneity, and computation of alpha (Shannon diversity index) and beta diversity measures (unweighted UniFrac distances) were performed with phyloseq in R. Alpha diversity measures species richness (number of taxa) within a single microbial ecosystem. Beta diversity can be represented by UniFrac distances which describe similarities and dissimilarities between bacterial communities using phylogenetic information, taking into account number of taxa and relative abundances within each taxon (Table 2). F-tests based on sequential sums of squares derived from 1,000 permutations on UniFrac distance matrices were performed with the null hypothesis that there is no difference in community structure between groups. Note that PCoA and the calculation of P values are measurements of clustering strength. Differences (and the P value) are derived from measuring differences between UniFrac distances. To find which taxa are most likely to explain the differences between our clinical groupings, taxa summaries and differential abundances were analyzed with DESeq2 with data from MICCA input into the R package DESeq2. This algorithm estimates variance-mean dependence in count data and tests for differential expression based on a model using the negative binomial distribution. Differentially abundant taxa that were statistically significant using an alpha of 0.01 and exceeded a log2-fold change of ±2 were visually represented on box plots. Comparisons between aorta and temporal artery microbiomes were performed using microbiome dataset collected from a parallel study of 47 temporal arteries with and without GCA [12]. However, both the aorta and temporal artery samples were sequenced at the same time to reduce batch effect.

Table 2. Measures of Microbial Diversity as Applied in this Study

* Genus or species

**OTU = operational taxonomic units

We predicted functional composition of microbiomes using the PICRUSt 1.0.0-dev bioinformatics package [24]. We filtered out all de novo OTUs and used this OTU table as our input into the PICRUSt algorithm, which calculates contributions of various OTUs to known biological pathways based on evolutionary modeling. Welch’s t test was used to calculate P values, and corresponding Storey q-values were used to control for the false discovery rates associated with multiple testing. These values were calculated using DESeq2 and LEfSe and visualized as a bar plot using STAMP [25].

Results

Patients

Twenty-six patients with aortitis-associated thoracic aneurysms, including 14 with GCA and 12 with CIA, and 23 controls with non-inflammatory aorta aneurysm, were enrolled in our study (Table 1). There were no differences in demographics, exposures, and relevant past history among the 3 groups, except for tobacco use (P = 0.02), a history of coronary artery disease (P = 0.05), and steroid use in 4/14 patients with GCA (Table 1). Importantly, there was no significant difference in concurrent systemic autoimmune disease across the groups (P = 0.31, Table 1). Nine of 14 patients with GCA had an ESR > 50 mm/hr compared to none in both the CIA group and controls (P = 0.01, Table 1). Other basic laboratory values did not differ across the 3 groups (P > 0.7, Table 1). Additionally, none of the patients were taking antibiotics at the time of their surgical procedures.

Histopathology

Of the 26 eligible patients with aortitis, 14 received a final diagnosis of biopsy-proven GCA revealing typical arteritis with mononuclear cell inflammatory infiltrates localized to the adventitia and media, giant cell granulomatous changes (9/14), fragmentation of the internal elastic lamina, and varying degrees of intimal proliferation and fibrinoid necrosis. The remaining 12 were diagnosed with CIA with giant cell granulomatous changes (8/12), a lymphoplasmacytic pattern of inflammation, a suppurative pattern of inflammation, or a mixed inflammatory pattern. Among the 23 control aortas, 21 had diagnoses that are considered to contribute to thoracic aorta aneurysm formation including bicuspid aortic valve (8), atherosclerosis with calcification of the aortic valve (5), cystic medial degeneration of uncertain etiology (3), hypertension (14), and Marfan syndrome (1).

Microbiome in Aortas From Patients With and Without Aortitis

Culture-independent, long-read genomic sequencing was used to characterize the entire microbial communities of aortas from the 49 research participants with aorta aneurysm. After sequencing and quality control, 10 samples, all from controls, were excluded from further analyses due to low read count, leaving all 26 samples with aortitis and 13 controls. Alpha diversity (Table 2) was reduced in aortitis samples compared to control aortas (P = 0.018, Figure 1A). Beta diversity (Table 2), as measured by UniFrac distances, differed in aortitis samples compared to controls (P = 0.024, Figure 1B). Importantly, there were no differences in alpha or beta diversity between GCA aorta samples and CIA aorta samples (P > 0.7, data not shown).

Figure 1. Differences in microbiome in aortas from patients with aortitis compared to those without aortitis. (A) Rarefaction (Shannon diversity) curves showing decreased alpha diversity in aortitis samples (red curve) compared to control aortas (blue) [P = 0.018]. (B) Principal component analysis (PCoA) of aorta microbiomes from patients with aortitis (red) compared to those from control patients (green) [P = 0.024].

Irrespective of diagnosis, there was borderline decreased beta diversity in aortas from patients with a history of hypertension compared to normotensive patients (P = 0.05), which disappeared when microbiomes were compared within the aortitis group (P = 0.2) and within the control group (P = 0.7). Of note, the histories of hypertension were not different between aortitis cases and controls (Table 1).

We then examined taxonomic classifications and relative abundances of bacterial species within aorta samples in order to determine whether differential taxon abundances exist in aortitis specimens compared to non-inflamed control aortas. At the phylum level, OTUs from Actinobacteria and Firmicutes were relatively less abundant in aortitis specimens compared to control aortas (Figure 2). An OTU from the genus Phascolarctobacterium had the greatest over-representation in aortitis samples compared to control aortas, followed by OTUs from the family Enterobacteriaceae and the genus Rothia (Figure 2). The most differentially less abundant OTU in aortitis compared to controls was the genus Prevotella, followed by the generaAcinetobacter, Klebsiella, Staphylococcus, and Cornybacterium. There were also 7 less abundant OTUs that were unclassified at the genus level from phylum Actinobacteria, 15 belonging to the genus Propionbacterium, 1 to the genus Onchrobactrum, and 9 to the class Clostridales (Figure 2).

Figure 2. Most differentially abundant taxa in aortas from patients with aortitis and from control patients. Bar plot representation showing the most (and statistically significant, all P < 0.05) over-represented (+) (right) and under-represented (-) (left) taxa in aortas from aortitis patients compared to those from control patients.

Predicting Functional Consequences of Differing Microbial Composition in Aortitis Specimens and Control Aortas

Metagenome functional content was predicted using PICRUSt, and group comparisons between the aortitis and control groups were performed with STAMP (Figure 3). With the different microbes in aortitis samples compared to control aortas as input, PICRUSt predicted relative upregulation of pathways involved in oxidative phosphorylation, porphyrin metabolism, TCA cycle, streptomycin biosynthesis, and glycine, serine, threonine, and histidine metabolism in aortitis. In contrast, predicted relative downregulated pathways included transcription factors, chromosome maintenance, the sulfur relay system, and selenocompound metabolism.

Figure 3. Predicted functional pathways differentially represented in aorta samples from patients with aortitis compared to those from controls. Representation of PICRUSt analysis yielding differentially regulated functional pathways in aortas from cases (blue bars) and controls (orange bars).

Comparison of Microbiomes in Aortas and Temporal Arteries

Gene sequencing data of 16S rRNA were compared between our 49 aorta samples and 47 temporal arteries, the latter from a parallel study [12], 23 without GCA and 24 with GCA. There was a significant difference in beta diversity between the microbiome of aorta and temporal arteries (P = 0.0002, Figure 4A). Microbiomes from GCA-affected temporal arteries [12] versus GCA-affected aortas were significantly different (P = 0.001), as was also noted for the non-inflammatory aortas and non-inflammatory temporal artery controls (P = 0.001). Of note, the separation by UniFrac distances by PCoA was qualitatively more marked in control temporal arteries versus control aortas (Figure 4B, almost complete separation) compared to those of temporal artery-GCA and aortitis-GCA (where there is overlap, Figure 4C).

Figure 4. Microbiome from aortas compared to temporal arteries. (A) PCoA of microbiomes from all aorta specimens (red) and all temporal arteries (green). (B) PCoA of microbiomes from control aortas (red) and control temporal arteries (green). (C) PCoA of microbiomes from aortitis specimens (red) and temporal arteritis samples (green).

Discussion

Our results suggest that, far from being sterile, aortas contain a microbiome in both inflammatory and non-inflammatory states. Interestingly, we found that the aortic microbiome of patients with inflammatory conditions differed substantially from our controls in both alpha diversity and beta diversity. However, we found no difference in diversity between the LVV subsets, GCA and CIA, within our aortitis group. Additionally, the observations that the histopathology of GCA and CIA are indistinguishable and that a significant minority of patients with clinically isolated aortitis phenotype may evolve into a giant cell arteritis phenotype [10], raise questions about shared critical pathogenetic mechanisms for these entities. Conversely, we do not know whether the inflammatory state and/or its consequences observed in these forms of LVV could impact the composition of the observed microbiome. Alternatively, these microbial differences could be influenced by variables not examined in these studies.

Prior studies of GCA pathogenesis have focused on temporal arteries.Mycoplasma pneumonia, Chlamydia pneumonia, andBurkholderia pseudomallei-like organisms [26-28] have been implicated in pathogenesis; however, none have been confirmed by others nor were they detected in our aorta dataset or in a separate study we performed on temporal arteries [12].

Not surprisingly, we found that the microbiomes of temporal arteries and aortas differed both between the GCA-affected temporal arteries and aortas, and the non-inflammatory aortas and non-inflammatory temporal artery controls. There are many precedents for the same stimuli having a range of effects in different tissues and even territories within a shared network (eg, the vasculature). Consider unique site predilections in atherosclerosis (eg, the infra-renal aorta) and even infection-mediated vasculitis (eg, hepatitis C, preferred targeting of skin, kidneys, and peripheral nerves). Disease targeting or sparing even within the same vessel (eg, proximal and distal aorta) may differ for many reasons including differences in embryogenic origins of smooth muscle cells within the aortic root and arch versus the more distal aorta, gene expression, and focal immune response [8].

The most thorough survey of the temporal arterial microbiomes so far examined the role of microbes in GCA [29]. Bhatt and colleagues utilized a whole genome sequencing approach to search for a pathogen in 17 formalin-fixed, paraffin-embedded temporal artery samples. More than 1000 non-human microbial species were detected; however, by far the most common were Propionibacterium acnes and Escherichia coli, which are commonly found within human skin flora. Moreover, these authors did not detect any significant differences between GCA-affected temporal arteries versus controls in regard to the relative abundance of bacterial species nor did they find any previously reported viral pathogens. Bhatt et al provided insight into the spectacular diversity of the vascular microbiome, but their study suffers from its use of non-sterile collection methods and the use of formalin-fixed, paraffin embedded tissue [29]. Consistent with these observations, we found that Proprionibacterium are common in the vasculature. Moreover, we found significant differences in the microbiomes of aortitis and control vessels. Our study had the advantage of using both aortic tissue and temporal arteries [12], as well as using fresh frozen samples collected aseptically, which we believe allows for a more accurate representation of constituents of the vascular microbiome, while minimizing microbial contaminants. Importantly, all samples and meta-data were blinded during preparation, sequencing and microbiome/bioinformatics analyses. That the vascular microbiome differs between the 2 conditions studied, inflammatory and non-inflammatory aorta control specimens, processed at the same time, also argues strongly against contamination.

In this study, we not only found distinct microbial communities in our aortitis and control groups, but also found that the aortitis specimens had significantly less phylogenic diversity than our control specimens. The finding that decreased microbial diversity may promote inflammation is well established. For example, in inflammatory bowel disease (IBD), a disease for which there is an extensive microbiome literature, among the most consistent findings is a reduction in overall microbial diversity, including species that are thought to have anti-inflammatory properties [30]. Decreased bacterial diversity has been associated with other rheumatologic diseases, even in areas not traditionally thought to contain a microbiome. A recent study of 30 patients with early rheumatoid arthritis (RA) or sarcoidosis, found that the lungs of patients with rheumatologic disease had distal airway dysbiosis with 40% fewer OTUs when compared to the microbiome of healthy controls [31]. Thus, 1 hypothesis that may be derived from our results is that rather than aortitis being incited by a single or a few pathogenic strains, bacterial diversity may actually confer protection against inflammation and loss of important commensals may promote disease susceptibility.

Our findings included a relative over-representation in the generaPhascolarctobacterium and Rothia and the family Enterobacteriaceae in the aortitis group. All of these are ubiquitous components of the human skin and gastrointestinal flora and known human pathogens with the potential to drive inflammation. Enterobacteriaceae has been implicated as a trigger of inflammation in both IBD and colon cancer and has been demonstrated to activate the inflammosome, promoting pro-inflammatory cytokine secretion [32]. Some of the over-represented genera have the potential to become pathogenic, but they do not dominate as would be expected with a frank bacterial infection; rather, they form a general pro-inflammatory milieu, which may contribute to the development of aortitis.

Interestingly, we found numerous OTUs were under-represented in our aortitis group when compared to the control aortas, consistent with the observed decrease in overall microbial diversity in the aortitis group. These included the genera Prevotella,Acinetobacter, Klebsiella, Staphylococcus, and Cornybacterium, many of which have been found to be protective against inflammation on mucosal surfaces. Additionally, Klebsiella and Staphylococcus have been previously isolated from atherosclerotic plaques and are known to form biofilms, allowing for adherence to vessel walls [11]. Prevotella, a GI commensal, is commonly observed to be elevated in several inflammatory conditions; however, it has also been found to be less frequent in inflammatory airway diseases such as asthma and chronic obstructive pulmonary disease [33, 34]. Similarly, 1 study found that increasing abundance of Prevotella was associated with lower levels of inflammatory cytokines IL-10, IL-23, and TNF-α in healthy lungs [34]. Acinetobacte has also been found to be protective against airway inflammation in asthma, through suppression of Th2 mediated T-cell responses to allergens [35]. Together, the reduction in these OTUs may have the effect of facilitating inflammation via loss of anti-inflammatory metabolites and host-microbe interactions with the immune system.

Prediction of up and downregulation of functional pathways based on the composition of the microbiome revealed a relative increase in oxidative phosphorylation and TCA cycle pathways in the aortitis group. One possible explanation for this is that it reflects an increase in aerobic and facultative aerobic microbes, including Enterobacteriaceae and Rothia. Porphyrins are a group of metabolites that have previously been linked to inflammation in acne via production of reactive oxygen species [36]. Several microbial species have been found to increase the secretion of porphyrins and upregulate enzymes responsible for porphyrin metabolism when exposed to hemin, an iron containing porphyrin found in hemoglobin [37, 38]. Consistent with these reports, we found that porphyrin metabolism was predicted to be upregulated in our aortitis group. It is tempting to speculate that these porphyrin metabolites may play a role in inducing inflammation in aortitis. Understanding the predicted downregulation in pathways responsible for chromosomal maintenance and transcription factors is more opaque, but it is possible that these pathways reflect decreased microbial replication in the setting of immune activation. However, further studies are required to clearly establish a firm relationship between the predicted up and downregulation of pathways and the microbes we detected in the inflammatory and non-inflammatory aorta samples.

The notion of a vascular microbiome is neither new nor novel [11]. Many studies have demonstrated the presence of both bacteria and viruses within the walls of blood vessels. The best evidence for the role of bacteria in vascular pathology comes from the isolation of microbes from atherosclerotic plaques, including C. pneumoniae, S. epidermidis, P. acnes, P. gingivalis, P. granulosum,Micrococcus sp., and A. viscosus [39, 40]. Moreover, the discovery of T cells specific for C. pneumoniae and P. gingivalis implicate microbes in the pathogenesis of atherosclerosis as does the finding that mice inoculated withP. gingivalis experienced accelerated atherosclerosis [41-43].

Recent studies have also lent support for the idea that apparently normal vessels may contain a stable microbiome. One of the most compelling of these studies applied a 16S-metagenomics sequencing approach to 56 fresh, sterile aortic aneurysm samples from patients with atherosclerosis and non-atherosclerotic disease. Bacterial DNA was isolated from about 90% of samples, and 10 samples were selected for speciation [39]. All but 1 specimen revealed the presence of multiple bacterial species. Similarly, a study of healthy femoral and coronary arteries in bypass grafts, using 16S sequencing, found 245 bacterial species and 90 genera, including Streptococcus, which was a major component of both the control and aortitis microbiome in our investigation [44].

The major limitations of our study are relatively small sample size and the inability to procure sterile healthy aortas as controls. For controls, we utilized samples from patients undergoing aorta reconstruction for aneurysms caused by non-inflammatory conditions such as hypertension or congenital/genetic anomalies. While we found that our non-inflammatory aorta aneurysm controls differed significantly from the aortitis-aneurysm samples, it is possible that our controls differ from healthy vessels in unforeseen ways. However, our study did have several important strengths. It is the first to be performed in a blinded fashion on surgically sterile aortas that were maintained under strict aseptic conditions throughout processing, reducing the possibility of contamination. Moreover, we used aortitis samples from 2 forms of LVV (GCA and CIA) that are histologically indistinguishable and may be pathogenically related, and compared their microbiomes to non-inflammatory aortic aneurysm controls.

In conclusion, our results emphasize that thoracic aorta aneurysms are not sterile but contain microbial populations regardless of diagnosis. These microbiomes differ between topographically distinct GCA-affected aortas and temporal arteries. While that could argue against specific microbial agents in these populations having an etiologic role, it is possible that these distinct populations could foster changes that have similar pathogenetic effects. Additionally, the microbiomes of patients that have aortitis have significantly less phylogenic diversity than our controls, and they cluster into a distinct group when compared to controls. Aortitis cases also had an increased prevalence of known pathogens that may contribute to a pro-inflammatory environment. While these findings imply possible roles for microbiota in LVV, further studies are needed to elucidate whether the observed differences in microbial populations in idiopathic aortitis versus non-inflammatory aneurysms are pathogenic or secondary to vessel substrate injury.

Acknowledgements

We thank the nurses and coordinators from the Cleveland Clinic Clinical Research Unit for help with study coordination and data collection, and especially, the patients who volunteered their samples to this study.

Funding

This work was supported, in part, by the Fasenmyer Clinical Immunology Center (to GSH, LC, and CE), and the National Center for Advancing Translational Sciences (NCATS) of the NIH (UL1TR000439 to GSH). CE is the Sondra J. and Stephen R. Hardis Endowed Chair in Cancer Genomic Medicine at the Cleveland Clinic, and an American Cancer Society Clinical Research Professor.

Conflict of Interest

Gary S. Hoffman and Leonard H. Calabrese are associate editors for Pathogens and Immunity. They did not participate in any aspect of the review or editorial process of this submitted manuscript.

References

1. Gornik HL, Creager MA. Aortitis. Circulation. 2008;117(23):3039-51. PubMed PMID: 18541754. Pubmed Central PMCID: 2759760. doi: 10.1161/CIRCULATIONAHA.107.760686

2. Belizna CC, Hamidou MA, Levesque H, Guillevin L, Shoenfeld Y. Infection and vasculitis. Rheumatology (Oxford). 2009;48(5):475-82. PubMed PMID: 19258377. doi: 10.1093/rheumatology/kep026

3. Stone JR, Bruneval P, Angelini A, Bartoloni G, Basso C, Batoroeva L, Buja LM, Butany J, d'Amati G, Fallon JT, Gittenberger-de Groot AC, Gouveia RH, Halushka MK, Kelly KL, Kholova I, Leone O, Litovsky SH, Maleszewski JJ, Miller DV, Mitchell RN, Preston SD, Pucci A, Radio SJ, Rodriguez ER, Sheppard MN, Suvarna SK, Tan CD, Thiene G, van der Wal AC, Veinot JP. Consensus statement on surgical pathology of the aorta from the Society for Cardiovascular Pathology and the Association for European Cardiovascular Pathology: I. Inflammatory diseases. Cardiovasc Pathol. 2015;24(5):267-78. PubMed PMID: 26051917. doi: 10.1016/j.carpath.2015.05.001

4. Cinar I, Wang H, Stone JR. Clinically isolated aortitis: pitfalls, progress, and possibilities. Cardiovasc Pathol. 2017;29:23-32. PubMed PMID: 28500877. doi: 10.1016/j.carpath.2017.04.003

5. Miller DV, Isotalo PA, Weyand CM, Edwards WD, Aubry MC, Tazelaar HD. Surgical pathology of noninfectious ascending aortitis: a study of 45 cases with emphasis on an isolated variant. Am J Surg Pathol. 2006;30(9):1150-8. PubMed PMID: 16931960. doi: 10.1097/01.pas.0000213293.04026.ec

6. Wang H, Smith RN, Spooner AE, Isselbacher EM, Cambria RP, MacGillivray TE, Stone JH, Stone JR. Giant cell aortitis of the ascending aorta without signs or symptoms of systemic vasculitis is associated with elevated risk of distal aortic events. Arthritis Rheum. 2012;64(1):317-9. PubMed PMID: 21953530. doi: 10.1002/art.33343

7. Clifford AH, Arafat A, Idrees JJ, Roselli EE, Tan CD, Rene Rodriguez E, Svensson LG, Blackstone E, Johnston D, Pettersson G, Soltesz E, Hoffman GS. Outcomes among 196 patients with non-infectious proximal aortitis. Arthritis Rheumatol. 2019. PubMed PMID: 30730604. doi: 10.1002/art.40855

8. Hoffman GS. Giant Cell Arteritis. Ann Intern Med. 2016;165(9):ITC65-ITC80. PubMed PMID: 27802475. doi: 10.7326/AITC201611010

9. Gordon LK, Goldman M, Sandusky H, Ziv N, Hoffman GS, Goodglick T, Goodglick L. Identification of candidate microbial sequences from inflammatory lesion of giant cell arteritis. Clin Immunol. 2004;111(3):286-96. PubMed PMID: 15183149. doi: 10.1016/j.clim.2003.12.016

10. Hoffman GS, Calabrese LH. Vasculitis: determinants of disease patterns. Nat Rev Rheumatol. 2014;10(8):454-62. PubMed PMID: 24934189. doi: 10.1038/nrrheum.2014.89

11. Clifford A, Hoffman GS. Evidence for a vascular microbiome and its role in vessel health and disease. Curr Opin Rheumatol. 2015;27(4):397-405. PubMed PMID: 26002032. doi: 10.1097/BOR.0000000000000184

12. Hoffman GS, Getz TM, Padmanabhan R, Villa-Forte A, Clifford AH, Funchain P, Sankunny M, Perry JD, Blandford A, Kosmorsky G, Lystad L, Calabrese LH, Eng C. The Microbiome of Temporal Arteries. Pathogens and Immunity. 2019;4(1):21-38. doi: 10.20411/pai.v4i1.270

13. Chakravarti R, Gupta K, Swain M, Willard B, Scholtz J, Svensson LG, Roselli EE, Pettersson G, Johnston DR, Soltesz EG, Yamashita M, Stuehr D, Daly TM, Hoffman GS. 14-3-3 in Thoracic Aortic Aneurysms: Identification of a Novel Autoantigen in Large Vessel Vasculitis. Arthritis Rheumatol. 2015;67(7):1913-21. PubMed PMID: 25917817. Pubmed Central PMCID: 4624269. doi: 10.1002/art.39130

14. Mukherjee PK, Wang H, Retuerto M, Zhang H, Burkey B, Ghannoum MA, Eng C. Bacteriome and mycobiome associations in oral tongue cancer. Oncotarget. 2017;8(57):97273-89. PubMed PMID: 29228609. Pubmed Central PMCID: 5722561. doi: 10.18632/oncotarget.21921

15. Albanese D, Fontana P, De Filippo C, Cavalieri D, Donati C. MICCA: a complete and accurate software for taxonomic profiling of metagenomic data. Sci Rep. 2015;5:9743. PubMed PMID: 25988396. Pubmed Central PMCID: 4649890. doi: 10.1038/srep09743

16. Magoc T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27(21):2957-63. PubMed PMID: 21903629. Pubmed Central PMCID: 3198573. doi: 10.1093/bioinformatics/btr507

17. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335-6. PubMed PMID: 20383131. Pubmed Central PMCID: 3156573. doi: 10.1038/nmeth.f.303

18. Rognes T, Flouri T, Nichols B, Quince C, Mahe F. VSEARCH: a versatile open source tool for metagenomics. PeerJ. 2016;4:e2584. PubMed PMID: 27781170. Pubmed Central PMCID: 5075697. doi: 10.7717/peerj.2584

19. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73(16):5261-7. PubMed PMID: 17586664. Pubmed Central PMCID: 1950982. doi: 10.1128/AEM.00062-07

20. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792-7. PubMed PMID: 15034147. Pubmed Central PMCID: 390337. doi: 10.1093/nar/gkh340

21. Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004;5:113. PubMed PMID: 15318951. Pubmed Central PMCID: PMC517706. doi: 10.1186/1471-2105-5-113

22. DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72(7):5069-72. PubMed PMID: 16820507. Pubmed Central PMCID: 1489311. doi: 10.1128/AEM.03006-05

23. Price MN, Dehal PS, Arkin AP. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009;26(7):1641-50. PubMed PMID: 19377059. Pubmed Central PMCID: 2693737. doi: 10.1093/molbev/msp077

24. Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, Clemente JC, Burkepile DE, Vega Thurber RL, Knight R, Beiko RG, Huttenhower C. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31(9):814-21. PubMed PMID: 23975157. Pubmed Central PMCID: 3819121. doi: 10.1038/nbt.2676

25. Parks DH, Beiko RG. Identifying biologically relevant differences between metagenomic communities. Bioinformatics. 2010;26(6):715-21. PubMed PMID: 20130030. doi: 10.1093/bioinformatics/btq041

26. Elling P, Olsson AT, Elling H. Synchronous variations of the incidence of temporal arteritis and polymyalgia rheumatica in different regions of Denmark; association with epidemics of Mycoplasma pneumoniae infection. J Rheumatol. 1996;23(1):112-9. PubMed PMID: 8838518.

27. Wagner AD, Gerard HC, Fresemann T, Schmidt WA, Gromnica-Ihle E, Hudson AP, Zeidler H. Detection of Chlamydia pneumoniae in giant cell vasculitis and correlation with the topographic arrangement of tissue-infiltrating dendritic cells. Arthritis Rheum. 2000;43(7):1543-51. PubMed PMID: 10902759. doi: 10.1002/1529-0131(200007)43:7<1543::AID-ANR19>3.0.CO;2-8

28. Koening CL, Katz BJ, Hernandez-Rodriguez J, Corbera-Bellata M, Cid MC, Schweizer HP, Li D, Kaplan J, Hoffman GS, De Domenico I. Identification of a Burkholderia-Like Strain From Temporal Arteries of Subjects with Giant Cell Arteritis. 2012 ACR/ARHP Annual Meeting, 2012.

29. Bhatt AS, Manzo VE, Pedamallu CS, Duke F, Cai D, Bienfang DC, Padera RF, Meyerson M, Docken WP. In search of a candidate pathogen for giant cell arteritis: sequencing-based characterization of the giant cell arteritis microbiome. Arthritis Rheumatol. 2014;66(7):1939-44. PubMed PMID: 24644069. Pubmed Central PMCID: 4113339. doi: 10.1002/art.38631

30. Nishida A, Inoue R, Inatomi O, Bamba S, Naito Y, Andoh A. Gut microbiota in the pathogenesis of inflammatory bowel disease. Clin J Gastroenterol. 2018;11(1):1-10. PubMed PMID: 29285689. doi: 10.1007/s12328-017-0813-5

31. Scher JU, Joshua V, Artacho A, Abdollahi-Roodsaz S, Ockinger J, Kullberg S, Skold M, Eklund A, Grunewald J, Clemente JC, Ubeda C, Segal LN, Catrina AI. The lung microbiota in early rheumatoid arthritis and autoimmunity. Microbiome. 2016;4(1):60. PubMed PMID: 27855721. Pubmed Central PMCID: 5114783. doi: 10.1186/s40168-016-0206-x

32. Allen-Vercoe E, Jobin C. Fusobacterium and Enterobacteriaceae: important players for CRC? Immunol Lett. 2014;162(2 Pt A):54-61. PubMed PMID: 24972311. Pubmed Central PMCID: 4259803. doi: 10.1016/j.imlet.2014.05.014

33. Larsen JM. The immune response to Prevotella bacteria in chronic inflammatory disease. Immunology. 2017;151(4):363-74. PubMed PMID: 28542929. Pubmed Central PMCID: 5506432. doi: 10.1111/imm.12760

34. Larsen JM, Musavian HS, Butt TM, Ingvorsen C, Thysen AH, Brix S. Chronic obstructive pulmonary disease and asthma-associated Proteobacteria, but not commensal Prevotella spp., promote Toll-like receptor 2-independent lung inflammation and pathology. Immunology. 2015;144(2):333-42. PubMed PMID: 25179236. Pubmed Central PMCID: 4298427. doi: 10.1111/imm.12376

35. Debarry J, Garn H, Hanuszkiewicz A, Dickgreber N, Blumer N, von Mutius E, Bufe A, Gatermann S, Renz H, Holst O, Heine H. Acinetobacter lwoffii and Lactococcus lactis strains isolated from farm cowsheds possess strong allergy-protective properties. J Allergy Clin Immunol. 2007;119(6):1514-21. PubMed PMID: 17481709. doi: 10.1016/j.jaci.2007.03.023

36. Schaller M, Loewenstein M, Borelli C, Jacob K, Vogeser M, Burgdorf WH, Plewig G. Induction of a chemoattractive proinflammatory cytokine response after stimulation of keratinocytes with Propionibacterium acnes and coproporphyrin III. Br J Dermatol. 2005;153(1):66-71. PubMed PMID: 16029328. doi: 10.1111/j.1365-2133.2005.06530.x

37. de la Fuente R, Schleifer KH, Götz F, Köst H-P. Accumulation of porphyrins and pyrrole pigments by Staphylococcus aureus ssp. anaerobius and its aerobic mutant. FEMS Microbiology Letters. 1986;35(2‐3):183-8. doi: 10.1111/j.1574-6968.1986.tb01524.x

38. McClary JS, Boehm AB. Transcriptional Response of Staphylococcus aureus to Sunlight in Oxic and Anoxic Conditions. Front Microbiol. 2018;9:249. PubMed PMID: 29599752. Pubmed Central PMCID: 5863498. doi: 10.3389/fmicb.2018.00249

39. Marques da Silva R, Caugant DA, Eribe ER, Aas JA, Lingaas PS, Geiran O, Tronstad L, Olsen I. Bacterial diversity in aortic aneurysms determined by 16S ribosomal RNA gene analysis. J Vasc Surg. 2006;44(5):1055-60. PubMed PMID: 17098542. doi: 10.1016/j.jvs.2006.07.021

40. Ramirez JA. Isolation of Chlamydia pneumoniae from the coronary artery of a patient with coronary atherosclerosis. The Chlamydia pneumoniae/Atherosclerosis Study Group. Ann Intern Med. 1996;125(12):979-82. PubMed PMID: 8967709.

41. Mosorin M, Surcel HM, Laurila A, Lehtinen M, Karttunen R, Juvonen J, Paavonen J, Morrison RP, Saikku P, Juvonen T. Detection of Chlamydia pneumoniae-reactive T lymphocytes in human atherosclerotic plaques of carotid artery. Arterioscler Thromb Vasc Biol. 2000;20(4):1061-7. PubMed PMID: 10764674.

42. Choi J, Chung SW, Kim SJ, Kim SJ. Establishment of Porphyromonas gingivalis-specific T-cell lines from atherosclerosis patients. Oral Microbiol Immunol. 2001;16(5):316-8. PubMed PMID: 11555310.

43. Lalla E, Lamster IB, Hofmann MA, Bucciarelli L, Jerud AP, Tucker S, Lu Y, Papapanou PN, Schmidt AM. Oral infection with a periodontal pathogen accelerates early atherosclerosis in apolipoprotein E-null mice. Arterioscler Thromb Vasc Biol. 2003;23(8):1405-11. PubMed PMID: 12816879. doi: 10.1161/01.ATV.0000082462.26258.FE

44. Mougeot JC, Stevens CB, Paster BJ, Brennan MT, Lockhart PB, Mougeot FK. Porphyromonas gingivalis is the most abundant species detected in coronary and femoral arteries. J Oral Microbiol. 2017;9(1):1281562. PubMed PMID: 28326156. Pubmed Central PMCID: 5328378. doi: 10.1080/20002297.2017.1281562

Footnotes