Association Between Gut Microbial Abundance and Sight-Threatening Diabetic Retinopathy

Abstract

Purpose: To study the association between gut microbial abundance and sight-threatening diabetic retinopathy among patients with a history of type 2 diabetes mellitus.

Methods: An observational case-control study was performed using a sample population of diabetics referred to a tertiary eye institute. Sample subjects were identified as cases if they were diagnosed with sight-threatening diabetic retinopathy and controls if they were not but had at least a 10-year history of diabetes. Fecal swabs for all patients were collected for enumeration and identification of sequenced gut microbes. Statistical analyses were performed to associate the clinically relevant Bacteroidetes to Firmicutes relative abundance ratio (B/F ratio) with sight-threatening diabetic retinopathy and an optimal cutoff value for the ratio was identified using Youden's J statistics.

Results: A sample size of 58 diabetic patients was selected (37 cases, 21 controls). No statistically significant difference in the relative abundance among the predominant phyla between the groups were found. In our univariate analysis, the B/F ratio was elevated in cases compared to controls (cases, 1.45; controls, 0.94; P = 0.049). However, this statistically significant difference was not seen in our multivariate regression model. Optimal cutoff value of 1.05 for the B/F ratio was identified, and significant clustering of cases above this value was noted in beta diversity plotting.

Conclusions: No difference in gut microbial abundance for any particular phylum was noted between the control and diseased population. Increased gut microbial B/F ratio can be a potential biomarker for the development of sight-threatening diabetic retinopathy among type 2 diabetic patients.

Diabetes mellitus (DM) is one of the world's fastest growing metabolic diseases.¹ In 2017, 451 million individuals were assessed to have diabetes, more than 90% of whom had type 2 DM, with an extended increment to 693 million by 2045.² Diabetic retinopathy (DR), a vision-threatening microvascular complications of DM, involves an interplay of angiogenic and inflammatory mediators in its pathogenesis.³ It progresses from mild nonproliferative diabetic retinopathy, characterized by increased vascular permeability, to moderate and severe nonproliferative diabetic retinopathy, characterized by vascular closure, and ultimately proliferative diabetic retinopathy (PDR).⁴ Sight-threatening diabetic retinopathy (STDR), which primarily comprises clinically significant macular edema (CSME), and PDR were observed in 14% and 23% of younger-onset patients, respectively. Although nonproliferative DR does not necessitate medical intervention besides close medical follow-up and strict glycemic management, these sight-threatening conditions that can present independently or together can result in severe visual loss if left untreated.^5,6 Proper management can prevent more than 90% of cases of visual loss in patients suffering from PDR, diabetic macular edema (DME) or both.⁷ The established risk factors for the occurrence and progression of DR includes longer duration of DM, obesity, poor glycemic control, altered lipid profile, and hypertension.⁸ DM causes various metabolic and physiological variations from the norm in the retina; however, which of these abnormalities contributes to perceived anatomical changes that are seen in DR is less clearly understood. 

Microbiome alludes to the multispecies network of microorganisms (bacteria, fungi, and viruses) that dwell within and over different parts of the human body varying in abundance and composition.^9,10 It has been noted that aberrant microbial composition relating to the relative abundance of the two most abundant gut bacteria namely Bacteroidetes and Firmicutes is associated with gut dysbiosis through a process of disruption in the gut/vascular barrier. This disruption consequently leads to pathological translocation of gut bacteria and its metabolic byproducts inducing systemic inflammation. A recently performed animal model study did not note any significant difference between healthy individuals and patients with ocular manifestations in the form of DR with respect to Firmicutes/Bacteroidetes ratio.¹¹ Further studies have also shown that type 2 diabetes, obesity, and high glycemic states are related with a modified gut microbiota resulting in some level of gut dysbiosis.^11–16 The likely cause of this association as described by previous studies is likely linked to the gut microbiota's involvement in physiological functions such as energy metabolism, metabolic signaling, and regulation of integrity of the gut barrier.^17–19 These variables also have a role in the presence and progression of DR, because DR accounts for more than 60% incidence in type 2 DM.²⁰ In humans, there is only one report by Moubayed et al.¹⁴ studying the influence of gut dysbiosis based on abundance restricted to a small cohort from Saudi Arabia; they observed higher Bacteroidetes ratio in diabetic group, and no significant difference of bacterial strains in subjects with or without DR. No similar such clinical studies have been conducted to identify the association of alterations in the gut microbiome with respect to DR among a South Asian cohort. 

In our present study, we studied the gut bacterial microbiomes of people with previously diagnosed type 2 DM from a primarily South Indian population to ascertain whether dysbiosis in the gut microbiome is associated with STDR. In addition, Bacteroidetes to Firmicutes ratio (B/F ratio) has been associated in previous literature as a potential diagnostic biomarker for diabetes,^11,14,21,22 and, hence, we wanted to assess its utility as a screening biomarker in STDR as well. 

Materials and Methods

Between April 2019 and October 2019, patients with pre-existing type 2 DM presenting to a tertiary referral eye care hospital to undergo ophthalmic evaluation were recruited for the study. All patients who met the above inclusion criteria underwent comprehensive ophthalmic evaluation including best-corrected visual acuity, slit-lamp biomicroscopy, and fundus examination and an exclusion criteria that included presence of any other vascular retinopathy other than DR, use of systemic antibiotic during the past six months, and any other systemic or ocular inflammatory or degenerative or neoplastic disorder that could independently affect the retina was applied. In total, 58 patients were recruited who were identified into two groups, the control group who were patients with at least a 10-year history of previously diagnosed type 2 DM on the basis of existing clinical criteria and no diagnosis of DR on eye examination and the cases group who were type 2 diabetics diagnosed with sight-threatening DR such CSME or PDR. 

The classification of sight-threatening DR was based on Early Treatment Diabetic Retinopathy Study Research Group and International clinical diabetic retinopathy and diabetic macular edema disease severity scales.^23,24 

• • CSME was defined as thickening of the retina at or within 500 µm of the center of the macula; hard exudate at or within 500 µm of the center of the macula associated with thickening of adjacent retina; or a zone of retinal thickening 1-disc area or larger, any part of which is within 1-disc diameter of the center of the macula
• • PDR is characterized by neovascularization of the disc, neovascularization of the retina, neovascularization of the iris, neovascularization of the angle, vitreous hemorrhage, or tractional retinal detachment. With regard to macular edema, it should be noted whether macular edema is present or absent.

Qualitative data regarding baseline demographic information including age, duration of DM, pre-existing comorbidities, antiglycemic medications, body mass index, HbA1c, fasting, and postprandial blood glucose levels, blood group type and dietary patterns were collected using patient questionnaire. 

The study protocol was approved by the Institutional Review Board at Vision Research Foundation, Sankara Nethralaya, Chennai. Written informed consent was obtained from all the studied subjects for sample collection and subsequent analyses, and it adhered to the tenets of the Declaration of Helsinki. 

Stool Collection and Next-Generation Sequencing

A clean catch of fecal sample was collected from each patient using fecal swabs (Norgen Biotek Corp, Thorold, Canada). Norgen's Microbiome DNA Isolation Kit was used to isolate the bacterial DNA from the stool sample. Purification was based on spin column chromatography. Initially, Lysis Additive A was added to the swab collection tube and incubated at 65˚C to efficiently and rapidly homogenize the sample. The sample was then spun in a centrifuge, and the supernatant is transferred to a DNAse-free microcentrifuge tube. Binding Buffer I was added, and the lysate was incubated for 10 minutes on ice. The lysate was then spun for two minutes to pellet any cell debris, the supernatant was collected, an equal volume of 70% ethanol was added to the lysate and the solution was loaded onto a spin-column. Norgen's spin column binds nucleic acids in a manner that depends on ionic concentrations, thus only the DNA will bind to the column while the proteins and other contaminants are removed in the flow through or retained on top of the resin. The bound DNA was then washed using the provided Binding Buffer B and Wash Solution A, and the purified DNA was eluted using the Elution Buffer B. 

Library Preparation and Next Generation Sequencing

The standard Illumina protocol was used to generate the microbiome library. In all sequencing runs V4 region of 16 rRNA gene were amplified. The primer pair used were 515F (5’-GTGCCAGCMGCCGCGGTAA-3′) and 806R (5’-GGACTACHVGGGTWTCTAAT-3′) targeting V4 region of the 16S rRNA gene. 

The 16s rRNA Bacterial Sequencing was conducted by Npedia technologies, in conjunction with National Centre for Cell Science, Pune on the stool samples and for this purpose Illumina next generation sequencing (NGS) Miseq platform was used with 250 × 2 paired end chemistry as their approved procedure. The raw output from Illumina MiSeq were demultiplexed to give FASTQ output that was used for downstream analysis. 

Data Analysis

Taxonomic Classification

The Illumina demultiplexed paired-end sequenced dataset was processed by the R package DADA2(1.14.1)²⁵ to correct for amplicon errors, to identify chimeras, and to merge paired-ends reads. The paired ends were filtered and trimmed at 220 for forward run and 160 for reverse run, on the basis of the quality profile of the reads. 

The end product of DADA2 yielded a total of 10,129 unique ASVs (amplicon sequence variants) and tallied as the number of times each exact ASV was observed for each sample. A phyloseq R object²⁶ comprised of all 10129 ASVs was generated, with a lookup table of taxonomy assignments to each ASVs obtained by using the SILVA v132 reference database, along with a phylogenetic tree and subject sample metadata. 

Statistical Analysis

A microbiome R package²⁷ and phyloseq R package were used to analyze the abundance, relative abundance, alpha diversity of samples. R 3.6.3 was used for all the analysis. All the R visualization was done using ggplot2 (v 3.3.2). Further statistical analyses were performed using a standard software package (Stata, version 16.1, StataCorp). Descriptive statistics on patients’ characteristics and summary statistics of gut microbiome relative abundance comparing controls and cases were performed. Two-sample t-testing was performed to compare the B/F abundance ratio between cases and controls. Multivariate linear regression analysis with adjusted cofactors decided a priori was performed to determine the significance of B/F ratio in differentiating cases from controls. Youden's J statistics method was used to identify an optimal B/F ratio cutoff point to differentiate between the two study groups. Bray-Curtis dissimilarity using principal co-ordinate analysis to differentiate clustering between cases and controls against this cutoff point was performed. 

Results

Baseline Characteristics of Cases and Controls

In our study cohort of 58 patients, based on eye examination 21 patients were found to be in the control group, and 37 were diagnosed with STDR. Of these 37 cases, 21 cases were found to be diagnosed as PDR with or without DME with a median disease duration of 0.50 (interquartile range, 0.50–1.00) years, and 16 cases were noted to be diagnosed as CSME alone with a median disease duration of 0.75 (interquartile range, 0.50–2.63) years. 

Table 1 illustrates the baseline characteristics of cases and controls. All of the variables—age, gender, vegetarian diet, body mass index, duration of DM, duration of hypertension, glycemic controls, Log MAR best-corrected visual acuity of the worse eye, use of metformin, and associated systemic diseases such and coronary artery disease and dyslipidemia—were similar in both of the groups; systemic hypertension, however, was present in higher proportion of subjects in cases than controls (70.27% vs. 42.86%; P = 0.04). 

Table 1.

 

Baseline Characteristics of Cases and Controls

 

Relative Abundance of Gut Microbiome in Cases and Controls

The distribution of various gut microbiome after genomic sequencing were identified into 17 different phyla in this study. When analyzing the microbial abundance for the study population as a whole, most sequences were affiliated predominantly with four phyla. Bacteroidetes occupied the largest portion (40.27% [35.3%–45.26%]) followed by Firmicutes (41.27% [37.38%–45.21%]), Proteobacteria (8.72% [5.92%–11.47%]), and Actinobacteria (8.00% [5.97%–10.03%]). 

Based on t-test statistics comparing control and cases (Table 2), a similar distribution within the individual groups was noted with the most common phyla being Bacteroidetes (cases, 43.38%; controls, 34.80%) and Firmicutes (cases, 39.75%; controls, 43.95%), and less common ones were Proteobacteria (cases, 8.53%; controls, 9.05%) and actinobacteria (cases, 6.63%; controls, 0.61%). However, no statistically significant difference in the distribution of these phyla between cases and controls were identified. 

Table 2.

 

Relative Abundance of Gut Microbiome in Cases and Controls

 

Table 3.

 

Multivariate Linear Regression Model Comparing the Association of B/F Ratio to Sight-Threatening Diabetic Retinopathy

 

Comparing B/F Ratio in Cases and Controls

When we compared the ratio of the most common identified gut phyla, Bacteroidetes (B) and Firmicutes (F), in cases and controls, the B/F ratio was significantly higher in cases than controls (cases, 1.45; controls, 0.94; P = 0.049) (Fig. 1a) was noted in our t-test analysis; however, no significant difference between the groups was seen in our linear regression model. Additional subgroup analysis exploring the diseased cohort individually as patients diagnosed with CSME alone or PDR (with or without DME) noted that both disease groups in univariate analysis had a higher B/F ratio compared to the control group. However, this difference was found to be statistically significant only in the CSME-alone group (CSME, 1.61; controls, 0.94; P = 0.014) and not statistically significant in the PDR group (PDR, 1.33; controls, 0.94; P = 0.180) (Fig. 1b). Multivariate linear regression models subgrouping the diseased cohort and analyzing this correlation once again did not reveal any significant difference against the control group (Table 3). 

Figure 1.

 

(a) Comparison of B/F ratio among cases and controls. (b) Comparison of B/F ratio among controls and cases subgroup as CSME or PDR.

B/F Ratio Optimum Cutoff Point

Youden's J statistics method, provided an optimal cutoff point of 1.05 (sensitivity, 62%; specificity, 57%), with controls being identified below this value and cases above. Difference in clustering within our sample population on the basis of taxonomic abundance profile using principle co-ordinate analysis, Bray-Curtis dissimilarity found a high percentage of cases clustered above this cutoff point (Fig. 2). 

Figure 2.

 

The Bray-Curtis principal coordinates analysis (PCoA) beta diversity plotting of cases versus controls against B/F ratio (1 = B/F ratio >1, 0 = B/F ratio < 1). The above PCoA plot groups microbial communities characterized by similar sets of 16S rRNA gene sequences. Cases have been labeled as triangles and the control group as circles along with blue coloration depicting patients with a B/F ratio >1 and red with a B/F ratio < 1. Within the PCoA plot, cases with a corresponding B/F ratio > 1 are grouped together primarily as depicted by the right circle. In similar fashion, majority of the control group with a corresponding B/F ratio < 1 tend to be grouped together based on gene sequencing as described by the left circle. This grouping based on gene sequencing is in line with our study findings of a positive correlation between elevated B/F ratio and patients with a diagnosis of STDR.

Discussion

Hippocrates, the father of modern medicine, suggested nearly more than 2000 years ago, “All diseases begin in the gut.” If not all diseases, at least many chronic metabolic diseases such as diabetes mellitus do.^13,28,29 The Human Microbiome Project (National Institute of Health) and Metagenomics of the Human Intestinal Tract (MetaHIT, the European Project) were initiated to characterize microbial communities from different parts of the body.^30–32 In our cross-sectional study, we observed that the most common phyla in the gut microbiome were Bacteroidetes (Gram-negative) and Firmicutes (Gram-positive), accounting for nearly 83% in cases and 79% in controls. The percentage value of Bacteroidetes abundance was noted to be higher, and the percentage value of Firmicutes abundance was noted to be lower in cases versus controls, respectively. Hence, we noted that the ratio of Bacteroidetes to Firmicutes was also found to be significantly higher in cases than controls. Although this difference did not provide statistical significance when accounting for multiple covariables, we speculate that this difference that we initially found could be due to the fact that gut dysbiosis, which is a multifactorial pathology, could be involved in chronic inflammation and retinal ischemia, both of which play an important role in the pathophysiology of STDR.^33–38 It was interesting to note that when subgroup analysis of different disease diagnosis within STDR (CSME vs. PDR) were performed, the same increased B/F ratio was found in both disease conditions against diabetics without retinopathy. However, statistical significance of this correlation was only noted among the CSME only group. Because these two STDR disease conditions have different pathologic mechanisms,^39,40 this variation in our findings might shed light into narrowing down which metabolic pathways might be contributing to the increased B/F ratio among those with STDR. 

Normally, the intestinal lining prevents the migration of microbes and their metabolites from the gut lumen to the bloodstream. However, a change in the intestinal milieu, intestinal dysbiosis, may deregulate the barrier effect of gut lining and cause leaky gut syndrome. Predominant gram-negative bacterial phyla such as Bacteroidetes releases bacterial endotoxin, lipopolysaccharides and triggers an innate or natural immunity and thereby proinflammatory pathways resulting in vascular dysfunction. A study done on a mouse with diabetes (db/db mice) showed its fecal bacterial composition (predominantly in Bacteroidetes and Firmicutes), presenting with impaired intestinal barrier function and replicating some of the features of diabetic retinopathy—acellular capillaries, activation of retina microglia, and infiltration of peripheral immune cells into the retina.⁴¹ 

Taking into consideration the previously researched significance of B/F ratio in systemic diseases such as diabetes, we wanted to explore its association with STDR in a clinical setting. In addition, on the basis of the observed potential association with STDR among diabetics from our study, we wanted to further explore the possibility of its role as a potential biomarker in the diagnosis of STDR.^11,14 We were able to formulate an optimal B/F ratio cutoff point of 1.05, above which we can predict with high accuracy the probability of STDR being present among type 2 diabetic patients. This was evidently noted in our beta diversity analysis where the plotted clustering above the 1.0 cutoff predominantly (86.67%) comprised patients diagnosed with STDR (15 patients in total, 13 cases, 2 controls), whereas the clustering below the 1.0 cutoff did not provide any significant discrimination between cases and controls. 

Although our study was able to identify a positive correlation between gut dysbiosis and developing STDR, further prospective studies comprising varied geographic populations and possibly a larger sample size would be required to strengthen this clinically relevant association. There are several questions yet to be answered with respect to gut microbiome, such as the influence of racial differences, diet, socioeconomic status, antibiotic use, role of antiglycemics, as well as prebiotic/probiotic supplementation. As our understanding of gut dysbiosis and its influences on health and disease expands, we will have a new direction to explore newer therapeutics to manage diabetes and its microvascular complications such as diabetic retinopathy. 

Conclusion

There is no significant difference in the relative abundance of individual phylum between patients with STDR and no retinopathy within patients with type 2 diabetes. Gut dysbiosis with respect to difference in B/F relative abundance ratio displays significant involvement when differentiating diabetic patients with or without STDR, and this B/F ratio can play a potential role as a biomarker in the diagnosis of STDR. 

Acknowledgments

The authors thank Vijay Vaidyanathan, Npedia Tech; Biohub Data Science Pvt. Ltd, Chennai, Tamil Nadu, India. 

Supported by Novartis health care Pvt. Ltd, India. 

Disclosure: R. Khan, None; A. Sharma, None; R. Ravikumar, None; A. Parekh, None; R. Srinivasan, None; R.J. George, None; R. Raman, None