THE NATURAL AND ORTHOGONAL INTERACTION (NOIA) MODELS FOR QUANTITATIVE TRAITS (QTs) AND COMPLEX DISEASES

My dissertation focuses on developing methods for gene-gene/environment interactions and imprinting effect detections for human complex diseases and quantitative traits. It includes three sections: (1) generalizing the Natural and Orthogonal interaction (NOIA) model for the coding technique originally developed for gene-gene (GxG) interaction and also to reduced models; (2) developing a novel statistical approach that allows for modeling gene-environment (GxE) interactions influencing disease risk, and (3) developing a statistical approach for modeling genetic variants displaying parent-of-origin effects (POEs), such as imprinting. In the past decade, genetic researchers have identified a large number of causal variants for human genetic diseases and traits by single-locus analysis, and interaction has now become a hot topic in the effort to search for the complex network between multiple genes or environmental exposures contributing to the outcome. Epistasis, also known as gene-gene interaction is the departure from additive genetic effects from several genes to a trait, which means that the same alleles of one gene could display different genetic effects under different genetic backgrounds. In this study, we propose to implement the NOIA model for association studies along with interaction for human complex traits and diseases. We compare the performance of the new statistical models we developed and the usual functional model by both simulation study and real data analysis. Both simulation and real data analysis revealed higher power of the NOIA GxG interaction model for detecting both main genetic effects and interaction effects. Through application on a melanoma dataset, we confirmed the previously identified significant regions for melanoma risk at 15q13.1, 16q24.3 and 9p21.3. We also identified potential interactions with these significant regions that contribute to melanoma risk. Based on the NOIA model, we developed a novel statistical approach that allows us to model effects from a genetic factor and binary environmental exposure that are jointly influencing disease risk. Both simulation and real data analyses revealed higher power of the NOIA model for detecting both main genetic effects and interaction effects for both quantitative and binary traits. We also found that estimates of the parameters from logistic regression for binary traits are no longer statistically uncorrelated under the alternative model when there is an association. Applying our novel approach to a lung cancer dataset, we confirmed four SNPs in 5p15 and 15q25 region to be significantly associated with lung cancer risk in Caucasians population: rs2736100, rs402710, rs16969968 and rs8034191. We also validated that rs16969968 and rs8034191 in 15q25 region are significantly interacting with smoking in Caucasian population. Our approach identified the potential interactions of SNP rs2256543 in 6p21 with smoking on contributing to lung cancer risk. Genetic imprinting is the most well-known cause for parent-of-origin effect (POE) whereby a gene is differentially expressed depending on the parental origin of the same alleles. Genetic imprinting affects several human disorders, including diabetes, breast cancer, alcoholism, and obesity. This phenomenon has been shown to be important for normal embryonic development in mammals. Traditional association approaches ignore this important genetic phenomenon. In this study, we propose a NOIA framework for a single locus association study that estimates both main allelic effects and POEs. We develop statistical (Stat-POE) and functional (Func-POE) models, and demonstrate conditions for orthogonality of the Stat-POE model. We conducted simulations for both quantitative and qualitative traits to evaluate the performance of the statistical and functional models with different levels of POEs. Our results showed that the newly proposed Stat-POE model, which ensures orthogonality of variance components if Hardy-Weinberg Equilibrium (HWE) or equal minor and major allele frequencies is satisfied, had greater power for detecting the main allelic additive effect than a Func-POE model, which codes according to allelic substitutions, for both quantitative and qualitative traits. The power for detecting the POE was the same for the Stat-POE and Func-POE models under HWE for quantitative traits.

The Role of Type I Insulin-Like Growth Factor Receptor Signaling in Breast Cancer Brain Metastasis

Brain metastasis is a common cause of mortality in cancer patients. Approximately 20-30% of breast cancer patients acquire brain metastasis, yet potential therapeutic targets remain largely unknown. The type I insulin-like growth factor receptor (IGF- IR) is known to play a role in the progression of breast cancer and is currently being investigated in the clinical setting for various types of cancer. The present study demonstrates that the IGF-IR signaling axis is constitutively active in brain-seeking sublines of breast cancer cells, driving an increase in in vitro metastatic properties. We demonstrate that IGF-IR signaling is activated in an autocrine manner as a result of IGFBP3 overexpression in brain-seeking cells. Transient and stable knockdown of IGF-IR results in a downregulation of IGF-IR downstream signaling through phospho-AKT, as well as decreased in vitro migration and invasion of MDA- MB-231Br brain-seeking cells. Using an in vivo experimental brain metastasis model, we show that IGF-IR ablation attenuates the establishment of brain metastases and prolongs survival. Finally, we demonstrate that the malignancy of brain-seeking cells is attenuated by pharmacological inhibition with picropodophyllin, an IGF-IR-specific tyrosine kinase inhibitor. Together, our data suggest that the IGF-IR is an important mediator of brain metastasis and its ablation delays the onset of brain metastases in our model system.

Risk factors for surgical site infections among patients undergoing major colon surgery in the United States 2007-2009

Background: Surgical site infections (SSIs) after abdominal surgeries account for approximately 26% of all reported SSIs. The Center for Disease Control and Prevention (CDC) defines 3 types of SSIs: superficial incisional, deep incisional, and organ/space. Preventing SSIs has become a national focus. This dissertation assesses several associations with the individual types of SSI in patients that have undergone colon surgery. ^ Methods: Data for this dissertation was obtained from the American College of Surgeons' National Surgical Quality Improvement Program (NSQIP); major colon surgeries were identified in the database that occurred between the time period of 2007 and 2009. NSQIP data includes more than 50 preoperative and 30 intraoperative factors; 40 collected postoperative occurrences are based on a follow-up period of 30 days from surgery. Initially, four individual logistic regressions were modeled to compare the associations between risk factors and each of the SSI groups: superficial, deep, organ/space and a composite of any single SSI. A second analysis used polytomous regression to assess simultaneously the associations between risk factors and the different types of SSIs, as well as, formally test the different effect estimates of 13 common risk factors for SSIs. The final analysis explored the association between venous thromboembolism (VTEs) and the different types of SSIs and risk factors. ^ Results: A total of 59,365 colon surgeries were included in the study. Overall, 13% of colon cases developed a single type of SSI; 8% of these were superficial SSIs, 1.4% was deep SSIs, and 3.8% were organ/space SSIs. The first article identifies the unique set of risk factors associated with each of the 4 SSI models. Distinct risk factors for superficial SSIs included factors, such as alcohol, chronic obstructive pulmonary disease, dyspnea and diabetes. Organ/space SSIs were uniquely associated with disseminated cancer, preoperative dialysis, preoperative radiation treatment, bleeding disorder and prior surgery. Risk factors that were significant in all models had different effect estimates. The second article assesses 13 common SSI risk factors simultaneously across the 3 different types of SSIs using polytomous regression. Then each risk factor was formally tested for the effect heterogeneity exhibited. If the test was significant the final model would allow for the effect estimations for that risk factor to vary across each type of SSI; if the test was not significant, the effect estimate would remain constant across the types of SSIs using the aggregate SSI value. The third article explored the relationship of venous thromboembolism (VTE) and the individual types of SSIs and risk factors. The overall incidence of VTEs after the 59,365 colon cases was 2.4%. All 3 types of SSIs and several risk factors were independently associated with the development of VTEs. ^ Conclusions: Risk factors associated with each type of SSI were different in patients that have undergone colon surgery. Each model had a unique cluster of risk factors. Several risk factors, including increased BMI, duration of surgery, wound class, and laparoscopic approach, were significant across all 4 models but no statistical inferences can be made about their different effect estimates. These results suggest that aggregating SSIs may misattribute and hide true associations with risk factors. Using polytomous regression to assess multiple risk factors with the multiple types of SSI, this study was able to identify several risk factors that had significant effect heterogeneity across the 3 types of SSI challenging the use of aggregate SSI outcomes. The third article recognizes the strong association between VTEs and the 3 types of SSIs. Clinicians understand the difference between superficial, deep and organ/space SSIs. Our results indicate that they should be considered individually in future studies.^

Zircons, isotope ratios and element analyses from tonalite and oxide gabbro of the mid-ocean ridge from ODP Hole 176-735B

A limiting factor in the accuracy and precision of U/Pb zircon dates is accurate correction for initial disequilibrium in the 238U and 235U decay chains. The longest-lived-and therefore most abundant-intermediate daughter product in the 235U isotopic decay chain is 231Pa (T1/2 = 32.71 ka), and the partitioning behavior of Pa in zircon is not well constrained. Here we report high-precision thermal ionization mass spectrometry (TIMS) U-Pb zircon data from two samples from Ocean Drilling Program (ODP) Hole 735B, which show evidence for incorporation of excess 231Pa during zircon crystallization. The most precise analyses from the two samples have consistent Th-corrected 206Pb/238U dates with weighted means of 11.9325 ± 0.0039 Ma (n = 9) and 11.920 ± 0.011 Ma (n = 4), but distinctly older 207Pb/235U dates that vary from 12.330 ± 0.048 Ma to 12.140 ± 0.044 Ma and 12.03 ± 0.24 to 12.40 ± 0.27 Ma, respectively. If the excess 207Pb is due to variable initial excess 231Pa, calculated initial (231Pa)/(235U) activity ratios for the two samples range from 5.6 ± 1.0 to 9.6 ± 1.1 and 3.5 ± 5.2 to 11.4 ± 5.8. The data from the more precisely dated sample yields estimated DPazircon/DUzircon from 2.2-3.8 and 5.6-9.6, assuming (231Pa)/(235U) of the melt equal to the global average of recently erupted mid-ocean ridge basaltic glasses or secular equilibrium, respectively. High precision ID-TIMS analyses from nine additional samples from Hole 735B and nearby Hole 1105A suggest similar partitioning. The lower range of DPazircon/DUzircon is consistent with ion microprobe measurements of 231Pa in zircons from Holocene and Pleistocene rhyolitic eruptions (Schmitt (2007; doi:10.2138/am.2007.2449) and Schmitt (2011; doi:10.1146/annurev-earth-040610-133330)). The data suggest that 231Pa is preferentially incorporated during zircon crystallization over a range of magmatic compositions, and excess initial 231Pa may be more common in zircons than acknowledged. The degree of initial disequilibrium in the 235U decay chain suggested by the data from this study, and other recent high precision datasets, leads to resolvable discordance in high precision dates of Cenozoic to Mesozoic zircons. Minor discordance in zircons of this age may therefore reflect initial excess 231Pa and does not require either inheritance or Pb loss.

Age models, accumulation rates, sand and carbonates at DSDP Leg 74 holes

Accumulation rates for the five sites drilled during Leg 74 of the Glomar Challenger are presented on a common timescale based on calibration of datum levels to paleomagnetic records in Leg 74 sediments for the Paleogene, and a new compilation by Berggren et al. (1985), for the Neogene, and using the seafloor-spreading magnetic anomaly timescale of Kent (1985). We present data on accumulation of total sediment, of foraminifers, of the noncarbonate portion, and of fish teeth that give a history of productivity, winnowing, carbonate dissolution, and nonbiogenic input to what was then a part of the South Atlantic at about 30 deg S.

Elevation Models of Antarctica and Greenland derived from CryoSat-2 in the period 2011 to 2013

This study focuses on the present-day surface elevation of the Greenland and Antarctic ice sheets. Based on 3 years of CryoSat-2 data acquisition we derived new elevation models (DEMs) as well as elevation change maps and volume change estimates for both ice sheets. Here we present the new DEMs and their corresponding error maps. The accuracy of the derived DEMs for Greenland and Antarctica is similar to those of previous DEMs obtained by satellite-based laser and radar altimeters. Comparisons with ICESat data show that 80% of the CryoSat-2 DEMs have an uncertainty of less than 3 m ± 15 m. The surface elevation change rates between January 2011 and January 2014 are presented for both ice sheets. We compared our results to elevation change rates obtained from ICESat data covering the time period from 2003 to 2009. The comparison reveals that in West Antarctica the volume loss has increased by a factor of 3. It also shows an anomalous thickening in Dronning Maud Land, East Antarctica which represents a known large-scale accumulation event. This anomaly partly compensates for the observed increased volume loss of the Antarctic Peninsula and West Antarctica. For Greenland we find a volume loss increased by a factor of 2.5 compared to the ICESat period with large negative elevation changes concentrated at the west and southeast coasts. The combined volume change of Greenland and Antarctica for the observation period is estimated to be -503 ± 107 km**3/yr. Greenland contributes nearly 75% to the total volume change with -375 ± 24 km**3/yr.

Age models of upper Pliocene samples of the North Atlantic

Abundance variations of six Pliocene species of discoasters have been analyzed over the time interval from 1.89 to 2.95 Ma at five contrasting sites in the North Atlantic: Deep Sea Drilling Project Sites 552 (56°N) and 607 (41°N) and Ocean Drilling Program 658 (20°N), 659 (18°N), and 662 (1°S). A sampling interval equivalent to approximately 3 k.y. was used. Total Discoaster abundance showed a reduction with increasing latitude and from the effects of upwelling. This phenomenon is most obvious in Discoaster brouweri, the only species that survived over the entire time interval studied. Prior to 2.38 Ma, Discoaster pentaradiatus and Discoaster surculus are important components of the Discoaster assemblage: Discoaster pentaradiatus increases slightly with latitude up to 41°N, and its abundance relative to D. brouweri increases up to 56°N; D. surculus increases in abundance with latitude and with upwelling conditions relative to both D. brouweri and D. pentaradiatus and is dominant to the latter species at upwelling Site 658 and at the highest latitude sites. Discoaster asymmetricus and Discoaster tamalis appear to increase in abundance with latitude relative to D. brouweri. Many of the abundance changes observed appear to be connected with the initiation of glaciation in the North Atlantic at 2.4 Ma. The long-term trend of decreasing Discoaster abundance probably reflects the fall of sea-surface temperatures. This trend of cooling is overprinted by short-term variations that are probably associated with orbital forcing. Evidence for the astronomical elements of eccentricity and obliquity periodicities were found at all sites; however, only at Sites 607, 659, and 662 were precessional periodicities detected. Furthermore, only at Site 659 was precession found to be dominant to obliquity. Abundance peaks of individual species were found to cross-correlate between sites. The distinct abundance fluctuations observed especially in the tropics suggest that temperature is not the only factor responsible for this variation. This study reveals for the first time the importance of productivity pressure on the suppression of Discoaster abundance.

Geochemistry of abyssal peridotites from ODP Hole 209-1274A

Aqueous dihydrogen (H2,aq) is produced in copious amounts when seawater interacts with peridotite and H2O oxidizes ferrous iron in olivine to ferric iron in secondary magnetite and serpentine. Poorly understood in this process is the partitioning of iron and its oxidation state in serpentine, although both impose an important control on dihydrogen production. We present results of detailed petrographic, mineral chemical, magnetic and Mößbauer analyses of partially to fully serpentinized peridotites from the Ocean Drilling Program (ODP) Leg 209, Mid-Atlantic Ridge (MAR) 15°N area. These results are used to constrain the fate of iron during serpentinization and are compared with phase equilibria considerations and peridotite-seawater reaction path models. In samples from Hole 1274A, mesh-rims reveal a distinct in-to-out zoning from brucite at the interface with primary olivine, followed by a zone of serpentine + brucite ± magnetite and finally serpentine + magnetite in the outermost mesh-rim. The compositions of coexisting serpentine (Mg# 95) and brucite (Mg# 80) vary little throughout the core. About 30-50% of the iron in serpentine/brucite mesh-rims is trivalent, irrespective of subbasement depth and protolith (harzburgite versus dunite). Model calculations suggest that both partitioning and oxidation state of iron are very sensitive to temperature and water-to-rock ratio during serpentinization. At temperatures above 330 °C the dissolution of olivine and coeval formation of serpentine, magnetite and dihydrogen depends on the availability of an external silica source. At these temperatures the extent of olivine serpentinization is insufficient to produce much hydrogen, hence conditions are not reducing enough to form awaruite. At T < 330 °C, hydrogen generation is facilitated by the formation of brucite, as dissolution of olivine to form serpentine, magnetite and brucite requires no addition of silica. The model calculations suggest that the iron distribution observed in serpentine and brucite is consistent with formation temperatures ranging from <150 to 250 °C and bulk water-to-rock ratios between 0.1 and 5. These conditions coincide with peak hydrogen fugacities during serpentinization and are conducive to awaruite formation during main stage serpentinization. The development of the common brucite rims around olivine is either due to an arrested reaction olivine -> brucite -> serpentine + brucite, or reflects metastable olivine-brucite equilibria developing in the strong gradient in silica activity between orthopyroxene (talc-serpentine) and olivine (serpentine-brucite).

Age models and sea-surface paleotemperature reconstruction of sediment cores from the eastern equatorial Atlantic

Based on the faunal record of planktonic foraminifers in three long gravity sediment cores from the eastern equatorial Atlantic, the sea-surface temperature history ove the last 750,000 years was studied at a resolution of 3,000 to 10,000 years. Detailed oxygen-isotope and paleomagnetic stratigraphy helped to identify the following major faunal events: Globorotaloides hexagonus and Globorotalia tumida flexuosa became extinct in the eastern tropical Atlantic at the isotope stage 4/5 boundary, now dated at 68,000 years B.P. The persistent occurrence of the pink variety of Globigerinoides ruber started during the late stage 12 at 410,000 years B.P. CARTUNE-age. This datum may provide an easily detectible faunal stratigraphic marker for the mid-Brunhes Chron. The updated scheme of the Ericson zones helped the recognition of a hiatus at the northwestern slope of the Sierra Leone Basin covering oxygen-isotope stages 10 to 12. Classifying the planktonic foraminifer counts into six faunal assemblages, according to the factor analysis derived model of Pflaumann (1985), the tropical and the tropical-upwelling communities account for 57 % at Site 16415, and 86 % at Site 13519, respectively of the variance of the faunal record. A largely continuous paleotemperature record for both winter and summer seasons was obtained from the top of the Sierra Leone Rise with the winter temperatures ranging between 20 and 25 °C, and the summer ones between 24 and 30 °C. The record of cores from greater water depths is frequently interrupted by samples with no-analogue faunal communities and/or poor preservation. Based on the seasonality signal, during cold periods the termal equator shifted to a geographically mnore asymmetrical northern position. Dissolution altering the faunal communities becomes stronger with greater water depth, the estimated mean minimum loss of specimens increases from 70 % to 80 % between 2,860 and 3,850 water depth although some species will be more susceptible than others. Enhanced dissolution occured during stage 4 but also during cold phases in the warm stage 7 and 9. Correlations between the Foraminiferal Dissolution Index and the estimated sea-surface temperatures are significant. Foraminiferal flux rates, negatively correlated to the flux rates of organic carbon and of diatoms, may be a result of enhanced dissolution during cold stages, destroying still more of the faunal signal than indicated by the calculated minimum loss. The fluctuations of the oxygen-isotope curves and the hibernal sea-surfave temperatures are fairly coherent. During warm oxygen-isotope stages the temperature maxima lag often by 5 to 15 ka behind the respective sotope minima. During cold stages, sea-surface temperature changes are partly out of phase and contain additional fluctuations.

Geochemistry and kinetic models od Cretaceous samples

The source rock potential of Cretaceous organic rich whole rock samples from deep sea drilling project (DSDP) wells offshore southwestern Africa was investigated using bulk and quantitative pyrolysis techniques. The sample material was taken from organic rich intervals of Aptian, Albian and Turonian aged core samples from DSDP site 364 offshore Angola, DSDP well 530A north of the Walvis Ridge offshore Namibia, and DSDP well 361 offshore South Africa. The analytical program included TOC, Rock-Eval, pyrolysis GC, bulk kinetics and micro-scale sealed vessel pyrolysis (MSSV) experiments. The results were used to determine differences in the source rock petroleum type organofacies, petroleum composition, gas/oil ratio (GOR) and pressure-volume-temperature (PVT) behavior of hydrocarbons generated from these black shales for petroleum system modeling purposes. The investigated Aptian and Albian organic rich shales proved to contain excellent quality marine kerogens. The highest source rock potential was identified in sapropelic shales in DSDP well 364, containing very homogeneous Type II and organic sulfur rich Type IIS kerogen. They generate P-N-A low wax oils and low GOR sulfur rich oils, whereas Type III kerogen rich silty sandstones of DSDP well 361 show a potential for gas/condensate generation. Bulk kinetic experiments on these samples indicate that the organic sulfur contents influence kerogen transformation rates, Type IIS kerogen being the least stable. South of the Walvis Ridge, the Turonian contains predominantly a Type III kerogen. North of the Walvis Ridge, the Turonian black shales contain Type II kerogen and have the potential to generate P-N-A low and high wax oils, the latter with a high GOR at high maturity. Our results provide the first compositional kinetic description of Cretaceous organic rich black shales, and demonstrate the excellent source rock potential, especially of the Aptian-aged source rock, that has been recognized in a number of the South Atlantic offshore basins.

Predicting bed form roughness: the influence of lee side angle

Flow transverse bedforms (ripples and dunes) are ubiquitous in rivers and coastal seas. Local hydrodynamics and transport conditions depend on the size and geometry of these bedforms, as they constitute roughness elements at the bed. Bedform influence on flow energy must be considered for the understanding of flow dynamics, and in the development and application of numerical models. Common estimations or predictors of form roughness (friction factors) are based mostly on data of steep bedforms (with angle-of-repose lee slopes), and described by highly simplified bedform dimensions (heights and lengths). However, natural bedforms often are not steep, and differ in form and hydraulic effect relative to idealised bedforms. Based on systematic numerical model experiments, this study shows how the hydraulic effect of bedforms depends on the flow structure behind bedforms, which is determined by the bedform lee side angle, aspect ratio and relative height. Simulations reveal that flow separation behind bedform crests and, thus, a hydraulic effect is induced at lee side angles steeper than 11 to 18° depending on relative height, and that a fully developed flow separation zone exists only over bedforms with a lee side angle steeper than 24°. Furthermore, the hydraulic effect of bedforms with varying lee side angle is evaluated and a reduction function to common friction factors is proposed. A function is also developed for the Nikuradse roughness (k s), and a new equation is proposed which directly relates k s to bedform relative height, aspect ratio and lee side angle.

Gridded climate data from 5 Global Climate Models (GCM) of the Last Glacial Maximum (LGM) downscaled to 30 arc seconds for Europe, with links to NetCDF files

Studies on the impact of historical, current and future global change require very high-resolution climate data (less or equal 1km) as a basis for modelled responses, meaning that data from digital climate models generally require substantial rescaling. Another shortcoming of available datasets on past climate is that the effects of sea level rise and fall are not considered. Without such information, the study of glacial refugia or early Holocene plant and animal migration are incomplete if not impossible. Sea level at the last glacial maximum (LGM) was approximately 125m lower, creating substantial additional terrestrial area for which no current baseline data exist. Here, we introduce the development of a novel, gridded climate dataset for LGM that is both very high resolution (1km) and extends to the LGM sea and land mask. We developed two methods to extend current terrestrial precipitation and temperature data to areas between the current and LGM coastlines. The absolute interpolation error is less than 1°C and 0.5 °C for 98.9% and 87.8% of all pixels for the first two 1 arc degree distance zones. We use the change factor method with these newly assembled baseline data to downscale five global circulation models of LGM climate to a resolution of 1km for Europe. As additional variables we calculate 19 'bioclimatic' variables, which are often used in climate change impact studies on biological diversity. The new LGM climate maps are well suited for analysing refugia and migration during Holocene warming following the LGM.

University-industry technology transfer models: an empirical analysis

The aim of this paper is to contribute to the understanding of the underlying factors in the process of transferring technology from university to industry. Findings point to strategic importance of critical factors as the definition of common objectives, cooperation, motivation, and the elimination of technical and legal barriers. These challenges must have implications in the incorporation of cooperative aspects of research projects in the design of public innovation policies.

Uniformly Reweighted Belief Propagation: A Factor Graph Approach

Tree-reweighted belief propagation is a message passing method that has certain advantages compared to traditional belief propagation (BP). However, it fails to outperform BP in a consistent manner, does not lend itself well to distributed implementation, and has not been applied to distributions with higher-order interactions. We propose a method called uniformly-reweighted belief propagation that mitigates these drawbacks. After having shown in previous works that this method can substantially outperform BP in distributed inference with pairwise interaction models, in this paper we extend it to higher-order interactions and apply it to LDPC decoding, leading performance gains over BP.