MorphoGraphX: A platform for quantifying morphogenesis in 4D

Morphogenesis emerges from complex multiscale interactions between genetic and mechanical processes. To understand these processes, the evolution of cell shape, proliferation and gene expression must be quantified. This quantification is usually performed either in full 3D, which is computationally expensive and technically challenging, or on 2D planar projections, which introduces geometrical artifacts on highly curved organs. Here we present MorphoGraphX (www.MorphoGraphX.org), a software that bridges this gap by working directly with curved surface images extracted from 3D data. In addition to traditional 3D image analysis, we have developed algorithms to operate on curved surfaces, such as cell segmentation, lineage tracking and fluorescence signal quantification. The software’s modular design makes it easy to include existing libraries, or to implement new algorithms. Cell geometries extracted with MorphoGraphX can be exported and used as templates for simulation models, providing a powerful platform to investigate the interactions between shape, genes and growth.DOI: http://dx.doi.org/10.7554/eLife.05864.001Author keywordsResearch organism

Hybrid Breakdown in Cichlid Fish

Studies from a wide diversity of taxa have shown a negative relationship between genetic compatibility and the divergence time of hybridizing genomes. Theory predicts the main breakdown of fitness to happen after the F1 hybrid generation, when heterosis subsides and recessive allelic (Dobzhansky-Muller) incompatibilities are increasingly unmasked. We measured the fitness of F2 hybrids of African haplochromine cichlid fish bred from species pairs spanning several thousand to several million years divergence time. F2 hybrids consistently showed the lowest viability compared to F1 hybrids and non-hybrid crosses (crosses within the grandparental species), in agreement with hybrid breakdown. Especially the short- and long-term survival (2 weeks to 6 months) of F2 hybrids was significantly reduced. Overall, F2 hybrids showed a fitness reduction of 21% compared to F1 hybrids, and a reduction of 43% compared to the grandparental, non-hybrid crosses. We further observed a decrease of F2 hybrid viability with the genetic distance between grandparental lineages, suggesting an important role for negative epistatic interactions in cichlid fish postzygotic isolation. The estimated time window for successful production of F2 hybrids resulting from our data is consistent with the estimated divergence time between the multiple ancestral lineages that presumably hybridized in three major adaptive radiations of African cichlids.

Untangling a Holocene pollen record with forest model simulations and independent climate data

Adaptation potential of forests to rapid climatic changes can be assessed from vegetation dynamics during past climatic changes as preserved in fossil pollen data. However, pollen data reflect the integrated effects of climate and biotic processes, such as establishment, survival, competition, and migration. To disentangle these processes, we compared an annually laminated late Würm and Holocene pollen record from the Central Swiss Plateau with simulations of a dynamic forest patch model. All input data used in the simulations were largely independent from pollen data; i.e. the presented analysis is non-circular. Temperature and precipitation scenarios were based on reconstructions from pollen-independent sources. The earliest arrival times of the species at the study site after the last glacial were inferred from pollen maps. We ran a series of simulations under different combinations of climate and immigration scenarios. In addition, the sensitivity of the simulated presence/absence of four major species to changes in the climate scenario was examined. The pattern of the pollen record could partly be explained by the used climate scenario, mostly by temperature. However, some features, in particular the absence of most species during the late Würm could only be simulated if the winter temperature anomalies of the used scenario were decreased considerably. Consequently, we had to assume in the simulations, that most species immigrated during or after the Younger Dryas (12 000 years BP), Abies and Fagus even later. Given the timing of tree species immigration, the vegetation was in equilibrium with climate during long periods, but responded with lags at the time-scale of centuries to millennia caused by a secondary succession after rapid climatic changes such as at the end of Younger Dryas, or immigration of dominant taxa. Climate influenced the tree taxa both directly and indirectly by changing inter-specific competition. We concluded, that also during the present fast climatic change, species migration might be an important process, particularly if geographic barriers, such as the Alps are in the migrational path.

Positron emission mammography in the diagnosis of breast cancer. Is maximum PEM uptake value a valuable threshold for malignant breast cancer detection?

AIM To evaluate the diagnostic value (sensitivity, specificity) of positron emission mammography (PEM) in a single site non-interventional study using the maximum PEM uptake value (PUVmax). PATIENTS, METHODS In a singlesite, non-interventional study, 108 patients (107 women, 1 man) with a total of 151 suspected lesions were scanned with a PEM Flex Solo II (Naviscan) at 90 min p.i. with 3.5 MBq 18F-FDG per kg of body weight. In this ROI(region of interest)-based analysis, maximum PEM uptake value (PUV) was determined in lesions, tumours (PUVmaxtumour), benign lesions (PUVmaxnormal breast) and also in healthy tissues on the contralateral side (PUVmaxcontralateral breast). These values were compared and contrasted. In addition, the ratios of PUVmaxtumour / PUVmaxcontralateral breast and PUVmaxnormal breast / PUVmaxcontralateral breast were compared. The image data were interpreted independently by two experienced nuclear medicine physicians and compared with histology in cases of suspected carcinoma. RESULTS Based on a criteria of PUV>1.9, 31 out of 151 lesions in the patient cohort were found to be malignant (21%). A mean PUVmaxtumour of 3.78 ± 2.47 was identified in malignant tumours, while a mean PUVmaxnormal breast of 1.17 ± 0.37 was reported in the glandular tissue of the healthy breast, with the difference being statistically significant (p < 0.001). Similarly, the mean ratio between tumour and healthy glandular tissue in breast cancer patients (3.15 ± 1.58) was found to be significantly higher than the ratio for benign lesions (1.17 ± 0.41, p < 0.001). CONCLUSION PEM is capable of differentiating breast tumours from benign lesions with 100% sensitivity along with a high specificity of 96%, when a threshold of PUVmax >1.9 is applied.

Investigating the potential of reported cattle mortality data in Switzerland for syndromic surveillance

Systems for the identification and registration of cattle have gradually been receiving attention for use in syndromic surveillance, a relatively recent approach for the early detection of infectious disease outbreaks. Real or near real-time monitoring of deaths or stillbirths reported to these systems offer an opportunity to detect temporal or spatial clusters of increased mortality that could be caused by an infectious disease epidemic. In Switzerland, such data are recorded in the "Tierverkehrsdatenbank" (TVD). To investigate the potential of the Swiss TVD for syndromic surveillance, 3 years of data (2009-2011) were assessed in terms of data quality, including timeliness of reporting and completeness of geographic data. Two time-series consisting of reported on-farm deaths and stillbirths were retrospectively analysed to define and quantify the temporal patterns that result from non-health related factors. Geographic data were almost always present in the TVD data; often at different spatial scales. On-farm deaths were reported to the database by farmers in a timely fashion; stillbirths were less timely. Timeliness and geographic coverage are two important features of disease surveillance systems, highlighting the suitability of the TVD for use in a syndromic surveillance system. Both time series exhibited different temporal patterns that were associated with non-health related factors. To avoid false positive signals, these patterns need to be removed from the data or accounted for in some way before applying aberration detection algorithms in real-time. Evaluating mortality data reported to systems for the identification and registration of cattle is of value for comparing national data systems and as a first step towards a European-wide early detection system for emerging and re-emerging cattle diseases.

Natural History of Geographic Atrophy Progression Secondary to Age-Related Macular Degeneration (Geographic Atrophy Progression Study).

PURPOSE The Geographic Atrophy Progression (GAP) study was designed to assess the rate of geographic atrophy (GA) progression and to identify prognostic factors by measuring the enlargement of the atrophic lesions using fundus autofluorescence (FAF) and color fundus photography (CFP). DESIGN Prospective, multicenter, noninterventional natural history study. PARTICIPANTS A total of 603 participants were enrolled in the study; 413 of those had gradable lesion data from FAF or CFP, and 321 had gradable lesion data from both FAF and CFP. METHODS Atrophic lesion areas were measured by FAF and CFP to assess lesion progression over time. Lesion size assessments and best-corrected visual acuity (BCVA) were conducted at screening/baseline (day 0) and at 3 follow-up visits: month 6, month 12, and month 18 (or early exit). MAIN OUTCOME MEASURES The GA lesion progression rate in disease subgroups and mean change from baseline visual acuity. RESULTS Mean (standard error) lesion size changes from baseline, determined by FAF and CFP, respectively, were 0.88 (0.1) and 0.78 (0.1) mm(2) at 6 months, 1.85 (0.1) and 1.57 (0.1) mm(2) at 12 months, and 3.14 (0.4) and 3.17 (0.5) mm(2) at 18 months. The mean change in lesion size from baseline to month 12 was significantly greater in participants who had eyes with multifocal atrophic spots compared with those with unifocal spots (P < 0.001) and those with extrafoveal lesions compared with those with foveal lesions (P = 0.001). The mean (standard deviation) decrease in visual acuity was 6.2 ± 15.6 letters for patients with image data available. Atrophic lesions with a diffuse (mean 0.95 mm(2)) or banded (mean 1.01 mm(2)) FAF pattern grew more rapidly by month 6 compared with those with the "none" (mean, 0.13 mm(2)) and focal (mean, 0.36 mm(2)) FAF patterns. CONCLUSIONS Although differences were observed in mean lesion size measurements using FAF imaging compared with CFP, the measurements were highly correlated with one another. Significant differences were found in lesion progression rates in participants stratified by hyperfluorescence pattern subtype. This large GA natural history study provides a strong foundation for future clinical trials.

A Novel Ultrasound-Based Registration for Image-Guided Laparoscopic Liver Ablation.

BACKGROUND Patient-to-image registration is a core process of image-guided surgery (IGS) systems. We present a novel registration approach for application in laparoscopic liver surgery, which reconstructs in real time an intraoperative volume of the underlying intrahepatic vessels through an ultrasound (US) sweep process. METHODS An existing IGS system for an open liver procedure was adapted, with suitable instrument tracking for laparoscopic equipment. Registration accuracy was evaluated on a realistic phantom by computing the target registration error (TRE) for 5 intrahepatic tumors. The registration work flow was evaluated by computing the time required for performing the registration. Additionally, a scheme for intraoperative accuracy assessment by visual overlay of the US image with preoperative image data was evaluated. RESULTS The proposed registration method achieved an average TRE of 7.2 mm in the left lobe and 9.7 mm in the right lobe. The average time required for performing the registration was 12 minutes. A positive correlation was found between the intraoperative accuracy assessment and the obtained TREs. CONCLUSIONS The registration accuracy of the proposed method is adequate for laparoscopic intrahepatic tumor targeting. The presented approach is feasible and fast and may, therefore, not be disruptive to the current surgical work flow.

Trends in mean and extreme precipitation in the Mount Kenya region from observations and reanalyses

Land and water management in semi-arid regions requires detailed information on precipitation distribution, including extremes, and changes therein. Such information is often lacking. This paper describes statistics of mean and extreme precipitation in a unique data set from the Mount Kenya region, encompassing around 50 stations with at least 30 years of data. We describe the data set, including quality control procedures and statistical break detection. Trends in mean precipitation and extreme indices calculated from these data for individual rainy seasons are compared with corresponding trends in reanalysis products. From 1979 to 2011, mean precipitation decreased at 75% of the stations during the ‘long rains’ (March to May) and increased at 70% of the stations during the ‘short rains’ (October to December). Corresponding trends are found in the number of heavy precipitation days, and maximum of consecutive 5-day precipitation. Conversely, an increase in consecutive dry days within both main rainy seasons is found. However, trends are only statistically significant in very few cases. Reanalysis data sets agree with observations with respect to interannual variability, while correlations are considerably lower for monthly deviations (ratios) from the mean annual cycle. While some products well reproduce the rainfall climatology and some the spatial trend pattern, no product reproduces both.

Origin and history of ureilitic material in the solar system: The view from asteroid 2008 TC3 and the Almahata Sitta meteorite

Asteroid 2008 TC3 (approximately 4m diameter) was tracked and studied in space for approximately 19h before it impacted Earth's atmosphere, shattering at 44-36km altitude. The recovered samples (>680 individual rocks) comprise the meteorite Almahata Sitta (AhS). Approximately 50-70% of these are ureilites (ultramafic achondrites). The rest are chondrites, mainly enstatite, ordinary, and Rumuruti types. The goal of this work is to understand how fragments of so many different types of parent bodies became mixed in the same asteroid. Almahata Sitta has been classified as a polymict ureilite with an anomalously high component of foreign clasts. However, we calculate that the mass of fallen material was 0.1% of the pre-atmospheric mass of the asteroid. Based on published data for the reflectance spectrum of the asteroid and laboratory spectra of the samples, we infer that the lost material was mostly ureilitic. Therefore, 2008 TC3 probably contained only a few percent nonureilitic materials, similar to other polymict ureilites except less well consolidated. From available data for the AhS meteorite fragments, we conclude that 2008 TC3 samples essentially the same range of types of ureilitic and nonureilitic materials as other polymict ureilites. We therefore suggest that the immediate parent of 2008 TC3 was the immediate parent of all ureilitic material sampled on Earth. We trace critical stages in the evolution of that material through solar system history. Based on various types of new modeling and re-evaluation of published data, we propose the following scenario. (1) The ureilite parent body (UPB) accreted 0.5-0.6Ma after formation of calcium-aluminum-rich inclusions (CAI), beyond the ice line (outer asteroid belt). Differentiation began approximately 1Ma after CAI. (2) The UPB was catastrophically disrupted by a major impact approximately 5Ma after CAI, with selective subsets of the fragments reassembling into daughter bodies. (3) Either the UPB (before breakup), or one of its daughters (after breakup), migrated to the inner belt due to scattering by massive embryos. (4) One daughter (after forming in or migrating to the inner belt) became the parent of 2008 TC3. It developed a regolith, mostly 3.8Ga ago. Clasts of enstatite, ordinary, and Rumuruti-type chondrites were implanted by low-velocity collisions. (5) Recently, the daughter was disrupted. Fragments were injected or drifted into Earth-crossing orbits. 2008 TC3 comes from outer layers of regolith, other polymict ureilites from deeper regolith, and main group ureilites from the interior of this body. In contrast to other models that have been proposed, this model invokes a stochastic history to explain the unique diversity of foreign materials in 2008 TC3 and other polymict ureilites.

Extracting more out of relocation data: building movement models as mixtures of random walks

We present a framework for fitting multiple random walks to animal movement paths consisting of ordered sets of step lengths and turning angles. Each step and turn is assigned to one of a number of random walks, each characteristic of a different behavioral state. Behavioral state assignments may be inferred purely from movement data or may include the habitat type in which the animals are located. Switching between different behavioral states may be modeled explicitly using a state transition matrix estimated directly from data, or switching probabilities may take into account the proximity of animals to landscape features. Model fitting is undertaken within a Bayesian framework using the WinBUGS software. These methods allow for identification of different movement states using several properties of observed paths and lead naturally to the formulation of movement models. Analysis of relocation data from elk released in east-central Ontario, Canada, suggests a biphasic movement behavior: elk are either in an "encamped" state in which step lengths are small and turning angles are high, or in an "exploratory" state, in which daily step lengths are several kilometers and turning angles are small. Animals encamp in open habitat (agricultural fields and opened forest), but the exploratory state is not associated with any particular habitat type.

Mixture modeling for joint analysis of survival, discrete, and continuous data

Mixture modeling is commonly used to model categorical latent variables that represent subpopulations in which population membership is unknown but can be inferred from the data. In relatively recent years, the potential of finite mixture models has been applied in time-to-event data. However, the commonly used survival mixture model assumes that the effects of the covariates involved in failure times differ across latent classes, but the covariate distribution is homogeneous. The aim of this dissertation is to develop a method to examine time-to-event data in the presence of unobserved heterogeneity under a framework of mixture modeling. A joint model is developed to incorporate the latent survival trajectory along with the observed information for the joint analysis of a time-to-event variable, its discrete and continuous covariates, and a latent class variable. It is assumed that the effects of covariates on survival times and the distribution of covariates vary across different latent classes. The unobservable survival trajectories are identified through estimating the probability that a subject belongs to a particular class based on observed information. We applied this method to a Hodgkin lymphoma study with long-term follow-up and observed four distinct latent classes in terms of long-term survival and distributions of prognostic factors. Our results from simulation studies and from the Hodgkin lymphoma study demonstrated the superiority of our joint model compared with the conventional survival model. This flexible inference method provides more accurate estimation and accommodates unobservable heterogeneity among individuals while taking involved interactions between covariates into consideration.^

n-Alkane analyses from an transect of marine surface sediments off southwest Africa

An isobathic transect of marine surface sediments from 1°N to 28°S off southwest Africa was used to further evaluate the potential of the chain length distribution and carbon stable isotope composition of higher plant n-alkanes as proxies for continental vegetation and climate conditions. We found a strong increase in the n-C29-33 weighted mean average d13C values from -33 per mil near the equator to around -26 per mil further south. Additionally, C25-35n-alkanes reveal a southward trend of increasing average chain length from 30.0 to 30.5. The data reflect the changing contribution of plants employing different photosynthetic pathways (C3 and C4) and/or being differently influenced by the environmental conditions of their habitat. The C4 plant proportions calculated from the data (ca. 20% for rivers draining the rainforest, to ca. 70% at higher latitude) correspond to the C4 plant abundance in continental catchment areas postulated by considering prevailing wind systems and river outflows. Furthermore, the C4 plant contribution to the sediments correlates with the mean annual precipitation and aridity at selected continental locations in the postulated catchment areas, suggesting that the C4 plant fraction in marine sediments can be used to assess these environmental parameters.