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.

Molecular mediators linking stroke and carotid artery disease

Carotid artery disease is the most prevalent etiologic precursor of ischemic stroke, which is a major health hazard and the second most common cause of death in the world. If a patient presents with a symptomatic high-grade (>70%) stenosis in the internal carotid artery, the treatment of choice is carotid endarterectomy. However, the natural course of radiologically equivalent carotid lesions may be clinically quite diverse, and the reason for that is unknown. It would be of utmost importance to develop molecular markers that predict the symptomatic phenotype of an atherosclerotic carotid plaque (CP) and help to differentiate vulnerable lesions from stable ones. The aim of this study was to investigate the morphologic and molecular factors that associate with stroke-prone CPs. In addition to immunohistochemistry, DNA microarrays were utilized to identify molecular markers that would differentiate between symptomatic and asymptomatic CPs. Endothelial adhesion molecule expression (ICAM-1, VCAM-1, P-selectin, and E-selectin) did not differ between symptomatic and asymptomatic patients. Denudation of endothelial cells was associated with symptom-generating carotid lesions, but in studies on the mechanism of decay of endothelial cells, markers of apoptosis (TUNEL, activated caspase 3) were found to be decreased in the endothelium of symptomatic lesions. Furthermore, markers of endothelial apoptosis were directly associated with those of cell proliferation (Ki-67) in all plaques. FasL expression was significantly increased on the endothelium of symptomatic CPs. DNA microarray analysis revealed prominent induction of specific genes in symptomatic CPs, including those subserving iron and heme metabolism, namely HO-1, and hemoglobin scavenger receptor CD163. HO-1 and CD163 proteins were also increased in symptomatic CPs and associated with intraplaque iron deposits, which, however, did not correlate with symptom status itself. ADRP, the gene for adipophilin, was also overexpressed in symptomatic CPs. Adipophilin expression was markedly increased in ulcerated CPs and colocalized with extravasated red blood cells and cholesterol crystals. Taken together, the phenotypic characteristics and the numerous possible molecular mediators of the destabilization of carotid plaques provide potential platforms for future research. The denudation of the endothelial lining observed in symptomatic CPs may lead to direct thromboembolism and maintain harmful oxidative and inflammatory processes, predispose to plaque microhemorrhages, and contribute to lipid accumulation into the plaque, thereby making it vulnerable to rupture.

The influence of OLR1 and PCSK9 gene polymorphisms on ischemic stroke: Evidence from a meta-analysis

It has been reported that both OLR1 and PCSK9 genes are related to various vascular diseases such as atherosclerosis, cardiovascular disease, peripheral artery disease and stroke, in particular ischemic stroke. The prevalence of PCSK9 rs505151 and OLR1 rs11053646 variants in ischemic stroke were 0.005 and 0.116, respectively. However, to date, association between OLR1 rs11053646 and PCSK9 rs505151 polymorphisms and the risk of ischemic stroke remains unclear and inconclusive. Therefore, this first meta-analysis was carried out to clarify the presumed influence of genetic polymorphisms on ischemic stroke, by analyzing the complete coverage of all relevant studies. All eligible case-control and cohort studies that met the search term were retrieved in multiple scientific databases. Data of interest such as demographic data and genotyping methods were extracted from each study, and the meta-analysis was performed using RevMan 5.3 and Metafor R 3.2.1. The pooled odd ratios (ORs) and 95% confidence intervals (CIs) were calculated using both fixed- and random-effect models. A total of seven case-control studies encompassing 1897 ischemic stroke cases and 2119 healthy controls were critically evaluated. Pooled results from the genetic models indicated that OLR1 rs11053646 dominant (OR=1.33. 95%CI:1.11-1.58) and co-dominant models (OR=1.24, 95%CI:1.02-1.51) were significantly associated with ischemic stroke. For PCSK9 rs505151 polymorphism, the OR of co-dominant model (OR=1.36, 95%CI:1.01-1.58) was found to be higher among ischemic stroke patients. In conclusion, the current meta-analysis highlighted that variant allele of OLR1 rs11053646 G>C and PCSK9 rs505151 A>G may contribute to the susceptibility risk of ischemic stroke.

Two-level k-means clustering algorithm for k-tau relationship establishment and linear-time classification

Partitional clustering algorithms, which partition the dataset into a pre-defined number of clusters, can be broadly classified into two types: algorithms which explicitly take the number of clusters as input and algorithms that take the expected size of a cluster as input. In this paper, we propose a variant of the k-means algorithm and prove that it is more efficient than standard k-means algorithms. An important contribution of this paper is the establishment of a relation between the number of clusters and the size of the clusters in a dataset through the analysis of our algorithm. We also demonstrate that the integration of this algorithm as a pre-processing step in classification algorithms reduces their running-time complexity.

Mutational analysis of active-site residues in the Mycobacterium leprae RecA intein, a LAGLIDADG homing endonuclease: Asp(122) and Asp(193) are crucial to the double-stranded DNA cleavage activity whereas Asp(218) is not

Mycobacterium leprae recA harbors an in-frame insertion sequence that encodes an intein homing endonuclease (PI-MleI). Most inteins (intein endonucleases) possess two conserved LAGLIDADG (DOD) motifs at their ctive center. A common feature of LAGLIDADG-type homing endonucleases is that they recognize and cleave the same or very similar DNA sequences. However, PI-MleI is distinctive from other members of the family of LAGLIDADG-type HEases for its modular structure with functionally separable domains for DNA-binding and cleavage, each with distinct sequence preferences. Sequence alignment analyses of PI-MleI revealed three putative LAGLIDADG motifs; however, there is conflicting bioinformatics data in regard to their identity and specific location within the intein polypeptide. To resolve this conflict and to determine the active-site residues essential for DNA target site recognition and double-stranded DNA cleavage, we performed site-directed mutagenesis of presumptive catalytic residues in the LAGLIDADG motifs. Analysis of target DNA recognition and kinetic parameters of the wild-type PI-MleI and its variants disclosed that the two amino acid residues, Asp(122) (in Block C) and Asp(193) (in functional Block E), are crucial to the double-stranded DNA endonuclease activity, whereas Asp(218) (in pseudo-Block E) is not. However, despite the reduced catalytic activity, the PI-MleI variants, like the wild-type PI-MleI, generated a footprint of the same length around the insertion site. The D122T variant showed significantly reduced catalytic activity, and D122A and D193A mutations although failed to affect their DNA-binding affinities, but abolished the double-stranded DNA cleavage activity. On the other hand, D122C variant showed approximately twofold higher double-stranded DNA cleavage activity, compared with the wild-type PI-MleI. These results provide compelling evidence that Asp(122) and Asp(193) in DOD motif I and II, respectively, are bona fide active-site residues essential for DNA cleavage activity. The implications of these results are discussed in this report.

Functional diversity of human protein kinase splice variants marks significant expansion of human kinome

Background: Protein kinases are involved in diverse spectrum of cellular processes. Availability of draft version of the human genomic data in the year 2001 enabled recognition of repertoire of protein kinases. However, over the years the human genomic data is being refined and the current release of human genomic data has helped us to recognize a larger repertoire of over 900 human protein kinases represented mainly by splice variants. Results: Many of these identified protein kinases are alternatively spliced products. Interestingly, some of the human kinase splice variants appear to be significantly diverged in terms of their functional properties as represented by incorporation or absence of one or more domains. Many sets of protein kinase splice variants have substantially different domain organization and in a few sets of splice variants kinase domains belong to different subfamilies of kinases suggesting potential participation in different signal transduction pathways. Conclusions: Addition or deletion of a domain between splice variants of multi-domain kinases appears to be a means of generating differences in the functional features of otherwise similar kinases. It is intriguing that marked sequence diversity within the catalytic regions of some of the splice variant kinases result in kinases belonging to different subfamilies. These human kinase splice variants with different functions might contribute to diversity of eukaryotic cellular signaling.

Information Ambiguity

We propose a quantity called information ambiguity that plays the same role in the worst-case information-theoretic nalyses as the well-known notion of information entropy performs in the corresponding average-case analyses. We prove various properties of information ambiguity and illustrate its usefulness in performing the worst-case analysis of a variant of distributed source coding problem.

SPECT-aided diagnostics and genetic risk factors in Parkinson's disease

Hypokinesia, rigidity, tremor, and postural instability are the cardinal symptoms of Parkinson s disease (PD). Since these symptoms are not specific to PD the diagnosis may be uncertain in early PD. Etiology and pathogenesis of PD remain unclear. There is no neuroprotective therapy. Genetic findings are expected to reveal metabolic routes in PD pathogenesis and thereby eventually lead to therapeutic innovations. In this thesis, we first aimed to study the usefulness and accuracy of 123I-b-CIT SPECT in the diagnosis of PD in a consecutive clinic-based material including various movement disorders. We subsequently a genetic project to identify genetic risk factors for sporadic PD using a candidate gene approach in a case-control setting including 147 sporadic PD patients and 137 spouse controls. Dopamine transporter imaging by 123I-b-CIT SPECT could distinguish PD from essential tremor, drug-induced parkinsonism, dystonia and psychogenic parkinsonism. However, b-CIT uptake in Parkinson plus syndromes (PSP and multiple system atrophy) and dementia with Lewy bodies was not significantly different from PD. 123I-b-CIT SPECT could not reliably differentiate PD from vascular parkinsonism. 123I-b-CIT SPECT was 100% sensitive and specific in the diagnosis of PD in patients younger than 55 years but less specific in older patients, due to differential distribution of the above conditions in the younger and older age groups. 123I-b-CIT SPECT correlated with symptoms and detected bilateral nigrostriatal defect in patients whose PD was still in unilateral stage. Thus, in addition to as a differential diagnostic aid, 123I-b-CIT SPECT may be used to detect PD early, even pre-symptomatically in at-risk individuals. 123I-b-CIT SPECT was used to aid in the collection of patients to the genetic studies. In the genetic part of this thesis we found an association between PD and a polymorphic CAG-repeat in POLG1 gene encoding the catalytic subunit of mitochondrial polymerase gamma. The CAG-repeat encodes a polyglutamine tract (polyQ), the two most common lengths of which are 10Q (86-90%) and 11Q. In our Finnish material, the rarer non-10Q or non-11Q length variants (6Q-9Q, 12Q-14Q, 4R+9Q) were more frequent in patients than in spouse controls (10% vs. 3.5 %, p=0.003), or population controls (p=0.001). Therefore, we performed a replication study in 652 North American PD patients and 292 controls. Non-10/11Q alleles were more common in the US PD patients compared to the controls but the difference did not reach statistical significance (p=0.07). This larger data suggested our original definition of variant length allele might need reconsideration. Most previous studies on phenotypic effects of POLG1 polyQ have defined 10Q as the only normal allele. Non-10Q alleles were significantly more common in patients compared to the controls (17.3% vs. 12.3 %, p= 0.005). This association between non-10Q length variants and PD remained significant when compared to a larger set of 1541 literature controls (p=0.00005). In conclusion, POLG1 polyQ alleles other than 10Q may predispose to PD. We did not find association between PD and parkin or DJ-1, genes underlying autosomal recessive parkinsonism. The functional Val158Met polymorphism, which affects the catalytic effect of COMT enzyme, and another coding polymorphism in COMT were not associated with PD in our patient material. The APOE e2/3/4 polymorphism modifies risk for Alzheimer s disease and prognosis of for example brain trauma. APOE promoter and enhancer polymorphisms 219G/T and +113G/C, and APOE e3 haplotypes, have also been shown to modify the risk of Alzheimer s disease but not reported in PD. No association was found between PD and APOE e2/3/4 polymorphism, the promoter or enhancer polymorphisms, or the e3 haplotypes.

Genotype-Phenotype Relationships in Long QT Syndrome: Role of Mental Stress, Adrenergic Activity and a Common KCNH2 Polymorphism

Long QT syndrome is a congenital or acquired arrhythmic disorder which manifests as a prolonged QT-interval on the electrocardiogram and as a tendency to develop ventricular arrhythmias which can lead to sudden death. Arrhythmias often occur during intense exercise and/or emotional stress. The two most common subtypes of LQTS are LQT1, caused by mutations in the KCNQ1 gene and LQT2, caused by mutations in the KCNH2 gene. LQT1 and LQT2 patients exhibit arrhythmias in different types of situations: in LQT1 the trigger is usually vigorous exercise whereas in LQT2 arrhythmia results from the patient being startled from rest. It is not clear why trigger factors and clinical outcome differ from each other in the different LQTS subtypes. It is possible that stress hormones such as catecholamines may show different effects depending on the exact nature of the genetic defect, or sensitivity to catecholamines varies from subject to subject. Furthermore, it is possible that subtle genetic variants of putative modifier genes, including those coding for ion channels and hormone receptors, play a role as determinants of individual sensitivity to life-threatening arrhythmias. The present study was designed to identify some of these risk modifiers. It was found that LQT1 and LQT2 patients show an abnormal QT-adaptation to both mental and physical stress. Furthermore, as studied with epinephrine infusion experiments while the heart was paced and action potentials were measured from the right ventricular septum, LQT1 patients showed repolarization abnormalities which were related to their propensity to develop arrhythmia during intense, prolonged sympathetic tone, such as exercise. In LQT2 patients, this repolarization abnormality was noted already at rest corresponding to their arrhythmic episodes as a result of intense, sudden surges in adrenergic tone, such as fright or rage. A common KCNH2 polymorphism was found to affect KCNH2 channel function as demonstrated by in vitro experiments utilizing mammalian cells transfected with the KCNH2 potassium channel as well as QT-dynamics in vivo. Finally, the present study identified a common β-1-adrenergic receptor genotype that is related a shorter QT-interval in LQT1 patients. Also, it was discovered that compound homozygosity for two common β-adrenergic polymorphisms was related to the occurrence of symptoms in the LQT1 type of long QT syndrome. The studies demonstrate important genotype-phenotype differences between different LQTS subtypes and suggest that common modifier gene polymorphisms may affect cardiac repolarization in LQTS. It will be important in the future to prospectively study whether variant gene polymorphisms will assist in clinical risk profiling of LQTS patients.

Genetic background of bipolar disorder and related cognitive impairments

Bipolar disorder (BP) is a complex psychiatric disorder characterized by episodes of mania and depression. BP affects approximately 1% of the world’s population and shows no difference in lifetime prevalence between males and females. BP arises from complex interactions among genetic, developmental and environmental factors, and it is likely that several predisposing genes are involved in BP. The genetic background of BP is still poorly understood, although intensive and long-lasting research has identified several chromosomal regions and genes involved in susceptibility to BP. This thesis work aims to identify the genetic variants that influence bipolar disorder in the Finnish population by candidate gene and genome-wide linkage analyses in families with many BP cases. In addition to diagnosis-based phenotypes, neuropsychological traits that can be seen as potential endophenotypes or intermediate traits for BP were analyzed. In the first part of the thesis, we examined the role of the allelic variants of the TSNAX/DISC1 gene cluster to psychotic and bipolar spectrum disorders and found association of distinct allelic haplotypes with these two groups of disorders. The haplotype at the 5’ end of the Disrupted-in-Schizophrenia-1 gene (DISC1) was over-transmitted to males with psychotic disorder (p = 0.008; for an extended haplotype p = 0.0007 with both genders), whereas haplotypes at the 3’ end of DISC1 associated with bipolar spectrum disorder (p = 0.0002; for an extended haplotype p = 0.0001). The variants of these haplotypes also showed association with different cognitive traits. The haplotypes at the 5’ end associated with perseverations and auditory attention, while the variants at the 3’ end associated with several cognitive traits including verbal fluency and psychomotor processing speed. Second, in our complete set of BP families with 723 individuals we studied six functional candidate genes from three distinct signalling systems: serotonin-related genes (SLC6A4 and TPH2), BDNF -related genes (BDNF, CREB1 and NTRK2) and one gene related to the inflammation and cytokine system (P2RX7). We replicated association of the functional variant Val66Met of BDNF with BP and better performance in retention. The variants at the 5’ end of SLC6A4 also showed some evidence of association among males (p = 0.004), but the widely studied functional variants did not yield any significant results. A protective four-variant haplotype on P2RX7 showed evidence of association with BP and executive functions: semantic and phonemic fluency (p = 0.006 and p = 0.0003, respectively). Third, we analyzed 23 bipolar families originating from the North-Eastern region of Finland. A genome-wide scan was performed using the 6K single nucleotide polymorphism (SNP) array. We identified susceptibility loci at chromosomes 7q31 with a LOD score of 3.20 and at 9p13.1 with a LOD score of 4.02. We followed up both linkage findings in the complete set of 179 Finnish bipolar families. The finding on chromosome 9p13 was supported (maximum LOD score of 3.02), but the susceptibility gene itself remains unclarified. In the fourth part of the thesis, we wanted to test the role of the allelic variants that have associated with bipolar disorder in recent genome-wide association studies (GWAS). We could confirm findings for the DFNB31, SORCS2, SCL39A3, and DGKH genes. The best signal in this study comes from DFNB31, which remained significant after multiple testing corrections. Two variants of SORCS2 were allelic replications and presented the same signal as the haplotype analysis. However, no association was detected with the PALB2 gene, which was the most significantly associated region in the previous GWAS. Our results indicate that BP is heterogeneous and its genetic background may accordingly vary in different populations. In order to fully understand the allelic heterogeneity that underlies common diseases such as BP, complete genome sequencing for many individuals with and without the disease is required. Identification of the specific risk variants will help us better understand the pathophysiology underlying BP and will lead to the development of treatments with specific biochemical targets. In addition, it will further facilitate the identification of environmental factors that alter risk, which will potentially provide improved occupational, social and psychological advice for individuals with high risk of BP.

Hematopoietic Stem Cell Development and Transcriptional Regulation

The continuous production of blood cells, a process termed hematopoiesis, is sustained throughout the lifetime of an individual by a relatively small population of cells known as hematopoietic stem cells (HSCs). HSCs are unique cells characterized by their ability to self-renew and give rise to all types of mature blood cells. Given their high proliferative potential, HSCs need to be tightly regulated on the cellular and molecular levels or could otherwise turn malignant. On the other hand, the tight regulatory control of HSC function also translates into difficulties in culturing and expanding HSCs in vitro. In fact, it is currently not possible to maintain or expand HSCs ex vivo without rapid loss of self-renewal. Increased knowledge of the unique features of important HSC niches and of key transcriptional regulatory programs that govern HSC behavior is thus needed. Additional insight in the mechanisms of stem cell formation could enable us to recapitulate the processes of HSC formation and self-renewal/expansion ex vivo with the ultimate goal of creating an unlimited supply of HSCs from e.g. human embryonic stem cells (hESCs) or induced pluripotent stem cells (iPS) to be used in therapy. We thus asked: How are hematopoietic stem cells formed and in what cellular niches does this happen (Papers I, II)? What are the molecular mechanisms that govern hematopoietic stem cell development and differentiation (Papers III, IV)? Importantly, we could show that placenta is a major fetal hematopoietic niche that harbors a large number of HSCs during midgestation (Paper I)(Gekas et al., 2005). In order to address whether the HSCs found in placenta were formed there we utilized the Runx1-LacZ knock-in and Ncx1 knockout mouse models (Paper II). Importantly, we could show that HSCs emerge de novo in the placental vasculature in the absence of circulation (Rhodes et al., 2008). Furthermore, we could identify defined microenvironmental niches within the placenta with distinct roles in hematopoiesis: the large vessels of the chorioallantoic mesenchyme serve as sites of HSC generation whereas the placental labyrinth is a niche supporting HSC expansion (Rhodes et al., 2008). Overall, these studies illustrate the importance of distinct milieus in the emergence and subsequent maturation of HSCs. To ensure proper function of HSCs several regulatory mechanisms are in place. The microenvironment in which HSCs reside provides soluble factors and cell-cell interactions. In the cell-nucleus, these cell-extrinsic cues are interpreted in the context of cell-intrinsic developmental programs which are governed by transcription factors. An essential transcription factor for initiation of hematopoiesis is Scl/Tal1 (stem cell leukemia gene/T-cell acute leukemia gene 1). Loss of Scl results in early embryonic death and total lack of all blood cells, yet deactivation of Scl in the adult does not affect HSC function (Mikkola et al., 2003b. In order to define the temporal window of Scl requirement during fetal hematopoietic development, we deactivated Scl in all hematopoietic lineages shortly after hematopoietic specification in the embryo . Interestingly, maturation, expansion and function of fetal HSCs was unaffected, and, as in the adult, red blood cell and platelet differentiation was impaired (Paper III)(Schlaeger et al., 2005). These findings highlight that, once specified, the hematopoietic fate is stable even in the absence of Scl and is maintained through mechanisms that are distinct from those required for the initial fate choice. As the critical downstream targets of Scl remain unknown, we sought to identify and characterize target genes of Scl (Paper IV). We could identify transcription factor Mef2C (myocyte enhancer factor 2 C) as a novel direct target gene of Scl specifically in the megakaryocyte lineage which largely explains the megakaryocyte defect observed in Scl deficient mice. In addition, we observed an Scl-independent requirement of Mef2C in the B-cell compartment, as loss of Mef2C leads to accelerated B-cell aging (Gekas et al. Submitted). Taken together, these studies identify key extracellular microenvironments and intracellular transcriptional regulators that dictate different stages of HSC development, from emergence to lineage choice to aging.

Timing of puberty : Genetic regulation

In the general population, the timing of puberty is normally distributed. This variation is determined by genetic and environmental factors, but the exact mechanisms underlying these influences remain elusive. The purpose of this study was to gain insight into genetic regulation of pubertal timing. Contributions of genetic versus environmental factors to the normal variation of pubertal timing were explored in twins. Familial occurrence and inheritance patterns of constitutional delay of growth and puberty, CDGP (a variant of normal pubertal timing), were studied in pedigrees of patients with this condition. To ultimately detect genes involved in the regulation of pubertal timing, genetic loci conferring susceptibility to CDGP were mapped by linkage analysis in the same family cohort. To subdivide the overall phenotypic variance of pubertal timing into genetic and environmental components, genetic modeling based on monozygous twins sharing 100% and dizygous twins sharing 50% of their genes was used in 2309 girls and 1828 boys from the FinnTwin 12-17 study. The timing of puberty was estimated from height growth, i.e. change in the relative height between the age when pubertal growth velocity peaks in the general population and adulthood. This reflects the percentage of adult height achieved at the average peak height velocity age, and thus, pubertal timing. Boys and girls diagnosed with CDGP were gathered through medical records from six pediatric clinics in Finland. First-degree relatives of the probands were invited to participate by letter; altogether, 286 families were recruited. When possible, families were extended to include also second-, third-, or fourth-degree relatives. The timing of puberty in all family members was primarily assessed from longitudinal growth data. Delayed puberty was defined by onset of pubertal growth spurt or peak height velocity taking place 1.5 (relaxed criterion) or 2 SD (strict criterion) beyond the mean. If growth data were unavailable, pubertal timing was based on interviews. In this case, CDGP criteria were set as having undergone pubertal development more than 2 (strict criterion) or 1.5 years (relaxed criterion) later than their peers, or menarche after 15 (strict criterion) or 14 years (relaxed criterion). Familial occurrence of strict CDGP was explored in families of 124 patients (95 males and 29 females) from two clinics in Southern Finland. In linkage analysis, we used relaxed CDGP criteria; 52 families with solely growth data-based CDGP diagnoses were selected from all clinics. Based on twin data, genetic factors explain 86% and 82% of the variance of pubertal timing in girls and boys, respectively. In families, 80% of male and 76% of female probands had affected first-degree relatives, in whom CDGP was 15 times more common than the expected (2.5%). In 74% (17 of 23) of the extended families with only one affected parent, familial patterns were consistent with autosomal dominant inheritance. By using 383 multiallelic markers and subsequently fine-mapping with 25 additional markers, significant linkage for CDGP was detected to the pericentromeric region of chromosome 2, to 2p13-2q13 (multipoint HLOD 4.44, α 0.41). The findings of the large twin study imply that the vast majority of the normal variation of pubertal timing is attributed to genetic effects. Moreover, the high frequency of dominant inheritance patterns and the large number of affected relatives of CDGP patients suggest that genetic factors also markedly contribute to constitutional delay of puberty. Detection of the locus 2p13-2q13 in the pericentromeric region of chromosome 2 associating with CDGP is one step towards unraveling the genes that determine pubertal timing.

Common and rare variants of the renin-angiotensin system and their relation to antihypertensive drug responses

Most of the diseases affecting public health, like hypertension, are multifactorial by etiology. Hypertension is influenced by genetic, life style and environmental factors. Estimation of the influence of genes to the risk of essential hypertension varies from 30 to 50%. It is plausible that in most of the cases susceptibility to hypertension is determined by the action of more than one gene. Although the exact molecular mechanism underlying essential hypertension remains obscure, several monogenic forms of hypertension have been identified. Since common genetic variations may predict, not only to susceptibility to hypertension, but also response to antihypertensive drug therapy, pharmacogenetic approaches may provide useful markers in finding relations between candidate genes and phenotypes of hypertension. The aim of this study was to identify genetic mutations and polymorphisms contributing to human hypertension, and examine their relationships to intermediate phenotypes of hypertension, such as blood pressure (BP) responses to antihypertensive drugs or biochemical laboratory values. Two groups of patients were investigated in the present study. The first group was collected from the database of patients investigated in the Hypertension Outpatient Ward, Helsinki University Central Hospital, and consisted of 399 subjects considered to have essential hypertension. Frequncies of the mutant or variant alleles were compared with those in two reference groups, healthy blood donors (n = 301) and normotensive males (n = 175). The second group of subjects with hypertension was collected prospectively. The study subjects (n=313) underwent a protocol lasting eight months, including four one-month drug treatment periods with antihypertensive medications (thiazide diuretic, β-blocker, calcium channel antagonist, and an angiotensin II receptor antagonist). BP responses and laboratory values were related to polymorphims of several candidate genes of the renin-angiotensin system (RAS). In addition, two patients with typical features of Liddle’s syndrome were screened for mutations in kidney epithelial sodium channel (ENaC) subunits. Two novel mutations causing Liddle’s syndrome were identified. The first mutation identified located in the beta-subunit of ENaC and the second mutation found located in the gamma-subunit, constituting the first identified Liddle mutation locating in the extracellular domain. This mutation showed 2-fold increase in channel activity in vitro. Three gene variants, of which two are novel, were identified in ENaC subunits. The prevalence of the variants was three times higher in hypertensive patients (9%) than in reference groups (3%). The variant carriers had increased daily urinary potassium excretion rate in relation to their renin levels compared with controls suggesting increased ENaC activity, although in vitro they did not show increased channel activity. Of the common polymorphisms of the RAS studied, angiotensin II receptor type I (AGTR1) 1166 A/C polymorphism was associated with modest changes in RAS activity. Thus, patients homozygous for the C allele tended to have increased aldosterone and decreased renin levels. In vitro functional studies using transfected HEK293 cells provided additional evidence that the AGTR1 1166 C allele may be associated with increased expression of the AGTR1. Common polymorphisms of the alpha-adducin and the RAS genes did not significantly predict BP responses to one-month monotherapies with hydroclorothiazide, bisoprolol, amlodipin, or losartan. In conclusion, two novel mutations of ENaC subunits causing Liddle’s syndrome were identified. In addition, three common ENaC polymorphisms were shown to be associated with occurrence of essential hypertension, but their exact functional and clinical consequences remain to be explored. The AGTR1 1166 C allele may modify the endocrine phenotype of hypertensive patients, when present in homozygous form. Certain widely studied polymorphisms of the ACE, angiotensinogen, AGTR1 and alpha-adducin genes did not significantly affect responses to a thiazide, β-blocker, calcium channel antagonist, and angiotensin II receptor antagonist.

Consistent quaternion interpolation for objective finite element approximation of geometrically exact beam

We explore an isoparametric interpolation of total quaternion for geometrically consistent, strain-objective and path-independent finite element solutions of the geometrically exact beam. This interpolation is a variant of the broader class known as slerp. The equivalence between the proposed interpolation and that of relative rotation is shown without any recourse to local bijection between quaternions and rotations. We show that, for a two-noded beam element, the use of relative rotation is not mandatory for attaining consistency cum objectivity and an appropriate interpolation of total rotation variables is sufficient. The interpolation of total quaternion, which is computationally more efficient than the one based on local rotations, converts nodal rotation vectors to quaternions and interpolates them in a manner consistent with the character of the rotation manifold. This interpolation, unlike the additive interpolation of total rotation, corresponds to a geodesic on the rotation manifold. For beam elements with more than two nodes, however, a consistent extension of the proposed quaternion interpolation is difficult. Alternatively, a quaternion-based procedure involving interpolation of relative rotations is proposed for such higher order elements. We also briefly discuss a strategy for the removal of possible singularity in the interpolation of quaternions, proposed in [I. Romero, The interpolation of rotations and its application to finite element models of geometrically exact rods, Comput. Mech. 34 (2004) 121–133]. The strain-objectivity and path-independence of solutions are justified theoretically and then demonstrated through numerical experiments. This study, being focused only on the interpolation of rotations, uses a standard finite element discretization, as adopted by Simo and Vu-Quoc [J.C. Simo, L. Vu-Quoc, A three-dimensional finite rod model part II: computational aspects, Comput. Methods Appl. Mech. Engrg. 58 (1986) 79–116]. The rotation update is achieved via quaternion multiplication followed by the extraction of the rotation vector. Nodal rotations are stored in terms of rotation vectors and no secondary storages are required.

Molecular Genetics of Age-related Macular Degeneration

Age-related macular degeneration (AMD; OMIM # 603075) is an eye disease of the elderly, signs of which appear after the age of 50. In the Western world it is a leading cause of permanent visual loss with a prevalence of 8.5% in persons under 54 years of age and of 37% in persons over 75 years of age. Early forms of AMD may be asymptomatic, but in the late forms usually a central scotoma in the visual field follows severely complicating daily tasks. Smoking, age, and genetic predisposition are known risk factors for AMD. Until recently no true susceptibility genes had been identified though the composition of drusen deposits, the hallmarks of AMD, has suggested that the complement system might play a role in the pathogenesis of AMD. When four groups reported in March 2005, that, on chromosome 1q32, a Y402H variant in the complement factor H (CFH) gene confers risk for AMD in independent Caucasian samples, a new period in the field of genetic research of AMD started. CFH is a key regulator of the complement system. Thus, it is logical to speculate, that it plays a role in the pathogenesis of AMD. We performed a case-control association study to analyse whether the CFH Y402H variant contain a risk for AMD in the Finnish population. Although the population of Finland represents a genetic isolate, the CFH Y402H polymorphism was associated with AMD also in our patient sample with similar risk allele frequencies as in the other Caucasian populations. We further evaluated the effects of this variant, but no association between lesion subtype (predominantly classic, minimally classic or occult lesion) or lesion size of neovascular AMD and the CFH Y402H variant was detected. Neither did the variant have an effect on the photodynamic therapy (PDT) outcome. The patients that respond to PDT carried the risk genotype as frequently as those who did not respond, and no difference was found in the number of PDT sessions needed in patients with or without the risk genotypes of CFH Y402H. Functional analyses, however, showed that the binding of C-reactive protein (CRP) to CFH was significantly reduced in patients with the risk genotype of Y402H. In the past two years, the LOC387715/ high-temperature requirement factor A1 (HTRA1) locus on 10q26 has also been repeatedly associated with AMD in several populations. The recent discovery of the LOC387715 protein on the mitochondrial outer membrane suggests that the LOC387715 gene, not HTRA1, is the true predisposing gene in this region, although its biological function is still unknown. In our Finnish patient material, patients with AMD carried the A69S risk genotype of LOC387715 more frequently than the controls. Also, for the first time, an interaction between the CFH Y402H and the LOC387715 A69S variants was found. The most recently detected susceptibilty gene of AMD, the complement component 3 (C3) gene, encodes the central component of the complement system, C3. In our Finnish sample, an additive gene effect for the C3 locus was detected, though weaker than the effects for the two main loci, CFH and LOC387715. Instead, the hemicentin-1 or the elongation of very long chain fatty acids-like 4 genes that have also been suggested as candidate genes for AMD did not carry a risk for AMD in the Finnish population. This was the first series of molecular genetic study of AMD in Finland. We showed that two common risk variants, CFH Y402H and LOC387715 A69S, represent a high risk of AMD also in the isolated Finnish population, and furthermore, that they had a statistical interaction. It was demonstrated that the CFH Y402H risk genotype affects the binding of CFH to CRP thus suggesting that complement indeed plays an important role in the pathogenesis of AMD.