Functional data analysis approaches for genotype-phenotype association studies from next-generation sequencing

Next-generation DNA sequencing platforms can effectively detect the entire spectrum of genomic variation and is emerging to be a major tool for systematic exploration of the universe of variants and interactions in the entire genome. However, the data produced by next-generation sequencing technologies will suffer from three basic problems: sequence errors, assembly errors, and missing data. Current statistical methods for genetic analysis are well suited for detecting the association of common variants, but are less suitable to rare variants. This raises great challenge for sequence-based genetic studies of complex diseases.^ This research dissertation utilized genome continuum model as a general principle, and stochastic calculus and functional data analysis as tools for developing novel and powerful statistical methods for next generation of association studies of both qualitative and quantitative traits in the context of sequencing data, which finally lead to shifting the paradigm of association analysis from the current locus-by-locus analysis to collectively analyzing genome regions.^ In this project, the functional principal component (FPC) methods coupled with high-dimensional data reduction techniques will be used to develop novel and powerful methods for testing the associations of the entire spectrum of genetic variation within a segment of genome or a gene regardless of whether the variants are common or rare.^ The classical quantitative genetics suffer from high type I error rates and low power for rare variants. To overcome these limitations for resequencing data, this project used functional linear models with scalar response to develop statistics for identifying quantitative trait loci (QTLs) for both common and rare variants. To illustrate their applications, the functional linear models were applied to five quantitative traits in Framingham heart studies. ^ This project proposed a novel concept of gene-gene co-association in which a gene or a genomic region is taken as a unit of association analysis and used stochastic calculus to develop a unified framework for testing the association of multiple genes or genomic regions for both common and rare alleles. The proposed methods were applied to gene-gene co-association analysis of psoriasis in two independent GWAS datasets which led to discovery of networks significantly associated with psoriasis.^

Cataract surgery in cancer patients and its impact on quality of life as measured by the National Eye Institute Visual Functioning Questionnaire - 25 (NEI-VFQ-25)

Purpose. The measurement of quality of life has become an important topic in healthcare and in the allocation of limited healthcare resources. Improving the quality of life (QOL) in cancer patients is paramount. Cataract removal and lens implantation appears to improve patient well-being of cancer patients, though a formal measurement has never been published in the US literature. In this current study, National Eye Institute Visual Functioning Questionnaire (NEI-VFQ-25), a validated vision quality of life metric, was used to study the change in vision-related quality of life in cancer patients who underwent cataract extraction with intraocular lens implantation. ^ Methods. Under an IRB approved protocol, cancer patients who underwent cataract surgery with intraocular lens implantation (by a single surgeon) from December 2008 to March 2011, and who had completed a pre- and postoperative NEI-VFQ-25 were retrospectively reviewed. Post-operative data was collected at their routine 4-6 week post-op visit. Patients' demographics, cancer history, their pre and postoperative ocular examinations, visual acuities, and NEI-VFQ-25 with twelve components were included in the evaluation. The responses were evaluated using the Student t test, Spearman correlation and Wilcoxon signed rank test. ^ Results. 63 cases of cataract surgery (from 54 patients) from the MD Anderson Cancer Center were included in the study. Cancer patients had a significant improvement in the visual acuity (P<0.0001) postoperatively, along with a significant increase in vision-related quality of life (P<0.0001). Patients also had a statistically significant improvement in ten of the twelve subcategories which are addressed in the NEI-VFQ-25. ^ Conclusions. In our study, cataract extraction and intraocular implantation showed a significant impact on the vision-related quality of life in cancer patients. Although this study includes a small sample size, it serves as a positive pilot study to evaluate and quantify the impact of a surgical intervention on QOL in cancer patients and may help to design a larger study to measure vision related QOL per dollar spent for health care cost in cancer patients.^

Exploration of genetic susceptibility to pancreatic cancer: A GWAS approach

Pancreatic cancer is the 4th most common cause for cancer death in the United States, accompanied by less than 5% five-year survival rate based on current treatments, particularly because it is usually detected at a late stage. Identifying a high-risk population to launch an effective preventive strategy and intervention to control this highly lethal disease is desperately needed. The genetic etiology of pancreatic cancer has not been well profiled. We hypothesized that unidentified genetic variants by previous genome-wide association study (GWAS) for pancreatic cancer, due to stringent statistical threshold or missing interaction analysis, may be unveiled using alternative approaches. To achieve this aim, we explored genetic susceptibility to pancreatic cancer in terms of marginal associations of pathway and genes, as well as their interactions with risk factors. We conducted pathway- and gene-based analysis using GWAS data from 3141 pancreatic cancer patients and 3367 controls with European ancestry. Using the gene set ridge regression in association studies (GRASS) method, we analyzed 197 pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Using the logistic kernel machine (LKM) test, we analyzed 17906 genes defined by University of California Santa Cruz (UCSC) database. Using the likelihood ratio test (LRT) in a logistic regression model, we analyzed 177 pathways and 17906 genes for interactions with risk factors in 2028 pancreatic cancer patients and 2109 controls with European ancestry. After adjusting for multiple comparisons, six pathways were marginally associated with risk of pancreatic cancer ( P < 0.00025): Fc epsilon RI signaling, maturity onset diabetes of the young, neuroactive ligand-receptor interaction, long-term depression (Ps < 0.0002), and the olfactory transduction and vascular smooth muscle contraction pathways (P = 0.0002; Nine genes were marginally associated with pancreatic cancer risk (P < 2.62 × 10−5), including five reported genes (ABO, HNF1A, CLPTM1L, SHH and MYC), as well as four novel genes (OR13C4, OR 13C3, KCNA6 and HNF4 G); three pathways significantly interacted with risk factors on modifying the risk of pancreatic cancer (P < 2.82 × 10−4): chemokine signaling pathway with obesity ( P < 1.43 × 10−4), calcium signaling pathway (P < 2.27 × 10−4) and MAPK signaling pathway with diabetes (P < 2.77 × 10−4). However, none of the 17906 genes tested for interactions survived the multiple comparisons corrections. In summary, our current GWAS study unveiled unidentified genetic susceptibility to pancreatic cancer using alternative methods. These novel findings provide new perspectives on genetic susceptibility to and molecular mechanisms of pancreatic cancer, once confirmed, will shed promising light on the prevention and treatment of this disease. ^

Geological and geochemical characteristics of the Tersang gold deposit, Malaysia: implications for genesis and gold exploration

The Central gold belt of peninsular Malaysia comprises a number of gold deposits located in the east of the N-S striking Bentong-Raub Suture Zone. The Tersang gold deposit is one of the gold deposits in the gold belt and hosted in sandstone, rhyolite and breccia units. The deposit has an inferred resource of 528,000 ounces of gold. The geochronology of the Tersang deposit has been newly constrained by LA ICP-MS U-Pb zircon dating. The maximum depositional age of the host sedimentary rocks ranges from Early Carboniferous to Early Permian (261.5 ± 4.9 Ma to 333.5 ± 2.5 Ma) for the host sandstone and Late Triassic for the rhyolite intrusion (218.8 ± 1.7 Ma). Textural characteristics of pyrite have revealed five types including (1) Euhedral to subhedral pyrite with internal fracturing and porous cores located in the sandstone layers (pyrite 1); (2) Anhedral pyrite overgrowths on pyrite 1 and disseminated in stage 1 vein (pyrite 2); (3) Fracture-filled or vein pyrite located in stages 1 and 2 vein (pyrite 3); (4) Euhedral pyrite with internal fractures also located in stage 2 vein (pyrite 4); and (5) Subhedral clean pyrite located in the rhyolite intrusion (pyrite 5). Based on pyrite mapping and spot analyses, two main stages of gold enrichment are documented from the Tersang gold deposit. Gold in sandstone-hosted pyrite 1 (mean 4.3 ppm) shows best correlation with Bi and Pb (as evidenced on pyrite maps). In addition, gold in pyrite 3 (mean 8 ppm) located in stage 2 vein shows a good correlation with As, Ag, Sb, Cu, Tl, and Pb. In terms of gold exploration, we suggest that elements such as As, Ag, Sb, Cu, Tl, Bi, and Pb associated with Au may serve as vectoring tools in gold exploration. Our new geological, structural, geochemical and isotopic data together with mineral paragenesis, pyrite chemistry and ore fluid characteristics indicate that the Tersang gold deposit is comparable to a sediment-hosted gold deposit. Our new genetic model suggests deposition of the Permo-Carboniferous sediments followed by intrusion of rhyolitic magma in the Late Triassic. At a later stage, gold mineralisation overprinted the rhyolite intrusion and the sandstone.

(Table 1) Stable isotope data of chemoherm carbonates at Hydrate Ridge on the Cascadia continental margin

Two active chemoherm build-ups growing freely up into the oceanic water column, the Pinnacle and the South East-Knoll Chemoherms, have been discovered at Hydrate Ridge on the Cascadia continental margin. These microbially-mediated carbonate formations rise above the seafloor by several tens of meters and display a pinnacle-shaped morphology with steep flanks. The recovered rocks are pure carbonates dominated by aragonite. Based on fabric and mineralogic composition different varieties of authigenic aragonite can be distinguished. Detailed visual and petrographic investigations unambiguously reveal the involvement of microbes during the formation of the carbonates. The fabric of the cryptocrystalline and fibrous aragonite can be described as thrombolitic. Fossilized microbial filaments in the microcrystalline aragonite indicate the intimate relationship between microbes and carbonates. The strongly 13C-depleted carbon isotope values of the samples (as low as -48.1 per mill PDB) are characteristic of methane as the major carbon source for the carbonate formation. The methane-rich fluids from which the carbonates are precipitated originate most probably from a gas reservoir below the bottom-simulating reflector (BSR) and rise through fault systems. The d18O values of the aragonitic chemoherm carbonates are substantially higher (as high as 5.0 per mill PDB) than the expected equilibrium value for an aragonite forming from ambient seawater (3.5 per mill PDB). As a first approximation this indicates formation from glacial ocean water but other factors are considered as well. A conceptual model is presented for the precipitation of these chemoherm carbonates based on in situ observations and the detailed petrographic investigation of the carbonates. This model explains the function of the consortium of archaea and sulfate-reducing bacteria that grows on the carbonates performing anaerobic oxidation of methane (AOM) and enabling the precipitation of the chemoherms above the seafloor surrounded by oxic seawater. Beggiatoa mats growing on the surface of the chemoherms oxidize the sulfide provided by sulfate-dependent anaerobic oxidation of methane within an oxic environment. The contact between Beggiatoa and the underlying microbial consortium represents the interface between the overlying oxic water column and an anoxic micro-environment where carbonate formation takes place.