A systematic method of parameterisation estimation using data assimilation

In numerical weather prediction, parameterisations are used to simulate missing physics in the model. These can be due to a lack of scientific understanding or a lack of computing power available to address all the known physical processes. Parameterisations are sources of large uncertainty in a model as parameter values used in these parameterisations cannot be measured directly and hence are often not well known; and the parameterisations themselves are also approximations of the processes present in the true atmosphere. Whilst there are many efficient and effective methods for combined state/parameter estimation in data assimilation (DA), such as state augmentation, these are not effective at estimating the structure of parameterisations. A new method of parameterisation estimation is proposed that uses sequential DA methods to estimate errors in the numerical models at each space-time point for each model equation. These errors are then fitted to pre-determined functional forms of missing physics or parameterisations that are based upon prior information. We applied the method to a one-dimensional advection model with additive model error, and it is shown that the method can accurately estimate parameterisations, with consistent error estimates. Furthermore, it is shown how the method depends on the quality of the DA results. The results indicate that this new method is a powerful tool in systematic model improvement.

Integration of steroid hormone initiated membrane action to genomic function in the brain

Estrogen is a ligand for the estrogen receptor (ER), which on binding 17beta-estradiol, functions as a ligand-activated transcription factor and regulates the transcription of target genes. This is the slow genomic mode of action. However, rapid non-genomic actions of estrogen also exist at the cell membrane. Using a novel two-pulse paradigm in which the first pulse rapidly initiates non-genomic actions using a membrane-limited estrogen conjugate (E-BSA), while the second pulse promotes genomic transcription from a consensus estrogen response element (ERE), we have demonstrated that rapid actions of estrogen potentiate the slower transcriptional response from an ERE-reporter in neuroblastoma cells. Since rapid actions of estrogen activate kinases, we used selective inhibitors in the two-pulse paradigm to determine the intracellular signaling cascades important in such potentiation. Inhibition of protein kinase A (PKA), PKC, mitogen activated protein kinase (MAPK) or phosphatidylinositol 3-OH kinase (PI-3K) in the first pulse decreases potentiation of transcription. Also, our data with both dominant negative and constitutive mutants of Galpha subunits show that Galpha(q) initiates the rapid signaling cascade at the membrane in SK-N-BE(2)C neuroblastoma cells. We discuss two models of multiple kinase activation at the membrane Pulses of estrogen induce lordosis behavior in female rats. Infusion of E-BSA into the ventromedial hypothalamus followed by 17beta-estradiol in the second pulse could induce lordosis behavior, demonstrating the applicability of this paradigm in vivo. A model where non-genomic actions of estrogen couple to genomic actions unites both aspects of hormone action.

Discovering biomarkers for antidepressant response: protocol from the Canadian biomarker integration network in depression (CAN-BIND) and clinical characteristics of the first patient cohort

Background Major Depressive Disorder (MDD) is among the most prevalent and disabling medical conditions worldwide. Identification of clinical and biological markers (“biomarkers”) of treatment response could personalize clinical decisions and lead to better outcomes. This paper describes the aims, design, and methods of a discovery study of biomarkers in antidepressant treatment response, conducted by the Canadian Biomarker Integration Network in Depression (CAN-BIND). The CAN-BIND research program investigates and identifies biomarkers that help to predict outcomes in patients with MDD treated with antidepressant medication. The primary objective of this initial study (known as CAN-BIND-1) is to identify individual and integrated neuroimaging, electrophysiological, molecular, and clinical predictors of response to sequential antidepressant monotherapy and adjunctive therapy in MDD. Methods CAN-BIND-1 is a multisite initiative involving 6 academic health centres working collaboratively with other universities and research centres. In the 16-week protocol, patients with MDD are treated with a first-line antidepressant (escitalopram 10–20 mg/d) that, if clinically warranted after eight weeks, is augmented with an evidence-based, add-on medication (aripiprazole 2–10 mg/d). Comprehensive datasets are obtained using clinical rating scales; behavioural, dimensional, and functioning/quality of life measures; neurocognitive testing; genomic, genetic, and proteomic profiling from blood samples; combined structural and functional magnetic resonance imaging; and electroencephalography. De-identified data from all sites are aggregated within a secure neuroinformatics platform for data integration, management, storage, and analyses. Statistical analyses will include multivariate and machine-learning techniques to identify predictors, moderators, and mediators of treatment response. Discussion From June 2013 to February 2015, a cohort of 134 participants (85 outpatients with MDD and 49 healthy participants) has been evaluated at baseline. The clinical characteristics of this cohort are similar to other studies of MDD. Recruitment at all sites is ongoing to a target sample of 290 participants. CAN-BIND will identify biomarkers of treatment response in MDD through extensive clinical, molecular, and imaging assessments, in order to improve treatment practice and clinical outcomes. It will also create an innovative, robust platform and database for future research.

Use of the DNA Checkerboard hybridization method for detection and quantitation of Candida species in oral microbiota

The DNA Checkerboard method enables the simultaneous identification of distinct microorganisms in a large number of samples and employs up to 45 whole genomic DNA probes to gram-negative and gram-positive bacterial species present in subgingival biofilms. Collectively, they account for 55%-60% of the bacteria in subgingival biofilms. In this study, we present the DNA Checkerboard hybridization as an alternative method for the detection and quantitation of Candida species in oral cavities. Our results reveal that DNA Checkerboard is sensitive enough and constitutes a powerful and appropriate method for detecting and quantifying Candida species found in the oral cavity.

Validation of an ELISA method for the serological diagnosis of canine brucellosis due to Brucella canis

In the present study, the validation of an enzyme-linked immunosorbent assay (ELISA) for serodiagnosis of canine brucellosis is described. Two different antigenic extracts, obtained by heat or ultrasonic homogenization of microbial antigens from a wild isolate of Brucella canis bacteria, were compared by ELISA and Western blot (WB). A total of 145 canine sera were used to define sensitivity, specificity and accuracy of the ELISA as follows: (1) sera from 34 animals with natural B. canis infection, confirmed by blood culture and PCR, as well as 51 sera samples from healthy dogs with negative results by the agar-gel immunodiffusion (ACID) test for canine brucellosis, were used as the control panel for B. cants infection; and (2) to scrutinize the possibility of cross reactions with other common dog infections in the same geographical area in Brazil, 60 sera samples from dogs harboring known infections by Leptospira sp., Ehrlichia canis, canine distemper virus (CDV), Neospora caninum, Babesia canis and Leishmania chagasi (10 in each group) were included in the study. The ELISA using heat soluble bacterial extract (HE-antigen) as antigen showed the best values of sensitivity (91.18%), specificity (100%) and accuracy (96.47%). In the WB analyses, the HE-antigen showed no cross-reactivity with sera from dogs with different infections, while the B. canis sonicate had various protein bands identified by those sera. The performance of the ELISA standardized with the heat soluble B. canis antigen indicates that this assay can be used as a reliable and practical method to confirm infection by this microorganism, as well as a tool for seroepidemiological studies. (C) 2010 Elsevier Ltd. All rights reserved.

A METHOD FOR THE STUDY OF ACCRETION DISK EMISSION IN CATACLYSMIC VARIABLES. I. THE MODEL

We have developed a spectrum synthesis method for modeling the ultraviolet (UV) emission from the accretion disk from cataclysmic variables (CVs). The disk is separated into concentric rings, with an internal structure from the Wade & Hubeny disk-atmosphere models. For each ring, a wind atmosphere is calculated in the comoving frame with a vertical velocity structure obtained from a solution of the Euler equation. Using simple assumptions, regarding rotation and the wind streamlines, these one-dimensional models are combined into a single 2.5-dimensional model for which we compute synthetic spectra. We find that the resulting line and continuum behavior as a function of the orbital inclination is consistent with the observations, and verify that the accretion rate affects the wind temperature, leading to corresponding trends in the intensity of UV lines. In general, we also find that the primary mass has a strong effect on the P Cygni absorption profiles, the synthetic emission line profiles are strongly sensitive to the wind temperature structure, and an increase in the mass-loss rate enhances the resonance line intensities. Synthetic spectra were compared with UV data for two high orbital inclination nova-like CVs-RW Tri and V347 Pup. We needed to include disk regions with arbitrary enhanced mass loss to reproduce reasonably well widths and line profiles. This fact and a lack of flux in some high ionization lines may be the signature of the presence of density-enhanced regions in the wind, or alternatively, may result from inadequacies in some of our simplifying assumptions.

SEM and Neurohistological Observations of Nerve Endings in the Middle Region of the Tongue of the Collared Peccary (Tayassu tajacu): A Silver Impregnation Method

The presence of lingual papillae and the nerve endings in the middle region of the tongue mucosa of collared peccary (Tayassu tajacu) were studied using scanning electron microscopy and light microscopy, based upon the silver impregnation method. The middle region of tongue mucosa revealed numerous filiform and fungiform papillae. The thick epithelial layer showed epithelial cells and a dense connective tissue layer containing nerve fibre bundles and capillaries. The sensory nerve endings, intensely stained by silver impregnation, were usually non-encapsulated and extended into the connective tissue of the filiform and fungiform papillae very close to the epithelial cells. In some regions, the sensory nerves fibres formed a dense and complex network of fine fibrils. The presence of these nerve fibrils may characterize the mechanisms of transmission of sensitive impulses to the tongue mucosa.

The geodiamolide H, derived from Brazilian sponge Geodia corticostylifera, regulates actin cytoskeleton, migration and invasion of breast cancer cells cultured in three-dimensional environment

We are investigating effects of the depsipeptide geodiamolide H, isolated from the Brazilian sponge Geodia corticostylifera, on cancer cell lines grown in 3D environment. As shown previously geodiamolide H disrupts actin cytoskeleton in both sea urchin eggs and breast cancer cell monolayers. We used a normal mammary epithelial cell line MCF 10A that in 3D assay results formation of polarized spheroids. We also used cell lines derived from breast tumors with different degrees of differentiation: MCF7 positive for estrogen receptor and the Hs578T, negative for hormone receptors. Cells were placed on top of Matrigel. Spheroids obtained from these cultures were treated with geodiamolide H. Control and treated samples were analyzed by light and confocal microscopy. Geodiamolide H dramatically affected the poorly differentiated and aggressive Hs578T cell line. The peptide reverted HsS78T malignant phenotype to polarized spheroid-like structures. MCF7 cells treated by geodiamolide H exhibited polarization compared to controls. Geodiamolide H induced striking phenotypic modifications in Hs578T cell line and disruption of actin cytoskeleton. We investigated effects of geodiamolide H on migration and invasion of Hs578T cells. Time-lapse microscopy showed that the peptide inhibited migration of these cells in a dose-dependent manner. Furthermore invasion assays revealed that geodiamolide H induced a 30% decrease on invasive behavior of Hs578T cells. Our results suggest that geodiamolide H inhibits migration and invasion of Hs578T cells probably through modifications in actin cytoskeleton. The fact that normal cell lines were not affected by treatment with geodiamolide H stimulates new studies towards therapeutic use for this peptide.

Numerical prediction of three-dimensional time-dependent viscoelastic extrudate swell using differential and algebraic models

This study investigates the numerical simulation of three-dimensional time-dependent viscoelastic free surface flows using the Upper-Convected Maxwell (UCM) constitutive equation and an algebraic explicit model. This investigation was carried out to develop a simplified approach that can be applied to the extrudate swell problem. The relevant physics of this flow phenomenon is discussed in the paper and an algebraic model to predict the extrudate swell problem is presented. It is based on an explicit algebraic representation of the non-Newtonian extra-stress through a kinematic tensor formed with the scaled dyadic product of the velocity field. The elasticity of the fluid is governed by a single transport equation for a scalar quantity which has dimension of strain rate. Mass and momentum conservations, and the constitutive equation (UCM and algebraic model) were solved by a three-dimensional time-dependent finite difference method. The free surface of the fluid was modeled using a marker-and-cell approach. The algebraic model was validated by comparing the numerical predictions with analytic solutions for pipe flow. In comparison with the classical UCM model, one advantage of this approach is that computational workload is substantially reduced: the UCM model employs six differential equations while the algebraic model uses only one. The results showed stable flows with very large extrudate growths beyond those usually obtained with standard differential viscoelastic models. (C) 2010 Elsevier Ltd. All rights reserved.

Numerical study of the atmospheric flow over a coastal cliff

This work presents a numerical method suitable for the study of the development of internal boundary layers (IBL) and their characteristics for flows over various types of coastal cliffs. The IBL is an important meteorological occurrence for flows with surface roughness and topographical step changes. A two-dimensional flow program was used for this study. The governing equations were written using the vorticity-velocity formulation. The spatial derivatives were discretized by high-order compact finite differences schemes. The time integration was performed with a low storage fourth-order Runge-Kutta scheme. The coastal cliff (step) was specified through an immersed boundary method. The validation of the code was done by comparison of the results with experimental and observational data. The numerical simulations were carried out for different coastal cliff heights and inclinations. The results show that the predominant factors for the height of the IBL and its characteristics are the upstream velocity, and the height and form (inclination) of the coastal cliff. Copyright (C) 2010 John Wiley & Sons, Ltd.

A discontinuous-Galerkin-based immersed boundary method with non-homogeneous boundary conditions and its application to elasticity

We propose a discontinuous-Galerkin-based immersed boundary method for elasticity problems. The resulting numerical scheme does not require boundary fitting meshes and avoids boundary locking by switching the elements intersected by the boundary to a discontinuous Galerkin approximation. Special emphasis is placed on the construction of a method that retains an optimal convergence rate in the presence of non-homogeneous essential and natural boundary conditions. The role of each one of the approximations introduced is illustrated by analyzing an analog problem in one spatial dimension. Finally, extensive two- and three-dimensional numerical experiments on linear and nonlinear elasticity problems verify that the proposed method leads to optimal convergence rates under combinations of essential and natural boundary conditions. (C) 2009 Elsevier B.V. All rights reserved.

Numerical solution of the PTT constitutive equation for unsteady three-dimensional free surface flows

This work deals with the development of a numerical technique for simulating three-dimensional viscoelastic free surface flows using the PTT (Phan-Thien-Tanner) nonlinear constitutive equation. In particular, we are interested in flows possessing moving free surfaces. The equations describing the numerical technique are solved by the finite difference method on a staggered grid. The fluid is modelled by a Marker-and-Cell type method and an accurate representation of the fluid surface is employed. The full free surface stress conditions are considered. The PTT equation is solved by a high order method, which requires the calculation of the extra-stress tensor on the mesh contours. To validate the numerical technique developed in this work flow predictions for fully developed pipe flow are compared with an analytic solution from the literature. Then, results of complex free surface flows using the FIT equation such as the transient extrudate swell problem and a jet flowing onto a rigid plate are presented. An investigation of the effects of the parameters epsilon and xi on the extrudate swell and jet buckling problems is reported. (C) 2010 Elsevier B.V. All rights reserved.

Numerical Solution of the Upper-Convected Maxwell Model for Three-Dimensional Free Surface Flows

This work presents a finite difference technique for simulating three-dimensional free surface flows governed by the Upper-Convected Maxwell (UCM) constitutive equation. A Marker-and-Cell approach is employed to represent the fluid free surface and formulations for calculating the non-Newtonian stress tensor on solid boundaries are developed. The complete free surface stress conditions are employed. The momentum equation is solved by an implicit technique while the UCM constitutive equation is integrated by the explicit Euler method. The resulting equations are solved by the finite difference method on a 3D-staggered grid. By using an exact solution for fully developed flow inside a pipe, validation and convergence results are provided. Numerical results include the simulation of the transient extrudate swell and the comparison between jet buckling of UCM and Newtonian fluids.

An automatic correction method for the heel effect in digitized mammography images

The most significant radiation field nonuniformity is the well-known Heel effect. This nonuniform beam effect has a negative influence on the results of computer-aided diagnosis of mammograms, which is frequently used for early cancer detection. This paper presents a method to correct all pixels in the mammography image according to the excess or lack on radiation to which these have been submitted as a result of the this effect. The current simulation method calculates the intensities at all points of the image plane. In the simulated image, the percentage of radiation received by all the points takes the center of the field as reference. In the digitized mammography, the percentages of the optical density of all the pixels of the analyzed image are also calculated. The Heel effect causes a Gaussian distribution around the anode-cathode axis and a logarithmic distribution parallel to this axis. Those characteristic distributions are used to determine the center of the radiation field as well as the cathode-anode axis, allowing for the automatic determination of the correlation between these two sets of data. The measurements obtained with our proposed method differs on average by 2.49 mm in the direction perpendicular to the anode-cathode axis and 2.02 mm parallel to the anode-cathode axis of commercial equipment. The method eliminates around 94% of the Heel effect in the radiological image and the objects will reflect their x-ray absorption. To evaluate this method, experimental data was taken from known objects, but could also be done with clinical and digital images.

Aerosol inorganic composition at a tropical site: Discrepancies between filter-based sampling and a semi-continuous method

The concentrations of the water-soluble inorganic aerosol species, ammonium (NH4+), nitrate (NO3-), chloride (Cl-), and sulfate (SO42-), were measured from September to November 2002 at a pasture site in the Amazon Basin (Rondnia, Brazil) (LBA-SMOCC). Measurements were conducted using a semi-continuous technique (Wet-annular denuder/Steam-Jet Aerosol Collector: WAD/SJAC) and three integrating filter-based methods, namely (1) a denuder-filter pack (DFP: Teflon and impregnated Whatman filters), (2) a stacked-filter unit (SFU: polycarbonate filters), and (3) a High Volume dichotomous sampler (HiVol: quartz fiber filters). Measurements covered the late dry season (biomass burning), a transition period, and the onset of the wet season (clean conditions). Analyses of the particles collected on filters were performed using ion chromatography (IC) and Particle-Induced X-ray Emission spectrometry (PIXE). Season-dependent discrepancies were observed between the WAD/SJAC system and the filter-based samplers. During the dry season, when PM2.5 (D-p <= 2.5 mu m) concentrations were similar to 100 mu g m(-3), aerosol NH4+ and SO42- measured by the filter-based samplers were on average two times higher than those determined by the WAD/SJAC. Concentrations of aerosol NO3- and Cl- measured with the HiVol during daytime, and with the DFP during day- and nighttime also exceeded those of the WAD/SJAC by a factor of two. In contrast, aerosol NO3- and Cl- measured with the SFU during the dry season were nearly two times lower than those measured by the WAD/SJAC. These differences declined markedly during the transition period and towards the cleaner conditions during the onset of the wet season (PM2.5 similar to 5 mu g m(-3)); when filter-based samplers measured on average 40-90% less than the WAD/SJAC. The differences were not due to consistent systematic biases of the analytical techniques, but were apparently a result of prevailing environmental conditions and different sampling procedures. For the transition period and wet season, the significance of our results is reduced by a low number of data points. We argue that the observed differences are mainly attributable to (a) positive and negative filter sampling artifacts, (b) presence of organic compounds and organosulfates on filter substrates, and (c) a SJAC sampling efficiency of less than 100%.