A novel risk model of three SUMOylation genes based on RNA expression for potential prognosis and treatment sensitivity prediction in kidney cancer.

Introduction: Kidney cancer is one of the most common and lethal urological malignancies. Discovering a biomarker that can predict prognosis and potential drug treatment sensitivity is necessary for managing patients with kidney cancer. SUMOylation is a type of posttranslational modification that could impact many tumor-related pathways through the mediation of SUMOylation substrates. In addition, enzymes that participate in the process of SUMOylation can also influence tumorigenesis and development. Methods: We analyzed the clinical and molecular data which were obtanied from three databases, The Cancer Genome Atlas (TCGA), the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium (CPTAC), and ArrayExpress. Results: Through analysis of differentially expressed RNA based on the total TCGA-KIRC cohort, it was found that 29 SUMOylation genes were abnormally expressed, of which 17 genes were upregulated and 12 genes were downregulated in kidney cancer tissues. A SUMOylation risk model was built based on the discovery TCGA cohort and then validated successfully in the validation TCGA cohort, total TCGA cohort, CPTAC cohort, and E-TMAB-1980 cohort. Furthermore, the SUMOylation risk score was analyzed as an independent risk factor in all five cohorts, and a nomogram was constructed. Tumor tissues in different SUMOylation risk groups showed different immune statuses and varying sensitivity to the targeted drug treatment. Discussion: In conclusion, we examined the RNA expression status of SUMOylation genes in kidney cancer tissues and developed and validated a prognostic model for predicting kidney cancer outcomes using three databases and five cohorts. Furthermore, the SUMOylation model can serve as a biomarker for selecting appropriate therapeutic drugs for kidney cancer patients based on their RNA expression.

Telomere-related gene risk model for prognosis and drug treatment efficiency prediction in kidney cancer.

Kidney cancer is one of the most common urological cancers worldwide, and kidney renal clear cell cancer (KIRC) is the major histologic subtype. Our previous study found that von-Hippel Lindau (VHL) gene mutation, the dominant reason for sporadic KIRC and hereditary kidney cancer-VHL syndrome, could affect VHL disease-related cancers development by inducing telomere shortening. However, the prognosis role of telomere-related genes in kidney cancer has not been well discussed. In this study, we obtained the telomere-related genes (TRGs) from TelNet. We obtained the clinical information and TRGs expression status of kidney cancer patients in The Cancer Genome Atlas (TCGA) database, The International Cancer Genome Consortium (ICGC) database, and the Clinical Proteomic Tumor Analysis Consortium (CPTAC) database. Totally 353 TRGs were differential between tumor and normal tissues in the TCGA-KIRC dataset. The total TCGA cohort was divided into discovery and validation TCGA cohorts and then using univariate cox regression, lasso regression, and multivariate cox regression method to conduct data analysis sequentially, ten TRGs (ISG15, RFC2, TRIM15, NEK6, PRKCQ, ATP1A1, ELOVL3, TUBB2B, PLCL1, NR1H3) risk model had been constructed finally. The kidney patients in the high TRGs risk group represented a worse outcome in the discovery TCGA cohort (p<0.001), and the result was validated by these four cohorts (validation TCGA cohort, total TCGA cohort, ICGC cohort, and CPTAC cohort). In addition, the TRGs risk score is an independent risk factor for kidney cancer in all these five cohorts. And the high TRGs risk group correlated with worse immune subtypes and higher tumor mutation burden in cancer tissues. In addition, the high TRGs risk group might benefit from receiving immune checkpoint inhibitors and targeted therapy agents. Moreover, the proteins NEK6, RF2, and ISG15 were upregulated in tumors both at the RNA and protein levels, while PLCL1 and PRKCQ were downregulated. The other five genes may display the contrary expression status at the RNA and protein levels. In conclusion, we have constructed a telomere-related genes risk model for predicting the outcomes of kidney cancer patients, and the model may be helpful in selecting treatment agents for kidney cancer patients.

Lupus susceptibility region containing CTLA4 rs17268364 functionally reduces CTLA4 expression by binding EWSR1 and correlates IFN-α signature.

BACKGROUND: Dysregulation of T cells mediated immune responses is a hallmark in the development of systemic lupus erythematosus (SLE). Recent genome wide association study (GWAS) revealed the genetic contribution of variants located in the cytotoxic T lymphocyte-associated protein-4 (CTLA4)-inducible T cell co-stimulator (ICOS) intergenic region to SLE susceptibility. Our aim is to find a functional variant in this region. METHODS: The genetic association results in the CTLA4-ICOS region from previous GWAS were adopted to select the potential variant which was further replicated in two independent cohorts (Henan cohort 2053 SLE patients and 1845 healthy controls, Beijing cohort 2303 SLE patients and 19,262 healthy). In order to explore the functional significance in SLE, bioinformatics with validation experiments (including electrophoretic mobility shift assay and luciferase reporter assay) and mRNA expression analysis were also performed. RESULTS: A variant located in the CTLA4-ICOS intergenic region, rs17268364, was associated with susceptibility to SLE patients in Chinese populations (risk allele, pmeta = 7.02×10-11, OR 1.19, 95%CI 1.13-1.26). The bioinformatics suggested that rs17268364 might affect the expression of CTLA4, not ICOS. The rs17268364 risk G allele containing sequence reduced the expression of the reporter gene by binding transcriptional repressor Ewing sarcoma breakpoint region 1 (EWSR1). Following genotype-mRNA expression, the analysis also showed the risk allele of rs17268364 was associated with low CTLA4 expression in lupus nephritis (LN) patients. Healthy individuals carrying rs17268364 risk G allele was significantly correlated with higher levels of IFN-α signature including increased lymphocyte antigen 6E (LY6E) (p=0.031), interferon-stimulated gene 15 (ISG15) (p=0.038), interferon regulatory factor 9 (IRF9) (p=0.028), and interferon regulatory factor 5 (IRF5) (p=0.040) mRNA expression. CONCLUSIONS: The present study confirmed the functional role of rs17268364 in the CTLA4-ICOS intergenic region that increased SLE susceptibility in the Chinese population.

Association of the PINX1 Variant rs6984094, Which Lengthens Telomeres, with Systemic Lupus Erythematosus Susceptibility in Chinese Populations.

A recent genome-wide association study (GWAS) of Asian ancestry reported that single nucleotide polymorphism (SNP) in TERT (telomerase reverse transcriptase) was associated with systemic lupus erythematosus (SLE). TERT has a critical role in maintaining the chromosomal stability and the length of telomere. Given that only a small portion of the genetic heritability of SLE has been explained so far, we aimed to identify novel loci in telomere-related genes responsible for SLE susceptibility in Chinese populations. We performed a comprehensive genetic association analysis of SLE with telomere-related genes. To identify functional significance, we analyzed the publicly available HaploReg v4.1 and RegulomeDB databases. Differential gene expression analysis was also performed using ArrayExpress. A novel signal of PINX1 rs6984094 was identified (P discovery = 4.13 × 10-2, OR = 0.58, 95% CI 0.35-0.98) and successfully replicated (P replication = 5.73 × 10-3, OR = 0.45, 95% CI 0.26-0.81). Multiple layers of functional analysis suggested that the PINX1 rs6984094 risk T allele exhibited increased nuclear protein binding. We also observed an increased expression of PINX1 mRNA in peripheral blood mononuclear cells from SLE patients compared with healthy controls. Overall, we observed a novel genetic association between PINX1 (encodes the PinX1 protein, an inhibitory telomerase enzyme that lengthens telomeres) and SLE susceptibility in Chinese populations.

The ZNF76 rs10947540 polymorphism associated with systemic lupus erythematosus risk in Chinese populations.

Systemic lupus erythematosus (SLE) is a typical autoimmune disease with a strong genetic disposition. Genetic studies have revealed that single-nucleotide polymorphisms (SNPs) in zinc finger protein (ZNF)-coding genes are associated with susceptibility to autoimmune diseases, including SLE. The objective of the current study was to evaluate the correlation between ZNF76 gene polymorphisms and SLE risk in Chinese populations. A total of 2801 individuals (1493 cases and 1308 controls) of Chinese Han origin were included in this two-stage genetic association study. The expression of ZNF76 was evaluated, and integrated bioinformatic analysis was also conducted. The results showed that 28 SNPs were associated with SLE susceptibility in the GWAS cohort, and the association of rs10947540 was successfully replicated in the independent replication cohort (Preplication = 1.60 × 10-2, OR 1.19, 95% CI 1.03-1.37). After meta-analysis, the association between rs10947540 and SLE was pronounced (Pmeta = 9.62 × 10-6, OR 1.29, 95% CI 1.15-1.44). Stratified analysis suggested that ZNF76 rs10947540 C carriers were more likely to develop relatively high levels of serum creatinine (Scr) than noncarriers (CC + CT vs. TT, p = 9.94 × 10-4). The bioinformatic analysis revealed that ZNF76 rs10947540 was annotated as an eQTL and that rs10947540 was correlated with decreased expression of ZNF76. Remarkably, significantly reduced expression of ZNF76 was confirmed by expression data from both our laboratory and an array-based expression database. Taken together, these results suggest that ZNF76 rs10947540 is a possible susceptibility factor associated with SLE susceptibility. The mechanism underlying the relationship between ZNF76 and SLE pathogenesis still requires further investigation.

Identification of a Hypoxia-Related Gene Model for Predicting the Prognosis and Formulating the Treatment Strategies in Kidney Renal Clear Cell Carcinoma.

PURPOSE: The present study aimed to establish a hypoxia related genes model to predict the prognosis of kidney clear cell carcinoma (KIRC) patients using data accessed from The Cancer Genome Atlas (TCGA) database and International Cancer Genome Consortium (ICGC) database. METHODS: Patients' data were downloaded from the TCGA and ICGC databases, and hypoxia related genes were accessed from the Molecular Signatures Database. The differentially expressed genes were evaluated and then the differential expressions hypoxia genes were screened. The TCGA cohort was randomly divided into a discovery TCGA cohort and a validation TCGA cohort. The discovery TCGA cohort was used for constructing the hypoxia genes risk model through Lasso regression, univariate and multivariate Cox regression analysis. Receiver operating characteristic (ROC) curves were used to assess the reliability and sensitivity of our model. Then, we established a nomogram to predict the probable one-, three-, and five-year overall survival rates. Lastly, the Tumor Immune Dysfunction and Exclusion (TIDE) score of patients was calculated. RESULTS: We established a six hypoxia-related gene prognostic model of KIRC patients in the TCGA database and validated in the ICGC database. The patients with high riskscore present poorer prognosis than those with low riskscore in the three TCGA cohorts and ICGC cohort. ROC curves show our six-gene model with a robust predictive capability in these four cohorts. In addition, we constructed a nomogram for KIRC patients in the TCGA database. Finally, the high risk-group had a high TIDE score than the patients with low riskscore. CONCLUSIONS: We established a six hypoxia-related gene risk model for independent prediction of the prognosis of KIRC patients was established and constructed a robust nomogram. The different riskscores might be a biomarker for immunotherapy strategy.

SCUBE3 Is Likely a Susceptibility Gene for Systemic Lupus Erythematosus for Chinese Populations.

BACKGROUND: Systemic lupus erythematosus (SLE) is a complex autoimmune disease with strong genetic disposition with more than 100 susceptibility genes identified until now. However, our knowledge on SLE genetic background is still limited. The present study was aimed at evaluating the role of single nucleotide polymorphisms (SNPs) in SCUBE3, a TGF-ß signaling activator, with SLE susceptibility in Chinese populations. METHODS: A total of 2801 individuals (490 cases and 493 controls from GWAS cohort and 1003 cases and 815 controls from our cohort) were enrolled, and SNPs located 10 kb up- and downstream of SCUBE3 (chr6:35182190-35218609) were included in the genetic association study. Multiple layers of bioinformatics were conducted, and the levels of SCUBE3 expression were confirmed. RESULTS: Of the 31 SNPs in SCUBE3 tested, 24 SNPs were significantly associated with SLE at p ≤ 0.05. The top locus was rs1888822 with p = 8.74∗10-6 in the discovery cohort and was confirmed by the replication cohort with p = 0.012. Additionally, the levels of SCUBE3 mRNA expression were significantly lower in patients with SLE comparing with healthy controls (p = 4.28∗10-4). Further expression data from ArrayExpress showed that the expression of SCUBE3 was also lower in CD3+ T cells and B cells from patients with SLE. CONCLUSIONS: Our research revealed that variants in SCUBE3, which encode SCUBE3 as a TGF-ß signaling activator, can be considered as a new genetic susceptibility factor for systemic lupus erythematosus. And the reduced mRNA expression of SCUBE3 was first reported in patients with SLE.

A Three Protein-Coding Gene Prognostic Model Predicts Overall Survival in Bladder Cancer Patients.

Bladder cancer (BLCA) is the most common urinary tract tumor and is the 11th most malignant cancer worldwide. With the development of in-depth multisystem sequencing, an increasing number of prognostic molecular markers have been identified. In this study, we focused on the role of protein-coding gene methylation in the prognosis of BLCA. We downloaded BLCA clinical and methylation data from The Cancer Genome Atlas (TCGA) database and used this information to identify differentially methylated genes and construct a survival model using lasso regression. We assessed 365 cases, with complete information regarding survival status, survival time longer than 30 days, age, gender, and tumor characteristics (grade, stage, T, M, N), in our study. We identified 353 differentially methylated genes, including 50 hypomethylated genes and 303 hypermethylated genes. After annotation, a total of 227 genes were differentially expressed. Of these, 165 were protein-coding genes. Three genes (zinc finger protein 382 (ZNF382), galanin receptor 1 (GALR1), and structural maintenance of chromosomes flexible hinge domain containing 1 (SMCHD1)) were selected for the final risk model. Patients with higher-risk scores represent poorer survival than patients with lower-risk scores in the training set (HR = 2.37, 95% CI 1.43-3.94, p = 0.001), in the testing group (HR = 1.85, 95% CI 1.16-2.94, p = 0.01), and in the total cohort (HR = 2.06, 95% CI 1.46-2.90, p < 0.001). Further univariate and multivariate analyses using the Cox regression method were conducted in these three groups, respectively. All the results indicated that risk score was an independent risk factor for BLCA. Our study screened the different methylation protein-coding genes in the BLCA tissues and constructed a robust risk model for predicting the outcome of BLCA patients. Moreover, these three genes may function in the mechanism of development and progression of BLCA, which should be fully clarified in the future.

Single Nucleotide Polymorphisms in PPARD Associated with Systemic Lupus Erythematosus in Chinese Populations.

BACKGROUND: Systemic lupus erythematosus (SLE) is a multifactorial autoimmune disease characterized by apoptotic clearance deficiency provoking autoimmune responses and leading to multiple organ damage. PPAR-δ, encoded by the PPARD gene, was induced in macrophages promoting the timely disposal of apoptotic cells. Biological studies had provided solid foundation of PPARD involvement in SLE; it is worthwhile to further explore the genetic contribution of PPARD to SLE. METHODS: We performed a discovery-replication genetic association study. The discovery study was based on previous reported GWAS data. And the replication study was conducted in 1003 SLE patients and 815 healthy controls from Henan, Middle East of China. Further, we analyzed the eQTL effect to identify possible functional significance. RESULTS: In the genetic association analysis, we observed significant association between the risk C allele of rs4713853 (p = 0.03, OR 1.167, 95% CI 1.015-1.341) and increased SLE susceptibility. Moreover, individuals with the risk C allele were associated with lower expression of PPARD and DEF6. Our clinical analysis showed that SLE patients with the risk C allele of rs4713853 were more likely to present a higher proportion of anti-Sm antibody presence (CC+CT vs. TT, 20.0% vs. 14.2%, p = 0.039) and higher level of Scr (median inter quarter range CC+CT vs. TT, 56 48-71 vs. 54 46-64 µmol/L, p = 0.002). CONCLUSIONS: In conclusion, our study identified a novel association between PPARD rs4713853 and SLE susceptibility in Chinese populations. By integrating multiple layers of analysis, we suggested that PPARD might be a main candidate in the pathogenesis of SLE.

Genotype and phenotype correlation in von Hippel-Lindau disease based on alteration of the HIF-α binding site in VHL protein.

PURPOSE: Von Hippel-Lindau (VHL) disease is a rare hereditary cancer syndrome that reduces life expectancy. We aimed to construct a more valuable genotype-phenotype correlation based on alterations in VHL protein (pVHL). METHODS: VHL patients (n = 339) were recruited and grouped based on mutation types: HIF-α binding site missense (HM) mutations, non-HIF-α binding site missense (nHM) mutations, and truncating (TR) mutations. Age-related risks of VHL-associated tumors and patient survival were compared. RESULTS: Missense mutations conferred an increased risk of pheochromocytoma (HR = 1.854, p = 0.047) compared with truncating mutations. The risk of pheochromocytoma was lower in the HM group than in the nHM group (HR = 0.298, p = 0.003) but was similar between HM and TR groups (HR = 0.901, p = 0.810). Patients in the nHM group had a higher risk of pheochromocytoma (HR = 3.447, p < 0.001) and lower risks of central nervous system hemangioblastoma (CHB) (HR = 0.700, p = 0.045), renal cell carcinoma (HR = 0.610, p = 0.024), and pancreatic tumor (HR = 0.382, p < 0.001) than those in the combined HM and TR (HMTR) group. Moreover, nHM mutations were independently associated with better overall survival (HR = 0.345, p = 0.005) and CHB-specific survival (HR = 0.129, p = 0.005) than HMTR mutations. CONCLUSION: The modified genotype-phenotype correlation links VHL gene mutation, substrate binding site, and phenotypic diversity (penetrance and survival), and provides more accurate information for genetic counseling and pathogenesis studies.

Risk factors for survival in patients with von Hippel-Lindau disease.

BACKGROUND: Historically, von Hippel-Lindau (VHL) disease is characterised by a poor survival. Although genotype-phenotype correlation has been described in many studies, the risk factors for VHL survival remain unclear. This study aims to evaluate the median survival of Chinese patients with VHL disease and explore whether VHL survival is influenced by genetic and clinical factors. METHODS: In this retrospective study, we recruited 340 patients from 127 VHL families. Kaplan-Meier plot and Cox regression model were used to evaluate the median survival and assess how survival was influenced by birth year, birth order, sex, family history, mutation type, onset age and first presenting symptom. RESULTS: The estimated median life expectancy for Chinese patients with VHL disease was 62 years. Patients with early-onset age, positive family history and truncating mutation types had poorer overall and VHL-related survival. Patients with haemangioblastoma as their first presenting symptom were related to a higher risk of death from central nervous system haemangioblastoma than those with abdominal lesions (HR 8.84, 95% CI 2.04 to 38.37, P=0.004). CONCLUSIONS: This largest VHL survival analysis indicates that onset age, family history, mutation type and first presenting symptom have an effect on the survival of patients with VHL disease, which is helpful to genetic counselling and clinical decision-making.

Shorter telomere length increases age-related tumor risks in von Hippel-Lindau disease patients.

Von Hippel-Lindau (VHL) disease is a rare autosomal dominant cancer syndrome caused by alterations of VHL gene. Patients are predisposed to develop pheochromocytomas and solid or cystic tumors of the central nervous system, kidney, pancreas, and retina. Remarkable phenotypic heterogeneity exits in organ involvement and tumor onset age between and within VHL families. However, no reliable markers have been found to predict the age-related tumor risks in VHL patients. A large Chinese cohort composed of 300 VHL patients and 92 healthy family controls was enrolled in our study. Blood relative telomere length was measured in 184 patients and all the controls available for genomic DNA samples. Age-related risks for the five major VHL-associated tumors were evaluated using Kaplan-Meier plots and Cox regression analysis. Differences in clinical phenotype were observed between Chinese cohort and the United Kingdom cohort. VHL patients showed significantly shorter telomere length than healthy family controls(P = 0.0183), and a positive correlation was found between telomere length and onset age of the five major tumors, respectively. Moreover, patients in the shorter telomere group (age-adjusted telomere length ≤ 0.44) suffered higher age-related risks for VHL-associated central nervous system hemangioblastomas (HR: 1.879, P = 0.004), renal cell carcinoma (HR: 2.126, P = 0.002) and pancreatic cyst and neuroendocrine tumors (HR: 2.093, P = 0.001). These results indicate that blood shorter telomere length is a new biomarker for age-related tumor risks in VHL patients, which will be crucial to genetic counseling and future research about the role of telomere shortening in the pathogenesis of VHL-associated tumors.

Higher programmed cell death 1 ligand 1 (PD-L1) mRNA level in clear cell renal cell carcinomas is associated with a favorable outcome due to the active immune responses in tumor tissues.

Renal cell carcinoma is one of the most common urological tumors. The role of programmed cell death 1 ligand 1 (PD-L1) in renal cell carcinomas in predicting outcome of the patients is yet unclear. We analyzed the clinical and RNA-seq data of 522 kidney clear cell cancer, 259 kidney papillary cell carcinoma and 66 kidney chromophobe patients from The Cancer Genome Atlas (TCGA) database. In kidney clear cell cancer patients with high PD-L1 mRNA level and low PD-L1 mRNA level in tumors, the median overall survival periods were 45.0 and 37.1 months respectively (p=0.002). Multivariate Cox regression tests found that PD-L1 mRNA level in tumor was an independent predictor for overall survival status in kidney clear cell cancer patients (HR=0.7, 95% CI 0.5-0.9, p=0.007). However, no significant difference in overall survival status was found between high and low PD-L1 groups in kidney papillary cell carcinoma and kidney chromophobe cohorts. Gene-set enrichment analysis on the data from databases of TCGA and GSE53757 dataset in Gene Expression Omnibus databases showed that several pathways relating to immunological functions were activated in kidney clear cell cancers with high PD-L1 mRNA expression, and glycolysis and epithelial-mesenchymal transition pathways relating to tumor progression and metastasis were increased in kidney clear cell cancers with low PD-L1 mRNA level. In conclusion, higher PD-L1 mRNA level in kidney clear cell cancer tissues was associated with a favorable outcome due to the higher immunological responses in tumor tissues.

Association between FBP1 and hypoxia-related gene expression in clear cell renal cell carcinoma.

Fructose-1,6-bisphosphatase 1 (FBP1) is a rate-limiting enzyme in gluconeogenesis. Recently, the catalytic activity-independent function of FBP1, hypoxia-induced factor (HIF) repression in the nucleus, was identified. The aim of the present study was to investigate the association between FBP1 and hypoxia-related gene expression in clear cell renal cell carcinoma (ccRCC). The protein expression levels of FBP1, HIF-1α, HIF-2α, erythropoietin (EPO) and carbonic anhydrase IX (CA9) were assessed by immunohistochemical staining of ccRCC paraffin blocks from 123 patients using the tissue microarray technique. The expression level of FBP1 was then correlated with various clinicopathological factors, and the protein expression levels of HIF-1α, HIF-2α, EPO and CA9. Clinicopathological factors, including age, gender, T stage and Fuhrman grade, were not significantly different between patients with low and high FBP1 expression in ccRCC (P>0.05). FBP1 protein expression level was significantly correlated with the expression levels of HIF-1α (P=0.005) and EPO (P=0.010), but not significantly correlated with the expression levels of HIF-2α (P=0.123) and CA9 (P=0.513) in ccRCC tissues. The current findings confirm the association between FBP1 and hypoxia-related gene expression, and may facilitate understanding of the mechanisms of ccRCC tumorigenesis.

UBE2T silencing suppresses proliferation and induces cell cycle arrest and apoptosis in bladder cancer cells.

Ubiquitin-conjugating enzyme E2T (UBE2T), a member of the ubiquitin-conjugating E2 family in the ubiquitin-proteasome pathway, has been reported to be overexpressed in certain tumor types and to have an important role in the Fanconi anemia pathway. In the present study, the expression of UBE2T and its association with bladder cancer were investigated; to the best of our knowledge, this has not been reported previously. Immunohistochemistry and western blot analysis demonstrated that UBE2T was significantly upregulated in bladder cancer tissues and cell lines compared with adjacent normal bladder tissues and a normal human urinary tract epithelial cell line, respectively. UBE2T was detectable in the nuclei and cytoplasm of cancer cells, exhibiting stronger expression in the nuclei. A UBE2T-siRNA-expressing lentivirus was constructed and used to infect human bladder cancer 5637 cells, in order to examine the role of UBE2T in bladder cancer cell growth in vitro. The knockdown of UBE2T significantly decreased bladder cancer cell proliferation and colony formation. Furthermore, UBE2T silencing induced cell cycle arrest at G2/M phase and increased cell apoptosis. Therefore, UBE2T serves an important role in the growth of bladder cancer cells, and may be considered as a potential biomarker and therapeutic target for bladder cancer.

Telomere shortening is associated with genetic anticipation in Chinese Von Hippel-Lindau disease families.

Von Hippel-Lindau (VHL) disease is a rare autosomal dominant cancer syndrome. A phenomenon known as genetic anticipation has been documented in some hereditary cancer syndromes, where it was proved to relate to telomere shortening. Because studies of this phenomenon in VHL disease have been relatively scarce, we investigated anticipation in 18 Chinese VHL disease families. We recruited 34 parent-child patient pairs (57 patients) from 18 families with VHL disease. Onset age was defined as the age when any symptom or sign of VHL disease first appeared. Anticipation of onset age was analyzed by paired t test and the other two special tests (HV and RY2). Relative telomere length of peripheral leukocytes was measured in 29 patients and 325 healthy controls. Onset age was younger in child than in parent in 31 of the 34 parent-child pairs. Patients in the first generation had older onset age with longer age-adjusted relative telomere length, and those in the next generation had younger onset age with shorter age-adjusted relative telomere length (P < 0.001) in the 10 parent-child pairs from eight families with VHL disease. In addition, relative telomere length was shorter in the 29 patients with VHL disease than in the normal controls (P = 0.003). The anticipation may relate to the shortening of telomere length in patients with VHL in successive generations. These findings indicate that anticipation is present in families with VHL disease and may be helpful for genetic counseling for families with VHL disease families and for further understanding the pathogenesis of VHL disease.