Bivariate genome-wide association study of genetically correlated neuroimaging phenotypes from DTI and MRI through a seemingly unrelated regression model

Large multisite efforts (e.g., the ENIGMA Consortium), have shown that neuroimaging traits including tract integrity (from DTI fractional anisotropy, FA) and subcortical volumes (from T1-weighted scans) are highly heritable and promising phenotypes for discovering genetic variants associated with brain structure. However, genetic correlations (rg) among measures from these different modalities for mapping the human genome to the brain remain unknown. Discovering these correlations can help map genetic and neuroanatomical pathways implicated in development and inherited risk for disease. We use structural equation models and a twin design to find rg between pairs of phenotypes extracted from DTI and MRI scans. When controlling for intracranial volume, the caudate as well as related measures from the limbic system - hippocampal volume - showed high rg with the cingulum FA. Using an unrelated sample and a Seemingly Unrelated Regression model for bivariate analysis of this connection, we show that a multivariate GWAS approach may be more promising for genetic discovery than a univariate approach applied to each trait separately.

Genetic influences on brain asymmetry: A DTI study of 374 twins and siblings

Brain asymmetry, or the structural and functional specialization of each brain hemisphere, has fascinated neuroscientists for over a century. Even so, genetic and environmental factors that influence brain asymmetry are largely unknown. Diffusion tensor imaging (DTI) now allows asymmetry to be studied at a microscopic scale by examining differences in fiber characteristics across hemispheres rather than differences in structure shapes and volumes. Here we analyzed 4. Tesla DTI scans from 374 healthy adults, including 60 monozygotic twin pairs, 45 same-sex dizygotic pairs, and 164 mixed-sex DZ twins and their siblings; mean age: 24.4 years ± 1.9 SD). All DTI scans were nonlinearly aligned to a geometrically-symmetric, population-based image template. We computed voxel-wise maps of significant asymmetries (left/right differences) for common diffusion measures that reflect fiber integrity (fractional and geodesic anisotropy; FA, GA and mean diffusivity, MD). In quantitative genetic models computed from all same-sex twin pairs (N=210 subjects), genetic factors accounted for 33% of the variance in asymmetry for the inferior fronto-occipital fasciculus, 37% for the anterior thalamic radiation, and 20% for the forceps major and uncinate fasciculus (all L > R). Shared environmental factors accounted for around 15% of the variance in asymmetry for the cortico-spinal tract (R > L) and about 10% for the forceps minor (L > R). Sex differences in asymmetry (men > women) were significant, and were greatest in regions with prominent FA asymmetries. These maps identify heritable DTI-derived features, and may empower genome-wide searches for genetic polymorphisms that influence brain asymmetry.

Genetic, functional, and histopathological evaluation of two C-terminal BRCA1 missense variants

Background: The vast majority of BRCA1 missense sequence variants remain uncharacterised for their possible effect on protein expression and function, and therefore are unclassified in terms of their pathogenicity. BRCA1 plays diverse cellular roles and it is unlikely that any single functional assay will accurately reflect the total cellular implications of missense mutations in this gene. Objective: To elucidate the effect of two BRCA1 variants, 5236G>C (G1706A) and 5242C>A (A1708E) on BRCA1 function, and to survey the relative usefulness of several assays to direct the characterisation of other unclassified variants in BRCA genes. Methods and Results: Data from a range of bioinformatic, genetic, and histopathological analyses, and in vitro functional assays indicated that the 1708E variant was associated with the disruption of different cellular functions of BRCA1. In transient transfection experiments in T47D and 293T cells, the 1708E product was mislocalised to the cytoplasm and induced centrosome amplification in 293T cells. The 1708E variant also failed to transactivate transcription of reporter constructs in mammalian transcriptional transactivation assays. In contrast, the 1706A variant displayed a phenotype comparable to wildtype BRCA1 in these assays. Consistent with functional data, tumours from 1708E carriers showed typical BRCA1 pathology, while tumour material from 1706A carriers displayed few histopathological features associated with BRCA1 related tumours. Conclusions: A comprehensive range of genetic, bioinformatic, and functional analyses have been combined for the characterisation of BRCA1 unclassified sequence variants. Consistent with the functional analyses, the combined odds of causality calculated for the 1706A variant after multifactorial likelihood analysis (1:142) indicates a definitive classification of this variant as "benign". In contrast, functional assays of the 1708E variant indicate that it is pathogenic, possibly through subcellular mislocalisation. However, the combined odds of 262:1 in favour of causality of this variant does not meet the minimal ratio of 1000:1 for classification as pathogenic, and A1708E remains formally designated as unclassified. Our findings highlight the importance of comprehensive genetic information, together with detailed functional analysis for the definitive categorisation of unclassified sequence variants. This combination of analyses may have direct application to the characterisation of other unclassified variants in BRCA1 and BRCA2.

Modeling strategies in longitudinal data analysis: Covariate, variance function and correlation structure selection

A modeling paradigm is proposed for covariate, variance and working correlation structure selection for longitudinal data analysis. Appropriate selection of covariates is pertinent to correct variance modeling and selecting the appropriate covariates and variance function is vital to correlation structure selection. This leads to a stepwise model selection procedure that deploys a combination of different model selection criteria. Although these criteria find a common theoretical root based on approximating the Kullback-Leibler distance, they are designed to address different aspects of model selection and have different merits and limitations. For example, the extended quasi-likelihood information criterion (EQIC) with a covariance penalty performs well for covariate selection even when the working variance function is misspecified, but EQIC contains little information on correlation structures. The proposed model selection strategies are outlined and a Monte Carlo assessment of their finite sample properties is reported. Two longitudinal studies are used for illustration.

An evaluation of genetic analyses, skull morphology and visual appearance for assessing dingo purity: implications for dingo conservation

The introgression of domestic dog genes into dingo populations threatens the genetic integrity of 'pure' dingoes. However, dingo conservation efforts are hampered by difficulties in distinguishing between dingoes and hybrids in the field. This study evaluates consistency in the status of hybridisation (i.e. dingo, hybrid or dog) assigned by genetic analyses, skull morphology and visual assessments. Of the 56 south-east Queensland animals sampled, 39 (69.6%) were assigned the same status by all three methods, 10 (17.9%) by genetic and skull methods, four (7.1%) by genetic and visual methods; and two (3.6%) by skull and visual methods. Pair-wise comparisons identified a significant relationship between genetic and skull methods, but not between either of these and visual methods. Results from surveying 13 experienced wild dog managers showed that hybrids were more easily identified by visual characters than were dingoes. A more reliable visual assessment can be developed through determining the relationship between (1) genetics and phenotype by sampling wild dog populations and (2) the expression of visual characteristics from different proportions and breeds of domestic dog genes by breeding trials. Culling obvious hybrids based on visual characteristics, such as sable and patchy coat colours, should slow the process of hybridisation.

Familial aggregation of type 1 diabetes and diabetic nephropathy in Finland

Type 1 diabetes (T1D) is a common, multifactorial disease with strong familial clustering. In Finland, the incidence of T1D among children aged 14 years or under is the highest in the world. The increase in incidence has been approximately 2.4% per year. Although most new T1D cases are sporadic the first-degree relatives are at an increased risk of developing the same disease. This study was designed to examine the familial aggregation of T1D and one of its serious complications, diabetic nephropathy (DN). More specifically the study aimed (1) to determine the concordance rates of T1D in monozygotic (MZ) and dizygotic (DZ) twins and to estimate the relative contributions of genetic and environmental factors to the variability in liability to T1D as well as to study the age at onset of diabetes in twins; (2) to obtain long-term empirical estimates of the risk of T1D among siblings of T1D patients and the factors related to this risk, especially the effect of age at onset of diabetes in the proband and the birth cohort effect; (3) to establish if DN is aggregating in a Finnish population-based cohort of families with multiple cases of T1D, and to assess its magnitude and particularly to find out whether the risk of DN in siblings is varying according to the severity of DN in the proband and/or the age at onset of T1D: (4) to assess the recurrence risk of T1D in the offspring of a Finnish population-based cohort of patients with childhood onset T1D, and to investigate potential sex-related effects in the transmission of T1D from the diabetic parents to their offspring as well as to study whether there is a temporal trend in the incidence. The study population comprised of the Finnish Young Twin Cohort (22,650 twin pairs), a population-based cohort of patients with T1D diagnosed at the age of 17 years or earlier between 1965 and 1979 (n=5,144) and all their siblings (n=10,168) and offspring (n=5,291). A polygenic, multifactorial liability model was fitted to the twin data. Kaplan-Meier analyses were used to provide the cumulative incidence for the development of T1D and DN. Cox s proportional hazards models were fitted to the data. Poisson regression analysis was used to evaluate temporal trends in incidence. Standardized incidence ratios (SIRs) between the first-degree relatives of T1D patients and background population were determined. The twin study showed that the vast majority of affected MZ twin pairs remained discordant. Pairwise concordance for T1D was 27.3% in MZ and 3.8% in DZ twins. The probandwise concordance estimates were 42.9% and 7.4%, respectively. The model with additive genetic and individual environmental effects was the best-fitting liability model to T1D, with 88% of the phenotypic variance due to genetic factors. The second paper showed that the 50-year cumulative incidence of T1D in the siblings of diabetic probands was 6.9%. A young age at diagnosis in the probands considerably increased the risk. If the proband was diagnosed at the age of 0-4, 5-9, 10-14, 15 or more, the corresponding 40-year cumulative risks were 13.2%, 7.8%, 4.7% and 3.4%. The cumulative incidence increased with increasing birth year. However, SIR among children aged 14 years or under was approximately 12 throughout the follow-up. The third paper showed that diabetic siblings of the probands with nephropathy had a 2.3 times higher risk of DN compared with siblings of probands free of nephropathy. The presence of end stage renal disease (ESRD) in the proband increases the risk three-fold for diabetic siblings. Being diagnosed with diabetes during puberty (10-14) or a few years before (5-9) increased the susceptibility for DN in the siblings. The fourth paper revealed that of the offspring of male probands, 7.8% were affected by the age of 20 compared with 5.3% of the offspring of female probands. Offspring of fathers with T1D have 1.7 times greater risk to be affected with T1D than the offspring of mothers with T1D. The excess risk in the offspring of male fathers manifested itself through the higher risk the younger the father was when diagnosed with T1D. Young age at onset of diabetes in fathers increased the risk of T1D greatly in the offspring, but no such pattern was seen in the offspring of diabetic mothers. The SIR among offspring aged 14 years or under remained fairly constant throughout the follow-up, approximately 10. The present study has provided new knowledge on T1D recurrence risk in the first-degree relatives and the risk factors modifying the risk. Twin data demonstrated high genetic liability for T1D and increased heritability. The vast majority of affected MZ twin pairs, however, remain discordant for T1D. This study confirmed the drastic impact of the young age at onset of diabetes in the probands on the increased risk of T1D in the first-degree relatives. The only exception was the absence of this pattern in the offspring of T1D mothers. Both the sibling and the offspring recurrence risk studies revealed dynamic changes in the cumulative incidence of T1D in the first-degree relatives. SIRs among the first-degree relatives of T1D patients seems to remain fairly constant. The study demonstrates that the penetrance of the susceptibility genes for T1D may be low, although strongly influenced by the environmental factors. Presence of familial aggregation of DN was confirmed for the first time in a population-based study. Although the majority of the sibling pairs with T1D were discordant for DN, its presence in one sibling doubles and presence of ESRD triples the risk of DN in the other diabetic sibling. An encouraging observation was that although the proportion of children to be diagnosed with T1D at the age of 4 or under is increasing, they seem to have a decreased risk of DN or at least delayed onset.

Male eating disorders and related traits : Genetic epidemiological study in Finnish twins

The aim of this study was to deepen the understanding of eating disorders, body image dissatisfaction and related traits in males by examining the epidemiology and genetic epidemiology of these conditions in representative population-based twin samples. The sample of Study I included adolescent twins from FinnTwin12 cohorts born 1983 87 and assessed by a questionnaire at ages 14 y (N=2070 boys, N=2062 girls) and 17 y (N=1857 boys, N=1984 girls). Samples of Studies II-V consisted of young adult twins born 1974-79 from FinnTwin16 cohorts (Study II N=1245 men, Study III N=724 men, Study IV N=2122 men, Study V N=2426 women and N=1962 men), who were assessed by a questionnaire at the age 22-28 y. In addition, 49 men and 526 women were assessed by a diagnostic interview. The overall response rates for both twin cohorts in all studies were 80-90%. In boys, mainly genetic factors (82%, 95% confidence interval [CI] 72-92) explained the covariation of self-esteem between the ages 14 y and 17 y, whereas in girls, environmental factors (69%, 95% CI 43-93) were the largest contributors. Of young men, 30% experienced high muscle dissatisfaction, while 12% used or had used muscle building supplements and/or anabolic steroids on a regular basis. Muscle dissatisfaction exhibited a robust association with the indicators of mental distress and a genetic component (42%, 95% CI 23-59) for its liability in this population was found. The variation of muscle-building substance use was primarily explained by the environmental factors. The incidence rate of anorexia nervosa in males for the age of 10-24 y was 15.7 (95% CI 6.6-37.8) per 100 000 person-years, and its lifetime prevalence by the young adulthood was 0.24% (95% CI 0.03-0.44). All detected probands with anorexia nervosa had recovered from eating disorders, but suffered from substantial psychiatric comorbidity, which manifested also in their co-twins. Additionally, male co-twins of the probands displayed significant dissatisfaction with body musculature, a male-specific feature of body dysmorphic disorder. All probands were from twin pairs discordant for eating disorders. Of the five male probands with anorexia nervosa, only one was from an opposite-sex twin pair. Among women from the opposite-sex pairs, the prevalence of DSM-IV or broad anorexia nervosa was no significantly different compared to that of the women from monozygotic pairs or from dizygotic same-sex pairs. The prevalence of DSM-IV or broad bulimia nervosa did not differ in opposite- versus same-sex female twin individuals either. In both sexes, the overall profile of indicators on eating disorders was rather similar between individuals from opposite-sex and same-sex pairs. In adolescence, development of self-esteem was differently regulated in boys compared to girls: this finding may have far-reaching implications on the etiology of sex discrepancy of internalizing and externalizing disorders. In young men, muscle dissatisfaction and muscle building supplement/steroid use were relatively common. Muscle dissatisfaction was associated with marked psychological distress such as symptoms of depression and disordered eating. Both genetic and environmental factors explained muscle dissatisfaction in the population, but environmental factors appeared to best explain the use of muscle-building substances. In this study, anorexia nervosa in boys and young men from the general population was more common, transient and accompanied by more substantial co-morbidity than previously thought. Co-twins of the probands with anorexia nervosa displayed significant psychopathology such as male specific symptoms of body dysmorphic disorder, but none of them had had an eating disorder: taken together, these traits are suggestive for an endophenotype of anorexia nervosa in males. Little evidence was found on that the risk for anorexia nervosa, bulimia nervosa, disordered eating or body dissatisfaction were associated with twin zygosity. Thus, it is unlikely that in utero femininization, masculinization or postnatal socialization according to the sex of the co-twin have a major influence on the later development of eating disorders or related traits.

Amyloid diseases at old age : a pathological, epidemiological, and genetic study

Along with the increased life span of individuals, the burden of old age-associated diseases has inevitably increased. Alzheimer s disease (AD), probably the most well known geriatric disease, belongs to the old age-associated amyloid diseases. The purpose of this study was to investigate the frequency, genetic and health-associated risk factors, mutual association, and amyloid proteins in two old age-associated amyloid disorders senile systemic amyloidosis (SSA) and cerebral amyloid angiopathy (CAA) as part of the prospective population-based Vantaa 85+ autopsy study on a Finnish population aged 85 years or more (Studies I-III), completed with a case report on a patient with advanced AGel amyloidosis (Study IV). The numbers of patients investigated in the studies (I-III) were 256, 74, and 63, respectively. The diagnosis and grading of amyloid were based upon histological examination of tissue samples obtained post mortem and stained with Congo red. The amyloid fibril and associated proteins were characterized by immunohistochemical staining methods. The genotype frequencies of 20 polymorphisms in 9 genes and information on health-associated risk factors in subjects with and without SSA and CAA were compared. In a Finnish population ≥ 95 years of age, SSA and CAA occurred in 36% and 49% of the subjects, respectively. In total, two-thirds of these very elderly individuals had SSA, CAA, or both. However, in only 14% of the population these two conditions co-occurred. In subjects 85 years or older, the prevalence of SSA was 25%. In this population, SSA was associated with age at the time of death (p=0.002), myocardial infarctions (MIs; p=0.004), the G/G (Val/Val) genotype of the exon 24 polymorphism in the alpha2-macroglobulin (α2M) gene (p=0.042) and with the H2 haplotype of the tau gene (p=0.016). In contrast, the presence of CAA was strongly associated with APOE e4 (p=0.0003), with histopathological AD (p=0.0005), and with clinical dementia (p=0.01) in both e4+ (p=0.02) and e4- (p=0.06) individuals. Apart from demonstrating the amyloid fibril proteins, complement proteins 3d (C3d) and 9 (C9) were detected in the amyloid deposits of CAA and AGel amyloidosis, and α2M protein was found in fibrous scar tissue close to SSA. In conclusion, this first population based study on SSA shows that both SSA and CAA are common in very elderly individuals. Old age, MIs, the exon 24 polymorphism of the α2M gene, and H1/H2 polymorphism of the tau gene associate with SSA while clinical dementia and APOE ε4 genotype associate with CAA. The high prevalence of CAA, combined with its association with clinical dementia independent of APOE genotype, neuropathological AD, or SSA, also highlights its clinical significance in the very aged, among which the serious end stage complications of CAA, namely multiple infarctions and hemorrhages, are rare. The report on a patient having advanced AGel amyloidosis added knowledge on the disease and showed that this generally benign condition occasionally may lead to death. Further studies are warranted to confirm the findings in other populations. Also, the role of α2M and tau in the pathogenesis of SSA and the involvement of complement in the process of amyloid beta (Aβ) protein elimination from the brain remain to be clarified. Finally, the high prevalence of SSA in the elderly raises the need for prospective clinical studies to define its clinical significance.

Whole-genome sequencing reveals untapped genetic potential in Africa’s indigenous cereal crop sorghum

Sorghum is a food and feed cereal crop adapted to heat and drought and a staple for 500 million of the world’s poorest people. Its small diploid genome and phenotypic diversity make it an ideal C4 grass model as a complement to C3 rice. Here we present high coverage (16-45 × ) resequenced genomes of 44 sorghum lines representing the primary gene pool and spanning dimensions of geographic origin, end-use and taxonomic group. We also report the first resequenced genome of S. propinquum, identifying 8 M high-quality SNPs, 1.9 M indels and specific gene loss and gain events in S. bicolor. We observe strong racial structure and a complex domestication history involving at least two distinct domestication events. These assembled genomes enable the leveraging of existing cereal functional genomics data against the novel diversity available in sorghum, providing an unmatched resource for the genetic improvement of sorghum and other grass species.

Behavioural, population, and genetic processes affecting metapopulation dynamics of the Glanville fritillary butterfly

In my thesis I have been studying the effects of population fragmentation and extinction-recolonization dynamics on genetic and evolutionary processes in the Glanville fritillary butterfly (Melitaea cinxia). By conducting crosses within and among newly-colonized populations and using several fitness measures, I found a strong decrease in fitness following colonization by a few related individuals, and a strong negative relationship between parental relatedness and offspring fitness. Thereafter, I was interested in determining the number and relatedness of individuals colonizing new populations, which I did using a set of microsatellites I had previously developed for this species. Additionally, I am interested in the evolution of key life-history traits. By following the lifetime reproductive success of males emerging at different times in a semi-natural setup, I demonstrated that protandry is adaptive in males, and I was able to rule out, for M. cinxia, alternative incidental hypotheses evoked to explain the evolution of protandry in insects. Finally, in work I did together with Prof. Hanna Kokko, I am proposing bet-hedging as a new mechanism that could explain the evolution of polyandry in M. cinxia.

Genetic studies on Finnish lupus erythematosus patients with cutaneous manifestations

Lupus erythematosus (LE) is a chronic, heterogeneous autoimmune disorder with abnormal immune responses, including production of autoantibodies and immune complexes. Clinical presentations of the disease range from mild cutaneous manifestations to a more generalised systemic involvement of internal organs. Cutaneous (CLE) forms are further subclassified into discoid LE (DLE), subacute cutaneous LE (SCLE) and acute cutaneous lupus erythematosus (ACLE), and may later progress to systemic disease (SLE). Both genetic and environmental factors contribute to the disease, although the precise aetiology is still elusive. Furthermore, complex gene-gene or gene-environment interactions may result in different subphenotypes of lupus. The genetic background of CLE is poorly known and only a few genes are confirmed, while the number of robust genetic associations in SLE exceeds 30. The aim of this thesis was to characterise the recruited patients clinically, and identify genetic variants conferring susceptibility to cutaneous variants of LE. Given that cutaneous and systemic disease may share underlying genetic factors, putative CLE candidate genes for genotyping were selected among those showing strong evidence of association in SLE. The correlation between relevant clinical manifestations and risk genotypes was investigated in order to find specific subphenotype associations. In addition, epistatic interactions in SLE were studied. Finally, the role of tissue degrading matrix metalloproteinases (MMP) in LE tissue injury was explored. These studies were conducted in Finnish case-control and family cohort, and Swedish case-control cohort. The clinical picture of the patients in terms of cutaneous, haematological and immunological findings resembled that described in the contemporary literature. However, the proportion of daily smokers was very high supporting the role of smoking in disease aetiology. The results confirmed that, even though clinically distinct entities, CLE and SLE share predisposing genetic factors. For the first time it was shown that known SLE susceptibility genes IRF5 and TYK2 also increase the risk of CLE. A tendency toward gene-gene interaction between these genes was found in SLE. As a remarkable novel finding, it was observed that ITGAM polymorphisms associated even more strongly to DLE than SLE, and the risk estimates were substantially higher than those reported for SLE. Several other recently identified SLE susceptibility genes showed signs of good or modest association especially in DLE. Subphenotype analyses indicated possible associations to clinical features, but marginally significant results reflected lack of sufficient power for these studies. Thorough immunohistochemical analyses of several MMPs demonstrated a role in epidermal changes and dermal tissue remodelling in diseased skin, and suggested that targeted action using selective MMP inhibitors may reduce lupus-induced damage in inflamed tissues. In conclusion, the results provide an insight into the genetics of CLE and demonstrate that genetic predisposition is at least in part shared between cutaneous and systemic variants of LE. This doctoral study has contributed IRF5, TYK2, ITGAM and several other novel genes to the so far short list of genes implicated in CLE susceptibility. Detailed examination of the function of these genes in CLE pathogenesis warrants further studies. Furthermore, the results support the need of subphenotype analysis with sample sizes large enough to reveal possible specific disease associations in order to better understand the heterogeneous nature and clinical specificities of the disease. Comprehensive analysis of clinical data suggests that smoking is an environmental triggering factor.

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.

Migraine Comorbidities : A Clinical and Molecular Genetic Study

Migraine is a highly prevalent disease, and despite several important breakthroughs there are still a many questions unanswered in the clinical, genetic and pathophysiological aspects of migraine research. Migraine has been linked to several other diseases such as epilepsy and stroke, but there are still unsolved issues concerning the true nature of these associations. Three genes predisposing to hemiplegic migraine and several loci associated to migraine have been identified, but so far no genes responsible for common forms of migraine have been recognized. Triptans have provided an important step in migraine treatment, but their usefulness in rare forms of migraine have been controversial. The Finnish Migraine Gene Project (FMGP) includes more than 1600 families and 7500 individuals. We evaluated comorbidity from 1000 consecutive subjects in the FMGP. To search for novel loci, we performed a genome-wide linkage scan in 36 families with high prevalences of migraine with visual aura. We collected 76 subjects from the FMGP who suffer from hemiplegic migraine and have used triptans. Finally, to study possible links between stroke and migraine we evaluated the prevalence of migraine in subjects with cervical artery dissection (CAD) and healthy controls. Migraine was associated with increased prevalence of allergy, hypotension and psychiatric diseases. Additionally, men suffering from migraine with aura had increased prevalence of epilepsy and stroke. Further evidence of association between migraine and epilepsy was found in our linkage study. The parametric two-point linkage analysis showed significant evidence of linkage between migraine aura and a locus on 9q21-q22. Interestingly, the same locus has been associated with occipitotemporal epilepsy. CAD seems to be a migraine risk factor, and therefore a link between stroke and migraine. Notably, CAD seems to alleviate migraine activity further indicating the association between these two conditions. Despite the contraindications of triptans, it seems that they are safe and effective in the abortive treatment of hemiplegic migraine.

Evolutionary and conservation biology of the Finnish house sparrow

Recently it has been recognized that evolutionary aspects play a major role in conservation issues of a species. In this thesis I have combined evolutionary research with conservation studies to provide new insight into these fields. The study object of this thesis is the house sparrow, a species that has features that makes it interesting for this type of study. The house sparrow has been ubiquitous almost all over the world. Even though being still abundant, several countries have reported major declines. These declines have taken place in a relatively short time covering both urban and rural habitats. In Finland this species has declined by more than two thirds in just over two decades. In addition, as the house sparrow lives only in human inhabited areas it can also raise public awareness to conservation issues. I used both an extensive museum collection of house sparrows collected in 1980s from all over Finland as well as samples collected in 2009 from 12 of the previously collected localities. I used molecular techniques to study neutral genetic variation within and genetic differentiation between the study populations. This knowledge I then combined with data gathered on morphometric measurements. In addition I analyzed eight heavy metals from the livers of house sparrows that lived in either rural or urban areas in the 1980s and evaluated the role of heavy metal pollution as a possible cause of the declines. Even though dispersal of house sparrows is limited I found that just as the declines started in 1980s the house sparrows formed a genetically panmictic population on the scale of the whole Finland. When compared to Norway, where neutral genetic divergence has been found even with small geographic distances, I concluded that this difference would be due to contrasting landscapes. In Finland the landscape is rather homogeneous facilitating the movements of these birds and maintaining gene flow even with the low dispersal. To see whether the declines have had an effect on the neutral genetic variation of the populations I did a comparison between the historical and contemporary genetic data. I showed that even though genetic diversity has not decreased due to the drastic declines the populations have indeed become more differentiated from each other. This shows that even in a still quite abundant species the declines can have an effect on the genetic variation. It is shown that genetic diversity and differentiation may approach their new equilibriums at different rates. This emphasizes the importance of studying both of them and if the latter has increased it should be taken as a warning sign of a possible loss of genetic diversity in the future. One of the factors suggested to be responsible for the house sparrow declines is heavy metal pollution. When studying the livers of house sparrows from 1980s I discovered higher levels of heavy metal concentrations in urban than rural habitats, but the levels of the metals were comparatively low and based on that heavy metal pollution does not seem to be a direct cause for the declines in Finland. However, heavy metals are known to decrease the amount of insects in urban areas and thus in the cities heavy metals may have an indirect effect on house sparrows. Although neutral genetic variation is an important tool for conservation genetics it does not tell the whole story. Since neutral genetic variation is not affected by selection, information can be one-sided. It is possible that even neutral genetic differentiation is low, there can be substantial variation in additive genetic traits indicating local adaptation. Therefore I performed a comparison between neutral genetic differentiation and phenotypic differentiation. I discovered that two traits out of seven are likely to be under directional selection, whereas the others could be affected by random genetic drift. Bergmann s rule may be behind the observed directional selection in wing length and body mass. These results highlight the importance of estimating both neutral and adaptive genetic variation.

Genetic variation in the LDL receptor-related protein 5 (LRP5) gene : Association with bone health and metabolic parameters

Bone mass accrual and maintenance are regulated by a complex interplay between genetic and environmental factors. Recent studies have revealed an important role for the low-density lipoprotein receptor-related protein 5 (LRP5) in this process. The aim of this thesis study was to identify novel variants in the LRP5 gene and to further elucidate the association of LRP5 and its variants with various bone health related clinical characteristics. The results of our studies show that loss-of-function mutations in LRP5 cause severe osteoporosis not only in homozygous subjects but also in the carriers of these mutations, who have significantly reduced bone mineral density (BMD) and increased susceptibility to fractures. In addition, we demonstrated for the first time that a common polymorphic LRP5 variant (p.A1330V) was associated with reduced peak bone mass, an important determinant of BMD and osteoporosis in later life. The results from these two studies are concordant with results seen in other studies on LRP5 mutations and in association studies linking genetic variation in LRP5 with BMD and osteoporosis. Several rare LRP5 variants were identified in children with recurrent fractures. Sequencing and multiplex ligation-dependent probe amplification (MLPA) analyses revealed no disease-causing mutations or whole-exon deletions. Our findings from clinical assessments and family-based genotype-phenotype studies suggested that the rare LRP5 variants identified are not the definite cause of fractures in these children. Clinical assessments of our study subjects with LPR5 mutations revealed an unexpectedly high prevalence of impaired glucose tolerance and dyslipidaemia. Moreover, in subsequent studies we discovered that common polymorphic LRP5 variants are associated with unfavorable metabolic characteristics. Changes in lipid profile were already apparent in pre-pubertal children. These results, together with the findings from other studies, suggest an important role for LRP5 also in glucose and lipid metabolism. Our results underscore the important role of LRP5 not only in bone mass accrual and maintenance of skeletal health but also in glucose and lipid metabolism. The role of LRP5 in bone metabolism has long been studied, but further studies with larger study cohorts are still needed to evaluate the specific role of LRP5 variants as metabolic risk factors.