Understanding and characterising the diffusive transport of capillary oxygen and nutrients in striated muscles is key to assessing angiogenesis and investigating the efficacy of experimental and therapeutic interventions for numerous pathological conditions, such as chronic ischaemia. In articular, the influence of both muscle tissue and microvascular heterogeneities on capillary oxygen supply is poorly understood. The objective of this thesis is to develop mathematical and computational modelling frameworks for the purpose of extending and generalising the current use of histology in estimating the regions of tissue supplied by individual capillaries to facilitate the exploration of functional capillary oxygen supply in striated muscles. In particular, we aim to investigate the balance between local capillary supply of oxygen and oxygen demand in the presence of various anatomical and functional heterogeneities, by capturing tissue details from histological imaging and estimating or predicting regions of capillary supply. Our computational method throughout is based on a finite element framework that captures the anatomical details of tissue cross sections. In Chapter 1 we introduce the problem. In Chapter 2 we develop a theoretical model to describe oxygen transport from capillaries to uniform muscle tissues (e.g. cardiac muscle). Transport is then explored in terms of oxygen levels and capillary supply regions. In Chapter 3 we extend this modelling framework to explore the influence of the surrounding tissue by accounting for the spatial anisotropies of fibre oxygen demand and diffusivity and the heterogeneity in fibre size and shape, as exemplified by mixed muscle tissues (e.g. skeletal muscle). We additionally explore the effects of diffusion through the interstitium, facilitated--diffusion by myoglobin, and Michaelis--Menten kinetics of tissue oxygen consumption. In Chapter 4, a further extension is pursued to account for intracellular heterogeneities in mitochondrial distribution and diffusive parameters. As a demonstration of the potential of the models derived in Chapters 2--4, in Chapter 5 we simulate oxygen transport in myocardial tissue biopsies from rats with either impaired angiogenesis or impaired arteriolar perfusion. Quantitative predictions are made to help explain and support experimental measurements of cardiac performance and metabolism. In the final chapter we summarize the main results and indicate directions for further work.

Anxiety and depression are highly prevalent diseases with common heritable elements, but the particular genetic mechanisms and biological pathways underlying them are poorly understood. Part of the challenge in understanding the genetic basis of these disorders is that they are polygenic and often context-dependent. In my thesis, I apply a series of modern statistical tools to ascertain some of the myriad genetic and environmental factors that underlie fear-related behaviours in nearly two thousand heterogeneous stock mice, which serve as animal models of anxiety and depression. Using a Bayesian method called Sparse Partitioning and a frequentist method called Bagphenotype, I identify gene-by-sex interactions that contribute to variation in fear-related behaviours, such as those displayed in the elevated plus maze and the open field test, although I demonstrate that the contributions are generally small. Also using Bagphenotype, I identify hundreds of gene-by-environment interactions related to these traits. The interacting environmental covariates are diverse, ranging from experimenter to season of the year. With gene expression data from a brain structure associated with anxiety called the hippocampus, I generate modules of co-expressed genes and map them to the genome. Two of these modules were enriched for key nervous system components — one for dendritic spines, another for oligodendrocyte markers — but I was unable to find significant correlations between them and fear-related behaviours. Finally, I employed another Bayesian technique, Sparse Instrumental Variables, which takes advantage of conditional probabilities to identify hippocampus genes whose expression appears not just to be associated with variation in fear-related behaviours, but cause variation in those phenotypes.

We investigate the biological constraints determined by the complex relationships between ecological and immunological processes of host-pathogen interactions, with emphasis on influenza viruses in human, which are responsible for a number of pandemics in the last 150 years. We begin by discussing prolegomenous reviews of historical perspectives on the use of theoretical modelling as a complementary tool in public health and epidemiology, current biological background motivating the objective of the thesis, and derivations of mathematical models of multi-locus-allele systems for infectious diseases with co-circulating serotypes. We provide detailed analysis of the multi-locus-allele model and its age-specific extension. In particular, we establish the necessary conditions for the local asymptotic stability of the steady states and the existence of oscillatory behaviours. For the age-structured model, results on the existence of a mild solution and stability conditions are presented. Numerical studies of various strain spaces show that the dynamic features are preserved. Specifically, we demonstrate that discrete antigenic forms of pathogens can exhibit three distinct dynamic features, where antigenic variants (i) fully self-organize and co-exist with no strain structure (NSS), (ii) sort themselves into discrete strain structure (DSS) with non-overlapping or minimally overlapping clusters under the principle of competitive exclusion, or (iii) exhibit cyclical strain structure (CSS) where dominant antigenic types are cyclically replaced with sharp epidemics dominated by (1) a single strain dominance with irregular emergence and re-emergence of certain pathogenic forms, (2) ordered alternating appearance of a single antigenic type in periodic or quasi-periodic form similar to periodic travelling waves, (3) erratic appearance and disappearance of synchrony between discrete antigenic types, and (4) phase-synchronization with uncorrelated amplitudes. These analyses allow us to gain insight into the age-specific immunological profile in order to untangle the effects of strain structures as captured by the clustering behaviours, and to provide public health implications. The age-structured model can be used to investigate the effect of age-specific targeting for public health purposes.

This biomedical engineering thesis explores the opportunities and challenges of 2D+t contrast echocardiography for left ventricle functional analysis, both clinically and within a computer vision atlas-based deformable template model framework. A database was created for the experiments in this thesis, with 21 studies of contrast Dobutamine Stress Echo, in all 4 principal planes. The database includes clinical variables, human expert hand-traced myocardial contours and visual scoring. First the problem is studied from a clinical perspective. Quantification of endocardial global and local function using standard measures shows expected values and agreement with human expert visual scoring, but the results are less reliable for myocardial thickening. Next, the problem of segmenting the endocardium with a computer is posed in a standard landmark and atlas-based deformable template model framework. The underlying assumption is that these models can emulate human experts in terms of integrating previous knowledge about the anatomy and physiology with three sources of information from the image: texture, geometry and kinetics. Probabilistic atlases of contrast echocardiography are computed, while noting from histograms at selected anatomical locations that modelling texture with just mean intensity values may be too naive. Intensity analysis together with the clinical results above suggest that lack of external boundary definition may preclude this imaging technique for appropriate measuring of myocardial thickening, while endocardial boundary definition is appropriate for evaluation of wall motion. Geometry is presented in a Principal Component Analysis (PCA) context, highlighting issues about Gaussianity, the correlation and covariance matrices with respect to physiology, and analysing different measures of dimensionality. A popular extension of deformable models ---Active Appearance Models (AAMs)--- is then studied in depth. Contrary to common wisdom, it is contended that using a PCA texture space instead of a fixed atlas is detrimental to segmentation, and that PCA models are not convenient for texture modelling. To integrate kinetics, a novel spatio-temporal model of cardiac contours is proposed. The new explicit model does not require frame interpolation, and it is compared to previous implicit models in terms of approximation error when the shape vector changes from frame to frame or remains constant throughout the cardiac cycle. Finally, the 2D+t atlas-based deformable model segmentation problem is formulated and solved with a gradient descent approach. Experiments using the similarity transformation suggest that segmentation of the whole cardiac volume outperforms segmentation of individual frames. A relatively new approach ---the inverse compositional algorithm--- is shown to decrease running times of the classic Lucas-Kanade algorithm by a factor of 20 to 25, to values that are within real-time processing reach.