Molecular Genetics of Schizophrenia and Related Intermediate Phenotypes in a Founder Population

Schizophrenia, affecting about 1% of population worldwide, is a severe mental disorder characterized by positive and negative symptoms, such as psychosis and anhedonia, as well as cognitive deficits. At present, schizophrenia is considered a complex disorder of neurodevelopmental origin with both genetic and environmental factors contributing to its onset. Although a number of candidate genes for schizophrenia have been highlighted, only very few schizophrenia patients are likely to share identical genetic liability. This study is based on the nation-wide schizophrenia family sample of the National Institute for Health and Welfare, and represents one of the largest and most well-characterized familial series in the world. In the first part of this study, we investigated the roles of the DTNBP1, NRG1, and AKT1 genes in the background of schizophrenia in Finland. Although these genes are associated with schizophrenia liability in several populations, any significant association with clinical diagnostic information of schizophrenia remained absent in our sample of 441 schizophrenia families. In the second part of this study, we first replicated schizophrenia linkage on the long arm of chromosome 7 in 352 schizophrenia families. In the following association analysis, we utilized additional clinical disorder features and intermediate phenotypes – endophenotypes - in addition to diagnostic information from altogether 290 neuropsychologically assessed schizophrenia families. An intragenic short tandem repeat allele of the regional RELN gene, supposed to play a role in the background of several neurodevelopmental disorders, showed significant association with poorer cognitive functioning and more severe schizophrenia symptoms. Additionally, this risk allele was significantly more prevalent among the individuals affected with schizophrenia spectrum disorders. We have previously identified linkage of schizophrenia and its cognitive endophenotypes on the long arms of chromosomes 2, 4, and 5. In the last part of this study, we selected altogether 104 functionally relevant candidate genes from the linked regions. We detected several promising associations, of which especially interesting are the ERBB4 gene, showing association with the severity of schizophrenia symptoms and impairments in traits related to verbal abilities, and the GRIA1 gene, showing association with the severity of schizophrenia symptoms. Our results extend the previous evidence that the genetic risk for schizophrenia is at least partially mediated via the effects of the candidate genes and their combinations on relevant brain systems, resulting in alterations in different disorder domains, such as the cognitive deficits.

α-Synuclein pathology in very elderly Finns : A population-based clinico-neuropathologic and molecular genetic study of Dementia with Lewy bodies

Populations in developed countries are ageing fast. The elderly have the greatest incidence of de-mentia, and thus the increase in the number of demented individuals, increases the immediate costs for the governments concerning healthcare and hospital treatment. Attention is being paid to disorders behind cognitive impairment with behavioural and psychological symptoms, which are enormous contributors to the hospital care required for the elderly. The highest dreams are in prevention; however, before discovering the tools for preventing dementia, the pathogenesis behind dementia disorders needs to be understood. Dementia with Lewy bodies (DLB), a relatively recently discovered dementia disorder compared to Alzheimer’s disease (AD), is estimated to account for up to one third of primary degenerative dementia, thus being the second most common cause of dementia in the elderly. Nevertheless, the impact of neuropathological and genetic findings on the clinical syndrome of DLB is not fully established. In this present series of studies, the frequency of neuropathological findings of DLB and its relation to the clinical findings was evaluated in a cohort of subjects with primary degenerative dementia and in a population-based prospective cohort study of individuals aged 85 years or older. α-synuclein (αS) immunoreactive pathology classifiable according to the DLB consensus criteria was found in one fourth of the primary degenerative dementia subjects. In the population-based study, the corresponding figure was one third of the population, 38% of the demented and one fifth of the non-demented very elderly Finns. However, in spite of the frequent discovery of αS pathology, its association with the clinical symptoms was quite poor. Indeed, the common clinical features of DLB, hypokinesia and visual hallucinations, associated better with the severe neurofibrillary AD-type pathology than with the extensive (diffuse neocortical) αS pathology when both types of pathology were taken into account. The severity of the neurofibrillary AD-type pathology (Braak stage) associated with the extent of αS pathology in the brain. In addition, the genetic study showed an interaction between tau and αS; common variation in the αS gene (SNCA) associated significantly with the severity of the neurofibrillary AD-type pathology and nominally significantly with the extensive αS pathology. Further, the relevance and temporal course of the substantia nigra (SN) degeneration and of the spinal cord αS pathology were studied in relation to αS pathology in the brain. The linear association between the extent of αS pathology in the brain and the neuron loss in SN suggests that in DLB the degeneration of SN proceeds as the αS pathology extends from SN to the neocortex instead of early destruction of SN seen in Parkinson’s disease (PD). Furthermore, the extent of αS pathology in the brain associated with the severity of αS pathology in the thoracic and sacral autonomic nuclei of the spinal cord. The thoracic αS pathology was more common and more severe compared to sacral cord, suggesting that the progress of αS pathology proceeds downwards from the brainstem towards the sacral spinal cord.

Insights into the molecular genetics of celiac disease : Applying family- and population-based strategies

Celiac disease, or gluten intolerance, is triggered by dietary glutens in genetically susceptible individuals and it affects approximately 1% of the Caucasian population. The best known genetic risk factors for celiac disease are HLA DQ2 and DQ8 heterodimers, which are necessary for the development of the disease. However, they alone are not sufficient for disease induction, other risk factors are required. This thesis investigated genetic factors for celiac disease, concentrating on susceptibility loci on chromosomes 5q31-q33, 19p13 and 2q12 previously reported in genome-wide linkage and association studies. In addition, a novel genotyping method for the detection of HLA DQ2 and DQ8 coding haplotypes was validated. This study was conducted using Finnish and Hungarian family materials, and Finnish, Hungarian and Italian case-control materials. Genetic linkage and association were analysed in these materials using candidate gene and fine-mapping approaches. The results confirmed linkage to celiac disease on the chromosomal regions 5q31-q33 and 19p13. Fine-mapping on chromosome 5q31-q33 revealed several modest associations in the region, and highlighted the need for further investigations to locate the causal risk variants. The MYO9B gene on chromosome 19p13 showed evidence for linkage and association particularly with dermatitis herpetiformis, the skin manifestation of celiac disease. This implies a potential difference in the genetic background of the intestinal and skin forms of the disease, although studies on larger samplesets are required. The IL18RAP locus on chromosome 2q12, shown to be associated with celiac disease in a previous genome-wide association study and a subsequent follow-up, showed association in the Hungarian population in this study. The expression of IL18RAP was further investigated in small intestinal tissue and in peripheral blood mononuclear cells. The results showed that IL18RAP is expressed in the relevant tissues. Two putative isoforms of IL18RAP were detected by Western blot analysis, and the results suggested that the ratios and total levels of these isoforms may contribute to the aetiology of celiac disease. A novel genotyping method for celiac disease-associated HLA haplotypes was also validated in this thesis. The method utilises single-nucleotide polymorphisms tagging these HLA haplotypes with high sensitivity and specificity. Our results suggest that this method is transferable between populations, and it is suitable for large-scale analysis. In conclusion, this doctorate study provides an insight into the roles of the 5q31-q33, MYO9B, IL18RAP and HLA loci in the susceptibility to celiac disease in the Finnish, Hungarian and Italian populations, highlighting the need for further studies at these genetic loci and examination of the function of the candidate genes.

Replication of the association of common rs9939609 variant of FTO with increased BMI in an Australian adult twin population but no evidence for gene by environment (G x E) interaction

OBJECTIVE: To further investigate a common variant (rs9939609) in the fat mass- and obesity-associated gene (FTO), which recent genome-wide association studies have shown to be associated with body mass index (BMI) and obesity. DESIGN: We examined the effect of this FTO variant on BMI in 3353 Australian adult male and female twins. RESULTS: The minor A allele of rs9939609 was associated with an increased BMI (P=0.0007). Each additional copy of the A allele was associated with a mean BMI increase of approximately 1.04 kg/m(2) (approximately 3.71 kg). Using variance components decomposition, we estimate that this single-nucleotide polymorphism accounts for approximately 3% of the genetic variance in BMI in our sample (approximately 2% of the total variance). By comparing intrapair variances of monozygotic twins of different genotypes we were able to perform a direct test of gene by environment (G x E) interaction in both sexes and gene by parity (G x P) interaction in women, but no evidence was found for either. CONCLUSIONS: In addition to supporting earlier findings that the rs9939609 variant in the FTO gene is associated with an increased BMI, our results indicate that the associated genetic effect does not interact with environment or parity.

Genetic Heterogeneity in Autism Spectrum Disorders in a Population Isolate

Positional cloning has enabled hypothesis-free, genome-wide scans for genetic factors contributing to disorders or traits. Traditionally linkage analysis has been used to identify regions of interest, followed by meticulous fine mapping and candidate gene screening using association methods and finally sequencing of regions of interest. More recently, genome-wide association analysis has enabled a more direct approach to identify specific genetic variants explaining a part of the variance of the phenotype of interest. Autism spectrum disorders (ASDs) are a group of childhood onset neuropsychiatric disorders with shared core symptoms but varying severity. Although a strong genetic component has been established in ASDs, genetic susceptibility factors have largely eluded characterization. Here, we have utilized modern molecular genetic methods combined with the advantages provided by the special population structure in Finland to identify genetic risk factors for ASDs. The results of this study show that numerous genetic risk factors exist for ASDs even within a population isolate. Stratification based on clinical phenotype resulted in encouraging results, as previously identified linkage to 3p14-p24 was replicated in an independent family set of families with Asperger syndrome, but no other ASDs. Fine-mapping of the previously identified linkage peak for ASDs at 3q25-q27 revealed association between autism and a subunit of the 5-hydroxytryptamine receptor 3C (HTR3C). We also used dense, genome-wide single nucleotide polymorphism (SNP) data to characterize the population structure of Finns. We observed significant population substructure which correlates with the known history of multiple consecutive bottle-necks experienced by the Finnish population. We used this information to ascertain a genetically homogenous subset of autism families to identify possible rare, enriched risk variants using genome-wide SNP data. No rare enriched genetic risk factors were identified in this dataset, although a subset of families could be genealogically linked to form two extended pedigrees. The lack of founder mutations in this isolated population suggests that the majority of genetic risk factors are rare, de novo mutations unique to individual nuclear families. The results of this study are consistent with others in the field. The underlying genetic architecture for this group of disorders appears highly heterogeneous, with common variants accounting for only a subset of genetic risk. The majority of identified risk factors have turned out to be exceedingly rare, and only explain a subset of the genetic risk in the general population in spite of their high penetrance within individual families. The results of this study, together with other results obtained in this field, indicate that family specific linkage, homozygosity mapping and resequencing efforts are needed to identify these rare genetic risk factors.

Helicobacter pylori in Vammala, Finland: seroepidemiological studies and a population-based ‘screen-and-treat’ programme

Helicobacter pylorin (helikobakteeri) tartunta saadaan yleensä lapsena ja tauti jää tavallisesti pysyväksi ilman täsmähoitoa. Onnistunut hoito parantaa pysyvästi helikobakteerista aiheutuvan mahan haavataudin ja näyttää ehkäisevän mahalaukun pahanlaatuisten muutosten kehittymistä. Aloitimme Vammalassa terveyskeskuksessa toteutetun kansainvälisesti ainutlaatuisen väestöpohjaisen helikobakteeritulehduksen seulonta- ja hoito-ohjelman pilottitutkimuksella 1994. 1996 kaikki 15-40-vuotiaat ja 1997-2000 15- ja 45-vuotiaat vammalalaiset kutsuttiin verinäyteseulontaan. Yhteensä 4626 henkilöä (75% kutsutuista) osallistui seulontaan. Vasta-ainepositiivisille tarjottiin helikobakteeritulehduksen lopettava lääkekuuri. Toiminnan seurauksena helikobakteeritulehduksen esiintyvyyden laskettiin vähentyneen 12%:sta 4%:iin 15-40-vuotiaiden ikäryhmässä. Tutkimme myös helikobakteerivasta-ainepositiivisten ja -negatiivisten eroja sekä helikobakteeritulehduksen riskitekijöitä kyselytutkimuksella. Lapsuudenkodin asumisahtauden, äidin matalan koulutusasteen, tupakoinnin, alkoholinkäytön, huonojen asunto-olojen ja ylävatsavaivoista johtuvien sairauslomien todettiin liittyvän helikobakteeritulehdukseen monimuuttuja-analyysissa. Tutkimme seulontaohjelmassa käyttämiemme IgG- ja IgA-luokan helikobakteeri-vasta-ainetestien luotettavuutta eri ikäryhmissä ottaen huomioon atrofisen gastriitin esiintyvyyden. 561 kliinisin perustein gastroskopoidun potilaan aineistossa IgG-testi osoittautui erittäin herkäksi kaikissa ikäryhmissä (99%). Tarkkuus oli myös vanhemmissa ikäryhmissä hyvä (97-93%), kun atrofista gastriittia sairastavat suljettiin pois. IgA- ja CagA-helikobakteerivasta-aineiden on todettu liittyvän lisääntyneeseen mahahaava- ja mahasyöpäriskiin. Analysoimme 560 henkilön pariseeruminäytteet, jotka oli otettu kahden vuosikymmenen välein, ja totesimme, että IgA-vasta-aineiden esiintyvyyden lisääntyyminen iän myötä johtuu paitsi syntymäajankohdasta ja uusista infektioista myös IgA-vasta-ainetasojen kohoamisesta helikobakteeritulehduksen aikana. Selvitimme myös CagA-vasta-ainetasojen muuttumista analysoimalla seeruminäytteet, jotka oli otettu kahden vuosikymmenen välein. Totesimme, että samanaikaisesti kun helikobakteerin esiintyvyys väestössä on alentunut, erityisesti CagA-positiiviset infektiot ovat vähentyneet. Tutkimuksemme osoittaa, että Suomessa terveyskeskuksen yhteydessä voidaan toteuttaa näin laajamittainen seulonta- ja hoito-ohjelma, johon suomalaiset osallistuvat aktiivisesti. Nähtäväksi jää, kuinka paljon ohjelma kykeni vähentämään helikobakteeritulehdukseen liittyviä myöhäisseuraamuksia.

Multiple host use by a sap-sucking membracid: population consequences of nymphal development on primary and secondary host plant species

Aconophora compressa is a gregarious, sap-sucking insect that uses multiple host plant species. Nymphal host plant species (and variety) significantly affected nymphal survival, nymphal development rate and the subsequent size and fecundity of adults, with fiddlewood ( Citharexylum spinosum ) being significantly best in all respects. Nymphs that developed on a relatively poor host ( Duranta erecta var “geisha girl”) and which were moved to fiddlewood as adults laid significantly fewer eggs (mean ± SE = 836 ± 130) than those that developed solely on fiddlewood (1,329 ± 105). Adults on geisha girl, regardless of having been reared as nymphs on fiddlewood or geisha girl, laid significantly fewer eggs (342 ± 83 and 317 ± 74, respectively) than adults on fiddlewood. A simple model that incorporates host plant related survival, development rate and fecundity suggests that the population dynamics of A. compressa are governed mainly by fiddlewood, the primary host. The results have general implications for understanding the population dynamics of herbivores that use multiple host plant species, and also for the way in which weed biological control host testing methods should be conducted.

Raw data for "Linkage disequilibrium estimation of effective population size with immigrants from divergent populations: a case study on Spanish mackerel (Scomberomorus commerson)."

Abstract of Macbeth, G. M., Broderick, D., Buckworth, R. & Ovenden, J. R. (In press, Feb 2013). Linkage disequilibrium estimation of effective population size with immigrants from divergent populations: a case study on Spanish mackerel (Scomberomorus commerson). G3: Genes, Genomes and Genetics. Estimates of genetic effective population size (Ne) using molecular markers are a potentially useful tool for the management of endangered through to commercial species. But, pitfalls are predicted when the effective size is large, as estimates require large numbers of samples from wild populations for statistical validity. Our simulations showed that linkage disequilibrium estimates of Ne up to 10,000 with finite confidence limits can be achieved with sample sizes around 5000. This was deduced from empirical allele frequencies of seven polymorphic microsatellite loci in a commercially harvested fisheries species, the narrow barred Spanish mackerel (Scomberomorus commerson). As expected, the smallest standard deviation of Ne estimates occurred when low frequency alleles were excluded. Additional simulations indicated that the linkage disequilibrium method was sensitive to small numbers of genotypes from cryptic species or conspecific immigrants. A correspondence analysis algorithm was developed to detect and remove outlier genotypes that could possibly be inadvertently sampled from cryptic species or non-breeding immigrants from genetically separate populations. Simulations demonstrated the value of this approach in Spanish mackerel data. When putative immigrants were removed from the empirical data, 95% of the Ne estimates from jacknife resampling were above 24,000.

Soil carbon modelling as a tool for carbon balance studies in forestry

Soils represent a remarkable stock of carbon, and forest soils are estimated to hold half of the global stock of soil carbon. Topical concern about the effects of climate change and forest management on soil carbon as well as practical reporting requirements set by climate conventions have created a need to assess soil carbon stock changes reliably and transparently. The large spatial variability of soil carbon commensurate with relatively slow changes in stocks hinders the assessment of soil carbon stocks and their changes by direct measurements. Models therefore widely serve to estimate carbon stocks and stock changes in soils. This dissertation aimed to develop the soil carbon model YASSO for upland forest soils. The model was aimed to take into account the most important processes controlling the decomposition in soils, yet remain simple enough to ensure its practical applicability in different applications. The model structure and assumptions were presented and the model parameters were defined with empirical measurements. The model was evaluated by studying the sensitivities of the model results to parameter values, by estimating the precision of the results with an uncertainty analysis, and by assessing the accuracy of the model by comparing the predictions against measured data and to the results of an alternative model. The model was applied to study the effects of intensified biomass extraction on the forest carbon balance and to estimate the effects of soil carbon deficit on net greenhouse gas emissions of energy use of forest residues. The model was also applied in an inventory based method to assess the national scale forest carbon balance for Finland’s forests from 1922 to 2004. YASSO managed to describe sufficiently the effects of both the variable litter and climatic conditions on decomposition. When combined with the stand models or other systems providing litter information, the dynamic approach of the model proved to be powerful for estimating changes in soil carbon stocks on different scales. The climate dependency of the model, the effects of nitrogen on decomposition and forest growth as well as the effects of soil texture on soil carbon stock dynamics are areas for development when considering the applicability of the model to different research questions, different land use types and wider geographic regions. Intensified biomass extraction affects soil carbon stocks, and these changes in stocks should be taken into account when considering the net effects of forest residue utilisation as energy. On a national scale, soil carbon stocks play an important role in forest carbon balances.

Modeling population growth and site specific control of the invasive Lantana camara L. (Verbenaceae) under differing fire regimes

It is at the population level that an invasion either fails or succeeds. Lantana camara L. (Verbenaceae) is a weed of great significance in Queensland Australia and globally but its whole life-history ecology is poorly known. Here we used 3 years of field data across four land use types (farm, hoop pine plantation and two open eucalyptus forests, including one with a triennial fire regime) to parameterise the weed’s vital rates and develop size-structured matrix models. Lantana camara in its re-colonization phase, as observed in the recently cleared hoop pine plantation, was projected to increase more rapidly (annual growth rate, λ = 3.80) than at the other three sites (λ 1.88–2.71). Elasticity analyses indicated that growth contributed more (64.6 %) to λ than fecundity (18.5 %) or survival (15.5 %), while across size groups, the contribution was of the order: juvenile (19–27 %) ≥ seed (17–28 %) ≥ seedling (16–25 %) > small adult (4–26 %) ≥ medium adult (7–20 %) > large adult (0–20 %). From a control perspective it is difficult to determine a single weak point in the life cycle of lantana that might be exploited to reduce growth below a sustaining rate. The triennial fire regime applied did not alter the population elasticity structure nor resulted in local control of the weed. However, simulations showed that, except for the farm population, periodic burning could work within 4–10 years for control of the weed, but fire frequency should increase to at least once every 2 years. For the farm, site-specific control may be achieved by 15 years if the biennial fire frequency is tempered with increased burning intensity.

The decline of northern malaria and population dynamics of Plamodium vivax

Europe was declared malaria free in 1975. The disappearance of malaria has traditionally been attributed to numerous deliberate actions like vector control, the screening of houses, more efficient medication etc. Malaria, however, disappeared from many countries like Finland before any counter measures had even started. The aim of this thesis is to study the population ecology of P. vivax and its interaction with the human host and the vector. By finding the factors that attributed to the extinction of vivax malaria it might be possible to improve the modern strategy against P. vivax. The parasite was studied with data from Finland, which provides the longest time series (1749-2008) of malaria statistics in the world. The malaria vectors, Anopheles messeae and A. beklemishevi are still common species in the country. The eradication of vivax malaria is difficult because the parasite has a dormant stage that can cause a relapse long after a primary infection. It was now shown that P. vivax is able to detect the presence of a potential vector. A dormant stage is triggered even from a bite of an uninfected Anopheles mosquito. This optimizes the chances for the Plasmodium to reach a mosquito vector for sexual reproduction. The longevity of the dormant stage could be shown to be at least nine years. The parasite spends several years in its human host and the behaviour of the human carrier had a profound impact on the decline of the disease in Finland. Malaria spring epidemics could be explained by a previous warm summer. Neither annual nor summer mean temperature had any impact on the long term malaria trend. Malaria disappeared slowly from Finland without mosquito control. The sociological change from extended families to nuclear families led to decreased household size. The decreased household size correlated strongly with the decline of malaria. That led to an increased isolation of the subpopulations of P. vivax. Their habitat consisted of the bedrooms in which human carriers slept together with the overwintering vectors. The isolation of the parasite ultimately led to the extinction of vivax malaria. Metapopulation models adapted to local conditions should therefore be implemented as a tool for settlement planning and socio-economic development and become an integrated part of the fight against malaria.

Migraine with aura and migraine without aura are not distinct entities: further evidence from a large Dutch population study

It is often debated whether migraine with aura (MA) and migraine without aura (MO) are etiologically distinct disorders. A previous study using latent class analysis (LCA) in Australian twins showed no evidence for separate subtypes of MO and MA. The aim of the present study was to replicate these results in a population of Dutch twins and their parents, siblings and partners (N = 10,144). Latent class analysis of International Headache Society (IHS)-based migraine symptoms resulted in the identification of 4 classes: a class of unaffected subjects (class 0), a mild form of nonmigrainous headache (class 1), a moderately severe type of migraine (class 2), typically without neurological symptoms or aura (8% reporting aura symptoms), and a severe type of migraine (class 3), typically with neurological symptoms, and aura symptoms in approximately half of the cases. Given the overlap of neurological symptoms and nonmutual exclusivity of aura symptoms, these results do not support the MO and MA subtypes as being etiologically distinct. The heritability in female twins of migraine based on LCA classification was estimated at .50 (95% confidence intervals [CI] .27 - .59), similar to IHS-based migraine diagnosis (h2 = .49, 95% CI .19-.57). However, using a dichotomous classification (affected-unaffected) decreased heritability for the IHS-based classification (h2 = .33, 95% CI .00-.60), but not the LCA-based classification (h2 = .51, 95% CI .23-.61). Importantly, use of the LCA-based classification increased the number of subjects classified as affected. The heritability of the screening question was similar to more detailed LCA and IHS classifications, suggesting that the screening procedure is an important determining factor in genetic studies of migraine.