Bayesian Inference for Retrospective Population Genetics Models Using Markov Chain Monte Carlo Methods

Genetics, the science of heredity and variation in living organisms, has a central role in medicine, in breeding crops and livestock, and in studying fundamental topics of biological sciences such as evolution and cell functioning. Currently the field of genetics is under a rapid development because of the recent advances in technologies by which molecular data can be obtained from living organisms. In order that most information from such data can be extracted, the analyses need to be carried out using statistical models that are tailored to take account of the particular genetic processes. In this thesis we formulate and analyze Bayesian models for genetic marker data of contemporary individuals. The major focus is on the modeling of the unobserved recent ancestry of the sampled individuals (say, for tens of generations or so), which is carried out by using explicit probabilistic reconstructions of the pedigree structures accompanied by the gene flows at the marker loci. For such a recent history, the recombination process is the major genetic force that shapes the genomes of the individuals, and it is included in the model by assuming that the recombination fractions between the adjacent markers are known. The posterior distribution of the unobserved history of the individuals is studied conditionally on the observed marker data by using a Markov chain Monte Carlo algorithm (MCMC). The example analyses consider estimation of the population structure, relatedness structure (both at the level of whole genomes as well as at each marker separately), and haplotype configurations. For situations where the pedigree structure is partially known, an algorithm to create an initial state for the MCMC algorithm is given. Furthermore, the thesis includes an extension of the model for the recent genetic history to situations where also a quantitative phenotype has been measured from the contemporary individuals. In that case the goal is to identify positions on the genome that affect the observed phenotypic values. This task is carried out within the Bayesian framework, where the number and the relative effects of the quantitative trait loci are treated as random variables whose posterior distribution is studied conditionally on the observed genetic and phenotypic data. In addition, the thesis contains an extension of a widely-used haplotyping method, the PHASE algorithm, to settings where genetic material from several individuals has been pooled together, and the allele frequencies of each pool are determined in a single genotyping.

Identification of loci for migraine with aura

Common migraine, i.e. migraine with (MA) or without aura (MO), is a chronic neurological disorder affecting about 10% of the Caucasian population. In MA, migraine headache is preceded by visual, sensoric and/or dysphasic reversible aura symptoms. Twin and family studies have suggested a multifactorial mode of inheritance for common migraine, and a stronger genetic component for MA than for MO. Since there is no biological or genetic marker to identify common migraine, aura symptoms provide a distinctive character to identify those suspected of suffering from migraine. The aim of this study was to identify MA susceptibility loci in well-phenotyped migraine samples with familial predisposition using different gene mapping methods. Genes coding for endothelin1 and its receptors EDNRA and ENDRB are potential candidate genes for cortical spreading depression (CSD), which is considered to be the underlying mechanism of migraine aura. The role of these genes in MA was studied in 850 Finnish migraine cases and 890 control individuals. Rare homozygous EDNRA SNPs showed nominal association with MA and with the age of onset trait (20 years). This result was also detected in the pooled analysis on 648 German MA cases and 651 control individuals when the test was adjusted for gender and sample origin. Evaluation of SNP genotyping reactions with two different DNA polymerase enzymes ensured that the genotype quality was high, and thus the discovered associations are considered reliable. The role of the 19p13 region was studied in a linkage analysis of 72 Finnish MA families. This region contains two migraine-associated genes: CACNA1A, which is associated with a predisposition to a rare Mendelian form of MA, familial hemiplegic migraine (FHM), and the insulin receptor gene (INSR) that is associated with common migraine. No evidence of linkage between the 19p13 and MA was detected. A novel visual aura locus was mapped to chromosome 9q21-q22 with significant evidence of linkage using a genome-wide linkage approach in 36 Finnish MA families. Five additional, potential loci were also detected. The 9q21-q22 region has previously been linked to occipitotemporal lobe epilepsy and MA, both of which involve prominent visual symptoms. Our result further supports a shared background for these episodic disorders.

Major histocompatibility complex and coronary artery disease: Special emphasis on Chlamydia pneumoniae, and periodontitis

Most of the genes in the MHC region are involveed in adaptive and innate immunity, with essential function in inflammatory reactions and in protection against infections. These genes might serve as a candidate region for infection and inflammation associated diseases. CAD is an inflammatory disease. The present set of studies was performed to assess whether the MHC region harbors genetic markers for CAD, and whether these genetic markers explain the CAD risk factors: e.g., C. pneumoniae, periodontitis, and periodontal pathogens. Study I was performed using two separate patient materials and age- and sex-matched healthy controls, categorizing them into two independent studies: the HTx and ACS studies. Both studies consistently showed the HLA-A3– B35– DR1 (35 ancestral haplotype) haplotype as a susceptible MHC genetic marker for CAD. HLA-DR1 alone was associated not only with CAD, but also with CAD risk factor diseases, e.g., diabetes mellitus, and hyperlipidemia. The ACS study further showed the HLA-B*07 and -DRB1*15 -related haplotype as a protective MHC haplotype for CAD. Study II showed that patients with CAD showed signs of chronic C. pneumoniae infection when compared to age- and sex-matched healthy controls. HLA-B*35 or -related haplotypes associated with the C. pneumoniae infection markers. Among these haplotype carriers, males and smokers associated with elevated C. pneumoniae infection markers. Study III showed that CAD patients with periodontitis had elevated serum markers of P. gingivalis and occurrence of the pathogen in saliva. LTA+496C strongly associated with periodontitis, while HLA-DRB1*01 with periodontitis and with the elevated serum antibodies of P. gingivalis. Study IV showed that the increased level of C3/C4 ratio was a new risk factor and was associated with recurrent cardiovascular end-points. The increased C3 and decreased C4 concentrations in serum explained the increased level of the C3/C4 ratio. Both the higher than cut-off value (4.53) and the highest quartile of the C3/C4 ratio were also associated with worst survival, increased end-points, and C4 null alleles. The presence of C4 null alleles associated with decreased serum C4 concentration, and increased C3/C4 ratio. In conclusion, the present studies show that the CAD susceptibility haplotype (HLA-A3− B35− DR1 -related haplotypes, Study I) partially explains the development of CAD in patients possessing several recognized and novel risk factors: diabetes mellitus, increased LDL, smoking, C4B*Q0, C. pneumnoiae, periodontitis, P. gingivalis, and complement C3/C4 ratio (Study II, III, and IV).

The DISC1 Pathway in the Genetic Etiology of Schizophrenia

Schizophrenia is a severe psychotic disorder affecting 0.5-1 % of the population. The disorder is characterized by hallucinations; delusions; disorganized behavior and speech; avolition; anhedonia; flattened affect and cognitive deficits. The etiology of the disorder is complex with evidence for multiple genes contributing to the onset of the disorder along with environmental factors. DISC1 is one of the most promising candidate genes for schizophrenia. It codes for a protein which takes part in numerous molecular interactions along several pathways. This network, termed as the DISC1 pathway, is evidently important for the development and maturation of the central nervous system from the embryo until young adulthood. Disruption at these pathways is thought to predispose schizophrenia. In the present study, we have studied the DISC1 pathway in the etiology of schizophrenia in the Finnish population. We have utilized large Finnish samples; the schizophrenia family sample where DISC1 was originally shown to associate with schizophrenia and the Northern Finland birth cohort 1966 (NFBC66). Several DISC1 binding partners displayed evidence for association in the family sample along with DISC1. Through a genome-wide linkage study, we found a significant linkage signal to a locus where a DISC1 binding partner NDE1 is located at the carriers of a certain DISC1 risk variant. In a follow-up study, genetic markers in NDE1 displayed significant evidence for association with schizophrenia. Further exploration of association between 11 genes of the DISC1 pathway and schizophrenia led to recognition of novel variants in NDEL1, PDE4B and PDE4D that significantly either increased or decreased the risk for schizophrenia. Further, we found evidence that DISC1 itself has a significant role in the human mental functioning even in the healthy population. Variants in DISC1 had a significant effect on anhedonia which is a trait present at everybody but is in its severe form one of the main symptoms of schizophrenia and correlates with the risk of developing the disorder. Further, utilizing genome-wide marker data, we recognized three genes; MIR620; CCDC141 and LCT; that are closely related to the DISC1 pathway but which effects on anhedonia were observable only at the individuals who carried these specific DISC1 variants. Our findings significantly add up to the previous evidence for the involvement of DISC1 and the DISC1 pathway in the etiology of schizophrenia and psychosis. Our results support the concept of a number of DISC1 pathway related genes contributing in the etiology of schizophrenia along with DISC1 and provide new candidates for the studies of schizophrenia. Our findings also significantly increase the importance of DISC1 itself as having a role in psychological functioning in the general population.

Molecular Genetic Background of Tumours in Hereditary Leiomyomatosis and Renal Cell Cancer Syndrome

Hereditary Leiomyomatosis and Renal Cell Cancer (HLRCC) is a hereditary tumour predisposition syndrome. Its phenotype includes benign cutaneous and uterine leiomyomas (CLM, ULM) with high penetrance and rarer renal cell cancer (RCC), most commonly of papillary type 2 subtype. Over 130 HLRCC families have been identified world-wide but the RCC phenotype seems to concentrate in families from Finland and North America for unknown reasons. HLRCC is caused by heterozygous germline mutations in the fumarate hydratase (FH) gene. FH encodes the enzyme fumarase from mitochondrial citric acid cycle. Fumarase enzyme activity or type or site of the FH mutation are unassociated with disease phenotype. The strongest evidence for tumourigenesis mechanism in HLRCC supports a hypoxia inducible factor driven process called pseudohypoxia resulting from accumulation of the fumarase substrate fumarate. In this study, to assess the importance of gene- or exon-level deletions or amplifications of FH in patients with HLRCC-associated phenotypes, multiplex ligation-dependent probe amplification (MLPA) method was used. One novel FH mutation, deletion of exon 1, was found in a Swedish male patient with an evident HLRCC phenotype with CLM, RCC, and a family history of ULM and RCC. Six other patients with CLM and 12 patients with only RCC or uterine leiomyosarcoma (ULMS) remained FH mutation-negative. These results suggest that copy number aberrations of FH or its exons are an infrequent cause of HLRCC and that only co-occurrence of benign tumour types justifies FH-mutation screening in RCC or ULMS patients. Determination of the genomic profile of 11 HLRCC-associated RCCs from Finnish patients was performed by array comparative genomic hybridization. The most common copy number aberrations were gains of 2, 7, and 17 and losses of 13q12.3-q21.1, 14, 18, and X. When compared to aberrations of sporadic papillary RCCs, HLRCC-associated RCCs harboured a distinct DNA copy number profile and lacked many of the changes characterizing the sporadic RCCs. The findings suggest a divergent molecular pathway for tumourigenesis of papillary RCCs in HLRCC. In order to find a genetic modifier of RCC risk in HLRCC, genome-wide linkage and identical by descent (IBD) analysis studies were performed in Finnish HLRCC families with microsatellite marker mapping and SNP-array platforms. The linkage analysis identified only one locus of interest, the FH gene locus in 1q43, but no mutations were found in the genes of the region. IBD analysis yielded no convincing haplotypes shared by RCC patients. Although these results do not exclude the existence of a genetic modifier for RCC risk in HLRCC, they emphasize the role of FH mutations in the malignant tumourigenesis of HLRCC. To study the benign tumours in HLRCC, genome-wide DNA copy number and gene expression profiles of sporadic and HLRCC ULMs were defined with modern SNP- and gene-expression array platforms. The gene expression array suggests novel genes involved in FH-deficient ULM tumourigenesis and novel genes with putative roles in propagation of sporadic ULM. Both the gene expression and copy number profiles of HLRCC ULMs differed from those of sporadic ULMs indicating distinct molecular basis of the FH-deficient HLRCC tumours.

Quantitative trait loci detection and benefits from marker-assisted selection in dairy cattle

Knowing the chromosomal areas or actual genes affecting the traits under selection would add more information to be used in the selection decisions which would potentially lead to higher genetic response. The first objective of this study was to map quantitative trait loci (QTL) affecting economically important traits in the Finnish Ayrshire population. The second objective was to investigate the effects of using QTL information in marker-assisted selection (MAS) on the genetic response and the linkage disequilibrium between the different parts of the genome. Whole genome scans were carried out on a grand-daughter design with 12 half-sib families and a total of 493 sons. Twelve different traits were studied: milk yield, protein yield, protein content, fat yield, fat content, somatic cell score (SCS), mastitis treatments, other veterinary treatments, days open, fertility treatments, non-return rate, and calf mortality. The average spacing of the typed markers was 20 cM with 2 to 14 markers per chromosome. Associations between markers and traits were analyzed with multiple marker regression. Significance was determined by permutation and genome-wise P-values obtained by Bonferroni correction. The benefits from MAS were investigated by simulation: a conventional progeny testing scheme was compared to a scheme where QTL information was used within families to select among full-sibs in the male path. Two QTL on different chromosomes were modelled. The effects of different starting frequencies of the favourable alleles and different size of the QTL effects were evaluated. A large number of QTL, 48 in total, were detected at 5% or higher chromosome-wise significance. QTL for milk production were found on 8 chromosomes, for SCS on 6, for mastitis treatments on 1, for other veterinary treatments on 5, for days open on 7, for fertility treatments on 7, for calf mortality on 6, and for non-return rate on 2 chromosomes. In the simulation study the total genetic response was faster with MAS than with conventional selection and the advantage of MAS persisted over the studied generations. The rate of response and the difference between the selection schemes reflected clearly the changes in allele frequencies of the favourable QTL. The disequilibrium between the polygenes and QTL was always negative and it was larger with larger QTL size. The disequilibrium between the two QTL was larger with QTL of large effect and it was somewhat larger with MAS for scenarios with starting frequencies below 0.5 for QTL of moderate size and below 0.3 for large QTL. In conclusion, several QTL affecting economically important traits of dairy cattle were detected. Further studies are needed to verify these QTL, check their presence in the present breeding population, look for pleiotropy and fine map the most interesting QTL regions. The results of the simulation studies show that using MAS together with embryo transfer to pre-select young bulls within families is a useful approach to increase the genetic merit of the AI-bulls compared to conventional selection.

The life cycle and genetic structure of the red alga Furcellaria lumbricalis on a salinity gradient

The life cycle and genetic diversity of the red alga Furcellaria lumbricalis (Hudson) Lamouroux were investigated in 15 populations in northern Europe. The occurrence of different life cycle phases and seasonality of reproduction were studied in four brackish populations in the northern Baltic Sea. Furthermore, a new method, based on genome screening with ISSR markers combined with a restriction-ligation method, was developed to discover microsatellite markers for population genetic analyses. The mitochondrial DNA cox2-3 spacer sequence and four microsatellite markers were used to examine the genetic diversity and differentiation of red algal populations in northern Europe. In addition, clonality and small-scale genetic structure of one Irish and four Baltic Sea populations were studied with microsatellite markers. It was discovered that at the low salinities of the northern Baltic Sea, only tetrasporophytes and males were present in the populations of F. lumbricalis and that winter was the main season for tetrasporangial production. Furthermore, the population occurring at the lowest salinity (3.6 practical salinity units, psu) did not produce spores. The size of the tetraspores was smaller in the Baltic Sea populations than that in the Irish population, and there were more deformed spores in the Baltic Sea populations than in the Irish populations. Studies with microsatellite markers indicated that clonality is a common phenomenon in the Baltic Sea populations of F. lumbricalis, although the proportion of clonal individuals varied among populations. Some genetic divergence occurred within locations both in Ireland and in the northern Baltic Sea. Even though no carpogonia were detected in the field samples during the study, the microsatellite data indicated that sexual reproduction occurs at least occasionally in the northern Baltic Sea. The genetic diversity of F. lumbricalis was highest in Brittany, France. Since no variation was discovered in the mtDNA cox2-3 spacer sequence, which is generally regarded as an informative phylogeographic marker in red algae, it can be assumed that the studied populations probably share the same origin.

Genetic aspects of outcome in kidney transplantation: cytokine and thrombosis associated candidate genes and gene expression biomarkers

Kidney transplantation (Tx) is the treatment of choice for end stage renal disease. Immunosuppressive medications are given to prevent an immunological rejection of the transplant. However, immunosuppressive drugs increase e.g. the risk of infection, cancer or nephrotoxicity. A major genetic contributors to immunological acceptance of the graft are human leukocyte antigen (HLA) genes. Also other non-HLA gene polymorphisms may predict the future risk of complications before Tx, possibly enabling individualised immunotherapy. Graft function after Tx is monitored using non-specific clinical symptoms and laboratory markers. The definitive diagnosis of graft rejection however relies on a biopsy of the graft. In the acute rejection (AR) diagnostics there is a need for an alternative to biopsy that would be an easily repeatable and simple method for regular use. Frequent surveillance of acute or subclinical rejection (SCR) may improve long-term function. In this thesis, associations between cytokine and thrombosis associated candidate genes and the outcome of kidney Tx were studied. Cytotoxic and co-stimulatory T lymphocyte molecule gene expression biomarkers for the diagnosis of the AR and the SCR were also investigated. We found that polymorphisms in the cytokine genes tumor necrosis factor and interleukin 10 (IL10) of the recipients were associated with AR. In addition, certain IL10 gene polymorphisms of the donors were associated with the incidence of cytomegalovirus infection and occurrence of later infection in a subpopulation of recipients. Further, polymorphisms in genes related to the risk of thrombosis and those of certain cytokines were not associated with the occurrence of thrombosis, infarction, AR or graft survival. In the study of biomarkers for AR, whole blood samples were prospectively collected from adult kidney Tx patients. With real-time quantitative PCR (RT-QPCR) gene expression quantities of CD154 and ICOS differentiated the patients with AR from those without, but not from the patients with other causes of graft dysfunction. Biomarkers for SCR were studied in paediatric kidney Tx patients. We used RT-QPCR to quantify the gene expression of immunological candidate genes in a low-density array format. In addition, we used RT-QPCR to validate the results of the microarray analysis. No gene marker differentiated patients with SCR from those without SCR. This research demonstrates the lack of robust markers among polymorphisms or biomarkers in investigated genes that could be included in routine analysis in a clinical laboratory. In genetic studies, kidney Tx can be regarded as a complex trait, i.e. several environmental and genetic factors may determine its outcome. A number of currently unknown genetic factors probably influence the results of Tx.

Genetic analysis of chromosomal regions 2q33, 7q32 and 19q13 in multiple sclerosis susceptibility

Multiple sclerosis (MS) is an immune-mediated demyelinating disorder of the central nervous system (CNS) affecting 0.1-0.2% of Northern European descent population. MS is considered to be a multifactorial disease, both environment and genetics play a role in its pathogenesis. Despite several decades of intense research, the etiological and pathogenic mechanisms underlying MS remain still largely unknown and no curative treatment exists. The genetic architecture underlying MS is complex with multiple genes involved. The strongest and the best characterized predisposing genetic factors for MS are located, as in other immune-mediated diseases, in the major histocompatibility complex (MHC) on chromosome 6. In humans MHC is called human leukocyte antigen (HLA). Alleles of the HLA locus have been found to associate strongly with MS and remained for many years the only consistently replicable genetic associations. However, recently other genes located outside the MHC region have been proposed as strong candidates for susceptibility to MS in several studies. In this thesis a new genetic locus located on chromosome 7q32, interferon regulatory factor 5 (IRF5), was identified in the susceptibility to MS. In particular, we found that common variation of the gene was associated with the disease in three different populations, Spanish, Swedish and Finnish. We also suggested a possible functional role for one of the risk alleles with impact on the expression of the IRF5 locus. Previous studies have pointed out a possible role played by chromosome 2q33 in the susceptibility to MS and other autoimmune disorders. The work described here also investigated the involvement of this chromosomal region in MS predisposition. After the detection of genetic association with 2q33 (article-1), we extended our analysis through fine-scale single nucleotide polymorphism (SNP) mapping to define further the contribution of this genomic area to disease pathogenesis (article-4). We found a trend (p=0.04) for association to MS with an intronic SNP located in the inducible T-cell co-stimulator (ICOS) gene, an important player in the co-stimulatory pathway of the immune system. Expression analysis of ICOS revealed a novel, previously uncharacterized, alternatively spliced isoform, lacking the extracellular domain that is needed for ligand binding. The stability of the newly-identified transcript variant and its subcellular localization were analyzed. These studies indicated that the novel isoform is stable and shows different subcellular localization as compared to full-length ICOS. The novel isoform might have a regulatory function, but further studies are required to elucidate its function. Chromosome 19q13 has been previously suggested as one of the genomic areas involved in MS predisposition. In several populations, suggestive linkage signals between MS predisposition and 19q13 have been obtained. Here, we analysed the role of allelic variation in 19q13 by family based association analysis in 782 MS families collected from Finland. In this dataset, we were not able to detect any statistically significant associations, although several previously suggested markers were included to the analysis. Replication of the previous findings on the basis of linkage disequilibrium between marker allele and disease/risk allele appears notoriously difficult because of limitations such as allelic heterogeneity. Re-sequencing based approaches may be required for elucidating the role of chromosome 19q13 with MS. This thesis has resulted in the identification of a new MS susceptibility locus (IRF5) previously associated with other inflammatory or autoimmune disorders, such as SLE. IRF5 is one of the mediators of interferons biological function. In addition to providing new insight in the possible pathogenetic pathway of the disease, this finding suggests that there might be common mechanisms between different immune-mediated disorders. Furthermore the work presented here has uncovered a novel isoform of ICOS, which may play a role in regulatory mechanisms of ICOS, an important mediator of lymphocyte activation. Further work is required to uncover its functions and possible involvement of the ICOS locus in MS susceptibility.

Genetic Variation and Clinical Manifestations in Inflammatory Bowel Disease

Crohn s disease (CD) and ulcerative colitis (UC), collectively known as inflammatory bowel disease (IBD), are characterised by chronic inflammation of the gastrointestinal tract. IBD prevalence in Finland is approximately 3-4 per 1000 inhabitants with a peak incidence in adolescence. The symptoms of IBD include diarrhoea, abdominal pain, fever, and weight loss. The precise aetiology of IBD is unknown but interplay of environmental risk factors and immunologic changes trigger the disease in a genetically susceptible individual. Twin and family studies have provided strong evidence for genetic factors in IBD susceptibility, and genetic factors may be more prominent in CD than UC. The first CD susceptibility gene was identified in 2001. Three common mutations R702W, G908R, and 1007fs of the CARD15/NOD2 gene are shown to associate independently with CD but the magnitude of association varies between different populations. The present study aimed at identifying mutations and genetic variations in IBD susceptibility and candidate genes. In addition, correlation to phenotype was also assessed. One of the main objectives of this study was to evaluate the role of CARD15 in a Finnish CD cohort. 271 CD patients were studied for the three common mutations and the results showed a lower mutation frequency than in other Caucasian populations. Only 16% of the patients carried one of the three mutations. Ileal location as well as stricturing and penetrating behaviour of the disease were associated with occurrence of the mutations. The whole protein coding region of CARD15 was screened for possible Finnish founder mutations. In addition to several sequence variants, five novel mutations (R38M, W355X, P727L, W907R, and R1019X) were identified in five patients. Functional consequences of these novel variants were studied in vitro, and these studies demonstrated a profound impairment of MDP response. Investigation of CARD15 mutation frequency in healthy people across three continents showed a large geographic fluctuation. No simple correlation between mutation frequency and disease incidence was seen in populations studied. The occurrence of double mutant carriers in healthy controls suggested that the penetrance of risk alleles is low. Other main objectives aimed at identifying other genetic variations that are involved in the susceptibility to IBD. We investigated the most plausible IBD candidate genes including TRAF6, SLC22A4, SLC22A5, DLG5, TLR4, TNFRSF1A, ABCB1/MDR1, IL23R, and ATG16L1. The marker for a chromosome 5 risk haplotype and the rare HLA-DRB1*0103 allele were also studied. The study cohort consisted of 699 IBD patients (240 CD and 459 UC), of which 23% had a first-degree relative with IBD. Of the several candidate genes studied, IL23R was associated with CD susceptibility, and TNFRSF1A as well as the HLA-DRB1*0103 allele with UC susceptibility. IL23R variants also showed association with the stricturing phenotype and longer disease duration in CD patients. In addition, TNFRSF1A variants were more common among familial UC and ileocolonic CD. In conclusion, the common CARD15 mutations were shown to account for 16% of CD cases in Finland. Novel CARD15 variants identified in the present study are most likely disease-causing mutations, as judged by the results of in vitro studies. The present study also confirms the IL23R association with CD susceptibility and, in addition, TNFRSF1A and HLA-DRB1*0103 allele association with UC of specific clinical phenotypes.

Work Stress and Early Atherosclerosis: Do Genetic Background and Pre-Employment Risk Factors Explain Conflicting Findings?

According to the models conceptualizing work stress, increased risk of health problems arise when high job demands co-occur with low job control (the demand-control model) or the efforts invested by the employee are disproportionately high compared to the rewards received (effort-reward imbalance model). This study examined the association between work stress and early atherosclerosis with particular attention to the role of pre-employment risk factors and genetic background in this association. The subjects were young healthy adults aged 24-39 who were participating in the 21-year follow-up of the ongoing prospective "Cardiovascular Risk in Young Finns" study in 2001-2002. Work stress was evaluated with questionnaires on demand-control model and on effort-reward model. Atherosclerosis was assessed with ultrasound of carotid artery intima-media thickness (IMT). In addition, risk for enhanced atherosclerotic process was assessed by measuring with heart rate variability and heart rate. Pre-employment risk factors, measured at age 12 to 18, included such as body mass index, blood lipids, family history of coronary heart disease, and parental socioeconomic position. Variants of the neuregulin-1 were determined using genomic DNA. The results showed that higher work stress was associated with higher IMT in men. This association was not attenuated by traditional risk factors of atherosclerosis and coronary heart disease or by pre-employment risk factors measured in adolescence. Neuregulin-1 gene moderated the association between work stress and IMT in men. A significant association between work stress and IMT was found only for the T/T genotype of the neuregulin-1 gene but not for other genotypes. Among women an association was found between higher work stress and lower heart rate variability, suggesting higher risk for developing atherosclerosis. These associations could not be explained by demographic characteristics or coronary risk factors. The present findings provide evidence for an association between work stress and atherosclerosis in relatively young population. This association seems to be modified by genetic influences but it does not appear to be confounded by pre-employment adolescent risk factors.

Genetic and Epigenetic Characterization of Colorectal and Endometrial Cancer

Colorectal cancer is one of the three most common cancers today, for both men and women. Approximately 90% of the cases are sporadic while the remaining 10% is hereditary. Among this 10% is hereditary nonpolyposis colorectal cancer (HNPCC), an autosomal dominant disease, accounting for up to 13% of these cases. HNPCC is associated with germline mutations in four mismatch repair (MMR) genes, MLH1, MSH2, MSH6, and PMS2, and is characterized by a familial accumulation of endometrial, gastric, urological, and ovarian tumors, in addition to colorectal cancer. An important etiological characteristic of HNPCC is the presence of microsatellite instability (MSI), caused by mutations of the MMR genes. Approximately 15% of sporadic cases share the MSI+ trait. Colon cancer is believed to be a consequence of an accumulation of mutations in tumor suppressor genes and oncogenes, eventually resulting in tumor development. This phenomena is accelerated in HNPCC due the presence of an inherited mutation in the MMR genes, accounting for one of the two hits proposed to be needed by Knudson (1971) in order for the manifestation of the MSI phenotype. MMR alterations alone, however, do not occur in the majority of sporadic colon cancers, prompting searches for other mechanisms. One such mechanism found to play a role in colon cancer development was DNA methylation, which is known to play a role in MLH1 inactivation. Our objective was clarification of mechanisms associated with tumor development in both HNPCC and sporadic colorectal cancer in relation to tumorigenic mechanisms. Of particular interest were underlying mechanisms of MSI in sporadic colorectal cancers, with attention to DNA methylation changes and their correlation to MSI. Of additional interest were the genetic and epigenetic events leading to the HNPCC tumor spectrum, chiefly colon and endometrial cancers, in regards to what extent the somatic changes in target tissue explained this phenomenon. We made a number of important findings pertaining to these questions. First, MSI tumor development differs epigenetically from stable tumor development, possibly underlying developmental pathway differences. Additionally, while epigenetic modification, principally DNA methylation, is a major mechanism in sporadic MSI colorectal cancer MLH1 inactivation it does not play a significant role in HNPCC tumors with germline MLH1 mutations. This is possibly an explanation for tumorigenic pathways and clinicopathological characteristic differences between sporadic and hereditary MSI colorectal cancers. Finally, despite indistinguishable genetic predisposition for endometrial and colorectal cancers, instability profiles highlighting organ-specific differences, may be important HNPCC tumor spectrum determinants.