Gene and Pathway-Based Analysis: Second Wave of Genome-wide Association Studies

Despite great success of GWAS in identification of common genetic variants associated with complex diseases, the current GWAS have focused on single SNP analysis. However, single SNP analysis often identifies a number of the most significant SNPs that account for only a small proportion of the genetic variants and offers limited understanding of complex diseases. To overcome these limitations, we propose gene and pathway-based association analysis as a new paradigm for GWAS. As a proof of concept, we performed a comprehensive gene and pathway-based association analysis for thirteen published GWAS. Our results showed that the proposed new paradigm for GWAS not only identified the genes that include significant SNPs found by single SNP analysis, but also detected new genes in which each single SNP conferred small disease risk, but their joint actions were implicated in the development of diseases. The results also demonstrated that the new paradigm for GWAS was able to identify biologically meaningful pathways associated with the diseases which were confirmed by gene-set rich analysis using gene expression data

Genetic variants in ELOVL2 and HSD17B12 predict melanoma‐specific survival

Fatty acids play a key role in cellular bioenergetics, membrane biosynthesis and intracellular signaling processes and thus may be involved in cancer development and progression. In the present study, we comprehensively assessed associations of 14,522 common single‐nucleotide polymorphisms (SNPs) in 149 genes of the fatty‐acid synthesis pathway with cutaneous melanoma disease‐specific survival (CMSS). The dataset of 858 cutaneous melanoma (CM) patients from a published genome‐wide association study (GWAS) by The University of Texas M.D. Anderson Cancer Center was used as the discovery dataset, and the identified significant SNPs were validated by a dataset of 409 CM patients from another GWAS from the Nurses’ Health and Health Professionals Follow‐up Studies. We found 40 noteworthy SNPs to be associated with CMSS in both discovery and validation datasets after multiple comparison correction by the false positive report probability method, because more than 85% of the SNPs were imputed. By performing functional prediction, linkage disequilibrium analysis, and stepwise Cox regression selection, we identified two independent SNPs of ELOVL2 rs3734398 T>C and HSD17B12 rs11037684 A>G that predicted CMSS, with an allelic hazards ratio of 0.66 (95% confidence interval = 0.51–0.84 and p = 8.34 × 10−4) and 2.29 (1.55–3.39 and p = 3.61 × 10−5), respectively. Finally, the ELOVL2 rs3734398 variant CC genotype was found to be associated with a significantly increased mRNA expression level. These SNPs may be potential markers for CM prognosis, if validated by additional larger and mechanistic studies

Prognostic relevance of acquired uniparental disomy in serous ovarian cancer

BACKGROUND: Acquired uniparental disomy (aUPD) can lead to homozygosity for tumor suppressor genes or oncogenes. Our purpose is to determine the frequency and profile aUPD regions in serous ovarian cancer (SOC) and investigated the association of aUPD with clinical features and patient outcomes.METHODS: We analyzed single nucleotide polymorphism (SNP) array-based genotyping data on 532 SOC specimens from The Cancer Genome Atlas database to identify aUPD regions. Cox univariate regression and Cox multivariate proportional hazards analyses were performed for survival analysis.RESULTS: We found that 94.7% of SOC samples harbored aUPD; the most common aUPD regions were in chromosomes 17q (76.7%), 17p (39.7%), and 13q (38.3%). In Cox univariate regression analysis, two independent regions of aUPD on chromosome 17q (A and C), and whole-chromosome aUPD were associated with shorter overall survival (OS), and five regions on chromosome 17q (A, D-G) and BRCA1 were associated with recurrence-free survival time. In Cox multivariable proportional hazards analysis, whole-chromosome aUPD was associated with shorter OS. One region of aUPD on chromosome 22q (B) was associated with unilateral disease. A statistically significant association was found between aUPD at TP53 loci and homozygous mutation of TP53 (p < 0.0001).CONCLUSIONS: aUPD is a common event and some recurrent loci are associated with a poor outcome for patients with serous ovarian cancer

Genetic variation in genes for the xenobiotic-metabolizing enzymes CYP1A1, EPHX1, GSTM1, GSTT1, and GSTP1 and susceptibility to colorectal cancer in Lynch syndrome.

Individuals with Lynch syndrome are predisposed to cancer due to an inherited DNA mismatch repair gene mutation. However, there is significant variability observed in disease expression likely due to the influence of other environmental, lifestyle, or genetic factors. Polymorphisms in genes encoding xenobiotic-metabolizing enzymes may modify cancer risk by influencing the metabolism and clearance of potential carcinogens from the body. In this retrospective analysis, we examined key candidate gene polymorphisms in CYP1A1, EPHX1, GSTT1, GSTM1, and GSTP1 as modifiers of age at onset of colorectal cancer among 257 individuals with Lynch syndrome. We found that subjects heterozygous for CYP1A1 I462V (c.1384A\u3eG) developed colorectal cancer 4 years earlier than those with the homozygous wild-type genotype (median ages, 39 and 43 years, respectively; log-rank test P = 0.018). Furthermore, being heterozygous for the CYP1A1 polymorphisms, I462V and Msp1 (g.6235T\u3eC), was associated with an increased risk for developing colorectal cancer [adjusted hazard ratio for AG relative to AA, 1.78; 95% confidence interval, 1.16-2.74; P = 0.008; hazard ratio for TC relative to TT, 1.53; 95% confidence interval, 1.06-2.22; P = 0.02]. Because homozygous variants for both CYP1A1 polymorphisms were rare, risk estimates were imprecise. None of the other gene polymorphisms examined were associated with an earlier onset age for colorectal cancer. Our results suggest that the I462V and Msp1 polymorphisms in CYP1A1 may be an additional susceptibility factor for disease expression in Lynch syndrome because they modify the age of colorectal cancer onset by up to 4 years

Functional photoacoustic microscopy of pH

pH is a tightly regulated indicator of metabolic activity. In mammalian systems, imbalance of pH regulation may result from or result in serious illness. Even though the regulation system of pH is very robust, tissue pH can be altered in many diseases such as cancer, osteoporosis and diabetes mellitus. Traditional high-resolution optical imaging techniques, such as confocal microscopy, routinely image pH in cells and tissues using pH sensitive fluorescent dyes, which change their fluorescence properties with the surrounding pH. Since strong optical scattering in biological tissue blurs images at greater depths, high-resolution pH imaging is limited to penetration depths of 1mm. Here, we report photoacoustic microscopy (PAM) of commercially available pH-sensitive fluorescent dye in tissue phantoms. Using both opticalresolution photoacoustic microscopy (OR-PAM), and acoustic resolution photoacoustic microscopy (AR-PAM), we explored the possibility of recovering the pH values in tissue phantoms. In this paper, we demonstrate that PAM was capable of recovering pH values up to a depth of 2 mm, greater than possible with other forms of optical microscopy

Variable Selection with False Discovery Control

Technological advances that allow routine identification of high-dimensional risk factors have led to high demand for statistical techniques that enable full utilization of these rich sources of information for genome-wide association studies (GWAS). Variable selection for censored outcome data as well as control of false discoveries (i.e. inclusion of irrelevant variables) in the presence of high-dimensional predictors present serious challenges. In the context of survival analysis with high-dimensional covariates, this paper develops a computationally feasible method for building general risk prediction models, while controlling false discoveries. We have proposed a high-dimensional variable selection method by incorporating stability selection to control false discovery. Comparisons between the proposed method and the commonly used univariate and Lasso approaches for variable selection reveal that the proposed method yields fewer false discoveries. The proposed method is applied to study the associations of 2,339 common single-nucleotide polymorphisms (SNPs) with overall survival among cutaneous melanoma (CM) patients. The results have confirmed that BRCA2 pathway SNPs are likely to be associated with overall survival, as reported by previous literature. Moreover, we have identified several new Fanconi anemia (FA) pathway SNPs that are likely to modulate survival of CM patients

Genetic variants in the vitamin D pathway genes VDBP and RXRA modulate cutaneous melanoma disease-specific survival

Single nucleotide polymorphisms (SNPs) in the vitamin D pathway genes have been implicated in cutaneous melanoma (CM) risk, but their role in CM disease-specific survival (DSS) remains obscure. We comprehensively analyzed the prognostic roles of 2669 common SNPs in the vitamin D pathway genes using data from a published genome-wide association study (GWAS) at The University of Texas M.D. Anderson Cancer Center (MDACC) and then validated the SNPs of interest in another GWAS from the Nurses' Health Study and Health Professionals Follow-up Study. Among the 2669 SNPs, 203 were significantly associated with DSS in MDACC dataset (P C was associated with a better DSS [combined hazards ratio (HR) = 0.66]; and the same for RXRA rs7850212 C > A (combined HR = 0.38), which were further confirmed by the Fine and Gray competing-risks regression model. Further bioinformatics analyses indicated that these loci may modulate corresponding gene methylation status

Genetic variants in the PIWI-piRNA pathway gene DCP1A predict melanoma disease-specific survival

The Piwi-piRNA pathway is important for germ cell maintenance, genome integrity, DNA methylation and retrotransposon control and thus may be involved in cancer development. In this study, we comprehensively analyzed prognostic roles of 3,116 common SNPs in PIWI-piRNA pathway genes in melanoma disease-specific survival. A published genome-wide association study (GWAS) by The University of Texas M.D. Anderson Cancer Center was used to identify associated SNPs, which were later validated by another GWAS from the Harvard Nurses' Health Study and Health Professionals Follow-up Study. After multiple testing correction, we found that there were 27 common SNPs in two genes (PIWIL4 and DCP1A) with false discovery rate < 0.2 in the discovery dataset. Three tagSNPs (i.e., rs7933369 and rs508485 in PIWIL4; rs11551405 in DCP1A) were replicated. The rs11551405 A allele, located at the 3' UTR microRNA binding site of DCP1A, was associated with an increased risk of melanoma disease-specific death in both discovery dataset [adjusted Hazards ratio (HR) = 1.66, 95% confidence interval (CI) = 1.21-2.27, p =1.50 × 10-3 ] and validation dataset (HR = 1.55, 95% CI = 1.03-2.34, p = 0.038), compared with the C allele, and their meta-analysis showed an HR of 1.62 (95% CI, 1.26-2.08, p =1.55 × 10-4 ). Using RNA-seq data from the 1000 Genomes Project, we found that DCP1A mRNA expression levels increased significantly with the A allele number of rs11551405. Additional large, prospective studies are needed to validate these findings