Genetic variations in the PI3K-AKT-mTOR pathway and bladder cancer susceptibility and clinical outcomes

Bladder cancer is the fourth most common cancer in men in the United States. There is compelling evidence supporting that genetic variations contribute to the risk and outcomes of bladder cancer. The PI3K-AKT-mTOR pathway is a major cellular pathway involved in proliferation, invasion, inflammation, tumorigenesis, and drug response. Somatic aberrations of PI3K-AKT-mTOR pathway are frequent events in several cancers including bladder cancer; however, no studies have investigated the role of germline genetic variations in this pathway in bladder cancer. In this project, we used a large case control study to evaluate the associations of a comprehensive catalogue of SNPs in this pathway with bladder cancer risk and outcomes. Three SNPs in RAPTOR were significantly associated with susceptibility: rs11653499 (OR: 1.79, 95%CI: 1.24–2.60), rs7211818 (OR: 2.13, 95%CI: 1.35–3.36), and rs7212142 (OR: 1.57, 95%CI: 1.19–2.07). Two haplotypes constructed from these 3 SNPs were also associated with bladder cancer risk. In combined analysis, a significant trend was observed for increased risk with an increase in the number of unfavorable genotypes (P for trend<0.001). Classification and regression tree analysis identified potential gene-environment interactions between RPS6KA5 rs11653499 and smoking. In superficial bladder cancer, we found that PTEN rs1234219 and rs11202600, TSC1 rs7040593, RAPTOR rs901065, and PIK3R1 rs251404 were significantly associated with recurrence in patients receiving BCG. In muscle invasive and metastatic bladder cancer, AKT2 rs3730050, PIK3R1 rs10515074, and RAPTOR rs9906827 were associated with survival. Survival tree analysis revealed potential gene-gene interactions: patients carrying the unfavorable genotypes of PTEN rs1234219 and TSC1 rs704059 exhibited a 5.24-fold (95% CI: 2.44–11.24) increased risk of recurrence. In combined analysis, with the increasing number of unfavorable genotypes, there was a significant trend of higher risk of recurrence and death (P for trend<0.001) in Cox proportional hazard regression analysis, and shorter event (recurrence and death) free survival in Kaplan-Meier estimates (P log rank<0.001). This study strongly suggests that genetic variations in PI3K-AKT-mTOR pathway play an important role in bladder cancer development. The identified SNPs, if validated in further studies, may become valuable biomarkers in assessing an individual's cancer risk, predicting prognosis and treatment response, and facilitating physicians to make individualized treatment decisions. ^

Genetic analysis of Pitx2 in left -right asymmetry and cardiovascular development

Pitx2, a paired-related homeobox gene that is mutated in human Rieger Syndrome, plays a key role in transferring the early asymmetric signals to individual organs. Pitx2 encodes three isoforms, Pitx2a, Pitx2b and Pitx2c. I found that Pitx2c was the Pitx2 isoform for regulating left-right asymmetry in heart, lung and the predominant isoform in guts. Previous studies suggested that the generation of left-right asymmetry within individual organs is an all or none, random event. Phenotypic analysis of various Pitx2 allelic combinations, that encode graded levels of Pitx2c, reveals an organ-intrinsic mechanism for regulating left-right asymmetric morphogenesis based on differential response to Pitx2c levels. The heart needs low Pitx2c levels, while the lungs and duodenum require higher doses of Pitx2c. In addition, the duodenal rotation is under strict control of Pitx2c activity. Left-right asymmetry development for aortic arch arteries involves complex vascular remodeling. Left-sided expression of Pitx2c in these developing vessels implied its potential function in this process. In order to determine if Pitx2c also can regulate the left-right asymmetry of the aortic arch arteries, a Pitx2c-specific loss of function mutation is generated. Although in wild type mice, the direction of the aortic arch is always oriented toward the left side, the directions of the aortic arches in the mutants were randomized, showing that Pitx2c also determined the left-right asymmetry of these vessels. I have further showed that the cardiac neural crest wasn't involved in this vascular remodeling process. In addition, all mutant embryos had Double Outlet Right Ventricle (DORV), a common congenital heart disease. This study provided insight into the mechanism of Pitx2c-mediated late stages of left-right asymmetry development and identified the roles of Pitx2c in regulation of aortic arch remodeling and heart development. ^

Genetic variants in the mammalian target of rapamycin (mTOR) signaling pathway as predictors of survival and clinical response in women with ovarian cancer

Background. The mTOR pathway is commonly altered in human tumors and promotes cell survival and proliferation. Preliminary evidence suggests this pathway's involvement in chemoresistance to platinum and taxanes, first line therapy for epithelial ovarian cancer. A pathway-based approach was used to identify individual germline single nucleotide polymorphisms (SNPs) and cumulative effects of multiple genetic variants in mTOR pathway genes and their association with clinical outcome in women with ovarian cancer. ^ Methods. The case-series was restricted to 319 non-Hispanic white women with high grade ovarian cancer treated with surgery and platinum-based chemotherapy. 135 SNPs in 20 representative genes in the mTOR pathway were genotyped. Hazard ratios (HRs) for death and Odds ratios (ORs) for failure to respond to primary therapy were estimated for each SNP using the multivariate Cox proportional hazards model and multivariate logistic regression model, respectively, while adjusting for age, stage, histology and treatment sequence. A survival tree analysis of SNPs with a statistically significant association (p<0.05) was performed to identify higher order gene-gene interactions and their association with overall survival. ^ Results. There was no statistically significant difference in survival by tumor histology or treatment regimen. The median survival for the cohort was 48.3 months. Seven SNPs were significantly associated with decreased survival. Compared to those with no unfavorable genotypes, the HR for death increased significantly with the increasing number of unfavorable genotypes and women in the highest risk category had HR of 4.06 (95% CI 2.29–7.21). The survival tree analysis also identified patients with different survival patterns based on their genetic profiles. 13 SNPs on five different genes were found to be significantly associated with a treatment response, defined as no evidence of disease after completion of primary therapy. Rare homozygous genotype of SNP rs6973428 showed a 5.5-fold increased risk compared to the wild type carrying genotypes. In the cumulative effect analysis, the highest risk group (individuals with ≥8 unfavorable genotypes) was significantly less likely to respond to chemotherapy (OR=8.40, 95% CI 3.10–22.75) compared to the low risk group (≤4 unfavorable genotypes). ^ Conclusions. A pathway-based approach can demonstrate cumulative effects of multiple genetic variants on clinical response to chemotherapy and survival. Therapy targeting the mTOR pathway may modify outcome in select patients.^

Detailed analysis of sequence changes occurring during vlsE antigenic variation in the mouse model of Borrelia burgdorferi infection.

Lyme disease Borrelia can infect humans and animals for months to years, despite the presence of an active host immune response. The vls antigenic variation system, which expresses the surface-exposed lipoprotein VlsE, plays a major role in B. burgdorferi immune evasion. Gene conversion between vls silent cassettes and the vlsE expression site occurs at high frequency during mammalian infection, resulting in sequence variation in the VlsE product. In this study, we examined vlsE sequence variation in B. burgdorferi B31 during mouse infection by analyzing 1,399 clones isolated from bladder, heart, joint, ear, and skin tissues of mice infected for 4 to 365 days. The median number of codon changes increased progressively in C3H/HeN mice from 4 to 28 days post infection, and no clones retained the parental vlsE sequence at 28 days. In contrast, the decrease in the number of clones with the parental vlsE sequence and the increase in the number of sequence changes occurred more gradually in severe combined immunodeficiency (SCID) mice. Clones containing a stop codon were isolated, indicating that continuous expression of full-length VlsE is not required for survival in vivo; also, these clones continued to undergo vlsE recombination. Analysis of clones with apparent single recombination events indicated that recombinations into vlsE are nonselective with regard to the silent cassette utilized, as well as the length and location of the recombination event. Sequence changes as small as one base pair were common. Fifteen percent of recovered vlsE variants contained "template-independent" sequence changes, which clustered in the variable regions of vlsE. We hypothesize that the increased frequency and complexity of vlsE sequence changes observed in clones recovered from immunocompetent mice (as compared with SCID mice) is due to rapid clearance of relatively invariant clones by variable region-specific anti-VlsE antibody responses.

Lobe specific Ca2+-calmodulin nano-domain in neuronal spines: a single molecule level analysis.

Calmodulin (CaM) is a ubiquitous Ca(2+) buffer and second messenger that affects cellular function as diverse as cardiac excitability, synaptic plasticity, and gene transcription. In CA1 pyramidal neurons, CaM regulates two opposing Ca(2+)-dependent processes that underlie memory formation: long-term potentiation (LTP) and long-term depression (LTD). Induction of LTP and LTD require activation of Ca(2+)-CaM-dependent enzymes: Ca(2+)/CaM-dependent kinase II (CaMKII) and calcineurin, respectively. Yet, it remains unclear as to how Ca(2+) and CaM produce these two opposing effects, LTP and LTD. CaM binds 4 Ca(2+) ions: two in its N-terminal lobe and two in its C-terminal lobe. Experimental studies have shown that the N- and C-terminal lobes of CaM have different binding kinetics toward Ca(2+) and its downstream targets. This may suggest that each lobe of CaM differentially responds to Ca(2+) signal patterns. Here, we use a novel event-driven particle-based Monte Carlo simulation and statistical point pattern analysis to explore the spatial and temporal dynamics of lobe-specific Ca(2+)-CaM interaction at the single molecule level. We show that the N-lobe of CaM, but not the C-lobe, exhibits a nano-scale domain of activation that is highly sensitive to the location of Ca(2+) channels, and to the microscopic injection rate of Ca(2+) ions. We also demonstrate that Ca(2+) saturation takes place via two different pathways depending on the Ca(2+) injection rate, one dominated by the N-terminal lobe, and the other one by the C-terminal lobe. Taken together, these results suggest that the two lobes of CaM function as distinct Ca(2+) sensors that can differentially transduce Ca(2+) influx to downstream targets. We discuss a possible role of the N-terminal lobe-specific Ca(2+)-CaM nano-domain in CaMKII activation required for the induction of synaptic plasticity.

Detailed analysis of sequence changes occurring during vlsE antigenic variation in the mouse model of Borrelia burgdorferi infection.

Lyme disease Borrelia can infect humans and animals for months to years, despite the presence of an active host immune response. The vls antigenic variation system, which expresses the surface-exposed lipoprotein VlsE, plays a major role in B. burgdorferi immune evasion. Gene conversion between vls silent cassettes and the vlsE expression site occurs at high frequency during mammalian infection, resulting in sequence variation in the VlsE product. In this study, we examined vlsE sequence variation in B. burgdorferi B31 during mouse infection by analyzing 1,399 clones isolated from bladder, heart, joint, ear, and skin tissues of mice infected for 4 to 365 days. The median number of codon changes increased progressively in C3H/HeN mice from 4 to 28 days post infection, and no clones retained the parental vlsE sequence at 28 days. In contrast, the decrease in the number of clones with the parental vlsE sequence and the increase in the number of sequence changes occurred more gradually in severe combined immunodeficiency (SCID) mice. Clones containing a stop codon were isolated, indicating that continuous expression of full-length VlsE is not required for survival in vivo; also, these clones continued to undergo vlsE recombination. Analysis of clones with apparent single recombination events indicated that recombinations into vlsE are nonselective with regard to the silent cassette utilized, as well as the length and location of the recombination event. Sequence changes as small as one base pair were common. Fifteen percent of recovered vlsE variants contained "template-independent" sequence changes, which clustered in the variable regions of vlsE. We hypothesize that the increased frequency and complexity of vlsE sequence changes observed in clones recovered from immunocompetent mice (as compared with SCID mice) is due to rapid clearance of relatively invariant clones by variable region-specific anti-VlsE antibody responses.

A molecular analysis of the conditional defect in MuSVts 110 viral RNA splicing

Cells infected with MuSVts110 express a viral RNA which contains an inherent conditional defect in RNA splicing. It has been shown previously that splicing of the MuSVts110 primary transcript is essential to morphological transformation of 6m2 cells in vitro. A growth temperature of 33$\sp\circ$C is permissive for viral RNA splicing,and, consequently, 6m2 cells appear morphologically transformed at this temperature. However, 6m2 cells appear phenotypically normal when incubated at 39$\sp\circ$C, the non-permissive temperature for viral RNA splicing.^ After a shift from 39$\sp\circ$C to 33$\sp\circ$C, the coordinate splicing of previously synthesized and newly transcribed MuSVts110 RNA was achieved. By S1 nuclease analysis of total RNA isolated at various times, 5$\sp\prime$ splice site cleavage of the MuSVts110 transcript appeared to occur 60 minutes after the shift to 33$\sp\circ$C, and 30 minutes prior to detectable exon ligation. In addition, consistent with the permissive temperatures and the kinetic timeframe of viral RNA splicing after a shift to 33$\sp\circ$C, four temperature sensitive blockades to primer extension were identified 26-75 bases upstream of the 3$\sp\prime$ splice site. These blockades likely reflect four branchpoint sequences utilized in the formation of MuSVts110 lariat splicing-intermediates.^ The 54-5A4 cell line is a spontaneous revertant of 6m2 cells and appears transformed at all growth temperatures. Primer extension sequence analysis has shown that a five base deletion occurred at the 3$\sp\prime$ splice site in MuSVts110 RNA allowing the expression of a viral transforming protein in 54-5A4 in the absence of RNA splicing, whereas in the parental 6m2 cell line, a splicing event is necessary to generate a similar transforming protein. As a consequence of this deletion, splicing cannot occur and the formation of the four MuSVts110 branched-intermediates were not observed at any temperature in 54-5A4 cells. However, 5$\sp\prime$ splice site cleavage was still detected at 33$\sp\circ$C.^ Finally, we have investigated the role of the 1488 bp deletion which occurred in the generation of MuSVts110 in the activation of temperature sensitive viral RNA splicing. This deletion appears solely responsible for splice site activation. Whether intron size is the crucial factor in MuSVts110 RNA splicing or whether inhibitory sequences were removed by the deletion is currently unknown. (Abstract shortened with permission of author.) ^

Analysis of pp60$\sp{{\rm c}{-}src}$ and pp62$\sp{{\rm c}{-}yes}$ intracellular protein tyrosine kinase function in colon tumor cell growth regulation using herbimycin A, a tyrosine kinase specific inhibitor

Two approaches were utilized to investigate the role of pp60c-src activation in growth control of model colon tumor cell lines. The first approach involved analysis of pp60c-src activity in response to growth factor treatment to determine if transient activation of the protein was associated with ligand induced mitogenic signal transduction as occurs in non-colonic cell types. Activation of pp60c-src was detected using colon tumor cell lysates after treatment with platelet derived growth factor (PDGF). Activation of pp60c-src was also detected in response to epidermal growth factor (EGF) treatment using cellular lysates and intact cells. In contrast, down-regulation of purified pp60c-src occurred after incubation with EGF-treated EGFr immune complexes in vitro suggesting additional cellular events were potentially required for the stimulatory response observed in intact cells. The results demonstrated activation of pp60c-src in colon tumor cells in response to PDGF and EGF which is consistent with the role of the protein in mitogenic signal transduction in non-colonic cell types.^ The second approach used to study the role of pp60c-src activation in colonic cell growth control focused on analysis of the role of constitutive activation of the protein, which occurs in approximately 80% of colon tumors and cell lines, in growth control. These studies involved analysis of the effects of the tyrosine kinase specific inhibitor Herbimycin A (HA) on monolayer growth and pp60c-src enzymatic activity using model colon tumor cell lines. HA induced dose-dependent growth inhibition of all colon tumor cell lines examined possessing elevated pp60c-src activity. In HT29 cells the dose-dependent growth inhibition induced by HA correlated with dose-dependent pp60c-src inactivation. Inactivation of pp60c-src was shown to be an early event in response to treatment with HA which preceded induction of HT29 colon tumor cell growth inhibition. The growth effects of HA towards the colon tumor cells examined did not appear to be associated with induction of differentiation or a cytotoxic mechanism of action as changes in morphology were not detected in treated cells and growth inhibition (and pp60c-src inactivation) were reversible upon release from treatment with the compound. The results suggested the constitutive activation of pp60c-src functioned as a proliferative signal in colon tumor cells. Correlation between pp60c-src inactivation and growth inhibition was also observed using HA chemical derivatives confirming the role of tyrosine kinase inactivation by these compounds in inhibition of mitogenic signalling. In contrast, in AS15 cells possessing specific antisense mRNA mediated inactivation of pp60c-src, HA-induced inactivation of the related pp62c-yes tyrosine kinase, which is also activated during colon tumor progression, was not associated with induction of monolayer growth inhibition. These results suggested a function for the constitutively activated pp62c-yes protein in colon tumor cell proliferation which was different from that of activated pp60c-src. (Abstract shortened by UMI.) ^

Demographic analysis of Certified Industrial Hygienists (CIH) with the use of survival analysis techniques

Objective. The goal of this study is to characterize the current workforce of CIHs, the lengths of professional practice careers of the past and current CIHs.^ Methods. This is a secondary data analysis of data compiled from all of the nearly 50 annual roster listings of the American Board of Industrial Hygiene (ABIH) for Certified Industrial Hygienists active in each year since 1960. Survival analysis was performed as a technique to measure the primary outcome of interest. The technique which was involved in this study was the Kaplan-Meier method for estimating the survival function.^ Study subjects: The population to be studied is all Certified Industrial Hygienists (CIHs). A CIH is defined by the ABIH as an individual who has achieved the minimum requirements for education, working experience and through examination, has demonstrated a minimum level of knowledge and competency in the prevention of occupational illnesses. ^ Results. A Cox-proportional hazards model analysis was performed by different start-time cohorts of CIHs. In this model we chose cohort 1 as the reference cohort. The estimated relative risk of the event (defined as retirement, or absent from 5 consecutive years of listing) occurred for CIHs for cohorts 2,3,4,5 relative to cohort 1 is 0.385, 0.214, 0.234, 0.299 relatively. The result show that cohort 2 (CIHs issued from 1970-1980) has the lowest hazard ratio which indicates the lowest retirement rate.^ Conclusion. The manpower of CIHs (still actively practicing up to the end of 2009) increased tremendously starting in 1980 and grew into a plateau in recent decades. This indicates that the supply and demand of the profession may have reached equilibrium. More demographic information and variables are needed to actually predict the future number of CIHs needed. ^

A descriptive and exploratory analysis of occupational injuries at a chemical manufacturing facility

A retrospective study of 1353 occupational injuries occurring at a chemical manufacturing facility in Houston, Texas from January, 1982 through May, 1988 was performed to investigate the etiology of the occupational injury process. Injury incidence rates were calculated for various sub-populations of workers to determine differences in the risk of injury for various groups. Linear modeling techniques were used to determine the association between certain collected independent variables and severity of an injury event. Finally, two sub-groups of the worker population, shiftworkers and injury recidivists, were examined. An injury recidivist as defined is any worker experiencing one or more injury per year. Overall, female shiftworkers evidenced the highest average injury incidence rate compared to all other worker groups analyzed. Although the female shiftworkers were younger and less experienced, the etiology of their increased risk of injury remains unclear, although the rigors of performing shiftwork itself or ergonomic factors are suspect. In general, females were injured more frequently than males, but they did not incur more severe injuries. For all workers, many injuries were caused by erroneous or foregone training, and risk taking behaviors. Injuries of these types are avoidable. The distribution of injuries by severity level was bimodal; either injuries were of minor or major severity with only a small number of cases falling in between. Of the variables collected, only the type of injury incurred and the worker's titlecode were statistically significantly associated with injury severity. Shiftworkers did not sustain more severe injuries than other worker groups. Injury to shiftworkers varied as a 24-hour pattern; the greatest number occurred between 1200-1230 hours, (p = 0.002) by Cosinor analysis. Recidivists made up 3.3% of the population (23 males and 10 females), yet suffered 17.8% of the injuries. Although past research suggests that injury recidivism is a random statistical event, analysis of the data by logistic regression implicates gender, area worked, age and job titlecode as being statistically significantly related to injury recidivism at this facility. ^

AN EPIDEMIOLOGIC DESCRIPTION AND ANALYSIS OF LOST WORKDAY OCCUPATIONAL BURN INJURIES

Background research consisted of a hospital case series of all adult burn patients (n = 162) admitted to John Sealy Hospital's burn unit from January 1978 to June 1979. Comparisons between occupationally and nonoccupationally burned adults demonstrated that occupationally burned adults were significantly more likely to have been active in the burn injury event and to have changed jobs during the prior year. They were significantly less likely to have physical or mental problems which contributed to sustaining the burn injury. Comparisons between occupational and nonoccupational burn injury events concluded that occupational burn injury events were significantly more likely to involve multiple sources of energy, sparks as the source of ignition and gases as the source of combustion. Other salient characteristics of occupational burn injuries indicated that subsequent research should focus upon lost workday occupational burns and other injuries sustained by blue-collar petrochemical workers employed in Galveston County, Texas.^ Subsequent research consisted of a historical cohort study of occupational injuries sustained in 1979 by a cohort of blue-collar petrochemical workers (n = 1771) who belonged to O.C.A.W. Local 4-449 in Texas City, Texas. Specific cohort injury rates included 15.08 occupational injuries per 100 person work-years, 11.98 lost workday occupational injuries per 100 person work-years, and 1.64 lost workday occupational burn injuries per 100 person work-years. Salient results from this study indicate that burn injuries are a very important type (in terms of both frequency and severity) of occupational injury sustained by blue-collar petrochemical workers, pipefitters are at greatest risk of lost workday injuries and lost workday burn injuries, company-specific experiences are comparable for lost workday occupational injuries, differences among company-specific nonlost workday occupational injury experiences may not be "safety-related", and minimal job-specific experience may not place employees at greater risk of lost workday burn injuries.^

Kawasaki syndrome: Cluster analysis of incident cases presenting to a pediatric tertiary care hospital in Texas from January 1, 2005 to December 31, 2009

Background: Despite almost 40 years of research into the etiology of Kawasaki Syndrome (KS), there is little research published on spatial and temporal clustering of KS cases. Previous analysis has found significant spatial and temporal clustering of cases, therefore cluster analyses were performed to substantiate these findings and provide insight into incident KS cases discharged from a pediatric tertiary care hospital. Identifying clusters from a single institution would allow for prospective analysis of risk factors and potential exposures for further insight into KS etiology. ^ Methods: A retrospective study was carried out to examine the epidemiology and distribution of patients presenting to Texas Children’s Hospital in Houston, Texas, with a diagnosis of Acute Febrile Mucocutaneous Lymph Node Syndrome (MCLS) upon discharge from January 1, 2005 to December 31, 2009. Spatial, temporal, and space-time cluster analyses were performed using the Bernoulli model with case and control event data. ^ Results: 397 of 102,761 total patients admitted to Texas Children’s Hospital had a principal or secondary diagnosis of Acute Febrile MCLS upon over the 5 year period. Demographic data for KS cases remained consistent with known disease epidemiology. Spatial, temporal, and space-time analyses of clustering using the Bernoulli model demonstrated no statistically significant clusters. ^ Discussion: Despite previous findings of spatial-temporal clustering of KS cases, there were no significant clusters of KS cases discharged from a single institution. This implicates the need for an expanded approach to conducting spatial-temporal cluster analysis and KS surveillance given the limitations of evaluating data from a single institution.^

Mixture modeling for joint analysis of survival, discrete, and continuous data

Mixture modeling is commonly used to model categorical latent variables that represent subpopulations in which population membership is unknown but can be inferred from the data. In relatively recent years, the potential of finite mixture models has been applied in time-to-event data. However, the commonly used survival mixture model assumes that the effects of the covariates involved in failure times differ across latent classes, but the covariate distribution is homogeneous. The aim of this dissertation is to develop a method to examine time-to-event data in the presence of unobserved heterogeneity under a framework of mixture modeling. A joint model is developed to incorporate the latent survival trajectory along with the observed information for the joint analysis of a time-to-event variable, its discrete and continuous covariates, and a latent class variable. It is assumed that the effects of covariates on survival times and the distribution of covariates vary across different latent classes. The unobservable survival trajectories are identified through estimating the probability that a subject belongs to a particular class based on observed information. We applied this method to a Hodgkin lymphoma study with long-term follow-up and observed four distinct latent classes in terms of long-term survival and distributions of prognostic factors. Our results from simulation studies and from the Hodgkin lymphoma study demonstrated the superiority of our joint model compared with the conventional survival model. This flexible inference method provides more accurate estimation and accommodates unobservable heterogeneity among individuals while taking involved interactions between covariates into consideration.^

Molecular and evolutionary genetics of the X -linked visual pigment genes in humans and New World monkeys

Normal humans have one red and at least one green visual pigment genes. These genes are tightly linked as tandem repeats on the X chromosome and each of them has six exons. There is only one X-linked visual pigment gene in New World monkeys (NWMs) but the locus has three polymorphic alleles encoding red, yellow and green visual pigments, respectively. The spectral properties of the squirrel monkey and the marmoset (both NWMs) have been studied and partial sequences of the three alleles are available. To study the evolutionary history of these X-linked opsin genes in humans and NWMs, coding and intron sequences of the three squirrel monkey alleles and the three marmoset alleles were amplified by PCR followed by subcloning and sequencing. Introns 2 and 4 of the human red and green pigment genes were also sequenced. The results obtained are as follows: (1) The sequences of introns 2 and 4 of the human red and green opsin genes are significantly more similar between the two genes than are coding sequences, contrary to the usual situation where coding regions are better conserved in evolution than are introns. The high similarities in the two introns are probably due to recent gene conversion events during evolution of the human lineage. (2) Phylogenetic analysis of both intron and exon sequences indicates that the phylogenetic tree of the available primate opsin genes is the same as the species tree. The two human genes were derived from a gene duplication event after the divergence of the human and NWM lineages. The three alleles in each of the two NWM species diverged after the split of the two NWMs but have persisted in the population for at least 5 million years. (3) Allelic gene conversion might have occurred between the three squirrel monkey alleles. (4) A model of additive effect of hydroxyl-bearing amino acids on spectral tuning is proposed by treating some unknown variables as groups. Under the assumption that some residues have no effect, it is found that at least five amino acid residues, at positions 178 (3 nm), 180 (5 nm), 230 ($-$4 nm), 277 (9 nm) and 285 (13 nm), have linear spectral tuning effects. (5) Adaptive evolution of the opsin genes to different spectral peaks was observed at four residues that are important for spectral tuning. ^

Analysis of interval -censored events in hypertension studies: Evaluation of treatment of coronary heart disease

Many statistical studies feature data with both exact-time and interval-censored events. While a number of methods currently exist to handle interval-censored events and multivariate exact-time events separately, few techniques exist to deal with their combination. This thesis develops a theoretical framework for analyzing a multivariate endpoint comprised of a single interval-censored event plus an arbitrary number of exact-time events. The approach fuses the exact-time events, modeled using the marginal method of Wei, Lin, and Weissfeld, with a piecewise-exponential interval-censored component. The resulting model incorporates more of the information in the data and also removes some of the biases associated with the exclusion of interval-censored events. A simulation study demonstrates that our approach produces reliable estimates for the model parameters and their variance-covariance matrix. As a real-world data example, we apply this technique to the Systolic Hypertension in the Elderly Program (SHEP) clinical trial, which features three correlated events: clinical non-fatal myocardial infarction, fatal myocardial infarction (two exact-time events), and silent myocardial infarction (one interval-censored event). ^