Mapping the regional influence of genetics on brain structure variability - A Tensor-Based Morphometry study

Genetic and environmental factors influence brain structure and function profoundly. The search for heritable anatomical features and their influencing genes would be accelerated with detailed 3D maps showing the degree to which brain morphometry is genetically determined. As part of an MRI study that will scan 1150 twins, we applied Tensor-Based Morphometry to compute morphometric differences in 23 pairs of identical twins and 23 pairs of same-sex fraternal twins (mean age: 23.8 ± 1.8 SD years). All 92 twins' 3D brain MRI scans were nonlinearly registered to a common space using a Riemannian fluid-based warping approach to compute volumetric differences across subjects. A multi-template method was used to improve volume quantification. Vector fields driving each subject's anatomy onto the common template were analyzed to create maps of local volumetric excesses and deficits relative to the standard template. Using a new structural equation modeling method, we computed the voxelwise proportion of variance in volumes attributable to additive (A) or dominant (D) genetic factors versus shared environmental (C) or unique environmental factors (E). The method was also applied to various anatomical regions of interest (ROIs). As hypothesized, the overall volumes of the brain, basal ganglia, thalamus, and each lobe were under strong genetic control; local white matter volumes were mostly controlled by common environment. After adjusting for individual differences in overall brain scale, genetic influences were still relatively high in the corpus callosum and in early-maturing brain regions such as the occipital lobes, while environmental influences were greater in frontal brain regions that have a more protracted maturational time-course.

Incorporating process variability into stormwater quality modelling

Process variability in pollutant build-up and wash-off generates inherent uncertainty that affects the outcomes of stormwater quality models. Poor characterisation of process variability constrains the accurate accounting of the uncertainty associated with pollutant processes. This acts as a significant limitation to effective decision making in relation to stormwater pollution mitigation. The study undertaken developed three theoretical scenarios based on research findings that variations in particle size fractions <150µm and >150µm during pollutant build-up and wash-off primarily determine the variability associated with these processes. These scenarios, which combine pollutant build-up and wash-off processes that takes place on a continuous timeline, are able to explain process variability under different field conditions. Given the variability characteristics of a specific build-up or wash-off event, the theoretical scenarios help to infer the variability characteristics of the associated pollutant process that follows. Mathematical formulation of the theoretical scenarios enables the incorporation of variability characteristics of pollutant build-up and wash-off processes in stormwater quality models. The research study outcomes will contribute to the quantitative assessment of uncertainty as an integral part of the interpretation of stormwater quality modelling outcomes.

Rainfall variability drives interannual variation in N2O emissions from a humid, subtropical pasture

Variations in interannual rainfall totals can lead to large uncertainties in annual N2O emission budget estimates from short term field studies. The interannual variation in nitrous oxide (N2O) emissions from a subtropical pasture in Queensland, Australia, was examined using continuous measurements of automated chambers over 2 consecutive years. Nitrous oxide emissions were highest during the summer months and were highly episodic, related more to the size and distribution of rain events than soil water content. Over 48% of the total N2O emitted was lost in just 16% of measurement days. Interannual variation in annual N2O estimates was high, with cumulative emissions increasing with decreasing rainfall. Cumulative emissions averaged 1826.7 ± 199.9 g N2O-N ha−1 yr−1 over the two year period, though emissions from 2008 (2148 ± 273 g N2O-N ha−1 yr−1) were 42% higher than 2007 (1504 ± 126 g N2O-N ha−1 yr−1). This increase in annual emissions coincided with almost half of the summer precipitation from 2007 to 2008. Emissions dynamics were chiefly driven by the distribution and size of rain events which varied on a seasonal and annual basis. Sampling frequency effects on cumulative N2O flux estimation were assessed using a jackknife technique to inform future manual sampling campaigns. Test subsets of the daily measured data were generated for the pasture and two adjacent land-uses (rainforest and lychee orchard) by selecting measured flux values at regular time intervals ranging from 1 to 30 days. Errors associated with weekly sampling were up to 34% of the sub-daily mean and were highly biased towards overestimation if strategically sampled following rain events. Sampling time of day also played a critical role. Morning sampling best represented the 24 hour mean in the pasture, whereas sampling at noon proved the most accurate in the shaded rainforest and lychee orchard.

Quantification of stability in an agility drill using linear and nonlinear measures of variability

This study implemented linear and nonlinear methods of measuring variability to determine differences in stability of two groups of skilled (n = 10) and unskilled (n = 10) participants performing 3m forward/backward shuttle agility drill. We also determined whether stability measures differed between the forward and backward segments of the drill. Finally, we sought to investigate whether local dynamic stability, measured using largest finite-time Lyapunov exponents, changed from distal to proximal lower extremity segments. Three-dimensional coordinates of five lower extremity markers data were recorded. Results revealed that the Lyapunov exponents were lower (P < 0.05) for skilled participants at all joint markers indicative of higher levels of local dynamic stability. Additionally, stability of motion did not differ between forward and backward segments of the drill (P > 0.05), signifying that almost the same control strategy was used in forward and backward directions by all participants, regardless of skill level. Furthermore, local dynamic stability increased from distal to proximal joints (P < 0.05) indicating that stability of proximal segments are prioritized by the neuromuscular control system. Finally, skilled participants displayed greater foot placement standard deviation values (P < 0.05), indicative of adaptation to task constraints. The results of this study provide new methods for sport scientists, coaches to characterize stability in agility drill performance.

Increasing functional variability in the preparatory phase of the takeoff improves elite springboard diving performance

PURPOSE: Previous research demonstrating that specific performance outcome goals can be achieved in different ways is functionally significant for springboard divers whose performance environment can vary extensively. This body of work raises questions about the traditional approach of balking (terminating the takeoff) by elite divers aiming to perform only identical, invariant movement patterns during practice. METHOD: A 12-week training program (2 times per day; 6.5 hr per day) was implemented with 4 elite female springboard divers to encourage them to adapt movement patterns under variable takeoff conditions and complete intended dives, rather than balk. RESULTS: Intraindividual analyses revealed small increases in variability in the board-work component of each diver's pretraining and posttraining program reverse-dive takeoffs. No topological differences were observed between movement patterns of dives completed pretraining and posttraining. Differences were noted in the amount of movement variability under different training conditions (evidenced by higher normalized root mean square error indexes posttraining). An increase in the number of completed dives (from 78.91%-86.84% to 95.59%-99.29%) and a decrease in the frequency of balked takeoffs (from 13.16%-19.41% to 0.63%-4.41%) showed that the elite athletes were able to adapt their behaviors during the training program. These findings coincided with greater consistency in the divers' performance during practice as scored by qualified judges. CONCLUSION: Results suggested that on completion of training, athletes were capable of successfully adapting their movement patterns under more varied takeoff conditions to achieve greater consistency and stability of performance outcomes.

Population of computational rabbit-specific ventricular action potential models for investigating sources of variability in cellular repolarisation

Variability is observed at all levels of cardiac electrophysiology. Yet, the underlying causes and importance of this variability are generally unknown, and difficult to investigate with current experimental techniques. The aim of the present study was to generate populations of computational ventricular action potential models that reproduce experimentally observed intercellular variability of repolarisation (represented by action potential duration) and to identify its potential causes. A systematic exploration of the effects of simultaneously varying the magnitude of six transmembrane current conductances (transient outward, rapid and slow delayed rectifier K(+), inward rectifying K(+), L-type Ca(2+), and Na(+)/K(+) pump currents) in two rabbit-specific ventricular action potential models (Shannon et al. and Mahajan et al.) at multiple cycle lengths (400, 600, 1,000 ms) was performed. This was accomplished with distributed computing software specialised for multi-dimensional parameter sweeps and grid execution. An initial population of 15,625 parameter sets was generated for both models at each cycle length. Action potential durations of these populations were compared to experimentally derived ranges for rabbit ventricular myocytes. 1,352 parameter sets for the Shannon model and 779 parameter sets for the Mahajan model yielded action potential duration within the experimental range, demonstrating that a wide array of ionic conductance values can be used to simulate a physiological rabbit ventricular action potential. Furthermore, by using clutter-based dimension reordering, a technique that allows visualisation of multi-dimensional spaces in two dimensions, the interaction of current conductances and their relative importance to the ventricular action potential at different cycle lengths were revealed. Overall, this work represents an important step towards a better understanding of the role that variability in current conductances may play in experimentally observed intercellular variability of rabbit ventricular action potential repolarisation.

Accounting for individual variability in the von Bertalanffy growth model

Estimation of von Bertalanffy growth parameters has received considerable attention in fisheries research. Since Sainsbury (1980, Can. J. Fish. Aquat. Sci. 37: 241-247) much of this research effort has centered on accounting for individual variability in the growth parameters. In this paper we demonstrate that, in analysis of tagging data, Sainsbury's method and its derivatives do not, in general, satisfactorily account for individual variability in growth, leading to inconsistent parameter estimates (the bias does not tend to zero as sample size increases to infinity). The bias arises because these methods do not use appropriate conditional expectations as a basis for estimation. This bias is found to be similar to that of the Fabens method. Such methods would be appropriate only under the assumption that the individual growth parameters that generate the growth increment were independent of the growth parameters that generated the initial length. However, such an assumption would be unrealistic. The results are derived analytically, and illustrated with a simulation study. Until techniques that take full account of the appropriate conditioning have been developed, the effect of individual variability on growth has yet to be fully understood.

Maximum likelihood estimation of mortality and growth with individual variability from multiple length-frequency data

We consider estimation of mortality rates and growth parameters from length-frequency data of a fish stock and derive the underlying length distribution of the population and the catch when there is individual variability in the von Bertalanffy growth parameter L-infinity. The model is flexible enough to accommodate 1) any recruitment pattern as a function of both time and length, 2) length-specific selectivity, and 3) varying fishing effort over time. The maximum likelihood method gives consistent estimates, provided the underlying distribution for individual variation in growth is correctly specified. Simulation results indicate that our method is reasonably robust to violations in the assumptions. The method is applied to tiger prawn data (Penaeus semisulcatus) to obtain estimates of natural and fishing mortality.

Effect of individual variability on estimation of population parameters from length-frequency data

We consider estimation of mortality rates and growth parameters from length-frequency data of a fish stock when there is individual variability in the von Bertalanffy growth parameter L-infinity and investigate the possible bias in the estimates when the individual variability is ignored. Three methods are examined: (i) the regression method based on the Beverton and Holt's (1956, Rapp. P.V. Reun. Cons. Int. Explor. Mer, 140: 67-83) equation; (ii) the moment method of Powell (1979, Rapp. PV. Reun. Int. Explor. Mer, 175: 167-169); and (iii) a generalization of Powell's method that estimates the individual variability to be incorporated into the estimation. It is found that the biases in the estimates from the existing methods are, in general, substantial, even when individual variability in growth is small and recruitment is uniform, and the generalized method performs better in terms of bias but is subject to a larger variation. There is a need to develop robust and flexible methods to deal with individual variability in the analysis of length-frequency data.

Sampling Methods for Exploring Between Subject Variability in Cardiac Electrophysiology Experiments

Between-subject and within-subject variability is ubiquitous in biology and physiology and understanding and dealing with this is one of the biggest challenges in medicine. At the same time it is difficult to investigate this variability by experiments alone. A recent modelling and simulation approach, known as population of models (POM), allows this exploration to take place by building a mathematical model consisting of multiple parameter sets calibrated against experimental data. However, finding such sets within a high-dimensional parameter space of complex electrophysiological models is computationally challenging. By placing the POM approach within a statistical framework, we develop a novel and efficient algorithm based on sequential Monte Carlo (SMC). We compare the SMC approach with Latin hypercube sampling (LHS), a method commonly adopted in the literature for obtaining the POM, in terms of efficiency and output variability in the presence of a drug block through an in-depth investigation via the Beeler-Reuter cardiac electrophysiological model. We show improved efficiency via SMC and that it produces similar responses to LHS when making out-of-sample predictions in the presence of a simulated drug block.

FTO genotype is associated with phenotypic variability of body mass index

There is evidence across several species for genetic control of phenotypic variation of complex traits1, 2, 3, 4, such that the variance among phenotypes is genotype dependent. Understanding genetic control of variability is important in evolutionary biology, agricultural selection programmes and human medicine, yet for complex traits, no individual genetic variants associated with variance, as opposed to the mean, have been identified. Here we perform a meta-analysis of genome-wide association studies of phenotypic variation using ~170,000 samples on height and body mass index (BMI) in human populations. We report evidence that the single nucleotide polymorphism (SNP) rs7202116 at the FTO gene locus, which is known to be associated with obesity (as measured by mean BMI for each rs7202116 genotype)5, 6, 7, is also associated with phenotypic variability. We show that the results are not due to scale effects or other artefacts, and find no other experiment-wise significant evidence for effects on variability, either at loci other than FTO for BMI or at any locus for height. The difference in variance for BMI among individuals with opposite homozygous genotypes at the FTO locus is approximately 7%, corresponding to a difference of ~0.5 kilograms in the standard deviation of weight. Our results indicate that genetic variants can be discovered that are associated with variability, and that between-person variability in obesity can partly be explained by the genotype at the FTO locus. The results are consistent with reported FTO by environment interactions for BMI8, possibly mediated by DNA methylation9, 10. Our BMI results for other SNPs and our height results for all SNPs suggest that most genetic variants, including those that influence mean height or mean BMI, are not associated with phenotypic variance, or that their effects on variability are too small to detect even with samples sizes greater than 100,000.

Variability in exposure to ambient ultrafine particles in urban schools: Comparative assessment between Australia and Spain

Ambient ultrafine particle number concentrations (PNC) have inhomogeneous spatio-temporal distributions and depend on a number of different urban factors, including background conditions and distant sources. This paper quantitatively compares exposure to ambient ultrafine particles at urban schools in two cities in developed countries, with high insolation climatic conditions, namely Brisbane (Australia) and Barcelona (Spain). The analysis used comprehensive indoor and outdoor air quality measurements at 25 schools in Brisbane and 39 schools in Barcelona. PNC modes were analysed with respect to ambient temperature, land use and urban characteristics, combined with the measured elemental carbon concentrations, NOx (Brisbane) and NO2 (Barcelona). The trends and modes of the quantified weekday average daily cycles of ambient PNC exhibited significant differences between the two cities. PNC increases were observed during traffic rush hours in both cases. However, the mid-day peak was dominant in Brisbane schools and had the highest contribution to total PNC for both indoors and outdoors. In Barcelona, the contribution from traffic was highest for ambient PNC, while the mid-day peak had a slightly higher contribution for indoor concentrations. Analysis of the relationships between PNC and land use characteristics in Barcelona schools showed a moderate correlation with the percentage of road network area and an anti-correlation with the percentage of green area. No statistically significant correlations were found for Brisbane. Overall, despite many similarities between the two cities, school-based exposure patterns were different. The main source of ambient PNC at schools was shown to be traffic in Barcelona and mid-day new particle formation in Brisbane. The mid-day PNC peak in Brisbane could have been driven by the combined effect of background and meteorological conditions, as well as other local/distant sources. The results have implications for urban development, especially in terms of air quality mitigation and management at schools.

Moral expansiveness: Examining variability in the extension of the moral world

The nature of our moral judgments—and the extent to which we treat others with care—depend in part on the distinctions we make between entities deemed worthy or unworthy of moral consideration— our moral boundaries. Philosophers, historians, and social scientists have noted that people’s moral boundaries have expanded over the last few centuries, but the notion of moral expansiveness has received limited empirical attention in psychology. This research explores variations in the size of individuals’ moral boundaries using the psychological construct of moral expansiveness and introduces the Moral Expansiveness Scale (MES), designed to capture this variation. Across 6 studies, we established the reliability, convergent validity, and predictive validity of the MES. Moral expansiveness was related (but not reducible) to existing moral constructs (moral foundations, moral identity, “moral” universalism values), predictors of moral standing (moral patiency and warmth), and other constructs associated with concern for others (empathy, identification with humanity, connectedness to nature, and social responsibility). Importantly, the MES uniquely predicted willingness to engage in prosocial intentions and behaviors at personal cost independently of these established constructs. Specifically, the MES uniquely predicted willingness to prioritize humanitarian and environmental concerns over personal and national self-interest, willingness to sacrifice one’s life to save others (ranging from human out-groups to animals and plants), and volunteering behavior. Results demonstrate that moral expansiveness is a distinct and important factor in understanding moral judgments and their consequences.