Next Article in Journal
Imaging Considerations before and after Liver-Directed Locoregional Treatments for Metastatic Colorectal Cancer
Next Article in Special Issue
Validation of Inter-Reader Agreement/Consistency for Quantification of Ellipsoid Zone Integrity and Sub-RPE Compartmental Features Across Retinal Diseases
Previous Article in Journal
Ambulatory Risk Stratification for Worsening Heart Failure in Patients with Reduced and Preserved Ejection Fraction Using Diagnostic Parameters Available in Implantable Cardiac Monitors
Previous Article in Special Issue
An Optical Coherence Tomography-Based Measure as an Independent Estimate of Retinal Function in Retinitis Pigmentosa
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Association of Genetic Risk for Age-Related Macular Degeneration with Morphological Features of the Retinal Microvascular Network

1
Chair and Clinical Department of Ophthalmology, Faculty of Medical Sciences in Zabrze, Medical University of Silesia, 40-752 Katowice, Poland
2
Faculty of Biomedical Engineering, Silesian University of Technology, 41-800 Zabrze, Poland
3
Institute of Biomedical Engineering and Informatics, Technische Universität Ilmenau, 98693 Ilmenau, Germany
4
Genomed S.A., 02-971 Warszawa, Poland
*
Author to whom correspondence should be addressed.
Diagnostics 2024, 14(7), 770; https://doi.org/10.3390/diagnostics14070770
Submission received: 13 March 2024 / Revised: 28 March 2024 / Accepted: 2 April 2024 / Published: 5 April 2024

Abstract

:
Background: Age-related macular degeneration (AMD) is a multifactorial disease encompassing a complex interaction between aging, environmental risk factors, and genetic susceptibility. The study aimed to determine whether there is a relationship between the polygenic risk score (PRS) in patients with AMD and the characteristics of the retinal vascular network visualized by optical coherence tomography angiography (OCTA). Methods: 235 patients with AMD and 97 healthy controls were included. We used data from a previous AMD PRS study with the same group. The vascular features from different retina layers were compared between the control group and the patients with AMD. The association between features and PRS was then analyzed using univariate and multivariate approaches. Results: Significant differences between the control group and AMD patients were found in the vessel diameter distribution (variance: p = 0.0193, skewness: p = 0.0457) and fractal dimension distribution (mean: p = 0.0024, variance: p = 0.0123). Both univariate and multivariate analyses showed no direct and significant association between the characteristics of the vascular network and AMD PRS. Conclusions: The vascular features of the retina do not constitute a biomarker of the risk of AMD. We have not identified a genotype–phenotype relationship, and the expression of AMD-related genes is perhaps not associated with the characteristics of the retinal vascular network.

1. Introduction

Age-related macular degeneration (AMD) is one of the leading causes of legal blindness in the industrialized world. AMD is characterized by the accumulation of extracellular deposits, namely drusen, along with the progressive degeneration of photoreceptors and adjacent tissues. AMD is a multifactorial disease that encompasses a complex interplay between aging, environmental risk factors, and genetic susceptibility. Chronic inflammation, lipid deposition, oxidative stress, and impaired maintenance of the extracellular matrix are strongly involved in the pathogenesis of AMD [1,2,3,4,5,6].
The retinal pigment epithelium (RPE) plays a crucial role in the degradation of photoreceptor cells in the retina. The RPE engulfs and breaks down the photoreceptor cells through phagocytosis, converting them into lipofuscin within its lysosomes. However, over time, lipofuscin accumulation in the RPE cells can lead to toxic reactions and, eventually, cell death [7].
Several genetic factors, lipid metabolism, oxidative stress, and aging are involved in the development of AMD. Dysregulation of genetics accounts for 46–71% of disease contribution, with CFH and ARMS2/HTRA1 identified as the two most significant risk loci among the 103 AMD-associated loci identified to date [8]. The genetic contribution to the development of AMD has been determined through linkage studies and genome-wide association approaches, which have revealed risk loci for AMD and several candidate genes. These include CFH, C3, C2-CFB, CFI, and a region on chromosome 10 that includes HTRA1/LOC387715/ARMS2, CETP, LIPC, VEGF-A, and APOE [9]. Although the exact pathophysiological mechanisms underlying AMD are still being investigated, there is general agreement that AMD is a multifactorial disease resulting from the synergistic effect of modifiable factors (e.g., smoking, eating habits, physical activity) and non-modifiable factors (age and genotype) [10]. The polygenic risk score (PRS) is used in numerous studies and is intended to estimate the risk of disease in various populations [11,12,13,14,15].
A key factor in the resulting retinal degeneration is reactive oxygen species (ROS) production. Increased levels of reactive oxygen species lead to damaged mitochondrial DNA (mtDNA), resulting in dysfunctional mitochondria, and, consequently, mtDNA causes further damage to retinal tissue [3]. These ROS are generated mainly when RPE cells are exposed to blue light within the range of 390 to 550 nanometers. This exposure triggers degenerative and apoptotic reactions in various retinal cells, including RPE cells, photoreceptors, pericytes, and endothelial cells [16,17].
The occurrence of oxidative stress triggers a cascade of events, including inflammation, which subsequently releases angiogenic factors [2]. In ocular diseases, neovascularization can be triggered by abnormal conditions such as ischemia, hypoxia, inflammation, and genetic predisposition. The presence of these conditions can lead to the stimulation of new blood vessel growth in the retina or choroid as a response to reduced oxygen supply to the affected area [18]. Neovascular age-related macular degeneration (nvAMD) is a condition that requires frequent and long-term intravitreal injections to preserve vision with a significant treatment burden to maintain sight [19,20,21].
nvAMD is a multifactorial disease with numerous risk factors that are thought to play a role in its development. Non-modifiable risk factors include increased age, female sex, genetic factors, Caucasian race, and light iris color. Modifiable risk factors include smoking, increased BMI, alcohol intake, and dietary habits [4]. Anti-VEGF therapies have emerged and are widely used, such as ranibizumab (Lucentis®, Novartis, Basel, Switzerland), aflibercept (Eylea®, Bayer, Berlin, Germany), brolucizumab (Beovu®, Novartis, Basel, Switzerland), faricimab (Vabysmo®, Roche, Basel, Switzerland), and the off-label use of bevacizumab (Avastin®, Roche, Basel, Switzerland) [21,22].
As the population continues to age, the incidence cases of AMD in general and neovascular AMD are projected to increase [5]. A study by Johanna M. Colijn et al. predicts that, by 2040, the number of individuals in Europe with early AMD will range between 14.9 and 21.5 million, and for late AMD, between 3.9 and 4.8 million [23].
Angio OCT, also known as optical coherence tomography angiography (OCTA), is a non-invasive imaging technique that allows the visualization of the microvasculature of the retina and choroid. It provides detailed information about blood flow in these structures without injecting any contrast agents [17]. Angio OCT utilizes the principles of OCT imaging based on low-coherence interferometry. It measures the interference pattern of light waves reflected from different depths within the tissue to create high-resolution cross-sectional images [24]. By detecting the dynamic movement of red blood cells, OCTA provides a unique perspective on blood flow in the retina and choroid. Angio OCT extends the capabilities of traditional OCT by incorporating motion contrast imaging. This innovative imaging technique enables a comprehensive analysis and visualization of the vascular structure of macular neovascularization (MNV), offering valuable insights into both physiological and pathological blood flow [25].
The retinal superficial capillary plexus (SCP) comprises the top three layers, housing arteries and veins responsible for the blood circulation in the deeper retina. Below these superficial layers, the deep capillary plexus (DCP) contains intermediate and deep capillaries. Notably, the SCP features larger, more prominently visible blood vessels due to its proximity to the imaging source, while the DCP displays finer vessels [24]. The choriocapillaris forms a densely interconnected network of capillaries, and a healthy choriocapillaris presents a distinct grainy texture in OCTA images. In contrast, the choroidal layer contains larger vessels in comparison to the retinal vasculature. Nevertheless, these larger choroidal vessels are not as well visualized in OCTA images, due to factors such as imaging depth limitations, signal attenuation, shadowing effects, the high reflectivity of the layers above, and the specific characteristics of blood flow [26,27].
The aim of the present study is to analyze the association between the characteristics of the retinal vascular network and the estimated risk of AMD based on genotypic analysis expressed as a PRS. To our knowledge, no articles have been published so far on OCTA images and the PRS to assess the relations between the vascular pattern in the macula and its deviations from the norm and genetic risk of AMD.

2. Materials and Methods

2.1. Dataset

Participants in the study were recruited by the Chair and Clinical Department of Ophthalmology at the Faculty of Medical Sciences in Zabrze, Medical University of Silesia in Katowice. The study was conducted in accordance with ethical guidelines and received approval from the Ethics Committee of the Medical University of Silesia (Resolutions No KNW/0022/KB1/105/13 and BNW/NWN/0052/KB1/97/I/22/23), following the principles outlined in The Declaration of Helsinki. Written informed consent was obtained from all participants after being informed of the purpose and protocol of the study.
The cohort comprised 332 patients, including 235 diagnosed with AMD and 97 healthy controls without symptoms of retinal degeneration. Individuals over 50 years of age were invited to participate, with AMD patients recruited during detailed follow-up visits at the ophthalmology clinic, including advanced forms of degeneration (neovascular AMD, geographic atrophy, and subretinal fibrosis) and early and intermediate forms. The control group included individuals recruited during follow-up visits after cataract surgery or planned check-ups, where fundus examination did not reveal any signs of AMD. The cohort included 20 eyes of patients treated for glaucoma (intraocular pressure was well controlled on one type of antiglaucoma drops): 8 of them were in the control group, and 12 were in the study group. The treatment received was similar in both groups.
Exclusion criteria for the study encompassed patients with macular diseases other than AMD (e.g., diabetic macular edema, macular dystrophy), limited vision that hinders fundus evaluation, prior retinal or choroidal inflammatory diseases, retinal detachment, intraocular procedures other than cataract surgery or posterior capsulotomy, and retinal laser therapy.
The study involved various ophthalmological examinations, such as best-corrected visual acuity (BCVA) using ETDRS charts, anterior segment examination with slit lamp biomicroscopy, and pupillary dilation with 1% tropicamide. A clinical examination of each eye using Volk superfield aspheric lenses 90D was performed to assess the posterior segments and further diagnostics to identify pathological changes in the macula. Digital fundus images, swept-source optical coherence tomography (SS-OCT) of the macula (radial and 3D wide scanning protocols), and optical coherence tomography angiography (OCTA) were obtained using the DRI OCT Triton tomograph from Topcon Healthcare, Tokyo, Japan. In this study, we used macular 6 mm × 6 mm blood flow OCTA images. Four en-face reconstructions were made based on segmentations of retinal layers in ranges defining the superficial capillary plexus (SCP), deep capillary plexus (DCP), outer retina (OR), and choriocapillaris (CC).

2.2. PRS Calculation

Ulańczyk et al. conducted research on 30 genes associated with critical retinal functions, including regulation of inflammation (e.g., TGFB1), immune response (e.g., C2, C3, CFB, CFH), lipid and protein synthesis (e.g., ELOVL4, HTRA1, RPL1) and the structure of the retinal layer (e.g., BEST1, C1QTNF5, GUCA1B). Other genes impact oxidative stress, extracellular matrix maintenance, transmembrane transport, transcription regulation, DNA repair, and AMD (e.g., ARMS2). Molecular Inverted Probes were used to enrich coding regions with flanking sequences. Subsequently, libraries were sequenced using the Illumina platform. The targeted enrichment of coding sequences from the 30 AMD-associated genes was achieved through this approach. The entire process took place at Genomed S.A., Warsaw [14].
The bioinformatic analysis included adapter trimming using Cutadapt v1.14 [28], reads mapping to the GRCh37.13 reference genome with Burrows–Wheeler Aligner v0.7.10, and deduplication based on Unique Molecular Identifiers using in-house scripts. For indel realignment and base recalibration, researchers used Genome Analysis Toolkit v3.5 (GATK) [29], a widely recognized and powerful software tool for best practices in these specific areas, ensuring the highest quality and accuracy of the results obtained.
The variant calling utilized HaplotypeCaller and UnifiedGenotyper from the GATK package to optimize single nucleotide variant and indel identification. Criteria included excluding variants missing in <95% of samples for data completeness and filtering out low-coverage variants (<10× in  80 %  of genotypes) to enhance reliability and accuracy.
In this study, data from a subgroup of subjects, previously considered by Wąsowska et al., were used in the PRS calculation using PLINK for QC of genotyped data and the additive model available via PRSice2 software (version 2.3.5). The PRS study was performed on the Polish population based on the targeted sequencing data involving the enrichment of genes known to be associated with AMD at the time of the study [30].

2.3. OCTA Image Selection

Figure 1 shows the flow chart of OCTA image selection based on objective exclusion criteria and experts’ evaluation. The best-quality images were selected from a database of 1255 reconstructed OCTA reports. Images with artifacts that distorted the vascular pattern in the OCTA image or showed poor scan quality were rejected. We avoided using images with poor resolution, high levels of noise, or significant artifacts that can detract from the clarity and interpretability of the content. Over-processing was also taken into consideration; it was important to avoid excessive image processing that may lead to the introduction of artifacts, loss of detail, or misrepresentation of the original data. Image processing was abandoned due to its substantial influence on quantitative features of the vessels’ network morphology.
Motion or doubling artifacts, blurry images, quality scores less than 40, or poor centration were the first stage exclusion criteria for the analysis. The following stages of image rejection were carried out based on the results of the segmentation algorithm, using a possible change in algorithm parameters to improve the results. Moreover, only those patients for whom correct OCTA images were found for the right and left eyes were included in further analysis. As a result of the selection, the accurate OCTA images of 132 patients were selected, allowing proper segmentation of the vessels in the SCP layer. In addition, a subgroup of 73 patients with nvAMD in whom MNV was evident on CC or OR layers of OCTA images was selected for detailed analysis.

2.4. Feature Extraction

Quantitative evaluation of the OCTA images obtained was performed based on automated analysis using OCTAVA (OCTA Vascular Analyzer) software (version 1.0) developed by Untrach et al. [31] and semi-automatic analysis in IMAGEnet® 6 software (version 1.26.16898, Topcon, Tokyo, Japan). The SCP images were preprocessed for noise detuning, shadow elimination, and enhancement of pixel position information using rolling ball background subtraction [32]. The processed images were then segmented and skeletonized in OCTAVA software. Vessel segmentation was carried out using the fuzzy means method and Frangi filtering with the maximum kernel value set to 2. After automatic skeletonization, the software assigns vessels and their connections to different categories. On the basis of this, quantitative indices of vessel area and length density (VAD and VLD, respectively), total number of nodes, branchpoint density (BD), fractal dimension (FD), and length, diameter, and tortuosity distributions are then calculated. Each feature describes a different aspect of morphology of the retina microvascular network. Feature distributions in a single image are presented in the form of descriptive statistics: mean, median, variance, skewness (measure of the asymmetry), and kurtosis (measure of the peakedness or flatness of the distribution). Based on the contours marked manually in IMAGEnet® the foveal avascular zone (FAZ) areas in the SCP and DCP layers were calculated. For a selected subset of MNV, the contour of neovascularization was marked in the CC and OR layers, for which the area was also calculated.

2.5. Statistical Analysis

The determined variables were used as a first step in a comparative analysis between control subjects and AMD patients. For this purpose, the t-test was used if the assumptions of normality of distributions and homogeneity of variance were met (Shapiro–Wilk and F tests, respectively). Otherwise, the Mann–Whitney U test was used. The values for both eyes of each subject were averaged beforehand so that one subject represented a single observation. Univariate testing of the correlation of individual vessel morphological features with PRS was performed for the averaged values using the Pearson correlation coefficient. The general association of features with PRS values in the multivariate approach was verified using the generalized estimating equation (GEE) linear model. This approach allows the features of each eye to be used independently through the ability to model repeated measures. The influence of covariates in the form of subjects’ age at examination, gender, pigment epithelial detachment (PED), MNV, and AMD were included in this analysis. The significance level was established at 0.05.

3. Results

Figure 2 shows the distributions of PRS values that represent the differences between the AMD patient group and the control group resulting from the first part of the ongoing study, based on the total data collection. The statistical analysis revealed a significant difference in PRS between groups ( p = 4.93 × 10 15 ) with a large effect size ( r = 0.58 ).
Table 1 presents the characteristics of the control group and patients with AMD from the selected dataset. There were no significant age differences between AMD patients and the control group. Excluding part of the data due to the established criteria for the OCTA SCP image introduced significant differences in the sex ratio between the groups. As expected, there was a significant deterioration in visual acuity and a reduction in choroidal thickness in the AMD group relative to the control group.
Image selection reduced the number of study participants and, consequently, the number of unique PRS values for individual subjects. Despite this, the significant difference of PRS between groups has been preserved ( p = 1.02 × 10 8 r = 0.43 ) (Figure 3a). In the subgroup of subjects diagnosed with MNV based on OCTA analysis, the PRS value is significantly higher than the group without MNV ( p = 0.013 r = 0.19 ). However, this difference is much smaller and may result from the disproportion in numbers between the AMD and control groups. Moreover, the range of PRS values in the MNV group is relatively wide and includes both subjects with low and high risk of AMD (Figure 3b). The results confirm the feasibility of analyzing the association of OCTA characteristics with PRS in the defined groups.
Figure 4 shows sample OCTA images from the dataset after selection according to established criteria, based on visual assessment. The presented results of OCTAVA software on images after preprocessing demonstrate the possibility of its effective use in the quantitative analysis of morphological features of the retinal microvascular network. The graphical presentation of the classification of individual vessels and structures into subclasses was the basis for the subsequent stages of selection of the analyzed dataset. The values of individual morphological features of the retinal microvascular network were determined based on the results of the classification.
The results of the statistical analysis of the difference in feature values between the group of patients with AMD and the control group are shown in Table 2. The mean FD is significantly larger in the subjects with AMD. The higher the FD, the more complex and developed the structure of the blood vessels is; we find the FD parameter as a general measure of the complexity of the blood vessel network. The standard deviation of FD is significantly higher in AMD. This indicates a greater variability of the complexity of the vessel structure within one image. The standard deviation of the vessels’ diameter is greater in AMD. Furthermore, another characteristic that differs significantly in patients with AMD is the skewness of the diameter distribution of the vessel. In both groups, it is negative, which means that there are more vessels in the image with a much smaller diameter compared to the diameter of the remaining vessels (many outliers with low values). In AMD, the skewness was greater (closer to zero), which means that fewer very thin vessels (much thinner than the rest) appeared in the images of AMD vessels than in the control group. In both cases, the skewness is greater than −1, which means that the distribution is close to normal, so the number of outliers is small.
Table 3 shows the results of the univariate and multivariate investigation of the relationship between individual morphological features of the retinal microvascular network and PRS including additional covariants. Using both correlation analysis and the GEE linear regression model, neither showed a direct and significant association of the characteristics of the vascular network with AMD PRS.
The general observation is that there is no association between the vascular network and PRS. The differences between FD and diameter characteristics in AMD and the control group are secondary to the occurrence of AMD. They do not constitute a biomarker of the risk of AMD. There is no genotype–phenotype relationship and the expression of AMD-related SNPs is not related to the characteristics of the retinal vascular network. It is difficult to clearly say whether these vessels are difficult to detect or whether they are absent. This may be due to retinal architecture disorders secondary to changes at the border of the RPE and outer retina.

4. Discussion

OCTA is a test whose interpretation may be problematic due to numerous artifacts occurring in the results. The rather long scanning time and difficulties in cooperation with the patient may pose difficulties. Problems with cooperation during examinations and eye movement artifacts during long examination time (longer than 10 s per eye) led to the rejection of many scans during our analysis. OCTA is a diagnostic tool that still requires the development of appropriate algorithms that integrate individual B-scans into a whole image [18]. Artifacts were also a problem with the segmentation algorithm, which sometimes marked an obvious line of artifact as a vessel. Rejecting such photos also narrowed the number of studies eligible for further analysis. The imperfection of this vascular imaging method was also highlighted in our study, where we struggled with many artifacts in the raw OCTA data. In three selection stages, we selected high-quality images for which the correct automatic segmentation of vessels was possible. The distribution of PRS values in the selected dataset correctly reflects the distribution in the entire dataset. However, the reduction resulted in a difference in the gender ratio between the group of patients with AMD and the control group. In the group of patients with AMD, the percentage of women was 55.05%, and the rate of men was 44.95%, while in the control group, a larger percentage was women (76.75%) and a smaller percentage was men (23.25%).
During the study, researchers repeatedly reported the need to manually mark the area of neovascularization and FAZ due to the imperfection of automatic algorithms [27]. When working with OCTA, it is important to consider the impact of low signal strength on image quality. Manufacturers recommend specific thresholds for sufficient signal strength, and, even if these measures are met, signal strength can still affect quantitative results [33,34]. For instance, improving signal strength may lead to an increase in retinal vessel densities. Therefore, ensuring adequate signal strength is crucial for accurate OCTA assessments [33,35]. In addition to signal strength, the positioning of scanning frames according to anatomical landmarks plays a significant role in ensuring comparability and reproducibility of quantitative OCTA images. Scanning frames are typically centered on the optic nerve head or the fovea. However, it is important to be aware of centration artifacts that can occur during OCTA analysis, with their frequency varying depending on whether the scanning frames are positioned manually or automatically by the OCTA device [18]. These artifacts are less relevant for smaller areas of interest, such as the foveal avascular zone (3 mm × 3 mm scans), but become critical when evaluating larger areas like the whole perifoveal macula (9 mm × 9 mm or 12 mm × 12 mm scans). The study conducted by Wicklein et al. has emphasized the importance of confirming accurate centration in every OCTA image to account for potential technical impacts on quantitative measures [36]. During our study, a panel of experts verified the correctness of the central positioning of the scanning frame in relation to the avascular zone or central neovascularization.
Furthermore, algorithmic failure in the segmentation of different retinal layers, such as the inner limiting membrane, the inner plexiform layer, or the retinal pigment epithelium, can lead to inaccurate mapping of blood flow signals to the retinal vascular plexus [3]. This issue is especially relevant for OCTA images incorporating retinal pathologies, as disruptions of the retinal architecture can severely affect segmentation and image quality [3,37]. It is essential to conduct a thorough assessment of accurate retinal layer segmentation within different B-scans to ensure correct segmentation, considering the potential impact of retinal pathologies on both retinal layer architecture and vasculature when interpreting OCTA images [38]. In our study, we also noticed similar problems: the proximity of superficial and deep plexuses and choriocapillaris, as well as their different thicknesses in different retinas, resulted in the interpenetration of these images, which led to incorrect automatic segmentation and necessitated the intervention of experts to correct errors in marking individual layers.
Moreover, the detection and differentiation of motion artifacts can be challenging, particularly for individuals with limited OCTA experience. Various scoring systems have been proposed to evaluate the severity of motion artifacts, distinguishing between minor artifacts that do not disrupt the retinal vasculature and major artifacts that can significantly affect the imaging of retinal vessels. Therefore, being mindful of these factors is crucial for obtaining accurate and reliable OCTA results [39].
PRS is of interest to researchers in relation to many diseases, and, due to the polygenic basis of many diseases, it is used as an important indicator in numerous scientific studies in many specializations [36,40,41,42]. PRS in the group presented is described in our previous publications, where the detailed genetic testing process is described [30]. Sekimitsu et al. conducted a rare study examining the correlation between retinal OCT images and PRS, with a focus on evaluating the risk of cognitive impairment [39]. The research employed a logistic regression model to evaluate the link between Alzheimer’s disease and Parkinson’s disease PRS and retinal layer thickness. The findings affirmed the connection between genetic predisposition to neurodegenerative diseases and retinal OCT measurements. Additionally, the study identified potential predictive biomarkers for cognitive decline, suggesting that incorporating retinal measurements enhances the predictive accuracy of Alzheimer’s and Parkinson’s disease.

5. Conclusions

Although the superficial vascular network is not involved in AMD directly, it can be significantly altered by the change of shape of the retina above MNV. OCTA is used most reliably to assess the superficial weave, and, to a slightly lesser extent, a deep weave. It was found that both are independent of the genetic changes related to AMD. However, we expected that the characteristics of the above weaves may change secondary to external retina and RPE disorders, especially in the form of neovascularization. It cannot be ruled out that analysis based on a larger group would show such changes. We take into account that the need to exclude a large group of patients is a limiting factor, because the scale of vascular lesions in this group is probably the largest, but we do not currently have the capability of checking the vascular network measurements in this group.

Author Contributions

Conceptualization, A.S., A.W.M., E.W. and S.T.; methodology, A.S., D.L. and S.T.; software, D.L., A.T. and J.N.; validation, A.T., A.B.-C., A.W. (Adam Wylęgała) and E.W.; formal analysis, D.L., A.T., J.N. and A.W. (Anna Wąsowska); investigation, A.S., A.W. (Anna Wąsowska) and A.W. (Adam Wylęgała); resources, A.B.-C., A.W.M. and S.T.; data curation, A.S., J.N. and A.W. (Anna Wąsowska); writing—original draft preparation, A.S., D.L., A.T., J.N. and A.W. (Anna Wąsowska); writing—review and editing, A.B.-C., A.W. (Adam Wylęgała), A.W.M., E.W. and S.T.; visualization, D.L.; supervision, A.B.-C., A.W.M., E.W. and S.T.; project administration, A.W.M. and S.T.; funding acquisition, E.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Centre for Research and Development grant STRATEGMED1/234261/2/NCBR/2014.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Ethics Committee of the Medical University od Silesia (Resolutions No KNW/0022/KB1/105/13, date of approval: 8 October 2013, and BNW/NWN/0052/KB1/97/I/ 22/23, date of approval: 4 April 2023).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The datasets generated and analyzed during the current study are available from the corresponding author on reasonable request.

Acknowledgments

This work is part of a PhD thesis by Adam Sendecki conducted at the Medical University of Silesia in cooperation with the Faculty of Biomedical Engineering, Silesian University of Technology and Genomed S.A. company.

Conflicts of Interest

Anna Wąsowska and Anna Boguszewska-Chachulska are employees of Genomed S.A. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results. Other authors declare no conflicts of interest.

References

  1. Fleckenstein, M.; Keenan, T.D.; Guymer, R.H.; Chakravarthy, U.; Schmitz-Valckenberg, S.; Klaver, C.C.; Wong, W.T.; Chew, E.Y. Age-related macular degeneration. Nat. Rev. Dis. Prim. 2021, 7, 31. [Google Scholar] [CrossRef] [PubMed]
  2. Pauleikhoff, D.; Harper, C.A.; Marshall, J.; Bird, A.C. Aging changes in Bruch’s membrane: A histochemical and morphologic study. Ophthalmology 1990, 97, 171–178. [Google Scholar] [CrossRef] [PubMed]
  3. Koller, A.; Lamina, C.; Brandl, C.; Zimmermann, M.E.; Stark, K.J.; Weissensteiner, H.; Würzner, R.; Heid, I.M.; Kronenberg, F. Systemic Evidence for Mitochondrial Dysfunction in Age-Related Macular Degeneration as Revealed by mtDNA Copy Number Measurements in Peripheral Blood. Int. J. Mol. Sci. 2023, 24, 16406. [Google Scholar] [CrossRef] [PubMed]
  4. Velilla, S.; García-Medina, J.J.; García-Layana, A.; Dolz-Marco, R.; Pons-Vázquez, S.; Pinazo-Durán, M.D.; Gómez-Ulla, F.; Arévalo, J.F.; Díaz-Llopis, M.; Gallego-Pinazo, R.; et al. Smoking and age-related macular degeneration: Review and update. J. Ophthalmol. 2013, 2013, 895147. [Google Scholar] [CrossRef] [PubMed]
  5. Brandl, C.; Finger, R.P.; Heid, I.M.; Mauschitz, M.M. Age-Related Macular Degeneration in an Ageing Society-Current Epidemiological Research. Klin. Monatsblatter Augenheilkd. 2023, 240, 1052–1059. [Google Scholar]
  6. Lüdtke, L.; Jürgens, C.; Ittermann, T.; Völzke, H.; Tost, F. Age-related macular degeneration and associated risk factors in the population-based study of health in pomerania (SHIP-Trend). Med Sci. Monit. Int. Med. J. Exp. Clin. Res. 2019, 25, 6383. [Google Scholar] [CrossRef] [PubMed]
  7. Bullock, J.; Polato, F.; Abu-Asab, M.; Bernardo-Colón, A.; Aflaki, E.; Agbaga, M.P.; Becerra, S.P. Degradation of photoreceptor outer segments by the retinal pigment epithelium requires pigment epithelium-derived factor receptor (PEDF-R). Investig. Ophthalmol. Vis. Sci. 2021, 62, 30. [Google Scholar] [CrossRef] [PubMed]
  8. Deng, Y.; Qiao, L.; Du, M.; Qu, C.; Wan, L.; Li, J.; Huang, L. Age-related macular degeneration: Epidemiology, genetics, pathophysiology, diagnosis, and targeted therapy. Genes Dis. 2022, 9, 62–79. [Google Scholar] [CrossRef] [PubMed]
  9. Winkler, T.W.; Grassmann, F.; Brandl, C.; Kiel, C.; Günther, F.; Strunz, T.; Weidner, L.; Zimmermann, M.E.; Korb, C.A.; Poplawski, A.; et al. Genome-wide association meta-analysis for early age-related macular degeneration highlights novel loci and insights for advanced disease. BMC Med. Genom. 2020, 13, 120. [Google Scholar] [CrossRef]
  10. Age-Related Eye Disease Study Research Group. Risk factors associated with age-related macular degeneration: A case-control study in the age-related eye disease study: Age-related eye disease study report number 3. Ophthalmology 2000, 107, 2224–2232. [Google Scholar] [CrossRef]
  11. Serra, R.; Rallo, V.; Pinna, A.; Steri, M.; Piras, M.G.; Marongiu, M.; Coscas, F.; Gorospe, M.; Schlessinger, D.; Fiorillo, E.; et al. Polygenic risk score and biochemical/environmental variables predict a low-risk profile of age-related macular degeneration in Sardinia. Graefe’s Arch. Clin. Exp. Ophthalmol. 2023, 261, 691–698. [Google Scholar] [CrossRef]
  12. Yu, C.; Robman, L.; He, W.; Woods, R.L.; Wolfe, R.; Phung, J.; Makeyeva, G.A.; Hodgson, L.A.; McNeil, J.J.; Guymer, R.H.; et al. Predictive performance of an updated polygenic risk score for age-related macular degeneration. Ophthalmology 2024, in press. [Google Scholar] [CrossRef]
  13. Shijo, T.; Sakurada, Y.; Yoneyama, S.; Kikushima, W.; Sugiyama, A.; Matsubara, M.; Fukuda, Y.; Mabuchi, F.; Kashiwagi, K. Association between polygenic risk score and one-year outcomes following as-needed aflibercept therapy for exudative age-related macular degeneration. Pharmaceuticals 2020, 13, 257. [Google Scholar] [CrossRef] [PubMed]
  14. Ulanczyk, Z.; Grabowicz, A.; Mozolewska-Piotrowska, K.; Safranow, K.; Kawa, M.P.; Pałucha, A.; Krawczyk, M.; Sikora, P.; Matczynska, E.; Machaliński, B.; et al. Genetic factors associated with age-related macular degeneration: Identification of a novel PRPH2 single nucleotide polymorphism associated with increased risk of the disease. Acta Ophthalmol. 2021, 99, 739–749. [Google Scholar] [CrossRef]
  15. Zekavat, S.M.; Sekimitsu, S.; Ye, Y.; Raghu, V.; Zhao, H.; Elze, T.; Segrè, A.V.; Wiggs, J.L.; Natarajan, P.; Del Priore, L.; et al. Photoreceptor layer thinning is an early biomarker for age-related macular degeneration: Epidemiologic and genetic evidence from UK Biobank OCT data. Ophthalmology 2022, 129, 694–707. [Google Scholar] [CrossRef] [PubMed]
  16. Tang, W.; Ma, J.; Gu, R.; Lei, B.; Ding, X.; Xu, G. Light-induced lipocalin 2 facilitates cellular apoptosis by positively regulating reactive oxygen species/Bim signaling in retinal degeneration. Investig. Ophthalmol. Vis. Sci. 2018, 59, 6014–6025. [Google Scholar] [CrossRef]
  17. Faatz, H.; Rothaus, K.; Ziegler, M.; Book, M.; Heimes-Bussmann, B.; Pauleikhoff, D.; Lommatzsch, A. Vascular analysis of type 1, 2, and 3 macular neovascularization in age-related macular degeneration using swept-source optical coherence tomography angiography shows new insights into differences of pathologic vasculature and may lead to a more personalized understanding. Biomedicines 2022, 10, 694. [Google Scholar] [CrossRef]
  18. Campochiaro, P.A. Ocular neovascularization. J. Mol. Med. 2013, 91, 311–321. [Google Scholar] [CrossRef]
  19. Prenner, J.L.; Halperin, L.S.; Rycroft, C.; Hogue, S.; Liu, Z.W.; Seibert, R. Disease burden in the treatment of age-related macular degeneration: Findings from a time-and-motion study. Am. J. Ophthalmol. 2015, 160, 725–731. [Google Scholar] [CrossRef]
  20. Steinmetz, J.D.; Bourne, R.R.; Briant, P.S.; Flaxman, S.R.; Taylor, H.R.; Jonas, J.B.; Abdoli, A.A.; Abrha, W.A.; Abualhasan, A.; Abu-Gharbieh, E.G.; et al. Causes of blindness and vision impairment in 2020 and trends over 30 years, and prevalence of avoidable blindness in relation to VISION 2020: The Right to Sight: An analysis for the Global Burden of Disease Study. Lancet Glob. Health 2021, 9, e144–e160. [Google Scholar] [CrossRef]
  21. Leroux, P.; Agard, E.; Billant, J.; Levron, A.; Bouvarel, H.; Badri, Y.; Douma, I.; Pradat, P.; Dot, C. Long Intervals between intravitreal injections using a Treat-and-Extend protocol in a Real-life context in AMD: The LIRE study. Ophthalmologica 2023, 247, 44–57. [Google Scholar] [CrossRef] [PubMed]
  22. Pugazhendhi, A.; Hubbell, M.; Jairam, P.; Ambati, B. Neovascular macular degeneration: A review of etiology, risk factors, and recent advances in research and therapy. Int. J. Mol. Sci. 2021, 22, 1170. [Google Scholar] [CrossRef] [PubMed]
  23. Colijn, J.M.; Buitendijk, G.H.; Prokofyeva, E.; Alves, D.; Cachulo, M.L.; Khawaja, A.P.; Cougnard-Gregoire, A.; Merle, B.M.; Korb, C.; Erke, M.G.; et al. Prevalence of age-related macular degeneration in Europe: The past and the future. Ophthalmology 2017, 124, 1753–1763. [Google Scholar] [CrossRef] [PubMed]
  24. Tun, Y.Z.; Aimmanee, P. A Complete Review of Automatic Detection, Segmentation, and Quantification of Neovascularization in Optical Coherence Tomography Angiography Images. Diagnostics 2023, 13, 3407. [Google Scholar] [CrossRef] [PubMed]
  25. Spaide, R.F.; Fujimoto, J.G.; Waheed, N.K. Image artifacts in optical coherence tomography angiography. Retina 2015, 35, 2163. [Google Scholar] [CrossRef] [PubMed]
  26. Zhu, Y.; Cui, Y.; Wang, J.C.; Lu, Y.; Zeng, R.; Katz, R.; Wu, D.M.; Eliott, D.; Vavvas, D.G.; Husain, D.; et al. Different scan protocols affect the detection rates of diabetic retinopathy lesions by wide-field swept-source optical coherence tomography angiography. Am. J. Ophthalmol. 2020, 215, 72–80. [Google Scholar] [CrossRef] [PubMed]
  27. Dabir, S.; Bhatt, V.; Bhatt, D.; Rajan, M.; Samant, P.; Munusamy, S.; Webers, C.; Berendschot, T. Need for manual segmentation in optical coherence tomography angiography of neovascular age-related macular degeneration. PLoS ONE 2020, 15, e0244828. [Google Scholar] [CrossRef] [PubMed]
  28. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011, 17, 10–12. [Google Scholar] [CrossRef]
  29. Van der Auwera, G.A.; Carneiro, M.O.; Hartl, C.; Poplin, R.; Del Angel, G.; Levy-Moonshine, A.; Jordan, T.; Shakir, K.; Roazen, D.; Thibault, J.; et al. From FastQ data to high-confidence variant calls: The genome analysis toolkit best practices pipeline. Curr. Protoc. Bioinform. 2013, 43, 11.10.1–11.10.33. [Google Scholar] [CrossRef]
  30. Wąsowska, A.; Sendecki, A.; Boguszewska-Chachulska, A.; Teper, S. Polygenic Risk Score and Rare Variant Burden Identified by Targeted Sequencing in a Group of Patients with Pigment Epithelial Detachment in Age-Related Macular Degeneration. Genes 2023, 14, 1707. [Google Scholar] [CrossRef]
  31. Untracht, G.R.; Matos, R.S.; Dikaios, N.; Bapir, M.; Durrani, A.K.; Butsabong, T.; Campagnolo, P.; Sampson, D.D.; Heiss, C.; Sampson, D.M. OCTAVA: An open-source toolbox for quantitative analysis of optical coherence tomography angiography images. PLoS ONE 2021, 16, e0261052. [Google Scholar] [CrossRef] [PubMed]
  32. Freedman, I.G.; Li, E.; Hui, L.; Adelman, R.A.; Nwanyanwu, K.; Wang, J.C. The Impact of Image Processing Algorithms on Optical Coherence Tomography Angiography Metrics and Study Conclusions in Diabetic Retinopathy. Transl. Vis. Sci. Technol. 2022, 11, 7. [Google Scholar] [CrossRef] [PubMed]
  33. Lim, H.B.; Kim, Y.W.; Kim, J.M.; Jo, Y.J.; Kim, J.Y. The importance of signal strength in quantitative assessment of retinal vessel density using optical coherence tomography angiography. Sci. Rep. 2018, 8, 12897. [Google Scholar] [CrossRef] [PubMed]
  34. Lee, J.J.; Lee, J.E.; Sadda, S.R.; Park, S.W.; Byon, I. Impact of signal strength on quantitative retinal and choriocapillaris flow measurement from optical coherence tomography angiography. Sci. Rep. 2022, 12, 4692. [Google Scholar] [CrossRef] [PubMed]
  35. Lim, H.B.; Kim, Y.W.; Nam, K.Y.; Ryu, C.K.; Jo, Y.J.; Kim, J.Y. Signal strength as an important factor in the analysis of peripapillary microvascular density using optical coherence tomography angiography. Sci. Rep. 2019, 9, 16299. [Google Scholar] [CrossRef] [PubMed]
  36. Oh, J.J.; Hong, S.K. Polygenic risk score in prostate cancer. Curr. Opin. Urol. 2022, 32, 466–471. [Google Scholar] [CrossRef] [PubMed]
  37. Bontzos, G.; Kabanarou, S.A.; Garnavou-Xirou, C.; Gkizis, I.; Kontou, E.; Triantafyllou, D.; Xirou, T. Segmentation errors and motion artifacts in OCT-A associated with epiretinal membranes. Can. J. Ophthalmol. 2020, 55, 293–300. [Google Scholar] [CrossRef] [PubMed]
  38. Wicklein, R.; Yam, C.; Noll, C.; Aly, L.; Banze, N.; Romahn, E.F.; Wolf, E.; Hemmer, B.; Oertel, F.C.; Zimmermann, H.; et al. The OSCAR-MP consensus criteria for quality assessment of retinal optical coherence tomography angiography. Neurol. Neuroimmunol. Neuroinflamm. 2023, 10, e200169. [Google Scholar] [CrossRef]
  39. Sekimitsu, S.; Shweikh, Y.; Shareef, S.; Zhao, Y.; Elze, T.; Segrè, A.; Wiggs, J.; Zebardast, N. Association of retinal optical coherence tomography metrics and polygenic risk scores with cognitive function and future cognitive decline. Br. J. Ophthalmol. 2023, 108, 599–606. [Google Scholar] [CrossRef]
  40. Saks, D.G.; Schulz, A.; Qassim, A.; Marshall, H.; Hewitt, A.W.; MacGregor, S.; Craig, J.E.; Graham, S.L. Genetic risk of glaucoma is associated with vascular and retinal nerve fibre wedge defects. Acta Ophthalmol. 2023, 102, e185–e194. [Google Scholar] [CrossRef]
  41. Wang, X.Y.; Wang, L.L.; Xu, L.; Liang, S.Z.; Yu, M.C.; Zhang, Q.Y.; Dong, Q.J. Evaluation of polygenic risk score for risk prediction of gastric cancer. World J. Gastrointest. Oncol. 2023, 15, 276. [Google Scholar] [CrossRef] [PubMed]
  42. Briggs, S.E.; Law, P.; East, J.E.; Wordsworth, S.; Dunlop, M.; Houlston, R.; Hippisley-Cox, J.; Tomlinson, I. Integrating genome-wide polygenic risk scores and non-genetic risk to predict colorectal cancer diagnosis using UK Biobank data: Population based cohort study. BMJ 2022, 379, e071707. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Flow chart of OCTA image selection into subgroups considered in the study.
Figure 1. Flow chart of OCTA image selection into subgroups considered in the study.
Diagnostics 14 00770 g001
Figure 2. Distribution of PRS values per individual as density plot.
Figure 2. Distribution of PRS values per individual as density plot.
Diagnostics 14 00770 g002
Figure 3. Distribution of PRS values for AMD patients and control group (a) and for MNV and non-MNV subjects (b) from the subset of participants with proper OCTA images selected.
Figure 3. Distribution of PRS values for AMD patients and control group (a) and for MNV and non-MNV subjects (b) from the subset of participants with proper OCTA images selected.
Diagnostics 14 00770 g003
Figure 4. Example of en-face OCTA images (top row) and corresponding classification results of different architectural components of the vessel network (bottom row) for subjects from the control group and with AMD, for different AMD PRS values.
Figure 4. Example of en-face OCTA images (top row) and corresponding classification results of different architectural components of the vessel network (bottom row) for subjects from the control group and with AMD, for different AMD PRS values.
Diagnostics 14 00770 g004
Table 1. Clinical characteristics of patients with AMD and control group—mean (standard deviation).
Table 1. Clinical characteristics of patients with AMD and control group—mean (standard deviation).
AMDControlp-Value
N8943
Age [years]72.47 (7.99)71.00 (7.04)0.307
Sex [male/female]40/4910/330.026
Visual acuity [logMAR]0.48 (0.45)0.11 (0.18)<0.001
Choroidal thickness [ μ m]237.23 (105.45)273.49 (82.03)0.001
Table 2. Comparison of OCTA features values between group of patients with AMD and control group.
Table 2. Comparison of OCTA features values between group of patients with AMD and control group.
AMDControl
Feature M SD M SD p-Value
Diameter [ μ m]Mean47.977.1046.671.470.4858
Median48.646.8747.421.640.9512
Variance35.3176.6925.214.230.0193
Skewness−0.210.64−0.380.440.0457
Kurtosis2.691.852.510.940.1231
Length [mm]Mean774.4479.50770.4446.240.5359
Median578.4176.37603.6263.090.0920
Total50.184.8750.463.990.9729
Variance515,754.41195,925.48449,083.87106,514.090.0521
Skewness1.970.551.850.480.2512
Kurtosis8.043.267.622.900.5247
TortuosityMean0.11390.01210.11550.01060.3614
Variance0.01900.00940.01810.00730.8270
Skewness3.040.723.120.800.8631
Kurtosis14.2195.67115.2677.0020.6392
Fractal DimensionMean1.420.031.400.010.0024
Variance0.5250.0240.5160.0050.0123
Nodes 43.769.1942.646.130.4683
VLD [%] 2.6760.2332.6280.2080.2574
VAD [%] 15.8534.14315.0811.4060.1373
BD [Nodes/mm] 0.8630.1290.8420.0790.3289
FAZSCP [mm2] 0.47370.38280.34460.00150.0616
FAZDCP [mm2] 0.51460.42610.40340.15550.6727
VLD—vessel length density; VAD—vessel area density; BD—branchpoint density; FAZ—foveal avascular zone; SCP—superficial capillary plexus; DCP—deep capillary plexus. Bold indicates statistically significant differences.
Table 3. Correlation and multivariate regression analysis result between the subjects’ characteristics and retinal morphological features with AMD PRS.
Table 3. Correlation and multivariate regression analysis result between the subjects’ characteristics and retinal morphological features with AMD PRS.
UnivariateMultivariate (GEE)
Feature Pearson’s r p -Value Coefficient p -Value
CovariantsAMD0.8041<0.001
Age−0.00130.915
Sex0.39240.014
PED0.31960.101
MNV0.02550.893
Diameter [ μ m]Mean0.0420.6301.3440.179
Median0.0400.650−1.2340.217
Variance0.0160.853−0.3840.701
Skewness−0.0490.574−0.0550.956
Kurtosis−0.1000.253−0.2980.765
Length [mm]Mean−0.0120.893−0.7090.479
Median−0.0780.3720.3450.730
Total−0.1240.157−0.3080.758
Variance0.0700.422−0.4310.667
Skewness0.0930.2881.8310.067
Kurtosis0.0560.525−1.5170.129
TortuosityMean0.0750.3961.3150.188
Variance0.0430.621−0.0080.994
Skewness−0.1420.105−1.5420.123
Kurtosis−0.1440.1011.1570.247
Fractal DimensionMean0.0270.755−0.5080.611
Variance0.0720.414−0.0780.938
Nodes −0.0270.7560.4300.667
VLD [%] −0.0980.2630.1240.901
VAD [%] −0.0110.8990.0590.953
BD [Nodes/mm] 0.0240.788−0.2960.768
FAZSCP [mm2] 0.1050.230−0.1330.894
FAZDCP [mm2] 0.1030.2410.4100.682
MNVCC area [mm2] −0.0540.692−0.7710.441
MNVOR area [mm2] 0.1040.6921.3980.162
VLD—vessel length density; VAD—vessel area density; BD—branchpoint density; FAZ—foveal avascular zone; SCP—superficial capillary plexus; DCP—deep capillary plexus; MNV—macular neovascularization; CC—choriocapillaris; OR—outer retina. Bold indicates statistically significant coefficients.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sendecki, A.; Ledwoń, D.; Tuszy, A.; Nycz, J.; Wąsowska, A.; Boguszewska-Chachulska, A.; Wylęgała, A.; Mitas, A.W.; Wylęgała, E.; Teper, S. Association of Genetic Risk for Age-Related Macular Degeneration with Morphological Features of the Retinal Microvascular Network. Diagnostics 2024, 14, 770. https://doi.org/10.3390/diagnostics14070770

AMA Style

Sendecki A, Ledwoń D, Tuszy A, Nycz J, Wąsowska A, Boguszewska-Chachulska A, Wylęgała A, Mitas AW, Wylęgała E, Teper S. Association of Genetic Risk for Age-Related Macular Degeneration with Morphological Features of the Retinal Microvascular Network. Diagnostics. 2024; 14(7):770. https://doi.org/10.3390/diagnostics14070770

Chicago/Turabian Style

Sendecki, Adam, Daniel Ledwoń, Aleksandra Tuszy, Julia Nycz, Anna Wąsowska, Anna Boguszewska-Chachulska, Adam Wylęgała, Andrzej W. Mitas, Edward Wylęgała, and Sławomir Teper. 2024. "Association of Genetic Risk for Age-Related Macular Degeneration with Morphological Features of the Retinal Microvascular Network" Diagnostics 14, no. 7: 770. https://doi.org/10.3390/diagnostics14070770

APA Style

Sendecki, A., Ledwoń, D., Tuszy, A., Nycz, J., Wąsowska, A., Boguszewska-Chachulska, A., Wylęgała, A., Mitas, A. W., Wylęgała, E., & Teper, S. (2024). Association of Genetic Risk for Age-Related Macular Degeneration with Morphological Features of the Retinal Microvascular Network. Diagnostics, 14(7), 770. https://doi.org/10.3390/diagnostics14070770

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop