Association between smoking history and optical coherence tomography angiography findings in diabetic patients without diabetic retinopathy

Purpose To investigate any associations between cigarette smoking and retinal microvascular changes in diabetic patients without visible retinopathy. Design Retrospective, cross-sectional study. Participants 1099 eyes from 1099 diabetic patients with no clinical evidence of diabetic retinopathy (DR) were included in this study. Methods Diabetic patients underwent optical coherence tomography angiography (OCTA) scanning at Zuckerberg San Francisco General Hospital and Trauma Center between April 2018 and September 2019. Patient demographic and clinical information was collected. Standard bivariate statistics and multivariate linear regression were performed. Main outcome measures OCTA parameters included metrics related to the foveal avascular zone (FAZ; area, perimeter, circularity), perfusion density (PD; full, center, inner), and vessel length density (VLD; full, center, inner). Results The study population included 750 non-smokers and 349 smokers. FAZ perimeter was the only OCTA parameter that was significantly different between the two groups on uncontrolled analysis (P = 0.033). Multivariate regression analyses revealed significant associations between lower VLD full (β = -0.31, P = 0.048), lower VLD inner (β = -0.35, P = 0.046) and a history of smoking. No significant associations between cigarette smoking and either FAZ or PD were detected. Conclusions Our results suggest that smoking is likely associated with deleterious changes in the retinal microvasculature of patients with a history of diabetes and no visible DR. Based on these findings, diabetic patients with a history of smoking may benefit from higher prioritization in terms of ophthalmic screening.


Introduction
Diabetes mellitus (DM) leads to microvascular and macrovascular complications that place a heavy burden on societies around the world [1]. One such complication is diabetic retinopathy (DR), which is becoming the leading cause of preventable blindness in working-age populations [2]. DR is a complex, chronic disease caused by interactions between genetic factors and environmental conditions. Implicated risk factors for DR development and progression include age, gender, race, duration of diabetes, hyperglycemia, hypertension, hyperlipidemia, body mass index, and insulin use [3][4][5]. However, the relationship between smoking and DR remains unclear despite the known association between smoking and other vasculopathic diseases. The epidemiologic literature has reported conflicting results on the specific association of smoking with DR, with some studies finding a higher risk of DR [6] and others showing no effect or even protection against DR [7][8][9][10].
For many years, fluorescein angiography (FA) has been considered the gold standard for assessing microvascular changes in retinal vascular diseases such as DR. However, optical coherence tomography angiography (OCTA) has recently emerged as a noninvasive alternative technology. It can visualize blood flow in retinal blood vessels and determine the size of capillaries by detecting the contrast of blood cells in motion. Compared with conventional FA, it is a faster and safer imaging technique with high contrast and resolution [11][12][13][14]. Recent work has found widefield swept-source OCTA to be comparable to ultra-widefield fluorescein angiography in the detection of DR lesions [15]. In addition, OCTA can detect preclinical retinal microvascular abnormalities before the appearance of visible microaneurysms [16].
Given the ability of OCTA to detect early diabetes-related vascular changes in the retina, this technology offers a unique opportunity to study associations between smoking and incipient DR development. Prior epidemiologic studies have depended upon the visualization of later-stage clinical manifestations [6][7][8][9][10]. In this setting, the cumulative contribution of other, more determinative, risk factors might have rendered a true association with smoking difficult to confirm. In contrast, a more sensitive OCTA-based approach in patients with no visible DR may be able to identify smaller sub-clinical changes.
At this time, there is a paucity of studies that have examined the relationship between smoking and OCTA-based vascular parameters. One report compared OCTA findings between a small group of smokers and non-smokers and found no differences in central macular perfusion; however, choroidal thickness was larger in the group of smokers [17]. In addition, studies that focus on acute changes immediately after smoking have shown mixed results [18,19]. As such, the relationship between smoking and microvascular changes detected by OCTA remains unclear. The evaluation of patients with DM and no visible DR using OCTA offers an opportunity to quantify discrete changes associated with smoking. Such changes may confer a higher risk of progression to DR, given that microvascular alterations on OCTA can be predictive of subsequent retinopathy risk [20]. For these reasons, we sought to test the hypothesis that smoking is associated with pathologic OCTA-based vascular changes in patients with DM and no visible DR.

Participants
This retrospective cross-sectional study involved patients evaluated through a DR screening program at Zuckerberg San Francisco General Hospital and Trauma Center from April 2018 to September 2019. This study was approved by the Human Research Protection Program at the University of California, San Francisco and followed the tenets of the Declaration of Helsinki.
Adult (18 years or older) patients with type 1 or type 2 DM presenting for DR screening underwent ultra-widefield fundus photography (Optos Daytona, Optos PLC, Dunfermline, UK) and OCTA imaging during the visit. All eyes without detectable DR on the ultra-widefield fundus images were considered for inclusion in the study. Eyes with a history of ocular trauma or any form of ocular disease, with the exception of cataract, were excluded.

Clinical data
Collected information included age, sex, race, duration of diabetes mellitus, most recent hemoglobin A1c level (HbA1c), insulin use, body mass index (BMI), hypertension, hyperlipidemia, and smoking history. Smoking history was self-reported by the patient through solicitation at the DR screening visit [21,22]. The smoking history was used to classify patients into two categories. Patients who had never smoked or had smoked fewer than 100 cigarettes in their lifetime [23] were labeled as non-smokers. Patients who had smoked greater than 100 cigarettes in their lifetime irrespective of current smoking status were labeled as smokers.

OCTA image acquisition
OCTA imaging was obtained using a Cirrus™ HD-OCT 5000 with AngioPlex OCT Angiography (Carl Zeiss Meditec, Dublin, CA, USA) machine. A 3 x 3 mm scan area centered on the fovea was obtained in the right eye and then the left eye of each patient. Each scan included 245 clusters of B-scans repeated at least twice for repeatability analysis, and each B-scan consisted of 245 A-scans. FastTrac motion-correction software was used during image acquisition to reduce artifact. Images were considered to be of satisfactory quality by the following criteria: signal strength (SS) of greater than or equal to 7, maintenance of the fovea within the center of the scanned image without significant decentration, minimal blinking, and no more than one motion artifact.

OCTA image processing
The following parameters were generated automatically by the built-in OCTA analysis software (CIRRUS 11.0; clearance by the US Food and Drug Administration pending): foveal avascular zone (FAZ), perfusion density (PD), and vessel length density (VLD) in the superficial capillary plexus (SCP). The SCP was defined by an inner boundary at the internal limiting membrane and an outer boundary at the junction of the inner plexiform layer and the inner nuclear layer. For each patient, the images for each eye were assessed for artifactual errors and the eye with fewer errors was selected for analysis. If such errors were similar, the eye with the higher SS was selected. If both quality metrics were similar, the right eye was selected.
The FAZ (Figs 1A, 1B and 2A) reported by the software includes three indicators: FAZ area, FAZ perimeter and FAZ circularity. The FAZ area was the two-dimensional area of the FAZ, measured in mm 2 ( Fig 2B). The FAZ perimeter was the length of the line around the FAZ, measured in mm ( Fig 2C). FAZ circularity was defined as 4π x FAZ area/(FAZ perimeter) 2 and quantifies the degree of irregularity of the shape of the FAZ on a scale from 0 to 1. PD (Fig 1C and 1D) was defined as the percentage of area comprised of blood vessels within the total selected area. VLD (Fig 1E and 1F) was defined as the total length of vasculature per unit area. Both PD and VLD were evaluated within specific regions: center (central circle with 1-mm diameter) (Fig 2E), inner (ring around the center with inner and outer radii of 0.5 and 1.5 mm) (Fig 2F), and full (both the central circle and the inner ring) ( Fig 2D). The data reported by the software did not include an outer ring since images were limited to a 3 x 3 mm scan area.

Statistical analysis
Qualitative data are expressed in frequencies or percentages. Quantitative data are presented as a mean ± standard deviation (SD). The normality of continuous variable distributions was tested with a Kolmogorov-Smirnov test. Quantitative variables were compared between the two groups using a t test or Mann-Whitney test if they were normally or non-normally distributed, respectively. Chi-squared testing was used to compare categorical variables. Multivariate linear regression was performed using OCTA parameters as dependent variables and the following clinical data factors as independent variables: age, sex, race, smoking history, diabetesrelated parameters (duration of disease, last hemoglobin A1c, insulin use), and systemic factors (hypertension, hyperlipidemia, BMI). A p-value < 0.05 was considered to be statistically significant. Corrections for multiplicity were not made due to the high degree of correlation among the OCTA parameters. All statistics were performed using GraphPad Prism software (version 8.3.1; GraphPad Inc., San Diego, California, USA).

Results
1099 eyes from 1099 patients with DM were included in this study and were comprised of 750 non-smokers and 349 smokers. The DM sub-type composition was as follows: type 2 DM (n = 1085), type 1 DM (n = 9), un-specified (n = 5). Table 1 summarizes the demographic and clinical characteristics of all patients stratified by smoking history. Gender (P < 0.0001) and race (P < 0.0001) compositions were significantly different between the two groups, with a larger proportion of women in the never smoker group and a larger proportion of non-Hispanic Black and other race patients in the current smoker group. Supporting information is available online in S1 Table. Results of uncontrolled analyses are reported in Table 2. FAZ perimeter was the only OCTA parameter that was significantly different between the two groups (P = 0.033). However, FAZ area was trending towards significance (P = 0.070).
Select significant results from multivariate linear regression models with PD full, PD center, PD inner, VLD full, VLD center, and VLD inner as dependent variables are shown in Table 3. There were significant associations between lower VLD full (β = -0.31, P = 0.048), lower VLD

PLOS ONE
Association of smoking with optical coherence tomography angiography in diabetic patients

PLOS ONE
Association of smoking with optical coherence tomography angiography in diabetic patients inner (β = -0.35, P = 0.046) and smoking history. Larger FAZ area was associated with female sex (β = -0.05, P < 0.0001) and race. Longer FAZ perimeter was also associated with female sex (β = -0.17, P < 0.0001) and race. Higher FAZ circularity was associated with older age (β = -0.24, P < 0.0001). No significant associations between cigarette smoking and either FAZ or PD were detected. In addition, no significant associations between duration of DM or most recent HbA1c and either FAZ or PD or VLD were detected.

Discussion
In this study, we evaluated the association of cigarette smoking with microvascular changes in the macular region of diabetic patients without visible retinopathy using OCTA. Smoking status is a categorical variable consisting of the following categories: Non-smoking and smoking.
The following independent variables were used in each multivariate linear regression model: Age, sex, race, smoking history, diabetes-related parameters (duration of disease, last hemoglobin A1c, insulin use), and systemic factors (hypertension, hyperlipidemia, BMI). � = P < 0.05; † = P < 0.01; ‡ = P < 0.001. https://doi.org/10.1371/journal.pone.0253928.t003 Microvascular abnormalities are considered to be major pathological changes in diabetic retinopathy [24]. In the early stages of diabetes, it is increasingly recognized that microvascular abnormalities can progressively affect retinal function [25]. In this study, we detected a negative association between a history of smoking and microvascular density, that is, diabetic patients with a history of smoking who had no retinopathy had lower full and inner ring region VLD values. Interestingly, no associations between smoking and the central region PD and FAZ parameters were detected. The explanation for this discrepancy is unclear; however, it may be that the extrafoveal microvasculature is more susceptible to damage from smoking. Nonetheless, these findings suggest that smoking may have a detrimental impact on the retinal microvasculature. By extension, smoking may accelerate DR-related vascular alterations at an early stage in the arc of retinopathy development and progression. These results are supported by the findings of Lee et al., [22] who found that current smoking was a risk factor for reduced vascular density on multivariate analysis in a small group of Korean patients who had DM for an average of three years and had no DR. The current study builds upon those findings in a large and diverse population who, on average, had diabetes for twice as long; the ability to discern the association between smoking and microvascular changes in the face of additional accumulated diabetes-related damage strengthens the attribution of these changes to smokingrelated mechanisms.
Uncontrolled analysis showed that only FAZ perimeter was significantly different between the two groups in our study; however, controlled analysis revealed a larger number of significant results. This discordance can be explained by inter-group demographic differences. For example, the majority of smokers were men. Previous studies have shown that women have lower macular vascular density [26,27] and a larger FAZ [28,29] compared to men in younger patient populations. However, in contrast, studies in older adults have demonstrated that men have a lower vascular perfusion density [30,31]. Further complicating the situation is the fact that yet other demographic or clinical factors are also likely associated with vascular parameter measurements determined by OCTA. Ultimately, these considerations underscore the importance of controlling for confounding in this setting.
Interestingly, no associations were found between diabetes-related variables, such as disease duration and most recent HbA1c, and any of the studied OCTA parameters. This result is somewhat counterintuitive since these variables have been implicated in diabetic retinopathy incidence and/or progression in prior studies [3][4][5][6][7]. However, this result may be explained by the pre-clinical phenotype of the patients that were included in the present study. In this population, such variables may not have a significant association with microvascular changes that are only detectable with sensitive imaging modalities such as OCTA.
Older age was also found to be associated with lower PD and VLD on controlled analysis, which has been reported in previous studies [32][33][34]. Iafe et al. [32] reported that SCP and deep capillary plexus (DCP) vessel density decreased with age while FAZ area increased with age in normal subjects. Of note, our study found that the beta coefficients for age were larger compared to smoking history in the VLD full but not the VLD inner regression models. Our study population included diabetic patients without retinopathy, who likely differ from nondiabetic subjects in important ways. Therefore, while the contribution of older age to macular microvascular changes may be relatively large in our study, the negative association between smoking and retinal microvasculature changes should not be ignored.
Of note, FAZ characteristics may not be as useful as PD and VLD in assessing retinal microvascular changes due to the significant variation in [35,36], and ambiguous segmentation of, boundaries between the SCP and DCP [37,38]. In fact, the large variation in FAZ measurements may have contributed to the lack of significant findings with regards to the center parameters for both PD and VLD.
A key component of the prevention efforts targeting DR is appropriate screening. Early screening enables timely intervention, which can delay the development of vision loss due to the sequelae of DR in many patients. However, in resource-limited contexts, it is often not possible to screen all diabetic patients for evidence of DR [39]. In this setting, risk stratification can inform optimal referral patterns. According to our results, a history of smoking may be associated with deleterious changes in the retinal microvasculature of diabetic patients. As such, this risk factor should be considered when triaging patients from primary care clinics for a formal ophthalmic evaluation.
Limitations of this study include its cross-sectional design, which precludes any inferences regarding the predictive value of smoking on DR progression. An additional limitation is the significant demographic differences between the three studied groups. It is not unexpected that different age, sex, and race categories have different rates of smoking, and this issue was addressed with multivariate regression; however, we cannot rule out the possibility of residual confounding affecting our results. In addition, it should be noted that the imaging data used in this study was derived from the SCP in a 3 x 3 mm scan area using a spectral domain device. The relationship between smoking and microvasculature changes may be different in deeper layers of the retina, the choroid, or the non-central macula. Future studies should include analyses of the DCP, choriocapillaris, and choroid derived from a larger scan area using a sweptsource device. Indeed, widefield swept-source OCTA can be particularly valuable in the detection of peripheral lesions in DR [40]. Lastly, smoking history was classified as a binary variable in our study. In future work, this exposure should be quantified (i.e. pack-years) rather than categorized in order to assess correlations including potential dose-response relationships.

Conclusions
In conclusion, our results show that there is an association between a history of smoking and a decrease in VLD in diabetic patients without visible retinopathy. Given these results, diabetic patients with a history of smoking may benefit from higher prioritization in terms of ophthalmic screening.
Supporting information S1