Stabilized Petrov-Galerkin methods for the convection-diffusion-reaction and the Helmholtz equations

We present two new stabilized high-resolution numerical methods for the convection–diffusion–reaction (CDR) and the Helmholtz equations respectively. The work embarks upon a priori analysis of some consistency recovery procedures for some stabilization methods belonging to the Petrov–Galerkin framework. It was found that the use of some standard practices (e.g. M-Matrices theory) for the design of essentially non-oscillatory numerical methods is not feasible when consistency recovery methods are employed. Hence, with respect to convective stabilization, such recovery methods are not preferred. Next, we present the design of a high-resolution Petrov–Galerkin (HRPG) method for the 1D CDR problem. The problem is studied from a fresh point of view, including practical implications on the formulation of the maximum principle, M-Matrices theory, monotonicity and total variation diminishing (TVD) finite volume schemes. The current method is next in line to earlier methods that may be viewed as an upwinding plus a discontinuity-capturing operator. Finally, some remarks are made on the extension of the HRPG method to multidimensions. Next, we present a new numerical scheme for the Helmholtz equation resulting in quasi-exact solutions. The focus is on the approximation of the solution to the Helmholtz equation in the interior of the domain using compact stencils. Piecewise linear/bilinear polynomial interpolation are considered on a structured mesh/grid. The only a priori requirement is to provide a mesh/grid resolution of at least eight elements per wavelength. No stabilization parameters are involved in the definition of the scheme. The scheme consists of taking the average of the equation stencils obtained by the standard Galerkin finite element method and the classical finite difference method. Dispersion analysis in 1D and 2D illustrate the quasi-exact properties of this scheme. Finally, some remarks are made on the extension of the scheme to unstructured meshes by designing a method within the Petrov–Galerkin framework.

Human urothelial cells grown on collagen adsorbed to surface-modified polymers.

OBJECTIVES: Tissue engineering methods can be applied to regenerate diseased, or congenitally missing, urinary tract tissues. Urinary tract tissue cell cultures must be established in vitro and adequate matrices, acting as cell carriers, must be developed. Although degradable and nondegradable polymer matrices offer adequate mechanical stability, they are not optimal for cell adherence and growth. To overcome this problem, extracellular matrix proteins, permitting cell adhesion and regulation of cell proliferation and differentiation, can be adsorbed to the surface-modified polymer. METHODS: In this study, nondegradable polymer films, poly(ethylene terephthalate), were used as an experimental model. Films were modified by graft polymerization of acrylic acid to subsequently allow collagen type I and III immobilization. The following adhesion, proliferation of human urothelial cells, and induction of their stratification were analyzed. RESULTS: Collagen adsorption on 0.2 microg/cm2 poly(acrylic acid)-grafted polymer films rendered the matrix apt for human urothelial cell adhesion and proliferation. Furthermore, stratification of urothelial cells was demonstrated on these surface-modified matrices. CONCLUSIONS: These results have shown that surface-modified polymer matrices can be used to act as cell carriers for cultured human urothelial cells. Such a cell-matrix construct could be applied in reparative surgery of the urinary tract.

The eukaryotic promoter database (EPD).

The Eukaryotic Promoter Database (EPD) is an annotated non-redundant collection of eukaryotic POL II promoters for which the transcription start site has been determined experimentally. Access to promoter sequences is provided by pointers to positions in nucleotide sequence entries. The annotation part of an entry includes a description of the initiation site mapping data, exhaustive cross-references to the EMBL nucleotide sequence database, SWISS-PROT, TRANSFAC and other databases, as well as bibliographic references. EPD is structured in a way that facilitates dynamic extraction of biologically meaningful promoter subsets for comparative sequence analysis. WWW-based interfaces have been developed that enable the user to view EPD entries in different formats, to select and extract promoter sequences according to a variety of criteria, and to navigate to related databases exploiting different cross-references. The EPD web site also features yearly updated base frequency matrices for major eukaryotic promoter elements. EPD can be accessed at http://www.epd.isb-sib.ch

Les TIC a l'ensenyament de la modelització matemàtica en Ciències Econòmiques i en Ciències Experimentals

Aquest projecte proposa materials didàctics per a un nou plantejament de les assignatures de Matemàtiques dels primers cursos de Ciències Empresarials i d'Enginyeria Tècnica, més acord amb el procés de convergència europea, basat en la realització de projectes que anomenem “Tallers de Modelització Matemàtica” (TMM) en els quals: (1) Els alumnes parteixen de situacions i problemes reals per als quals han de construir per sí mateixos els models matemàtics més adients i, a partir de la manipulació adequada d’aquests models, poden obtenir la informació necessària per donar-los resposta. (2) El treball de construcció, experimentació i avaluació dels models es realitza amb el suport de la calculadora simbòlica Wiris i del full de càlcul Excel com a instruments “normalitzats” del treball matemàtic d’estudiants i professors. (3) S’adapten els programes de les assignatures de matemàtiques de primer curs per tal de poder-les associar a un petit nombre de Tallers que parteixen de situacions adaptades a cada titulació. L’assignatura de Matemàtiques per a les Ciències Empresarials s’articula entorn de dos tallers independents: “Matrius de transició” pel que fa a l’àlgebra lineal i “Previsió de vendes” per a la modelització funcional en una variable. L’assignatura de Matemàtiques per a l’Enginyeria s’articula entorn d’un únic taller, “Models de poblacions”, que abasta la majoria de continguts del curs: successions i models funcionals en una variable, àlgebra lineal i equacions diferencials. Un conjunt d’exercicis interactius basats en la calculadora simbòlica WIRIS (Wiris-player) serveix de suport per al treball tècnic imprescindible per al desenvolupament de les dues assignatures. L’experimentació d’aquests tallers durant 2 cursos consecutius (2006/07 i 2007/08) en dues universitats catalanes (URL i UAB) ha posat en evidència tant els innegables avantatges del nou dispositiu docent per a l’aprenentatge dels estudiants, així com les restriccions institucionals que actualment dificulten la seva gestió i difusió.

Naissances vaudoises 1993-1994: caractéristiques et facteurs de risque pour une affection chronique: rapport EDEN no 1

Une cohorte de 6477 nouveau-nés de mères résidant dans le Canton du Vaud a été recrutée pendant une année (1993-1994) dans les 18 maternités vaudoises et celle de Châtel-St-Denis. L'objectif de l'étude EDEN (Etude du DEveloppement des Nouveau-nés) est de calculer l'incidence et la prévalence des affections chroniques de toute étiologie et pour toutes les catégories de poids de naissance, à 18 mois et à 4 ans. Ce rapport présente la méthode de l'étude et l'état de santé à la naissance. Cinq critères de sélection non exclusifs ont permis de cibler un groupe de nouveau-nés à haut risque de développer une affection chronique (12% des nouveau-nés, n=760): (1) le petit poids de naissance (n=408, 6.5% des naissances vivantes); (2) une malformation congénitale ou une maladie génétique (n=157, 2.4% des naissances vivantes); (3) une affection susceptible de devenir chronique liée à une utilisation importante des services de soins au cours de la petite enfance (n=61, 0.9% des naissances vivantes); (4) le transfert aux soins intensifs (n=287, 4.4% des naissances vivantes); (5) des difficultés sociales importantes (n=105, 1.6% des naissances vivantes). Le taux d'acceptation de l'étude par les parents est bon (90%). En tout 5.9% des enfants étaient prématurés et 2.2 pour mille sont décédés à < ou = à 7 jours de vie. Selon les indicateurs à disposition, le réseau vaudois répond efficacement aux besoins en soins obstétricaux et néonatals La durée moyenne du séjour hospitalier était de 7 jours, avec des variations importantes. L'influence néfaste du tabagisme pendant la grossesse se manifeste par un doublement du risque de poids de naissance < ou = à 2500g chez les fumeuses; 24% des femmes ont fumé pendant leur grossesse, pour les trois quarts jusqu'à l'accouchement. Un grand potentiel de prévention subsiste dans ce domaine. L'examen des enfants à 18 mois, terminé fin mai 1996, ainsi que celui des 4 ans, permettront de valider les critères de sélection à la naissance comme indicateurs précoces de problèmes de santé chroniques dans la petite enfance. Les nouveaux cas d'affection chronique seront alors signalés par les pédiatres et les médecins spécialistes. [Auteurs, p. 9]

Linear Aggregation In The Social Accounting Matrix Framework

In economic literature, information deficiencies and computational complexities have traditionally been solved through the aggregation of agents and institutions. In inputoutput modelling, researchers have been interested in the aggregation problem since the beginning of 1950s. Extending the conventional input-output aggregation approach to the social accounting matrix (SAM) models may help to identify the effects caused by the information problems and data deficiencies that usually appear in the SAM framework. This paper develops the theory of aggregation and applies it to the social accounting matrix model of multipliers. First, we define the concept of linear aggregation in a SAM database context. Second, we define the aggregated partitioned matrices of multipliers which are characteristic of the SAM approach. Third, we extend the analysis to other related concepts, such as aggregation bias and consistency in aggregation. Finally, we provide an illustrative example that shows the effects of aggregating a social accounting matrix model.

Simulating the Impact of Sectorial Productivity Gains on Two Regional Economies: Key Sectors from a Supply Side Perspective

In this paper we simulate and analyse the economic impact that sectorial productivity gains have on two regional Spanish economies (Catalonia and Extremadura). In particular we study the quantitative effect that each sector’s productivity gain has on household welfare (real disposable income and equivalent variation), on the consumption price indices and factor relative prices, on real production (GDP) and on the government’s net income (net taxation revenues of social transfers to households). The analytical approach consists of a computable general equilibrium model, in which we assume perfect competition and cleared markets, including factor markets. All the parameters and exogenous variables of the model are calibrated by means of two social accounting matrices, one for each region under study. The results allow us to identify those sectors with the greatest impact on consumer welfare as the key sectors in the regional economies. Keywords: Productivity gains, key sectors, computable general equilibrium

Analyse des cannabinoïdes par spectrométrie de masse en mode tandem = Analysis of cannabinoids by tandem mass spectrometry

Depuis quelques années, la spectrométrie de masse en tandem (MS/MS) ne cesse de gagner du terrain comme méthode d'analyse en toxicologie forensique, notamment pour le dosage des cannabinoïdes. Couplée à la chromatographie liquide (LC) ou gazeuse (GC), elle permet l'identification fiable et le dosage rapide du THC, de son précurseur acide, et de ses principaux métabolites, y compris les glucuronides. Au cours de ces dix dernières années, un nombre significatif de publications sont parues sur ce sujet. L'objectif de cet article est de passer en revue les analyses par spectrométrie de masse en tandem des cannabinoïdes dans diverses matrices biologiques. In recent years, tandem mass spectrometry (MS/MS) is gaining ground as a reference method of analysis in clinical and forensic toxicology, especially for the determination of cannabinoids. Coupled to liquid chromatography (LC) or gas chromatography (GC), it allows the definitive identification and rapid determination of THC, its acid precursor, and its major metabolites, including the glucuronides. During the past decade, several methods of analysis of cannabinoids in different matrices have appeared on this subject. The aim of this paper is to review the analysis of cannabinoids by tandem mass spectrometry methods in various biological matrices

Identifying genetic signatures of selection in a non-model species, alpine gentian (Gentiana nivalis L.), using a landscape genetic approach

It is generally accepted that most plant populations are locally adapted. Yet, understanding how environmental forces give rise to adaptive genetic variation is a challenge in conservation genetics and crucial to the preservation of species under rapidly changing climatic conditions. Environmental variation, phylogeographic history, and population demographic processes all contribute to spatially structured genetic variation, however few current models attempt to separate these confounding effects. To illustrate the benefits of using a spatially-explicit model for identifying potentially adaptive loci, we compared outlier locus detection methods with a recently-developed landscape genetic approach. We analyzed 157 loci from samples of the alpine herb Gentiana nivalis collected across the European Alps. Principle coordinates of neighbor matrices (PCNM), eigenvectors that quantify multi-scale spatial variation present in a data set, were incorporated into a landscape genetic approach relating AFLP frequencies with 23 environmental variables. Four major findings emerged. 1) Fifteen loci were significantly correlated with at least one predictor variable (R (adj) (2) > 0.5). 2) Models including PCNM variables identified eight more potentially adaptive loci than models run without spatial variables. 3) When compared to outlier detection methods, the landscape genetic approach detected four of the same loci plus 11 additional loci. 4) Temperature, precipitation, and solar radiation were the three major environmental factors driving potentially adaptive genetic variation in G. nivalis. Techniques presented in this paper offer an efficient method for identifying potentially adaptive genetic variation and associated environmental forces of selection, providing an important step forward for the conservation of non-model species under global change.

Prenatally engineered autologous amniotic fluid stem cell-based heart valves in the fetal circulation.

Prenatal heart valve interventions aiming at the early and systematic correction of congenital cardiac malformations represent a promising treatment option in maternal-fetal care. However, definite fetal valve replacements require growing implants adaptive to fetal and postnatal development. The presented study investigates the fetal implantation of prenatally engineered living autologous cell-based heart valves. Autologous amniotic fluid cells (AFCs) were isolated from pregnant sheep between 122 and 128 days of gestation via transuterine sonographic sampling. Stented trileaflet heart valves were fabricated from biodegradable PGA-P4HB composite matrices (n = 9) and seeded with AFCs in vitro. Within the same intervention, tissue engineered heart valves (TEHVs) and unseeded controls were implanted orthotopically into the pulmonary position using an in-utero closed-heart hybrid approach. The transapical valve deployments were successful in all animals with acute survival of 77.8% of fetuses. TEHV in-vivo functionality was assessed using echocardiography as well as angiography. Fetuses were harvested up to 1 week after implantation representing a birth-relevant gestational age. TEHVs showed in vivo functionality with intact valvular integrity and absence of thrombus formation. The presented approach may serve as an experimental basis for future human prenatal cardiac interventions using fully biodegradable autologous cell-based living materials.

Multivariate QST-FST comparisons: a neutrality test for the evolution of the g matrix in structured populations.

Neutrality tests in quantitative genetics provide a statistical framework for the detection of selection on polygenic traits in wild populations. However, the existing method based on comparisons of divergence at neutral markers and quantitative traits (Q(st)-F(st)) suffers from several limitations that hinder a clear interpretation of the results with typical empirical designs. In this article, we propose a multivariate extension of this neutrality test based on empirical estimates of the among-populations (D) and within-populations (G) covariance matrices by MANOVA. A simple pattern is expected under neutrality: D = 2F(st)/(1 - F(st))G, so that neutrality implies both proportionality of the two matrices and a specific value of the proportionality coefficient. This pattern is tested using Flury's framework for matrix comparison [common principal-component (CPC) analysis], a well-known tool in G matrix evolution studies. We show the importance of using a Bartlett adjustment of the test for the small sample sizes typically found in empirical studies. We propose a dual test: (i) that the proportionality coefficient is not different from its neutral expectation [2F(st)/(1 - F(st))] and (ii) that the MANOVA estimates of mean square matrices between and among populations are proportional. These two tests combined provide a more stringent test for neutrality than the classic Q(st)-F(st) comparison and avoid several statistical problems. Extensive simulations of realistic empirical designs suggest that these tests correctly detect the expected pattern under neutrality and have enough power to efficiently detect mild to strong selection (homogeneous, heterogeneous, or mixed) when it is occurring on a set of traits. This method also provides a rigorous and quantitative framework for disentangling the effects of different selection regimes and of drift on the evolution of the G matrix. We discuss practical requirements for the proper application of our test in empirical studies and potential extensions.

Interest of qualitative and quantitative profiling of endocannabinoids for evaluation of their implication in drug addiction process

Introduction: In the middle of the 90's, the discovery of endogenous ligands for cannabinoid receptors opened a new era in this research field. Amides and esters of arachidonic acid have been identified as these endogenous ligands. Arachidonoylethanolamide (anandamide or AEA) and 2-Arachidonoylglycerol (2-AG) seem to be the most important of these lipid messengers. In addition, virodhamine (VA), noladin ether (2-AGE), and N-arachidonoyl dopamine (NADA) have been shown to bind to CB receptors with varying affinities. During recent years, it has become more evident that the EC system is part of fundamental regulatory mechanisms in many physiological processes such as stress and anxiety responses, depression, anorexia and bulimia, schizophrenia disorders, neuroprotection, Parkinson disease, anti-proliferative effects on cancer cells, drug addiction, and atherosclerosis. Aims: This work presents the problematic of EC analysis and the input of Information Dependant Acquisition based on hybrid triple quadrupole linear ion trap (QqQLIT) system for the profiling of these lipid mediators. Methods: The method was developed on a LC Ultimate 3000 series (Dionex, Sunnyvale, CA, USA) coupled to a QTrap 4000 system (Applied biosystems, Concord, ON, Canada). The ECs were separated on an XTerra C18 MS column (50 × 3.0 mm i.d., 3.5 μm) with a 5 min gradient elution. For confirmatory analysis, an information-dependant acquisition experiment was performed with selected reaction monitoring (SRM) as survey scan and enhanced produced ion (EPI) as dependant scan. Results: The assay was found to be linear in the concentration range of 0.1-5 ng/mL for AEA, 0.3-5 ng/mL for VA, 2-AGE, and NADA and 1-20 ng/mL for 2-AG using 0.5 mL of plasma. Repeatability and intermediate precision were found less than 15% over the tested concentration ranges. Under non-pathophysiological conditions, only AEA and 2-AG were actually detected in plasma with concentration ranges going from 104 to 537 pg/mL and from 2160 to 3990 pg/mL respectively. We have particularly focused our scopes on the evaluation of EC level changes in biological matrices through drug addiction and atherosclerosis processes. We will present preliminary data obtained during pilot study after administration of cannabis on human patients. Conclusion: ECs have been shown to play a key role in regulation of many pathophysiological processes. Medical research in these different fields continues to growth in order to understand and to highlight the predominant role of EC in the CNS and peripheral tissues signalisation. The profiling of these lipids needs to develop rapid, highly sensitive and selective analytical methods.

Connectome alterations in schizophrenia

Introduction : DTI has proven to be an exquisite biomarker of tissue microstructure integrity. This technique has been successfully applied to schizophrenia in showing that fractional anisotropy (FA, a marker of white matter integrity) is diminished in several areas of the brain (Kyriakopoulos M et al (2008)). New ways of representing diffusion data emerged recently and achieved to create structural connectivity maps in healthy brains (Hagmann P et al. (2008)). These maps have the capacity to study alterations over the entire brain at the connection and network level. This is of high interest in complex disconnection diseases like schizophrenia. We report on the specific network alterations of schizophrenic patients. Methods : 13 patients with chronic schizophrenia were recruited from in-patient, day treatment, out-patient clinics. Comparison subjects were recruited and group-matched to patients on age, sex, handedness, and parental social economic-status. This study was approved by the local IRB and subjects had to give informed written consent. They were scanned with a 3T clinical MRI scanner. DTI and high-resolution anatomical T1w imaging were performed during the same session. The path from diffusion MRI to a multi-resolution structural connection matrices of the entire brain is a five steps process that was performed in a similar way as described in Hagmann P et al. (2008). (1) DTI and T1w MRI of the brain, (2) segmentation of white and gray matter, (3) white matter tractography, (4) segmentation of the cortex into 242 ROIs of equal surface area covering the entire cortex (Fig 1), (5) the connection network was constructed by measuring for each ROI to ROI connection the related average FA along the corresponding tract. Results : For every connection between 2 ROIs of the network we tested the hypothesis H0: "average FA along fiber pathway is larger or equal in patients than in controls". H0 was rejected for connections where average FA in a connection was significantly lower in patients than in controls. Threshold p-value was 0.01 corrected for multiple comparisons with false discovery rate. We identified consistently that temporal, occipito-temporal, precuneo-temporal as well as frontal inferior and precuneo-cingulate connections were altered (Fig 2: significant connections in yellow). This is in agreement with the known literature, which showed across several studies that FA is diminished in several areas of the brain. More precisely, abnormalities were reported in the prefrontal and temporal white matter and to some extent also in the parietal and occipital regions. The alterations reported in the literature specifically included the corpus callosum, the arcuate fasciculus and the cingulum bundle, which was the case here as well. In addition small world indexes are significantly reduced in patients (p<0.01) (Fig. 3). Conclusions : Using connectome mapping to characterize differences in structural connectivity between healthy and diseased subjects we were able to show widespread connectional alterations in schizophrenia patients and systematic small worldness decrease, which is a marker of network desorganization. More generally, we described a method that has the capacity to sensitively identify structure alterations in complex disconnection syndromes where lesions are widespread throughout the connectional network.

The additive genetic variance after bottlenecks is affected by the number of loci involved in epistatic interactions.

We investigated the role of the number of loci coding for a neutral trait on the release of additive variance for this trait after population bottlenecks. Different bottleneck sizes and durations were tested for various matrices of genotypic values, with initial conditions covering the allele frequency space. We used three different types of matrices. First, we extended Cheverud and Routman's model by defining matrices of "pure" epistasis for three and four independent loci; second, we used genotypic values drawn randomly from uniform, normal, and exponential distributions; and third we used two models of simple metabolic pathways leading to physiological epistasis. For all these matrices of genotypic values except the dominant metabolic pathway, we find that, as the number of loci increases from two to three and four, an increase in the release of additive variance is occurring. The amount of additive variance released for a given set of genotypic values is a function of the inbreeding coefficient, independently of the size and duration of the bottleneck. The level of inbreeding necessary to achieve maximum release in additive variance increases with the number of loci. We find that additive-by-additive epistasis is the type of epistasis most easily converted into additive variance. For a wide range of models, our results show that epistasis, rather than dominance, plays a significant role in the increase of additive variance following bottlenecks.

Ultrastructural description of rabies virus infection in cultured sensory neurons

Primary cultures were made from adult mouse spinal ganglia for depicting an ultrastructural description of rabies virus (RABV) infection in adult mouse sensory neuron cultures; they were infected with rabies virus for 24, 36, and 48 h. The monolayers were processed for transmission electron microscopy and immunochemistry studies at the end of each period. As previously reported, sensory neurons showed great susceptibility to infection by RABV; however, in none of the periods evaluated were assembled virions observed in the cytoplasm or seen to be associated with the cytoplasmic membrane. Instead, fibril matrices of aggregated ribonucleoprotein were detected in the cytoplasm. When infected culture lysate were inoculated into normal animals via intra-cerebral route it was observed that these animals developed clinical symptoms characteristic of infection and transmission electron microscopy revealed assembled virions in the cerebral cortex and other areas of the brain. Sensory neurons infected in vitro by RABV produced a large amount of unassembled viral ribonucleoprotein. However, this intracellular material was able to produce infection and virions on being intra-cerebrally inoculated. It can thus be suggested that the lack of intracellular assembly in sensory neurons forms part of an efficient dissemination strategy.