Sociodemographic, behavioral and genetic determinants of allostatic load in a Swiss population-based study.

Allostatic load (AL) is a marker of physiological dysregulation which reflects exposure to chronic stress. High AL has been related to poorer health outcomes including mortality. We examine here the association of socioeconomic and lifestyle factors with AL. Additionally, we investigate the extent to which AL is genetically determined. We included 803 participants (52% women, mean age 48±16years) from a population and family-based Swiss study. We computed an AL index aggregating 14 markers from cardiovascular, metabolic, lipidic, oxidative, hypothalamus-pituitary-adrenal and inflammatory homeostatic axes. Education and occupational position were used as indicators of socioeconomic status. Marital status, stress, alcohol intake, smoking, dietary patterns and physical activity were considered as lifestyle factors. Heritability of AL was estimated by maximum likelihood. Women with a low occupational position had higher AL (low vs. high OR=3.99, 95%CI [1.22;13.05]), while the opposite was observed for men (middle vs. high OR=0.48, 95%CI [0.23;0.99]). Education tended to be inversely associated with AL in both sexes(low vs. high OR=3.54, 95%CI [1.69;7.4]/OR=1.59, 95%CI [0.88;2.90] in women/men). Heavy drinking men as well as women abstaining from alcohol had higher AL than moderate drinkers. Physical activity was protective against AL while high salt intake was related to increased AL risk. The heritability of AL was estimated to be 29.5% ±7.9%. Our results suggest that generalized physiological dysregulation, as measured by AL, is determined by both environmental and genetic factors. The genetic contribution to AL remains modest when compared to the environmental component, which explains approximately 70% of the phenotypic variance.

Bayesian models for the analysis of genetic structure when populations are correlated

Motivation: Population allele frequencies are correlated when populations have a shared history or when they exchange genes. Unfortunately, most models for allele frequency and inference about population structure ignore this correlation. Recent analytical results show that among populations, correlations can be very high, which could affect estimates of population genetic structure. In this study, we propose a mixture beta model to characterize the allele frequency distribution among populations. This formulation incorporates the correlation among populations as well as extending the model to data with different clusters of populations. Results: Using simulated data, we show that in general, the mixture model provides a good approximation of the among-population allele frequency distribution and a good estimate of correlation among populations. Results from fitting the mixture model to a dataset of genotypes at 377 autosomal microsatellite loci from human populations indicate high correlation among populations, which may not be appropriate to neglect. Traditional measures of population structure tend to over-estimate the amount of genetic differentiation when correlation is neglected. Inference is performed in a Bayesian framework.

Nonlinear tests for genetic studies of complex disease

Linkage and association studies are major analytical tools to search for susceptibility genes for complex diseases. With the availability of large collection of single nucleotide polymorphisms (SNPs) and the rapid progresses for high throughput genotyping technologies, together with the ambitious goals of the International HapMap Project, genetic markers covering the whole genome will be available for genome-wide linkage and association studies. In order not to inflate the type I error rate in performing genome-wide linkage and association studies, multiple adjustment for the significant level for each independent linkage and/or association test is required, and this has led to the suggestion of genome-wide significant cut-off as low as 5 × 10 −7. Almost no linkage and/or association study can meet such a stringent threshold by the standard statistical methods. Developing new statistics with high power is urgently needed to tackle this problem. This dissertation proposes and explores a class of novel test statistics that can be used in both population-based and family-based genetic data by employing a completely new strategy, which uses nonlinear transformation of the sample means to construct test statistics for linkage and association studies. Extensive simulation studies are used to illustrate the properties of the nonlinear test statistics. Power calculations are performed using both analytical and empirical methods. Finally, real data sets are analyzed with the nonlinear test statistics. Results show that the nonlinear test statistics have correct type I error rates, and most of the studied nonlinear test statistics have higher power than the standard chi-square test. This dissertation introduces a new idea to design novel test statistics with high power and might open new ways to mapping susceptibility genes for complex diseases. ^

Immune and genetic risk factors in glioma

Despite extensive research, the etiology of adult glioma remains largely unknown. We sought to further explore the role of immune and genetic factors in glioma etiology using data from the Harris County Brain Tumor Study and the first U.S. genome-wide association study of glioma. First, using a case-control study design, we examined the association between adult glioma risk and surrogates of the timing and frequency of common early childhood infections, birth order and sibship size, respectively. We found that each one-unit increase in birth order was associated with a 12% decreased risk of glioma development in adulthood (OR=0.88, 95% CI=0.81-0.96); however, sibship size was not associated with adult glioma risk (OR=0.96, 95% CI=0.91-1.02). Second, we used a multi-strategic approach to explore the relationships between glioma risk, history of asthma/allergies, and 23 functional SNPs in 11 inflammation genes. We found three inflammation gene SNPs to be significantly associated with glioma risk (COX2/PTGS2 rs20417 [OR=1.41]; IL10 rs1800896 [OR=1.57]; and IL13 rs20541 [OR=0.39]). Joint effects analysis of the risk-conferring alleles of these three SNPs revealed a trend of increasing risk with increasing number of adverse alleles among those without asthma/allergies (p<0.0001). Finally, we conducted a case-only study to explore pairwise SNP-SNP interactions in immune-related pathways among a population of 1304 non-Hispanic white glioma cases. After correction for multiple comparisons, we found 279 significant SNP-SNP interactions among polymorphisms of immune-related genes, many of which have not been previously examined. Our results, cumulatively, suggest an important role for immune and genetic factors in glioma etiology and provide several new hypotheses for future studies.^

Genetic epidemiology of gallbladder disease in Mexican Americans and cholesterol 7$\alpha$-hydroxylase gene sequence variation

Among Mexican Americans, the second largest minority group in the United States, the prevalence of gallbladder disease is markedly elevated. Previous data from both genetic admixture and family studies indicate that there is a genetic component to the occurrence of gallbladder disease in Mexican Americans. However, prior to this thesis no formal genetic analysis of gallbladder disease had been carried out nor had any contributing genes been identified.^ The results of complex segregation analysis in a sample of 232 Mexican American pedigrees documented the existence of a major gene having two alleles with age- and gender-specific effects influencing the occurrence of gallbladder disease. The estimated frequency of the allele increasing susceptibility was 0.39. The lifetime probabilities that an individual will be affected by gallbladder disease were 1.0, 0.54, and 0.00 for females of genotypes "AA", "Aa", and "aa", respectively, and 0.68, 0.30, and 0.00 for males, respectively. This analysis provided the first conclusive evidence for the existence of a common single gene having a large effect on the occurrence of gallbladder disease.^ Human cholesterol 7$\alpha$-hydroxylase is the rate-limiting enzyme in bile acid synthesis. The results of an association study in both a random sample and a matched case/control sample showed that there is a significant association between cholesterol 7$\alpha$-hydroxylase gene variation and the occurrence of gallbladder disease in Mexican Americans males but not in females. These data have implicated a specific gene, 7$\alpha$-hydroxylase, in the etiology of gallbladder disease in this population.^ Finally, I asked whether the inferred major gene from complex segregation analysis is genetically linked to the cholesterol 7$\alpha$-hydroxylase gene. Three pedigrees predicted to be informative for linkage analysis by virtue of supporting the major gene hypothesis and having parents with informative genotypes and multiple offspring were selected for this linkage analysis. In each of these pedigrees, the recombination fractions maximized at 0 with a positive, albeit low, LOD score. The results of this linkage analysis provide preliminary and suggestive evidence that the cholesterol 7$\alpha$-hydroxylase gene and the inferred gallbladder disease susceptibility gene are genetically linked. ^

Genetic variants associated with metastatic prostate cancer

Introduction. Distant metastasis remains the leading cause of death among prostate cancer patients. Several genetic susceptibility loci associated with Prostate cancer have been identified by the Genome Wide Association Studies (GWAS). To date, few studies have explored the ability of these SNPs to identify metastatic prostate cancer. Based on the identification of genetic variants as predictors of aggressive disease, a case comparison study of prostate cancer patients was designed to explore the association of 96 GWAS single nucleotide polymorphisms (SNPs) with metastatic disease. ^ Method. 1242 histologically confirmed prostate cancer patients, with and without metastatic disease, were enrolled into the study. Data were collected from personal interviews, hospital database and abstraction of medical records. Ninety six SNPs identified from GWAS studies based on their associations with prostate cancer risk were genotyped in the study population. Univariate and multivariate logistic regression analyses were used to explore the relationships of the variants with metastatic prostate cancer in Whites and African American men. ^ Results. Four SNPs showed independent associations with metastatic prostate cancer (rs721048 in EHBP1 (2p15), rs3025039 in VEGF (6p12), rs11228565 in Intergenic(11q13.2) and rs2735839 in KLK3(19q13.33)) in the White population. For SNP rs2735839 in KLK3, genotype GA was 1.71 times as likely to be associated with metastatic prostate cancer diagnosis as genotype AA after adjusting for other significant SNPs and covariates (95% CI, 1.12-2.60; p=0.012). In men of African descent, three SNPs: rs1512268 in NKX3-1(8p21.2), rs12155172 in intergenic (7p15.3) & rs10486567 in JAZF1 (7p15.2) were positively associated with metastatic disease in the multivariate analysis. The strongest SNP was rs1512268 heterozygous genotype AG in NKX3-1(8p21.2) which was associated with 3.97-fold increased risk of metastatic prostate cancer diagnosis (95% CI, 1.69-9.34; p =0.002). ^ Conclusion. Genetic variants associated with metastatic prostate cancer were different in Whites and African American men. Given the high mortality rate recorded in men diagnosed with metastatic prostate tumor, further studies are needed to validate associations and establish their clinical application.^

Literature review of estimation of prevalence rate of HIV/AIDS in India using population survey and sentinel surveillance based data

HIV/AIDS is a treatable although incurable disease that presents immense challenges to those infected including physical, social and psychological effects. As of 2009, an estimated 2.4 million people were living with HIV or AIDS in India, 0.3% of the country's population. In India, it is difficult to not only treat but also to track because it is associated with socio-economic factors such as illiteracy, social biases, poor sanitation, malnutrition and social class. Nevertheless, it is important to know the prevalence of HIV/AIDS for several reasons. At the individual level, the quality of life of people living with HIV/AIDS is markedly lower than their counterparts without the disease and is associated with challenges. At the community level, it is important to identify high risk groups, monitor prevention efforts, and allocate appropriate resources to target programs for the reduction of transmission of HIV. ^

Exploration of genetic susceptibility to pancreatic cancer: A GWAS approach

Pancreatic cancer is the 4th most common cause for cancer death in the United States, accompanied by less than 5% five-year survival rate based on current treatments, particularly because it is usually detected at a late stage. Identifying a high-risk population to launch an effective preventive strategy and intervention to control this highly lethal disease is desperately needed. The genetic etiology of pancreatic cancer has not been well profiled. We hypothesized that unidentified genetic variants by previous genome-wide association study (GWAS) for pancreatic cancer, due to stringent statistical threshold or missing interaction analysis, may be unveiled using alternative approaches. To achieve this aim, we explored genetic susceptibility to pancreatic cancer in terms of marginal associations of pathway and genes, as well as their interactions with risk factors. We conducted pathway- and gene-based analysis using GWAS data from 3141 pancreatic cancer patients and 3367 controls with European ancestry. Using the gene set ridge regression in association studies (GRASS) method, we analyzed 197 pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Using the logistic kernel machine (LKM) test, we analyzed 17906 genes defined by University of California Santa Cruz (UCSC) database. Using the likelihood ratio test (LRT) in a logistic regression model, we analyzed 177 pathways and 17906 genes for interactions with risk factors in 2028 pancreatic cancer patients and 2109 controls with European ancestry. After adjusting for multiple comparisons, six pathways were marginally associated with risk of pancreatic cancer ( P < 0.00025): Fc epsilon RI signaling, maturity onset diabetes of the young, neuroactive ligand-receptor interaction, long-term depression (Ps < 0.0002), and the olfactory transduction and vascular smooth muscle contraction pathways (P = 0.0002; Nine genes were marginally associated with pancreatic cancer risk (P < 2.62 × 10−5), including five reported genes (ABO, HNF1A, CLPTM1L, SHH and MYC), as well as four novel genes (OR13C4, OR 13C3, KCNA6 and HNF4 G); three pathways significantly interacted with risk factors on modifying the risk of pancreatic cancer (P < 2.82 × 10−4): chemokine signaling pathway with obesity ( P < 1.43 × 10−4), calcium signaling pathway (P < 2.27 × 10−4) and MAPK signaling pathway with diabetes (P < 2.77 × 10−4). However, none of the 17906 genes tested for interactions survived the multiple comparisons corrections. In summary, our current GWAS study unveiled unidentified genetic susceptibility to pancreatic cancer using alternative methods. These novel findings provide new perspectives on genetic susceptibility to and molecular mechanisms of pancreatic cancer, once confirmed, will shed promising light on the prevention and treatment of this disease. ^

Decomposition of $\rm R\times C$ contingency table chi -square: Application to binned DNA fragment size data and population structure analysis

The purpose of this research is to develop a new statistical method to determine the minimum set of rows (R) in a R x C contingency table of discrete data that explains the dependence of observations. The statistical power of the method will be empirically determined by computer simulation to judge its efficiency over the presently existing methods. The method will be applied to data on DNA fragment length variation at six VNTR loci in over 72 populations from five major racial groups of human (total sample size is over 15,000 individuals; each sample having at least 50 individuals). DNA fragment lengths grouped in bins will form the basis of studying inter-population DNA variation within the racial groups are significant, will provide a rigorous re-binning procedure for forensic computation of DNA profile frequencies that takes into account intra-racial DNA variation among populations. ^

Data from PhD thesis

The ocean quahog, Arctica islandica is the longest-lived non-colonial animal known to science. A maximum individual age of this bivalve of 405 years has been found in a population off the north western coast of Iceland. Conspicuously shorter maximum lifespan potentials (MLSPs) were recorded from other populations of A. islandica in European waters (e.g. Kiel Bay: 30 years, German Bight: 150 years) which experience wider temperature and salinity fluctuations than the clams from Iceland. The aim of my thesis was to identify possible life-prolonging physiological strategies in A. islandica and to examine the modulating effects of extrinsic factors (e.g. seawater temperature, food availability) and intrinsic factors (e.g. species-specific behavior) on these strategies. Burrowing behavior and metabolic rate depression (MRD), tissue-specific antioxidant and anaerobic capacities as well as cell-turnover (= apoptosis and proliferation) rates were investigated in A. islandica from Iceland and the German Bight. An inter-species comparison of the quahog with the epibenthic scallop Aequipecten opercularis (MLSP = 8-10 years) was carried out in order to determine whether bivalves with short lifespans and different lifestyles also feature a different pattern in cellular maintenance and repair. The combined effects of a low-metabolic lifestyle, low oxidative damage accumulation, and constant investment into cellular protection and tissue maintenance, appear to slow-down the process of physiological aging in A. islandica and to afford the extraordinarily long MLSP in this species. Standard metabolic rates were lower in A. islandica when compared to the shorter-lived A. opercularis. Furthermore, A. islandica regulate mantle cavity water PO2 to mean values < 5 kPa, a PO2 at which the formation of reactive oxygen species (ROS) in isolated gill tissues of the clams was found to be 10 times lower than at normoxic conditions (21 kPa). Burrowing and metabolic rate depression (MRD) in Icelandic specimens were more pronounced in winter, possibly supported by low seawater temperature and food availability, and seem to be key energy-saving and life-prolonging parameters in A. islandica. The signaling molecule nitric oxide (NO) may play an important role during the onset of MRD in the ocean quahog by directly inhibiting cytochome-c-oxidase at low internal oxygenation upon shell closure. In laboratory experiments, respiration of isolated A. islandica gills was completely inhibited by chemically produced NO at low experimental PO2 <= 10 kPa. During shell closure, mantle cavity water PO2 decreased to 0 kPa for longer than 24 h, a state in which ROS production is supposed to subside. Compared to other mollusk species, onset of anaerobic metabolism is late in A. islandica in the metabolically reduced state. Increased accumulation of the anaerobic metabolite succinate was initially detected in the adductor muscle of the clams after 3.5 days under anoxic incubation or in burrowed specimens. A ROS-burst was absent in isolated gill tissue of the clams following hypoxia (5 kPa)-reoxygenation (21 kPa). Accordingly, neither the activity of antioxidant enzymes superoxide dismutase (SOD) and catalase (CAT), nor the specific content of the ROS-scavenger glutathione (GSH) was enhanced in different tissues of the ocean quahog after 3.5 days of self-induced or forced hypoxia/anoxia to prepare for an oxidative burst. While reduced ROS formation compared to routine levels lowers oxidative stress during MRD and also during surfacing, the general preservation of high cellular defense and the efficient removal and replacement of damaged cells over lifetime seem to be of crucial importance in decelerating the senescent decline in tissues of A. islandica. Along with stable antioxidant protection over 200 years of age, proliferation rates and apoptosis intensities in most investigated tissues of the ocean quahog were low, but constant over 140 years of age. Accordingly, age-dependent accumulations of protein and lipid oxidation products are lower in A. islandica tissues when compared to the shorter-lived bivalve A. opercularis. The short-lived swimming scallop is a model bivalve species representing the opposite life and aging strategy to A. islandica. In this species permanently high energy throughput, reduced investment into antioxidant defense with age, and higher accumulation of oxidation products are met by higher cell turnover rates than in the ocean quahog. The only symptoms of physiological change over age ever found in A. islandica were decreasing cell turnover rates in the heart muscle over a lifetime of 140 years. This may either indicate higher damage levels and possibly ongoing loss of functioning in the heart of aging clams, or, the opposite, lower rates of cell damage and a reduced need for cell renewal in the heart tissue of A. islandica over lifetime. Basic physiological capacities of different A. islandica populations, measured at controlled laboratory conditions, could not explain considerable discrepancies in population specific MLSPs. For example, levels of tissue-specific antioxidant capacities and cell turnover rates were similarly high in individuals from the German Bight and from Iceland. Rather than genetic differences, the local impacts of environmental conditions on behavioral and physiological traits in the ocean quahog seem to be responsible for differences in population-specific MLSPs.

AFLP data based on presence (1) and absence (0) of AFLP loci (A4-F171) in Gonyostomum semen

Many recent studies have found genetically differentiated populations in microorganisms despite potentially high dispersal. We designed a study to specifically examine the importance of physical dispersal barriers, i.e. geographic distance and lack of hydrological connectivity, in restricting gene flow and enhancing divergence in limnic microorganisms. We focused on the nuisance microalga Gonyostomum semen, which has recently expanded in northern Europe and differentiated into genetically distinct populations. Gonyostomum semen was sampled from six lakes distributed in two adjacent watersheds, which thereby comprised, both connected and non-connected lakes. The individual isolates were genotyped by Amplified Fragment Length Polymorphism. Several lake populations were differentiated from each other, but connectivity within watersheds could not explain the observed population genetic pattern. However, isolation by distance was moderate and might limit the gene flow among distant populations. In addition, we found low, but significant linkage disequilibrium, which indicates regular sexual recombination in this species, despite its high degree of asexual reproduction. Therefore, we conclude that the genetic properties of microalgae with occasional sexual reproduction essentially mirror regularly recombining species. Furthermore, the data indicated bottlenecks supporting the hypothesized recent range expansion of this species.

PumaPlex 1.0: Data generated during the development of the SNP markers and SNP genotypes of pumas

Pumas are one of the most studied terrestrial mammals because of their widespread distribution, substantial ecological impacts, and conflicts with humans. Extensive efforts, often employing genetic methods, are undertaken to manage this species. However, the comparison of population genetic data is difficult because few of the microsatellite loci chosen are shared across research programs. Here, we describe the development of PumaPlex, a high-throughput assay to genotype 25 single nucleotide polymorphisms in pumas. We validated PumaPlex in more than 700 North American pumas (Puma concolor couguar), and demonstrated its ability to generate reproducible genotypes and accurately identify individuals. Furthermore, we compared PumaPlex with traditional genotyping of 12 microsatellite loci in fecal DNA samples and found that PumaPlex produced significantly more genotypes with fewer false alleles. PumaPlex promotes the cross-laboratory comparison of genotypes, is easily expandable in the future, and is a valuable tool for the genetic monitoring and management of North American puma populations.

Experimental data on photophysiology and growth rates of Pararotalia calcariformata

The eastern Mediterranean is a hotspot of biological invasions. Numerous species of Indo-pacific origin have colonized the Mediterranean in recent times, including tropical symbiont-bearing foraminifera. Among these is the species Pararotalia calcariformata. Unlike other invasive foraminifera, this species has been discovered only two decades ago and is restricted to the eastern Mediterranean coast. Combining ecological, genetic and physiological observations, we attempt to explain the recent invasion of this species in the Mediterranean Sea. Using morphological and genetic data, we confirm the species attribution to P. calcariformata McCulloch 1977 and identify its symbionts as a consortium of diatom species dominated by Minutocellus polymorphus. We document photosynthetic activity of its endosymbionts using Pulse Amplitude Modulated Fluorometry and test the effects of elevated temperatures on growth rates of asexual offspring. The culturing of asexual offspring for 120 days shows a 30-day period of rapid growth followed by a period of slower growth. A subsequent 48-day temperature sensitivity experiment indicates a similar developmental pathway and high growth rate at 28°C, whereas an almost complete inhibition of growth was observed at 20°C and 35°C. This indicates that the offspring of this species may have lower tolerance to cold temperatures than what would be expected for species native to the Mediterranean. We expand this hypothesis by applying a Species Distribution Model (SDM) based on modern occurrences in the Mediterranean using three environmental variables: irradiance, turbidity and yearly minimum temperature. The model reproduces the observed restricted distribution and indicates that the range of the species will drastically expand westwards under future global change scenarios. We conclude that P. calcariformata established a population in the Levant because of the recent warming in the region. In line with observations from other groups of organisms, our results indicate that continued warming of the eastern Mediterranean will facilitate the invasion of more tropical marine taxa into the Mediterranean, disturbing local biodiversity and ecosystem structure.

Genetic data from Staurosira (Fragilariophyceae) sedimentary DNA

Genetic investigations on eukaryotic plankton confirmed the existence of modern biogeographic patterns, but analyses of palaeoecological data exploring the temporal variability of these patterns have rarely been presented. Ancient sedimentary DNA proved suitable for investigations of past assemblage turnover in the course of environmental change, but genetic relatedness of the identified lineages has not yet been undertaken. Here, we investigate the relatedness of diatom lineages in Siberian lakes along environmental gradients (i.e. across treeline transects), over geographic distance and through time (i.e. the last 7000 years) using modern and ancient sedimentary DNA. Our results indicate that closely-related Staurosira lineages occur in similar environments and less-related lineages in dissimilar environments, in our case different vegetation and co-varying climatic and limnic variables across treeline transects. Thus our study reveals that environmental conditions rather than geographic distance is reflected by diatom-relatedness patterns in space and time. We tentatively speculate that the detected relatedness pattern in Staurosira across the treeline could be a result of adaptation to diverse environmental conditions across the arctic boreal treeline, however, a geographically-driven divergence and subsequent repopulation of ecologically different habitats might also be a potential explanation for the observed pattern.