Identification of Key Genes and Pathways in Persistent Hyperplastic Primary Vitreous of the Eye Using Bioinformatic Analysis

Abstract

Background: The failure of the embryonic hyaloid vascular system to regress naturally causes persistent hyperplastic primary vitreous (PHPV), a congenital eye disease. PHPVs molecular pathway, candidate genes, and drug targets are unknown. The current paper describes a comprehensive analysis using bioinformatics to identify the key genes and molecular pathways associated with PHPV, and to evaluate potential therapeutic agents for disease management.

Methods: The genes associated with PHPV were identified using the pubmed2ensembl text mining platform. GeneCodis was employed to evaluate the Gene Ontology (GO) biological process terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Search Tool for the Retrieval of Interacting Genes (STRING) constructed a protein-protein interaction (PPI) network from the text mining genes (TMGs) in Cytoscape. The significant modules were clustered using Molecular Complex Detection (MCODE), and the GO and KEGG analysis for the hub genes were analyzed with the Database of Annotation, Visualization and Integrated Discovery (DAVID) tool. ClueGO, CluePedia, and ShinyGo were used to illustrate the functions and pathways of the clustered hub genes in a significant module. The Drug-Gene Interaction database (DGIdb) was used to evaluate drug–gene interactions of the hub genes to identify potential PHPV drug candidates.

Results: A total of 50 genes associated with PHPV were identified. Overall, 35 enriched GO terms and 15 KEGG pathways were discovered by the gene functional enrichment analysis. Two gene modules were obtained from the PPI network constructed with 31 nodes with 42 edges using MCODE. We selected 14 hub genes as core candidate genes: TP53, VEGFA, SMAD2, CDKN2A, FOXC, FZD4, LRP5, KDR, FZD5, PAX6, MYCN, NDP, PITX2, and PAX2, primarily associated with camera-type eye morphogenesis, pancreatic cancer, the apoptotic process involved in morphogenesis, and the VEGF receptor signaling pathway. We discovered that 26 Food and Drug Administration (FDA)-approved drugs could target 7 of the 14 hub genes.

Conclusions: In conclusion, the results revealed a total of 14 potential genes, 4 major pathways, 7 drug gene targets, and 26 candidate drugs that could provide the basis of novel targeted therapies for targeted treatment and management of PHPV.

Schematic diagram for evaluating and identifying the key genes, molecular pathways, predicted drug targets of PHPV.

Introduction

Persistent hyperplastic primary vitreous (PHPV) is a rare vitreoretinal disorder which accounts for up to 5% of blindness (1). PHPV's pathophysiology occurs during the embryonic stage, with vessel development occurring in the third week of pregnancy. The hyaloid artery system expands and extends to the anterior part of the eye forming the iridohyaloid or capsulopupillary artery during this period. At this point of development, the posterior tunica vasculosa lentis, which is an anastomosis of vessels at the back of the lens begins to develop and nourishes the lens. Secondary vitreous begins to develop in place of primary vitreous during the second trimester of pregnancy. The pathological persistence of fetal intraocular vessels including the hyaloid artery in embryonic vitreous causes this congenital eye disease (2, 3). Apoptosis or macrophage activation causes hyaloidal artery regression which is accompanied by vasa hyaloidal propria, iridohyaloid, and tunica vasculosa lentis (4). White retrolental tissue, an anteriorly swollen lens, centrally dragging ciliary structures, and varying degrees of lenticular opacification are the most prominent clinical symptoms (5). These vascular remnants can hinder the normal retinal development leading to retinal detachment and optic nerve or macula anomalies, and it can also appear in anterior, lateral, or combined forms in various patients (5, 6). PHPV is usually detected in infants within the first 3 months of life due to leukocoria, microphthalmos, and strabismus (7). Figure 1 represents the typical morphology of a PHPV subject. PHPV is also known as persistent fetal vasculature (PFV) (5). Bilateral PHPV is rare and sporadic compared to unilateral PHPV; however, it is an autosomal dominant or recessive trait that may be inherited (8, 9).

Figure 1

A large number of genes are involved in the development and regression of the hyaloid artery. PHPV traits have been observed in human and animal models. In humans, PHPV incidence was found to be an autosomal dominant inheritance pattern in an Egyptian family (10). Mutations in the NDP gene and the COX15 gene on chromosome 10 have been found in cases of bilateral PHPV (11–14). The ZNF408 gene, which had previously been found in retinitis pigmentosa and autosomal dominant familial exudative vitreoretinopathy (ADFEVR) was also identified in PHPV cases of microcornea, posterior megalolenticonus, and coloboma syndrome (MPPC syndrome) (15). FZD4 (frizzled-type receptor 4) was reported to be associated in certain PHPV cases as well as a gene linked to familial exudative vitreoretinopathy (FEVR) (16). In animal models, various signaling pathways have been implicated in the pathogenesis of PHPV including protooncogene ski, p53, tumor suppressor gene Arf, ephrin-B2, βA3/A1-crystallin, LRP5, ang-2, Bax and Bak, FZD4, and ephrin-A5. FEVR, incontinentia pigmenti, retinoblastoma, and retinopathy of prematurity (RoP) are some of the conditions that mimic PHPV-like symptoms (17–23). However, the regulatory mechanisms responsible and genes involved in the process of fetal vascular regression continue to be unclear, as does the underlying cause of failure of regression.

The current surgical management of PHPV is primarily based on the pathological presence of individual cases. Depending on the ocular pathology of PHPV, the limbal and pars plicata incisions are the two most frequent surgical incision methods (24). The most common criteria for surgical intervention are severe media opacities due to cataract or retrolental membranes, progressive anterior chamber shallowing due to cataract, uncontrolled glaucoma or secondary ocular hypotony related to ciliary process dragging, vitreous hemorrhage, and retinal detachment following vitreoretinal traction (5, 25). In cases with advanced pathology, such as acute optic nerve hypoplasia, severe retinal detachment, or microphthalmia, surgery is not a preferred choice since post-operative vision is often low (24). Non-surgical management is currently used in non-progressive conditions and patients with non-central opacity that does not cause any visual impairment. If a non-surgical alternative is used, diligent follow-up should be carried out to detect any potential risks, such as cataract progression or glaucoma (26). The disease's heterogeneity continues to render PHPV diagnosis and treatment challenging.

Since PHPV is a rare disease, understanding the mechanisms that constitute a group of phenotypes is often restricted by small sampling sizes. Therefore, comprehending the molecular mechanisms underlying the expression of the mutated gene, which leads to improper vascular remodeling and the formation of PHPV is often individual-specific and critical for diagnosis, prevention, and therapeutic management. The assessment and analysis of molecular pathways and genetic variant analysis using conventional variant detection approaches such as Sanger sequencing, next generation sequencing, FISH, aCGH, and GTG banding can be time-consuming, expensive, and results in complicated data analysis for unspecified variants (27–31). Text mining is an effective tool for generating a hypothesis since it can reveal novel correlations between genes and the disease pathologies (32). Integration of text mining with biological knowledge and a bioinformatic approach provides new insights into the potential to reconfigure existing drugs (33). By integrating biological databases and in silico tools, the present paper aims to explore possible molecular mechanisms (if any) and classify the causative genes responsible for the heterogenic disease PHPV, thus discovering new drug targets for the treatment of the disease.

Methods

Selection of Key Genes Using Text Mining Analysis

To identify genes related to PHPV, text mining analysis was performed using pubmed2ensembl (http://pubmed2ensembl.ls.manchester.ac.uk) which revealed associations between genes and the literature for data extraction. It is a freely accessible database that connects over 2,000,000 articles in PubMed publications to 150,000 Ensembl genes from 50 species (34, 35). To create a list of key genes, we used search terminology “Persistent Hyperplastic Primary Vitreous” and “Persistent Fetal Vasculature” from 100,000 relevant document IDs. The search terms used were confined to avoid overlapping genes related with other ocular disorders. The species dataset was set to “Homo sapiens (GRCh37)” and the query result was constrained using “filter on MEDLINE: PubMed ID”. The unduplicated genes were extracted and the TMGs were recovered as the intersection of gene hits from the two sets. Figure 2 represents the methodology flowchart and summary of the study design.

Figure 2

Pathway Enrichment and Biological Process Analysis

The TMGs obtained from text mining were analyzed for biological process annotations. The tool GeneCodis, a web-based server, was used (http://genecodis.cnb.csic.es/) to execute an enrichment analysis of the TMGs. GeneCodis assesses functional analysis of gene lists that integrates different sources of information which includes Gene Ontology (GO) [a collection of terminology that describe gene products in terms of Biological Process (BP), Molecular Function (MF), and Cellular Component (CC)], KEGG pathways (offers evidence on biological metabolic pathways that are well-known), and Inter Pro motifs (36). The organism chosen for the analysis was set as Homo sapiens. The TMGs were used as the input set, and genes with significantly enriched biological processes relevant to eye development and vasculogenesis were chosen using an adjusted P-value and analyzed using the GO and BP categories. Using the same method, the genes from the selected annotations were used for KEGG pathway analysis and the genes obtained by the KEGG pathway analysis were further analyzed (28). GeneMania (version 3.5.2), a Cytoscape plugin (version 3.8.2), was used to construct a gene-gene functional interaction network from the TMGs. The advanced statistical options used were max resultant genes = 20, max resultant attributes = 10, and the automatically selected network weighting function. The resulting network comprised functional annotations from GO as well as genes most closely related to the original list.

Construction of Protein-Protein Interaction Network and Module Analysis

STRING (version 1.6.0) was used to construct the PPI network of 35 enriched genes based on GO. STRING is a web-based database comprising nearly 24.6 million proteins and over 3.1 billion interactions from 5,090 distinct species [https://string-db.org/cgi/input.pl; (37)]. The fundamental metrics of nodes in network theory are connectivity degree (k), Betweenness Centrality (BC), Closeness Centrality (CC), Eigenvector Centrality (EC), and eccentricity. However, the main advantage of PPI network analysis is to accommodate a wide range of biological processes including inputs pathway information, providing confidence scores based on evidence from conserved genomic neighborhoods, gene-fusion events, co-occurrence events, co-expression data, experimental data, database information, text mining, and homology. In the PPI network, nodes with a high degree known, as hub proteins, are critical proteins because they may correlate to disease-causing genes while nodes with a high BC, known as bottlenecks, prefer to signify important genes because they can be compared to highly used intersections on major highways or bridges. The confidence score of 0.900 was specified as the minimum criterion. The molecular interaction network was then visualized and hub genes were identified using the Cytoscape software which visually presents the integration of gene expression, biological network, and genotype (38). In this study, the hub nodes were classified by a high score based on the network's scale-free property and was used for centrality analysis by analyzing the network topology (39) and considered the sub-network of these key proteins as the backbone which was worth exploring further in the signaling pathways involved in eye development. Further, a built in Cytoscape plugin Molecular Complex Detection (MCODE, version 2.0.0) was used to distinguish the significant gene modules (clusters) and hub genes from the PPI network (40). The cutoff parameters were “degree cutoff = 2,” “node score cutoff = 0.2,” “k-core = 2,” and “max depth = 100” (41).

Drug-Gene Interactions

The Drug-Gene Interaction Database (DGIdb) (www.dgidb.org) is an online resource that consolidates data from various sources to illustrate drug–gene interactions and gene druggability (42). We investigated drug-gene interactions used in significant module genes as the potential targets for existing drugs or compounds using DGIdb (Version 3.0). The PubChem database was used to obtain the chemical structure of the identified drugs (https://pubchem.ncbi.nlm.nih.gov). It has over 25 million specific chemical structures and 90 million bioactivity outcomes linked to thousands of macromolecular targets.

Results

Identification of Candidate Genes

We obtained 50 unique genes in Homo sapiens associated with PHPV using the TMG approach. Figure 3 depicts the network, genetic interactions, co-expression analysis, and pathways of the 50 TMGs assessed by GeneMania. From these, 35 genes were selected as candidate genes for enrichment analysis based on their GO and molecular pathways.

Figure 3

Enrichment Analysis of TMGs

The most enriched terminology directly linked to the pathology of vasculature morphogenesis of the eye contributing to PHPV was identified using GeneCodis (with P = 1.00E-07), GO, biological process (BP), and KEGG. The GO and BP annotations analysis identified 35 significantly enriched genes. The 10 most enriched functions were “extracellular matrix cell signaling” (P = 2.33E-08), “multicellular organism development” (P = 2.51E-08), “regulation of transcription by RNA polymerase ii” (P = 3.12 E-08), “Norrin signaling pathway” (P = 3.50E-08), “regulation of transcription, DNA templated” (P = 3.77E-08), “negative regulation of gene expression” (P = 1.27E-07), “retina vasculature morphogenesis in camera-type eye” (P = 1.53E-07), “heart development” (P = 2.04E-06), “vascular endothelial growth factor (VEGF) signaling pathway” (P = 2.17E-06), and “post-embryonic camera type eye development” (P = 4.52E-05). Overall, 15 major pathways involving 12 TMGs were discovered by KEGG enrichment analysis. The five most significantly enriched pathways were “pathways in cancer” (P = 11.24E-12), “proteoglycans in cancer” (P = 5.54E-12), “hepatocellular carcinoma” (P = 3.16E-11), “gastric cancer” (P = 3.94E-11), and “pancreatic cancer” (P = 8.69E-10), involving 9, 7, 6, 6, and 5 text mining genes, respectively. Table 1 displays 15 enriched GO terms and Table 2 exhibits the KEGG analysis of 10 enriched molecular pathways of the TMGs.

PPI Network Construction, Modular Analysis, and Key Genes Identification

STRING was used to construct a PPI network for the 35 target genes with a high confidence score >0.900. There were 31 nodes and 42 edges in the network (Figure 4A). Using a cluster analysis of filtering nodes, 14 hub node genes were identified among 31 nodes (Table 3). The hub genes identified were TP53, VEGFA, SMAD2, CDKN2A, FOXC, FZD4, LRP5, KDR, FZD5, PAX6, MYCN, NDP, PITX2, and PAX2. The REVIGO analysis of the hub genes revealed five clusters based on GO similarity which were primarily related to eye development, Wnt signaling pathway, cell proliferation, regulation of cell migration, and regulation of angiogenesis (Figure 5). The modular analysis performed using MCODE yielded two modules. The PPI network relies on a total of 9 genes, as module 1 (FZD4, FZD5, LRP5, and NDP) contained 4 genes with 10 edges and module 2 (TP53, KDR, VEGFA, CDK2NA, and SMAD2) contained 5 genes with 6 edges (Figures 4B,C). According to the pathway enrichment analysis using KEGG and the ShinyGo platform, the genes in module 1 were associated with VEGF signaling pathway, regulation of execution process of apoptosis, and cell migration involved in sprouting angiogenesis. The module 2 genes were significantly associated with eye development, retinal vasculature development, and Wnt signaling pathway (Figure 6). Overall, the enrichment analysis revealed that these genes were substantially enriched in cell proliferation, anatomical structure morphogenesis, and regulation of developmental process which play a crucial role in vasculature formation of the lens causing PHPV.

Figure 4

Figure 5

Figure 6

Drug-Gene Interaction Analysis of Core Genes

In the drug-gene interaction study, we selected 14 hub genes as potential drug targets (Table 4). Overall, 7 of the 14 are potential gene targets and 26 FDA-approved drugs are expected to have drug-gene interactions. FOXC1, FZD4, LRP5, FZD5, PAX6, NDP, and PITX2 were the exceptions. The major interactions among drugs, genes, and pathways are depicted in Table 4. Table 5 represents the chemical structure and formula of the identified drugs.

Table 5

Drug	Gene target	Chemical structure	Chemical formula/protein formula
Pegaptanib sodium	VEGFA		C22H44N3O10P
Minocycline	VEGFA		C23H27N3O7
Carvedilol	VEGFA		C24H26N2O4
Sonidegib	MYCN		C26H26F3N3O3
Amifostine	SMAD2		C5H15N2O3PS
Hydrocortisone	SMAD2		C21H30O5
Bleomycin	SMAD2		C110H168N34O46S7
Tretinoin	SMAD2		C20H28O2
Progesterone	PAX2		C21H30O2
Vandetanib	KDR		C22H24BrFN4O2
Lenvatinib	KDR		21H19ClN4O4
Sunitinib	KDR		C22H27FN4O2
Sorafenib	KDR		C21H16ClF3N4O3
Axinib	KDR		C22H18N4OS
Regorafenib	KDR		C21H15ClF4NO3
Sorafenib	KDR		C21H16ClF3N4O3
Abemaciclib	CDKN2A		C27H32F2N8
Palbociclib	CDKN2A		C24H29N7O2
Alectinib	CDKN2A		C30H34N4O2
Zinc chloride	TP53		Cl2Zn
Trifluridine	TP53		C10H11F3N2O5
Bortezomib	TP53		C19H25BN4O4
Ranibizumab	VEGFA	Monoclonal antibody	C2158H3282N562O681S12
Aflibercept	VEGFA	Monoclonal antibody	C4318H6788N1164S32
Dinutuximab	MYCN	Monoclonal antibody	C6422H9982N1722O2008S48
Ramucirumab	KDR	Monoclonal antibody	C6374H9864N1692O1996S46

Chemical structure of the potential drugs that target the seven candidate genes.

Discussion

PHPV is a disease that leads to blindness or severe vision loss, although there are currently few therapeutic choices (24, 25). On the other hand, PHPV patients are more likely to develop cataracts and closed-angle glaucoma early on in life. Terminal glaucoma, uveitis, retinal detachment, and intra-ocular hemorrhage can be inevitable for these patients (35). As a consequence, the molecular mechanisms that contribute to PHPV must be established. Our analysis discloses that the molecular mechanism of PHPV overlaps with various other signaling pathways contributing to a broader range of therapeutic targets and prognostic biomarkers. The present paper reports 35 genes that might be involved in the development of the eye's vasculature process in the PHPV condition. The enriched GO and BP terms assigned to these genes were associated mainly with extracellular matrix-cell signaling, multicellular organism development, regulation of transcription by RNA polymerase ii, Norrin signaling pathway, and retina vasculature morphogenesis in camera-type eyes. The PPI network and enrichment analysis identified 14 hub genes, TP53, VEGFA, SMAD2, CDKN2A, FOXC, FZD4, LRP5, KDR, FZD5, PAX6, MYCN, NDP, PITX2, and PAX2 that were involved in camera-type eye morphogenesis, pancreatic cancer, the apoptotic process involved in morphogenesis, and the VEGF receptor signaling pathway (Figure 7). The functional analysis and pathways of the key genes in module 1 and module 2 illustrated using ClueGO are displayed in Figure 7A. Figure 7B displays the distribution of functions and pathways among core genes, while Figure 7C reveals KEGG pathways and enriched GO terms, with colors allocated to each pathway.

Figure 7

Based on the evaluations, four genes such as FZD4, LRP5, FZD5, NDP were involved in the process of eye development, retina vasculature development, retinal blood vessel morphogenesis, and the Wnt signaling pathway [Figure 6B; (43, 44)]. The architecture of the retinal vasculature is dependent on highly organized signaling between various cell types of retina, combining internal metabolic conditions with external influences such as oxygen and nutrient supply. In various organs, including the eye, the Wnt signaling pathway is essential for vascular morphogenesis. During eye development, and in vascular eye disorders, Wnt ligands and receptors are key regulators of ocular angiogenesis and also control the development of structured layers of vasculature in retinas as well as the regression of hyaloid vessels (45). FEVR (an inherited disease in which the peripheral retina is hypovascularized to varying degrees) has been attributed to mutations in Wnt pathway components FZD4, LRP5, and the secreted cysteine-knot protein Norrin (46, 47). Norrin is a non-Wnt ligand with a high affinity FZD4 receptor located in the retina and activates the Wnt/β-catenin pathway. Norrie disease, retinopathy of prematurity, and Coats disease are vascular retinopathies caused by defects in the Norrin gene (48). In humans, mutations in NDP and FZD4 have been identified in a limited number of unilateral and bilateral PHPV patients (14, 49–51). ATOH7 mutation (N46H-homozygous) in a family of autosomal recessive PHPV disease traits linked to 10q21 has been identified (52). These variations include deletions, insertions, and missense and nonsense mutations. However, individuals with X-linked FEVR, autosomal dominant FEVR, retinopathy of prematurity, and Norrie disease have also been reported to have mutations in NDP and FZD4 genes (53). According to the GO analysis, five genes TP53, VEGFA, SMAD2, CDKN2A, and KDR (Figure 6A) are involved in the process of regulation of cell migration by the VEGF signaling pathway, angiogenesis, regulation of muscle cell apoptotic process, and embryonic organ development process. Apoptosis is another crucial process in eye development involving extensive programmed cell death associated with morphogenesis (54).

Previous research on transgenic mice models supports our in silico analysis of PHPV to validate the function of these hub genes in hyaloid vasculature regression such as knockouts of the Arf tumor suppressor gene (23, 55, 56), p53 (21, 57), and Frizzled-5 (57) which were associated with PHPV-like phenotypes in mouse models. During mouse eye development, the arf tumor suppressor gene promoted hyaloid vasculature regression and its deficiency may cause a retrolental membrane with persistent hyaloid vessels (9, 23). In Atoh7 knockout mice, hyaloid vessels persist in the vitreous and proliferate to supply the retina which lacks intrinsic vasculature (58, 59). The ephirin-A5 family of receptor tyrosine has been demonstrated to be significant in the regression of the primary vitreous in mouse models (60). Furthermore, in mice lacking LRP5, a Wnt receptor displayed hyaloid vasculature that lasted throughout their lives (61, 62). Given the correlation between transgenic mouse PHPV phenotypes and the hub genes in human congenital defects affecting the eye morphogenesis or retinal vasculature and molecular signaling pathways in module 1 and module 2, it suggests that the pathogenesis of PHPV is regulated by genes in modules 1 and 2.

Twenty-six drugs identified by the drug-gene interaction analysis were classified as anti-neoplastic agents, ocular vascular disorder agents, kinase inhibitors, immune system functioning agents, or corticosteroids. Among them, four potential drugs such as Ranibizumab, Dinutuximab, Pegaptanib Sodium, and Sonidegib were identified based on their high drug-gene interaction score (Table 4). Ranibizumab is a recombinant humanized monoclonal antibody fragment that binds to human vascular endothelial growth factor A (VEGF-A) and thereby prevents it from binding to its receptor and blocking the development of new blood vessels (63). Pegaptanib Sodium is an anti-angiogenic drug used to treat neovascular diseases. It specifically binds to the 165 isoform of VEGF, a protein that is involved in angiogenesis and increased blood vessel leakage (64). Ranibizumab and Pegaptanib Sodium are typically used to treat wet age-related macular degeneration, a type of eye disease (61, 62, 65, 66). They are also used to treat macular edema after retinal vein occlusion, diabetic macular edema, and diabetic retinopathy. Dinutuximab is a GD2-binding human/mouse chimeric monoclonal antibody. It has been proven that the action of pro and anti-angiogenic factors regulates angiogenesis in the development of new capillaries from a pre-existing capillary network (67). Dinutuximab binds to GD2 on the cell surface and induces GD2 expressing cells to lyse by antibody-dependent cell-mediated cytotoxicity and complement-dependent cytotoxicity (68, 69). Sonidegib is an anticancer drug that inhibits the hedgehog (Hh) pathway which is involved in cell differentiation, tissue polarity, and stem cell maintenance during embryonic growth. Hh is essential for the development of the hyaloid loop on the lens's ventral surface by promoting VEGF-mediated angiogenesis. In a zebrafish model, the loss of Hh signaling induced excess sprouting of blood vessels in the dorsal eye and impaired the growth of blood vessels in the ventral eye (70). Regulation of the Hh signaling pathway has been associated with the growth and progression of cancers such as basal cell carcinoma, medulloblastoma (71), and periocular basal cell carcinoma (72).

In PHPV, the ocular fetal vasculature does not go through normal developmental regression. The reasons could be due to presumed loss of apoptosis in PHPV; these natural apoptotic pathways could be pathologically disrupted (22, 73). Apoptosis is a process of cell death that is regulated by a number of gene families. In mice models, the macrophage has been established as a key mediator in studies examining the mechanisms of regression (74, 75). In silico drug-gene analysis using the hub genes of PHPV revealed high interaction with anticancer compounds in the present study. It is understood that vascular quiescence can be regulated by a combination of pro and anti-angiogenic factors. Previous reports have demonstrated that the equilibrium of angiogenic factors such as VEGF and placental growth factor is crucial in vascular regression and mice lacking angiopoietin 2 which regulates angiogenesis by binding to the Tie2 receptor, maintaining fetal vessels in the eyes (18, 76). In addition, the identified drugs can be used for pharmacological screening in mice and zebrafish models to identify compounds affecting vasculature development that could be of therapeutic importance. The results of the study could lead to a better understanding of the potential molecular pathway and possible hyaloid vasculature mechanism, as well as the development of novel therapeutics to prevent or cure this blinding disease, PHPV. Since the current paper focuses on the appropriate path for understanding molecular pathways and therapeutic options for PHPV through in silico analysis, further experimental analysis using animal models is highly recommended to confirm the significance of the candidate genes and pro- and antiangiogenic factors in hyaloid vasculature development and physiology. This continues to be a limitation of the study.

Conclusion and Future Perspectives

To conclude, no specific candidate genes, molecular pathways, or drug targets have been associated with PHPV until now. TP53, VEGFA, SMAD2, CDKN2A, FOXC, FZD4, LRP5, KDR, FZD5, PAX6, MYCN, NDP, PITX2, and PAX2 were identified for the first time as hub genes using in silico tools that may be involved in the development of retinal vasculature and dysfunction of these genes, leading to PHPV. These genes appear to be predominantly associated with functions related to eye morphogenesis, cancer, apoptosis, and VEGF receptor signaling pathways. Previous reports in knockout TP53, VEGFA, FZD4, and NDP transgenic mouse models confirmed the failure of regression of hyaloid vessels and abnormalities in the retinal vasculature. In the future, these in silico analyses will be validated by mutation screening of the hub genes in PHPV patients in order to identify pathogenic variants and gene product expressivity. Clearly, more research is warranted on animals and in human patients as the phenotypic differences will differ from individual to individual based on the expressivity of the gene product. In addition, we identified four genes that may be potential drug targets. Precision medicine for a fetal ocular condition like PHPV presents new challenges but with a possibility. Since PHPV is a rare and often autosomal recessive condition, the present paper is useful when there is little pathological knowledge about the disease or where there is substantial pathway heterogeneity, underlying the clinical phenotype. As a result, a combination of therapeutic methods such as surgical intervention and candidate gene identification may be used not only to analyze biological pathways unique to specific cases, but also to propose potential drug combinations based on gene products annotated to the disease associated with PHPV. This research sheds light on the potential of personalized intervention in the treatment of PHPV indicating a substantial advancement in management strategy.

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.

References

• 1.

  KhaliqSHameedAIsmailMAnwarKLeroyBPayneAMet al. Locus for autosomal recessive nonsyndromic persistent hyperplastic primary vitreous. Investigat Ophthalmol Vis Sci. (2001) 42:2225–8.

  • Pubmed Abstract
  • Google Scholar

• 2.

  ReeseAB. Persistent hyperplastic primary vitreous: the Jackson Memorial Lecture. Am J Ophthalmol. (1955) 40:317–31. 10.1016/0002-9394(55)91866-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3.

  MaqsoodHYounusSFatimaMSaimMQaziS. Bilateral persistent hyperplastic primary vitreous: a case report and review of the literature. Cureus. (2021) 13:e13105. 10.7759/cureus.13105

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4.

  ItoMYoshiokaM. Regression of the hyaloid vessels and pupillary membrane of the mouse. Anatomy Embryol. (1999) 200:403–11. 10.1007/s004290050289

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5.

  GoldbergMF. Persistent fetal vasculature (PFV): an integrated interpretation of signs and symptoms associated with persistent hyperplastic primary vitreous (PHPV) LIV Edward Jackson Memorial Lecture. Am J Ophthalmol. (1997) 124:587–626. 10.1016/S0002-9394(14)70899-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6.

  Müllner-EidenböckAAmonMMoserEKlebermassN. Persistent fetal vasculature and minimal fetal vascular remnants: a frequent cause of unilateral congenital cataracts. Ophthalmology. (2004) 111:906–13. 10.1016/j.ophtha.2003.07.019

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7.

  AntoszykJHAntoszykA. Persistent fetal vasculature 743.51. In: RoyFHFraunfelderFWFraunfelderFT editors. Roy and Fraunfelder's Current Ocular Therapy.Elsevier Health Sciences (2008). p. 660–2.

  • Google Scholar

• 8.

  JainTP. Bilateral persistent hyperplastic primary vitreous. Indian J Ophthalmol. (2009) 57:53. 10.4103/0301-4738.44487

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9.

  ShastryBS. Persistent hyperplastic primary vitreous: congenital malformation of the eye. Clin Exp Ophthalmol. (2009) 37:884–90. 10.1111/j.1442-9071.2009.02150.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10.

  GalalAHKotouryAIAzzabAA. Bilateral persistent hyperplastic primary vitreous: an Egyptian family supporting a rare autosomal dominant inheritance. Genet Counsel. (2006) 17:441–7.

  • Pubmed Abstract
  • Google Scholar

• 11.

  HasbrookMYonekawaYvan LaereLShahARCaponeA. Bilateral persistent fetal vasculature and a chromosome 10 mutation including COX15. Can J Ophthalmol. (2017) 52:e203–5. 10.1016/j.jcjo.2017.04.019

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12.

  PayabvashSAndersonJSNasceneDR. Bilateral persistent fetal vasculature due to a mutation in the Norrie disease protein gene. Neuroradiol J. (2015) 28:623–7. 10.1177/1971400915609350

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13.

  DhingraSShearsDJBlakeVStewartHPatelCK. Advanced bilateral persistent fetal vasculature associated with a novel mutation in the Norrie gene. Br J Ophthalmol. (2006) 90:1324–5. 10.1136/bjo.2005.088625

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14.

  PendergastSDTreseMTLiuXShastryBS. Study of the Norrie disease gene in 2 patients with bilateral persistent hyperplastic primary vitreous. Archiv Ophthalmol. (1998) 116:381–2.

  • Pubmed Abstract
  • Google Scholar

• 15.

  WeinerGANudlemanE. Microcornea, posterior megalolenticonus, persistent fetal vasculature, and coloboma syndrome associated with a new mutation in ZNF408. Ophthal Surg Lasers Imaging Retina. (2019) 50:253–6. 10.3928/23258160-20190401-10

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16.

  RobitailleJMWallaceKZhengBBeisMJSamuelsMHoskin-MottAet al. Phenotypic overlap of familial exudative vitreoretinopathy (FEVR) with persistent fetal vasculature (PFV) caused by FZD4 mutations in two distinct pedigrees. Ophthal Genet. (2009) 30:23–30. 10.1080/13816810802464312

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17.

  McGannonPMiyazakiYGuptaPCTraboulsiEIColmenaresC. Ocular abnormalities in mice lacking the Ski proto-oncogene. Investig Ophthalmol Visual Sci. (2006) 47:4231–7. 10.1167/iovs.05-1543

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18.

  HackettSFWiegandSYancopoulosGCampochiaroPA. Angiopoietin-2 plays an important role in retinal angiogenesis. J Cell Physiol. (2002) 192:182–7. 10.1002/jcp.10128

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19.

  Freeman-AndersonNEZhengYMcCalla-MartinACTreanorLMZhaoYDGarfinPMet al. Expression of the Arf tumor suppressor gene is controlled by Tgfβ2 during development. Development. (2009) 136:2081–9. 10.1242/dev.033548

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20.

  SalvucciOOhnukiHMaricDHouXLiXYoonSOet al. EphrinB2 controls vessel pruning through STAT1-JNK3 signalling. Nat Commun. (2015) 6:6576. 10.1038/ncomms7576

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21.

  ReichelMBAliRRD'EspositoFClarkeARLuthertPJBhattacharyaSSet al. High frequency of persistent hyperplastic primary vitreous and cataracts in p53-deficient mice. Cell Death Different. (1998) 5:156–62. 10.1038/sj.cdd.4400326

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22.

  HahnPLindstenTTolentinoMThompsonCBBennettJDunaiefJL. Persistent fetal ocular vasculature in mice deficient in bax and bak. Archiv Ophthalmol. (2005) 123:797–802. 10.1001/archopht.123.6.797

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23.

  McKellerRNFowlerJLCunninghamJJWarnerNSmeyneRJZindyFet al. The Arf tumor suppressor gene promotes hyaloid vascular regression during mouse eye development. Proc Natl. Acad Sci USA. (2002) 99:3848–53. 10.1073/pnas.052484199

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24.

  DassABTreseMT. Surgical results of persistent hyperplastic primary vitreous. Ophthalmology. (1999) 106:280–4. 10.1016/S0161-6420(99)90066-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25.

  SmithRE. Persistent hyperplastic primary vitreous: results of surgery. Trans Am Acad Ophthalmol Otolaryngol. (1974) 78:911–25.

  • Pubmed Abstract
  • Google Scholar

• 26.

  PrakhunhungsitSBerrocalAM. Diagnostic and management strategies in patients with persistent fetal vasculature: current insights. Clin Ophthalmol. (2020) 14:4325–35. 10.2147/OPTH.S236117

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27.

  FeukLCarsonARSchererSW. Structural variation in the human genome. Nat Rev Genet. (2006) 7:85–97. 10.1038/nrg1767

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28.

  ConradDFHurlesME. The population genetics of structural variation. Nat Genet. (2007) 39:S30–6. 10.1038/ng2042

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29.

  BarkholtLFloryEJekerleVLucas-SamuelSAhnertPBissetLet al. Risk of tumorigenicity in mesenchymal stromal cell–based therapies—bridging scientific observations and regulatory viewpoints. Cytotherapy. (2013) 15:753–9. 10.1016/j.jcyt.2013.03.005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30.

  BishopR. Applications of fluorescence in situ hybridization (FISH) in detecting genetic aberrations of medical significance. Bioscience Horizons. (2010) 3:85–95. 10.1093/biohorizons/hzq009

  • CrossRef
  • Google Scholar

• 31.

  BejjaniBAShafferLG. Application of array-based comparative genomic hybridization to clinical diagnostics. J Mol Diagnost. (2006) 8:528–33. 10.2353/jmoldx.2006.060029

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32.

  MoscaEBertoliGPiscitelliEVilardoLReinboldRAZucchiIet al. Identification of functionally related genes using data mining and data integration: a breast cancer case study. BMC Bioinformatics. (2009) 10(Suppl. 12):S8. 10.1186/1471-2105-10-S12-S8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33.

  ButcherECBergELKunkelEJ. Systems biology in drug discovery. Nat Biotechnol. (2004) 22:1253–9. 10.1038/nbt1017

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34.

  YuSTrancheventL-CDe MoorBMoreauY. Gene prioritization and clustering by multi-view text mining. BMC Bioinformatics. (2010) 11:28. 10.1186/1471-2105-11-28

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35.

  HuDJiangJLinZZhangCMoonasarNQianS. Identification of key genes and pathways in scleral extracellular matrix remodeling in glaucoma: potential therapeutic agents discovered using bioinformatics analysis. Int J Med Sci. (2021) 18:1554–65. 10.7150/ijms.52846

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36.

  Nogales-CadenasRCarmona-SaezPVazquezMVicenteCYangXTiradoFet al. GeneCodis: interpreting gene lists through enrichment analysis and integration of diverse biological information. Nucleic Acids Res. (2009) 37:W317–22. 10.1093/nar/gkp416

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37.

  von MeringCHuynenMJaeggiDSchmidtSBorkPSnelB. STRING: a database of predicted functional associations between proteins. Nucleic Acids Res. (2003) 31:258–61. 10.1093/nar/gkg034

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38.

  SmootMEOnoKRuscheinskiJWangP-LIdekerT. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. (2011) 27:431–2. 10.1093/bioinformatics/btq675

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39.

  TangYLiMWangJPanYWuF-X. CytoNCA: a cytoscape plugin for centrality analysis and evaluation of protein interaction networks. Biosystems. (2015) 127:67–72. 10.1016/j.biosystems.2014.11.005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40.

  WuWFangCZhangCHuNWangL. Identification of 10 important genes with poor prognosis in non-small cell lung cancer through bioinformatical analysis. J Biomed Res Rev. (2020) 3:41–50. 10.21203/rs.3.rs-21044/v1

  • CrossRef
  • Google Scholar

• 41.

  BandettiniWPKellmanPManciniCBookerOJVasuSLeungSWet al. MultiContrast Delayed Enhancement (MCODE) improves detection of subendocardial myocardial infarction by late gadolinium enhancement cardiovascular magnetic resonance: a clinical validation study. J Cardiovasc Magnet Reson. (2012) 14:1–10. 10.1186/1532-429X-14-83

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42.

  CottoKCWagnerAHFengY-YKiwalaSCoffmanACSpiesGet al. DGIdb 3.0: a redesign and expansion of the drug–gene interaction database. Nucleic Acids Res. (2018) 46:D1068–73. 10.1093/nar/gkx1143

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43.

  HuttlinELTingLBrucknerRJGebreabFGygiMPSzpytJet al. The BioPlex network: a systematic exploration of the human interactome. Cell. (2015) 162:425–40. 10.1016/j.cell.2015.06.043

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44.

  ChenT-CLinK-TChenC-HLeeS-ALeeP-YLiuY-Wet al. Using an in situ proximity ligation assay to systematically profile endogenous protein–protein interactions in a pathway network. J Proteome Res. (2014) 13:5339–46. 10.1021/pr5002737

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45.

  WangZLiuC-HHuangSChenJ. Wnt signaling in vascular eye diseases. Progr Retinal Eye Res. (2019) 70:110–33. 10.1016/j.preteyeres.2018.11.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46.

  GongYSleeRBFukaiNRawadiGRoman-RomanSReginatoAMet al. LDL receptor-related protein 5 (LRP5) affects bone accrual and eye development. Cell. (2001) 107:513–23. 10.1016/S0092-8674(01)00571-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47.

  XuQWangYDabdoubASmallwoodPMWilliamsJWoodsCet al. Vascular development in the retina and inner ear: control by Norrin and Frizzled-4, a high-affinity ligand-receptor pair. Cell. (2004) 116:883–95. 10.1016/s0092-8674(04)00216-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48.

  FuhrmannS. Wnt signaling in eye organogenesis. Organogenesis. (2008) 4:60–7. 10.4161/org.4.2.5850

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49.

  WuW-CDrenserKTreseMCaponeADaileyW. Retinal phenotype–genotype correlation of pediatric patients expressing mutations in the Norrie disease gene. Archiv Ophthalmol. (2007) 125:225–30. 10.1001/archopht.125.2.225

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50.

  DhingraSShearsDJBlakeVStewartHPatelCK. Advanced bilateral persistent fetal vasculature. Am J Ophthalmol. (1981) 91:312–7.

  • Pubmed Abstract
  • Google Scholar

• 51.

  AponteEPPulidoJSEllisonJWQuiramPAMohneyBG. A novel NDP mutation in an infant with unilateral persistent fetal vasculature and retinal vasculopathy. Ophthalmic Genet. (2009) 30:99–102. 10.1080/13816810802705755

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52.

  PrasovLMasudTKhaliqSMehdiSQAbidAOliverERet al. ATOH7 mutations cause autosomal recessive persistent hyperplasia of the primary vitreous. Hum Mol Genet. (2012) 21:3681–94. 10.1093/hmg/dds197

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53.

  ShastryBSHiraokaM. Molecular genetics of familial exudative vitreoretinopathy and Norrie disease. Curr Genomics. (2000) 1:259–69. 10.2174/1389202003351445

  • CrossRef
  • Google Scholar

• 54.

  LangRA. Apoptosis in mammalian eye development: lens morphogenesis, vascular regression and immune privilege. Cell Death Different. (1997) 4:12–20. 10.1038/sj.cdd.4400211

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55.

  MartinACThorntonJDLiuJWangXZuoJJablonskiMMet al. Pathogenesis of persistent hyperplastic primary vitreous in mice lacking the arf tumor suppressor gene. Investig Ophthalmol Visual Sci. (2004) 45:3387–96. 10.1167/iovs.04-0349

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56.

  ThorntonJDSwansonDJMaryMNPeiDMartinACPoundsSet al. Persistent hyperplastic primary vitreous due to somatic mosaic deletion of the arf tumor suppressor. Investig Ophthalmol Visual Sci. (2007) 48:491–9. 10.1167/iovs.06-0765

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57.

  ZhangJFuhrmannSVetterML. A nonautonomous role for retinal frizzled-5 in regulating hyaloid vitreous vasculature development. Investig Ophthalmol Visual Sci. (2008) 49:5561–7. 10.1167/iovs.08-2226

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58.

  BrzezinskiJASchulzSMCrawfordSWroblewskiEBrownNLGlaserT. Math5 null mice have abnormal retinal and persistent hyaloid vasculatures. In: Developmental Biology.San Diego, CA; Academic Press Inc.; Elsevier Science (2003). p. 553.

  • Google Scholar

• 59.

  EdwardsMMMcLeodDSLiRGrebeRBhuttoIMuXet al. The deletion of Math5 disrupts retinal blood vessel and glial development in mice. Exp Eye Res. (2012) 96:147–56. 10.1016/j.exer.2011.12.005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60.

  SonAIShelegMCooperMASunYKleimanNJZhouR. Formation of persistent hyperplastic primary vitreous in ephrin-A5^−/− mice. Investig Ophthalmol Visual Sci. (2014) 55:1594–606. 10.1167/iovs.13-12706

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61.

  GaudreaultJFeiDBeyerJCRyanARangellLShiuVet al. Pharmacokinetics and retinal distribution of ranibizumab, a humanized antibody fragment directed against VEGF-A, following intravitreal administration in rabbits. Retina. (2007) 27:1260–6. 10.1097/IAE.0b013e318134eecd

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 62.

  BakriSJSnyderMRReidJMPulidoJSEzzatMKSinghRJ. Pharmacokinetics of intravitreal ranibizumab (Lucentis). Ophthalmology. (2007) 114:2179–82. 10.1016/j.ophtha.2007.09.012

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 63.

  KlettnerARoiderJ. Comparison of bevacizumab, ranibizumab, and pegaptanib in vitro: efficiency and possible additional pathways. Investig Ophthalmol Visual Sci. (2008) 49:4523–7. 10.1167/iovs.08-2055

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64.

  ShuklaDNamperumalsamyPGoldbaumMCunninghamETJr. Pegaptanib sodium for ocular vascular disease. Indian J Ophthalmol. (2007) 55:427. 10.4103/0301-4738.36476

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 65.

  ChongV. Ranibizumab for the treatment of wet AMD: a summary of real-world studies. Eye. (2016) 30:270–86. 10.1038/eye.2015.217

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 66.

  MiyakeMYamashiroKAkagi-KurashigeYKumagaiKNakataINakanishiHet al. Vascular endothelial growth factor gene and the response to anti-vascular endothelial growth factor treatment for choroidal neovascularization in high myopia. Ophthalmology. (2014) 121:225–33. 10.1016/j.ophtha.2013.06.043

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 67.

  HamanoYKalluriR. Tumstatin, the NC1 domain of α3 chain of type IV collagen, is an endogenous inhibitor of pathological angiogenesis and suppresses tumor growth. Biochem Biophys Res Commun. (2005) 333:292–8. 10.1016/j.bbrc.2005.05.130

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 68.

  MujooKChereshDAYangHMReisfeldRA. Disialoganglioside GD2 on human neuroblastoma cells: target antigen for monoclonal antibody-mediated cytolysis and suppression of tumor growth. Cancer Res. (1987) 47:1098–104.

  • Pubmed Abstract
  • Google Scholar

• 69.

  BarkerEReisfeldRA. A mechanism for neutrophil-mediated lysis of human neuroblastoma cells. Cancer Res. (1993) 53:362–7.

  • Pubmed Abstract
  • Google Scholar

• 70.

  WeissOKaufmanRMishaniEInbalA. Ocular vessel patterning in zebrafish is indirectly regulated by Hedgehog signaling. Int J Develop Biol. (2017) 61:277–84. 10.1387/ijdb.160273ai

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 71.

  KoolMJonesDTWJägerNNorthcottPAPughTJHovestadtVet al. Genome sequencing of SHH medulloblastoma predicts genotype-related response to smoothened inhibition. Cancer Cell. (2014) 25:393–405.

  • Pubmed Abstract
  • Google Scholar

• 72.

  JainSSongRXieJ. Sonidegib: mechanism of action, pharmacology, and clinical utility for advanced basal cell carcinomas. Onco Targets Ther. (2017) 10:1645. 10.2147/OTT.S130910

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73.

  PapermasterDSTillyJL. Apoptosis of tissues of the eye during development and disease. In: When Cells Die: A Comprehensive Evaluation of Apoptosis and Programmed Cell Death.New York, NY: Wiley-Liss (1998). p. 321–46.

  • Google Scholar

• 74.

  LangRABishopJM. Macrophages are required for cell death and tissue remodeling in the developing mouse eye. Cell. (1993) 74:453–62. 10.1016/0092-8674(93)80047-I

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 75.

  BalazsEATothLZOzanicsV. Cytological studies on the developing vitreous as related to the hyaloid vessel system. Albrecht von Graefes Archiv für klinische und experimentelle Ophthalmologie. (1980) 213:71–85. 10.1007/BF00413534

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 76.

  FeeneySASimpsonDACGardinerTABoyleCJamisonPStittAW. Role of vascular endothelial growth factor and placental growth factors during retinal vascular development and hyaloid regression. Investig Ophthalmol Visual Sci. (2003) 44:839–47. 10.1167/iovs.02-0040

  • Pubmed Abstract
  • CrossRef
  • Google Scholar