Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

Multi-task learning for activity detection in neovascular age-related macular degeneration

View ORCID ProfileMurat Seçkin Ayhan, Hanna Faber, Laura Kühlewein, Werner Inhoffen, Gulnar Aliyeva, Focke Ziemssen, Philipp Berens
doi: https://doi.org/10.1101/2022.06.13.22276315
Murat Seçkin Ayhan
1Institute for Ophthalmic Research, University of Tübingen, Tübingen, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Murat Seçkin Ayhan
Hanna Faber
1Institute for Ophthalmic Research, University of Tübingen, Tübingen, Germany
2University Eye Clinic, University of Tübingen, Tübingen, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Laura Kühlewein
1Institute for Ophthalmic Research, University of Tübingen, Tübingen, Germany
2University Eye Clinic, University of Tübingen, Tübingen, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Werner Inhoffen
2University Eye Clinic, University of Tübingen, Tübingen, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gulnar Aliyeva
2University Eye Clinic, University of Tübingen, Tübingen, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Focke Ziemssen
2University Eye Clinic, University of Tübingen, Tübingen, Germany
3University Eye Clinic, University of Leipzig, Leipzig, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Philipp Berens
1Institute for Ophthalmic Research, University of Tübingen, Tübingen, Germany
4Tübingen AI Center, Tübingen, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: philipp.berens{at}uni-tuebingen.de
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Data/Code
  • Preview PDF
Loading

Abstract

Purpose Comparison of performance and explainability of a multi-task convolutional deep neuronal network to single-task networks for activity detection in neovascular age-dependent macular degeneration.

Methods From n = 70 patients (46 female, 24 male) who attended the University Eye Hospital Tübingen 3762 optical coherence tomography B-scans (right eye: 2011, left eye: 1751) were acquired with Heidelberg Spectralis, Heidelberg, Germany. B-scans were graded by a retina specialist and an ophthalmology resident, and then used to develop a multi-task deep learning model to predict disease activity in neovascular age-related macular degeneration along with the presence of sub- and intraretinal fluid. We used performance metrics for comparison to single-task networks and visualized the DNN-based decision with t-distributed stochastic neighbor embedding and clinically validated saliency mapping techniques.

Results The multi-task model surpassed single-task networks in accuracy for activity detection (94.2). Further-more, compared to single-task networks, visualizations via t-distributed stochastic neighbor embedding and saliency maps highlighted that multi-task networks’ decisions for activity detection in neovascular age-related macular degeneration were highly consistent with the presence of both sub- and intraretinal fluid.

Conclusions Multi-task learning increases the performance of neuronal networks for predicting disease activity, while providing clinicians with an easily accessible decision control, which resembles human reasoning.

Translational Relevance By improving nAMD activity detection performance and transparency of automated decisions, multi-task DNNs can support the translation of machine learning research into clinical decision support systems for nAMD activity detection.

1 Introduction

Neovascular age-related macular degeneration (nAMD) is a sight-threatening disease and a common cause of vision loss worldwide.1–3 Among the basic features of nAMD are subretinal fluid (SRF) and intraretinal fluid (IRF), which serve as surrogate markers of nAMD activity and can be monitored using optical coherence tomography (OCT)4,5 (Fig. 1).

Figure 1:
  • Download figure
  • Open in new tab
Figure 1:

Exemplary retinal images (B-scans) with neovascular age-related macular degeneration (nAMD). Solid and dotted arrows indicate subretinal and intraretinal fluid, respectively. (a): no nAMD activity. (b): nAMD activity due to subretinal fluid (SRF). (c): nAMD activity due to intraretinal fluid (IRF). (d): nAMD activity due to both SRF and IRF.

In nAMD, increased levels of vascular endothelial growth factor (VEGF) lead to formation of new vessels from the choroidal and/or retinal vasculature. If leakage from these vessels exceeds local clearance rates, fluid builds up, leading to IRF and SRF.4 IRF is assumed to originate from vascular leakage from intraretinal neovasculaturisation and/or retinal vasculature or from diffusion through the outer retina due to changes within the external limiting membrane.4 In contrast, SRF formation likely results from malfunction of the retinal pigment epithelium with reduced removal rates.4 Due to the partially different pathophysiology, IRF and SRF can occur both simultaneously and independently from each other.4,6 In addition, the characterisation of the lesion based on IRF and SRF could help to determine visual outcome.7

Treatment with intravitreal anti-VEGF agents efficiently restores the balance between fluid formation and retinal removal and is standard of care, when IRF or SRF in nAMD is detected via OCT.5 Prompt treatment initiation is necessary to prevent vision loss.8–10 Additionally, this chronic disease demands high-frequent therapy monitoring, which has put considerable burden on patients, their families and ophthalmological care since its initial approval in 2006.11–14 Since the number of patients suffering from AMD is thought to rise from 196 million in 2020 to 288 million in 2040, the care needed will also rise.2 Hence, automated solutions making the diagnostic processes more efficient have considerable appeal. For example, deep neural networks (DNNs) have been used for automatic referral decisions15 and predicting disease conversion to nAMD.16 Automated algorithms could detect both SRF and IRF more reliably than retinal specialists especially in less conspicuous cases.17 Ideally, such automated tools serve to support retinal specialists in their decision making. In collaboration, a retina specialist assisted by an artificial intelligence (AI) tool can outperform the model alone, e.g., for the task of diabetic retinopathy grading.18 To this end, computational tools need to explain their decisions and communicate their uncertainty to the treating ophthalmologist.19, 20

Here, we develop a convolutional deep learning model based on the concept of multi-task learning.21, 22 Multi-task learning is a generalization of the widely used single-task learning, where models are trained for multiple input-output mappings simultaneously (Fig. 2). For instance, multi-task models can be used to capture different characteristics of dry AMD, such as drusen area, geographic atrophy, increased pig-ment, and depigmentation, to combine these outputs into final AMD diagnosis w.r.t. a 9-step severity scale.23 Multi-task learning has also shown prognostic value when applied to survival analysis via two simultaneous prediction tasks: drusen and pigmentation grading.24 In a similar vein, our multi-task model detects SRF, IRF and nAMD activity in parallel. However, it generates distinct outputs for each of these tasks and offers well-calibrated uncertainty estimates for each of them, which is unique to our study. As the fluid compartment plays a decisive role in the treatment outcome25–27 with the simultaneous presence of IRF and SRF being associated with the worst prognosis,9 we visualize the representation driving the DNN-based decisions using t-distributed stochastic neighbor embedding (t-SNE)28,29 and investigate the model’s decisions using clinically validated saliency mapping techniques.30 Thus, together with well-calibrated uncertainty reports, our work provides an interpretable tool for the ophthalmologist to rapidly access the neural network’s decision process on both population-based and individual-patient levels as a prerequisite for clinical application.

Figure 2:
  • Download figure
  • Open in new tab
Figure 2:

A deep neural network for simultaneous detection of subretinal and intraretinal fluid as well as the nAMD activity from OCT B-scans. Given a B-scan, convolutional stack of the InceptionV3 architecture extracts 2048 feature maps. These are average and max pooled, and fed into a fully connected (dense) layer with 1024 units for shared representation. Then, task-specific heads specialize into individual tasks and single units with sigmoid function achieve binary classification based on 256 task-specific features.

2 Methods

2.1 Data Collection

This study included 70 patients (46 females, 24 males) with nAMD at least in one eye, seen by an ophthalmologist resident (AG) in the Macula clinic at the University Eye Hospital Tübingen. Exclusion criteria were any other cause of neovascularisation, any coexisting retinal pathology (e.g. epiretinal membrane, macular hole, diabetic retinopathy), glaucoma and media opacity preventing sufficient image quality.

3762 B-scans (2011 right eye, 1751 left eye) of 440 × 512 pixels with Heidelberg Spectralis OCT (Heidelberg Engineering, Heidelberg, Germany) were included in the study. A retina specialist of the same hospital (IW) assessed disease activitiy and presence of IRF and SRF on each individual B-Scan. (Fig. 1). Disease activity was also graded by a resident (AG). B-scans were assigned to a training, validation or test set (Table 1). All images of one patient were assigned to one set to avoid information leakage. The study was conducted in accordance with the tenets of the Declaration of Helsinki and approved by the local institutional ethics committee of the University of Tübingen, which waived the requirement for patient consent due to the study’s retrospective character.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 1:

OCT Data distribution of subretinal fluid (SRF), intraretinal fluid (IRF) and active nAMD in B-Scans in training, validation and test sets, respectively. Absolute and relative numbers are shown.

2.2 Diagnostic Tasks, Network Architecture and Model Development

We developed a multi-task DNN to detect the presence of SRF and IRF as well as the nAMD activity from OCT B-scans (Fig. 2). As backbone, we used the InceptionV3 architecture31 via Keras,32 which was pre-trained on ImageNet33 for 1000-way classification via a softmax function. We used the InceptionV3 DNN’s convolutional stack as is but linked max pooling and average pooling layers to the end of convolutional stack and concatenated their outputs to obtain 4096-dimensional feature vectors. These were followed by a dense layer, which yielded a shared representation with 1024 features. To this, we added task-specific heads with 256 units, which specialized into their respective tasks. Then, task-specific binary decisions were achieved by single units equipped with sigmoid functions. For training our DNNs in both single and multi-task scenarios, we resorted to the retina specialist’s set of labels.

We trained our networks with equally weighted cross-entropy losses for all tasks on the training images: Embedded Image, where yn was a vector of binary labels indicating nAMD activity and the presence of IRF or SRF in an image xn. Parameterized by θ, a DNN fθ(·) was optimized with respect to the total cross-entropy on the training data: Embedded Image Embedded Image pn,t was a probability estimated via the sigmoid function for a task indicated by t, and T was the total number of tasks. For T = 1, multi-task learning was reduced to single-task learning based on the same architecture but with only one task head. We also developed a 2-task model to perform the SRF and IRF detection tasks (T = 2), while eliminating the redundancy of the nAMD activity detection task, which is, in principle, a function of the former two.

To address the class imbalance (Table 1), we used random oversampling (see Section 2.2.2 for details). We trained the DNN using Stochastic Gradient Descent (SGD) with Nesterov’s Accelerated Gradients (NAG),34,35 minibatch size of eight, a momentum coefficient of 0.9, an initial learning rate of 5·10−4, a decay rate of 10−6 and a regularization constant of 10−5 for 120 or 150 epochs (see Section 2.2.1 for longer training). During the first five epochs, the convolutional stack was frozen and only dense layers were trained. Then, all layers were fine-tuned to all tasks. The best models were selected based on total validation loss after each epoch and used for inference on the test set.

2.2.1 Data augmentation and preprocessing

We used mixup36 for data augmentation during training. Mixup generates artificial examples through the convex combinations of randomly sampled data points. We adapted mixup to our multi-task learning scenario as follows: Embedded Image Mixing was controlled by λ ∼Beta(α, α), where α ∈ (0, ∞). For α = 0, λ is either 0 or 1, and there is no mixing. We used 0, 0.05 0.1, and 0.2 for α and trained networks for 120 epochs when not mixing and 150 epochs when mixing. Also, to allow for a warm-up period when mixing,36 we set α = 0 for the first five epochs.

In addition, we applied common data augmentation operations such as adjustment of brightness within ±10%, horizontal and vertical flipping, up and down scaling within ±10%, translation of pixels horizontally and vertically within ±30 positions and random rotation within ±45 degrees. After all data augmentation operations, we used an appropriate preprocessing function1 from the Keras API.32

2.2.2 Quantification of uncertainty via mixup and Deep Ensembles

DNNs often do not generate well-calibrated and reliable uncertainty estimates for their decision.37–41 How-ever, quantification of diagnostic uncertainty is crucial for treatment decisions since proper management can minimize diagnostic errors, delays or excess healthcare utilization.42 mixup36 improves the calibration of DNN outputs by smoothing labels through their convex combinations (Eq. 3).43 In addition, we used Deep Ensembles39 consisting of multiple DNNs with different random initializations.39,44 This can improve upon the single network performance both in accuracy and calibration, even with small numbers of DNNs.39,44–46 We used ensembles with three DNNs, for which we enforced diversity by a specialized oversampling strategy: for each DNN, we oversampled training images with respect to one of the task’s labels. This enabled DNNs to train on a balanced dataset while also learning about other tasks, even though the data were not balanced for these. We then used the ensemble’s mean output for predictions and quantified uncertainty in terms of entropy, given the average predictive probabilities.

2.3 Low-dimensional embedding of images

We used t-SNE28 to obtain further insights into the decision-making process of our ensemble model. t-SNE is a non-linear dimensionality reduction method, that embeds high-dimensional data points into a low-dimensional space. We concatenated features from ensemble members’ predetermined read-out layers and performed t-SNE based on them, embedding each B-scan into the two-dimensional plane. We used openTSNE47 with PCA initialization to better preserve the global structure of the data and improve the reproducibility.29 A perplexity of 200 for 1500 iterations with an early exaggeration coefficient of 12 for the first 500 iterations was used according to best-practice strategies.29 Similarities between data points were measured by Euclidean distance in the feature space.

2.4 Saliency Maps

We used Layer-wise Relevance Propagation (LRP)48 to compute saliency maps highlighting the regions in the OCT images which contributed to the DNN decisions, as it provides most clinically relevant tasks.30 We created three saliency maps for each OCT slice: subretinal (cyan), intraretinal (magenta) and diesease activity in nAMD (yellow) (Fig. 5). To improve the visualization of the salient regions, saliency maps were postprocessed.30 Saliency maps were only shown for predictions with an estimated probability greater than 0.5 since previous work has shown, that especially in absence of disease, saliency maps can lead physicans to overdiagnosis.18

3 Results

We developed an ensemble of three multi-task DNNs to simultaneously detect SRF, IRF and activity of nAMD on OCT B-scans (Fig. 1). Each DNN consisted of a shared convolutional core combined with pooling operations and a fully connected (dense) layer (Fig. 2). The resulting shared representation served as the basis for the decisions of the three task-specific heads. The idea behind this approach is that the DNN can benefit from the shared representation induced by combining information from different tasks. We compared the performance of the multi-task model with more specialized single-task models, where we constructed three DNNs for each task, which did not share any representation but were trained independently. In addition, we also used a 2-task model that simultaneously detected only SRF and IRF, without the nAMD activity detection head.

All DNNs were trained on the same dataset (see Table 1 and Methods), which was graded according to the nAMD activity by a retina specialist (IW) and and an ophthalmologist resident (AG) with high inter-grader agreement on disease activity (Cohen’s kappa = 0.86). In a second step, the retina specialist further examined the data for the presence of IRF and SRF. The two retinal fluid types occurred largely independently, while there was natural overlap of both with the active AMD label (Table 2).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 2:

Agreement of task-specific labels across training, validation and test sets, measured via Cohen’s kappa statistic, which is essentially a number between -1 and 1. While 1 indicates a full agreement, lower scores mean less agreement. Negative scores indicate disagreement.

We selected the 3-task model with the best accuracy for the activity detection task on the validation set and report accuracy values computed on an independent test set (Table 3). The 3-task model was well calibrated on the test set (Adaptive expected calibration error41 of 0.0147 for SRF, 0.0104 for IRF and 0.0263 for active nAMD).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 3:

Accuracy of ensembles for various degrees of mixing (indicated by α). Gray row indicates the ensemble of choice for further analysis based on the validation performance for the activity detection task. In the 2-task scenario, the average validation accuracy of SRF and IRF detection tasks was used for model selection.

We found that the performance of the 3-task model surpassed the single-task model performance in disease activity detection, reaching an accuracy of 94.2 % for the multi-task model vs. 91.4% for the single task model (Table 3, Fig. 3). This 3-task model optimized for AMD activity detection performed slightly worse than the single-task models for SRF and IRF detection (SRF: accuracy of 0.917 vs. 0.924 for multi-task vs. single-task; IRF: 0.937 vs. 0.950). For the 2-task scenario, we selected the model with the highest average validation accuracy across the SRF and IRF detection tasks. Interestingly, the 2-task model performed worse than the single-task and 3-task models. This highlights the importance of the explicit nAMD activity detection head in the 3-task model.

Figure 3:
  • Download figure
  • Open in new tab
Figure 3:

Performance curves of the selected models on the test images. Area under the curve (AUC) values given for models also summarize the overall performance into one number (higher is better). (a) Receiver Operating Characteristics (ROC) curves. (b) Precision-recall curves.

We then further studied the representations learned by the models to gain insight into their decision making-process. To this end, we extracted the representations of individual OCT scans from both single-task and multi-task models and created two-dimensional embeddings of these via t-SNE (Fig. 4. In these visualizations, each point represents an individual OCT scan. Scans which are similar to each other according to the learned representation are mapped to nearby points. Of note, distances and in particular the size of white space between clusters in t-SNE plots should be carefully interpreted.29,49

Figure 4:
  • Download figure
  • Open in new tab
Figure 4:

Visualization of data via t-SNE of ensemble-based representations. Only the test data are shown. (a) Low dimensional embedding of images based on the 1024-dimensional features from the prepenultimate layers of single-task networks. Colored with respect to the task-specific labels. (b) Same as in (a) but w.r.t. 1024 features from the shared representation layer of multi-task networks. (c) Same map as in (b) but colored w.r.t. correct and wrong predictions. (d) Same map as in (b) but colored w.r.t. uncertainty min-max normalized to [0, 1].

We labeled individual points according the evidence for SRF or IRF and overall AMD activity. In the single-taks DNNs, well-separated clusters were found, indicating only the learned task-label (Fig. 4a). For example, OCT scans with SRF present formed a single cluster, clearly distinct from the OCT scans without this label. In contrast, in the multi-task network subclusters within the active nAMD data points were observed (Fig. 4a, b): OCT scans labeled with SRF formed a well-separated cluster at the bottom right, as did scans with IRF labels at the top right (Fig. 4b). Interestingly, there was a small cluster in between these two which contained scans labeled with both. This suggests that multi-task DNNs learned a representation which could differentiate between the two fluid types. The few incorrectly classified OCT scans could be found within their clusters to be placed close towards other clusters (Fig. 4c) in areas where we also found examples with high classifier uncertainty (Fig. 4d).

We next studied how the multi-task representations emerged through processing in the network (Appendix, Fig. 7). While in the initial layers data points representing active nAMD were still uniformly distributed (Fig. 7, a-c), a clear separation of active nAMD cases developed gradually in later layers of the DNN (Fig. 7, d-g), leading to best separation in the shared representation (Fig. 7, h). The decision head for active AMD refined this representation only very little (Fig. 7, i).

We finally analyzed the saliency maps of the multi-task DNNs and asked whether the saliency maps for the subtasks of SRF and IRF detection obtained from the multi-task model allowed reasoning about evidence specific to these tasks. We generated saliency maps on four exemplary OCT scans using LRP48 (Figure 5). For an OCT scan with clearly active AMD and both SRF and IRF present (Figure 5a), we found that the active AMD saliency map focused on intraretinal fluids, which were also clearly visible in the task-specific saliency map, and faintly highlighted regions with SRF. The SRF saliency map, however, clearly highlighted SRF. In two further example scans with either IRF or SRF, respectively, active AMD saliency maps clearly corresponded to the individual task maps (Figure 5b,c). We also identified a rare failure case of the obtained saliency maps (Fig. 5d), where an OCT scan was falsely classified positive for SRF with a confidence of 0.614 due to the misclassification of IRF to SRF. We hypothesize that the DNN misclassified the superior border of the IRF as photoreceptor layer detached from the retinal pigment epithelium. The assumption that the DNN primarily recognizes contrast-rich interfaces such as SRF and IRF is further supported by the false labeling of cystoid spaces within choroid in Fig. 5b and d, while in a smoother, lower-contrast choroid saliency maps do not highlight any structures (Fig. 5. Comparision with salinecy maps from the single-task DNNs (Fig.6) to those generated from the multi-task models shows that those single-task saliency maps appear slightly more defined, but generally highlight similar areas.

Figure 5:
  • Download figure
  • Open in new tab
Figure 5:

Exemplary saliency maps for four optical coherence tomography (OCT) images. The first column displays the OCT B-scan with the corresponding labeling of a retinal specialist. Second to fourth column show saliency maps and the network’s confidence for active nAMD (yellow), subretinal fluid (SRF) (cyan) and intraretinal fluid (IRF) (magenta). Note, that saliency maps are only shown in case of confidence > 0.5.

Figure 6:
  • Download figure
  • Open in new tab
Figure 6:

Exemplary saliency maps as in Fig. 5 but results were obtained from single-task models.

4 Discussion

In this study, we developed a multi-task learning model to simultaneously detect SRF, IRF as well as disease activity in OCT B-scans of nAMD patients. We showed that a 3-task model, which takes the presence of IRF and SRF into account to detect disease activity in nAMD, surpassed a single task model regarding accuracy in the activity detection task. Furthermore, our visualization of the multi-task model’s decision-making process via t-SNE showed that inactive and active nAMD B-scans formed different clusters. Among active AMD B-scans, three distinct clusters were observed, which contained OCT B-scans with either SRF, IRF or both fluid types. This separation could not be seen in the single-task models. Saliency maps of exemplary B-scans further corroborate that task-relevant information can be extracted from the multi-task networks.

Thus, multi-task DNN could serve as a basis for an explainable clinical decision support system for nAMD activity, providing support for clinicians in detecting active AMD, but would also allow clinicians to identify evidence in the relevant sub-tasks of finding SRF and IRF. A recent meta-analysis has provided evidence of varying influences of SRF and IRF on the visual outcome in nAMD patients.50 Stable SRF might not affect visual outcome, while fluctuations in IRF during treatment seem to negatively influence visual acuity.50 For this reason, treatment decisions in nAMD solely on a yes/no basis may not meet future treatment guidelines, which might rather require a sophisticated decision depending on the present fluid type and its variation in volume for or against an anti-VEGF injection.

Ophthalmology has recently seen a development ofvarious artificial intelligence systems, yet their use in clinical routine remains rare with only few systems available on the market.51, 52 One big barrier is potential harm of the patient-physician relationship going hand in hand with the lack of trust in those systems.53 Here, we combined multi-task DNNs with different visualization methods to give an insight into the DNNs’ reasoning and increase transparency. First, we used t-SNE as visualization method for high-dimensional data28,29 (Fig. 4) to present the decision-making process of the model. This form of visualization provides an intuitively interpretable rationale for how OCT B-scans were graded by visualizing which other B-scans are similar. The results visualization may also increase an ophthalmologist’s confidence in the model since it illustrates shows that model’s decision making reasoning resembles their own (Fig. 7). In the future, the multi-task system could be extended for other signs of active nAMD such as hard exsudates, pigment epithelial detachment, subretinal hyperreflective materia or hyperreflective foci.4

We further analyzed the multi-task model’s decision on saliency maps ofindividual OCT-scans. Saliency maps highlight critical regions for the model’s decision and thus allow a quick visual control of its reasoning. This may be important in cases of advanced AMD, where fluid is due to degeneration rather than exudation to avoid overtreatment. However, different methods of saliency map agree to differing degrees with clinical annotations30,54,55 and saliency maps can lead to overdiagnosis.18 Therefore we used the saliency map technique with the best clinical relevance for AMD activity30 and displayed saliency maps in case of a confidence of the algorithm > 0.5. Compared to saliency maps of single task DNNs, multi-task saliency maps seem to draw slightly less sharp contours, however, we found good overlap between regions used for active AMD detection and those for SRF and IRF.

Future studies will need to assess how well these multi-task learning results transfer from this data sample acquired at a tertiary center in Germany. It would be desirable to perform similar analysis with larger and more diverse data sets, to thest also the generalization to other populations, different recording qualities as well as OCT devices (including mobile devices). Further, performance could be potentially increased by combining the multi-task network with a a segmentation layer,15 which could reduce false positive cases. Additionally, in clinical routine, activity decision is made on a whole volume not a single B-scan, which could technically be implemented by combining the results from individual B-scans, e.g. by majority voting or uncertainty propagation.

While the approval of anti-VEGF has decreased economic and overall treatment burden of nAMD measured in disability-adjusted life,56,57 a large number of patients still discontinues treatment.58 Patients named the need for assistance, either in the form of a travel companion or a family member, as the main reason for discontinuation.14 Additionally, recurrence of quiescent disease requiring prompt treatment is common, making life-long monitoring necessary.59 For these reasons, automated solutions allowing monitoring close home or even at home are promising technologies:60,61 They provide easier access and reduce the disease burden on the individual.62 Automated solutions for fluid detection have further gained popularity during the Covid-19 pandemic, which showed the devastating effects of delay or interruption of nAMD treatment on visual function.9,59 Despite promising results in laboratory settings, real-world data revealed significantly lower performance rates of home-based OCT with in particular SRF being overlooked by the system.63 This shows the necessity of further developments on the machine learning side to guarantee reliable use, with multi-task learning as suggested in this study being a viable option.

Data Availability

The optical coherence tomography scans were obtained from the University Eye Clinic and their use was permitted by the Institutional Ethics Committee of the University of Tuebingen.

Author contribution statement

MSA, HF and PB designed the research. MSA performed the experiments. GA, WI, FZ, LK were involved in data acuqisition. FZ, HF, GA, LK and WI provided medical advice. MSA, HF and PB wrote the manuscript with input from all authors. All authors approved the final version of the manuscript and agreed on being accountable for the work.

5 Acknowledgments

We thank the German Ministry of Science and Education (BMBF) for funding through the Tübingen AI Center (FKZ 01IS18039A) and the German Science Foundation for funding through a Heisenberg Professorship (BE5601/4-2) and the Excellence Cluster “Machine Learning — New Perspectives for Science” (EXC 2064, project number 390727645). H. Faber thanks the Faculty of Medicine, Eberhard Karls University of Tuebingen, Germany (application number 463–0–0) for additionally funding her research through the Junior Clinician Scientist Program (application number 463–0–0). We further thank Novartis AG for funding part of the research. The funding bodies did not have any influence in the study planning and design.

Appendix

Figure 7:
  • Download figure
  • Open in new tab
Figure 7:

Layer-wise visualization of test data via t-SNE. Starting just before the first Inception module (a) and reading out feature representations yielded by every other module (b-f) along with the last Inception module (g), the shared representation layer (h) and the nAMD activity detection head’s penultimate layer (i), we performed t-SNE with the aforementioned settings. Useful representations emerged towards the end of convolutional stack and the task-specific representation allowed the best separation of nAMD active cases from those inactive. Exact read-out locations can be found in Fig. 8.

Figure 8:
  • Download figure
  • Open in new tab
Figure 8:

Read-out locations within the convolutional stack of the InceptionV3 architecture (indicated by big black arrows). In addition to these, we used the shared representation layer and task-specific layers of our multi-task networks (see Fig. 2). Base figure was obtained from https://cloud.google.com/tpu/docs/inception-v3-advanced.

Footnotes

  • Meeting Presentation: The manuscript was presented as poster presentation (ID PDo04-01) at the 120. Congress of the DOG (Deutsche Ophthalmologische Gesellschaft), 29.09.–02.10.2022, Berlin, Germany

  • Financial Support: Financial support was provided by German Ministry of Science and Education (BMBF) through the Tübingen AI Center (FKZ 01IS18039A) and the German Science Foundation for funding through a Heisenberg Professorship (BE5601/4-2) and the Excellence Cluster “Machine Learning — New Perspectives for Science” (EXC 2064, project number 390727645), the Junior Clinician Scientist Program of the Faculty of Medicine, Eberhard Karls University of Tübingen, Germany (application number 463–0–0) (HF) and the Novartis AG. The sponsor or funding organization had no role in the design or conduct of this research.

  • Conflict of Interest: FZ has received consulting fees from Allergan, Bayer HealthCare, Boehringer-Ingelheim, Novo Nordisk, MSD and Novartis and speaker fees from Alimera, Allergan, Bayer HealthCare and Novartis. FZ was involved in research funded by grants from Bayer Healthcare, Biogen, Clearside, Ionis, Kodiak, Novartis, Ophtea, Regeneron and Roche/Genentech. LK receives, via third-party accounts of the University Eye Hospital, research funding and honoraria from Novartis and research funding from the Tistou and Charlotte Kerstan Foundation. HF received medical training event costs from Novartis. MSA, WI and GA declare no competing interest.

  • Address for preprints: Prof. Dr. rer. nat. P. Berens, Werner Reichardt Centre for Integrative Neuroscience (CIN) Institute for Ophthalmic Research, University of Tübingen, Otfried-Müller-Str. 25, D-72076 Tübingen, Germany. Phone: +49 (0)7071 29-88833, philipp.berens{at}uni-tuebingen.de

  • We updated the manuscript based on peer-review, added new results and discussion points.

  • ↵1 keras.applications.inception_v3.preprocess_input

References

  1. 1.↵
    W. H. Organization, et al. (2019).
  2. 2.↵
    W. L. Wong, et al., The Lancet Global Health 2, e106 (2014).
    OpenUrl
  3. 3.↵
    F. L. Ferris, S. L. Fine, L. Hyman, Archives of ophthalmology 102, 1640 (1984).
    OpenUrlCrossRefPubMedWeb of Science
  4. 4.↵
    R. F. Spaide, et al., Ophthalmology 127, 616 (2020).
    OpenUrlCrossRefPubMed
  5. 5.↵
    P. J. Rosenfeld, Investigative ophthalmology & visual science 57, OCT14 (2016).
    OpenUrlPubMed
  6. 6.↵
    S. Riedl, W.-D. Vogl, S. M. Waldstein, U. Schmidt-Erfurth, H. Bogunovi, Ophthalmology Retina 6, 291 (2022).
    OpenUrl
  7. 7.↵
    S. M. Waldstein, et al., Ophthalmology 123, 1521 (2016).
    OpenUrl
  8. 8.↵
    F. G. Holz, et al., British Journal of Ophthalmology 100, 1623 (2016).
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    A. Valverde-Megías, et al., Journal of Clinical Medicine 10, 3551 (2021).
    OpenUrl
  10. 10.↵
    C. Arruabarrena, et al., Journal of clinical medicine 10, 3281 (2021).
    OpenUrl
  11. 11.↵
    S. Day, K. Acquah, P. P. Lee, P. Mruthyunjaya, F. A. Sloan, American journal of ophthalmology 152, 1014 (2011).
    OpenUrlCrossRefPubMed
  12. 12.
    A. Almony, et al., Journal of managed care & specialty pharmacy 27, 1260 (2021).
    OpenUrl
  13. 13.
    J. M. Ruiz-Moreno, L. Arias, M. J. Abraldes, J. Montero, P. Udaondo, International ophthalmology 41, 3427 (2021).
    OpenUrlCrossRef
  14. 14.↵
    B. Sobolewska, M. Sabsabi, F. Ziemssen, Clinical Ophthalmology (Auckland, NZ) 15, 4317 (2021).
    OpenUrl
  15. 15.↵
    J. De Fauw, et al., Nature medicine 24, 1342 (2018).
    OpenUrlPubMed
  16. 16.↵
    J. Yim, et al., Nature Medicine 26, 892 (2020).
    OpenUrlCrossRefPubMed
  17. 17.↵
    T. D. Keenan, et al., Ophthalmology 128, 100 (2021).
    OpenUrl
  18. 18.↵
    R. Sayres, et al., Ophthalmology 126, 552 (2019).
    OpenUrlCrossRefPubMed
  19. 19.↵
    T. Grote, P. Berens, Journal of medical ethics 46, 205 (2020).
    OpenUrlAbstract/FREE Full Text
  20. 20.↵
    T. Grote, P. Berens, Bioethics 36, 134 (2022).
    OpenUrl
  21. 21.↵
    R. Caruana, Machine learning 28, 41 (1997).
    OpenUrl
  22. 22.↵
    Y. Wang, et al., CoRR abs/2008.05808 (2020).
  23. 23.↵
    Q. Chen, et al., AMIA Summits on Translational Science Proceedings 2019, 505 (2019).
    OpenUrl
  24. 24.↵
    G. Ghahramani, et al., AMIA Annual Symposium Proceedings (American Medical Informatics Associa-tion, 2021), vol. 2021, p. 506.
  25. 25.↵
    S. Sharma, et al., Ophthalmology 123, 865 (2016).
    OpenUrlPubMed
  26. 26.
    S. Llorente-González, et al., Acta Ophthalmologica (2021).
  27. 27.↵
    M. Saenz-de Viteri, et al., Acta Ophthalmologica 99, 861 (2021).
    OpenUrl
  28. 28.↵
    L. Van der Maaten, G. Hinton, Journal of machine learning research 9 (2008).
  29. 29.↵
    D. Kobak, P. Berens, Nature communications 10, 1 (2019).
    OpenUrl
  30. 30.↵
    M. S. Ayhan, et al., Medical Image Analysis p. 102364 (2022).
  31. 31.↵
    C. Szegedy, V. Vanhoucke, S. Ioffe, J. Shlens, Z. Wojna, Proceedings of the IEEE conference on computer vision and pattern recognition (2016), pp. 2818–2826.
  32. 32.↵
    F. Chollet, et al., Keras (2015).
  33. 33.↵
    O. Russakovsky, et al., International journal of computer vision 115, 211 (2015).
    OpenUrl
  34. 34.↵
    Y. E. Nesterov, Dokl. akad. nauk Sssr (1983), vol. 269, pp. 543–547.
    OpenUrl
  35. 35.↵
    I. Sutskever, J. Martens, G. E. Dahl, G. E. Hinton, ICML (3) 28, 5 (2013).
    OpenUrl
  36. 36.↵
    H. Zhang, M. Cisse, Y. N. Dauphin, D. Lopez-Paz, 6th International Conference on Learning Representations (ICLR) (2018), pp. 1–13.
  37. 37.↵
    C. Guo, G. Pleiss, Y. Sun, K. Q. Weinberger, Proceedings of the 34th International Conference on Machine Learning-Volume 70 (JMLR. org, 2017), pp. 1321–1330.
    OpenUrl
  38. 38.
    A. Kendall, Y. Gal, Advances in Neural Information Processing Systems (2017), pp. 5580–5590.
  39. 39.↵
    B. Lakshminarayanan, A. Pritzel, C. Blundell, Advances in Neural Information Processing Systems (2017), pp. 6405–6416.
  40. 40.
    A. Malinin, M. Gales, Advances in Neural Information Processing Systems (2018), pp. 7047–7058.
  41. 41.↵
    Y. Ding, J. Liu, J. Xiong, Y. Shi, arXiv preprint arXiv:1903.02050 (2019).
  42. 42.↵
    V. Bhise, et al., Journal of general internal medicine 33, 103 (2018).
    OpenUrlCrossRefPubMed
  43. 43.↵
    S. Thulasidasan, G. Chennupati, J. A. Bilmes, T. Bhattacharya, S. Michalak, Advances in Neural Information Processing Systems 32 (2019).
  44. 44.↵
    S. Fort, H. Hu, B. Lakshminarayanan, arXiv preprint arXiv:1912.02757 (2019).
  45. 45.
    Y. Ovadia, et al., Advances in Neural Information Processing Systems (2019), pp. 13991–14002.
  46. 46.↵
    M. S. Ayhan, et al., Medical Image Analysis p. 101724 (2020).
  47. 47.↵
    P. G. Poliar, M. Stražar, B. Zupan, bioRxiv (2019).
  48. 48.↵
    S. Bach, et al., PloS one 10, e0130140 (2015).
    OpenUrlCrossRefPubMed
  49. 49.↵
    J. N. Böhm, P. Berens, D. Kobak, Journal of Machine Learning Research 23, 1 (2022).
    OpenUrl
  50. 50.↵
    U. Chakravarthy, et al., Eye 35, 2983 (2021).
    OpenUrlPubMed
  51. 51.↵
    M. D. Abràmoff, P. T. Lavin, M. Birch, N. Shah, J. C. Folk, NPJ digital medicine 1, 1 (2018).
    OpenUrlCrossRef
  52. 52.↵
    A. Olvera-Barrios, et al., British Journal of Ophthalmology 105, 265 (2021).
    OpenUrlAbstract/FREE Full Text
  53. 53.↵
    J. J. Hatherley, Journal of Medical Ethics 46, 478 (2020).
    OpenUrlAbstract/FREE Full Text
  54. 54.↵
    1. H. Fu,
    2. M. K. Garvin,
    3. T. MacGillivray,
    4. Y. Xu,
    5. Y. Zheng, eds
    A. Singh, et al., Ophthalmic Medical Image Analysis, H. Fu, M. K. Garvin, T. MacGillivray, Y. Xu, Y. Zheng, eds. (Springer International Publishing, Cham, 2020), pp. 21–31.
  55. 55.↵
    T. Van Craenendonck, B. Elen, N. Gerrits, P. De Boever, Translational vision science & technology 9, 64 (2020).
    OpenUrl
  56. 56.↵
    X. Xu, et al., BMC Public Health 20, 1 (2020).
    OpenUrlCrossRefPubMed
  57. 57.↵
    K. Mulligan, et al., JAMA ophthalmology 138, 40 (2020).
    OpenUrl
  58. 58.↵
    I. Westborg, et al., Acta Ophthalmologica 95, 787 (2017).
    OpenUrlCrossRefPubMed
  59. 59.↵
    S. Aslanis, U. Amrén, C. Lindberg, D. Epstein, Ophthalmology Retina 6, 15 (2022).
    OpenUrl
  60. 60.↵
    F. A. Sloan, G. Picone, D. S. Brown, P. P. Lee, Journal of the American Geriatrics Society 53, 1867 (2005).
    OpenUrlCrossRefPubMedWeb of Science
  61. 61.↵
    E. Y. Chew, et al., Contemporary clinical trials 37, 294 (2014).
    OpenUrl
  62. 62.↵
    K. Nahen, G. Benyamini, A. Loewenstein, Klinische Monatsblätter für Augenheilkunde 237, 1410 (2020).
    OpenUrl
  63. 63.↵
    Y. Liu, N. L. Holekamp, J. S. Heier, Ophthalmology Retina (2022).
Back to top
PreviousNext
Posted November 23, 2022.
Download PDF
Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
Multi-task learning for activity detection in neovascular age-related macular degeneration
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Multi-task learning for activity detection in neovascular age-related macular degeneration
Murat Seçkin Ayhan, Hanna Faber, Laura Kühlewein, Werner Inhoffen, Gulnar Aliyeva, Focke Ziemssen, Philipp Berens
medRxiv 2022.06.13.22276315; doi: https://doi.org/10.1101/2022.06.13.22276315
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools
Multi-task learning for activity detection in neovascular age-related macular degeneration
Murat Seçkin Ayhan, Hanna Faber, Laura Kühlewein, Werner Inhoffen, Gulnar Aliyeva, Focke Ziemssen, Philipp Berens
medRxiv 2022.06.13.22276315; doi: https://doi.org/10.1101/2022.06.13.22276315

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Ophthalmology
Subject Areas
All Articles
  • Addiction Medicine (349)
  • Allergy and Immunology (668)
  • Allergy and Immunology (668)
  • Anesthesia (181)
  • Cardiovascular Medicine (2648)
  • Dentistry and Oral Medicine (316)
  • Dermatology (223)
  • Emergency Medicine (399)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (942)
  • Epidemiology (12228)
  • Forensic Medicine (10)
  • Gastroenterology (759)
  • Genetic and Genomic Medicine (4103)
  • Geriatric Medicine (387)
  • Health Economics (680)
  • Health Informatics (2657)
  • Health Policy (1005)
  • Health Systems and Quality Improvement (985)
  • Hematology (363)
  • HIV/AIDS (851)
  • Infectious Diseases (except HIV/AIDS) (13695)
  • Intensive Care and Critical Care Medicine (797)
  • Medical Education (399)
  • Medical Ethics (109)
  • Nephrology (436)
  • Neurology (3882)
  • Nursing (209)
  • Nutrition (577)
  • Obstetrics and Gynecology (739)
  • Occupational and Environmental Health (695)
  • Oncology (2030)
  • Ophthalmology (585)
  • Orthopedics (240)
  • Otolaryngology (306)
  • Pain Medicine (250)
  • Palliative Medicine (75)
  • Pathology (473)
  • Pediatrics (1115)
  • Pharmacology and Therapeutics (466)
  • Primary Care Research (452)
  • Psychiatry and Clinical Psychology (3432)
  • Public and Global Health (6527)
  • Radiology and Imaging (1403)
  • Rehabilitation Medicine and Physical Therapy (814)
  • Respiratory Medicine (871)
  • Rheumatology (409)
  • Sexual and Reproductive Health (410)
  • Sports Medicine (342)
  • Surgery (448)
  • Toxicology (53)
  • Transplantation (185)
  • Urology (165)