STUDYING AGGREGATE FORMATION BY AMYOTROPHIC LATERAL SCLEROSIS-ASSOCIATED MUTANT SOD1 PROTEIN IN DROSOPHILA MODEL

A common pathological hallmark of most neurodegenerative disorders is the presence of protein aggregates in the brain. Understanding the regulation of aggregate formation is thus important for elucidating disease pathogenic mechanisms and finding effective preventive avenues and cures. Amyotrophic Lateral Sclerosis (ALS), also known as Lou Gehrig’s disease, is a selective neurodegenerative disorder predominantly affecting motor neurons. The majority of ALS cases are sporadic, however, mutations in superoxide dismutase 1 (SOD1) are responsible for about 20% of familial ALS (fALS). Mutated SOD1 proteins are prone to misfold and form protein aggregates, thus representing a good candidate for studying aggregate formation. The long-term goal of this project is to identify regulators of aggregate formation by mutant SOD1 and other ALS-associated disease proteins. The specific aim of this thesis project is to assess the possibility of using the well-established Drosophila model system to study aggregation by human SOD1 (hSOD1) mutants. To this end, using wild type and the three mutant hSOD1 (A4V, G85R and G93A) most commonly found among fALS, I have generated 16 different SOD1 constructs containing either eGFP or mCherry in-frame fluorescent reporters, established and tested both cell- and animal-based Drosophila hSOD1 models. The experimental strategy allows for clear visualization of ectopic hSOD1 expression as well as versatile co-expression schemes to fully investigate protein aggregation specifically by mutant hSOD1. I have performed pilot cell-transfection experiments and verified induced expression of hSOD1 proteins. Using several tissue- or cell type-specific Gal4 lines, I have confirmed the proper expression of hSOD1 from established transgenic fly lines. Interestingly, in both Drosophila S2 cells and different fly tissues including the eye and motor neurons, robust aggregate formation by either wild type or mutant hSOD1 proteins was not observed. These preliminary observations suggest that Drosophila might not be a good experimental organism to study aggregation and toxicity of mutant hSOD1 protein. Nevertheless this preliminary conclusion implies the potential existence of a potent protective mechanism against mutant hSOD1 aggregation and toxicity in Drosophila. Thus, results from my SOD1-ALS project in Drosophila will help future studies on how to best employ this classic model organism to study ALS and other human brain degenerative diseases.

Amyotrophic lateral sclerosis and Parkinsonism-dementia complex of Guam: Descriptive epidemiology, secular trends, and birth cohort effects on incidence, 1950--1989

This study describes the patterns of occurrence of amyotrophic lateral sclerosis (ALS) and parkinsonism-dementia complex (PDC) of Guam during 1950-1989. Both ALS and PDC occur with high frequency among the indigenous Chamorro population, first recognized in the early 1950's. Reports in the early 1980's indicated that both ALS and PDC were disappearing, due to a purported reduction in exposure to harmful environmental factors as a result of the dramatic changes in lifestyle that took place after World War II. However, this study provides compelling evidence that ALS and PDC have not disappeared on Guam and that rates for both are higher during 1980-1989 than previously reported.^ The patterns of occurrence for both ALS and PDC overlap in most respects: (1) incidence and mortality are decreasing; (2) median age at onset is increasing; (3) males are at increased risk for developing disease; (4) risk is higher for those residing in the south compared to the non-south; and (5) age-specific incidence is decreasing over time except in the oldest age groups.^ Age-specific incidence of ALS and PDC, separately and together, is generally higher for cohorts born before 1920 than for those born after 1920. A significant birth cohort effect on the incidence of PDC for the 1906-1915 birth cohort was found, but not for ALS and for ALS and PDC together. Whether or not a cohort effect, period effect, or both are associated with incidence of ALS and PDC cannot be determined from the data currently available and will require additional follow-up of individuals born after 1920.^ The epidemiological data amassed over this 40-year period provide evidence that supports an environmental exposure model for disease occurrence as opposed to a simple genetic or infectious disease model. Whether neurodegenerative disease in this population occurs as a consequence of a single exposure or is explained by a multifactorial model such as a genetic predisposition with some environmental interaction is yet to be determined. However, descriptive studies such as this can provide clues concerning timing and location of potential adverse exposures but cannot determine etiology, underscoring the urgent need for analytic studies of ALS and PDC to further investigate existing etiologic hypotheses and to test new hypotheses. ^

Prognosis of amyotrophic lateral sclerosis

Two cohorts of amyotrophic lateral sclerosis (ALS) patients were identified. One incidence-based cohort from Harris County, Texas with 97 cases, and the other a clinic referral series from an ALS clinic in Houston, Texas with 439 cases were followed-up to evaluate the prognosis of ALS. The overall Kaplan-Meier 3-year survival after diagnosis was similar, 0.287 for the incidence cohort and 0.313 for the referral cohort. However, the 5-year survival was much lower for the incidence cohort than the referral cohort (0.037 vs. 0.206). The large difference in 5-year survival was thought to be the results of a stronger unfavorable effect of the prognostic factors in the incidence cohort than in the referral cohort.^ Cohort-specific Weibull regression models were derived to evaluate the cohort-specific prognostic factors and survival probability with adjustment of certain prognostic factors.^ The major prognostic factors were: age at diagnosis, bulbar onset, black ethnicity, and positive family history of ALS in both cohorts. Female gender, simultaneous upper and lower extremities onset were specifically unfavorable factors in the incidence cohort. In the incidence cohort the prognosis was relatively favorable for cases with duration from onset to diagnosis longer than 4 months, however in the referral cohort the relatively favorable prognosis only occurred in cases with duration from onset to diagnosis 1 year or longer and was strongest in cases with duration 5 years and longer. Age at diagnosis modified the effect of bulbar onset in the incidence cohort but not in the referral cohort. The estimated survival with presence of an unfavorable prognostic factor identified in the incidence cohort was higher for the referral cohort than for the incidence cohort. Future studies are indicated to investigate the disease heterogeneity issue of ALS based on survival distribution of ALS. ^

AMPA RECEPTOR ALLOSTERISM: MEASUREMENT OF THE CONFORMATIONAL CHANGES IN THE LIGAND BINDING DOMAIN OF A FUNCTIONAL RECEPTOR

Ionotropic glutamate receptors are important excitatory neurotransmitter receptors in the mammalian central nervous system that have been implicated in a number of neuropathologies such as epilepsy, ischemia, and amyotrophic lateral sclerosis. Glutamate binding to an extracellular ligand binding domain initiates a series of structural changes that leads to the formation of a cation selective transmembrane channel, which consequently closes due to desensitization of the receptor. The crystal structures of the AMPA subtype of the glutamate receptor have been particularly useful in providing initial insight into the conformational changes in the ligand binding domain; however, these structures are limited by crystallographic constraint. To gain a clear picture of how agonist binding is coupled to channel activation and desensitization, it is essential to study changes in the ligand binding domain in a dynamic, physiological state. In this dissertation, a technique called Luminescence Resonance Energy Transfer was used to determine the conformational changes associated with activation and desensitization in a functional AMPA receptor (ÄN*-AMPA) that contains the ligand binding domain and transmembrane segments; ÄN*-AMPA has been modified such that fluorophores can be introduced at specific sites to serve as a readout of cleft closure or to establish intersubunit distances. Previous structural studies of cleft closure of the isolated ligand binding domain in conjunction with functional studies of the full receptor suggest that extent of cleft closure correlates with extent of activation. Here, LRET has been used to show that a similar relationship between cleft closure and activation is observed in the “full length” receptor showing that the isolated ligand binding domain is a good model of the domain in the full length receptor for changes within a subunit. Similar LRET investigations were used to study intersubunit distances specifically to probe conformational changes between subunits within a dimer in the tetrameric receptor. These studies show that the dimer interface is coupled in the open state, and decoupled in the desensitized state, similar to the isolated ligand binding domain crystal structure studies. However, we show that the apo state dimer interface is not pre-formed as in the crystal structure, hence suggesting a mechanism for functional transitions within the receptor based on LRET distances obtained.

CONFORMATIONAL CHANGES IN THE EXTRACELLULAR DOMAIN OF GLUTAMATE RECEPTORS

The family of membrane protein called glutamate receptors play an important role in the central nervous system in mediating signaling between neurons. Glutamate receptors are involved in the elaborate game that nerve cells play with each other in order to control movement, memory, and learning. Neurons achieve this communication by rapidly converting electrical signals into chemical signals and then converting them back into electrical signals. To propagate an electrical impulse, neurons in the brain launch bursts of neurotransmitter molecules like glutamate at the junction between neurons, called the synapse. Glutamate receptors are found lodged in the membranes of the post-synaptic neuron. They receive the burst of neurotransmitters and respond by fielding the neurotransmitters and opening ion channels. Glutamate receptors have been implicated in a number of neuropathologies like ischemia, stroke and amyotrophic lateral sclerosis. Specifically, the NMDA subtype of glutamate receptors has been linked to the onset of Alzheimer’s disease and the subsequent degeneration of neuronal cells. While crystal structures of AMPA and kainate subtypes of glutamate receptors have provided valuable information regarding the assembly and mechanism of activation; little is known about the NMDA receptors. Even the basic question of receptor assembly still remains unanswered. Therefore, to gain a clear understanding of how the receptors are assembled and how agonist binding gets translated to channel opening, I have used a technique called Luminescence Resonance Energy Transfer (LRET). LRET offers the unique advantage of tracking large scale conformational changes associated with receptor activation and desensitization. In this dissertation, LRET, in combination with biochemical and electrophysiological studies, were performed on the NMDA receptors to draw a correlation between structure and function. NMDA receptor subtypes GluN1 and GluN2A were modified such that fluorophores could be introduced at specific sites to determine their pattern of assembly. The results indicated that the GluN1 subunits assembled across each other in a diagonal manner to form a functional receptor. Once the subunit arrangement was established, this was used as a model to further examine the mechanism of activation in this subtype of glutamate receptor. Using LRET, the correlation between cleft closure and activation was tested for both the GluN1 and GluN2A subunit of the NMDA receptor in response to agonists of varying efficacies. These investigations revealed that cleft closure plays a major role in the mechanism of activation in the NMDA receptor, similar to the AMPA and kainate subtypes. Therefore, suggesting that the mechanism of activation is conserved across the different subtypes of glutamate receptors.

BAYESIAN STATISTICAL METHODS IN GENE-ENVIRONMENT AND GENE-GENE INTERACTION STUDIES

Complex diseases such as cancer result from multiple genetic changes and environmental exposures. Due to the rapid development of genotyping and sequencing technologies, we are now able to more accurately assess causal effects of many genetic and environmental factors. Genome-wide association studies have been able to localize many causal genetic variants predisposing to certain diseases. However, these studies only explain a small portion of variations in the heritability of diseases. More advanced statistical models are urgently needed to identify and characterize some additional genetic and environmental factors and their interactions, which will enable us to better understand the causes of complex diseases. In the past decade, thanks to the increasing computational capabilities and novel statistical developments, Bayesian methods have been widely applied in the genetics/genomics researches and demonstrating superiority over some regular approaches in certain research areas. Gene-environment and gene-gene interaction studies are among the areas where Bayesian methods may fully exert its functionalities and advantages. This dissertation focuses on developing new Bayesian statistical methods for data analysis with complex gene-environment and gene-gene interactions, as well as extending some existing methods for gene-environment interactions to other related areas. It includes three sections: (1) Deriving the Bayesian variable selection framework for the hierarchical gene-environment and gene-gene interactions; (2) Developing the Bayesian Natural and Orthogonal Interaction (NOIA) models for gene-environment interactions; and (3) extending the applications of two Bayesian statistical methods which were developed for gene-environment interaction studies, to other related types of studies such as adaptive borrowing historical data. We propose a Bayesian hierarchical mixture model framework that allows us to investigate the genetic and environmental effects, gene by gene interactions (epistasis) and gene by environment interactions in the same model. It is well known that, in many practical situations, there exists a natural hierarchical structure between the main effects and interactions in the linear model. Here we propose a model that incorporates this hierarchical structure into the Bayesian mixture model, such that the irrelevant interaction effects can be removed more efficiently, resulting in more robust, parsimonious and powerful models. We evaluate both of the 'strong hierarchical' and 'weak hierarchical' models, which specify that both or one of the main effects between interacting factors must be present for the interactions to be included in the model. The extensive simulation results show that the proposed strong and weak hierarchical mixture models control the proportion of false positive discoveries and yield a powerful approach to identify the predisposing main effects and interactions in the studies with complex gene-environment and gene-gene interactions. We also compare these two models with the 'independent' model that does not impose this hierarchical constraint and observe their superior performances in most of the considered situations. The proposed models are implemented in the real data analysis of gene and environment interactions in the cases of lung cancer and cutaneous melanoma case-control studies. The Bayesian statistical models enjoy the properties of being allowed to incorporate useful prior information in the modeling process. Moreover, the Bayesian mixture model outperforms the multivariate logistic model in terms of the performances on the parameter estimation and variable selection in most cases. Our proposed models hold the hierarchical constraints, that further improve the Bayesian mixture model by reducing the proportion of false positive findings among the identified interactions and successfully identifying the reported associations. This is practically appealing for the study of investigating the causal factors from a moderate number of candidate genetic and environmental factors along with a relatively large number of interactions. The natural and orthogonal interaction (NOIA) models of genetic effects have previously been developed to provide an analysis framework, by which the estimates of effects for a quantitative trait are statistically orthogonal regardless of the existence of Hardy-Weinberg Equilibrium (HWE) within loci. Ma et al. (2012) recently developed a NOIA model for the gene-environment interaction studies and have shown the advantages of using the model for detecting the true main effects and interactions, compared with the usual functional model. In this project, we propose a novel Bayesian statistical model that combines the Bayesian hierarchical mixture model with the NOIA statistical model and the usual functional model. The proposed Bayesian NOIA model demonstrates more power at detecting the non-null effects with higher marginal posterior probabilities. Also, we review two Bayesian statistical models (Bayesian empirical shrinkage-type estimator and Bayesian model averaging), which were developed for the gene-environment interaction studies. Inspired by these Bayesian models, we develop two novel statistical methods that are able to handle the related problems such as borrowing data from historical studies. The proposed methods are analogous to the methods for the gene-environment interactions on behalf of the success on balancing the statistical efficiency and bias in a unified model. By extensive simulation studies, we compare the operating characteristics of the proposed models with the existing models including the hierarchical meta-analysis model. The results show that the proposed approaches adaptively borrow the historical data in a data-driven way. These novel models may have a broad range of statistical applications in both of genetic/genomic and clinical studies.

Trends in utilization of computer tomography (CT) scans for victims of motor vehicle collision (MVCs) in a level I trauma center from 1996 to 2006

Purpose. To evaluate trends in the utilization of head, abdominal, thoracic and other body regions CTs in the management of victims of MVC at a level I trauma center from 1996 to 2006.^ Method. From the trauma registry, I identified patients involved in MVC's in a level I trauma center and categorized them into three age groups of 13-18, 19-55 and ≥56. I used International Classification of Disease (ICD-9-CM) codes to find the type and number of CTs examinations performed for each patient. I plotted the mean number of CTs per patient against year of admission to find the crude estimate of change in utilization pattern for each type of CT. I used logistic regression to assess whether repetitive CTs (≥ 2) for head, abdomen, thorax and other body regions were associated with age group and year of admission for MVC patients. I adjusted the estimates for gender, ethnicity, insurance status, mechanism and severity of injury, intensive care unit admission status, patient disposition (dead or alive) and year of admission.^ Results. Utilization of head, abdominal, thoracic and other body regions CTs significantly increased over 11-year period. Utilization of head CT was greatest in the 13-18 age group, and increased from 0.58 CT/patient in 1996 to 1.37 CT/patient in 2006. Abdominal CTs were more common in the ≥56+ age group, and increased from 0.33 CT/patient in 1996 to 0.72 CT/patient in 2006. Utilization of thoracic CTs was higher in the 56+ age group, and increased from 0.01 CT/patient in 1996 to 0.42 CT/patient in 2006. Utilization of other CTs did not change materially during the study period for adolescents, adults or older adults. In the multivariable analysis, after adjustment for potential confounders, repetitive head CTs significantly increased in the 13-18 age group (95% CI: 1.29-1.87, p=<0.001) relative to the 19-55 age group. Repetitive thoracic CT use was lower in adolescents (95% CI: 0.22-0.70, p=<0.001) relative to the 19-55 age group.^ Conclusion. There has been a substantial increase in the utilization of head, abdominal, thoracic and other CTs in the management of MVC patients. Future studies need to identify if increased utilization of CTs have resulted in better health outcome for these patients. ^

HPV vaccine coverage in Texas counties: An application of multilevel, small area estimation

Health departments, research institutions, policy-makers, and healthcare providers are often interested in knowing the health status of their clients/constituents. Without the resources, financially or administratively, to go out into the community and conduct health assessments directly, these entities frequently rely on data from population-based surveys to supply the information they need. Unfortunately, these surveys are ill-equipped for the job due to sample size and privacy concerns. Small area estimation (SAE) techniques have excellent potential in such circumstances, but have been underutilized in public health due to lack of awareness and confidence in applying its methods. The goal of this research is to make model-based SAE accessible to a broad readership using clear, example-based learning. Specifically, we applied the principles of multilevel, unit-level SAE to describe the geographic distribution of HPV vaccine coverage among females aged 11-26 in Texas.^ Multilevel (3 level: individual, county, public health region) random-intercept logit models of HPV vaccination (receipt of ≥ 1 dose Gardasil® ) were fit to data from the 2008 Behavioral Risk Factor Surveillance System (outcome and level 1 covariates) and a number of secondary sources (group-level covariates). Sampling weights were scaled (level 1) or constructed (levels 2 & 3), and incorporated at every level. Using the regression coefficients (and standard errors) from the final models, I simulated 10,000 datasets for each regression coefficient from the normal distribution and applied them to the logit model to estimate HPV vaccine coverage in each county and respective demographic subgroup. For simplicity, I only provide coverage estimates (and 95% confidence intervals) for counties.^ County-level coverage among females aged 11-17 varied from 6.8-29.0%. For females aged 18-26, coverage varied from 1.9%-23.8%. Aggregated to the state level, these values translate to indirect state estimates of 15.5% and 11.4%, respectively; both of which fall within the confidence intervals for the direct estimates of HPV vaccine coverage in Texas (Females 11-17: 17.7%, 95% CI: 13.6, 21.9; Females 18-26: 12.0%, 95% CI: 6.2, 17.7).^ Small area estimation has great potential for informing policy, program development and evaluation, and the provision of health services. Harnessing the flexibility of multilevel, unit-level SAE to estimate HPV vaccine coverage among females aged 11-26 in Texas counties, I have provided (1) practical guidance on how to conceptualize and conduct modelbased SAE, (2) a robust framework that can be applied to other health outcomes or geographic levels of aggregation, and (3) HPV vaccine coverage data that may inform the development of health education programs, the provision of health services, the planning of additional research studies, and the creation of local health policies.^

AN ESTIMATION OF IN UTERO DRUG EXPOSURE BY DETERMINANTS OF FREQUENCY AND DISTRIBUTION OF MATERNAL DRUG CONSUMPTION

Data derived from 1,194 gravidas presenting at the observation unit of a city/county hospital between October 11, 1979 through December 7, 1979 were evaluated with respect to the proportion ingesting drugs during pregnancy. The mean age of the mother at the time of the interview was 22.0 years; 43.0 percent were Black; 34.0 percent Latin-American, 21.0 percent White and 2.0 percent other; mean gravida was 2.5 pregnancies; mean parity was 1.0; and mean number of previous abortions was 0.34. Completed interview data was available for 1,119 gravida, corresponding urinalyses for 997 subjects. Ninety and one-tenth percent (90.1 percent) of the subjects reported ingestion of one or more drug preparation(s) (prescription, OTC, or substances used for recreational purposes) during pregnancy with a range of 0 to 11 substances and a mean of 2.7. Dietary supplements (vitamins and minerals) were most frequently reported followed by non-narcotic analgesics. Seventy-six and one tenth percent (76.1 percent) of the population reported consumption of prescription medication, 42.5 percent reported consumption of over-the-counter medications, 45.7 percent reported consumption of a substance for recreational purposes and 4.3 percent reported illicit consumption of a substance. For selected substances, no measurable difference was found between obtaining the information from the interview method or from a urinalysis assay. ^