News from compositions, the R package

The R-package “compositions”is a tool for advanced compositional analysis. Its basicfunctionality has seen some conceptual improvement, containing now some facilitiesto work with and represent ilr bases built from balances, and an elaborated subsys-tem for dealing with several kinds of irregular data: (rounded or structural) zeroes,incomplete observations and outliers. The general approach to these irregularities isbased on subcompositions: for an irregular datum, one can distinguish a “regular” sub-composition (where all parts are actually observed and the datum behaves typically)and a “problematic” subcomposition (with those unobserved, zero or rounded parts, orelse where the datum shows an erratic or atypical behaviour). Systematic classificationschemes are proposed for both outliers and missing values (including zeros) focusing onthe nature of irregularities in the datum subcomposition(s).To compute statistics with values missing at random and structural zeros, a projectionapproach is implemented: a given datum contributes to the estimation of the desiredparameters only on the subcompositon where it was observed. For data sets withvalues below the detection limit, two different approaches are provided: the well-knownimputation technique, and also the projection approach.To compute statistics in the presence of outliers, robust statistics are adapted to thecharacteristics of compositional data, based on the minimum covariance determinantapproach. The outlier classification is based on four different models of outlier occur-rence and Monte-Carlo-based tests for their characterization. Furthermore the packageprovides special plots helping to understand the nature of outliers in the dataset.Keywords: coda-dendrogram, lost values, MAR, missing data, MCD estimator,robustness, rounded zeros

Compositional data and Simpson's paradox

Simpson's paradox, also known as amalgamation or aggregation paradox, appears whendealing with proportions. Proportions are by construction parts of a whole, which canbe interpreted as compositions assuming they only carry relative information. TheAitchison inner product space structure of the simplex, the sample space of compositions, explains the appearance of the paradox, given that amalgamation is a nonlinearoperation within that structure. Here we propose to use balances, which are specificelements of this structure, to analyse situations where the paradox might appear. Withthe proposed approach we obtain that the centre of the tables analysed is a naturalway to compare them, which avoids by construction the possibility of a paradox.Key words: Aitchison geometry, geometric mean, orthogonal projection

The cover basement contact beneath the Rawil axial depression (Western Alps) - True amplitude seismic processing, petrophysical properties, and modeling

Since 1986, several near-vertical seismic reflection profiles have been recorded in Switzerland in order to map the deep geologic structure of the Alps. One objective of this endeavour has been to determine the geometries of the autochthonous basement and of the external crystalline massifs, important elements for understanding the geodynamics of the Alpine orogeny. The PNR-20 seismic line W1, located in the Rawil depression of the western Swiss Alps, provides important information on this subject. It extends northward from the `'Penninic front'' across the Helvetic nappes to the Prealps. The crystalline massifs do not outcrop along this profile. Thus, the interpretation of `'near-basement'' reflections has to be constrained by down-dip projections of surface geology, `'true amplitude'' processing, rock physical property studies and modelling. 3-D seismic modelling has been used to evaluate the seismic response of two alternative down-dip projection models. To constrain the interpretation in the southern part of the profile, `'true amplitude'' processing has provided information on the strength of the reflections. Density and velocity measurements on core samples collected up-dip from the region of the seismic line have been used to evaluate reflection coefficients of typical lithologic boundaries in the region. The cover-basement contact itself is not a source of strong reflections, but strong reflections arise from within the overlaying metasedimentary cover sequence, allowing the geometry of the top of the basement to be determined on the basis of `'near-basement'' reflections. The front of the external crystalline massifs is shown to extend beneath the Prealps, about 6 km north of the expected position. A 2-D model whose seismic response shows reflection patterns very similar to the observed is proposed.

Implication of the endocannabinoid system in the locomotor hyperactivity associated with congenital hypothyroidism

Alterations in motor functions are well-characterized features observed in humans and experimental animals subjected to thyroid hormone dysfunctions during development. Here we show that congenitally hypothyroid rats display hyperactivity in the adult life. This phenotype was associated with a decreased content of cannabinoid receptor type 1 (CB(1)) mRNA in the striatum and a reduction in the number of binding sites in both striatum and projection areas. These findings suggest that hyperactivity may be the consequence of a thyroid hormone deficiency-induced removal of the endocannabinoid tone, normally acting as a brake for hyperactivity at the basal ganglia. In agreement with the decrease in CB(1) receptor gene expression, a lower cannabinoid response, measured by biochemical, genetic and behavioral parameters, was observed in the hypothyroid animals. Finally, both CB(1) receptor gene expression and the biochemical and behavioral dysfunctions found in the hypothyroid animals were improved after a thyroid hormone replacement treatment. Thus, the present study suggests that impairment in the endocannabinoid system can underlay the hyperactive phenotype associated with hypothyroidism.

Calcium imaging of odor-evoked responses in the Drosophila antennal lobe.

The antennal lobe is the primary olfactory center in the insect brain and represents the anatomical and functional equivalent of the vertebrate olfactory bulb. Olfactory information in the external world is transmitted to the antennal lobe by olfactory sensory neurons (OSNs), which segregate to distinct regions of neuropil called glomeruli according to the specific olfactory receptor they express. Here, OSN axons synapse with both local interneurons (LNs), whose processes can innervate many different glomeruli, and projection neurons (PNs), which convey olfactory information to higher olfactory brain regions. Optical imaging of the activity of OSNs, LNs and PNs in the antennal lobe - traditionally using synthetic calcium indicators (e.g. calcium green, FURA-2) or voltage-sensitive dyes (e.g. RH414) - has long been an important technique to understand how olfactory stimuli are represented as spatial and temporal patterns of glomerular activity in many species of insects. Development of genetically-encoded neural activity reporters, such as the fluorescent calcium indicators G-CaMP and Cameleon, the bioluminescent calcium indicator GFP-aequorin, or a reporter of synaptic transmission, synapto-pHluorin has made the olfactory system of the fruitfly, Drosophila melanogaster, particularly accessible to neurophysiological imaging, complementing its comprehensively-described molecular, electrophysiological and neuroanatomical properties. These reporters can be selectively expressed via binary transcriptional control systems (e.g. GAL4/UAS, LexA/LexAop, Q system) in defined populations of neurons within the olfactory circuitry to dissect with high spatial and temporal resolution how odor-evoked neural activity is represented, modulated and transformed. Here we describe the preparation and analysis methods to measure odor-evoked responses in the Drosophila antennal lobe using G-CaMP. The animal preparation is minimally invasive and can be adapted to imaging using wide-field fluorescence, confocal and two-photon microscopes.

Mate and fuse: how yeast cells do it.

Many cells are able to orient themselves in a non-uniform environment by responding to localized cues. This leads to a polarized cellular response, where the cell can either grow or move towards the cue source. Fungal haploid cells secrete pheromones to signal mating, and respond by growing a mating projection towards a potential mate. Upon contact of the two partner cells, these fuse to form a diploid zygote. In this review, we present our current knowledge on the processes of mating signalling, pheromone-dependent polarized growth and cell fusion in Saccharomyces cerevisiae and Schizosaccharomyces pombe, two highly divergent ascomycete yeast models. While the global architecture of the mating response is very similar between these two species, they differ significantly both in their mating physiologies and in the molecular connections between pheromone perception and downstream responses. The use of both yeast models helps enlighten both conserved solutions and species-specific adaptations to a general biological problem.

Ultrastructural description of Ceratomyxa microlepis sp. nov. (Phylum Myxozoa): a parasite infecting the gall bladder of Hemiodus microlepis, a freshwater teleost from the Amazon River

A new ceratomyxid parasite was examined for taxonomic identification, upon being found infecting the gall bladder of Hemiodus microlepis (Teleostei: Hemiodontidae), a freshwater teleost collected from the Amazon River, Brazil. Light and transmission electron microscopy revealed elongated crescent-shaped spores constituted by two asymmetrical shell valves united along a straight sutural line, each possessing a lateral projection. The spores body measured 5.2 ± 0.4 µm (n = 25) in length and 35.5 ± 0.9 µm (n = 25) in total thickness. The lateral projections were asymmetric, one measuring 18.1 ± 0.5 µm (n = 25) in thickness and the other measuring 17.5 ± 0.5 µm (n = 25) in thickness. Two equal-sized subspherical polar capsules measuring 2.2 ± 0.3 µm in diameter were located at the same level, each possessing a polar filament with 5-6 coils. The sporoplasm was binucleate. Considering the morphometric data analyzed from the microscopic observations, as well as the host species and its geographical location, this paper describes a new myxosporean species, herein named Ceratomyxa microlepis sp. nov.; therefore representing the first description of a freshwater ceratomyxid from the South American region.

Mechanical and elemental characterization of solder joints and welds using a gold-palladium alloy

Purpose. This study was conducted to determine whether newer infrared or laser welding technologies created joints superior to traditional furnace or torch soldering methods of joining metals. It was designed to assess the mechanical resistance, the characteristics of the fractured surfaces, and the elemental diffusion of joints obtained by four different techniques: (1) preceramic soldering with a propane-oxygen torch, (2) postceramic soldering with a porcelain furnace, (3) preceramic and (4) postceramic soldering with an infrared heat source, and (5) laser welding. Material and methods. Mechanical resistance was determined by measuring the ultimate tensile strength of the joint and by determining their resistance to fatigue loading. Elemental diffusion to and from the joint was assessed with microprobe tracings. Scanning electron microscopy micrographs of the fractured surface were also obtained and evaluated. Results. Under monotonic tensile stress, three groups emerged: The laser welds were the strongest, the preceramic joints ranged second, and the postceramic joints were the weakest. Under fatigue stress, the order was as follows: first, the preceramic joints, and second, a group that comprised both postceramic joints and the laser welds. Inspection of the fractographs revealed several fracture modes but no consistent pattern emerged. Microprobe analyses demonstrated minor diffusion processes in the preceramic joints, whereas significant diffusion was observed in the postceramic joints. Clinical Implications. The mechanical resistance data conflicted as to the strength that could be expected of laser welded joints. On the basis of fatigue resistance of the joints, neither infrared solder joints nor laser welds were stronger than torch or furnace soldered joints.

Neurofilament triplet proteins are restricted to a subset of neurons in the rat neocortex.

The cellular localisation of neurofilament triplet subunits was investigated in the rat neocortex. A subset of mainly pyramidal neurons showed colocalisation of subunit immunolabelling throughout the neocortex, including labelling with the antibody SMI32, which has been used extensively in other studies of the primate cortex as a selective cellular marker. Neurofilament-labelled neurons were principally localised to two or three cell layers in most cortical regions, but dramatically reduced labelling was present in areas such as the perirhinal cortex, anterior cingulate and a strip of cortex extending from caudal motor regions through the medial parietal region to secondary visual areas. However, quantitative analysis demonstrated a similar proportion (10-20%) of cells with neurofilament triplet labelling in regions of high or low labelling. Combining retrograde tracing with immunolabelling showed that cellular content of the neurofilament proteins was not correlated with the length of projection. Double labelling immunohistochemistry demonstrated that neurofilament content in axons was closely associated with myelination. Analysis of SMI32 labelling in development indicated that content of this epitope within cell bodies was associated with relatively late maturation, between postnatal days 14 and 21. This study is further evidence of a cell type-specific regulation of neurofilament proteins within neocortical neurons. Neurofilament triplet content may be more closely related to the degree of myelination, rather than the absolute length, of the projecting axon.

Trabecular bone score (TBS) the new parameter of 2D texture analysis for the evaluation of 3D bone micro architecture status

X-ray is a technology that is used for numerous applications in the medical field. The process of X-ray projection gives a 2-dimension (2D) grey-level texture from a 3- dimension (3D) object. Until now no clear demonstration or correlation has positioned the 2D texture analysis as a valid indirect evaluation of the 3D microarchitecture. TBS is a new texture parameter based on the measure of the experimental variogram. TBS evaluates the variation between 2D image grey-levels. The aim of this study was to evaluate existing correlations between 3D bone microarchitecture parameters - evaluated from μCT reconstructions - and the TBS value, calculated on 2D projected images. 30 dried human cadaveric vertebrae were acquired on a micro-scanner (eXplorer Locus, GE) at isotropic resolution of 93 μm. 3D vertebral body models were used. The following 3D microarchitecture parameters were used: Bone volume fraction (BV/TV), Trabecular thickness (TbTh), trabecular space (TbSp), trabecular number (TbN) and connectivity density (ConnD). 3D/2D projections has been done by taking into account the Beer-Lambert Law at X-ray energy of 50, 100, 150 KeV. TBS was assessed on 2D projected images. Correlations between TBS and the 3D microarchitecture parameters were evaluated using a linear regression analysis. Paired T-test is used to assess the X-ray energy effects on TBS. Multiple linear regressions (backward) were used to evaluate relationships between TBS and 3D microarchitecture parameters using a bootstrap process. BV/TV of the sample ranged from 18.5 to 37.6% with an average value at 28.8%. Correlations' analysis showedthat TBSwere strongly correlatedwith ConnD(0.856≤r≤0.862; p<0.001),with TbN (0.805≤r≤0.810; p<0.001) and negatively with TbSp (−0.714≤r≤−0.726; p<0.001), regardless X-ray energy. Results show that lower TBS values are related to "degraded" microarchitecture, with low ConnD, low TbN and a high TbSp. The opposite is also true. X-ray energy has no effect onTBS neither on the correlations betweenTBS and the 3Dmicroarchitecture parameters. In this study, we demonstrated that TBS was significantly correlated with 3D microarchitecture parameters ConnD and TbN, and negatively with TbSp, no matter what X-ray energy has been used. This article is part of a Special Issue entitled ECTS 2011. Disclosure of interest: None declared.

Intracortical connectivity of layer VI pyramidal neurons in the somatosensory cortex of normal and barrelless mice.

In mice, barrels in layer IV of the somatosensory cortex correspond to the columnar representations of whisker follicles. In barrelless (BRL) mice, barrels are absent, but functionally, a columnar organization persists. Previously we characterized the aberrant geometry of thalamic projection of BRL mice using axonal reconstructions of individual neurons. Here we proceeded with the analysis of the intracortical projections from layer VI pyramidal neurons, to assess their contribution to the columnar organization. From series of tangential sections we reconstructed the axon collaterals of individual layer VI pyramidal neurons in the C2 barrel column that were labelled with biocytin [controls from normal (NOR) strain, 19 cells; BRL strain, nine cells]. Using six morphological parameters in a cluster analysis, we showed that layer VI neurons in NOR mice are distributed into four clusters distinguished by the radial and tangential extent of their intracortical projections. These clusters correlated with the cortical or subcortical projection of the main axon. In BRL mice, neurons were distributed within the same four clusters, but their projections to the granular and supragranular layers were significantly smaller and their tangential projection was less columnar than in NOR mice. However, in both strains the intracortical projections had a preference for the appropriate barrel column (C2), indicating that layer VI pyramidal cells could participate in the functional columnar organization of the barrel cortex. Correlative light and electron microscopy analyses provided morphometric data on the intracortical synaptic boutons and synapses of layer VI pyramidal neurons and revealed that projections to layer IV preferentially target excitatory dendritic spines and shafts.

Low-dose multidetector computed tomography of the cervical spine: optimization of iterative reconstruction strength levels.

BACKGROUND: Iterative reconstruction (IR) techniques reduce image noise in multidetector computed tomography (MDCT) imaging. They can therefore be used to reduce radiation dose while maintaining diagnostic image quality nearly constant. However, CT manufacturers offer several strength levels of IR to choose from. PURPOSE: To determine the optimal strength level of IR in low-dose MDCT of the cervical spine. MATERIAL AND METHODS: Thirty consecutive patients investigated by low-dose cervical spine MDCT were prospectively studied. Raw data were reconstructed using filtered back-projection and sinogram-affirmed IR (SAFIRE, strength levels 1 to 5) techniques. Image noise, signal-to-noise ratio (SNR), and contrast-to-noise ratio (CNR) were measured at C3-C4 and C6-C7 levels. Two radiologists independently and blindly evaluated various anatomical structures (both dense and soft tissues) using a 4-point scale. They also rated the overall diagnostic image quality using a 10-point scale. RESULTS: As IR strength levels increased, image noise decreased linearly, while SNR and CNR both increased linearly at C3-C4 and C6-C7 levels (P < 0.001). For the intervertebral discs, the content of neural foramina and dural sac, and for the ligaments, subjective image quality scores increased linearly with increasing IR strength level (P ≤ 0.03). Conversely, for the soft tissues and trabecular bone, the scores decreased linearly with increasing IR strength level (P < 0.001). Finally, the overall diagnostic image quality scores increased linearly with increasing IR strength level (P < 0.001). CONCLUSION: The optimal strength level of IR in low-dose cervical spine MDCT depends on the anatomical structure to be analyzed. For the intervertebral discs and the content of neural foramina, high strength levels of IR are recommended.

Classification of sags according to their origin based on the waveform similarity

A statistical method for classification of sags their origin downstream or upstream from the recording point is proposed in this work. The goal is to obtain a statistical model using the sag waveforms useful to characterise one type of sags and to discriminate them from the other type. This model is built on the basis of multi-way principal component analysis an later used to project the available registers in a new space with lower dimension. Thus, a case base of diagnosed sags is built in the projection space. Finally classification is done by comparing new sags against the existing in the case base. Similarity is defined in the projection space using a combination of distances to recover the nearest neighbours to the new sag. Finally the method assigns the origin of the new sag according to the origin of their neighbours

The spatiotemporal trajectory of a dengue epidemic in a medium-sized city

Understanding the transmission dynamics of infectious diseases is important to allow for improvements of control measures. To investigate the spatiotemporal pattern of an epidemic dengue occurred at a medium-sized city in the Northeast Region of Brazil in 2009, we conducted an ecological study of the notified dengue cases georeferenced according to epidemiological week (EW) and home address. Kernel density estimation and space-time interaction were analysed using the Knox method. The evolution of the epidemic was analysed using an animated projection technique. The dengue incidence was 6.918.7/100,000 inhabitants; the peak of the epidemic occurred from 8 February-1 March, EWs 6-9 (828.7/100,000 inhabitants). There were cases throughout the city and was identified space-time interaction. Three epicenters were responsible for spreading the disease in an expansion and relocation diffusion pattern. If the health services could detect in real time the epicenters and apply nimbly control measures, may possibly reduce the magnitude of dengue epidemics.

Coronary MR imaging using free-breathing 3D steady-state free precession with radial k-space sampling.

PURPOSE: To investigate the potential of free-breathing 3D steady-state free precession (SSFP) imaging with radial k-space sampling for coronary MR-angiography (MRA), coronary projection MR-angiography and coronary vessel wall imaging. MATERIALS AND METHODS: A navigator-gated free-breathing T2-prepared 3D SSFP sequence (TR = 6.1 ms, TE = 3.0 ms, flip angle = 120 degrees, field-of-view = 360 mm(2)) with radial k-space sampling (384 radials) was implemented for coronary MRA. For projection coronary MRA, this sequence was combined with a 2D selective aortic spin tagging pulse. Coronary vessel wall imaging was performed using a high-resolution inversion-recovery black-blood 3D radial SSFP sequence (384 radials, TR = 5.3 ms, TE = 2.7 ms, flip angle = 55 degrees, reconstructed resolution 0.35 x 0.35 x 1.2 mm(3)) and a local re-inversion pulse. Six healthy volunteers (two for each sequence) were investigated. Motion artifact level was assessed by two radiologists. Results: In coronary MRA, the coronary lumen was displayed with a high signal and high contrast to the surrounding lumen. Projection coronary MRA demonstrated selective visualization of the coronary lumen while surrounding tissue was almost completely suppressed. In coronary vessel wall imaging, the vessel wall was displayed with a high signal when compared to the blood pool and the surrounding tissue. No visible motion artifacts were seen. Conclusion: 3D radial SSFP imaging enables coronary MRA, coronary projection MRA and coronary vessel wall imaging with a low motion artifact level.