The Genotype-Phenotype Association of Von Hipple Lindau Disease Based on Mutation Locations: A Retrospective Study of 577 Cases in a Chinese Population

Purpose: Von Hippel-Lindau (VHL) disease is a hereditary kidney cancer syndrome, with which patients are more likely to get affected by renal cell carcinoma (RCC), pancreatic cyst or tumor (PCT), central nervous system hemangioblastoma (CHB), retinal angiomas (RA), and pheochromocytoma (PHEO). Mutations of VHL gene located in 3p25 may impair the function of the VHL protein and lead to the disease. It’s unclear why obvious phenotype varieties exist among VHL patients. Here we aimed to ascertain whether the mutation types and locations affect the phenotype.

Methods: We enrolled 577 Chinese VHL patients from 211 families and divided them into three groups and six subgroups according to their mutation types and locations. Cox survival analysis and Kaplan-Meier analysis were used to compare intergroup age-related tumor risks.

Results: Patients with nonsense or frameshift mutations that were located before residues 117 of VHL protein (NoF1 subgroup) hold lower age-related risks of VHL associated tumors (HR = 0.638, 95%CI 0.461–0.883, p = 0.007), CHB (HR = 0.596, 95%CI 0.409–0.868, p = 0.007) or PCT (HR = 0.595, 95%CI 0.368–0.961, p = 0.034) than patients whose mutations were located after residues 117 (NoF2 subgroup). Patients in NoF1 subgroup still had lower age-related risks of CHB (HR = 0.652, 95%CI 0.476–0.893, p = 0.008) and PCT (HR = 0.605, 95%CI 0.398–0.918, p = 0.018) compared with those in combined NoF2 subgroup and other truncating mutation patients. NoF1 subgroup correspondingly had a longer estimated median lifespan (64 vs. 55 year, p = 0.037) than NoF2 subgroup. Among patients with missense mutations of VHL, only a small minority (23 of 286 missense mutations carriers) carried mutations involving neither HIF-α binding region nor elongin C binding region, who were grouped in MO subgroup. MO subgroup seemed to have a higher age-related risk of PHEO. In the whole cohort (n = 577), PHEO was an independent protective factor for CHB (p = 0.001) and survival (p = 0.005). RA and CHB failed to predict the age-related risk of each other.

Conclusion: The mutation types and locations of VHL gene are associated with phenotypes. Genetic counselors could predict phenotypes more accurately based on more detailed genotype-phenotype correlations. Further genotype-phenotype studies should focus on the prediction of tumor recurrence, progression, and metastasis. The deep molecular mechanism of genotype-phenotype correlation is worth further exploring.

Introduction

Von Hippel–Lindau (VHL) disease (OMIM 193300) is a hereditary kidney cancer syndrome that is an autosomal dominant disease caused by mutations of VHL gene, characterized by multiorgan and multicenter tumors, such as renal cell carcinoma (RCC), pancreatic cyst or tumor (PCT), central nervous system hemangioblastoma (CHB), retinal angiomas (RA), pheochromocytoma (PHEO), papillary cystadenoma of the broad ligament or epididymis, and endolymphatic sac tumor (Ong et al., 2007). VHL patients were diagnosed as type 1 when pheochromocytomas were absent and as type 2 when pheochromocytomas were proved (Haddad et al., 2013). CHB, RCC, RA, PCT, and PHEO are the most common symptoms of VHL disease (Nielsen et al., 2016; Carlo et al., 2019; Hong et al., 2019).

The VHL protein (pVHL), coded by the VHL gene located on 3p25, functions as an E3 ubiquitin protein ligase. The protein has 213 amino acids and consists of 2 domains: the α domain and the β domain. The best-characterized and well-known downstream pathway is that pVHL, elongin B and C, Cul2, and Rbx1 compose the pVHL complex that can bind and degrade transcription factor hypoxia-inducible factor α (actually HIF-1α, HIF-2α, HIF-3α; collectively HIF-α), which is an oxygen-dependent pathway. The hypoxia-inducible factor (HIF) transcription factor is composed by two subunits: the oxygen-dependent α subunit (1α, 2α, and 3α) and the constitutive expressing HIF1β subunit. Under the hypoxia condition or in the presence of the mutant VHL gene, the pVHL complex can’t recognize HIF-α so that HIF-α accumulates. Then it heterodimerizes with HIF1β and affects a group of downstream genes’ expression, including erythropoietin (EPO), vascular endothelial growth factor (VEGF), and transforming growth factor (TGF-α), to promote oncogenesis. Most of the pathogenic missense mutations are located in two regions of the VHL protein: the elongin C-binding site (residues 158–184 of the α domain) (Feldman et al., 1999; Ohh et al., 2000) and the HIF-α binding site (residues 65–117 of the β domain) (Hon et al., 2002; Lee et al., 2016).

There have been several studies focusing on genotype-phenotype correlations of VHL disease, providing some valuable tumor risk stratification methods mainly based on mutation types or locations of VHL gene. According to classical theories, type 2 VHL disease was due to missense mutations of VHL and truncating mutations were responsible for type 1 VHL disease (Ong et al., 2007). This theory was to be improved because even among patients with the same mutation types of the VHL gene, there existed phenotype diversities not explained clearly.

Several previous studies reported that significant phenotypic heterogeneity existed among missense mutations carriers (Hong et al., 2019). Mettu et al. (2010) collected phenotypic characters of 412 VHL patients with documented germline missense mutations and concluded that mutations in the α-domain of VHL protein were associated with a higher incidence of RA, compared with mutations in the β-domain. Liu et al. (2018) found that missense mutations located in HIF-binding region of VHL gene were associated with an increased risk of CHB and a decreased risk of PHEO. Forman et al. (2009) held the opinion that carriers of mutations involving the elongin C binding site were more likely to be affected by PEHO through structural analysis, compared with other missense mutations carriers. Only a few studies focused on internal phenotypic heterogeneity in truncating mutations. Ong et al. (2007) found nonsense or frameshift mutations were associated with earlier age of onset of RCC, higher age-related risks of RCC and RA than missense mutations or deletions. Large VHL gene deletions involving C3orf10 (HSPC300) were discovered to be related to a lower lifetime risk of RCC than deletions that do not involve C3orf10 (McNeill et al., 2009).

In this study, we regrouped VHL patients according to their mutation types: missense mutations group (M group), nonsense or frameshift mutations group (NoF group), and other truncating mutations that affect more than one nucleotide of VHL gene (oTR group). M group was further divided into HIF-α binding site mutations subgroup (MH subgroup), elongin C binding site mutations subgroup (ME subgroup), and mutations at other sites subgroup (MO subgroup). NoF group was seen as NoF2 subgroup if the mutations are located after residue 117 of VHL protein, otherwise as NoF1 subgroup. We aim to discover more detailed genotype-phenotype correlations through comparing the difference of age-related penetrance of VHL tumors and survival among subgroups.

Materials and Methods

Medical Ethics

This project had gained the approval of the Institutional Ethics Committee of Peking University First Hospital. Informed consents were signed by all patients or their legal guardians.

Patient Recruitment and Assessment

A total of 577 Chinese VHL patients from 211 unrelated families diagnosed at Peking University First Hospital were recruited for this study. An individual was diagnosed with VHL disease when he or she carried a germline VHL mutation or met the clinical criteria for VHL disease (Shuin et al., 2006). Medical history and previous imaging data of all enrolled patients were collected. The age of onset of five main symptoms, which were CHB, RA, RCC, PCT, PHEO, were reviewed according to medical records, disease history, and imaging data of all VHL patients. Follow-up was carried out from the time genetic or clinically diagnosis of VHL disease for patients until death or January 1, 2019. Patients whose information was unavailable or inaccurate were removed from this cohort.

Genetic Test of VHL Mutations

Peripheral blood was drawn from each enrolled VHL patient unless he or she refused or was dead. DNA was extracted from all collected peripheral blood samples. The three exons of VHL and their flanking intronic sequences were amplified through polymerase chain reaction (PCR) and sequenced by Sanger sequencing. If a suspected patient was negative for PCR and Sanger sequencing, multiplex ligation-dependent probe amplification (MLPA, P016-C2 kit, MRC-Holland, Amsterdam, The Netherlands) followed by real-time quantitative PCR was used to detect large-scale mutations. Patients with point missense mutations were divided into MH subgroup (mutations in residues 65–117), ME subgroup (mutations in residues 158–184), and MO subgroup (mutations in other residues) (Figure 1A). Generally, truncating mutations were those mutations that can truncate or shorten the protein, such as nonsense mutations, frameshift mutations, splice-site mutations, and large-scale mutations. Patients with nonsense or frameshift mutations that were located at or before residues 117 of VHL protein were regrouped in NoF1 group, and others were regrouped in NoF2 subgroup. Other truncating mutations carriers were categorized in oTR group.

FIGURE 1

Figure 1. (A) The structure of the VHL gene. (B) The number of cases and families with missense mutations on each residues of pVHL.

Statistical Analysis

The Mann-Whitney U-test was used to compare the age at onset of each symptom. The survival time from birth to the end of follow-up or to death was used for survival analysis. The age-related risks of the five symptoms of VHL disease and survival analysis were calculated by Kaplan-Meier plot and log-rank analysis. And the Cox regression model was used for multivariate survival analysis to assess the effect of genotype on risk of the five symptoms. The correlation between two symptoms was analyzed by Chi-square (χ²) test. It’s considered to be statistically significant when p-value was less than 0.05. Data were analyzed by SPSS version 23.0 software (IBM, Chicago, IL).

Results

VHL Disease Patients’ Diagnosis and Enrollment

A total of 577 patients from 211 unrelated families diagnosed with the VHL disease in our center were enrolled in our retrospective study. One hundred families (286 patients), 60 families (165 patients), and 52 families (126 patients) were assigned to M group, NoF group, and oTR group, respectively (Table 1). For the 100 families in the M group, they carried a total of 50 mutations involving 42 gene loci or 31 residues of pVHL (Figure 1B). The median age at onset was 30 year (from 2 to 74 years), and the median lifespan of 131 dead patients was 42 years (from 11 to 71 year). Fifty-seven of 577 patients were diagnosed only by genetic tests but had not yet displayed any symptoms. Their age range was from 3 to 60 year, with a median of 15 year. Only five patients with pathogenic VHL mutations over the age of 40 year were asymptomatic. The five patients’ basic information is shown as Supplementary Table 1. Interestingly, some of the five patients’ direct relatives showed symptoms before their 40 year.

TABLE 1

Table 1. The number of families and patients in each subgroup.

Differences in Age of Onset Among Subgroups

As mentioned before, among 520 patients who have shown clinical symptoms, the age at onset of VHL disease in our cohort ranged from 2 to 74 year, with a median of 30 year. CHB was the first symptom of more than half of the patients (311 of 520, 59.8%). Even in each subgroup, CHB was still the most common first symptom. The median age of onset of CHB, RA, RCC, PCT, PHEO is 31 year (3–66 year), 26 year (2–66 year), 37 year (14–74 year), 33 year (14–67 year), 37 year (9–68 year). The number of patients affected by each symptom was shown as Supplementary Table 2. There was no significant distinction in the age of onset between males and females (30 vs. 30.5 year, p = 0.748). NoF1 subgroup had a later age of onset than NoF2 (median age of onset, 28 vs. 32 year, p = 0.018). Significant differences in age of onset were not discovered among subgroup MH, ME, and MO (p = 0.42). There was no significant difference between group NoF and oTR or between group M and those in the combined NoF and oTR together, either (p = 0.983 and p = 0.065, respectively).

Pathologically, RA and CHB were all hemangioblastomas. But in the whole cohort (n = 577) the onset of the two symptoms did not appear to be related (χ² = 0.138, p = 0.711). Patients with CHB didn’t show higher age-related risks of RA (p = 0.494, Supplementary Figure 1A), and patients with RA weren’t associated with higher age-related risks of CHB (p = 0.055, Supplementary Figure 1B).

Comparison of Age-Related Tumor Risk Among Subgroup

Because of the significant difference in age of onset between subgroup NoF1 and NoF2, we compared these two subgroups further. NoF1 subgroup took a lower age-related risk of VHL-associated tumors than subgroup NoF2 (p = 0.004). NoF1 subgroup had lower age-related risks of PCT (p = 0.033) and CHB (p = 0.006) compared with NoF2, but risks of RA (p = 0.717), RCC (p = 0.171), and PHEO (p = 0.711) were similar between the two subgroups (Figure 2). Accordingly, multivariate Cox regression analysis revealed that NoF1 mutations were independently associated with a lower risk of VHL-associated tumors (HR = 0.689, 95%CI 0.497–0.953, p = 0.025), CHB (HR = 0.596, 95%CI 0.409–0.868, p = 0.007), and PCT (HR = 0.595, 95%CI 0.368–0.961, p = 0.034) compared with NoF2 mutations (Table 2).

FIGURE 2

Figure 2. Comparison of age-related risks in NoF1 and NoF2 subgroup. (A) All the five VHL-associated tumors, (B) CHB, (C) RA, (D) RCC, (E) PCT, (F) PHEO. NoF1, nonsense or frameshift mutations before residue 117 of VHL protein; NoF2, nonsense or frameshift mutations after residue 117 of VHL protein; CHB, central nervous system hemangioblastoma; RA, retinal angiomas; RCC, renal cell carcinoma; PCT, pancreatic cyst or tumor; PHEO, pheochromocytoma.

TABLE 2

Table 2. Univariate and multivariate Cox regression analysis of age-related VHL tumor risks of NoF1 vs. NoF2.

Nonsense mutations, frameshift mutations were often regarded as truncating mutations in previous studies (Ong et al., 2007; Forman et al., 2009; Liu et al., 2018). Through comparisons in the age of onset among subgroups, it was proven that there were some differences in truncating mutations. Age-related risks of overall VHL-associated tumors or each VHL tumor were similar between NoF2 subgroup and oTR group. So, we combined patients in NoF2 subgroup and oTR group into one group (N2TR) and compared the group with NoF1 subgroup. NoF1 subgroup was related to lower risks of overall VHL-associated tumors (p = 0.038), CHB (p = 0.006), PCT (p = 0.022) compared with N2TR group. Meanwhile, the risks of RA, RCC, or PHEO didn’t bring out significant differences between the two groups (Supplementary Figure 2). Correspondingly, multivariate Cox regression analysis revealed that compared with the N2TR group, NoF1 subgroup was the protective factor of CHB (HR = 0.652, 95% CI 0.476–0.893, p = 0.008) and PCT (HR = 0.605, 95% CI 0.398–0.918, p = 0.018) (Table 3). It’s worth noting that males were associated with a lower risk of PCT (HR = 0.654, 95% CI 0.465–0.919, p = 0.015) in patients with truncating mutations.

TABLE 3

Table 3. Univariate and multivariate Cox regression analysis of age-related VHL tumor risks of NoF1 vs. N2TR.

On the basis of predecessors, we classified patients with missense mutations (M group) by mutation locations into three subgroups, whose missense mutations were located in the HIF-binding region (residues 65–117, MH subgroup), Elongin C binding region (residues 158–184, ME subgroup), and others region (MO subgroup), respectively (Ong et al., 2007; Lee et al., 2016; Liu et al., 2018; Salama et al., 2019). It’s a remarkable fact that the number of MO group was far less than the other two groups in our cohort. Compared with MH subgroup and ME subgroup, patients in MO subgroup seemed to be more likely to have PHEO (p < 0.001, compared with MH subgroup, Figure 3A; p = 0.067, compared with ME subgroup, Figure 3B). Though age-associated risks of PHEO between MO and ME subgroup didn’t show significant differences (p = 0.067), it may be due to the limited sample size of MO subgroup. If we regarded ME and MO as a whole subgroup, this subgroup had a higher penetrance of PHEO than MH subgroup (p < 0.001, Figure 3C), or than those in combined all other patients (p < 0.001, Figure 3D). No significant difference in age-related risks of other symptoms was discovered among the three M subgroups.

FIGURE 3

Figure 3. Comparison of age-related risks of PHEO between MO subgroup and MH subgroup (A), MO subgroup and ME subgroup (B), MH subgroup and those in combination with MO and ME subgroups (C), those in combined with MO and ME subgroups and all other VHL patients (D). PHEO, pheochromocytoma; MH, missense mutations in HIF-α binding site; ME, missense mutations in elongin C binding site; MO, missense mutations subgroup in other sites.

We diagnosed patients affected by PHEO before death or the end of follow-up as type 2 VHL disease and those not affected by PHEO as type 1 VHL disease. Finally, 87 of 577 patients were diagnosed as type 2 VHL disease. Type 1 patients were associated with higher risks of CHB (p = 0.001), RA (p = 0.051), and lower risks of RCC (p = 0.053) than those with type 2 disease. The risks of VHL-associated tumors (p = 0.068) or PCT (p = 0.083) between the two types didn’t show a significant difference. We could then come to the conclusion that appearance of PHEO was a protective factor for CHB, and its effect on RA, RCC was still unclear (Supplementary Figure 3).

It’s generally acknowledged that residues prior to 113 of VHL protein were coded by exon 1, residues 115–153 by exon 2, and residues after 155 by exon 3 (Gnarra et al., 1996; Hansen et al., 2002). Patients in M group with mutations in exon 2 vs. exons 1 or 3 had a higher age-related risk of PCT on univariate (HR = 2.045, 95%CI 1.249–3.348, p = 0.004) and multivariate analyses (HR = 2.053, 95%CI 1.253–3.362, p = 0.004; Table 4).

TABLE 4

Table 4. The univariate and multivariate cox analyses of age-related risks of PCT.

Difference in Survival Between Groups

The estimated median lifespan of all 577 patients was 62 year. Because MO subgroup was too small and all of them were alive, we didn’t have access to their survival data. The differences between the median survival of MH subgroup and ME subgroup weren’t statistically significant (62 vs. 63 year, p = 0.259). Missense mutations carriers seemed to have longer estimated median survival than others (63 vs. 58 year), but the difference was not statistically significant (p = 0.090). The estimated median survival of patients in NoF1 subgroup was 9 years longer than those in NoF2 subgroup (64 vs. 55 year, p = 0.037, Figure 4). A NoF2 mutation was an independent risk factor for VHL patients’ survival in comparison with NoF1 mutations on univariate (HR = 1.900, 1.024–3.525, p = 0.042) and multivariate (HR = 1.888, 95%CI 1.017–3.506, p = 0.044) Cox regression analysis (Supplementary Table 3). Type 2 patients had significantly better survival than type 1 patients (p = 0.005, Figure 4B). Type 2 disease was an independent protective factor of survival, compared with type 1 patients (HR = 0.361, 95%CI 0.177–0.740, p = 0.005, Supplementary Table 3).

FIGURE 4

Figure 4. Comparison of survival in patients with different mutation groups. (A) Comparison of survival in patients from NoF1 and NoF2 subgroups. (B) Comparison of survival in type 1 and type 2 patients. (C) Comparison of CHB-specific survival in patients with MH and ME mutations. (D) Comparison of RCC-specific survival in patients with MH and ME mutations. (E) Comparison of CHB-specific survival in patients with NoF1 and NoF2 mutations. (F) Comparison of RCC-specific survival in patients with NoF1 and NoF2 mutations. NoF1, nonsense or frameshift mutations before residue 117 of VHL protein; NoF2, nonsense or frameshift mutations after residue 117 of VHL protein; CHB, central nervous system hemangioblastoma; MH, missense mutations in HIF-α binding site; ME, missense mutations in elongin C binding site; RCC, renal cell carcinoma.

CHB and RCC Specific Survival Among Subgroups

A total of 131 of 577 patients had died of some causes, which were categorized as CHB related, RCC related, and other causes of death. CHB and RCC were the most common death causes of VHL patients, which were shown as Supplementary Table 4. A total of 90.8% (119 of 131) of deaths were VHL related. Their death age ranged from 11 to 71 with a median of 42. As mentioned before, we didn’t have access to survival data of patients in the MO subgroup. We failed to find differences between MH subgroup and ME subgroup for CHB- or RCC-specific survival (p = 0.067 and p = 0.724, respectively). There was no difference for CHB or RCC specific survival between NoF1 and NoF2 subgroups (p = 0.413 and p = 0.062, respectively, Figure 4).

Discussion

The VHL disease is characterized by lifelong multiorgan and multisite tumors risks. Genotype-phenotype correlations are crucial to clinical genetic counseling. Several previous studies have focused on genotype-phenotype correlations of the VHL disease, but there was still a lack of authoritative consensus (Ong et al., 2007; Mettu et al., 2010; Lee et al., 2016; Lomte et al., 2018). The typical classification methods divide VHL patients into type 1 and type 2 based on the penetrance of PHEO. In our study, type 1 patients were associated with higher lifetime risks of CHB and worse survival compared with type 2 patients. These results suggested that type 1 or 2 VHL disease classification was an easy but rough way of risk stratification in clinicians’ daily work.

Our study suggests that MO subgroup should be treated as a separate subgroup independent of ME and MH subgroup, and NoF1 subgroup should be seen as a separate subgroup with a relatively better prognosis independent of the N2TR subgroup. These discoveries have the potential to improve and detail clinical counseling of VHL disease and reveal the underlying molecular mechanisms.

In our cohort, we subdivided the M group into MH, ME, and MO subgroups. The age-related risks of PHEO in MO and ME subgroups were higher than in other patients. MO subgroup seemed to be more likely to be affected by PHEO than MH subgroup, but no significant difference was discovered. The small number of people in MO subgroup and not long enough follow-up time limited more detailed genotype-phenotype research. There wasn’t significant difference in lifespan, CHB, or RCC specific survival between MH and ME subgroup. This conclusion wasn’t consistent with some previous papers, which held the opinions that MH subgroup had worse CHB specific survival compared with other missense mutations carriers (Lee et al., 2016; Liu et al., 2018). It might explain this divergence that we defined the new MO subgroup from non-MH missense mutations carriers and we expanded the sample size. However, there is a lack of evidence at molecular levels to support the rationality of this classification method, which is mainly based on VHL-HIF pathway. VHL protein is able to interact with several other molecular pathways involving a diverse array of genes’ transcription, deposition of the extracellular matrix, and the regulation of the microtubule cytoskeleton (Czyzyk-Krzeska and Meller, 2004; Frew and Krek, 2008; Xue et al., 2012). In vitro and in vivo studies revealed that some mutations involved in type 2 VHL disease retained the ability to regulate HIF-α subunits, suggesting that there must be more underlying mechanisms that are responsible for VHL disease (Clifford et al., 2001; Hoffman et al., 2001). In other words, there are several pathways besides the VHL-HIF pathway through which the mutant VHL gene leads to the disease (Roe et al., 2006; Forman et al., 2009). So, the classification of mutations based on the VHL-HIF pathway is not suitable for all VHL-associated phenotypes. Mutations at different locations are likely to change the intrinsic structure and interfaces of the VHL protein (Minervini et al., 2019), resulting in deregulation of the corresponding pathways and tumorigenesis. Further studies should classify mutations based on their effects on the structure and functions of the VHL protein rather than merely on the basis of mutation locations.

Ong et al. (2007) had made nonsense and frameshift mutations and deletions as two independent types of mutations and found nonsense and frameshift mutations were associated with higher lifetime risks of RCC. In our cohort, the NoF group was divided into NoF1 and NoF2 subgroups to explore phenotype heterogeneity in NoF group. NoF1 subgroup had lower age-related risks of overall VHL-associated tumors, CHB, and PCT than N2TR subgroup and longer estimated lifespans than NoF2 subgroup.

Although RA and CHB all belonged to hemangioblastomas, one’s appearance failed to remind higher lifetime risks of the other.

According to our conclusion, it’s consistent with previous studies that CHB is the most common symptom and CHB and RCC are the lead causes of death (Binderup et al., 2017). On this account, the management of VHL disease should focus on the control of CHB and RCC. In our cohort, the estimated median lifespan of all patients is 62 year. There isn’t access to survival information of MO subgroup.

For PCT, patients in M group with mutations in exon 2 vs. exons 1 or 3 were under higher age-related risks. This conclusion is inconsistent with Tirosh et al. (2018) whose work supported mutations in exon 3 were more likely affected by pancreatic neuroendocrine tumors (PNETs). Perhaps it’s because the definition of PCT is different from the definition of PNETs. In oTR group, males had a lower risk of PCT than females. In view of the divergences and uncertainties, more research is needed to explore the risk factors of PCT in VHL patients.

At present, the main management measures of VHL disease are active surveillance and surgical intervention when necessary. Clinicians need to weigh the risks of tumor progression and metastasis against the risks of surgical intervention and complication. There is a lack of scientific and authoritative surgical indications for each of five main symptoms to refer to. Existing standards often focus mainly on clinical parameters but ignore mutation types of patients. VHL associated endolymphatic sac tumors (ELSTs) are potential to cause audio vestibular morbidity even when an ELST mass can’t be identified on imaging examinations. So, excision is suggested as soon as an ELST or ELST associated morbidity is discovered. VHL associated CHB are usually multisite, of which natural history is still unpredictable. According to the current mainstream view, only symptomatic tumors need to be resected (Dornbos et al., 2018). It’s generally recognized that RCC of VHL patients whose longest diameter is below 3 cm don’t need surgical interventions (Duffey et al., 2004; Carlo et al., 2019). Peng et al. (2019) held the opinion that the criterion should be relaxed up to 4 cm. PNETs of VHL patients were considered not at risk of metastasis with a diameter below 1.2 cm and high risks of malignancy and metastasis with a diameter greater than 3 cm. For patients with a PNET diameter between 1.2 and 3 cm, only those with VHL missense mutations were observed to develop metastasis (Tirosh et al., 2018). However, these results aren’t enough to make an accepted criterion for PNET surgical intervention. In the following studies, the VHL genotype should be taken into consideration to determine whether surgical treatments are needed or not, to predict the risk of tumor recurrence, progression, and metastasis, and to predict therapeutic effects of non-surgical interventions.

In this study, we describe a modified genotype-phenotype correlation of VHL disease on the basis of previous studies, according to the VHL-HIF pathway. The correlation tries to subgroup truncating mutations and missense mutations, predicting tumor penetrance and survival of VHL patients. Our results will work as a VHL disease risk stratification method and provide more detailed information for genetic counseling.

Data Availability Statement

The datasets for this article are not publicly available because we have a responsibility to protect patients’ privacy. Requests to access the datasets should be directed to KG, gongkan_pku@126.com.

Ethics Statement

The studies involving human participants were reviewed and approved by the Institutional Ethics Committee of the Peking University First Hospital. Written informed consent to participate in this study was provided by all patients or their legal guardians.

Author Contributions

JQ: conceptualization. KM and JZ: data curation. JQ and KZ: formal analysis and original draft writing. KG: funding acquisition. KG, YG, and LC: methodology and draft review and revision. KG and LC: project administration and supervision. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China (No. 81872081) and the Fundamental Research Funds for the Central Universities (Grant BMU2018JI002).

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.

Acknowledgments

We would like to give our sincere appreciation to all the patients in our cohort and their families for their support and cooperation.

Supplementary Material

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

References

Binderup, M. L. M., Jensen, A. M., Budtz-Jørgensen, E., and Bisgaard, M. L. (2017). Survival and causes of death in patients with von Hippel-Lindau disease. J. Med. Genet. 54, 11–18. doi: 10.1136/jmedgenet-2016-104058

PubMed Abstract | CrossRef Full Text | Google Scholar

Carlo, M. I., Hakimi, A. A., Stewart, G. D., Bratslavsky, G., Brugarolas, J., Chen, Y.-B., et al. (2019). Familial kidney cancer: implications of new syndromes and molecular insights. Eur. Urol. 76, 754–764. doi: 10.1016/j.eururo.2019.06.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Clifford, S. C., Cockman, M. E., Smallwood, A. C., Mole, D. R., Woodward, E. R., Maxwell, P. H., et al. (2001). Contrasting effects on HIF-1alpha regulation by disease-causing pVHL mutations correlate with patterns of tumourigenesis in von Hippel-Lindau disease. Hum. Mol. Genet. 10, 1029–1038. doi: 10.1093/hmg/10.10.1029

PubMed Abstract | CrossRef Full Text | Google Scholar

Czyzyk-Krzeska, M. F., and Meller, J. (2004). von Hippel-Lindau tumor suppressor: not only HIF’s executioner. Trends Mol. Med. 10, 146–149. doi: 10.1016/j.molmed.2004.02.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Dornbos, D., Kim, H. J., Butman, J. A., and Lonser, R. R. (2018). Review of the neurological implications of von hippel-lindau disease. JAMA Neurol. 75, 620–627. doi: 10.1001/jamaneurol.2017.4469

PubMed Abstract | CrossRef Full Text | Google Scholar

Duffey, B. G., Choyke, P. L., Glenn, G., Grubb, R. L., Venzon, D., Linehan, W. M., et al. (2004). The relationship between renal tumor size and metastases in patients with von hippel-lindau disease. J. Urol. 172, 63–65. doi: 10.1097/01.ju.0000132127.79974.3f

CrossRef Full Text | Google Scholar

Feldman, D. E., Thulasiraman, V., Ferreyra, R. G., and Frydman, J. (1999). Formation of the VHL-elongin BC tumor suppressor complex is mediated by the chaperonin TRiC. Mol. Cell 4, 1051–1061. doi: 10.1016/s1097-2765(00)80233-6

CrossRef Full Text | Google Scholar

Forman, J. R., Worth, C. L., Bickerton, G. R. J., Eisen, T. G., and Blundell, T. L. (2009). Structural bioinformatics mutation analysis reveals genotype-phenotype correlations in von Hippel-Lindau disease and suggests molecular mechanisms of tumorigenesis. Proteins 77, 84–96. doi: 10.1002/prot.22419

PubMed Abstract | CrossRef Full Text | Google Scholar

Frew, I. J., and Krek, W. (2008). pVHL: a multipurpose adaptor protein. Sci. Signal. 1:e30. doi: 10.1126/scisignal.124pe30

PubMed Abstract | CrossRef Full Text | Google Scholar

Gnarra, J. R., Duan, D. R., Weng, Y., Humphrey, J. S., Chen, D. Y., Lee, S., et al. (1996). Molecular cloning of the von Hippel-Lindau tumor suppressor gene and its role in renal carcinoma. Biochim. Biophys. Acta 1242, 201–210. doi: 10.1016/0304-419x(95)00012-5

CrossRef Full Text | Google Scholar

Haddad, N. M. N., Cavallerano, J. D., and Silva, P. S. (2013). Von hippel-lindau disease: a genetic and clinical review. Semin. Ophthalmol. 28, 377–386. doi: 10.3109/08820538.2013.825281

PubMed Abstract | CrossRef Full Text | Google Scholar

Hansen, W. J., Ohh, M., Moslehi, J., Kondo, K., Kaelin, W. G., and Welch, W. J. (2002). Diverse effects of mutations in exon II of the von Hippel-Lindau (VHL) tumor suppressor gene on the interaction of pVHL with the cytosolic chaperonin and pVHL-dependent ubiquitin ligase activity. Mol. Cell. Biol. 22, 1947–1960. doi: 10.1128/mcb.22.6.1947-1960.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoffman, M. A., Ohh, M., Yang, H., Klco, J. M., Ivan, M., and Kaelin, W. G. (2001). von Hippel-Lindau protein mutants linked to type 2C VHL disease preserve the ability to downregulate HIF. Hum. Mol. Genet. 10, 1019–1027. doi: 10.1093/hmg/10.10.1019

PubMed Abstract | CrossRef Full Text | Google Scholar

Hon, W.-C., Wilson, M. I., Harlos, K., Claridge, T. D. W., Schofield, C. J., Pugh, C. W., et al. (2002). Structural basis for the recognition of hydroxyproline in HIF-1 alpha by pVHL. Nature 417, 975–978. doi: 10.1038/nature00767

PubMed Abstract | CrossRef Full Text | Google Scholar

Hong, B., Ma, K., Zhou, J., Zhang, J., Wang, J., Liu, S., et al. (2019). Frequent mutations of VHL gene and the clinical phenotypes in the largest chinese cohort with von hippel-lindau disease. Front. Genet. 10:867. doi: 10.3389/fgene.2019.00867

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, J.-S., Lee, J.-H., Lee, K. E., Kim, J. H., Hong, J. M., Ra, E. K., et al. (2016). Genotype-phenotype analysis of von Hippel-Lindau syndrome in Korean families: HIF-α binding site missense mutations elevate age-specific risk for CNS hemangioblastoma. BMC Med. Genet. 17:48. doi: 10.1186/s12881-016-0306-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, S.-J., Wang, J.-Y., Peng, S.-H., Li, T., Ning, X.-H., Hong, B.-A., et al. (2018). Genotype and phenotype correlation in von Hippel-Lindau disease based on alteration of the HIF-α binding site in VHL protein. Genet. Med. Off. J. Am. Coll. Med. Genet. 20, 1266–1273. doi: 10.1038/gim.2017.261

PubMed Abstract | CrossRef Full Text | Google Scholar

Lomte, N., Kumar, S., Sarathi, V., Pandit, R., Goroshi, M., Jadhav, S., et al. (2018). Genotype phenotype correlation in Asian Indian von Hippel-Lindau (VHL) syndrome patients with pheochromocytoma/paraganglioma. Fam. Cancer 17, 441–449. doi: 10.1007/s10689-017-0058-y

PubMed Abstract | CrossRef Full Text | Google Scholar

McNeill, A., Rattenberry, E., Barber, R., Killick, P., MacDonald, F., and Maher, E. R. (2009). Genotype-phenotype correlations in VHL exon deletions. Am. J. Med. Genet. A 149A, 2147–2151. doi: 10.1002/ajmg.a.33023

PubMed Abstract | CrossRef Full Text | Google Scholar

Mettu, P., Agrón, E., Samtani, S., Chew, E. Y., and Wong, W. T. (2010). Genotype-phenotype correlation in ocular von Hippel-Lindau (VHL) disease: the effect of missense mutation position on ocular VHL phenotype. Invest. Ophthalmol. Vis. Sci. 51, 4464–4470. doi: 10.1167/iovs.10-5223

PubMed Abstract | CrossRef Full Text | Google Scholar

Minervini, G., Quaglia, F., Tabaro, F., and Tosatto, S. C. E. (2019). Genotype-phenotype relations of the von Hippel-Lindau tumor suppressor inferred from a large-scale analysis of disease mutations and interactors. PLoS Comput. Biol. 15:e1006478. doi: 10.1371/journal.pcbi.1006478

PubMed Abstract | CrossRef Full Text | Google Scholar

Nielsen, S. M., Rhodes, L., Blanco, I., Chung, W. K., Eng, C., Maher, E. R., et al. (2016). Von hippel-lindau disease: genetics and role of genetic counseling in a multiple neoplasia syndrome. J. Clin. Oncol. Off. J. Am. Soc. Clin. Oncol. 34, 2172–2181. doi: 10.1200/JCO.2015.65.6140

PubMed Abstract | CrossRef Full Text | Google Scholar

Ohh, M., Park, C. W., Ivan, M., Hoffman, M. A., Kim, T. Y., Huang, L. E., et al. (2000). Ubiquitination of hypoxia-inducible factor requires direct binding to the beta-domain of the von Hippel-Lindau protein. Nat. Cell Biol. 2, 423–427. doi: 10.1038/35017054

PubMed Abstract | CrossRef Full Text | Google Scholar

Ong, K. R., Woodward, E. R., Killick, P., Lim, C., Macdonald, F., and Maher, E. R. (2007). Genotype-phenotype correlations in von Hippel-Lindau disease. Hum. Mutat. 28, 143–149. doi: 10.1002/humu.20385

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, X., Chen, J., Wang, J., Peng, S., Liu, S., Ma, K., et al. (2019). Natural history of renal tumours in von Hippel-Lindau disease: a large retrospective study of Chinese patients. J. Med. Genet. 56, 380–387. doi: 10.1136/jmedgenet-2018-105567

PubMed Abstract | CrossRef Full Text | Google Scholar

Roe, J.-S., Kim, H., Lee, S.-M., Kim, S.-T., Cho, E.-J., and Youn, H.-D. (2006). p53 stabilization and transactivation by a von Hippel-Lindau protein. Mol. Cell 22, 395–405. doi: 10.1016/j.molcel.2006.04.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Salama, Y., Albanyan, S., Szybowska, M., Bullivant, G., Gallinger, B., Giles, R. H., et al. (2019). Comprehensive characterization of a Canadian cohort of von Hippel-Lindau disease patients. Clin. Genet. 96, 461–467. doi: 10.1111/cge.13613

PubMed Abstract | CrossRef Full Text | Google Scholar

Shuin, T., Yamasaki, I., Tamura, K., Okuda, H., Furihata, M., and Ashida, S. (2006). Von Hippel-Lindau disease: molecular pathological basis, clinical criteria, genetic testing, clinical features of tumors and treatment. Jpn. J. Clin. Oncol. 36, 337–343. doi: 10.1093/jjco/hyl052

PubMed Abstract | CrossRef Full Text | Google Scholar

Tirosh, A., Sadowski, S. M., Linehan, W. M., Libutti, S. K., Patel, D., Nilubol, N., et al. (2018). Association of VHL genotype with pancreatic neuroendocrine tumor phenotype in patients with von hippel–lindau disease. JAMA Oncol. 4, 124–126. doi: 10.1001/jamaoncol.2017.3428

PubMed Abstract | CrossRef Full Text | Google Scholar

Xue, J., Lv, D., Jiao, S., Zhao, W., Li, X., Sun, H., et al. (2012). pVHL mediates K63-linked ubiquitination of nCLU. PLoS One 7:e0035848. doi: 10.1371/journal.pone.0035848

PubMed Abstract | CrossRef Full Text | Google Scholar