Comparison of multiple comparison methods for identifying differential gene expression in simulated and real papillary thyroid cancer microarray data

The difficulty of detecting differential gene expression in microarray data has existed for many years. Several correction procedures try to avoid the family-wise error rate in multiple comparison process, including the Bonferroni and Sidak single-step p-value adjustments, Holm's step-down correction method, and Benjamini and Hochberg's false discovery rate (FDR) correction procedure. Each multiple comparison technique has its advantages and weaknesses. We studied each multiple comparison method through numerical studies (simulations) and applied the methods to the real exploratory DNA microarray data, which detect of molecular signatures in papillary thyroid cancer (PTC) patients. According to our results of simulation studies, Benjamini and Hochberg step-up FDR controlling procedure is the best process among these multiple comparison methods and we discovered 1277 potential biomarkers among 54675 probe sets after applying the Benjamini and Hochberg's method to PTC microarray data.^

REGULATION OF TOXIN SYNTHESIS BY CLOSTRIDIUM DIFFICILE

Clostridium difficile is the leading definable cause of nosocomial diarrhea worldwide due to its virulence, multi-drug resistance, spore-forming ability, and environmental persistence. The incidence of C. difficile infection (CDI) has been increasing exponentially in the last decade. Virulent strains of C. difficile produce either toxin A and/or toxin B, which are essential for the pathogenesis of this bacterium. Current methods for diagnosing CDI are mostly qualitative tests that detect the bacterium, the toxins, or the toxin genes. These methods do not differentiate virulent C. difficile strains that produce active toxins from non-virulent strains that do not produce toxins or produce inactive toxins. Based on the knowledge that C. difficile toxins A and B cleave a substrate that is stereochemically similar to the native substrate of the toxins, uridine diphosphoglucose, a quantitative, cost-efficient assay, the Cdifftox activity assay, was developed to measure C. difficile toxin activity. The concept behind the activity assay was modified to develop a novel, rapid, sensitive, and specific assay for C. difficile toxins in the form of a selective and differential agar plate culture medium, the Cdifftox Plate assay (CDPA). This assay combines in a single step the specific identification of C. difficile strains and the detection of active toxin(s). The CDPA was determined to be extremely accurate (99.8% effective) at detecting toxin-producing strains based on the analysis of 528 C. difficile isolates selected from 50 tissue culture cytotoxicity assay-positive clinical stool samples. This new assay advances and improves the culture methodology in that only C. difficile strains will grow on this medium and virulent strains producing active toxins can be differentiated from non-virulent strains. This new method reduces the time and effort required to isolate and confirm toxin-producing C. difficile strains and provides a clinical isolate for antibiotic susceptibility testing and strain typing. The Cdifftox activity assay was used to screen for inhibitors of toxin activity. Physiological levels of the common human conjugated bile salt, taurocholate, was found to inhibit toxin A and B in vitro activities. When co-incubated ex vivo with purified toxin B, taurocholate protected Caco-2 colonic epithelial cells from the damaging effects of the toxin. Furthermore, using a caspase-3 detection assay, taurocholate reduced the extent of toxin B-induced Caco-2 cell apoptosis. These results suggest that bile salts can be effective in protecting the gut epithelium from C. difficile toxin damage, thus, the delivery of physiologic amounts of taurocholate to the colon, where it is normally in low concentration, could be useful in CDI treatment. These findings may help to explain why bile rich small intestine is spared damage in CDI, while the bile salt poor colon is vulnerable in CDI. Toxin synthesis in C. difficile occurs during the stationary phase, but little is known about the regulation of these toxins. It was hypothesized that C. difficile toxin synthesis is regulated by a quorum sensing mechanism. Two lines of evidence supported this hypothesis. First, a small (KDa), diffusible, heat-stable toxin-inducing activity accumulates in the medium of high-density C. difficile cells. This conditioned medium when incubated with low-density log-phase cells causes them to produce toxin early (2-4 hrs instead of 12-16 hrs) and at elevated levels when compared with cells grown in fresh medium. These data suggested that C. difficile cells extracellularly release an inducing molecule during growth that is able to activate toxin synthesis prematurely and demonstrates for the first time that toxin synthesis in C. difficile is regulated by quorum signaling. Second, this toxin-inducing activity was partially purified from high-density stationary-phase culture supernatant fluid by HPLC and confirmed to induce early toxin synthesis, even in C. difficile virulent strains that over-produce the toxins. Mass spectrometry analysis of the purified toxin-inducing fraction from HPLC revealed a cyclic compound with a mass of 655.8 Da. It is anticipated that identification of this toxin-inducing compound will advance our understanding of the mechanism involved in the quorum-dependent regulation of C. difficile toxin synthesis. This finding should lead to the development of even more sensitive tests to diagnose CDI and may lead to the discovery of promising novel therapeutic targets that could be harnessed for the treatment C. difficile infections.

Single nucleotide polymorphism (SNP) selection for genotype-phenotype association studies

With hundreds of single nucleotide polymorphisms (SNPs) in a candidate gene and millions of SNPs across the genome, selecting an informative subset of SNPs to maximize the ability to detect genotype-phenotype association is of great interest and importance. In addition, with a large number of SNPs, analytic methods are needed that allow investigators to control the false positive rate resulting from large numbers of SNP genotype-phenotype analyses. This dissertation uses simulated data to explore methods for selecting SNPs for genotype-phenotype association studies. I examined the pattern of linkage disequilibrium (LD) across a candidate gene region and used this pattern to aid in localizing a disease-influencing mutation. The results indicate that the r2 measure of linkage disequilibrium is preferred over the common D′ measure for use in genotype-phenotype association studies. Using step-wise linear regression, the best predictor of the quantitative trait was not usually the single functional mutation. Rather it was a SNP that was in high linkage disequilibrium with the functional mutation. Next, I compared three strategies for selecting SNPs for application to phenotype association studies: based on measures of linkage disequilibrium, based on a measure of haplotype diversity, and random selection. The results demonstrate that SNPs selected based on maximum haplotype diversity are more informative and yield higher power than randomly selected SNPs or SNPs selected based on low pair-wise LD. The data also indicate that for genes with small contribution to the phenotype, it is more prudent for investigators to increase their sample size than to continuously increase the number of SNPs in order to improve statistical power. When typing large numbers of SNPs, researchers are faced with the challenge of utilizing an appropriate statistical method that controls the type I error rate while maintaining adequate power. We show that an empirical genotype based multi-locus global test that uses permutation testing to investigate the null distribution of the maximum test statistic maintains a desired overall type I error rate while not overly sacrificing statistical power. The results also show that when the penetrance model is simple the multi-locus global test does as well or better than the haplotype analysis. However, for more complex models, haplotype analyses offer advantages. The results of this dissertation will be of utility to human geneticists designing large-scale multi-locus genotype-phenotype association studies. ^