The far-infrared/submillimeter properties of galaxies located behind the Bullet cluster

The Herschel Lensing Survey (HLS) takes advantage of gravitational lensing by massive galaxy clusters to sample a population of high-redshift galaxies which are too faint to be detected above the confusion limit of current far-infrared/submillimeter telescopes. Measurements from 100-500 μm bracket the peaks of the far-infrared spectral energy distributions of these galaxies, characterizing their infrared luminosities and star formation rates. We introduce initial results from our science demonstration phase observations, directed toward the Bullet cluster (1E0657-56). By combining our observations with LABOCA 870 μm and AzTEC 1.1 mm data we fully constrain the spectral energy distributions of 19 MIPS 24 μm-selected galaxies which are located behind the cluster. We find that their colors are best fit using templates based on local galaxies with systematically lower infrared luminosities. This suggests that our sources are not like local ultra-luminous infrared galaxies in which vigorous star formation is contained in a compact highly dust-obscured region. Instead, they appear to be scaled up versions of lower luminosity local galaxies with star formation occurring on larger physical scales.

Paradigm shift in determining Neoproterozoic atmospheric oxygen

ACKNOWLEDGMENTS We thank the Geological Survey of Australia for permission to sample the Empress 1A and Lancer 1 cores, the Natural Sciences and Engineering Research Council of Canada for financial support (grant #7961–15) of U. Brand, and the National Natural Science Foundation of China for support of F. Meng and P. Ni (grants 41473039 and 4151101015). We thank M. Lozon (Brock University) for drafting and constructing the figures. We thank the editor, Brendan Murphy, as well as three reviewers (Steve Kesler, Erik Sperling, and an anonymous reviewer), for improving the manuscript into its final form ©The Authors Gold Open Access: This paper is published under the terms of the CC-BY license.

Radiotherapy Treatment Assessment using DCE-MRI

Abstract

The goal of modern radiotherapy is to precisely deliver a prescribed radiation dose to delineated target volumes that contain a significant amount of tumor cells while sparing the surrounding healthy tissues/organs. Precise delineation of treatment and avoidance volumes is the key for the precision radiation therapy. In recent years, considerable clinical and research efforts have been devoted to integrate MRI into radiotherapy workflow motivated by the superior soft tissue contrast and functional imaging possibility. Dynamic contrast-enhanced MRI (DCE-MRI) is a noninvasive technique that measures properties of tissue microvasculature. Its sensitivity to radiation-induced vascular pharmacokinetic (PK) changes has been preliminary demonstrated. In spite of its great potential, two major challenges have limited DCE-MRI’s clinical application in radiotherapy assessment: the technical limitations of accurate DCE-MRI imaging implementation and the need of novel DCE-MRI data analysis methods for richer functional heterogeneity information.

This study aims at improving current DCE-MRI techniques and developing new DCE-MRI analysis methods for particular radiotherapy assessment. Thus, the study is naturally divided into two parts. The first part focuses on DCE-MRI temporal resolution as one of the key DCE-MRI technical factors, and some improvements regarding DCE-MRI temporal resolution are proposed; the second part explores the potential value of image heterogeneity analysis and multiple PK model combination for therapeutic response assessment, and several novel DCE-MRI data analysis methods are developed.

I. Improvement of DCE-MRI temporal resolution. First, the feasibility of improving DCE-MRI temporal resolution via image undersampling was studied. Specifically, a novel MR image iterative reconstruction algorithm was studied for DCE-MRI reconstruction. This algorithm was built on the recently developed compress sensing (CS) theory. By utilizing a limited k-space acquisition with shorter imaging time, images can be reconstructed in an iterative fashion under the regularization of a newly proposed total generalized variation (TGV) penalty term. In the retrospective study of brain radiosurgery patient DCE-MRI scans under IRB-approval, the clinically obtained image data was selected as reference data, and the simulated accelerated k-space acquisition was generated via undersampling the reference image full k-space with designed sampling grids. Two undersampling strategies were proposed: 1) a radial multi-ray grid with a special angular distribution was adopted to sample each slice of the full k-space; 2) a Cartesian random sampling grid series with spatiotemporal constraints from adjacent frames was adopted to sample the dynamic k-space series at a slice location. Two sets of PK parameters’ maps were generated from the undersampled data and from the fully-sampled data, respectively. Multiple quantitative measurements and statistical studies were performed to evaluate the accuracy of PK maps generated from the undersampled data in reference to the PK maps generated from the fully-sampled data. Results showed that at a simulated acceleration factor of four, PK maps could be faithfully calculated from the DCE images that were reconstructed using undersampled data, and no statistically significant differences were found between the regional PK mean values from undersampled and fully-sampled data sets. DCE-MRI acceleration using the investigated image reconstruction method has been suggested as feasible and promising.

Second, for high temporal resolution DCE-MRI, a new PK model fitting method was developed to solve PK parameters for better calculation accuracy and efficiency. This method is based on a derivative-based deformation of the commonly used Tofts PK model, which is presented as an integrative expression. This method also includes an advanced Kolmogorov-Zurbenko (KZ) filter to remove the potential noise effect in data and solve the PK parameter as a linear problem in matrix format. In the computer simulation study, PK parameters representing typical intracranial values were selected as references to simulated DCE-MRI data for different temporal resolution and different data noise level. Results showed that at both high temporal resolutions (<1s) and clinically feasible temporal resolution (~5s), this new method was able to calculate PK parameters more accurate than the current calculation methods at clinically relevant noise levels; at high temporal resolutions, the calculation efficiency of this new method was superior to current methods in an order of 102. In a retrospective of clinical brain DCE-MRI scans, the PK maps derived from the proposed method were comparable with the results from current methods. Based on these results, it can be concluded that this new method can be used for accurate and efficient PK model fitting for high temporal resolution DCE-MRI.

II. Development of DCE-MRI analysis methods for therapeutic response assessment. This part aims at methodology developments in two approaches. The first one is to develop model-free analysis method for DCE-MRI functional heterogeneity evaluation. This approach is inspired by the rationale that radiotherapy-induced functional change could be heterogeneous across the treatment area. The first effort was spent on a translational investigation of classic fractal dimension theory for DCE-MRI therapeutic response assessment. In a small-animal anti-angiogenesis drug therapy experiment, the randomly assigned treatment/control groups received multiple fraction treatments with one pre-treatment and multiple post-treatment high spatiotemporal DCE-MRI scans. In the post-treatment scan two weeks after the start, the investigated Rényi dimensions of the classic PK rate constant map demonstrated significant differences between the treatment and the control groups; when Rényi dimensions were adopted for treatment/control group classification, the achieved accuracy was higher than the accuracy from using conventional PK parameter statistics. Following this pilot work, two novel texture analysis methods were proposed. First, a new technique called Gray Level Local Power Matrix (GLLPM) was developed. It intends to solve the lack of temporal information and poor calculation efficiency of the commonly used Gray Level Co-Occurrence Matrix (GLCOM) techniques. In the same small animal experiment, the dynamic curves of Haralick texture features derived from the GLLPM had an overall better performance than the corresponding curves derived from current GLCOM techniques in treatment/control separation and classification. The second developed method is dynamic Fractal Signature Dissimilarity (FSD) analysis. Inspired by the classic fractal dimension theory, this method measures the dynamics of tumor heterogeneity during the contrast agent uptake in a quantitative fashion on DCE images. In the small animal experiment mentioned before, the selected parameters from dynamic FSD analysis showed significant differences between treatment/control groups as early as after 1 treatment fraction; in contrast, metrics from conventional PK analysis showed significant differences only after 3 treatment fractions. When using dynamic FSD parameters, the treatment/control group classification after 1st treatment fraction was improved than using conventional PK statistics. These results suggest the promising application of this novel method for capturing early therapeutic response.

The second approach of developing novel DCE-MRI methods is to combine PK information from multiple PK models. Currently, the classic Tofts model or its alternative version has been widely adopted for DCE-MRI analysis as a gold-standard approach for therapeutic response assessment. Previously, a shutter-speed (SS) model was proposed to incorporate transcytolemmal water exchange effect into contrast agent concentration quantification. In spite of richer biological assumption, its application in therapeutic response assessment is limited. It might be intriguing to combine the information from the SS model and from the classic Tofts model to explore potential new biological information for treatment assessment. The feasibility of this idea was investigated in the same small animal experiment. The SS model was compared against the Tofts model for therapeutic response assessment using PK parameter regional mean value comparison. Based on the modeled transcytolemmal water exchange rate, a biological subvolume was proposed and was automatically identified using histogram analysis. Within the biological subvolume, the PK rate constant derived from the SS model were proved to be superior to the one from Tofts model in treatment/control separation and classification. Furthermore, novel biomarkers were designed to integrate PK rate constants from these two models. When being evaluated in the biological subvolume, this biomarker was able to reflect significant treatment/control difference in both post-treatment evaluation. These results confirm the potential value of SS model as well as its combination with Tofts model for therapeutic response assessment.

In summary, this study addressed two problems of DCE-MRI application in radiotherapy assessment. In the first part, a method of accelerating DCE-MRI acquisition for better temporal resolution was investigated, and a novel PK model fitting algorithm was proposed for high temporal resolution DCE-MRI. In the second part, two model-free texture analysis methods and a multiple-model analysis method were developed for DCE-MRI therapeutic response assessment. The presented works could benefit the future DCE-MRI routine clinical application in radiotherapy assessment.

Wear and its effects on dental topography measures in howling monkeys (Alouatta palliata).

OBJECTIVES: Three dental topography measurements: Dirichlet Normal Energy (DNE), Relief Index (RFI), and Orientation Patch Count Rotated (OPCR) are examined for their interaction with measures of wear, within and between upper and lower molars in Alouatta palliata. Potential inferences of the "dental sculpting" phenomenon are explored. MATERIALS AND METHODS: Fifteen occluding pairs of howling monkey first molars (15 upper, 15 lower) opportunistically collected from La Pacifica, Costa Rica, were selected to sample wear stages ranging from unworn to heavily worn as measured by the Dentine Exposure Ratio (DER). DNE, RFI, and OPCR were measured from three-dimensional surface reconstructions (PLY files) derived from high-resolution CT scans. Relationships among the variables were tested with regression analyses. RESULTS: Upper molars have more cutting edges, exhibiting significantly higher DNE, but have significantly lower RFI values. However, the relationships among the measures are concordant across both sets of molars. DER and EDJL are curvilinearly related. DER is positively correlated with DNE, negatively correlated with RFI, and uncorrelated with OPCR. EDJL is not correlated with DNE, or RFI, but is positively correlated with OPCR among lower molars only. DISCUSSION: The relationships among these metrics suggest that howling monkey teeth adaptively engage macrowear. DNE increases with wear in this sample presumably improving food breakdown. RFI is initially high but declines with wear, suggesting that the initially high RFI safeguards against dental senescence. OPCR values in howling monkey teeth do not show a clear relationship with wear changes.

Le décours temporel de l'utilisation des fréquences spatiales dans les troubles du spectre autistique

Notre système visuel extrait d'ordinaire l'information en basses fréquences spatiales (FS) avant celles en hautes FS. L'information globale extraite tôt peut ainsi activer des hypothèses sur l'identité de l'objet et guider l'extraction d'information plus fine spécifique par la suite. Dans les troubles du spectre autistique (TSA), toutefois, la perception des FS est atypique. De plus, la perception des individus atteints de TSA semble être moins influencée par leurs a priori et connaissances antérieures. Dans l'étude décrite dans le corps de ce mémoire, nous avions pour but de vérifier si l'a priori de traiter l'information des basses aux hautes FS était présent chez les individus atteints de TSA. Nous avons comparé le décours temporel de l'utilisation des FS chez des sujets neurotypiques et atteints de TSA en échantillonnant aléatoirement et exhaustivement l'espace temps x FS. Les sujets neurotypiques extrayaient les basses FS avant les plus hautes: nous avons ainsi pu répliquer le résultat de plusieurs études antérieures, tout en le caractérisant avec plus de précision que jamais auparavant. Les sujets atteints de TSA, quant à eux, extrayaient toutes les FS utiles, basses et hautes, dès le début, indiquant qu'ils ne possédaient pas l'a priori présent chez les neurotypiques. Il semblerait ainsi que les individus atteints de TSA extraient les FS de manière purement ascendante, l'extraction n'étant pas guidée par l'activation d'hypothèses.

Impact of ocean acidification and warming on the larval development of the spider crab Hyas araneus from different latitudes (54° vs 79°N), 2010

The combined impacts of future scenarios of ocean acidification and global warming on the larvae of a cold-eurythermal spider crab, Hyas araneus L., were investigated in one of its southernmost populations (living around Helgoland, southern North Sea, 54°N) and one of the northernmost populations (Svalbard, North Atlantic, 79°N). Larvae were exposed at temperatures of 3, 9 and 15°C to present day normocapnia (380 ppm CO2) and to CO2 conditions expected for the near or medium-term future (710 ppm by 2100 and 3000 ppm CO2 by 2300 and beyond). Larval development time and biochemical composition were studied in the larval stages Zoea I, II, and Megalopa. Permanent differences in instar duration between both populations were detected in all stages, likely as a result of evolutionary temperature adaptation. With the exception of Zoea II at 3°C and under all CO2 conditions, development in all instars from Svalbard was delayed compared to those from Helgoland, under all conditions. Most prominently, development was much longer and fewer specimens morphosed to the first crab instar in the Megalopa from Svalbard than from Helgoland. Enhanced CO2 levels (710 and particularly 3000 ppm), caused extended duration of larval development and reduced larval growth (measured as dry mass) and fitness (decreasing C/N ratio, a proxy of the lipid content). Such effects were strongest in the zoeal stages in Svalbard larvae, and during the Megalopa instar in Helgoland larvae.

Annotated record of the detailed examination of Mn deposits from DSDP Hole 2-8

The Atlantic Advisory Panel proposed that Site 8 should be drilled on the rise between the Hatteras and Sohm Abyssal Plains (lat 35° 2l'N., long 67° 3l'W.) This location was considered to offer the best opportunity for realizing two primary objectives. The first of these objectives was to sample and date the oldest available rock in a region adjacent to the North American continent and as far as possible from the Mid-Atlantic Ridge. The second objective relates to the potential paleobiological and paleoecological information that could be derived from the sedimentary column in this general area.

Pore water chemistry of samples near the sediment-water interface

Rhizon samplers were originally designed as micro-tensiometers for soil science to sample seepage water in the unsaturated zone. This study shows applications of Rhizons for porewater sampling from sediments in aquatic systems and presents a newly developed Rhizon in situ sampler (RISS). With the inexpensive Rhizon sampling technique, porewater profiles can be sampled with minimum disturbance of both the sediment structure and possible flow fields. Field experiments, tracer studies, and numerical modeling were combined to assess the suitability of Rhizons for porewater sampling. It is shown that the low effort and simple application makes Rhizons a powerful tool for porewater sampling and an alternative to classical methods. Our investigations show that Rhizons are well suited for sampling porewater on board a ship, in a laboratory, and also for in situ sampling. The results revealed that horizontally aligned Rhizons can sample porewater with a vertical resolution of 1 cm. Combined with an in situ benthic chamber system, the RISS allows studies of benthic fluxes and porewater profiles at the same location on the seafloor with negligible effect on the incubated sediment water interface. Results derived by porewater sampling of sediment cores from the Southern Ocean (Atlantic sector) and by in situ sampling of tidal flat sediments of the Wadden Sea (Sahlenburg/Cuxhaven, Germany) are presented.

Molecular hydrogen (H2) mixing ratio and stable isotopic composition (dD) at the Cabauw tall tower in the Netherlands (2008-2012)

Measurements of the stable isotopic composition (dD(H2) or dD) of atmospheric molecular hydrogen (H2) are a useful addition to mixing ratio (X(H2)) measurements for understanding the atmospheric H2 cycle. dD datasets published so far consist mostly of observations at background locations. We complement these with observations from the Cabauw tall tower at the CESAR site, situated in a densely populated region of the Netherlands. Our measurements show a large anthropogenic influence on the local H2 cycle, with frequently occurring pollution events that are characterized by X(H2) values that reach up to 1 ppm and low dD values. An isotopic source signature analysis yields an apparent source signature below -400 per mil, which is much more D-depleted than the fossil fuel combustion source signature commonly used in H2 budget studies. Two diurnal cycles that were sampled at a suburban site near London also show a more D-depleted source signature (-340 per mil), though not as extremely depleted as at Cabauw. The source signature of the Northwest European vehicle fleet may have shifted to somewhat lower values due to changes in vehicle technology and driving conditions. Even so, the surprisingly depleted apparent source signature at Cabauw requires additional explanation; microbial H2 production seems the most likely cause. The Cabauw tower site also allowed us to sample vertical profiles. We found no decrease in (H2) at lower sampling levels (20 and 60m) with respect to higher sampling levels (120 and 200m). There was a significant shift to lower median dD values at the lower levels. This confirms the limited role of soil uptake around Cabauw, and again points to microbial H2 production during an extended growing season, as well as to possible differences in average fossil fuel combustion source signature between the different footprint areas of the sampling levels. So, although knowledge of the background cycle of H2 has improved over the last decade, surprising features come to light when a non-background location is studied, revealing remaining gaps in our understanding.

Annotated record of the detailed examination of Mn deposits from DSDP Site 42, Leg 5 (Cores 5-42-5, 5-42-6, 5-42-9)

Site 42 is one of the series of sites selected by the Pacific Advisory Panel along the 140th meridian to sample the longitudinal variation in sediment composition in the eastern Pacific. The site is located in an area of abyssal hills between the Clarion and Clipperton Fracture Zones, and is at the northern margin of the thick development of acoustically transparent sediment extending along the equator.

Annotated record of the detailed examination of Mn deposits from DSDP Site 25, Leg 4 (Core 4-25A-1)

At the Leg 4 pre-cruise meeting in Miami, Florida, in January 1969, the Ridge site was proposed as an additional site to be investigated, if possible. Since the water depth and sediment thickness at this site were both minimal (2000 meters and less than 100 meters, respectively) and the location was not far from the track between initially planned sites, the shipboard party decided, upon completion of work at Site 24, to make a 24-hour drilling effort on the North Brazilian Ridge. The drilling objectives at this location were to sample and date the sediments atop the Ridge crest and to determine the nature, age and origin of the consolidated material beneath the thin mantle of sediment.

On the multi-dimensionality and sampling of air transport networks

Complex network theory is a framework increasingly used in the study of air transport networks, thanks to its ability to describe the structures created by networks of flights, and their influence in dynamical processes such as delay propagation. While many works consider only a fraction of the network, created by major airports or airlines, for example, it is not clear if and how such sampling process bias the observed structures and processes. In this contribution, we tackle this problem by studying how some observed topological metrics depend on the way the network is reconstructed, i.e. on the rules used to sample nodes and connections. Both structural and simple dynamical properties are considered, for eight major air networks and different source datasets. Results indicate that using a subset of airports strongly distorts our perception of the network, even when just small ones are discarded; at the same time, considering a subset of airlines yields a better and more stable representation. This allows us to provide some general guidelines on the way airports and connections should be sampled.

Type Ia supernovae from violent mergers of carbon-oxygen white dwarfs: polarization signatures

The violent merger of two carbon-oxygen white dwarfs has been proposed as a viable progenitor for some Type Ia supernovae. However, it has been argued that the strong ejecta asymmetries produced by this model might be inconsistent with the low degree of polarization typically observed in Type Ia supernova explosions. Here, we test this claim by carrying out a spectropolarimetric analysis for the model proposed by Pakmor et al. for an explosion triggered during the merger of a 1.1 and 0.9 M⊙ carbon-oxygen white dwarf binary system. Owing to the asymmetries of the ejecta, the polarization signal varies significantly with viewing angle. We find that polarization levels for observers in the equatorial plane are modest (≲1 per cent) and show clear evidence for a dominant axis, as a consequence of the ejecta symmetry about the orbital plane. In contrast, orientations out of the plane are associated with higher degrees of polarization and departures from a dominant axis. While the particular model studied here gives a good match to highly polarized events such as SN 2004dt, it has difficulties in reproducing the low polarization levels commonly observed in normal Type Ia supernovae. Specifically, we find that significant asymmetries in the element distribution result in a wealth of strong polarization features that are not observed in the majority of currently available spectropolarimetric data of Type Ia supernovae. Future studies will map out the parameter space of the merger scenario to investigate if alternative models can provide better agreement with observations.

The peculiar Type Ia supernova iPTF14atg: Chandrasekhar-mass explosion or violent merger?

iPTF14atg, a subluminous peculiar Type Ia supernova (SN Ia) similar to SN 2002es, is the first SN Ia for which a strong UV flash was observed in the early-time light curves. This has been interpreted as evidence for a single-degenerate (SD) progenitor system, where such a signal is expected from interactions between the SN ejecta and the non-degenerate companion star. Here, we compare synthetic observables of multidimensional state-of-the-art explosion models for different progenitor scenarios to the light curves and spectra of iPTF14atg. From our models, we have difficulties explaining the spectral evolution of iPTF14atg within the SD progenitor channel. In contrast, we find that a violent merger of two carbon-oxygen white dwarfs with 0.9 and 0.76 M⊙, respectively, provides an excellent match to the spectral evolution of iPTF14atg from 10 d before to several weeks after maximum light. Our merger model does not naturally explain the initial UV flash of iPTF14atg. We discuss several possibilities like interactions of the SN ejecta with the circumstellar medium and surface radioactivity from an He-ignited merger that may be able to account for the early UV emission in violent merger models.

Building a 'Repository of Science': the importance of integrating Biobanks within molecular pathology programmes

Repositories containing high quality human biospecimens linked with robust and relevant clinical and pathological information are required for the discovery and validation of biomarkers for disease diagnosis, progression and response to treatment. Current molecular based discovery projects using either low or high throughput technologies rely heavily on ready access to such sample collections. It is imperative that modern biobanks align with molecular diagnostic pathology practices not only to provide the type of samples needed for discovery projects but also to ensure requirements for ongoing sample collections and the future needs of researchers are adequately addressed. Biobanks within comprehensive molecular pathology programmes are perfectly positioned to offer more than just tumour derived biospecimens; for example, they have the ability to facilitate researchers gaining access to sample metadata such as digitised scans of tissue samples annotated prior to macrodissection for molecular diagnostics or pseudoanonymised clinical outcome data or research results retrieved from other users utilising the same or overlapping cohorts of samples. Furthermore, biobanks can work with molecular diagnostic laboratories to develop standardized methodologies for the acquisition and storage of samples required for new approaches to research such as ‘liquid biopsies’ which will ultimately feed into the test validations required in large prospective clinical studies in order to implement liquid biopsy approaches for routine clinical practice. We draw on our experience in Northern Ireland to discuss how this harmonised approach of biobanks working synergistically with molecular pathology programmes is key for the future success of precision medicine.