Bayesian Inference for Retrospective Population Genetics Models Using Markov Chain Monte Carlo Methods

Genetics, the science of heredity and variation in living organisms, has a central role in medicine, in breeding crops and livestock, and in studying fundamental topics of biological sciences such as evolution and cell functioning. Currently the field of genetics is under a rapid development because of the recent advances in technologies by which molecular data can be obtained from living organisms. In order that most information from such data can be extracted, the analyses need to be carried out using statistical models that are tailored to take account of the particular genetic processes. In this thesis we formulate and analyze Bayesian models for genetic marker data of contemporary individuals. The major focus is on the modeling of the unobserved recent ancestry of the sampled individuals (say, for tens of generations or so), which is carried out by using explicit probabilistic reconstructions of the pedigree structures accompanied by the gene flows at the marker loci. For such a recent history, the recombination process is the major genetic force that shapes the genomes of the individuals, and it is included in the model by assuming that the recombination fractions between the adjacent markers are known. The posterior distribution of the unobserved history of the individuals is studied conditionally on the observed marker data by using a Markov chain Monte Carlo algorithm (MCMC). The example analyses consider estimation of the population structure, relatedness structure (both at the level of whole genomes as well as at each marker separately), and haplotype configurations. For situations where the pedigree structure is partially known, an algorithm to create an initial state for the MCMC algorithm is given. Furthermore, the thesis includes an extension of the model for the recent genetic history to situations where also a quantitative phenotype has been measured from the contemporary individuals. In that case the goal is to identify positions on the genome that affect the observed phenotypic values. This task is carried out within the Bayesian framework, where the number and the relative effects of the quantitative trait loci are treated as random variables whose posterior distribution is studied conditionally on the observed genetic and phenotypic data. In addition, the thesis contains an extension of a widely-used haplotyping method, the PHASE algorithm, to settings where genetic material from several individuals has been pooled together, and the allele frequencies of each pool are determined in a single genotyping.

AndroMedia - Towards a Context-aware Mobile Music Recommender

Portable music players have made it possible to listen to a personal collection of music in almost every situation, and they are often used during some activity to provide a stimulating audio environment. Studies have demonstrated the effects of music on the human body and mind, indicating that selecting music according to situation can, besides making the situation more enjoyable, also make humans perform better. For example, music can boost performance during physical exercises, alleviate stress and positively affect learning. We believe that people intuitively select different types of music for different situations. Based on this hypothesis, we propose a portable music player, AndroMedia, designed to provide personalised music recommendations using the user’s current context and listening habits together with other user’s situational listening patterns. We have developed a prototype that consists of a central server and a PDA client. The client uses Bluetooth sensors to acquire context information and logs user interaction to infer implicit user feedback. The user interface also allows the user to give explicit feedback. Large user interface elements facilitate touch-based usage in busy environments. The prototype provides the necessary framework for using the collected information together with other user’s listening history in a context- enhanced collaborative filtering algorithm to generate context-sensitive recommendations. The current implementation is limited to using traditional collaborative filtering algorithms. We outline the techniques required to create context-aware recommendations and present a survey on mobile context-aware music recommenders found in literature. As opposed to the explored systems, AndroMedia utilises other users’ listening habits when suggesting tunes, and does not require any laborious set up processes.

Bringing nature to schools to promote children’s physical activity

Physical activity (PA) is essential for human health and wellbeing across all age, socioeconomic and ethnic groups. Engagement with the natural world is a new defining criterion for enhancing the benefits of PA particularly for children and young people. Interacting with nature benefits children’s social and emotional wellbeing, develops resilience and reduces the risk of obesity and type 2 diabetes across all population groups. Governments around the world are now recognising the importance of children spending more active time outdoors. However, children’s outdoor activities, free play and nature-related exploration are often structured and supervised by adults due to safety concerns and risks. In this context schools become more accessible and safe options for children to engage in PA outdoors with the presence of nature features. Research on school designs involving young children has revealed that children prefer nature-related features in school environments. Affordances in nature may increase children’s interest in physically active behaviours. Given that present school campuses are designed for operational efficiency and economic reasons there is a need to re-design schools responding to the positive role of nature on human health. If schools were re-designed to incorporate diverse natural features children’s PA and consequent health and wellbeing would likely improve markedly.

Doctors’ perspectives on law and life-sustaining treatment: Survey design and recruitment strategies for a challenging cohort

- Background Palliative medicine and other specialists play significant legal roles in decisions to withhold and withdraw life-sustaining treatment at the end of life. Yet little is known about their knowledge of or attitudes to the law, and the role they think it should play in medical practice. Consideration of doctors’ views is critical to optimizing patient outcomes at the end of life. However, doctors are difficult to engage as participants in empirical research, presenting challenges for researchers seeking to understand doctors’ experiences and perspectives. - Aims To determine how to engage doctors involved in end-of-life care in empirical research about knowledge of the law and the role it plays in medical practice at the end of life. - Methods Postal survey of all specialists in palliative medicine, emergency medicine, geriatric medicine, intensive care, medical oncology, renal medicine, and respiratory medicine in three Australian states: New South Wales, Victoria, and Queensland. The survey was sent in hard copy with two reminders and a follow up reminder letter was also sent to the directors of hospital emergency departments. Awareness was further promoted through engagement with the relevant medical colleges and publications in professional journals; various incentives to respond were also used. The key measure is the response rate of doctors to the survey. - Results Thirty-two percent of doctors in the main study completed their survey with response rate by specialty ranging from 52% (palliative care) to 24% (medical oncology). This overall response rate was twice that of the reweighted pilot study (16%). - Conclusions Doctors remain a difficult cohort to engage in survey research but strategic recruitment efforts can be effective in increasing response rate. Collaboration with doctors and their professional bodies in both the development of the survey instrument and recruitment of participants is essential.

Characterisation, cloning and production of industrially interesting enzymes: gluconolactone oxidase of Penicillium cyaneo-fulvum and gluconate 5-dehydrogenase of Gluconobacter suboxydans

The work covered in this thesis is focused on the development of technology for bioconversion of glucose into D-erythorbic acid (D-EA) and 5-ketogluconic acid (5-KGA). The task was to show on proof-of-concept level the functionality of the enzymatic conversion or one-step bioconversion of glucose to these acids. The feasibility of both studies to be further developed for production processes was also evaluated. The glucose - D-EA bioconversion study was based on the use of a cloned gene encoding a D-EA forming soluble flavoprotein, D-gluconolactone oxidase (GLO). GLO was purified from Penicillium cyaneo-fulvum and partially sequenced. The peptide sequences obtained were used to isolate a cDNA clone encoding the enzyme. The cloned gene (GenBank accession no. AY576053) is homologous to the other known eukaryotic lactone oxidases and also to some putative prokaryotic lactone oxidases. Analysis of the deduced protein sequence of GLO indicated the presence of a typical secretion signal sequence at the N-terminus of the enzyme. No other targeting/anchoring signals were found, suggesting that GLO is the first known lactone oxidase that is secreted rather than targeted to the membranes of the endoplasmic reticulum or mitochondria. Experimental evidence supports this analysis, as near complete secretion of GLO was observed in two different yeast expression systems. Highest expression levels of GLO were obtained using Pichia pastoris as an expression host. Recombinant GLO was characterised and the suitability of purified GLO for the production of D-EA was studied. Immobilised GLO was found to be rapidly inactivated during D-EA production. The feasibility of in vivo glucose - D-EA conversion using a P. pastoris strain co-expressing the genes of GLO and glucose oxidase (GOD, E.C. 1.1.3.4) of A. niger was demonstrated. The glucose - 5-KGA bioconversion study followed a similar strategy to that used in the D-EA production research. The rationale was based on the use of a cloned gene encoding a membrane-bound pyrroloquinoline quinone (PQQ)-dependent gluconate 5-dehydrogenase (GA 5-DH). GA 5-DH was purified to homogeneity from the only source of this enzyme known in literature, Gluconobacter suboxydans, and partially sequenced. Using the amino acid sequence information, the GA 5-DH gene was cloned from a genomic library of G. suboxydans. The cloned gene was sequenced (GenBank accession no. AJ577472) and found to be an operon of two adjacent genes encoding two subunits of GA 5-DH. It turned out that GA 5-DH is a rather close homologue of a sorbitol dehydrogenase from another G. suboxydans strain. It was also found that GA 5-DH has significant polyol dehydrogenase activity. The G. suboxydans GA 5-DH gene was poorly expressed in E. coli. Under optimised conditions maximum expression levels of GA 5-DH did not exceed the levels found in wild-type G. suboxydans. Attempts to increase expression levels resulted in repression of growth and extensive cell lysis. However, the expression levels were sufficient to demonstrate the possibility of bioconversion of glucose and gluconate into 5-KGA using recombinant strains of E. coli. An uncharacterised homologue of GA 5-DH was identified in Xanthomonas campestris using in silico screening. This enzyme encoded by chromosomal locus NP_636946 was found by a sequencing project of X. campestris and named as a hypothetical glucose dehydrogenase. The gene encoding this uncharacterised enzyme was cloned, expressed in E. coli and found to encode a gluconate/polyol dehydrogenase without glucose dehydrogenase activity. Moreover, the X. campestris GA 5-DH gene was expressed in E. coli at nearly 30 times higher levels than the G. suboxydans GA 5-DH gene. Good expressability of the X. campestris GA-5DH gene makes it a valuable tool not only for 5-KGA production in the tartaric acid (TA) bioprocess, but possibly also for other bioprocesses (e.g. oxidation of sorbitol into L-sorbose). In addition to glucose - 5-KGA bioconversion, a preliminary study of the feasibility of enzymatic conversion of 5-KGA into TA was carried out. Here, the efficacy of the first step of a prospective two-step conversion route including a transketolase and a dehydrogenase was confirmed. It was found that transketolase convert 5-KGA into TA semialdehyde. A candidate for the second step was suggested to be succinic dehydrogenase, but this was not tested. The analysis of the two subprojects indicated that bioconversion of glucose to TA using X. campestris GA 5-DH should be prioritised first and the process development efforts in future should be focused on development of more efficient GA 5-DH production strains by screening a more suitable production host and by protein engineering.

A constraint Jacobian based approach for static analysis of pantograph masts

This paper presents a constraint Jacobian matrix based approach to obtain the stiffness matrix of widely used deployable pantograph masts with scissor-like elements (SLE). The stiffness matrix is obtained in symbolic form and the results obtained agree with those obtained with the force and displacement methods available in literature. Additional advantages of this approach are that the mobility of a mast can be evaluated, redundant links and joints in the mast can be identified and practical masts with revolute joints can be analysed. Simulations for a hexagonal mast and an assembly with four hexagonal masts is presented as illustrations.

Modelling heterogeneity of concrete using 2D lattice network for concrete fracture and comparison with AE study

In this paper, numerical modelling of fracture in concrete using two-dimensional lattice model is presented and also a few issues related to lattice modelling technique applicable to concrete fracture are reviewed. A comparison is made with acoustic emission (AE) events with the number of fractured elements. To implement the heterogeneity of the plain concrete, two methods namely, by generating grain structure of the concrete using Fuller's distribution and the concrete material properties are randomly distributed following Gaussian distribution are used. In the first method, the modelling of the concrete at meso level is carried out following the existing methods available in literature. The shape of the aggregates present in the concrete are assumed as perfect spheres and shape of the same in two-dimensional lattice network is circular. A three-point bend (TPB) specimen is tested in the experiment under crack mouth opening displacement (CMOD) control at a rate of 0.0004 mm/sec and the fracture process in the same TPB specimen is modelled using regular triangular 2D lattice network. Load versus crack mouth opening isplacement (CMOD) plots thus obtained by using both the methods are compared with experimental results. It was observed that the number of fractured elements increases near the peak load and beyond the peak load. That is once the crack starts to propagate. AE hits also increase rapidly beyond the peak load. It is compulsory here to mention that although the lattice modelling of concrete fracture used in this present study is very similar to those already available in literature, the present work brings out certain finer details which are not available explicitly in the earlier works.

Reconstruction of the ocular surface using biomaterial templates

This chapter discusses the effect on vision of a large group of pathological conditions, known as ocular surface disorders (OSDs), and presents the therapeutic strategies to reconstruct the abnormal ocular surface. If left untreated, most of the OSDs will lead to partial or total loss of eyesight, especially when limbal stem cell deficiency is involved. An overview of various treatment strategies is presented, with the emphasis on the development of the ex vivo expansion of corneal limbal epithelial cells (presumed to be progenitor or stem cells) and the creation of transplantable epithelial constructs. The use of naturally derived biomaterials (collagen, fibrin, amnion, etc.) or synthetic polymers (polylactides, thermoresponsive polymers, etc.) as substrata in these constructs is critically analyzed. Emphasis is placed on the templates from silk proteins, which are being developed by the authors.

Biomaterial templates for the culture and transplantation of retinal pigment epithelial cells: A critical review

The idea of retinal cell transplantation as a potential treatment for age-related retinal degeneration, a leading cause of blindness in the Western world, has been around for a number of decades. To date, however, it has not been entirely successful; one of the main reasons for this is the lack of an ideal substratum for the retinal cells, specifically for the growth of retinal pigment epithelial cells prior to transplantation. This chapter reviews the reasoning behind this potential treatment, the development of animal transplantation models for human trials, the prerequisites of an ideal substratum, the past and current research on substratum materials, and the potential for future developments in this area.

Characterization of anomalous relaxation using the time-fractional Bloch equation and multiple echo T2 *-weighted magnetic resonance imaging at 7T

PURPOSE To study the utility of fractional calculus in modeling gradient-recalled echo MRI signal decay in the normal human brain. METHODS We solved analytically the extended time-fractional Bloch equations resulting in five model parameters, namely, the amplitude, relaxation rate, order of the time-fractional derivative, frequency shift, and constant offset. Voxel-level temporal fitting of the MRI signal was performed using the classical monoexponential model, a previously developed anomalous relaxation model, and using our extended time-fractional relaxation model. Nine brain regions segmented from multiple echo gradient-recalled echo 7 Tesla MRI data acquired from five participants were then used to investigate the characteristics of the extended time-fractional model parameters. RESULTS We found that the extended time-fractional model is able to fit the experimental data with smaller mean squared error than the classical monoexponential relaxation model and the anomalous relaxation model, which do not account for frequency shift. CONCLUSIONS We were able to fit multiple echo time MRI data with high accuracy using the developed model. Parameters of the model likely capture information on microstructural and susceptibility-induced changes in the human brain.

Testing and Modeling of MR Damper and Its Application to SDOF Systems Using Integral Backstepping Technique

Magnetorheological dampers are intrinsically nonlinear devices, which make the modeling and design of a suitable control algorithm an interesting and challenging task. To evaluate the potential of magnetorheological (MR) dampers in control applications and to take full advantages of its unique features, a mathematical model to accurately reproduce its dynamic behavior has to be developed and then a proper control strategy has to be taken that is implementable and can fully utilize their capabilities as a semi-active control device. The present paper focuses on both the aspects. First, the paper reports the testing of a magnetorheological damper with an universal testing machine, for a set of frequency, amplitude, and current. A modified Bouc-Wen model considering the amplitude and input current dependence of the damper parameters has been proposed. It has been shown that the damper response can be satisfactorily predicted with this model. Second, a backstepping based nonlinear current monitoring of magnetorheological dampers for semi-active control of structures under earthquakes has been developed. It provides a stable nonlinear magnetorheological damper current monitoring directly based on system feedback such that current change in magnetorheological damper is gradual. Unlike other MR damper control techniques available in literature, the main advantage of the proposed technique lies in its current input prediction directly based on system feedback and smooth update of input current. Furthermore, while developing the proposed semi-active algorithm, the dynamics of the supplied and commanded current to the damper has been considered. The efficiency of the proposed technique has been shown taking a base isolated three story building under a set of seismic excitation. Comparison with widely used clipped-optimal strategy has also been shown.

Impact of exercise selection on hamstring muscle activation

Objective: To determine the extent to which different strength training exercises selectively activate the commonly injured biceps femoris long head (BFLH) muscle. Methods: This two-part observational study recruited 24 recreationally active males. Part 1 explored the amplitudes and the ratios of lateral to medial hamstring (BF/MH) normalised electromyography (nEMG) during the concentric and eccentric phases of 10 common strength training exercises. Part 2 used functional magnetic resonance imaging (fMRI) to determine the spatial patterns of hamstring activation during two exercises which i) most selectively, and ii) least selectively activated the BF in part 1. Results: Eccentrically, the largest BF/MH nEMG ratio was observed in the 45° hip extension exercise and the lowest was observed in the Nordic hamstring (NHE) and bent-knee bridge exercises. Concentrically, the highest BF/MH nEMG ratio was observed during the lunge and 45° hip extension and the lowest was observed for the leg curl and bent-knee bridge. fMRI revealed a greater BFLH to semitendinosus activation ratio in the 45° hip extension than the NHE (p<0.001). The T2 increase after hip extension for BFLH, semitendinosus and semimembranosus muscles were greater than that for BFSH (p<0.001). During the NHE, the T2 increase was greater for the semitendinosus than for the other hamstrings (p≤0.002). Conclusion: This investigation highlights the non-uniformity of hamstring activation patterns in different tasks and suggests that hip extension exercise more selectively activates the BFLH while the NHE preferentially recruits the semitendinosus. These findings have implications for strength training interventions aimed at preventing hamstring injury.

Fidgeting is associated with lower mortality risk

Context Evidence from prospective cohort studies has suggested that high volumes of reported daily sitting time is associated with mortality.1 ,2 However, not all have observed the same association.3 Fidgeting (small movements associated with nervousness or impatience), could provide additional energy expenditure when sitting, although the relationship with sitting and health outcomes had yet to be examined. Hagger-Johnson et al examined data from nearly 13 000 women to determine whether fidgeting modified the association between sitting time and mortality. Methods This study featured prospective data from 12 778 participants (aged 37–78 years) in the Women's Cohort Study (UK). Average daily sitting time was reported for weekdays and weekend days, and combined …

The grumpy stump - ultrasound and stump appendicitis

BACKGROUND: Stump appendicitis, the inflammation of an incompletely removed appendix, is a rare clinical presentation. Sonography can be useful in the diagnosis of the condition; by either directly visualising the inflamed stump, or by identifying signs of peri-caecal inflammation that can raise suspicion of the condition. SUMMARY OF WORK: A potential case of stump appendicitis was identified. This prompted a review of literature focused on the incidence of stump appendicitis, utility of ultrasound to identify an inflamed appendiceal stump, and surgical techniques used in appendectomy. SUMMARY OF RESULTS: Stump appendicitis is rare, with as few as 61 cases identified in literature during the last 60 years1. Of the two common techniques of appendectomy, which are ligation and invagination, the former can sometimes leave a residual stump that acts as a potential lumen for the pathophysiological process of appendicitis to recur. Established sonographic criteria for appendicitis also apply to the residual stump. Sonographic secondary signs that suggest the presence of acute appendicitis2 have also been demonstrated in cases of stump appendicitis3, even in the absence of an identifiable stump. DISCUSSION AND CONCLUSIONS: Appendicitis is usually dismissed in patients with a history of appendectomy. Though uncommon, sonographers should be aware of stump appendicitis in post-appendectomy patients that present with right iliac fossa pain from months to decades later. REFERENCE(S) 1. Subramanian A, Liang MK. A 60-year literature review of stump appendicitis: The need for a critical view. Am J Surg; 2012;203(4):503–7. 2. Reddan T, Corness J, Mengersen K, Harden F. Ultrasound of paediatric appendicitis and its secondary sonographic signs: providing a more meaningful finding. J Med Radiat Sci. 2015;DOI: 10.1002/jmrs.154 3. Martínez Chamorro E, Merina Castilla A, Muñoz Fraile B, Koren Fernández L, Borruel Nacenta S. Stump appendicitis: Preoperative imaging findings in four cases. Abdom Imaging. 2013;38(6):1214–9.

Chromatographic separation, fractionation and oxidation of selected beta-lactoglobulin peptides

The literature review elucidates the mechanism of oxidation in proteins and amino acids and gives an overview of the detection and analysis of protein oxidation products as well as information about ?-lactoglobulin and studies carried out on modifications of this protein under certain conditions. The experimental research included the fractionation of the tryptic peptides of ?-lactoglobulin using preparative-HPLC-MS and monitoring the oxidation process of these peptides via reverse phase-HPLC-UV. Peptides chosen to be oxidized were selected with respect to their amino acid content which were susceptible to oxidation and fractionated according to their m/z values. These peptides were: IPAVFK (m/z 674), ALPMHIR (m/z 838), LIVTQTMK (m/z 934) and VLVLDTDYK (m/z 1066). Even though it was not possible to solely isolate the target peptides due to co-elution of various fractions, the percentages of target peptides in the samples were satisfactory to carry out the oxidation procedure. IPAVFK and VLVLDTDYK fractions were found to yield the oxidation products reviewed in literature, however, unoxidized peptides were still present in high amounts after 21 days of oxidation. The UV data at 260 and 280 nm enabled to monitor both the main peptides and the oxidation products due to the absorbance of aromatic side-chains these peptides possess. ALPMHIR and LIVTQTMK fractions were oxidatively consumed rapidly and oxidation products of these peptides were observed even on day 0. High rates of depletion of these peptides were acredited to the presence of His (H) and sulfur-containing side-chains of Met (M). In conclusion, selected peptides hold the potential to be utilized as marker peptides in ?-lactoglobulin oxidation.