Monte Carlo Methods in Parameter Estimation of Nonlinear Models

Yksi keskeisimmistä tehtävistä matemaattisten mallien tilastollisessa analyysissä on mallien tuntemattomien parametrien estimointi. Tässä diplomityössä ollaan kiinnostuneita tuntemattomien parametrien jakaumista ja niiden muodostamiseen sopivista numeerisista menetelmistä, etenkin tapauksissa, joissa malli on epälineaarinen parametrien suhteen. Erilaisten numeeristen menetelmien osalta pääpaino on Markovin ketju Monte Carlo -menetelmissä (MCMC). Nämä laskentaintensiiviset menetelmät ovat viime aikoina kasvattaneet suosiotaan lähinnä kasvaneen laskentatehon vuoksi. Sekä Markovin ketjujen että Monte Carlo -simuloinnin teoriaa on esitelty työssä siinä määrin, että menetelmien toimivuus saadaan perusteltua. Viime aikoina kehitetyistä menetelmistä tarkastellaan etenkin adaptiivisia MCMC menetelmiä. Työn lähestymistapa on käytännönläheinen ja erilaisia MCMC -menetelmien toteutukseen liittyviä asioita korostetaan. Työn empiirisessä osuudessa tarkastellaan viiden esimerkkimallin tuntemattomien parametrien jakaumaa käyttäen hyväksi teoriaosassa esitettyjä menetelmiä. Mallit kuvaavat kemiallisia reaktioita ja kuvataan tavallisina differentiaaliyhtälöryhminä. Mallit on kerätty kemisteiltä Lappeenrannan teknillisestä yliopistosta ja Åbo Akademista, Turusta.

Identification and characterization of a resistance gene analog (RGA) from the Caricaceae Dumort family

The majority of cloned resistance (R) genes characterized so far contain a nucleotide-binding site (NBS) and a leucine-rich repeat (LRR) domain, where highly conserved motifs are found. Resistance genes analogs (RGAs) are genetic markers obtained by a PCR-based strategy using degenerated oligonucleotide primers drawn from these highly conserved "motifs". This strategy has the advantage of the high degree of structural and amino acid sequence conservation that is observed in R genes. The objective of the present study was to search for RGAs in Carica papaya L. and Vasconcellea cauliflora Jacq. A. DC. Out of three combinations of primers tested, only one resulted in amplification. The amplified product was cloned in pCR2.1TOPO and than sequenced using M13 forward and reverse primers. Forty-eight clones were sequenced from each species. The 96 sequences generated for each species were cleaned of vector sequences and clustered using CAP3 assembler. From the GENEBANK, one RGA was identified in C. papaya showing a BlastX e-value of 2x10-61 to the gb|AAP45165.1| putative disease resistant protein RGA3 (Solanum bulbocastanum). To the extent of our knowledge this is the first report of a RGA in the Caricaceae Dumort family. Preliminary structural studies were performed to further characterize this putative NBS-LRR type protein. Efforts to search for other RGAs in papaya should continue, mostly to provide basis for the development of transgenic papaya with resistance to diseases.

Autosomal dominant dopa-responsive parkinsonism in a multigenerational Swiss family.

AIM: To describe a large family with autosomal dominant parkinsonism. BACKGROUND: Seven genes are directly implicated in autosomally inherited parkinsonism. However, there are several multigenerational large families known with no identifiable mutation. MATERIAL AND METHODS: Family members were evaluated clinically, by history and chart review. Genetic investigation included SCA2, SCA3, UCHL1, SNCA, LRRK2, PINK1, PRKN, PGRN, FMR1 premutation, and MAPT. The proband underwent brain fluorodopa PET (FD-PET) scan, and one autopsy was available. RESULTS: Eleven patients had a diagnosis of Parkinson's disease (PD), nine women. Mean age of onset was 52 with tremor-predominant dopa-responsive parkinsonism. Disease progression was slow but severe motor fluctuations occurred. One patient required subthalamic nucleus deep-brain stimulation with a good motor outcome. One patient had mental retardation, schizophrenia and became demented, and another patient was demented. Three patients and also two unaffected subjects had mild learning difficulties. All genetic tests yielded negative results. FD-PET showed marked asymmetric striatal tracer uptake deficiency, consistent with PD. Pathological examination demonstrated no Lewy bodies and immunostaining was negative for alpha-synuclein. CONCLUSION: Apart from a younger age of onset and a female predominance, the phenotype was indistinguishable from sporadic tremor-predominant PD, including FD-PET scan results. As known genetic causes of autosomal dominant PD were excluded, this family harbors a novel genetic defect.

The profile of general practitioners (GPs) who publish in selected family practice journals

Background: Providing support for research is one of the key issues in the ongoing attempts to improve Primary Care. However, when patient care takes up a significant part of a GP's time, conducting research is difficult. In this study we examine the working conditions and profile of GPs who publish in three leading medical journals and propose possible remedial policy actions. Findings: The authors of all articles published in 2006 and 2007 in three international Family Medicine journals - Annals of Family Medicine, Family Practice, and Journal of Family Practice - were contacted by E-mail. They were asked to complete a questionnaire investigating the following variables: availability of specific time for research, time devoted to research, number of patients attended, and university affiliation. Only GPs were included in the study. Three hundred and ten relevant articles published between 2006 and 2007 were identified and the authors contacted using a survey tool. 124 researchers responded to our questionnaire; 45% of respondents who were not GPs were excluded. On average GPs spent 2.52 days per week and 6.9 hours per day on patient care, seeing 45 patients per week. Seventy-five per cent of GPs had specific time assigned to research, on average 13 hours per week; 79% were affiliated to a university and 69% held teaching positions. Conclusions: Most GPs who publish original articles in leading journals have time specifically assigned to research as part of their normal working schedule. They see a relatively small number of patients. Improving the working conditions of family physicians who intend to investigate is likely to lead to better research results.

Is blood still thicker than water? A life-course perspective on the transformation of family and friends' roles in personal networks

This thesis addresses the issue of the moving boundaries between family and friends' roles in personal networks, adopting a life-course perspective and using Switzerland as a case study. In a period of major changes in personal life happening in contemporary Western societies, understanding the organization of personal networks intertwined with the unfolding of individual life courses is of prime importance in facing new challenges with regard to social integration. The data stem from a representative national survey carried out in 2011 named Family tiMes, including 803 individuals born either in 1950-1955 or in 1970-1975. An innovative research design was adopted, combing cross-sectional ego-centered network data and retrospective longitudinal life-course data. The results show continuing boundaries between family and friends' roles and that family keeps a prominent role in personal networks despite the notable importance of friendship ties. One relationship stands out above all, that with the partner, followed quite a few steps behind by those with children. Regarding life courses, de-standardization tendencies were found in family formation and also a persistent gendering of occupational trajectories. Two kinds of life trajectories are particularly intertwined with personal networks, co-residence and partnership trajectories, both related to the unfolding of family life. In particular, transition to parenthood functions as a turning point in individuals' lives, deeply transforming their sociability. Finally, a twofold pluralization process was identified, affecting simultaneously the organization of personal networks and the unfolding of individual life courses. This thesis contributes to the literature on the sociology of family and personal life, and to fruitful interlinkage between the network approach and the life-course perspective.

The AMPK family member Snf1 protects Saccharomyces cerevisiae cells upon glutatione oxidation

The AMPK/Snf1 kinase has a central role in carbon metabolism homeostasis in Saccharomyces cerevisiae. In this study, we show that Snf1 activity, which requires phosphorylation of the Thr210 residue, is needed for protection against selenite toxicity. Such protection involves the Elm1 kinase, which acts upstream of Snf1 to activate it. Basal Snf1 activity is sufficient for the defense against selenite, although Snf1 Thr210 phosphorylation levels become increased at advanced treatment times, probably by inhibition of the Snf1 dephosphorylation function of the Reg1 phosphatase. Contrary to glucose deprivation, Snf1 remains cytosolic during selenite treatment, and the protective function of the kinase does not require its known nuclear effectors. Upon selenite treatment, a null snf1 mutant displays higher levels of oxidized versus reduced glutathione compared to wild type cells, and its hypersensitivity to the agent is rescued by overexpression of the glutathione reductase gene GLR1. In the presence of agents such as diethyl maleate or diamide, which cause alterations in glutathione redox homeostasis by increasing the levels of oxidized glutathione, yeast cells also require Snf1 in an Elm1-dependent manner for growth. These observations demonstrate a role of Snf1 to protect yeast cells in situations where glutathione-dependent redox homeostasis is altered to a more oxidant intracellular environment and associates AMPK to responses against oxidative stress.

Implementing evidence-based patient and family education on oral anticoagulation therapy: a community-based participatory project.

AIMS AND OBJECTIVES: This study aimed at developing and implementing evidence-based patient and family education on oral anticoagulation therapy. BACKGROUND: The number of persons with chronic diseases who live at home is increasing. They have to manage multiple diseases and complex treatments. One such treatment is oral anticoagulation therapy, a high risk variable dose medication. Adherence to oral anticoagulation therapy is jeopardised by limited information about the medications, their risk and complications, the impact of individual daily routine and the limited inclusion of family members in education. Hence, improved and tailored education is essential for patients and families to manage oral anticoagulation therapy at home. DESIGN AND METHODS: A community-based participatory research design combined with the Precede-Proceed model was used including a systematic literature review, posteducation analysis, an online nurse survey, a documentation analysis and patient/family interviews. The study was conducted between April 2010-December 2012 at a department of general internal medicine in a teaching hospital in Switzerland. Participants were the department's nursing and medical professionals including the patients and their families. RESULTS: The evidence-based patient and family education on oral anticoagulation therapy emerged comprising a learning assessment, teaching units, clarification of responsibilities of nurse professionals and documentation guidelines. CONCLUSION AND CLINICAL RELEVANCE: The inclusion of the whole department has contributed to the development and implementation of this evidence-based patient family education on oral anticoagulation therapy, which encompasses local characteristics and patient preferences. This education is now being used throughout the department.

Tie Strength and Family Formation : which Personal Relationships are Influential ?

Based on the analysis of qualitative interviews in western Germany, we argue that personal relationships have a strong impact on individuals' family formation processes and fertility intentions. We identify persons who influence individuals' fertility choices. Strong ties, such as those among core family members (i.e., parents and siblings), are an important contributing factor, but we are also able to show that weak ties, such as those among colleagues, and acquaintances need to be considered when examining social influence on family formation processes. Apart from single network partners, we have identified influential groups of persons. Such groups serve as a comparative standard regarding the timing of having one's first child and subsequent children.

OmpA family proteins and Pmp-like autotransporter: new adhesins of Waddlia chondrophila.

Waddlia chondrophila is a obligate intracellular bacterium belonging to the Chlamydiales order, a clade that also includes the well-known classical Chlamydia responsible for a number of severe human and animal diseases. Waddlia is an emerging pathogen associated with adverse pregnancy outcomes in humans and abortion in ruminants. Adhesion to the host cell is an essential prerequisite for survival of every strict intracellular bacteria and, in classical Chlamydia, this step is partially mediated by polymorphic outer membrane proteins (Pmps), a family of highly diverse autotransporters that represent about 15% of the bacterial coding capacity. Waddlia chondrophila genome however only encodes one putative Pmp-like protein. Using a proteomic approach, we identified several bacterial proteins potentially implicated in the adhesion process and we characterized their expression during the replication cycle of the bacteria. In addition, we demonstrated that the Waddlia Pmp-like autotransporter as well as OmpA2 and OmpA3, two members of the extended Waddlia OmpA protein family, exhibit adhesive properties on epithelial cells. We hypothesize that the large diversity of the OmpA protein family is linked to the wide host range of these bacteria that are able to enter and multiply in various host cells ranging from protozoa to mammalian and fish cells.

Bcl-2 family proteins and cytoskeleton changes involved in DM-1 cytotoxic effect on melanoma cells

Melanoma is one of the most aggressive types of skin cancer and its incidence rate is still increasing. All existing treatments are minimally effective. Consequently, new therapeutic agents for melanoma treatment should be developed. The DM-1 compound is a curcumin analog that possesses several curcumin characteristics, such as antiproliferative, antitumor, and anti-metastatic properties. The aim of this study was to evaluate the different signaling pathways involved in the cytotoxic effect of DM-1 on melanoma cells. The apoptotic process and cytoskeletal changes were evaluated by immunoblotting and immunofluorescence, respectively, in melanoma cells. After DM-1 treatment, SK-MEL-5 melanoma cells showed actin filament disorganization with spicule formation throughout the cytoskeleton and significant reduction of focal adhesion as well as they were present only at cell extremities, conferring a poor connection between the cell and the substrate. Besides this, there was significant filopodium retraction and loss of typical cytoskeleton scaffold. These modifications contributed to cell detachment followed by cell death. Furthermore, DM-1-induced apoptosis was triggered by multiple Bcl-2 proteins involved in both the extrinsic and the intrinsic apoptotic pathways. SK-MEL-5 cells showed a death mechanism mainly by Bcl-2/Bax ratio decrease, whereas A375 cells presented apoptosis induction by Mcl-1 and Bcl-xL downregulation. In SK-MEL-5 and A375 melanoma cells, there was a significant increase in the active form of caspase 9, and the inactive form of the effector caspase 3 was decreased in both cell lines. Expression of cleaved poly ADP ribose polymerase was increased after DM-1 treatment in these melanoma cell lines, demonstrating that the apoptotic process occurred. Altogether, these data elucidate the cellular and molecular mechanisms involved in the cytotoxicity induced by the antitumor agent DM-1 in melanoma cells.

Compound Heterozygous VSX2 Mutation Causing Bilateral Anophthalmia in a Consanguineous Egyptian Family

Purpose: To report the clinical and genetic study of a child with bilateral anophthalmia. Methods: A 14-year-old Egyptian boy, born from consanguineous parents, underwent a general and a full ophthalmological examination. Mutation screen of the A/M genes with recessive inheritance was done stepwise and DNA was analyzed by Sanger sequencing. Results: Bilateral anophthalmia, arachnodactyly of the feet and high arched palate were observed on general examination. The parents were first cousins and healthy. Sequencing analysis revealed a novel compound heterozygous mutation in one of the copy of exon 2 of VSX2 and a possible deletion of at least exon 2 on the other allele. Conclusions: A compound heterozygous VSX2 mutation associated with anophthalmia was identified in a patient from an Egyptian consanguineous family. This report brings the number of VSX2 mutation in anophthalmia/microphthalmia (A/M) to 13. Functional consequences of the reported changes still need to be characterized, as well as the percentage of A/M caused by mutations in the VSX2 gene. This family also shows that despite consanguinity, heterozygous mutations can also happen and one should not restrict the molecular analysis to homozygous mutations.

High-density linkage maps fail to detect any genetic component to sex determination in a Rana temporaria family.

Sex chromosome differentiation in Rana temporaria varies strikingly among populations or families: whereas some males display well-differentiated Y haplotypes at microsatellite markers on linkage group 2 (LG2 ), others are genetically undistinguishable from females. We analysed with RADseq markers one family from a Swiss lowland population with no differentiated sex chromosomes, and where sibship analyses had failed to detect any association between the phenotypic sex of progeny and parental haplotypes. Offspring were reared in a common tank in outdoor conditions and sexed at the froglet stage. We could map a total of 2177 SNPs (1123 in the mother, 1054 in the father), recovering in both adults 13 linkage groups (= chromosome pairs) that were strongly syntenic to Xenopus tropicalis despite > 200 My divergence. Sexes differed strikingly in the localization of crossovers, which were uniformly distributed in the female but limited to chromosome ends in the male. None of the 2177 markers showed significant association with offspring sex. Considering the very high power of our analysis, we conclude that sex determination was not genetic in this family; which factors determined sex remain to be investigated.

CYP17A1 Enzyme activity is linked to ambulatory blood pressure in a family-based population study.

BACKGROUND: Genome-wide association studies have linked CYP17A1 coding for the steroid hormone synthesizing enzyme 17α-hydroxylase (CYP17A1) to blood pressure (BP). We hypothesized that the genetic signal may translate into a correlation of ambulatory BP (ABP) with apparent CYP17A1 activity in a family-based population study and estimated the heritability of CYP17A1 activity. METHODS: In the Swiss Kidney Project on Genes in Hypertension, day and night urinary excretions of steroid hormone metabolites were measured in 518 participants (220 men, 298 women), randomly selected from the general population. CYP17A1 activity was assessed by 2 ratios of urinary steroid metabolites: one estimating the combined 17α-hydroxylase/17,20-lyase activity (ratio 1) and the other predominantly 17α-hydroxylase activity (ratio 2). A mixed linear model was used to investigate the association of ABP with log-transformed CYP17A1 activities exploring effect modification by urinary sodium excretion. RESULTS: Daytime ABP was positively associated with ratio 1 under conditions of high, but not low urinary sodium excretion (P interaction <0.05). Ratio 2 was not associated with ABP. Heritability estimates (SE) for day and night CYP17A1 activities were 0.39 (0.10) and 0.40 (0.09) for ratio 1, and 0.71 (0.09) and 0.55 (0.09) for ratio 2 (P values <0.001). CYP17A1 activities, assessed with ratio 1, were lower in older participants. CONCLUSIONS: Low apparent CYP17A1 activity (assessed with ratio 1) is associated with elevated daytime ABP when salt intake is high. CYP17A1 activity is heritable and diminished in the elderly. These observations highlight the modifying effect of salt intake on the association of CYP17A1 with BP.

Radiographic adenoid evaluation: proposal of an objective parameter

Objective The objective of the present study was to evaluate current radiographic parameters designed to investigate adenoid hypertrophy and nasopharyngeal obstruction, and to present an alternative radiographic assessment method. Materials and Methods In order to do so, children (4 to14 years old) who presented with nasal obstruction or oral breathing complaints were submitted to cavum radiographic examination. One hundred and twenty records were evaluated according to quantitative radiographic parameters, and data were correlated with a gold-standard videonasopharyngoscopic study, in relation to the percentage of choanal obstruction. Subsequently, a regression analysis was performed in order to create an original model so the percentage of the choanal obstruction could be predicted. Results The quantitative parameters demonstrated moderate, if not weak correlation with the real percentage of choanal obstruction. The regression model (110.119*A/N) demonstrated a satisfactory ability to “predict” the actual percentage of choanal obstruction. Conclusion Since current adenoid quantitative radiographic parameters present limitations, the model presented by the present study might be considered as an alternative assessment method in cases where videonasopharyngoscopic evaluation is unavailable.

A nonparametric method for the measurement of size diversity with emphasis on data standardization

The most suitable method for estimation of size diversity is investigated. Size diversity is computed on the basis of the Shannon diversity expression adapted for continuous variables, such as size. It takes the form of an integral involving the probability density function (pdf) of the size of the individuals. Different approaches for the estimation of pdf are compared: parametric methods, assuming that data come from a determinate family of pdfs, and nonparametric methods, where pdf is estimated using some kind of local evaluation. Exponential, generalized Pareto, normal, and log-normal distributions have been used to generate simulated samples using estimated parameters from real samples. Nonparametric methods include discrete computation of data histograms based on size intervals and continuous kernel estimation of pdf. Kernel approach gives accurate estimation of size diversity, whilst parametric methods are only useful when the reference distribution have similar shape to the real one. Special attention is given for data standardization. The division of data by the sample geometric mean is proposedas the most suitable standardization method, which shows additional advantages: the same size diversity value is obtained when using original size or log-transformed data, and size measurements with different dimensionality (longitudes, areas, volumes or biomasses) may be immediately compared with the simple addition of ln k where kis the dimensionality (1, 2, or 3, respectively). Thus, the kernel estimation, after data standardization by division of sample geometric mean, arises as the most reliable and generalizable method of size diversity evaluation