The population health performance of medical students at two medical schools

Public Health and medicine are complimentary disciplines dedicated to the health and well-being of humankind. Worldwide, medical school accreditation bodies require the inclusion of population health in medical education. In 2003, the Institutes of Medicine (IOM) recommended that all medical students receive basic public health training in population-based prevention. The purpose of this study was to (1) examine the public health clinical performance of third-year medical students at two independent medical schools, (2) compare the public health clinical practice performance of the schools, and (3) identify underlying predictors of high and low public health clinical performance at one of the medical schools. ^ This study is unique in its analysis and report of observed medical student public health clinical practices. The cohort consisted of 751 third-year medical students who completed a required clinical performance exam using trained standardized patients. Medical student performance scores on 24 consensus public health items derived from nine patient cases were analyzed.^ The analysis showed nearly identical results for both medical schools at the 60%, 65%, and 70% pass rate. Students performed poorly on items associated with prevention, behavioral science, and surveillance. Factors associated with high student performance included being from an underrepresented minority, matching to a primary care residency, and high class ranking. A review of medical school curriculum at both schools revealed a lack of training in four public health domains. Nationally, 32% of medical students reported inadequate training in public health in the year 2006.^ These findings suggest more dedicated teaching time for public health domains is needed at the medical schools represented in this study. Finally, more research is needed to assess attainment of public health knowledge and skills for medical students nationwide if we are to meet the recommendations of the IOM. ^

Characterization of the ebp(fm) pilus-encoding operon of Enterococcus faecium and its role in biofilm formation and virulence in a murine model of urinary tract infection.

We recently identified 15 genes encoding putative surface proteins with features of MSCRAMMs and/or pili in the Enterococcus faecium TX0016 (DO) genome, including four predicted pilus-encoding gene clusters; we also demonstrated that one of these, ebpABC(fm), is transcribed as an operon, that its putative major pilus subunit, EbpC(fm) (also called pilB), is polymerized into high molecular weight complexes, and that it is enriched among clinical E. faecium isolates. Here, we created a deletion of the ebpABC(fm) operon in an endocarditis-derived E. faecium strain (TX82) and showed, by a combination of whole-cell ELISA, flow cytometry, immunoblot and immunogold electron microscopy, that this deletion abolished EbpC(fm) expression and eliminated EbpC(fm)-containing pili from the cell surface. However, transcription of the downstream sortase, bps(fm), was not affected. Importantly, the ebpABC(fm) deletion resulted in significantly reduced biofilm formation (p < 0.0001) and initial adherence (p < 0.0001) versus the wild-type; both were restored by complementing ebpABC(fm) in trans, which also restored cell surface expression of EbpC(fm) and pilus production. Furthermore, the deletion mutant was significantly attenuated in two independent mixed infection mouse urinary tract experiments, i.e., outnumbered by the wild-type in kidneys (p = 0.0003 and < 0.0001, respectively) and urinary bladders (p = 0.0003 and = 0.002). In conclusion, we have shown that the ebpABC(fm) locus encodes pili on the E. faecium TX82 cell surface and provide the first evidence that pili of this emerging pathogen are important for its ability to form biofilm and to cause infection in an ascending UTI model.

Dominance of the proximal coordinate frame in determining the locations of hippocampal place cell activity during navigation.

The place-specific activity of hippocampal cells provides downstream structures with information regarding an animal's position within an environment and, perhaps, the location of goals within that environment. In rodents, recent research has suggested that distal cues primarily set the orientation of the spatial representation, whereas the boundaries of the behavioral apparatus determine the locations of place activity. The current study was designed to address possible biases in some previous research that may have minimized the likelihood of observing place activity bound to distal cues. Hippocampal single-unit activity was recorded from six freely moving rats as they were trained to perform a tone-initiated place-preference task on an open-field platform. To investigate whether place activity was bound to the room- or platform-based coordinate frame (or both), the platform was translated within the room at an "early" and at a "late" phase of task acquisition (Shift 1 and Shift 2). At both time points, CA1 and CA3 place cells demonstrated room-associated and/or platform-associated activity, or remapped in response to the platform shift. Shift 1 revealed place activity that reflected an interaction between a dominant platform-based (proximal) coordinate frame and a weaker room-based (distal) frame because many CA1 and CA3 place fields shifted to a location intermediate to the two reference frames. Shift 2 resulted in place activity that became more strongly bound to either the platform- or room-based coordinate frame, suggesting the emergence of two independent spatial frames of reference (with many more cells participating in platform-based than in room-based representations).

DUAL ENERGY COMPUTED TOMOGRAPHY FOR PROTON THERAPY TREATMENT PLANNING

Proton radiation therapy is gaining popularity because of the unique characteristics of its dose distribution, e.g., high dose-gradient at the distal end of the percentage-depth-dose curve (known as the Bragg peak). The high dose-gradient offers the possibility of delivering high dose to the target while still sparing critical organs distal to the target. However, the high dose-gradient is a double-edged sword: a small shift of the highly conformal high-dose area can cause the target to be substantially under-dosed or the critical organs to be substantially over-dosed. Because of that, large margins are required in treatment planning to ensure adequate dose coverage of the target, which prevents us from realizing the full potential of proton beams. Therefore, it is critical to reduce uncertainties in the proton radiation therapy. One major uncertainty in a proton treatment is the range uncertainty related to the estimation of proton stopping power ratio (SPR) distribution inside a patient. The SPR distribution inside a patient is required to account for tissue heterogeneities when calculating dose distribution inside the patient. In current clinical practice, the SPR distribution inside a patient is estimated from the patient’s treatment planning computed tomography (CT) images based on the CT number-to-SPR calibration curve. The SPR derived from a single CT number carries large uncertainties in the presence of human tissue composition variations, which is the major drawback of the current SPR estimation method. We propose to solve this problem by using dual energy CT (DECT) and hypothesize that the range uncertainty can be reduced by a factor of two from currently used value of 3.5%. A MATLAB program was developed to calculate the electron density ratio (EDR) and effective atomic number (EAN) from two CT measurements of the same object. An empirical relationship was discovered between mean excitation energies and EANs existing in human body tissues. With the MATLAB program and the empirical relationship, a DECT-based method was successfully developed to derive SPRs for human body tissues (the DECT method). The DECT method is more robust against the uncertainties in human tissues compositions than the current single-CT-based method, because the DECT method incorporated both density and elemental composition information in the SPR estimation. Furthermore, we studied practical limitations of the DECT method. We found that the accuracy of the DECT method using conventional kV-kV x-ray pair is susceptible to CT number variations, which compromises the theoretical advantage of the DECT method. Our solution to this problem is to use a different x-ray pair for the DECT. The accuracy of the DECT method using different combinations of x-ray energies, i.e., the kV-kV, kV-MV and MV-MV pair, was compared using the measured imaging uncertainties for each case. The kV-MV DECT was found to be the most robust against CT number variations. In addition, we studied how uncertainties propagate through the DECT calculation, and found general principles of selecting x-ray pairs for the DECT method to minimize its sensitivity to CT number variations. The uncertainties in SPRs estimated using the kV-MV DECT were analyzed further and compared to those using the stoichiometric method. The uncertainties in SPR estimation can be divided into five categories according to their origins: the inherent uncertainty, the DECT modeling uncertainty, the CT imaging uncertainty, the uncertainty in the mean excitation energy, and SPR variation with proton energy. Additionally, human body tissues were divided into three tissue groups – low density (lung) tissues, soft tissues and bone tissues. The uncertainties were estimated separately because their uncertainties were different under each condition. An estimate of the composite range uncertainty (2s) was determined for three tumor sites – prostate, lung, and head-and-neck, by combining the uncertainty estimates of all three tissue groups, weighted by their proportions along typical beam path for each treatment site. In conclusion, the DECT method holds theoretical advantages in estimating SPRs for human tissues over the current single-CT-based method. Using existing imaging techniques, the kV-MV DECT approach was capable of reducing the range uncertainty from the currently used value of 3.5% to 1.9%-2.3%, but it is short to reach our original goal of reducing the range uncertainty by a factor of two. The dominant source of uncertainties in the kV-MV DECT was the uncertainties in CT imaging, especially in MV CT imaging. Further reduction in beam hardening effect, the impact of scatter, out-of-field object etc. would reduce the Hounsfeld Unit variations in CT imaging. The kV-MV DECT still has the potential to reduce the range uncertainty further.

Intracellular signalling by the insulin receptor

The insulin receptor transduces insulin's biological signal through the tyrosine kinase present in the receptor's B subunit. The activated insulin receptor kinase then phosphorylates a series of intracellular substrate including insulin receptor substrate 1 (IRS-1), which has been shown to be the pivotal substrate for insulin receptor signal transduction. The phosphorylated tyrosine residues in IRS-1 can bind and activate the downstream effectors, many of which are SH2 domain containing proteins such as phosphotidylinositol 3-kinase, growth factor binding protein 2, and SH2 phosphotyrosine phosphatase 2. Phosphorylated synthetic IRS-1 peptides with the corresponding sequences of the IRS-1 have been shown to associate and activate their respective SH2 domain containing proteins. Another important event happening during insulin binding with the insulin receptor is that the insulin receptor rapidly undergoes internalization. However, the insulin receptor signalling and the receptor endocytosis have been studied as two independent processes. The hypothesis of the present thesis is that the insulin receptor endocytosis is involved in insulin receptor signalling and signal termination. The results of the present investigation demonstrate that insulin receptors in the earliest stage of endocytosis contain significantly greater kinase activity towards IRS-1 peptides than the receptors localized at the plasma membrane, indicating that they are potentially more capable of transducing signals. On the other hand, insulin receptors in the middle and late stage of endocytosis lose their kinase activity, suggesting that insulin receptor kinase activity inactivation and signal termination might take place in the late phase of the insulin receptor internalization. In addition, this study also found that the increased insulin receptor kinase activity in the endosomes is related to the tyrosyl phosphorylation of the specific domains of the receptor's $\beta$ subunit. ^

The story of RP10: A positional cloning approach to the identification of an autosomal dominant retinitis pigmentosa gene

Retinitis pigmentosa (RP) is an inherited retinal degenerative disease that is the leading cause of inherited blindness worldwide. Characteristic features of the disease include night blindness, progressive loss of visual fields, and deposition of pigment on the retina in a bone spicule-like pattern. RP is marked by extreme genetic heterogeneity with at least 19 autosomal dominant, autosomal recessive and X-linked loci identified. RP10, which maps to chromosome 7q, was the fifth autosomal dominant RP locus identified, and accounts for the early-onset disease in two independent families. Extensive linkage and haplotype analyses have been performed in these two families which have allowed the assignment of the disease locus to a 5-cM region on chromosome 7q31.3. In collaboration with Dr. Eric Green (National Center for Human Genome Research, National Institutes of Health), a well-characterized physical map of the region was constructed which includes YAC, BAC and cosmid coverage. The entire RP10 critical region resides within a 9-Mb well-characterized YAC contig. These physical maps not only provided the resources to undertake the CAIGES (cDNA amplification for identification of genomic expressed sequences) procedure for identification of retinal candidate genes within the critical region, but also identified a number of candidate genes, including transducin-$\gamma$ and blue cone pigment genes. All candidate genes examined were excluded. In addition, a number of ESTs were mapped within the critical region. EST20241, which was isolated from an eye library, corresponded to the 3$\sp\prime$ region of the ADP-ribosylation factor (ARF) 5 gene. ARF5, with its role in vesicle transport and possible participation in the regulation of the visual transduction pathway, became an extremely interesting candidate gene. Using a primer walking approach, the entire 3.2 kb genomic sequence of the ARF5 gene was generated and developed intronic primers to screen for coding region mutations in affected family members. No mutations were found in the ARF5 gene, however, a number of additional ESTs have been mapped to the critical region, and, as the large-scale sequencing projects get underway, megabases of raw sequence data from the RP10 region are becoming available. These resources will hasten the isolation and characterization of the RP10 gene. ^

Structure-function analyses of {\it PAX6\/} and the molecular bases of aniridia

PAX6 is a transcription activator that regulates eye development in animals ranging from Drosophila to human. The C-terminal region of PAX6 is proline/serine/threonine-rich (PST) and functions as a potent transactivation domain when attached to a heterologous DNA-binding domain of the yeast transcription factor, GAL4. The PST region comprises 152 amino acids encoded by four exons. The transactivation function of the PST region has not been defined and characterized in detail by in vitro mutagenesis. I dissected the PST domain in two independent systems, a heterologous system using a GAL4 DNA-binding site and the native system of PAX6. In both systems, the results show consistently that all four constituent exons of the PST domain are responsible for the transactivation function. The four exon fragments act cooperatively to stimulate transcription, although none of them can function individually as an independent transactivation domain. Combinations of two or more exon fragments can reconstitute substantial transactivation activity when fused to the DNA-binding domain of GAL4, but they surprisingly do not produce much activity in the context of native PAX6 even though the mutant PAX6 proteins are stable and their DNA-binding function remains unaffected. I conclude that the PAX6 protein contains an unusually large transactivation domain that is evolutionarily conserved to a high degree, and that its full transactivation activity relies on the cooperative action of the four exon fragments.^ Most PAX6 mutations detected in patients with aniridia result in truncations of the protein. Some of the truncation mutations occur in the PST region of PAX6, resulting in mutant proteins that retain their DNA-binding ability but have no significant transactivation activity. It is not clear whether such mutants are true loss-of-function or dominant-negative mutants. I show that these mutants are dominant-negative if they are coexpressed with wild-type PAX6 in cultured cells and that the dominant-negative effects result from enhanced DNA-binding ability of these mutants due to removal of the PST domain. These mutants are able to repress the wild-type PAX6 activity not only at target genes with paired domain binding sites but also at target genes with homeodomain binding sites.^ Mutations in the human PAX6 gene produce various phenotypes, including aniridia, Peters' anomaly, autosomal dominant keratitis, and familial foveal dysplasia. The various phenotypes may arise from different mutations in the same gene. To test this theory, I performed a functional analysis of two missense mutations in the paired domain: the R26G mutation reported in a case of Peters' anomaly, and the I87R mutation identified in a patient with aniridia. While both the R26 and the I87 positions are conserved in the paired boxes of all known PAX genes, X-ray crystallography has shown that only R26 makes contact with DNA. I found that the R26G mutant failed to bind a subset of paired domain binding sites but, surprisingly, bound other sites and successfully transactivated promoters containing those sites. In contrast, the I87R mutant had lost the ability to bind DNA at all tested sites and failed to transactivate promoters. My data support the haploinsufficiency hypothesis of aniridia, and the hypothesis that R26G is a hypomorphic allele. ^

Delineating the molecular basis of human genetic diseases: Epigenetic and functional studies

Identifying and characterizing the genes responsible for inherited human diseases will ultimately lead to a more holistic understanding of disease pathogenesis, catalyze new diagnostic and treatment modalities, and provide insights into basic biological processes. This dissertation presents research aimed at delineating the genetic and molecular basis of human diseases through epigenetic and functional studies and can be divided into two independent areas of research. The first area of research describes the development of two high-throughput melting curve based methods to assay DNA methylation, referred to as McMSP and McCOBRA. The goal of this project was to develop DNA methylation methods that can be used to rapidly determine the DNA methylation status at a specific locus in a large number of samples. McMSP and McCOBRA provide several advantages over existing methods, as they are simple, accurate, robust, and high-throughput making them applicable to large-scale DNA methylation studies. McMSP and McCOBRA were then used in an epigenetic study of the complex disease Ankylosing spondylitis (AS). Specifically, I tested the hypothesis that aberrant patterns of DNA methylation in five AS candidate genes contribute to disease susceptibility. While no statistically significant methylation differences were observed between cases and controls, this is the first study to investigate the hypothesis that epigenetic variation contributes to AS susceptibility and therefore provides the conceptual framework for future studies. ^ In the second area of research, I performed experiments to better delimit the function of aryl hydrocarbon receptor-interacting protein-like 1 (AIPL1), which when mutated causes various forms of inherited blindness such as Leber congenital amaurosis. A yeast two-hybrid screen was performed to identify putative AIPL1-interacting proteins. After screening 2 × 106 bovine retinal cDNA library clones, 6 unique putative AIPL1-interacting proteins were identified. While these 6 AIPL1 protein-protein interactions must be confirmed, their identification is an important step in understanding the functional role of AIPL1 within the retina and will provide insight into the molecular mechanisms underlying inherited blindness. ^

Dysregulation of alternative splicing of coagulation factor V results in bleeding disorder, east Texas type and other disorders of hemostasis

The studies completed herein explore different phenotypes related to the genetic defects that predispose individuals to a disruption of normal hemostasis. In the first study, a novel autosomal dominant bleeding disorder, which is characterized by excessive bleeding with trauma or surgery and menorrhagia in affected women, was studied in a large family (16 affected individuals) from east Texas. Affected members had a prolongation of their PT and/or aPTT, but normal clinical coagulation studies. Previous linkage analysis by Kuang et. al. (2001) mapped the defective gene to 1g23-24 (LODmax 7.22), which contains the gene for coagulation factor V (FV). I identified an alteration (A2440G) in the FV gene in exon 13 that segregated with the disease and was not present in 62 controls. Interestingly, this alteration resulted in a 22-fold up-regulation of a novel alternative splicing variant in patients' RNA versus controls. This translated into a similar fold increase in a 250-kDa isoform of FV seen in patients' plasma versus controls. A recombinant of this splicing event exhibited an increased sensitivity to cleavage by activated protein C (APC) that was more striking in the presence of PS. In addition, this novel isoform had increased APC cofactor activity, thus increasing the degradation of FVIIIa. These data indicated that A2440G up-regulates an alternatively spliced transcript of FV, and increases a FV isoform that hinders coagulation as opposed to promoting it like its wild-type counterpart. ^ The second study reports the largest screening to date of African Americans in two independent cohorts for a rare prothrombin variant, C20209T, which is suspected to be associated with thrombotic disease. The Texas Medical Center Genetics Resource (TexGen) Stroke DNA repository revealed 1.67% (Fisher p=0.27) of African American stroke patients were heterozygous for the 20209*T allele. Screening of the Atherosclerosis Risk in Communities Study (ARIC) cohort (n=3470) for the 20209*T allele revealed a population prevalence of 0.58% in individuals of African American descent; however, all associations with thrombotic disease were negative. Analysis of these two independent cohorts revealed that, unlike its neighbor G20210A, the C20209T variant does not increase the risk of thrombotic events in the African American population. ^

The Texas driver responsibility program: A preliminary analysis of the impact on impaired driving and trauma system funding

Objective. In 2003, the State of Texas instituted the Driver Responsibility Program (TDRP), a program consisting of a driving infraction point system coupled with a series of graded fines and annual surcharges for specific traffic violations such as driving while intoxicated (DWI). Approximately half of the revenues generated are earmarked to be disbursed to the state's trauma system to cover uncompensated trauma care costs. This study examined initial program implementation, the impact of trauma system funding, and initial impact on impaired driving knowledge, attitudes and behaviors. A model for targeted media campaigns to improve the program's deterrence effects was developed. ^ Methods. Data from two independent driver survey samples (conducted in 1999 and 2005), department of public safety records, state health department data and a state auditor's report were used to evaluate the program's initial implementation, impact and outcome with respect to drivers' impaired driving knowledge, attitudes and behavior (based on constructs of social cognitive theory) and hospital uncompensated trauma care funding. Survey results were used to develop a regression model of high risk drivers who should be targeted to improve program outcome with respect to deterring impaired driving. ^ Results. Low driver compliance with fee payment (28%) and program implementation problems were associated with lower surcharge revenues in the first two years ($59.5 million versus $525 million predicted). Program revenue distribution to trauma hospitals was associated with a 16% increase in designated trauma centers. Survey data demonstrated that only 28% of drivers are aware of the TDRP and that there has been no initial impact on impaired driving behavior. Logistical regression modeling suggested that target media campaigns highlighting the likelihood of DWI detection by law enforcement and the increased surcharges associated with the TDRP are required to deter impaired driving. ^ Conclusions. Although the TDRP raised nearly $60 million in surcharge revenue for the Texas trauma system over the first two years, this study did not find evidence of a change in impaired driving knowledge, attitudes or behaviors from 1999 to 2005. Further research is required to measure whether the program is associated with decreased alcohol-related traffic fatalities. ^

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 role of environmental tobacco smoke in the development and exacerbation of asthma among children in developing countries: A systematic review

Environmental tobacco smoke (ETS) is an important indoor air pollutant associated with adverse effects on the respiratory health of the general population, especially people with asthma. ETS consists mainly of sidestream smoke from burning cigarettes and a smaller quantity of mainstream smoke which is exhaled by the smoker. At least one out of every three children is frequently exposed to ETS. ^ This paper reviewed the literature for studies on the role of ETS in the development and exacerbation of asthma among children in developing countries, specifically the low and middle income countries from the year 1980 to the present. The databases searched in this systematic review were: Ovid Medline; PubMed (National Library of Medicine); and Cumulative Index to Nursing and Allied Health (CINAHL) (EBSCOhost). Out of a total of 197 articles initially identified, only four studies (two from China, one from Macedonia and one from Brazil) were rated by two independent raters as being of high quality, and were selected for final abstraction, synthesis and evidence weighting. Results from these four studies suggests that, in developing countries, ETS exposure is associated with childhood asthma, and that asthma prevalence increases with an increase in the amount and duration of exposure to ETS. Similarly, exposure to ETS is associated with persistent cough, current night dry cough, and exacerbation of asthma symptoms. ^ Therefore, as is the case in developed nations, there is suggestive evidence in the literature that ETS exposure plays substantial role in the development and/or exacerbation of asthma among children in developing countries. To decrease the likelihood of new asthma development, enhance asthma control, and reduce the rate of medical service utilization in children exposed to ETS, smoking should be eliminated at home and in public places.^

Bisphosphonates for metastatic bone disease in lung cancer

Bisphosphonates have proven effectiveness in preventing skeletal-related events (SREs) in advanced breast cancer, prostate cancer and multiple myeloma. The purpose of this study was to assess efficacy of bisphosphonates in preventing SREs, in controlling pain, and in increasing life expectancy in lung cancer patients with bone metastases.^ We performed an electronic search in MEDLINE, EMBASE, Web of Science, and Cochrane library databases up to April 4, 2010. Hand searching and searching in clinicaltrials.gov were also performed. Two independent reviewers selected all clinical trials that included lung cancer patients with bone metastases treated with bisphosphonates. We excluded articles that involved cancers other than lung, patients without bone metastasis and treatment other than bisphosphonates. Outcome questions answered were efficacy measured as overall pain control, overall improvement in survival and reduction in skeletal-related events or SREs (fracture, cord compression, radiation or surgery to the bone, hypercalcemia of malignancy). The quality of each study was evaluated using the Cochrane Back Review group questionnaire to assess risk of bias (0-worst to 11-best). Data extraction and quality assessments were independently performed by two assessors. Meta-analyses were performed where more than one study with similar outcomes were found.^ We identified eight trials that met our inclusion criteria. Three studies evaluated zoledronic acid, three pamidronate, three clodronate and two ibandronate. Two were placebocontrol trials while two had multi-group comparisons (radiotherapy, radionucleotides, and chemotherapy) and two had different bisphosphonate as active controls. Quality scores ranged from 1-4 out of 11 suggesting high risk of bias. Studies failed to report adequate explanation of randomization procedures, concealment of randomization and blinding. Metaanalysis showed that patients treated with zoledronic acid alone had lower rates of developing SREs compared to placebo at 21 months (RR=0.80, 95% CI=0.66-0.97, p=0.02). Meta-analyses also showed increased pain control when a bisphosphonate was added to the existing treatment modality like chemotherapy or radiation (RR=1.17, 95% CI=1.03-1.34, p=0.02). However, pain control was not statistically significantly different among various bisphosphonates when other treatment modalities were not present. Despite improvement in SRE and pain control, bisphosphonates failed to show improvement in overall survival (Difference in means=109.1 days, 95% CI= -51.52 – 269.71, p=0.183).^ Adding biphosphonates to standard care improved pain control and reduced SREs. Biphosphonates did not improve overall survival. Further larger studies with higher quality are required to stengthen the evidence.^ Keywords/MeSH terms Bisphosphonates/diphosphonates: generic, chemical and trade names.^

A hybrid method in combining treatment effects from matched and unmatched studies

In the biomedical studies, the general data structures have been the matched (paired) and unmatched designs. Recently, many researchers are interested in Meta-Analysis to obtain a better understanding from several clinical data of a medical treatment. The hybrid design, which is combined two data structures, may create the fundamental question for statistical methods and the challenges for statistical inferences. The applied methods are depending on the underlying distribution. If the outcomes are normally distributed, we would use the classic paired and two independent sample T-tests on the matched and unmatched cases. If not, we can apply Wilcoxon signed rank and rank sum test on each case. ^ To assess an overall treatment effect on a hybrid design, we can apply the inverse variance weight method used in Meta-Analysis. On the nonparametric case, we can use a test statistic which is combined on two Wilcoxon test statistics. However, these two test statistics are not in same scale. We propose the Hybrid Test Statistic based on the Hodges-Lehmann estimates of the treatment effects, which are medians in the same scale.^ To compare the proposed method, we use the classic meta-analysis T-test statistic on the combined the estimates of the treatment effects from two T-test statistics. Theoretically, the efficiency of two unbiased estimators of a parameter is the ratio of their variances. With the concept of Asymptotic Relative Efficiency (ARE) developed by Pitman, we show ARE of the hybrid test statistic relative to classic meta-analysis T-test statistic using the Hodges-Lemann estimators associated with two test statistics.^ From several simulation studies, we calculate the empirical type I error rate and power of the test statistics. The proposed statistic would provide effective tool to evaluate and understand the treatment effect in various public health studies as well as clinical trials.^

Gene-gene interaction and gene-pathway analysis of genome-wide association for schizophrenia

Schizophrenia (SZ) is a complex disorder with high heritability and variable phenotypes that has limited success in finding causal genes associated with the disease development. Pathway-based analysis is an effective approach in investigating the molecular mechanism of susceptible genes associated with complex diseases. The etiology of complex diseases could be a network of genetic factors and within the genes, interaction may occur. In this work we argue that some genes might be of small effect that by itself are neither sufficient nor necessary to cause the disease however, their effect may induce slight changes to the gene expression or affect the protein function, therefore, analyzing the gene-gene interaction mechanism within the disease pathway would play crucial role in dissecting the genetic architecture of complex diseases, making the pathway-based analysis a complementary approach to GWAS technique. ^ In this study, we implemented three novel linkage disequilibrium based statistics, the linear combination, the quadratic, and the decorrelation test statistics, to investigate the interaction between linked and unlinked genes in two independent case-control GWAS datasets for SZ including participants of European (EA) and African (AA) ancestries. The EA population included 1,173 cases and 1,378 controls with 729,454 genotyped SNPs, while the AA population included 219 cases and 288 controls with 845,814 genotyped SNPs. We identified 17,186 interacting gene-sets at significant level in EA dataset, and 12,691 gene-sets in AA dataset using the gene-gene interaction method. We also identified 18,846 genes in EA dataset and 19,431 genes in AA dataset that were in the disease pathways. However, few genes were reported of significant association to SZ. ^ Our research determined the pathways characteristics for schizophrenia through the gene-gene interaction and gene-pathway based approaches. Our findings suggest insightful inferences of our methods in studying the molecular mechanisms of common complex diseases.^