Molecular epidemiology of human rhinoviruses

The first part of this work investigates the molecular epidemiology of a human enterovirus (HEV), echovirus 30 (E-30). This project is part of a series of studies performed in our research team analyzing the molecular epidemiology of HEV-B viruses. A total of 129 virus strains had been isolated in different parts of Europe. The sequence analysis was performed in three different genomic regions: 420 nucleotides (nt) in the VP4/VP2 capsid protein coding region, the entire VP1 capsid protein coding gene of 876 nt, and 150 nt in the VP1/2A junction region. The analysis revealed a succession of dominant sublineages within a major genotype. The temporally earlier genotypes had been replaced by a genetically homogenous lineage that has been circulating in Europe since the late 1970s. The same genotype was found by other research groups in North America and Australia. Globally, other cocirculating genetic lineages also exist. The prevalence of a dominant genotype makes E-30 different from other previously studied HEVs, such as polioviruses and coxsackieviruses B4 and B5, for which several coexisting genetic lineages have been reported. The second part of this work deals with molecular epidemiology of human rhinoviruses (HRVs). A total of 61 field isolates were studied in the 420-nt stretch in the capsid coding region of VP4/VP2. The isolates were collected from children under two years of age in Tampere, Finland. Sequences from the clinical isolates clustered in the two previously known phylogenetic clades. Seasonal clustering was found. Also, several distinct serotype-like clusters were found to co-circulate during the same epidemic season. Reappearance of a cluster after disappearing for a season was observed. The molecular epidemiology of the analyzed strains turned out to be complex, and we decided to continue our studies of HRV. Only five previously published complete genome sequences of HRV prototype strains were available for analysis. Therefore, all designated HRV prototype strains (n=102) were sequenced in the VP4/VP2 region, and the possibility of genetic typing of HRV was evaluated. Seventy-six of the 102 prototype strains clustered in HRV genetic group A (HRV-A) and 25 in group B (HRV-B). Serotype 87 clustered separately from other HRVs with HEV species D. The field strains of HRV represented as many as 19 different genotypes, as judged with an approximate demarcation of a 20% nt difference in the VP4/VP2 region. The interserotypic differences of HRV were generally similar to those reported between different HEV serotypes (i.e. about 20%), but smaller differences, less than 10%, were also observed. Because some HRV serotypes are genetically so closely related, we suggest that the genetic typing be performed using the criterion "the closest prototype strain". This study is the first systematic genetic characterization of all known HRV prototype strains, providing a further taxonomic proposal for classification of HRV. We proposed to divide the genus Human rhinoviruses into HRV-A and HRV-B. The final part of the work comprises a phylogenetic analysis of a subset (48) of HRV prototype strains and field isolates (12) in the nonstructural part of the genome coding for the RNA-dependent RNA polymerase (3D). The proposed division of the HRV strains in the species HRV-A and HRV-B was also supported by 3D region. HRV-B clustered closer to HEV species B, C, and also to polioviruses than to HRV-A. Intraspecies variation within both HRV-A and HRV-B was greater in the 3D coding region than in the VP4/VP2 coding region, in contrast to HEV. Moreover, the diversity of HRV in 3D exceeded that of HEV. One group of HRV-A, designated HRV-A', formed a separate cluster outside other HRV-A in the 3D region. It formed a cluster also in the capsid region, but located within HRV-A. This may reflect a different evolutionary history of distinct genomic regions among HRV-A. Furthermore, the tree topology within HRV-A in the 3D region differed from that in the VP4/VP2, suggesting possible recombination events in the evolution of the strains. No conflicting phylogenies were observed in any of the 12 field isolates. Possible recombination was further studied using the Similarity and Bootscanning analyses of the complete genome sequences of HRV available in public databases. Evidence for recombination among HRV-A was found, as HRV2 and HRV39 showed higher similarity in the nonstructural part of the genome. Whether HRV2 and HRV39 strains - and perhaps also some other HRV-A strains not yet completely sequenced - are recombinants remains to be determined.

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.

Clinical Phenotype and Genetic Epidemiology of Schizophrenia in a Finnish Isolate

This study is part of the joint project "The Genetic Epidemiology and Molecular Genetics of schizophrenia in Finland" between the Departments of Mental Health and Alcohol Research, and Molecular Medicine at the National Public Health Institute. In the study, we utilized three nationwide health care registers: 1) the Hospital Discharge Register, 2) the Free Medication Register, and 3) the Disability Pension Register, plus the National Population Register, in order to identify all patients with schizophrenia born from 1940 to 1976 (N=33,731) in Finland, and their first degree-relatives. 658 patients with at least one parent born in a homogeneous isolate in northeastern Finland were identified, as well as 4904 familial schizophrenia patients with at least two affected siblings from the whole country. The comparison group was derived from the Health 2000 Study. We collected case records and reassessed the register diagnosis. Were contacted the isolate patients and a random sample of patients from the whole country to make diagnostic clinical interviews and to assess the negative and positive symptoms and signs of schizophrenia. In addition to these patients, we interviewed siblings who were initially healthy according to the Hospital Discharge Register. Of those with a register diagnosis of schizophrenia, schizoaffective or schizophreniform disorder, 69% received a record-based consensus diagnosis and 63% an interview-based diagnosis of schizophrenia. Patients with schizophrenia having first-degree relatives with psychotic disorder had more severe affective flattening and alogia than those who were the only affected individuals in their family. The novel findings were: 1) The prevalence of schizophrenia in the isolate was relatively high based on register (1.5%), case record (0.9-1.3%), and interview (0.7-1.2%) data. 2) Isolate patients, regardless of their familial loading for schizophrenia, had less delusions and hallucinations than the whole country familial patients, which may be related to the genetic homogeneity in the isolate. This phenotype encourages the use of endophenotypes in genetic analyses instead of diagnoses alone. 3) The absence of register diagnosis does not confirm that siblings are healthy, because 7.7% of siblings had psychotic symptoms already before the register diagnoses were identified in 1991. For genetic research, the register diagnosis should therefore be reassessed using either a structured interview or a best- estimate case note consensus diagnosis. Structured clinical interview methods need be considered also in clinical practice.

Intrahepatic cholestasis of pregnancy : genetic background, epidemiology and hepatobiliary consequences

Intrahepatic cholestasis of pregnancy (ICP) is the most common cholestatic liver disease during pregnancy. The reported incidence varies from 0.4 to 15% of full-term pregnancies. The etiology is heterogeneous but familial clustering is known to occur. Here we have studied the genetic background, epidemiology, and long-term hepatobiliary consequences of ICP. In a register-based nation-wide study (n=1 080 310) the incidence of ICP was 0.94% during 1987-2004. A slightly higher incidence, 1.3%, was found in a hospital-based series (n=5304) among women attending the University Hospital of Helsinki in 1992-1993. Of these 16% (11/69) were familial and showed a higher (92%) recurrence rate than the sporadic (40%) cases. In the register-based epidemiological study, advanced maternal age and, to a lesser degree, parity were identified as new risk factors for ICP. The risk was 3-fold higher in women >39 years of age compared to women <30 years. Multiple pregnancy also associated with an elevated risk. In a genetic study we found no association of ICP with the genes regulating bile salt transport (ABCB4, ABCB11 and ATP8B1). The livers of postmenopausal women with a history of ICP tolerated well the short-term exposure to oral and transdermal estradiol, although the doses used were higher than those in routine clinical use. The response of serum levels of sex hormone-binding globulin (SHBG) to oral estradiol was slightly reduced in the ICP group. Transdermal estradiol had no effect on C-reactive protein (CRP) or SHBG. A number of liver and biliary diseases were found to be associated with ICP. Women with a history of ICP showed elevated risks for non-alcoholic liver cirrhosis (8.2 CI 1.9-36), cholelithiasis and cholecystitis (3.7 CI 3.2-4.2), hepatitis C (3.5 CI 1.6-7.6) and non-alcoholic pancreatitis (3.2 CI 1.7-5.7). In conclusion, ICP complicates around 1% of all full-term pregnancies in Finland and its incidence has remained unchanged since 1987. It is familial in 16% of cases with a higher recurrence rate. Although the cause remains unknown, several risk factors, namely advanced maternal age, parity and multiple pregnancies, can be identified. Both oral and transdermal regimens of postmenopausal hormone therapy (HT) are safe for women with a history of ICP when liver function is considered. Some ICP patients are at risk of other liver and biliary diseases and, contrary to what has been thought, a follow-up is warranted.

Data integration from two microarray platforms identifies bi-allelic genetic inactivation of RIC8A in a breast cancer cell line

Background: Using array comparative genomic hybridization (aCGH), a large number of deleted genomic regions have been identified in human cancers. However, subsequent efforts to identify target genes selected for inactivation in these regions have often been challenging. Methods: We integrated here genome-wide copy number data with gene expression data and non-sense mediated mRNA decay rates in breast cancer cell lines to prioritize gene candidates that are likely to be tumour suppressor genes inactivated by bi-allelic genetic events. The candidates were sequenced to identify potential mutations. Results: This integrated genomic approach led to the identification of RIC8A at 11p15 as a putative candidate target gene for the genomic deletion in the ZR-75-1 breast cancer cell line. We identified a truncating mutation in this cell line, leading to loss of expression and rapid decay of the transcript. We screened 127 breast cancers for RIC8A mutations, but did not find any pathogenic mutations. No promoter hypermethylation in these tumours was detected either. However, analysis of gene expression data from breast tumours identified a small group of aggressive tumours that displayed low levels of RIC8A transcripts. qRT-PCR analysis of 38 breast tumours showed a strong association between low RIC8A expression and the presence of TP53 mutations (P = 0.006). Conclusion: We demonstrate a data integration strategy leading to the identification of RIC8A as a gene undergoing a classical double-hit genetic inactivation in a breast cancer cell line, as well as in vivo evidence of loss of RIC8A expression in a subgroup of aggressive TP53 mutant breast cancers.

Cognitive Functioning in Alcohol and Other Substance Use Disorders in Young Adulthood : A Genetic Epidemiological Study

Alcohol and other substance use disorders (SUDs) result in great costs and suffering for individuals and families and constitute a notable public health burden. A multitude of factors, ranging from biological to societal, are associated with elevated risk of SUDs, but at the level of individuals, one of the best predictors is a family history of SUDs. Genetically informative twin and family studies have consistently indicated this familial risk to be mainly genetic. In addition, behavioral and temperamental factors such as early initiation of substance use and aggressiveness are associated with the development of SUDs. These familial, behavioral and temperamental risk factors often co-occur, but their relative importance is not well known. People with SUDs have also been found to differ from healthy controls in various domains of cognitive functioning, with poorer verbal ability being among the most consistent findings. However, representative population-based samples have rarely been used in neuropsychological studies of SUDs. In addition, both SUDs and cognitive abilities are influenced by genetic factors, but whether the co-variation of these traits might be partly explained by overlapping genetic influences has not been studied. Problematic substance use also often co-occurs with low educational level, but it is not known whether these outcomes share part of their underlying genetic influences. In addition, educational level may moderate the genetic etiology of alcohol problems, but gene-environment interactions between these phenomena have also not been widely studied. The incidence of SUDs peaks in young adulthood rendering epidemiological studies in this age group informative. This thesis investigated cognitive functioning and other correlates of SUDs in young adulthood in two representative population-based samples of young Finnish adults, one of which consisted of monozygotic and dizygotic twin pairs enabling genetically informative analyses. Using data from the population-based Mental Health in Early Adulthood in Finland (MEAF) study (n=605), the lifetime prevalence of DSM-IV any substance dependence or abuse among persons aged 21—35 years was found to be approximately 14%, with a majority of the diagnoses being alcohol use disorders. Several correlates representing the domains of behavioral and affective factors, parental factors, early initiation of substance use, and educational factors were individually associated with SUDs. The associations between behavioral and affective factors (attention or behavior problems at school, aggression, anxiousness) and SUDs were found to be largely independent of factors from other domains, whereas daily smoking and low education were still associated with SUDs after adjustment for behavioral and affective factors. Using a wide array of neuropsychological tests in the MEAF sample and in a subsample (n=602) of the population-based FinnTwin16 (FT16) study, consistent evidence of poorer verbal cognitive ability related to SUDs was found. In addition, participants with SUDs performed worse than those without disorders in a task assessing psychomotor processing speed in the MEAF sample, whereas no evidence of more specific cognitive deficits was found in either sample. Biometrical structural equation models of the twin data suggested that both alcohol problems and verbal ability had moderate heritabilities (0.54—0.72), and that their covariation could be explained by correlated genetic influences (genetic correlations -0.20 to -0.31). The relationship between educational level and alcohol problems, studied in the full epidemiological FT16 sample (n=4,858), was found to reflect both genetic correlation and gene-environment interaction. The co-occurrence of low education and alcohol problems was influenced by overlapping genetic factors. In addition, higher educational level was associated with increased relative importance of genetic influences on alcohol problems, whereas environmental influences played a more important role in young adults with lower education. In conclusion, SUDs, especially alcohol abuse and dependence, are common among young Finnish adults. Behavioral and affective factors are robustly related to SUDs independently of many other factors, and compared to healthy peers, young adults who have had SUDs during their life exhibit significantly poorer verbal cognitive ability, and possibly less efficient psychomotor processing. Genetic differences between individuals explain a notable proportion of individual differences in risk of alcohol dependence, verbal ability, and educational level, and the co-occurrence of alcohol problems with poorer verbal cognition and low education is influenced by shared genetic backgrounds. Finally, various environmental factors related to educational level in young adulthood moderate the relative importance of genetic factors influencing the risk of alcohol problems, possibly reflecting differences in social control mechanisms related to educational level.

Dissection of the genetic background of childhood onset progressive myoclonic epilepsies

The progressive myoclonic epilepsies (PMEs) are a clinically and etiologically heterogeneous group of symptomatic epilepsies characterized by myoclonus, tonic-clonic seizures, psychomotor regression and ataxia. Different disorders have been classified as PMEs. Of these, the group of neuronal ceroid lipofuscinoses (NCLs) comprise an entity that has onset in childhood, being the most common cause of neurodegeneration in children. The primary aim of this thesis was to dissect the molecular genetic background of patients with childhood onset PME by studying candidate genes and attempting to identify novel PME-associated genes. Another specific aim was to study the primary protein properties of the most recently identified member of the NCL-causing proteins, MFSD8. To dissect the genetic background of a cohort of Turkish patients with childhood onset PME, a screen of the NCL-associated genes PPT1, TPP1, CLN3, CLN5, CLN6, MFSD8, CLN8 and CTSD was performed. Altogether 49 novel mutations were identified, which together with 56 mutations found by collaborators raised the total number of known NCL mutations to 364. Fourteen of the novel mutations affect the recently identified MFSD8 gene, which had originally been identified in a subset of mainly Turkish patients as the underlying cause of CLN7 disease. To investigate the distribution of MFSD8 defects, a total of 211 patients of different ethnic origins were evaluated for mutations in the gene. Altogether 45 patients from nine different countries were provided with a CLN7 molecular diagnosis, denoting the wide geographical occurrence of MFSD8 defects. The mutations are private with only one having been established by a founder-effect in the Roma population from the former Czechoslovakia. All mutations identified except one are associated with the typical clinical picture of variant late-infantile NCL. To address the trafficking properties of MFSD8, lysosomal targeting of the protein was confirmed in both neuronal and non-neuronal cells. The major determinant for this lysosomal sorting was identified to be an N-terminal dileucine based signal (9-EQEPLL-14), recognized by heterotetrameric AP-1 adaptor proteins, suggesting that MFSD8 takes the direct trafficking pathway en route to the lysosomes. Expression studies revealed the neurons as the primary cell-type and the hippocampus and cerebellar granular cell layer as the predominant regions in which MFSD8 is expressed. To identify novel genes associated with childhood onset PME, a single nucleotide polymorphism (SNP) genomewide scan was performed in three small families and 18 sporadic patients followed by homozygosity mapping to determine the candidate loci. One of the families and a sporadic patient were positive for mutations in PLA2G6, a gene that had previously been shown to cause infantile neuroaxonal dystrophy. Application of next-generation sequencing of candidate regions in the remaining two families led to identification of a homozygous missense mutation in USP19 for the first and TXNDC6 for the second family. Analysis of the 18 sporadic cases mapped the best candidate interval in a 1.5 Mb region on chromosome 7q21. Screening of the positional candidate KCTD7 revealed six mutations in seven unrelated families. All patients with mutations in KCTD7 were reported to have early onset PME, rapid disease progression leading to dementia and no pathologic hallmarks. The identification of KCTD7 mutations in nine patients and the clinical delineation of their phenotype establish KCTD7 as a gene for early onset PME. The findings presented in this thesis denote MFSD8 and KCTD7 as genes commonly associated with childhood onset symptomatic epilepsy. The disease-associated role of TXNDC6 awaits verification through identification of additional mutations in patients with similar phenotypes. Completion of the genetic spectrum underlying childhood onset PMEs and understanding of the gene products functions will comprise important steps towards understanding the underlying pathogenetic mechanisms, and will possibly shed light on the general processes of neurodegeneration and nervous system regulation, facilitating the diagnosis, classification and possibly treatment of the affected cases.