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.^

The use of second-generation antipsychotics and the changes in physical growth in children and adolescents with perinatally acquired HIV.

Second-generation antipsychotics (SGAs) are increasingly prescribed to treat psychiatric symptoms in pediatric patients infected with HIV. We examined the relationship between prescribed SGAs and physical growth in a cohort of youth with perinatally acquired HIV-1 infection. Pediatric AIDS Clinical Trials Group (PACTG), Protocol 219C (P219C), a multicenter, longitudinal observational study of children and adolescents perinatally exposed to HIV, was conducted from September 2000 until May 2007. The analysis included P219C participants who were perinatally HIV-infected, 3-18 years old, prescribed first SGA for at least 1 month, and had available baseline data prior to starting first SGA. Each participant prescribed an SGA was matched (based on gender, age, Tanner stage, baseline body mass index [BMI] z score) with 1-3 controls without antipsychotic prescriptions. The main outcomes were short-term (approximately 6 months) and long-term (approximately 2 years) changes in BMI z scores from baseline. There were 236 participants in the short-term and 198 in the long-term analysis. In linear regression models, youth with SGA prescriptions had increased BMI z scores relative to youth without antipsychotic prescriptions, for all SGAs (short-term increase = 0.192, p = 0.003; long-term increase = 0.350, p < 0.001), and for risperidone alone (short-term = 0.239, p = 0.002; long-term = 0.360, p = 0.001). Participants receiving both protease inhibitors (PIs) and SGAs showed especially large increases. These findings suggest that growth should be carefully monitored in youth with perinatally acquired HIV who are prescribed SGAs. Future research should investigate the interaction between PIs and SGAs in children and adolescents with perinatally acquired HIV infection.

Biochemical and genetic characterization of the Saccharomyces cerevisiae Hsp110 molecular chaperone Sse1

The Ssel/Hsp110 molecular chaperones are a poorly understood subgroup of the Hsp70 chaperone family. Hsp70 can refold denatured polypeptides via a carboxyl-terminal peptide binding domain (PBD), which is regulated by nucleotide cycling in an amino-terminal ATPase domain. However, unlike Hsp70, both Sse1 and mammalian Hsp110 bind unfolded peptide substrates but cannot refold them. To test the in vivo requirement for interdomain communication, SSE1 alleles carrying amino acid substitutions in the ATPase domain were assayed for their ability to complement sse1Δ phenotypes. Surprisingly, all mutants predicted to abolish ATP hydrolysis complemented the temperature sensitivity of sse1Δ, whereas mutations in predicted ATP binding residues were non-functional. Remarkably, the two domains of Ssel when expressed in trans functionally complement the sse1Δ growth phenotype and interact by coimmunoprecipitation analysis, indicative of a novel type of interdomain communication. ^ Relatively little is known regarding the interactions and cellular functions of Ssel. Through co-immunoprecipitation analysis, we found that Ssel forms heterodimeric complexes with the abundant cytosolic Hsp70s Ssa and Ssb in vivo. Furthermore, these complexes can be efficiently reconstituted in vitro using purified proteins. The ATPase domains of Ssel and the Hsp70s were found to be critical for interaction as inactivating point mutations severely reduced interaction efficiency. Ssel stimulated Ssal ATPase activity synergistically with the co-chaperone Ydj1 via a novel nucleotide exchange activity. Furthermore, FES1, another Ssa nucleotide exchange factor, can functionally substitute for SSE1/2 when overexpressed, suggesting that Hsp70 nucleotide exchange is the fundamental role of the Sse proteins in yeast, and by extension, the Hsp110 homologs in mammals. ^ Cells lacking SSE1 were found to accumulate prepro-α-factor, but not the cotranslationally imported protein Kar2, similar to mutants in the Ssa chaperones. This indicates that the interaction between Ssel and Ssa is functionally significant in vivo. In addition, sse10 cells are compromised for cell wall strength, likely a result of decreased Hsp90 chaperone activity with the cell integrity MAP kinase SIC. Taken together, this work established that the Hsp110 family must be considered an essential component of Hsp70 chaperone biology in the eukaryotic cell.^

Bifurcation and singularity analysis of a molecular network for the induction of long-term memory.

Withdrawal reflexes of the mollusk Aplysia exhibit sensitization, a simple form of long-term memory (LTM). Sensitization is due, in part, to long-term facilitation (LTF) of sensorimotor neuron synapses. LTF is induced by the modulatory actions of serotonin (5-HT). Pettigrew et al. developed a computational model of the nonlinear intracellular signaling and gene network that underlies the induction of 5-HT-induced LTF. The model simulated empirical observations that repeated applications of 5-HT induce persistent activation of protein kinase A (PKA) and that this persistent activation requires a suprathreshold exposure of 5-HT. This study extends the analysis of the Pettigrew model by applying bifurcation analysis, singularity theory, and numerical simulation. Using singularity theory, classification diagrams of parameter space were constructed, identifying regions with qualitatively different steady-state behaviors. The graphical representation of these regions illustrates the robustness of these regions to changes in model parameters. Because persistent protein kinase A (PKA) activity correlates with Aplysia LTM, the analysis focuses on a positive feedback loop in the model that tends to maintain PKA activity. In this loop, PKA phosphorylates a transcription factor (TF-1), thereby increasing the expression of an ubiquitin hydrolase (Ap-Uch). Ap-Uch then acts to increase PKA activity, closing the loop. This positive feedback loop manifests multiple, coexisting steady states, or multiplicity, which provides a mechanism for a bistable switch in PKA activity. After the removal of 5-HT, the PKA activity either returns to its basal level (reversible switch) or remains at a high level (irreversible switch). Such an irreversible switch might be a mechanism that contributes to the persistence of LTM. The classification diagrams also identify parameters and processes that might be manipulated, perhaps pharmacologically, to enhance the induction of memory. Rational drug design, to affect complex processes such as memory formation, can benefit from this type of analysis.

SOCIAL AND BIOLOGICAL DETERMINANTS OF SHORT-TERM GROWTH VELOCITY AMONG PRESCHOOLERS OF RURAL GUATEMALA: A PATH ANALYSIS APPROACH (WEIGHT-GAIN, IRON DEFICIENCY, GASTROINTESTINAL SYMPTOMS)

This dissertation develops and tests through path analysis a theoretical model to explain how socioeconomic, socioenvironmental, and biologic risk factors simultaneously influence each other to further produce short-term, depressed growth in preschoolers. Three areas of risk factors were identified: child's proximal environment, maturational stage, and biological vulnerability. The theoretical model represented both the conceptual framework and the nature and direction of the hypotheses. Original research completed in 1978-80 and in 1982 provided the background data. It was analyzed first by nested-analysis of variance, followed by path analysis. The study provided evidence of mild iron deficiency and gastrointestinal symptomatology in the etiology of depressed, short-term weight gain. Also, there was evidence suggesting that family resources for material and social survival significantly contribute to the variability of short-term, age-adjusted growth velocity. These results challenge current views of unifocal intervention, whether for prevention or control. For policy formulations, though, the mechanisms underlying any set of interlaced relationships must be decoded. Theoretical formulations here proposed should be reassessed under a more extensive research design. It is suggested that studies should be undertaken where social changes are actually in progress; otherwise, nutritional epidemiology in developing countries operates somewhere between social reality and research concepts, with little grasp of its real potential. The study stresses that there is a connection between substantive theory, empirical observation, and policy issues. ^