Independent genomewide screens identify the tumor suppressor VTRNA2-1 as a human epiallele responsive to periconceptional environment

Background: Interindividual epigenetic variation that occurs systemically must be established prior to gastrulation in the very early embryo and, because it is systemic, can be assessed in easily biopsiable tissues. We employ two independent genome-wide approaches to search for such variants.

Results: First, we screen for metastable epialleles by performing genomewide bisulfite sequencing in peripheral blood lymphocyte (PBL) and hair follicle DNA from two Caucasian adults. Second, we conduct a genomewide screen for genomic regions at which PBL DNA methylation is affected by season of conception in rural Gambia. Remarkably, both approaches identify the genomically imprinted VTRNA2-1 as a top environmentally responsive epiallele. We demonstrate systemic and stochastic interindividual variation in DNA methylation at the VTRNA2-1 differentially methylated region in healthy Caucasian and Asian adults and show, in rural Gambians, that periconceptional environment affects offspring VTRNA2-1 epigenotype, which is stable over at least 10 years. This unbiased screen also identifies over 100 additional candidate metastable epialleles, and shows that these are associated with cis genomic features including transposable elements.

Conclusions: The non-coding VTRNA2-1 transcript (also called nc886) is a putative tumor suppressor and modulator of innate immunity. Thus, these data indicating environmentally induced loss of imprinting at VTRNA2-1 constitute a plausible causal pathway linking early embryonic environment, epigenetic alteration, and human disease. More broadly, the list of candidate metastable epialleles provides a resource for future studies of epigenetic variation and human disease.

KATNAL1 regulation of sertoli cell microtubule dynamics is essential for spermiogenesis and male fertility

Spermatogenesis is a complex process reliant upon interactions between germ cells (GC) and supporting somatic cells. Testicular Sertoli cells (SC) support GCs during maturation through physical attachment, the provision of nutrients, and protection from immunological attack. This role is facilitated by an active cytoskeleton of parallel microtubule arrays that permit transport of nutrients to GCs, as well as translocation of spermatids through the seminiferous epithelium during maturation. It is well established that chemical perturbation of SC microtubule remodelling leads to premature GC exfoliation demonstrating that microtubule remodelling is an essential component of male fertility, yet the genes responsible for this process remain unknown. Using a random ENU mutagenesis approach, we have identified a novel mouse line displaying male-specific infertility, due to a point mutation in the highly conserved ATPase domain of the novel KATANIN p60-related microtubule severing protein Katanin p60 subunit A-like1 (KATNAL1). We demonstrate that Katnal1 is expressed in testicular Sertoli cells (SC) from 15.5 days post-coitum (dpc) and that, consistent with chemical disruption models, loss of function of KATNAL1 leads to male-specific infertility through disruption of SC microtubule dynamics and premature exfoliation of spermatids from the seminiferous epithelium. The identification of KATNAL1 as an essential regulator of male fertility provides a significant novel entry point into advancing our understanding of how SC microtubule dynamics promotes male fertility. Such information will have resonance both for future treatment of male fertility and the development of non-hormonal male contraceptives.

Factors influencing success of clinical genome sequencing across a broad spectrum of disorders

To assess factors influencing the success of whole-genome sequencing for mainstream clinical diagnosis, we sequenced 217 individuals from 156 independent cases or families across a broad spectrum of disorders in whom previous screening had identified no pathogenic variants. We quantified the number of candidate variants identified using different strategies for variant calling, filtering, annotation and prioritization. We found that jointly calling variants across samples, filtering against both local and external databases, deploying multiple annotation tools and using familial transmission above biological plausibility contributed to accuracy. Overall, we identified disease-causing variants in 21% of cases, with the proportion increasing to 34% (23/68) for mendelian disorders and 57% (8/14) in family trios. We also discovered 32 potentially clinically actionable variants in 18 genes unrelated to the referral disorder, although only 4 were ultimately considered reportable. Our results demonstrate the value of genome sequencing for routine clinical diagnosis but also highlight many outstanding challenges.

Variation in DNA repair genes XRCC3, XRCC4, XRCC5 and susceptibility to myeloma

Cytogenetic analysis in myeloma reveals marked chromosomal instability. Both widespread genomic alterations and evidence of aberrant class switch recombination, the physiological process that regulates maturation of the antibody response, implicate the DNA repair pathway in disease pathogenesis. We therefore assessed 27 SNPs in three genes (XRCC3, XRCC4 and XRCC5) central to DNA repair in patients with myeloma and controls from the EpiLymph study and from an Irish hospital registry (n = 306 cases, 263 controls). For the haplotype-tagging SNP (htSNP) rs963248 in XRCC4, Allele A was significantly more frequent in cases than in controls (86.4 versus 80.8%; odds ratio 1.51; 95% confidence interval 1.10-2.08; P = 0.0133), as was the AA genotype (74 versus 65%) (P = 0.026). Haplotype analysis was performed using Unphased for rs963248 in combination with additional SNPs in XRCC4. The strongest evidence of association came from the A-T haplotype from rs963248-rs2891980 (P = 0.008). For XRCC5, the genotype GG from rs1051685 was detected in 10 cases from different national populations but in only one control (P = 0.015). This SNP is located in the 3'-UTR of XRCC5. Overall, these data provide support for the hypothesis that common variation in the genes encoding DNA repair proteins contributes to susceptibility to myeloma.

Shared common variants in prostate cancer and blood lipids

BACKGROUND: Epidemiological and clinical studies suggest comorbidity between prostate cancer (PCA) and cardiovascular disease (CVD) risk factors. However, the relationship between these two phenotypes is still not well understood. Here we sought to identify shared genetic loci between PCA and CVD risk factors.

METHODS: We applied a genetic epidemiology method based on conjunction false discovery rate (FDR) that combines summary statistics from different genome-wide association studies (GWAS), and allows identification of genetic overlap between two phenotypes. We evaluated summary statistics from large, multi-centre GWA studies of PCA (n=50 000) and CVD risk factors (n=200 000) [triglycerides (TG), low-density lipoprotein (LDL) cholesterol and high-density lipoprotein (HDL) cholesterol, systolic blood pressure, body mass index, waist-hip ratio and type 2 diabetes (T2D)]. Enrichment of single nucleotide polymorphisms (SNPs) associated with PCA and CVD risk factors was assessed with conditional quantile-quantile plots and the Anderson-Darling test. Moreover, we pinpointed shared loci using conjunction FDR.

RESULTS: We found the strongest enrichment of P-values in PCA was conditional on LDL and conditional on TG. In contrast, we found only weak enrichment conditional on HDL or conditional on the other traits investigated. Conjunction FDR identified altogether 17 loci; 10 loci were associated with PCA and LDL, 3 loci were associated with PCA and TG and additionally 4 loci were associated with PCA, LDL and TG jointly (conjunction FDR <0.01). For T2D, we detected one locus adjacent to HNF1B.

CONCLUSIONS: We found polygenic overlap between PCA predisposition and blood lipids, in particular LDL and TG, and identified 17 pleiotropic gene loci between PCA and LDL, and PCA and TG, respectively. These findings provide novel pathobiological insights and may have implications for trials using targeting lipid-lowering agents in a prevention or cancer setting.

Exome sequencing of prostate cancer supports the hypothesis of independent tumour origins

BACKGROUND: Prostate cancer (PCa) is a clinically and pathologically heterogeneous disease. The rapid development of sequencing technology has the potential to deliver new biomarkers with emphasis on aggressive disease and to revolutionise personalised cancer treatment. However, a prostate harbouring cancer commonly contains multiple separate tumour foci, with the potential to aggravate tumour sampling. The level of intraprostatic tumour heterogeneity remains to be determined.

OBJECTIVE: To determine the level of intraprostatic tumour heterogeneity through genome-wide, high-resolution profiling of multiple tumour samples from the same individual.

DESIGN, SETTINGS, AND PARTICIPANTS: Multiple tumour samples were obtained from four individuals following radical prostatectomy. One individual (SWE-1) contained >70% cancer cells in all tumour samples, whereas the other three (SWE-2 to SWE-4) required the use of laser capture microdissection for tumour cell enrichment. Subsequently, DNA was extracted from all tissue samples, and exome sequencing was performed. All tumour foci of SWE-1 were also profiled using a high-resolution array for the identification of copy number alterations (CNA).

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Shared somatic high-frequency single nucleotide variants (SNV) and CNAs were used to infer the level of intraprostatic tumour heterogeneity.

RESULTS AND LIMITATIONS: No high-frequency mutations, common for the three tumour samples of SWE-1, were identified. Ten randomly chosen positions were validated with Sanger sequencing in all foci, which verified the exome data. The high level of intraprostatic heterogeneity was consistent in all individuals. In total, three out of four individuals harboured tumours without an apparent common somatic denominator. Although we cannot exclude the presence of common structural rearrangements, a high-density array was used for the detection of deletions and amplifications in SWE-1, which agreed with the exome data.

CONCLUSIONS: We present evidence for the presence of somatically independent tumours within the same prostate. This finding will have implications for personalised cancer treatment and biomarker discovery.

Integrin genetic variants and stage-specific tumor recurrence in patients with stage II and III colon cancer

Integrins (ITGs) are key elements in cancer biology, regulating tumor growth, angiogenesis and lymphangiogenesis through interactions of the tumor cells with the microenvironment. Moving from the hypothesis that ITGs could have different effects in stage II and III colon cancer, we tested whether a comprehensive panel of germline single-nucleotide polymorphisms (SNPs) in ITG genes could predict stage-specific time to tumor recurrence (TTR). A total of 234 patients treated with 5-fluorouracil-based chemotherapy at the University of Southern California were included in this study. Whole-blood samples were analyzed for germline SNPs in ITG genes using PCR-restriction fragment length polymorphism or direct DNA sequencing. In the multivariable analysis, stage II colon cancer patients with at least one G allele for ITGB3 rs4642 had higher risk of recurrence (hazard ratio (HR)=4.027, 95% confidence interval (95% CI) 1.556-10.421, P=0.004). This association was also significant in the combined stage II-III cohort (HR=1.975, 95% CI 1.194-3.269, P=0.008). The predominant role of ITGB3 rs4642 in stage II diseases was confirmed using recursive partitioning, showing that ITGB3 rs4642 was the most important factor in stage II diseases. In contrast, in stage III diseases the combined analysis of ITGB1 rs2298141 and ITGA4 rs7562325 allowed to identify three distinct prognostic subgroups (P=0.009). The interaction between stage and the combined ITGB1 rs2298141 and ITGA4 rs7562325 on TTR was significant (P=0.025). This study identifies germline polymorphisms in ITG genes as independent stage-specific prognostic markers for stage II and III colon cancer. These data may help to select subgroups of patients who may benefit from ITG-targeted treatments.

Plastin Polymorphisms Predict Gender- and Stage-Specific Colon Cancer Recurrence after Adjuvant Chemotherapy

Tumor recurrence after curative resection remains a major problem in patients with locally advanced colorectal cancer treated with adjuvant chemotherapy. Genetic single-nucleotide polymorphisms (SNP) may serve as useful molecular markers to predict clinical outcomes in these patients and identify targets for future drug development. Recent in vitro and in vivo studies have demonstrated that the plastin genes PLS3 and LCP1 are overexpressed in colon cancer cells and play an important role in tumor cell invasion, adhesion, and migration. Hence, we hypothesized that functional genetic variations of plastin may have direct effects on the progression and prognosis of locally advanced colorectal cancer. We tested whether functional tagging polymorphisms of PLS3 and LCP1 predict time to tumor recurrence (TTR) in 732 patients (training set, 234; validation set, 498) with stage II/III colorectal cancer. The PLS3 rs11342 and LCP1 rs4941543 polymorphisms were associated with a significantly increased risk for recurrence in the training set. PLS3 rs6643869 showed a consistent association with TTR in the training and validation set, when stratified by gender and tumor location. Female patients with the PLS3 rs6643869 AA genotype had the shortest median TTR compared with those with any G allele in the training set [1.7 vs. 9.4 years; HR, 2.84; 95% confidence interval (CI), 1.32-6.1; P = 0.005] and validation set (3.3 vs. 13.7 years; HR, 2.07; 95% CI, 1.09-3.91; P = 0.021). Our findings suggest that several SNPs of the PLS3 and LCP1 genes could serve as gender- and/or stage-specific molecular predictors of tumor recurrence in stage II/III patients with colorectal cancer as well as potential therapeutic targets.

Association of common gene variants in the WNT/β-catenin pathway with colon cancer recurrence

Wnt/β-catenin signaling has a central role in the development and progression of most colon cancers (CCs). Germline variants in Wnt/β-catenin pathway genes may result in altered gene function and/or activity, thereby causing inter-individual differences in relation to tumor recurrence capacity and chemoresistance. We investigated germline polymorphisms in a comprehensive panel of Wnt/β-catenin pathway genes to predict time to tumor recurrence (TTR) in patients with stage III and high-risk stage II CC. A total of 234 patients treated with 5-fluorouracil-based chemotherapy were included in this study. Whole-blood samples were analyzed for putative functional germline polymorphisms in SFRP3, SFRP4, DKK2, DKK3, Axin2, APC, TCF7L2, WNT5B, CXXC4, NOTCH2 and GLI1 genes by PCR-based restriction fragment-length polymorphism or direct DNA sequencing. Polymorphisms with statistical significance were validated in an independent study cohort. The minor allele of WNT5B rs2010851 T>G was significantly associated with a shorter TTR (10.7 vs 4.9 years; hazard ratio: 2.48; 95% CI, 0.96-6.38; P=0.04) in high-risk stage II CC patients. This result remained significant in multivariate Cox's regression analysis. This study shows that the WNT5B germline variant rs2010851 was significantly identified as a stage-dependent prognostic marker for CC patients after 5-fluorouracil-based adjuvant therapy.

Gender-specific profiling in SCN1A polymorphisms and time-to-recurrence in patients with stage II/III colorectal cancer treated with adjuvant 5-fluoruracil chemotherapy

This study was designed to analyze the gender-related association between SCN1A polymorphisms (voltage-gated sodium channels; α-subunit) and time-to-recurrence (TTR) in patients with colorectal cancer (CRC) treated with 5-fluoruracil (5-FU)-based adjuvant chemotherapy. We enrolled from a prospective database patients with stage II and III CRC treated with adjuvant 5-FU-based chemotherapy. Genotypes for SCN1A rs3812718 and rs229877 were determined by direct DNA sequencing. One hundred twenty-seven males and 107 females were included in the study. In the univariate and multivariate analysis, the shortest TTR was associated with female patients carrying the rs3812718-TT genotype (hazard ratio (HR): 2.26 (95% confidence interval (CI): 0.89, 5.70), P=0.039) but with male patients carrying the rs3812718-CC genotype (HR: 0.49 (95% CI: 0.18, 1.38), P=0.048). For rs229877 the CT genotype was associated with a trend for shorter TTR in both gender populations. The study validated gender-dependent association between genomic SCN1A rs3812718 polymorphism and TTR in CRC patients treated with adjuvant 5-FU-based chemotherapy. This study confirms that voltage-gated Na+ channels may be a potential therapeutic target and a useful predictive biomarker before 5-FU infusion.

Pharmacogenetic angiogenesis profiling for first-line Bevacizumab plus oxaliplatin-based chemotherapy in patients with metastatic colorectal cancer

PURPOSE: There is substantial germline genetic variability within angiogenesis pathway genes, thereby causing interindividual differences in angiogenic capacity and resistance to antiangiogenesis therapy. We investigated germline polymorphisms in genes involved in VEGF-dependent and -independent angiogenesis pathways to predict clinical outcome and tumor response in metastatic colorectal cancer (mCRC) patients treated with bevacizumab and oxaliplatin-based chemotherapy.

EXPERIMENTAL DESIGN: A total of 132 patients treated with first-line bevacizumab and FOLFOX or XELOX were included in this study. Genomic DNA was isolated from whole-blood samples by PCR-RFLP or direct DNA sequencing. The endpoints of the study were progression-free survival (PFS), overall survival (OS), and response rate (RR).

RESULTS: The minor alleles of EGF rs444903 A>G and IGF-1 rs6220 A>G were associated with increased OS and remained significant in multivariate Cox regression analysis (HR: 0.52; 95% CI: 0.31-0.87; adjusted P = 0.012 and HR: 0.60; 95% CI: 0.36-0.99; adjusted P = 0.046, respectively). The minor allele of HIF1α rs11549465 C>T was significantly associated with increased PFS but lost its significance in multivariate analysis. CXCR1 rs2234671 G>C, CXCR2 rs2230054 T>C, EGFR rs2227983 G>A, and VEGFR-2 rs2305948 C>T predicted tumor response, with CXCR1 rs2234671 G>C remaining significant in multiple testing (P(act) = 0.003).

CONCLUSION: In this study, we identified common germline variants in VEGF-dependent and -independent angiogenesis genes predicting clinical outcome and tumor response in patients with mCRC receiving first-line bevacizumab and oxaliplatin-based chemotherapy.