Skip to main content

ORIGINAL RESEARCH article

Front. Neuroanat., 13 December 2021
This article is part of the Research Topic How Animals See: Structure and Function of Light Sensory Tissues Along Evolution View all 12 articles

Characterization of the Canine Retinal Vasculature With Optical Coherence Tomography Angiography: Comparisons With Histology and Fluorescein Angiography

  • 1Division of Experimental Retinal Therapies, Department of Clinical Sciences & Advanced Medicine, School of Veterinary Medicine, University of Pennsylvania, Philadelphia, PA, United States
  • 2Department of Pathobiology, School of Veterinary Medicine, University of Pennsylvania, Philadelphia, PA, United States
  • 3Department of Ophthalmology, Scheie Eye Institute, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, United States
  • 4Department of Physiology, Genetics and Microbiology, University of Alicante, Alicante, Spain

Purpose: To present a methodology for quantification of the canine retinal vasculature imaged by optical coherence tomography angiography (OCTA) and validate this approach by comparison with fluorescein angiography (FA) and confocal imaging of retinal wholemounts labelled by immunohistochemistry (IHC).

Methods: Six normal adult dogs underwent retinal OCTA imaging in both eyes. The images extracted from the different microvascular plexuses at eight retinal locations spanning the central and mid-peripheral fundus were analyzed using the AngioTool software. FA was performed in one eye and was compared to the OCTA images. Six eyes from three dogs were processed by IHC to examine the retinal vasculature.

Results: A total of four retinal plexuses were identified by OCTA in the canine retina, and their density and topographical pattern varied with eccentricity. OCTA offered improved resolution over FA with the advantage of allowing imaging of the individual plexuses. Detection by OCTA of small vessels within the deep capillary plexus was possible and approached the level of resolution achieved with ex vivo imaging of the retinal vasculature by confocal microscopy/IHC. The plexuses herein described are analogous to human retinal vasculature.

Conclusion: OCTA can be used to image and quantify non-invasively the vascular retinal networks of the canine retina. We provide normative data in eight different retinal locations that can be imaged non-invasively with this technology. This could support analysis of retinal vascular changes associated with disease and following therapeutic intervention.

Introduction

The retina is a complex and highly metabolic extracranial part of the central nervous system that requires a continuous and self-regulated blood supply (Harris et al., 1998; Yu and Cringle, 2001). In vascularised retinas, oxygenation is ensured by a well-organised retinal and choroidal vascular network. The choriocapillaris, a network of capillary vessels located under the retinal pigmented epithelium, is the main source of oxygen for the outer retina, and retinal vessels are the main supply for the inner retina (Michaelson, 1954). While choroidal vasculature has been maintained throughout evolution in vertebrates, the retinal vasculature pattern differs widely between species (Rochon-Duvigneaud, 1943). Based on the presence and distribution of this vasculature, four patterns have been described: euangiotic/holangiotic, merangiotic, paurangiotic and anangiotic (Leber, 1875; Schaepdrijver et al., 1989). The canine retina is classified as holangiotic as its blood vessels extend from the optic nerve head to the far periphery. In the temporally located area centralis, they converge toward the highly specialized fovea-like area (Peichl, 1992; Beltran et al., 2014).

However, the vasculature of the canine retina has not been thoroughly studied, and the location of the different vascular networks, as well as their nomenclature, remains unclear (Leber, 1875; Barrett, 1886; Mutlu and Leopold, 1964). In dogs, previous studies have imprecisely identified the presence of radially distributed blood vessels surrounding the optic nerve papilla, parallel to the nerve fibre layer (NFL) (Parry, 1953; Mutlu and Leopold, 1964). A second vascular network composed of major arterioles and venules has also been described and identified within the ganglion cell layer (GCL) (Parry, 1953; Mutlu and Leopold, 1964). Additionally, two capillary beds have been identified in the inner plexiform layer (IPL) (Schaepdrijver et al., 1989) and inner half of the outer plexiform layer (OPL) (Barrett, 1886; Parry, 1953; Mutlu and Leopold, 1964; Engerman et al., 1966; Schaepdrijver et al., 1989). However, the techniques used to describe the canine retinal microvasculature to date, such as vascular corrosion casting (Fischer and Slatter, 1978; Schaepdrijver et al., 1989), do not include structural analysis of the exact location of the vasculature in relation to the retinal layers. Furthermore, there are discrepancies between studies in the number of plexuses described (Michaelson, 1954). In contrast, four vascular plexuses have been identified and extensively described in the human retina (Chan et al., 2012).

Evaluation and quantification of the retinal vasculature is a valuable method to diagnose and monitor the progression of many retinal diseases that can result in blindness. This information serves as an essential diagnostic, prognostic, and therapeutic tool (Yao et al., 2020). Fluorescein angiography (FA) is currently the most widely applied technique in human and veterinary medicine for evaluation of retinal vessels in vivo (Gelatt et al., 1976; Cavallerano, 1996). Conventional FA has several disadvantages such as: suboptimal resolution in visualization of deeper retinal capillaries, and the need to use an intravenous contrast agent which may cause adverse systemic reactions (Kwiterovich et al., 1991; López-Sáez et al., 1998). Nonetheless, FA can detect increased vessel permeability evidenced by dye leakage, which is pathognomonic of blood-retinal barrier breakdown.

Optical coherence tomography angiography (OCTA) is a novel non-invasive technology that allows detailed visualization of retinal circulation, enabling the study of the vascular plexuses separately using their intrinsic motion decorrelation. The OCTA probabilistic full-spectrum amplitude decorrelation algorithm generates images based on the motion signal between a series of OCT b-scans acquired at the same position. The differences between the images acquired are due to erythrocyte movement in the blood vessels. By focusing on the movement, the algorithm is able eliminate all static sections of the scan area in order to generate a resulting image. This image is then further processed to produce the final en face angiogram (Coscas et al., 2016). Hence, it is able to detect vessel networks using their intrinsic motion decorrelation, thus avoiding the need for an intravenous dye (Jia et al., 2012).

While there are many different techniques and parameters utilized to quantify retinal vasculature, vessel density (VD) is the most widely used across studies (Pechauer et al., 2015). VD is generally defined as the percent area of an image that is occupied by blood vessels. However, the absolute values can vary widely between studies and there is a need for validation of novel OCTA-based approaches in canine retinas. Currently, the most recognized method of validating OCTA images is via direct comparison with retinal wholemounts processed for detection of the retinal vasculature by immunohistochemistry (IHC) (Yu et al., 2021). Ultimately, in order to be considered a reliable imaging technique, OCTA must recognize the majority of the vessels that are visualized by IHC.

The main objective of this study was to establish a methodology and normative data for retinal vasculature quantification by OCTA in canine eyes. We also aimed to validate our proposed methodology by comparison with FA and IHC. In order to achieve these goals, we investigated the anatomic features of the different retinal plexuses, their interconnections, and how their specific topographic patterns vary as a function of eccentricity from the optic nerve head.

Materials and Methods

Study Animals

The animals enrolled in the study were part of the Retinal Diseases Studies Facility research colony at the University of Pennsylvania. All the procedures were approved by the Institutional Animal Care and Use Committee of the University of Pennsylvania (IACUC# 804956).

The study population consisted of 6 normal, mesocephalic, mongrel dogs (12 eyes) with a mean age of 4.5 ± 1.3 years of age, and both genders were equally represented. All the animals underwent OCTA imaging in both eyes (OCTA module, Spectralis HRA + OCT2, version 6.9.4.0, Heidelberg Engineering Inc., Franklin, MA, United States) prior to FA being performed in one eye at the end of the procedure. Three animals, with both genders represented, where humanely euthanized with an overdose of a pentobarbital-based euthanasia solution (Euthasol, Virbac, Westlake, TX, United States) administered via intravenous injection, and their retinas were processed for IHC analysis, as previously described (Beltran et al., 2014).

For the in vivo imaging, the pupils were dilated with atropine sulfate 1%, tropicamide 1% (both drugs from Akorn Inc., Lake Forest, IL, United States) and phenylephrine 10% (Paragon Biotech, Portland, OR, United States). The dogs were anesthetised by intravenous propofol induction (2–6 mg/kg, Zoetis, Kalamazoo, MI, United States) and maintained with gas inhalation (isoflurane 2–3%, Akorn Inc., Lake Forest, IL, United States). Once positioned in sternal recumbency, the eyelids were kept open with an eye speculum and two stay sutures (Vicryl 4-0, Ethicon Inc., Somerville, NJ, United States). The stay sutures were placed in the conjunctiva 1–2 mm away from the limbus at the 2 and 10 o’clock position in order to gently manipulate the orientation of the visual axis as required to image the entire retina. The ocular surface was lubricated frequently with saline (Sodium Chloride 0.9%, ICU Medical, Inc., Lake Forest, IL, United States) and the anaesthetic plane was closely monitored via heart and respiratory rate, electrocardiography, blood pressure and oxygen saturation measurements.

Optical Coherence Tomography Angiography Characterisation of the Retinal Vasculature

The retinal vascular characteristics were studied and quantified by obtaining 10-degree square OCTA images along the superior, inferior, temporal, and nasal retinal quadrants. For study purposes, 8 areas were preselected within the central and mid-peripheral retina (see Table 1). The mid-periphery scans corresponded to the most eccentric locations that could be reliably imaged with the OCTA unit. Indeed, the position of the globes in the orbit, and the shape of the canine skull precluded imaging of the peripheral retina. An additional area centered on the optic nerve head (ONH) was included for qualitative analysis only. The first scan was acquired in the area centralis and centered on the canine fovea-like area whose exact location was confirmed by anatomical features previously described including ONL thinning and a centrally converging vascular pattern (Supplementary Figure S1) (Beltran et al., 2014). Additionally, the orientation of the NFL on en face OCT shows a particular pattern bending around the very center of this cone-rich region (Supplementary Figure S1A), and this newly identified feature was used to accurately position the OCTA scanning window.

TABLE 1
www.frontiersin.org

Table 1. Location of the areas of interest and their distance from the optic nerve.

Each imaged area was composed of volumetric scans that included 512 a-scans and 512 b-scans with 6 μm distance in between b-scans. This provided a lateral resolution of 5.75 μm/pixel and a depth resolution of 3.87 μm/pixel. In order to account for curvature bias, we calculated the retinal magnification factor based on the individual axial globe length and the following formula RMF = 2*π*(0.59*AGL)/360, as previously described (Beltran et al., 2014). The axial globe length was measured in all 12 eyes after each imaging session by ultrasound (Sonomed A-scan A1500, Sonomed-Escalon, Lake Success, NY, United States), and the mean ± SD was 21.2 ± 0.5 mm. By using this approach, we calculated that the 10 × 10 degree OCTA images corresponded in these dogs to a 2.2 ± 0.05 mm2 area. The exact area for each eye, was used to calibrate the AngioTool software (see below) and obtain accurate quantifications.

In order to individually image the vascular networks, the boundaries of the plexuses were identified and selected by adjusting the location where the algorithm is applied on the OCT b-scan, with the designated software (Heidelberg Eye Explorer, HEYEX). In order to define the most superficial individual retinal plexus the first boundary of the slab was placed on the ILM and the second segmentation line was manually displaced in depth until no vessels were captured on the en face OCTA. The upper boundary was then displaced to the endpoint of the previous one, and the same process was repeated to define the following plexus. In this way, we not only identified different plexuses but also described their location based on the retinal lamination seen by OCT.

Quantitative analysis was performed in all the areas of interest, with the exception of the ONH scan (Table 1). The OCT angiograms were extracted as .tiff files and then imported into the semiautomated vessel-analysis program AngioTool (Zudaire et al., 2011). This is a validated software, accessible in the public domain1 (0.6a version, National Cancer Institute, Bethesda, MD, United States). AngioTool processes images in order to obtain a set of quantitative parameters and has previously been used to quantify retinal vasculature in other species (Zudaire et al., 2011; Giannakaki-Zimmermann et al., 2016; Liu et al., 2017; Munk et al., 2017; Sun et al., 2021; Told et al., 2021). Once the angiogram is uploaded into AngioTool, the program allows adjustments to the vessel diameter and intensity to ensure that the resulting image matches the correct vessel profile. Once the analysis is completed, the program extracts the results in a working spreadsheet. For each analyzed slab, we selected among the AngioTool parameters: VD, Junction Density (JD), Average Vessels Length (AVL), and Lacunarity (LC). For more details, see Table 2.

TABLE 2
www.frontiersin.org

Table 2. AngioTool parameters used to characterise the retinal vasculature.

Validation of Optical Coherence Tomography Angiography Quantifications by Comparison With Fluorescein Angiography and Immunohistochemistry

Fluorescein angiography images were obtained from the last eye evaluated (left eye) in all dogs. The FA images were obtained via fluorescein injection into the cephalic vein (20 mg/kg of Fluorescein sodium 10%, AK-FLUOR, Akorn, Lane Forest, IL, United States). Several cSLO fundus photographs were acquired with a short wavelength laser (488 nm, Spectralis Blue Autofluorescence mode) at the completion of the venous phase defined by a homogeneous and complete filling of the retinal venules (Schaepdrijver et al., 1996). The images were acquired with a 30 degrees lens that provided a lateral resolution of 5.60 μm/pixel.

For wholemount preparation, the retinas were prepared as previously described (Beltran et al., 2014). Briefly, after 15 min of fixation in 4% paraformaldehyde, the neuroretina was separated from the RPE and the fixation continued for 4 h at 4°C. The tissues were incubated for 4 days at 4°C in a solution containing primary antibodies directed against Collagen IV (1:100, Millipore Sigma AB769, Burlington, MA, United States). Fluorochrome-labelled secondary antibodies (Alexa Fluor 488 and 568 dyes; 1:200, Thermo Scientific) and DAPI nuclear stain were applied for 24 h at 4°C. The retinas were imaged using a 2-photon confocal microscope (Leica SP8 Multiphoton, Leica Microsystems, Wetzlar, Germany) with a 20x (HCX APO L, 1.0 NA) water immersion objective lens. The captured images studied were composed of individual 0.5 mm2 fields, acquired from z stacks spanning from the ILM to the ONL. Z-stacks were acquired at 1.5 μm Z-steps in 1024 × 1024 pixel format at 400 Hz with a line average of 2. These settings offered a lateral resolution of 0.541 μm/pixel and axial resolution of 1.5 μm/pixel. The images were acquired, merged, and analyzed with Leica Application Suite (LAS X, 3D Visualisation Module).

The evaluation of the anatomical features of the plexuses by IHC was performed in three eyes by stitching tiled z-stacks, covering an area that extended from the ONH to the ora serrata in all four retinal quadrants. Using the 3D representation module of the Leica Application Suite, we measured the distance from the optic nerve to the area where each of the retinal vascular plexuses was no longer present, to identify the topographical distribution of each vascular network.

For quantitative analysis, the OCTA images were matched with those obtained with FA and IHC using a graphic editing software (Adobe Photoshop CC, Adobe Systems Inc., CA, United States) and by matching vascular landmarks in the areas of interest. For this purpose, the individual 0.5 mm2 fields on IHC were stitched to cover the same squared area that was previously imaged on OCTA in the eight areas of interest. From this merged three-dimensional stack, the plexuses were extracted by cropping the subset of data that contained the same vascular networks defined by OCTA. Once the individual stacks were extracted, the mass projection image was exported as a .tif file and used for quantification using AngioTool, as previously described.

Statistical Analysis

All statistical analyses were performed using SAS v9.4 (SAS Institute, Cary, NC, United States). For all tests, two-sided p ≤ 0.05 was considered statistically significant.

The mean (standard deviation, SD) of each OCTA measurements from the different locations within each quadrant was calculated, and 95% normal limits were derived as mean ± 1.96*SD. For the statistical comparison of each OCTA measures between mid-periphery vs. central areas, generalised linear models were fitted and generalised estimating equations (GEE) were used to account for the inter-eye correlation for measures from two eyes of the same dog and correlation of measures at different locations from the same eye. In the GEE, the unstructured covariance structure was used.

For comparison of measures between imaging techniques (OCTA vs. FA and OCTA vs. IHC) that were measured in one eye per dog, paired t-tests were used to test their statistical significance.

Results

Optical Coherence Tomography Angiography Enables Identification of Four Vascular Plexuses in the Canine Retina

Manual segmentation in the OCTA module of the HEYEX software enabled precise localisation of individual retinal vascular plexuses within the retinal layers by concurrent visualisation of both the OCTA b-scan and the en face view of the angiogram (Figure 1). The most inner vascular network found within the NFL extended only into the peripapillary region and did not reach central areas. This plexus shared similarities with the human Radial Peripapillary Capillary Network (RPCP; Figures 1A1,2). The second vascular network spanned from the ILM to the GCL/IPL border, was formed by large arterioles, venules and their connecting capillaries, and branched from the ONH to the mid-periphery (most eccentric areas scanned by OCTA). It was analogous to the Superficial Vascular Plexus described in human retinas (SVP; Figures 1B1,2). We did not find an avascular zone in the canine SVP at the corresponding location of the canine fovea-like area, as is commonly seen in species with a foveal pit. The third vascular network was located within the IPL in proximity with the INL and contained capillaries which formed a homogeneous mesh. This vascular network was defined as the Intermediate Capillary Plexus (ICP; Figures 1C1,2). The SVP and ICP were seen in very close proximity, and they were frequently extracted together (Figures 1C1,2). The fourth and outermost plexus was located within the inner half of the OPL and formed a dense network of capillaries. This plexus was similar to the human Deep Capillary Plexus (DCP; Figures 1E1,2).

FIGURE 1
www.frontiersin.org

Figure 1. Illustration of the vascular plexuses of the canine retina imaged by OCTA. (A1,2) Radial peripapillary capillary plexus (RPCP). This plexus is only seen within the NFL (red arrow) around the optic nerve head. (B1,2) Superficial vascular plexus (SVP) in the area centralis. (C1,2) Intermediate capillary plexus (ICP) in the same area. (D1,2) SVP and ICP extracted together. (E1,2) Deep capillary plexus (DCP). (F1,2) OCT angiogram containing all the retinal vascular networks ( = whole retina thickness). The red dotted lines in panels (A1–F1) delineate the inner and outer boundaries of the slab selected to extract the angiograms shown in panels (A2–F2) (green dashed lines). NFL, nerve fibre layer; OS, left eye; ONH, optic nerve head; GCL, ganglion cell layer; IPL, inner plexiform layer; OPL, outer plexiform layer; ILM, inner limiting membrane; RPE, retinal pigment epithelium.

A previously undescribed feature seen in 9 out of 12 eyes consisted of a small focal area with a lower density of vessels in the DCP at the centre of the canine fovea-like area (Supplementary Figure S2A). We confirmed by IHC that this same anatomical feature was visible on retinal wholemounts in the DCP but not in the SVP (Supplementary Figure S2B).

Vessel Density and Other Vascular Optical Coherence Tomography Angiography Parameters Vary With Retinal Eccentricity

Due to the presence of large vessels in the SVP that occasionally invade the ICP, individual differentiation and segmentation of these two plexuses was challenging. Consequently, both plexuses were extracted together in the same imaging slab to be quantified as a single image (Figures 1D1,2). The DCP was extracted next, and then the last slab, whose limitations were set at the ILM and RPE (whole retina, WR; Figures 1F1,2). Representative images from each area and slabs used for quantification are compiled in Figure 2.

FIGURE 2
www.frontiersin.org

Figure 2. Representative OCTA images of the areas and slabs used for quantification. SVP, superficial vascular plexus; ICP, intermediate capillary plexus; DCP, deep capillary plexus.

The 95% confidence intervals of OCTA parameters for each slab (SVP + ICP, DCP and WR) were calculated, providing a range of reference values for this adult dog cohort (Table 3).

TABLE 3
www.frontiersin.org

Table 3. Normal range of the vascular parameters in each retinal area.

Along each quadrant, the VD was consistently lower in the mid-periphery in comparison to more central areas. This was observed for both the SVP + ICP and the DCP, although not always reaching statistical significance (Figures 3A1,B1,C1). JD values followed a similar topographical pattern, being higher in the more central areas than in the mid-periphery. This reached statistical significance in all four quadrants of the DCP (Figures 3A2,B2,C2). AVL results also showed a similar topographical trend, with mid-peripheral vessels being shorter due to decreased tortuosity (Figures 3A3,B3,C3). Finally, LC values were found to be higher in mid-peripheral than in central areas, suggesting increased vascular heterogeneity toward the periphery (Figures 3A4,B4,C4).

FIGURE 3
www.frontiersin.org

Figure 3. Topographical quantification of vascular plexuses by OCTA in the canine retina. The vessel density (A1,B1,C1), junction density (A2,B2,C2), average vessel length (A3,B3,C3) and lacunarity (A4,B4,C4) values are provided for the superficial vascular plexus and intermediate capillary plexus (SVP+ICP, A1–4), deep capillary plexus (DCP, B1–4) and whole retina slabs (WR, C1–4). A.C., area centralis; T, temporal; N, nasal; S: superior; I, inferior; mid-per, mid-periphery. Results are shown as a mean ± SD (n = 12). Paired t-test: *p ≤ 0.05, **p ≤ 0.01, ***p ≤ 0.001.

Analysis of the quality of the OCTA scans using signal to noise ratio values provided by the Spectralis software (HEYEX) showed that the best images were acquired in the area centralis and the poorest in the inferior mid-periphery. For more details see Supplementary Table S1.

Optical Coherence Tomography Angiography Supersedes Fundus Fluorescein Angiography at Detecting Small Vascular Networks

Fluorescein angiography images were acquired at the late venous phase, which was seen at 18 ± 4.7 s after contrast agent injection. Since FA does not allow visualisation of individual vascular networks, the images were compared to the OCTA image extracted from the combination of all the slabs (WR). Qualitatively, OCT angiograms had a greater ability to resolve smaller vascular branches and capillaries than FA (Supplementary Figure S3A) and this was confirmed quantitatively when values of vascular parameters (VD, JD, AVL, and LC) were compared between the two techniques (Supplementary Figure S3B).

Immunohistochemistry Validates Optical Coherence Tomography Angiography as a Reliable Technique to Image Retinal Vessels

To confirm the presence of the four plexuses identified by OCTA, retinal wholemounts from three previously OCTA-imaged eyes were processed for IHC with an antibody (Collagen IV) that labels vessels and imaged by 2-photon/confocal microscopy (Figure 4). A pan-retinal three-dimensional representation of the plexuses and their inter-connections was evaluated, in regions that extended from the ONH to the periphery (ora serrata) in each quadrant (Figure 5A and Supplementary Video S1). The SVP contained large arterioles and venules that radiated from the ONH and diverged off, creating tangential branches above and below to form the closely related ICP and RPCP (Figure 5A1). The vessels in the ICP were also seen to branch out and cross through the INL forming the DCP (Figure 5A1,2, white arrows). The DCP is a dense plexus of capillaries and the closest vascular supply to the photoreceptors somatas in the ONL (Figure 4A5). The DCP approaches the SVP towards the periphery (Figure 5A3) to finally merge with it near the ora serrata (Figure 5A4). In the two dense capillary networks (ICP and DCP), intercapillary bridging cells were identified and positively labelled with Collagen IV (Supplementary Video S1). No vascular structures were found beyond the ONL (Figure 4A6).

FIGURE 4
www.frontiersin.org

Figure 4. Canine retinal vascular plexuses imaged by 2-photon/confocal microscopy/IHC. (A1,2) En face and cross-sectional visualisation of the superficial vascular plexus (SVP), (A3) the intermediate capillary plexus (ICP), and (A5) the deep capillary plexus (DCP). (A4) The INL contains connecting capillaries from the ICP to the DCP. (A6) The ONL is avascular.

FIGURE 5
www.frontiersin.org

Figure 5. Topographical analysis by IHC of the vascular networks on canine retinal wholemounts. (A) Pseudocolored three-dimensional representation of the vascular plexuses imaged by 2 photon/confocal microscopy at four locations extending from the peripapillary region to the ora serrata. In the peripapillary region (A1), all four plexuses are identified (RPCP, radial peripapillary capillary plexus; SVP, superficial vascular plexus; ICP, intermediate capillary plexus; DCP, deep capillary plexus). In the central and mid-peripheral retina (A2), the RPCP is no longer seen. At the periphery (A3) only the SVP and DCP remain, and both merge at the far periphery close to the ora serrata (A4). The right-side portion of each panel shows localisation of the Collagen IV vascular marker (in red) and DAPI (in blue) for improved visualisation of the retinal layers in which the plexuses are located. (B) Graphic representation of the topographical extension from the ONH of each vascular plexus measured in three individual retinal wholemounts (B1) and averaged as percentage of retinal length (ONH to ora serrata distance) in (B2). Numerical values shown on the retinal diagrams (B1) represent the distance from the centre of the ONH to the point where a change in the number of plexuses in seen (B2).

It is noteworthy that, these four plexuses were not identified in all locations of the retina (Figure 5B1). Peripapillary, the RPCP was only found to extend within 4–5.3 mm from the center of the ONH (Figure 5B2 and Table 4). Although easily identifiable in the central retina, the ICP becomes closer to the SVP in the periphery, precluding its independent visualisation. It was found to merge with the SVP within 9–15.6 mm from the centre of the ONH (Figure 5B and Table 4). In the far periphery, the DCP merged with the SVP (Figure 5B).

TABLE 4
www.frontiersin.org

Table 4. Variation in the number of plexuses per quadrant, imaged by 2-photon confocal microscopy.

Qualitative comparison of both imaging techniques showed that small capillaries that were well resolved by confocal microscopy/IHC could also be detected by OCTA, albeit with a lower resolution (Figure 6). Images of the SVP + ICP obtained by IHC showed more precisely the smaller capillaries within the ICP, thus resulting in higher VD and JD values in all retinal quadrants, although statistical significance was not always reached, likely due to the small sample size (Supplementary Figure S4). Since the homogeneous ICP was better imaged by IHC than OCTA, LC values were lower in IHC images of the SVP + ICP (Supplementary Figure S4). It is worth noting that small vessels within the DCP were overall well detected with OCTA, and VD values were not found to be significantly different when comparing the two imaging approaches (Supplementary Figure S4). Images from the “Whole retina” slab followed the same trend described for SVP + ICP (Supplementary Figure S4) with a higher VD detected by IHC.

FIGURE 6
www.frontiersin.org

Figure 6. Qualitative comparison of retinal vascular network imaging in dogs by OCTA and IHC. In vivo OCTA (A1,B1,C1) approaches the level of resolution seen with confocal microscopy/IHC on retinal wholemounts for detection of small capillaries (A2,B2,C2). SVP, superficial vascular plexus; ICP, intermediate capillary plexus; DCP, deep capillary plexus; OS, left eye.

Discussion

Optical Coherence Tomography Angiography Reveals Four Vascular Plexuses in the Canine Retina That Share Similarities With Vascular Networks Reported in Human

The retinal distribution of the main arterioles, venules and capillary networks has been previously, albeit imprecisely, described in dogs through the use of vascular corrosion casting combined with light and electron microscopy (Engerman et al., 1966; Flower et al., 1985; Schaepdrijver et al., 1989). However, the number of plexuses, as well as their nomenclature, differs between canine studies (Engerman et al., 1966; Schaepdrijver et al., 1989). With the use of modern imaging techniques, this work has confirmed the presence of four distinct canine retinal plexuses. Since the topographical pattern and location of these plexuses within the retinal layers is analogous to that found in human retinas, we propose that the same nomenclature (RPCP, SVP, ICP, and DCP) outlined by Campbell et al. (2017) be used in all future description of canine retinas (Cuenca et al., 2020).

Similar to the human retinal vasculature, the holangiotic canine retina has a dual blood supply, with oxygenation of the inner retinal layers provided by the four vascular plexuses, and the outer retina relying on the choriocapillaris (Michaelson, 1954). This contrasts with the vascularisation of other vertebrates in which the entire retina relies quasi exclusively on choroidal blood supply (Rochon-Duvigneaud, 1943). Although no studies have reported measurements of intraretinal oxygen tension levels in the dog, it is likely that it shares a similar profile as that of other holangiotic species (Braun et al., 1995). In holangiotic retinas, the outer retinal oxygen levels gradually decrease as a function of distance from the choriocapillaris, and in the inner retina higher levels are found in close proximity to the retinal plexuses, reflecting the dual blood supply. In contrast, in anangiotic or avascular areas of merangiotic retinas, the oxygen values decrease rapidly as the distance from the choriocapillaris increases, reaching very low levels across the majority of the inner retina (Yu and Cringle, 2001). Our detailed analysis of the canine retinal vasculature, particularly the description of the exact location of each plexus within retinal layers, will help improve further understanding of the physiologic and pathologic intraretinal oxygen tension distribution as well as the energy demand and metabolic pathways throughout the various retinal layers in this species.

It is worth noting that while dogs possess a fovea-like region that contains a density of cones comparable to that of primates (Beltran et al., 2014), we did not identify a well-defined foveal avascular zone (FAZ) as described in human eyes (Chui et al., 2012). Nevertheless, in our small cohort of dogs, a poorly defined area with a lower density of vessels was seen in the DCP. Interestingly, this closely resembles the vascular pattern seen in people with a fovea plana. In this condition, the lack of a FAZ and the absence of foveal pit do not affect vision (Marmor et al., 2008; Dolz-Marco et al., 2016; Cicinelli et al., 2017; Chatzistergiou et al., 2021). Hence, the lack of a FAZ in dogs is likely associated with the anatomy of the canine fovea-like area that lacks a foveal pit.

Optical Coherence Tomography Angiography Can Be Used to Quantify Vascular Parameters in Dogs

Although many studies have reported quantitative data for OCTA images in human retinas, no canine vascular parameters have been described. This study provides normative data for the SVP + ICP, DCP, and WR, successfully establishing a reference range that can be consulted and compared in future studies.

In human eyes, the number of retinal plexuses and their densities diminish toward the periphery (Campbell et al., 2017; Cuenca et al., 2020; Park et al., 2021), and this is similar to what Engerman et al. (1966) previously described in dogs. Our work not only confirms this finding, but now provides quantitative data for four parameters that are frequently used to characterise vascular networks.

Interestingly, we found the VD, JD, and homogeneity (lowest LC) of the DCP to be highest within the area centralis. This is very similar to what occurs in the human fovea and this finding indicates that this region of the central canine retina shares closer homology to the primate fovea than previously described (Chan et al., 2012).

Optical coherence tomography angiography signal strength affects the detection of smaller vessels, and Yu et al. (2019) confirmed that a low signal strength can underestimate VD. This highlights the importance of including signal strength or quality parameters when reporting quantitative data (Al-Sheikh et al., 2017). Straightforward ways to optimise quality include, frequently lubricating the ocular surface as well as minimising eye movements through adjustment of the anaesthetic plane to minimise projection and motion artefacts.

Optical Coherence Tomography Angiography Is a Valid Method for Non-invasive in vivo Characterisation of Retinal Vascular Networks in Dogs

We validated the OCTA results by comparison with two other imaging techniques, FA, and IHC (acquired with two-photon/confocal microscopy). The lateral resolution differed between techniques: 5.75 μm/pixel for OCTA (Rocholz et al., 2018), 5.6 μm/pixel for FA and 0.541 μm/pixel for confocal microscopy. The axial resolution of OCTA and confocal microscopy images was 3.87 μm/pixel and 1.5 μm/pixel respectively (Rocholz et al., 2018). All three techniques had sufficient resolution to enable detection of the smallest capillaries, which in our samples had a diameter of approximately 8 μm, similar to that reported in human retinas (Tan et al., 2012). Limitations inherent to each technique explained the differences in image quality but did not affect the quantitative comparison of vascular parameters measured by OCTA and IHC/confocal microscopy.

In our evaluation, OCTA angiograms were of higher resolution when compared to FA images at the same location. Similar to what is reported in humans (Spaide et al., 2015), OCTA in dogs enabled identification of capillary beds (ICP and DCP) that were not identified with FA. However, FA provided a larger field of view and the artefacts that were found in some of the OCTA scans (motion artefacts and decorrelation abnormalities due to projection artefact) were not seen in FA images.

When compared to the images obtained by IHC on retinal wholemounts, our study confirms that OCTA provides good visualisation of the SVP and the DCP. We also found that there was an underrepresentation of small calibre vessels on OCTA, especially those located in highly reflective layers (ICP). These findings have also been reported when evaluating OCTA in ex vivo pig eyes (Yu et al., 2021). When compared to images acquired at the same locations by confocal microscopy/IHC, our results suggest that OCTA is a valuable technique for visualising and quantifying retinal vasculature in dogs, especially for analysis of VD in the DCP. Additionally, by IHC we found that the ICP merges with the SVP but not with the DCP as it occurs in human retinas (Chan et al., 2012).

There are some limitations to be considered, such as the small number of eyes that were used for quantitative validation of the OCTA by histology analysis. Additionally, there was minimal misalignment between the images obtained by OCTA and IHC, likely due to shrinkage/expansion due to fixation and distortion of the retinal tissue during the flatmounting process.

Conclusion

Our study confirmed the feasibility of using OCTA in dogs, providing depth-resolved images from different segmented retinal layers that allow evaluation of individual plexuses. This paves the way for further in vivo analyses of canine retinal vasculature in a wide number of retinal pathologies with a vascular phenotype.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Ethics Statement

The animal study was reviewed and approved by IACUC of the University of Pennsylvania.

Author Contributions

AR-G conducted all the in vivo imaging in the animals, performed the immunohistochemical examination, conceived the project, and wrote the original draft. GR acquired the confocal microscopy images of immunolabeled cryosections. G-SY and YC performed the statistical analysis. NC conceived the project. GA acquired the funding and provided the animal resources. WB conceived the project, acquired the funding, provided the animal resources, supervised the project, reviewed and edited the original graft. All the authors contributed to the article and approved the submitted version.

Funding

This study was supported by NIH grants U24EY029890, RO1EY017549, RO1EY06855, P30EY001583, S10 OD021633-01, Foundation Fighting Blindness, Fighting Blindness Canada.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors would like to thank Lydia Melnyk for research coordination. Also, Terry Jordan, Jacqueline Wivel, Nanci Newsom, Mayra Quiroz, and the staff of the RDSF for anaesthesia and animal care support. And finally, Tahira Akala and Dr. Alexa P. Gray for help with manuscript editing.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnana.2021.785249/full#supplementary-material

Footnotes

  1. ^ https://ccrod.cancer.gov/confluence/display/ROB2/Downloads

References

Al-Sheikh, M., Ghasemi Falavarjani, K., Akil, H., and Sadda, S. R. (2017). Impact of image quality on OCT angiography based quantitative measurements. Int. J. Retina Vitreous 3:13. doi: 10.1186/s40942-017-0068-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrett, J. W. (1886). The distribution of blood vessels in the retinae of mammals. J. Physiol. 7, 230–239. doi: 10.1113/jphysiol.1886.sp000218

PubMed Abstract | CrossRef Full Text | Google Scholar

Beltran, W. A., Cideciyan, A. V., Guziewicz, K. E., Iwabe, S., Swider, M., Scott, E. M., et al. (2014). Canine retina has a primate fovea-like bouquet of cone photoreceptors which is affected by inherited macular degenerations. PLoS One 9:e90390. doi: 10.1371/journal.pone.0090390

PubMed Abstract | CrossRef Full Text | Google Scholar

Braun, R. D., Linsenmeier, R. A., and Goldstick, T. K. (1995). Oxygen consumption in the inner and outer retina of the cat. Invest. Ophthalmol. Vis. Sci. 36, 542–554.

Google Scholar

Campbell, J. P., Zhang, M., Hwang, T. S., Bailey, S. T., Wilson, D. J., Jia, Y., et al. (2017). Detailed vascular anatomy of the human retina by projection-resolved optical coherence tomography angiography. Sci. Rep. 7:42201. doi: 10.1038/srep42201

PubMed Abstract | CrossRef Full Text | Google Scholar

Cavallerano, A. A. (1996). Ophthalmic fluorescein angiography. Optom. Clin. 5, 1–23.

Google Scholar

Chan, G., Balaratnasingam, C., Yu, P. K., Morgan, W. H., Mcallister, I. L., Cringle, S. J., et al. (2012). Quantitative morphometry of perifoveal capillary networks in the human retina. Invest. Ophthalmol. Vis. Sci. 53, 5502–5514. doi: 10.1167/iovs.12-10265

PubMed Abstract | CrossRef Full Text | Google Scholar

Chatzistergiou, V., Cilliers, H., Pournaras, J. A., and Ambresin, A. (2021). Foveaplana on optical coherence tomography angiography: new perspectives. Retina 41, 1541–1546. doi: 10.1097/IAE.0000000000003046

PubMed Abstract | CrossRef Full Text | Google Scholar

Chui, T. Y., Zhong, Z., Song, H., and Burns, S. A. (2012). Foveal avascular zone and its relationship to foveal pit shape. Optom. Vis. Sci. 89, 602–610. doi: 10.1097/OPX.0b013e3182504227

PubMed Abstract | CrossRef Full Text | Google Scholar

Cicinelli, M. V., Carnevali, A., Rabiolo, A., Querques, L., Zucchiatti, I., Scorcia, V., et al. (2017). Clinical Spectrum of Macular-Foveal Capillaries Evaluated with Optical Coherence Tomography Angiography. Retina 37, 436–443. doi: 10.1097/IAE.0000000000001199

PubMed Abstract | CrossRef Full Text | Google Scholar

Coscas, G., Lupidi, M., and Coscas, F. (2016). “Heidelberg Spectralis optical coherence tomography angiography: technical aspects,” in OCT Angiography in Retinal and Macular Diseases. Developments in Ophthalmology, eds F. Bandello, E. H. Souied, and G. Querques (Basel: Karger), 1–5. doi: 10.1159/000442768

PubMed Abstract | CrossRef Full Text | Google Scholar

Cuenca, N., Ortuno-Lizaran, I., Sanchez-Saez, X., Kutsyr, O., Albertos-Arranz, H., Fernandez-Sanchez, L., et al. (2020). Interpretation of OCT and OCTA images from a histological approach: clinical and experimental implications. Prog. Retin. Eye Res. 77:100828. doi: 10.1016/j.preteyeres.2019.100828

PubMed Abstract | CrossRef Full Text | Google Scholar

Dolz-Marco, R., Phasukkijwatana, N., Sarraf, D., and Freund, K. B. (2016). Optical coherence tomography angiography in fovea plana. Ophthalmic Surg. Lasers Imaging Retina 47, 670–673. doi: 10.3928/23258160-20160707-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Engerman, R. L., Molitor, D. L., and Bloodworth, J. M. Jr. (1966). Vascular system of the dog retina: light and electron microscopic studies. Exp. Eye Res. 5, 296–301. doi: 10.1016/s0014-4835(66)80039-8

CrossRef Full Text | Google Scholar

Fischer, M. W., and Slatter, D. H. (1978). Preparation and orientation of canine Retinal vasculature: a modified trypsin digestion technique. Aust. J. Ophthalmol. 6, 46–50. doi: 10.1111/j.1442-9071.1978.tb00252.x

CrossRef Full Text | Google Scholar

Flower, R. W., Mcleod, D. S., Lutty, G. A., Goldberg, B., and Wajer, S. D. (1985). Postnatal retinal vascular development of the puppy. Invest. Ophthalmol. Vis. Sci. 26, 957–968.

Google Scholar

Gelatt, K. N., Henderson, J. D. Jr., and Steffen, G. R. (1976). Fluorescein angiography of the normal and diseased ocular fundi of the laboratory dog. J. Am. Vet. Med. Assoc. 169, 980–984.

Google Scholar

Giannakaki-Zimmermann, H., Kokona, D., Wolf, S., Ebneter, A., and Zinkernagel, M. S. (2016). Optical Coherence Tomography Angiography in Mice: comparison with Confocal Scanning Laser Microscopy and Fluorescein Angiography. Transl. Vis. Sci. Technol. 5:11. doi: 10.1167/tvst.5.4.11

PubMed Abstract | CrossRef Full Text | Google Scholar

Harris, A., Ciulla, T. A., Chung, H. S., and Martin, B. (1998). Regulation of retinal and optic nerve blood flow. Arch. Ophthalmol. 116, 1491–1495. doi: 10.1001/archopht.116.11.1491

PubMed Abstract | CrossRef Full Text | Google Scholar

Jia, Y., Tan, O., Tokayer, J., Potsaid, B., Wang, Y., Liu, J. J., et al. (2012). Split-spectrum amplitude-decorrelation angiography with optical coherence tomography. Opt. Express 20, 4710–4725. doi: 10.1364/oe.20.004710

PubMed Abstract | CrossRef Full Text | Google Scholar

Kwiterovich, K. A., Maguire, M. G., Murphy, R. P., Schachat, A. P., Bressler, N. M., Bressler, S. B., et al. (1991). Frequency of Adverse Systemic Reactions after Fluorescein Angiography. Ophthalmology 98, 1139–1142. doi: 10.1016/s0161-6420(91)32165-1

CrossRef Full Text | Google Scholar

Leber, T. (1875). Die Circulations-und Ernährungsverhältnisse des Auges. Canada: AbeBooks Inc

Google Scholar

Liu, W., Luisi, J., Liu, H., Motamedi, M., and Zhang, W. (2017). OCT-Angiography for non-Invasive monitoring of neuronal and Vascular structure in mouse retina: implication for characterization of retinal neurovascular coupling. EC Ophthalmol. 5, 89–98.

Google Scholar

López-Sáez, M. P., Ordoqui, E., Tornero, P., Baeza, A., Sainza, T., Zubeldia, J. M., et al. (1998). Fluorescein-induced allergic reaction. Ann. Allergy Asthma Immunol. 81, 428–430. doi: 10.1016/s1081-1206(10)63140-7

CrossRef Full Text | Google Scholar

Marmor, M. F., Choi, S. S., Zawadzki, R. J., and Werner, J. S. (2008). Visual insignificance of the foveal pit: reassessment of foveal hypoplasia as fovea plana. Arch. Ophthalmol. 126, 907–913. doi: 10.1001/archopht.126.7.907

PubMed Abstract | CrossRef Full Text | Google Scholar

Michaelson, I. C. (1954). Retinal Circulation in Man and Animals. Springfield, Illinois: Charles C Thomas Publisher.

Google Scholar

Munk, M. R., Giannakaki-Zimmermann, H., Berger, L., Huf, W., Ebneter, A., Wolf, S., et al. (2017). OCT-angiography: a qualitative and quantitative comparison of 4 OCT-A devices. PLoS One 12:e0177059. doi: 10.1371/journal.pone.0177059

PubMed Abstract | CrossRef Full Text | Google Scholar

Mutlu, F., and Leopold, I. H. (1964). Structure of the retinal vascular system of the dog, monkey, rat, mouse and cow. Am. J. Ophthalmol. 58, 261–270. doi: 10.1016/0002-9394(64)91575-2

CrossRef Full Text | Google Scholar

Park, M. M., Young, B. K., Shen, L. L., Adelman, R. A., and Del Priore, L. V. (2021). Topographic Variation of Retinal Vascular Density in Normal Eyes Using Optical Coherence Tomography Angiography. Transl. Vis. Sci. Technol. 10, 15–15. doi: 10.1167/tvst.10.12.15

PubMed Abstract | CrossRef Full Text | Google Scholar

Parry, H. B. (1953). Degenerations of the dog retina. I. Structure and development of the retina of the normal dog. Br. J. Ophthalmol. 37, 385–404. doi: 10.1136/bjo.37.7.385

PubMed Abstract | CrossRef Full Text | Google Scholar

Pechauer, A. D., Jia, Y., Liu, L., Gao, S. S., Jiang, C., and Huang, D. (2015). Optical coherence tomography angiography of peripapillary retinal blood flow response to hyperoxia. Invest. Ophthalmol. Vis. Sci. 56, 3287–3291. doi: 10.1167/iovs.15-16655

PubMed Abstract | CrossRef Full Text | Google Scholar

Peichl, L. (1992). Topography of ganglion cells in the dog and wolf retina. J. Comp. Neurol. 324, 603–620. doi: 10.1002/cne.903240412

PubMed Abstract | CrossRef Full Text | Google Scholar

Rocholz, R., Teussink, M., Dolz-Marco, R., Holzhey, C., Dechent, J., Tafreshi, A., et al. (2018). SPECTRALIS optical coherence tomography angiography (OCTA): principles and clinical applications. Germany: Heidelberg Engineering Academy, 1–10.

Google Scholar

Rochon-Duvigneaud, A. (1943). Les yeux et la vision des vertébrés. Paris: Masson.

Google Scholar

Schaepdrijver, L. D., Simoens, P., and Lauwers, H. (1996). Fluorescein angiography of the canine retina. Vet. Comp. Ophthalmol. 6, 111–119.

Google Scholar

Schaepdrijver, L. D., Simoens, P., Lauwers, H., and Geest, J. P. D. (1989). Retinal vascular patterns in domestic animals. Res. Vet. Sci. 47, 34–42. doi: 10.1016/s0034-5288(18)31228-1

CrossRef Full Text | Google Scholar

Spaide, R. F., Klancnik, J. M. Jr., and Cooney, M. J. (2015). Retinal vascular layers imaged by fluorescein angiography and optical coherence tomography angiography. JAMA Ophthalmol. 133, 45–50. doi: 10.1001/jamaophthalmol.2014.3616

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, M. T., Huang, S., Chan, W., Craig, J. E., Knight, L. S. W., Sanders, P., et al. (2021). Impact of cardiometabolic factors on retinal vasculature: A 3 x 3, 6 x 6 and 8 x 8-mm ocular coherence tomography angiography study. Clin. Exp. Ophthalmol. 49, 260–269. doi: 10.1111/ceo.13913

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, P. E., Yu, P. K., Balaratnasingam, C., Cringle, S. J., Morgan, W. H., Mcallister, I. L., et al. (2012). Quantitative confocal imaging of the retinal microvasculature in the human retina. Invest. Ophthalmol. Vis. Sci. 53, 5728–5736. doi: 10.1167/iovs.12-10017

PubMed Abstract | CrossRef Full Text | Google Scholar

Told, R., Reiter, G. S., Schranz, M., Reumueller, A., Hacker, V., Mittermueller, T. J., et al. (2021). Correlation of retinal thickness and swept-source optical coherence tomography angiography derived vascular changes in patients with neovascular age-related macular degeneration. Curr. Eye Res. 46, 1002–1009. doi: 10.1080/02713683.2020.1849734

PubMed Abstract | CrossRef Full Text | Google Scholar

Yao, X., Alam, M. N., Le, D., and Toslak, D. (2020). Quantitative optical coherence tomography angiography: a review. Exp. Biol. Med. 245, 301–312. doi: 10.1177/1535370219899893

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, D. Y., and Cringle, S. J. (2001). Oxygen distribution and consumption within the retina in vascularised and avascular retinas and in animal models of retinal disease. Prog. Retin. Eye Res. 20, 175–208. doi: 10.1016/s1350-9462(00)00027-6

CrossRef Full Text | Google Scholar

Yu, J. J., Camino, A., Liu, L., Zhang, X., Wang, J., Gao, S. S., et al. (2019). Signal strength reduction effects in OCT angiography. Ophthalmol. Retina 3, 835–842. doi: 10.1016/j.oret.2019.04.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, P. K., Mehnert, A., Athwal, A., Sarunic, M. V., and Yu, D. Y. (2021). Use of the retinal vascular histology to validate an optical coherence tomography angiography technique. Transl. Vis. Sci. Technol. 10:29. doi: 10.1167/tvst.10.1.29

PubMed Abstract | CrossRef Full Text | Google Scholar

Zudaire, E., Gambardella, L., Kurcz, C., and Vermeren, S. (2011). A computational tool for quantitative analysis of vascular networks. PLoS One 6:e27385. doi: 10.1371/journal.pone.0027385

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: OCTA, canine retinal vasculature, vasculature quantification, vessel density (VD), vascular plexuses

Citation: Ripolles-Garcia A, Ruthel G, Ying G-S, Chen Y, Cuenca N, Aguirre GD and Beltran WA (2021) Characterization of the Canine Retinal Vasculature With Optical Coherence Tomography Angiography: Comparisons With Histology and Fluorescein Angiography. Front. Neuroanat. 15:785249. doi: 10.3389/fnana.2021.785249

Received: 28 September 2021; Accepted: 23 November 2021;
Published: 13 December 2021.

Edited by:

Enrica Strettoi, Institute of Neuroscience, Consiglio Nazionale delle Ricerche (CNR), Italy

Reviewed by:

Miruna G. Ghinia-Tegla, City College of New York (CUNY), United States
Michael Kalloniatis, University of New South Wales, Australia

Copyright © 2021 Ripolles-Garcia, Ruthel, Ying, Chen, Cuenca, Aguirre and Beltran. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: William A. Beltran, wbeltran@vet.upenn.edu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.