Genetic association studies under the population stratification, family pedigree and application to genome-wide association studies

This dissertation has three separate parts: the first part deals with the general pedigree association testing incorporating continuous covariates; the second part deals with the association tests under population stratification using the conditional likelihood tests; the third part deals with the genome-wide association studies based on the real rheumatoid arthritis (RA) disease data sets from Genetic Analysis Workshop 16 (GAW16) problem 1. Many statistical tests are developed to test the linkage and association using either case-control status or phenotype covariates for family data structure, separately. Those univariate analyses might not use all the information coming from the family members in practical studies. On the other hand, the human complex disease do not have a clear inheritance pattern, there might exist the gene interactions or act independently. In part I, the new proposed approach MPDT is focused on how to use both the case control information as well as the phenotype covariates. This approach can be applied to detect multiple marker effects. Based on the two existing popular statistics in family studies for case-control and quantitative traits respectively, the new approach could be used in the simple family structure data set as well as general pedigree structure. The combined statistics are calculated using the two statistics; A permutation procedure is applied for assessing the p-value with adjustment from the Bonferroni for the multiple markers. We use simulation studies to evaluate the type I error rates and the powers of the proposed approach. Our results show that the combined test using both case-control information and phenotype covariates not only has the correct type I error rates but also is more powerful than the other existing methods. For multiple marker interactions, our proposed method is also very powerful. Selective genotyping is an economical strategy in detecting and mapping quantitative trait loci in the genetic dissection of complex disease. When the samples arise from different ethnic groups or an admixture population, all the existing selective genotyping methods may result in spurious association due to different ancestry distributions. The problem can be more serious when the sample size is large, a general requirement to obtain sufficient power to detect modest genetic effects for most complex traits. In part II, I describe a useful strategy in selective genotyping while population stratification is present. Our procedure used a principal component based approach to eliminate any effect of population stratification. The paper evaluates the performance of our procedure using both simulated data from an early study data sets and also the HapMap data sets in a variety of population admixture models generated from empirical data. There are one binary trait and two continuous traits in the rheumatoid arthritis dataset of Problem 1 in the Genetic Analysis Workshop 16 (GAW16): RA status, AntiCCP and IgM. To allow multiple traits, we suggest a set of SNP-level F statistics by the concept of multiple-correlation to measure the genetic association between multiple trait values and SNP-specific genotypic scores and obtain their null distributions. Hereby, we perform 6 genome-wide association analyses using the novel one- and two-stage approaches which are based on single, double and triple traits. Incorporating all these 6 analyses, we successfully validate the SNPs which have been identified to be responsible for rheumatoid arthritis in the literature and detect more disease susceptibility SNPs for follow-up studies in the future. Except for chromosome 13 and 18, each of the others is found to harbour susceptible genetic regions for rheumatoid arthritis or related diseases, i.e., lupus erythematosus. This topic is discussed in part III.

Characterization of thermal and mechanical properties of polypropylene-based composites for fuel cell bipolar plates and development of educational tools in hydrogen and fuel cell technologies

In this project we developed conductive thermoplastic resins by adding varying amounts of three different carbon fillers: carbon black (CB), synthetic graphite (SG) and multi-walled carbon nanotubes (CNT) to a polypropylene matrix for application as fuel cell bipolar plates. This component of fuel cells provides mechanical support to the stack, circulates the gases that participate in the electrochemical reaction within the fuel cell and allows for removal of the excess heat from the system. The materials fabricated in this work were tested to determine their mechanical and thermal properties. These materials were produced by adding varying amounts of single carbon fillers to a polypropylene matrix (2.5 to 15 wt.% Ketjenblack EC-600 JD carbon black, 10 to 80 wt.% Asbury Carbon's Thermocarb TC-300 synthetic graphite, and 2.5 to 15 wt.% of Hyperion Catalysis International's FIBRILTM multi-walled carbon nanotubes) In addition, composite materials containing combinations of these three fillers were produced. The thermal conductivity results showed an increase in both through-plane and in-plane thermal conductivities, with the largest increase observed for synthetic graphite. The Department of Energy (DOE) had previously set a thermal conductivity goal of 20 W/m·K, which was surpassed by formulations containing 75 wt.% and 80 wt.% SG, yielding in-plane thermal conductivity values of 24.4 W/m·K and 33.6 W/m·K, respectively. In addition, composites containing 2.5 wt.% CB, 65 wt.% SG, and 6 wt.% CNT in PP had an in–plane thermal conductivity of 37 W/m·K. Flexural and tensile tests were conducted. All composite formulations exceeded the flexural strength target of 25 MPa set by DOE. The tensile and flexural modulus of the composites increased with higher concentration of carbon fillers. Carbon black and synthetic graphite caused a decrease in the tensile and flexural strengths of the composites. However, carbon nanotubes increased the composite tensile and flexural strengths. Mathematical models were applied to estimate through-plane and in-plane thermal conductivities of single and multiple filler formulations, and tensile modulus of single-filler formulations. For thermal conductivity, Nielsen's model yielded accurate thermal conductivity values when compared to experimental results obtained through the Flash method. For prediction of tensile modulus Nielsen's model yielded the smallest error between the predicted and experimental values. The second part of this project consisted of the development of a curriculum in Fuel Cell and Hydrogen Technologies to address different educational barriers identified by the Department of Energy. By the creation of new courses and enterprise programs in the areas of fuel cells and the use of hydrogen as an energy carrier, we introduced engineering students to the new technologies, policies and challenges present with this alternative energy. Feedback provided by students participating in these courses and enterprise programs indicate positive acceptance of the different educational tools. Results obtained from a survey applied to students after participating in these courses showed an increase in the knowledge and awareness of energy fundamentals, which indicates the modules developed in this project are effective in introducing students to alternative energy sources.

Development of a physically-based method for delineation of hydrologically homogeneous regions and flood quantile estimation in ungauged basins via the index flood method

Regional flood frequency techniques are commonly used to estimate flood quantiles when flood data is unavailable or the record length at an individual gauging station is insufficient for reliable analyses. These methods compensate for limited or unavailable data by pooling data from nearby gauged sites. This requires the delineation of hydrologically homogeneous regions in which the flood regime is sufficiently similar to allow the spatial transfer of information. It is generally accepted that hydrologic similarity results from similar physiographic characteristics, and thus these characteristics can be used to delineate regions and classify ungauged sites. However, as currently practiced, the delineation is highly subjective and dependent on the similarity measures and classification techniques employed. A standardized procedure for delineation of hydrologically homogeneous regions is presented herein. Key aspects are a new statistical metric to identify physically discordant sites, and the identification of an appropriate set of physically based measures of extreme hydrological similarity. A combination of multivariate statistical techniques applied to multiple flood statistics and basin characteristics for gauging stations in the Southeastern U.S. revealed that basin slope, elevation, and soil drainage largely determine the extreme hydrological behavior of a watershed. Use of these characteristics as similarity measures in the standardized approach for region delineation yields regions which are more homogeneous and more efficient for quantile estimation at ungauged sites than those delineated using alternative physically-based procedures typically employed in practice. The proposed methods and key physical characteristics are also shown to be efficient for region delineation and quantile development in alternative areas composed of watersheds with statistically different physical composition. In addition, the use of aggregated values of key watershed characteristics was found to be sufficient for the regionalization of flood data; the added time and computational effort required to derive spatially distributed watershed variables does not increase the accuracy of quantile estimators for ungauged sites. This dissertation also presents a methodology by which flood quantile estimates in Haiti can be derived using relationships developed for data rich regions of the U.S. As currently practiced, regional flood frequency techniques can only be applied within the predefined area used for model development. However, results presented herein demonstrate that the regional flood distribution can successfully be extrapolated to areas of similar physical composition located beyond the extent of that used for model development provided differences in precipitation are accounted for and the site in question can be appropriately classified within a delineated region.

Computational study of internal and external condensing flows and experimental synthesis to investigate their attainability and stability in ground-based and space-based environments

This doctoral thesis presents the computational work and synthesis with experiments for internal (tube and channel geometries) as well as external (flow of a pure vapor over a horizontal plate) condensing flows. The computational work obtains accurate numerical simulations of the full two dimensional governing equations for steady and unsteady condensing flows in gravity/0g environments. This doctoral work investigates flow features, flow regimes, attainability issues, stability issues, and responses to boundary fluctuations for condensing flows in different flow situations. This research finds new features of unsteady solutions of condensing flows; reveals interesting differences in gravity and shear driven situations; and discovers novel boundary condition sensitivities of shear driven internal condensing flows. Synthesis of computational and experimental results presented here for gravity driven in-tube flows lays framework for the future two-phase component analysis in any thermal system. It is shown for both gravity and shear driven internal condensing flows that steady governing equations have unique solutions for given inlet pressure, given inlet vapor mass flow rate, and fixed cooling method for condensing surface. But unsteady equations of shear driven internal condensing flows can yield different “quasi-steady” solutions based on different specifications of exit pressure (equivalently exit mass flow rate) concurrent to the inlet pressure specification. This thesis presents a novel categorization of internal condensing flows based on their sensitivity to concurrently applied boundary (inlet and exit) conditions. The computational investigations of an external shear driven flow of vapor condensing over a horizontal plate show limits of applicability of the analytical solution. Simulations for this external condensing flow discuss its stability issues and throw light on flow regime transitions because of ever-present bottom wall vibrations. It is identified that laminar to turbulent transition for these flows can get affected by ever present bottom wall vibrations. Detailed investigations of dynamic stability analysis of this shear driven external condensing flow result in the introduction of a new variable, which characterizes the ratio of strength of the underlying stabilizing attractor to that of destabilizing vibrations. Besides development of CFD tools and computational algorithms, direct application of research done for this thesis is in effective prediction and design of two-phase components in thermal systems used in different applications. Some of the important internal condensing flow results about sensitivities to boundary fluctuations are also expected to be applicable to flow boiling phenomenon. Novel flow sensitivities discovered through this research, if employed effectively after system level analysis, will result in the development of better control strategies in ground and space based two-phase thermal systems.

Development of Optically Based pH Sensing Hydrogel and Controlled Nitric Oxide Release Polymer

Nitric oxide has the potential to greatly improve intravascular measurements by locally inhibiting thrombus formation and dilating blood vessels. pH, the partial pressure of oxygen, and the partial pressure of carbon dioxide are three arterial blood parameters that are of interest to clinicians in the intensive care unit that can benefit from an intravascular sensor. This work explores fabrication of absorbance and fluorescence based pH sensing chemistry, the sensing chemistries' compatibility with nitric oxide, and a controllable nitric oxide releasing polymer. The pH sensing chemistries utilized various substrates, dyes, and methods of immobilization. Absorbance sensing chemistries used sol-gels, fumed silica particles, mesoporous silicon oxide, bromocresol purple, phenol red, bromocresol green, physical entrapment, molecular interactions, and covalent linking. Covalently linking the dyes to fumed silica particles and mesoporous silicon oxide eliminated leaching in the absorbance sensing chemistries. The structures of the absorbance dyes investigated were similar and bromocresol green in a sol-gel was tested for compatibility with nitric oxide. Nitric oxide did not interfere with the use of bromocresol green in a pH sensor. Investigated fluorescence sensing chemistries utilized silica optical fibers, poly(allylamine) hydrogel, SNARF-1, molecular interactions, and covalent linking. SNARF-1 covalently linked to a modified poly(allylamine) hydrogel was tested in the presence of nitric oxide and showed no interference from the nitric oxide. Nitric oxide release was controlled through the modulation of a light source that cleaved the bond between the nitric oxide and a sulfur atom in the donor. The nitric oxide donor in this work is S-nitroso-N-acetyl-D-penicillamine which was covalently linked to a silicone rubber made from polydimethylsiloxane. It is shown that the surface flux of nitric oxide released from the polymer films can be increased and decreased by increasing and decreasing the output power of the LED light source. In summary, an optical pH sensing chemistry was developed that eliminated the chronic problem of leaching of the indicator dye and showed no reactivity to nitric oxide released, thereby facilitating the development of a functional, reliable intravascular sensor.

DESIGN AND DEVELOPMENT OF BODIPY-BASED FLUORESCENT PROBES FOR SENSING AND IMAGING OF CYANIDE, Zn (II) IONS, LYSOSOMAL pH AND CANCER CELLS

BODIPY (4,4-Difluoro-3a,4a-diaza-s-indacene) dyes have gained lots of attention in application of fluorescence sensing and imaging in recent years because they possess many distinctive and desirable properties such as high extinction coefficient, narrow absorption and emission bands, high quantum yield and low photobleaching effect. However, most of BODIPY-based fluorescent probes have very poor solubilities in aqueous solution, emit less than 650 nm fluorescence that can cause cell and tissue photodamages compared with bio-desirable near infrared (650-900 nm) light. These undesirable properties extremely limit the applications of BODIPY-based fluorescent probes in sensing and imaging applications. In order to overcome these drawbacks, we have developed a very effective strategy to prepare a series of neutral highly water- soluble BODIPY dyes by enhancing the water solubilities of BODIPY dyes via incorporation of tri(ethylene glycol)methyl ether (TEG) and branched oligo(ethylene glycol)methyl ether (BEG) residues onto BODIPY dyes at 1,7-, 2,6-, 3,5-, 4- and meso- positions. We also have effectively tuned absorptions and emissions of BOIDPY dyes to red, deep red and near infrared regions via significant extension of π-conjugation of BODIPY dyes by condensation reactions of aromatic aldehydes with 2,6-diformyl BODIPY dyes at 1,3,5,7-positions. Based on the foundation that we built for enhancing water solubility and tuning wavelength, we have designed and developed a series of water-soluble, BODIPY-based fluorescent probes for sensitive and selective sensing and imaging of cyanide, Zn (II) ions, lysosomal pH and cancer cells. We have developed three BODIPY-based fluorescent probes for sensing of cyanide ions by incorporating indolium moieties onto the 6-position of TEG- or BEG-modified BOIDPY dyes. Two of them are highly water-soluble. These fluorescent probes showed selective and fast ratiometric fluorescent responses to cyanide ions with a dramatic fluorescence color change from red to green accompanying a significant increase in fluorescent intensity. The detection limit was measured as 0.5 mM of cyanide ions. We also have prepared three highly water-soluble fluorescent probes for sensing of Zn (II) ions by introducing dipicoylamine (DPA, Zn ion chelator) onto 2- and/or 6-positions of BEG-modified BODIPY dyes. These probes showed selective and sensitive responses to Zn (II) ion in the range from 0.5 mM to 24 mM in aqueous solution at pH 7.0. Particularly, one of the probes displayed ratiometric responses to Zn (II) ions with fluorescence quenching at 661 nm and fluorescence enhancement at 521 nm. This probe has been successfully applied to the detection of intracellular Zn (II) ions inside the living cells. Then, we have further developed three acidotropic, near infrared emissive BODIPY- based fluorescent probes for detection of lysosomal pH by incorporating piperazine moiety at 3,5-positions of TEG- or BEG-modified BODIPY dyes as parts of conjugation. The probes have low auto-fluorescence at physiological neutral condition while their fluorescence intensities will significant increase at 715 nm when pH shift to acidic condition. These three probes have been successfully applied to the in vitro imaging of lysosomes inside two types of living cells. At the end, we have synthesized one water- soluble, near infrared emissive cancer cell targetable BODIPY-based fluorescent polymer bearing cancer homing peptide (cRGD) residues for cancer cell imaging applications. This polymer exhibited excellent water-solubility, near infrared emission (712 nm), good biocompatibility. It also showed low nonspecific interactions to normal endothelial cells and can effectively detect breast tumor cells.

Development and characterization of fluorescent pH sensors based on porous silica and hydrogel support matrices

The hydrogen ion activity (pH) is a very important parameter in environment monitoring, biomedical research and other applications. Optical pH sensors have several advantages over traditional potentiometric pH measurement, such as high sensitivity, no need of constant calibration, easy for miniaturization and possibility for remote sensing. Several pH indicators has been successfully immobilized in three different solid porous materials to use as pH sensing probes. The fluorescent pH indicator fluorescein-5-isothiocyanate (FITC) was covalently bound onto the internal surface of porous silica (pore size ~10 nm) and retained its pH sensitivity. The excited state pK* a of FITC in porous silica (5.58) was slightly smaller than in solution (5.68) due to the free silanol groups (Si-OH) on the silica surface. The pH sensitive range for this probe is pH 4.5 - 7.0 with an error less than 0.1 pH units. The probe response was reproducible and stable for at least four month, stored in DI water, but exhibit a long equilibrium of up to 100 minutes. Sol-gel based pH sensors were developed with immobilization of two fluorescent pH indicators fluorescein-5-(and-6)-sulfonic acid, trisodium salt (FS) and 8-hydroxypyrene- 1,3,6-trisulfonic acid (HPTS) through physical entrapment. Prior to immobilization, the indicators were ion-paired with a common surfactant hexadecyltrimethylammonium bromide (CTAB) in order to prevent leaching. The sol-gel films were synthesized through the hydrolysis of two different precursors, ethyltriethoxysilane (ETEOS) and 3- glycidoxypropyltrimethoxysilane (GPTMS) and deposited on a quartz slide through spin coating. The pK a of the indicators immobilized in sol-gel films was much smaller than in solutions due to silanol groups on the inner surface of the sol-gel films and ammonium groups from the surrounding surfactants. Unlike in solution, the apparent pK a of the indicators in sol-gel films increased with increasing ionic strength. The equilibrium time for these sensors was within 5 minutes (with film thickness of ~470 nm). Polyethylene glycol (PEG) hydrogel was of interest for optical pH sensor development because it is highly proton permeable, transparent and easy to synthesize. pH indicators can be immobilized in hydrogel through physical entrapment and copolymerization. FS and HPTS ion-pairs were physically entrapped in hydrogel matrix synthesized via free radical initiation. For covalent immobilization, three indicators, 6,8-dihydroxypyrene-1,3- disulfonic acid (DHPDS), 2,7-dihydroxynaphthalene-3,6-disulfonic acid (DHNDS) and cresol red were first reacted with methacrylic anhydride (MA) to form methacryloylanalogs for copolymerization. These hydrogels were synthesized in aqueous solution with a redox initiation system. The thickness of the hydrogel film is controlled as ~ 0.5 cm and the porosity can be adjusted with the percentage of polyethylene glycol in the precursor solutions. The pK a of the indicators immobilized in the hydrogel both physically and covalently were higher than in solution due to the medium effect. The sensors are stable and reproducible with a short equilibrium time (less than 4 minutes). In addition, the color change of cresol red immobilized hydrogel is vivid from yellow (acidic condition) to purple (basic condition). Due to covalently binding, cresol red was not leaching out from the hydrogel, making it a good candidate of reusable "pH paper".