Mass conservative three-dimensional water tracer distribution from Markov chain Monte Carlo inversion of time-lapse ground-penetrating radar data

Time-lapse geophysical measurements are widely used to monitor the movement of water and solutes through the subsurface. Yet commonly used deterministic least squares inversions typically suffer from relatively poor mass recovery, spread overestimation, and limited ability to appropriately estimate nonlinear model uncertainty. We describe herein a novel inversion methodology designed to reconstruct the three-dimensional distribution of a tracer anomaly from geophysical data and provide consistent uncertainty estimates using Markov chain Monte Carlo simulation. Posterior sampling is made tractable by using a lower-dimensional model space related both to the Legendre moments of the plume and to predefined morphological constraints. Benchmark results using cross-hole ground-penetrating radar travel times measurements during two synthetic water tracer application experiments involving increasingly complex plume geometries show that the proposed method not only conserves mass but also provides better estimates of plume morphology and posterior model uncertainty than deterministic inversion results.

Improved segmentation of deep brain grey matter structures using magnetization transfer (MT) parameter maps.

Basal ganglia and brain stem nuclei are involved in the pathophysiology of various neurological and neuropsychiatric disorders. Currently available structural T1-weighted (T1w) magnetic resonance images do not provide sufficient contrast for reliable automated segmentation of various subcortical grey matter structures. We use a novel, semi-quantitative magnetization transfer (MT) imaging protocol that overcomes limitations in T1w images, which are mainly due to their sensitivity to the high iron content in subcortical grey matter. We demonstrate improved automated segmentation of putamen, pallidum, pulvinar and substantia nigra using MT images. A comparison with segmentation of high-quality T1w images was performed in 49 healthy subjects. Our results show that MT maps are highly suitable for automated segmentation, and so for multi-subject morphometric studies with a focus on subcortical structures.

A new locus for congenital cataract, microcornea, microphthalmia, and atypical iris coloboma maps to chromosome 2.

OBJECTIVE: To report a novel phenotype of autosomal dominant atypical congenital cataract associated with variable expression of microcornea, microphthalmia, and iris coloboma linked to chromosome 2. Molecular analysis of this phenotype may improve our understanding of anterior segment development. DESIGN: Observational case study, genome linkage analysis, and gene mutation screening. PARTICIPANTS: Three families, 1 Egyptian and 2 Belgians, with a total of 31 affected were studied. METHODS: Twenty-one affected subjects and 9 first-degree relatives underwent complete ophthalmic examination. In the Egyptian family, exclusion of PAX6, CRYAA, and MAF genes was demonstrated by haplotype analysis using microsatellite markers on chromosomes 11, 16, and 21. Genome-wide linkage analysis was then performed using 385 microsatellite markers on this family. In the 2 Belgian families, the PAX6 gene was screened for mutations by direct sequencing of all exons. MAIN OUTCOME MEASURES: Phenotype description, genome-wide linkage of the phenotype, linkage to the PAX6, CRYAA, and MAF genes, and mutation detection in the PAX6 gene. RESULTS: Affected members of the 3 families had bilateral congenital cataracts inherited in an autosomal dominant pattern. A novel form of hexagonal nuclear cataract with cortical riders was expressed. Among affected subjects with available data, 95% had microcornea, 39% had microphthalmia, and 38% had iris coloboma. Seventy-five percent of the colobomata were atypical, showing a nasal superior location in 56%. A positive lod score of 4.86 was obtained at theta = 0 for D2S2309 on chromosome 2, a 4.9-Mb common haplotype flanked by D2S2309 and D2S2358 was obtained in the Egyptian family, and linkage to the PAX6, CRYAA, or MAF gene was excluded. In the 2 Belgian families, sequencing of the junctions and all coding exons of PAX6 did not reveal any molecular change. CONCLUSIONS: We describe a novel phenotype that includes the combination of a novel form of congenital hexagonal cataract, with variably expressed microcornea, microphthalmia, and atypical iris coloboma, not caused by PAX6 and mapping to chromosome 2. FINANCIAL DISCLOSURE(S): The authors have no proprietary or commercial interest in any materials discussed in this article.

Bayesian Markov-chain-Monte-Carlo inversion of time-lapse cross hole ground-penetrating radar data to characterize the vadose zone at the Arrenaes field site, Denmark

The ground-penetrating radar (GPR) geophysical method has the potential to provide valuable information on the hydraulic properties of the vadose zone because of its strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR traveltime data can allow for a significant reduction in uncertainty regarding subsurface van Genuchten-Mualem (VGM) parameters. Much of the previous work on the stochastic estimation of VGM parameters from crosshole GPR data has considered the case of steady-state infiltration conditions, which represent only a small fraction of practically relevant scenarios. We explored in detail the dynamic infiltration case, specifically examining to what extent time-lapse crosshole GPR traveltimes, measured during a forced infiltration experiment at the Arreneas field site in Denmark, could help to quantify VGM parameters and their uncertainties in a layered medium, as well as the corresponding soil hydraulic properties. We used a Bayesian Markov-chain-Monte-Carlo inversion approach. We first explored the advantages and limitations of this approach with regard to a realistic synthetic example before applying it to field measurements. In our analysis, we also considered different degrees of prior information. Our findings indicate that the stochastic inversion of the time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions compared with the corresponding priors, which in turn significantly improves knowledge of soil hydraulic properties. Overall, the results obtained clearly demonstrate the value of the information contained in time-lapse GPR data for characterizing vadose zone dynamics.

On the convergence of EM-like algorithms for image segmentation using Markov random fields.

Inference of Markov random field images segmentation models is usually performed using iterative methods which adapt the well-known expectation-maximization (EM) algorithm for independent mixture models. However, some of these adaptations are ad hoc and may turn out numerically unstable. In this paper, we review three EM-like variants for Markov random field segmentation and compare their convergence properties both at the theoretical and practical levels. We specifically advocate a numerical scheme involving asynchronous voxel updating, for which general convergence results can be established. Our experiments on brain tissue classification in magnetic resonance images provide evidence that this algorithm may achieve significantly faster convergence than its competitors while yielding at least as good segmentation results.

Ecological-niche factor analysis: how to compute habitat-suitability maps without absence data?

We propose a multivariate approach to the study of geographic species distribution which does not require absence data. Building on Hutchinson's concept of the ecological niche, this factor analysis compares, in the multidimensional space of ecological variables, the distribution of the localities where the focal species was observed to a reference set describing the whole study area. The first factor extracted maximizes the marginality of the focal species, defined as the ecological distance between the species optimum and the mean habitat within the reference area. The other factors maximize the specialization of this focal species, defined as the ratio of the ecological variance in mean habitat to that observed for the focal species. Eigenvectors and eigenvalues are readily interpreted and can be used to build habitat-suitability maps. This approach is recommended in Situations where absence data are not available (many data banks), unreliable (most cryptic or rare species), or meaningless (invaders). We provide an illustration and validation of the method for the alpine ibex, a species reintroduced in Switzerland which presumably has not yet recolonized its entire range.

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.

High-density sex-specific linkage maps of a European tree frog (Hyla arborea) identify the sex chromosome without information on offspring sex.

Identifying homology between sex chromosomes of different species is essential to understanding the evolution of sex determination. Here, we show that the identity of a homomorphic sex chromosome pair can be established using a linkage map, without information on offspring sex. By comparing sex-specific maps of the European tree frog Hyla arborea, we find that the sex chromosome (linkage group 1) shows a threefold difference in marker number between the male and female maps. In contrast, the number of markers on each autosome is similar between the two maps. We also find strongly conserved synteny between H. arborea and Xenopus tropicalis across 200 million years of evolution, suggesting that the rate of chromosomal rearrangement in anurans is low. Finally, we show that recombination in males is greatly reduced at the centers of large chromosomes, consistent with previous cytogenetic findings. Our research shows the importance of high-density linkage maps for studies of recombination, chromosomal rearrangement and the genetic architecture of ecologically or economically important traits.