Risk factors associated with corneal nerve alteration in type 1 diabetes in the absence of neuropathy: A longitudinal in vivo corneal confocal microscopy study

Purpose The aim of this study was to determine alterations to the corneal subbasal nerve plexus (SNP) over four years using in vivo corneal confocal microscopy (IVCM) in participants with type 1 diabetes and to identify significant risk factors associated with these alterations. Methods A cohort of 108 individuals with type 1 diabetes and no evidence of peripheral neuropathy at enrollment underwent laser-scanning IVCM, ocular screening, and health and metabolic assessment at baseline and the examinations continued for four subsequent annual visits. At each annual visit, eight central corneal images of the SNP were selected and analyzed to quantify corneal nerve fiber density (CNFD), branch density (CNBD) and fiber length (CNFL). Linear mixed model approaches were fitted to examine the relationship between risk factors and corneal nerve parameters. Results A total of 96 participants completed the final visit and 91 participants completed all visits. No significant relationships were found between corneal nerve parameters and time, sex, duration of diabetes, smoking, alcohol consumption, blood pressure or BMI. However, CNFD was negatively associated with HbA1c (β=-0.76, P<0.01) and age (β=-0.13, P<0.01) and positively related to high density lipids (HDL) (β=2.01, P=0.03). Higher HbA1c (β=-1.58, P=0.04) and age (β=-0.23, P<0.01) also negatively impacted CNBD. CNFL was only affected by higher age (β=-0.06, P<0.01). Conclusions Glycemic control, HDL and age have significant effects on SNP structure. These findings highlight the importance of diabetic management to prevent corneal nerve damage as well as the capability of IVCM for monitoring subclinical alterations in the corneal SNP in diabetes.

The inferior whorl for detecting diabetic peripheral neuropathy using corneal confocal microscopy

PURPOSE: In vivo corneal confocal microscopy (CCM) is increasingly used as a surrogate endpoint in studies of diabetic polyneuropathy (DPN). However, it is not clear whether imaging the central cornea provides optimal diagnostic utility for DPN. Therefore, we compared nerve morphology in the central cornea and the inferior whorl, a more distal and densely innervated area located inferior and nasal to the central cornea. METHODS: A total of 53 subjects with type 1/type 2 diabetes and 15 age-matched control subjects underwent detailed assessment of neuropathic symptoms (NPS), deficits (neuropathy disability score [NDS]), quantitative sensory testing (vibration perception threshold [VPT], cold and warm threshold [CT/WT], and cold- and heat-induced pain [CIP/HIP]), and electrophysiology (sural and peroneal nerve conduction velocity [SSNCV/PMNCV], and sural and peroneal nerve amplitude [SSNA/PMNA]) to diagnose patients with (DPN+) and without (DPN-) neuropathy. Corneal nerve fiber density (CNFD) and length (CNFL) in the central cornea, and inferior whorl length (IWL) were quantified. RESULTS: Comparing control subjects to DPN- and DPN+ patients, there was a significant increase in NDS (0 vs. 2.6 ± 2.3 vs. 3.3 ± 2.7, P < 0.01), VPT (V; 5.4 ± 3.0 vs. 10.6 ± 10.3 vs. 17.7 ± 11.8, P < 0.01), WT (°C; 37.7 ± 3.5 vs. 39.1 ± 5.1 vs. 41.7 ± 4.7, P < 0.05), and a significant decrease in SSNCV (m/s; 50.2 ± 5.4 vs. 48.4 ± 5.0 vs. 39.5 ± 10.6, P < 0.05), CNFD (fibers/mm2; 37.8 ± 4.9 vs. 29.7 ± 7.7 vs. 27.1 ± 9.9, P < 0.01), CNFL (mm/mm2; 27.5 ± 3.6 vs. 24.4 ± 7.8 vs. 20.7 ± 7.1, P < 0.01), and IWL (mm/mm2; 35.1 ± 6.5 vs. 26.2 ± 10.5 vs. 23.6 ± 11.4, P < 0.05). For the diagnosis of DPN, CNFD, CNFL, and IWL achieved an area under the curve (AUC) of 0.75, 0.74, and 0.70, respectively, and a combination of IWL-CNFD achieved an AUC of 0.76. CONCLUSIONS: The parameters of CNFD, CNFL, and IWL have a comparable ability to diagnose patients with DPN. However, IWL detects an abnormality even in patients without DPN. Combining IWL with CNFD may improve the diagnostic performance of CCM.

Bayesian estimation of the dynamics of pandemic (H1N1) 2009 influenza transmission in Queensland: A space–time SIR-based model

Background A pandemic strain of influenza A spread rapidly around the world in 2009, now referred to as pandemic (H1N1) 2009. This study aimed to examine the spatiotemporal variation in the transmission rate of pandemic (H1N1) 2009 associated with changes in local socio-environmental conditions from May 7–December 31, 2009, at a postal area level in Queensland, Australia. Method We used the data on laboratory-confirmed H1N1 cases to examine the spatiotemporal dynamics of transmission using a flexible Bayesian, space–time, Susceptible-Infected-Recovered (SIR) modelling approach. The model incorporated parameters describing spatiotemporal variation in H1N1 infection and local socio-environmental factors. Results The weekly transmission rate of pandemic (H1N1) 2009 was negatively associated with the weekly area-mean maximum temperature at a lag of 1 week (LMXT) (posterior mean: −0.341; 95% credible interval (CI): −0.370–−0.311) and the socio-economic index for area (SEIFA) (posterior mean: −0.003; 95% CI: −0.004–−0.001), and was positively associated with the product of LMXT and the weekly area-mean vapour pressure at a lag of 1 week (LVAP) (posterior mean: 0.008; 95% CI: 0.007–0.009). There was substantial spatiotemporal variation in transmission rate of pandemic (H1N1) 2009 across Queensland over the epidemic period. High random effects of estimated transmission rates were apparent in remote areas and some postal areas with higher proportion of indigenous populations and smaller overall populations. Conclusions Local SEIFA and local atmospheric conditions were associated with the transmission rate of pandemic (H1N1) 2009. The more populated regions displayed consistent and synchronized epidemics with low average transmission rates. The less populated regions had high average transmission rates with more variations during the H1N1 epidemic period.

Lost and found: recovery of the holotype of the ocellated angelshark, Squatina tergocellatoides Chen, 1963 (Squatinidae), with comments on western Pacific squatinids

The ocellated angelshark, Squatina tergocellatoides, Chen, 1963 is redescribed from the holotype, which was thought to be lost. Its recent recovery has allowed for a revised description, including new data, and comparison to other Western Pacific squatinids. Squatina tergocellatoides can be distinguished from its congeners by three pairs of prominent large black spots, each with a diameter greater than eye length; two on each pectoral fin at anterior and posterior angles and one on each side near the tail base; another three pairs of lesser defined spots, one large spot on base of each dorsal fin and one located laterally on each side of tail located below first dorsal fin. Ventral surface is uniformly white to cream coloured, and margins of pectoral fins and tail similar in colour to dorsal side. Pectoral fins with angular lateral apices and rounded posterior lobe, pelvic fin tips not reaching origin of first dorsal fin, strongly fringed nasal barbels, small inter-orbital space, head and mouth lengths, broad internarial width and pelvic fin base, a very small pelvic girdle width, and a caudal fin with triangular ventral lobe greater in length than dorsal lobe. Comments on additional specimens are provided, as well as observations on biogeography. A review of western Pacific squatinids is also provided.

Assessing the combined effect of asbestos exposure and smoking on lung cancer: a Bayesian approach

We review the literature on the combined association between lung cancer and two environmental exposures, asbestos exposure and smoking, and explore a Bayesian approach to assess evidence of interaction between the exposures. The meta-analysis combines separate indices of additive and multiplicative relationships and multivariate relative risk estimates. By making inferences on posterior probabilities we can explore both the form and strength of interaction. This analysis may be more informative than providing evidence to support one relation over another on the basis of statistical significance. Overall, we find evidence for a more than additive and less than multiplicative relation.

Morphological changes in the conjunctiva, episclera and sclera following short-term miniscleral contact lens wear in rigid lens neophytes

PURPOSE To quantify the influence of short-term wear of miniscleral contact lenses on the morphology of the corneo-scleral limbus, the conjunctiva, episclera and sclera. METHODS OCT images of the anterior eye were captured before, immediately following 3h of wear and then 3h after removal of a miniscleral contact lens for 10 young (27±5 years) healthy participants (neophyte rigid lens wearers). The region of analysis encompassed 1mm anterior, to 3.5mm posterior to the scleral spur. Natural diurnal variations in thickness were measured on a separate day and compensated for in subsequent analyses. RESULTS Following 3h of lens wear, statistically significant tissue thinning was observed across all quadrants, with a mean decrease in thickness of -24.1±3.6μm (p<0.001), which diminished, but did not return to baseline 3h after lens removal (-16.9±1.9μm, p<0.001). The largest tissue compression was observed in the superior quadrant (-49.9±8.5μm, p<0.01) and in the annular zone 1.5mm from the scleral spur (-48.2±5.7μm), corresponding to the approximate edge of the lens landing zone. Compression of the conjunctiva/episclera accounted for about 70% of the changes. CONCLUSIONS Optimal fitting miniscleral contact lenses worn for three hours resulted in significant tissue compression in young healthy eyes, with the greatest thinning observed superiorly, potentially due to the additional force of the eyelid, with a partial recovery of compression 3h after lens removal. Most of the morphological changes occur in the conjunctiva/episclera layers.

The distribution and abundance of electrosensory pores in two benthic sharks: a comparison of the wobbegong shark, Orectolobus maculatus, and the angel shark, Squatina australis

Electroreception is an ancient sense found in many aquatic animals, including sharks, which may be used in the detection of prey, predators and mates. Wobbegong sharks (Orectolobidae) and angel sharks (Squatinidae) represent two distantly related families that have independently evolved a similar dorso-ventrally compressed body form to complement their benthic ambush feeding strategy. Consequently, these groups represent useful models in which to investigate the specific morphological and physiological adaptations that are driven by the adoption of a benthic lifestyle. In this study, we compared the distribution and abundance of electrosensory pores in the spotted wobbegong shark (Orectolobus maculatus) with the Australian angel shark (Squatina australis) to determine whether both species display a similar pattern of clustering of sub-dermal electroreceptors and to further understand the functional importance of electroreception in the feeding behaviour of these benthic sharks. Orectolobus maculatus has a more complex electrosensory system than S. australis, with a higher abundance of pores and an additional cluster of electroreceptors positioned in the snout (the superficial ophthalmic cluster). Interestingly, both species possess a cluster of pores (the hyoid cluster, positioned slightly posterior to the first gill slit) more commonly found in rays, but which may be present in all benthic elasmobranchs to assist in the detection of approaching predators.

On probability-based inference under data missing by design

Whether a statistician wants to complement a probability model for observed data with a prior distribution and carry out fully probabilistic inference, or base the inference only on the likelihood function, may be a fundamental question in theory, but in practice it may well be of less importance if the likelihood contains much more information than the prior. Maximum likelihood inference can be justified as a Gaussian approximation at the posterior mode, using flat priors. However, in situations where parametric assumptions in standard statistical models would be too rigid, more flexible model formulation, combined with fully probabilistic inference, can be achieved using hierarchical Bayesian parametrization. This work includes five articles, all of which apply probability modeling under various problems involving incomplete observation. Three of the papers apply maximum likelihood estimation and two of them hierarchical Bayesian modeling. Because maximum likelihood may be presented as a special case of Bayesian inference, but not the other way round, in the introductory part of this work we present a framework for probability-based inference using only Bayesian concepts. We also re-derive some results presented in the original articles using the toolbox equipped herein, to show that they are also justifiable under this more general framework. Here the assumption of exchangeability and de Finetti's representation theorem are applied repeatedly for justifying the use of standard parametric probability models with conditionally independent likelihood contributions. It is argued that this same reasoning can be applied also under sampling from a finite population. The main emphasis here is in probability-based inference under incomplete observation due to study design. This is illustrated using a generic two-phase cohort sampling design as an example. The alternative approaches presented for analysis of such a design are full likelihood, which utilizes all observed information, and conditional likelihood, which is restricted to a completely observed set, conditioning on the rule that generated that set. Conditional likelihood inference is also applied for a joint analysis of prevalence and incidence data, a situation subject to both left censoring and left truncation. Other topics covered are model uncertainty and causal inference using posterior predictive distributions. We formulate a non-parametric monotonic regression model for one or more covariates and a Bayesian estimation procedure, and apply the model in the context of optimal sequential treatment regimes, demonstrating that inference based on posterior predictive distributions is feasible also in this case.

Bayesian fisheries stock assessment: integrating and updating knowledge

In this thesis the use of the Bayesian approach to statistical inference in fisheries stock assessment is studied. The work was conducted in collaboration of the Finnish Game and Fisheries Research Institute by using the problem of monitoring and prediction of the juvenile salmon population in the River Tornionjoki as an example application. The River Tornionjoki is the largest salmon river flowing into the Baltic Sea. This thesis tackles the issues of model formulation and model checking as well as computational problems related to Bayesian modelling in the context of fisheries stock assessment. Each article of the thesis provides a novel method either for extracting information from data obtained via a particular type of sampling system or for integrating the information about the fish stock from multiple sources in terms of a population dynamics model. Mark-recapture and removal sampling schemes and a random catch sampling method are covered for the estimation of the population size. In addition, a method for estimating the stock composition of a salmon catch based on DNA samples is also presented. For most of the articles, Markov chain Monte Carlo (MCMC) simulation has been used as a tool to approximate the posterior distribution. Problems arising from the sampling method are also briefly discussed and potential solutions for these problems are proposed. Special emphasis in the discussion is given to the philosophical foundation of the Bayesian approach in the context of fisheries stock assessment. It is argued that the role of subjective prior knowledge needed in practically all parts of a Bayesian model should be recognized and consequently fully utilised in the process of model formulation.

Mathematical modelling of human corneal oxygenation and cell kinetics in corneal epithelial wound healing

Healthy transparent cornea depends upon the regulation of fluid, nutrient and oxygen transport through the tissue to sustain cell metabolism and other critical processes for normal functioning. This research considers the corneal geometry and investigates oxygen distribution using a two-dimensional Monod kinetic model, showing that previous studies make assumptions that lead to predictions of near-anoxic levels of oxygen tension in the limbal regions of the cornea. It also considers the comparison of experimental spatial and temporal data with the predictions of novel mathematical models with respect to distributed mitotic rates during corneal epithelial wound healing.

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.

Evidence from molecules and morphology expands Podonomopsis Brundin (Diptera: Chironomidae: Podonominae) to include ‘genus Chile'

The informal taxon ‘genus Chile’ of Brundin, based solely on pupal exuviae of a podonomine Chironomidae, has remained inadequately known for half a century. New collections reveal life associations, and provide molecular data to hypothesise a precise phylogenetic placement in the austral Podonominae. A densely sampled molecular phylogeny based on two nuclear and one mitochondrial DNA markers shows ‘genus Chile’ to be the sister group to Podonomopsis Brundin, 1966. Within Podonomopsis a clade of South American species is sister to all Australian species. We discuss how to rank such a sister group taxon and treat ‘genus Chile’ as a new subgenus Araucanopsis, subg. nov. with the new species, Podonomopsis (Araucanopsis) avelasse, sp. nov. from Chile and Argentina as genotype of the monotypic subgenus. We describe P. (A.) avelasse in all stages and provide an expanded diagnosis and description of Podonomopsis to include Araucanopsis. A dated biogeographic hypothesis (chronogram) infers the most recent common ancestor (tmcra) of expanded Podonomopsis at 95 million years ago (Mya) (68–122 Mya 95% highest posterior density), ‘core’ Podonomopsis at 83 Mya (58–108) and Australian Podonomopsis at 65 Mya (44–87). All dates are before the South America–Australia geological separation through Antarctica, supporting previous conclusions that the taxon distribution is ‘Gondwanan’ in origin. Podonomopsis, even as expanded here, remains unknown from New Zealand or elsewhere on extant Zealandia.

Computational Techniques for Haplotype Inference and for Local Alignment Significance

This thesis which consists of an introduction and four peer-reviewed original publications studies the problems of haplotype inference (haplotyping) and local alignment significance. The problems studied here belong to the broad area of bioinformatics and computational biology. The presented solutions are computationally fast and accurate, which makes them practical in high-throughput sequence data analysis. Haplotype inference is a computational problem where the goal is to estimate haplotypes from a sample of genotypes as accurately as possible. This problem is important as the direct measurement of haplotypes is difficult, whereas the genotypes are easier to quantify. Haplotypes are the key-players when studying for example the genetic causes of diseases. In this thesis, three methods are presented for the haplotype inference problem referred to as HaploParser, HIT, and BACH. HaploParser is based on a combinatorial mosaic model and hierarchical parsing that together mimic recombinations and point-mutations in a biologically plausible way. In this mosaic model, the current population is assumed to be evolved from a small founder population. Thus, the haplotypes of the current population are recombinations of the (implicit) founder haplotypes with some point--mutations. HIT (Haplotype Inference Technique) uses a hidden Markov model for haplotypes and efficient algorithms are presented to learn this model from genotype data. The model structure of HIT is analogous to the mosaic model of HaploParser with founder haplotypes. Therefore, it can be seen as a probabilistic model of recombinations and point-mutations. BACH (Bayesian Context-based Haplotyping) utilizes a context tree weighting algorithm to efficiently sum over all variable-length Markov chains to evaluate the posterior probability of a haplotype configuration. Algorithms are presented that find haplotype configurations with high posterior probability. BACH is the most accurate method presented in this thesis and has comparable performance to the best available software for haplotype inference. Local alignment significance is a computational problem where one is interested in whether the local similarities in two sequences are due to the fact that the sequences are related or just by chance. Similarity of sequences is measured by their best local alignment score and from that, a p-value is computed. This p-value is the probability of picking two sequences from the null model that have as good or better best local alignment score. Local alignment significance is used routinely for example in homology searches. In this thesis, a general framework is sketched that allows one to compute a tight upper bound for the p-value of a local pairwise alignment score. Unlike the previous methods, the presented framework is not affeced by so-called edge-effects and can handle gaps (deletions and insertions) without troublesome sampling and curve fitting.

The neuroscience of inspirational leadership: the importance of collective-oriented language and shared group membership

Effective leaders are believed to inspire followers by providing inclusive visions of the future that followers can identify with. In the present study, we examined the neural mechanisms underlying this process, testing key hypotheses derived from transformational and social identity approaches to leadership. While undergoing functional MRI, supporters from the two major Australian political parties (Liberal vs. Labor) were presented with inspirational collective-oriented and noninspirational personal-oriented statements made by in-group and out-group leaders. Imaging data revealed that inspirational (rather than noninspirational) statements from in-group leaders were associated with increased activation in the bilateral rostral inferior parietal lobule, pars opercularis, and posterior midcingulate cortex: brain areas that are typically implicated in controlling semantic information processing. In contrast, for out-group leaders, greater activation in these areas was associated with noninspirational statements. In addition, noninspirational statements by in-group (but not out-group) leaders resulted in increased activation in the medial prefrontal cortex, an area typically associated with reasoning about a person’s mental state. These results show that followers processed identical statements qualitatively differently as a function of leaders’ group membership, thus demonstrating that shared identity acts as an amplifier for inspirational leadership communication.

Tissue thickness calculation in ocular optical coherence tomography

Thickness measurements derived from optical coherence tomography (OCT) images of the eye are a fundamental clinical and research metric, since they provide valuable information regarding the eye’s anatomical and physiological characteristics, and can assist in the diagnosis and monitoring of numerous ocular conditions. Despite the importance of these measurements, limited attention has been given to the methods used to estimate thickness in OCT images of the eye. Most current studies employing OCT use an axial thickness metric, but there is evidence that axial thickness measures may be biased by tilt and curvature of the image. In this paper, standard axial thickness calculations are compared with a variety of alternative metrics for estimating tissue thickness. These methods were tested on a data set of wide-field chorio-retinal OCT scans (field of view (FOV) 60° x 25°) to examine their performance across a wide region of interest and to demonstrate the potential effect of curvature of the posterior segment of the eye on the thickness estimates. Similarly, the effect of image tilt was systematically examined with the same range of proposed metrics. The results demonstrate that image tilt and curvature of the posterior segment can affect axial tissue thickness calculations, while alternative metrics, which are not biased by these effects, should be considered. This study demonstrates the need to consider alternative methods to calculate tissue thickness in order to avoid measurement error due to image tilt and curvature.