Space Competition and Time Delays in Human Range Expansions. Application to the Neolithic Transition

Space competition effects are well-known in many microbiological and ecological systems. Here we analyze such an effectin human populations. The Neolithic transition (change from foraging to farming) was mainly the outcome of a demographic process that spread gradually throughout Europe from the Near East. In Northern Europe, archaeological data show a slowdown on the Neolithic rate of spread that can be related to a high indigenous (Mesolithic) population density hindering the advance as a result of the space competition between the two populations. We measure this slowdown from a database of 902 Early Neolithic sites and develop a time-delayed reaction-diffusion model with space competition between Neolithic and Mesolithic populations, to predict the observed speeds. The comparison of the predicted speed with the observations and with a previous non-delayed model show that both effects, the time delay effect due to the generation lag and the space competition between populations, are crucial in order to understand the observations

Phenotypic plasticity in salmonid embryos : characterizing pathogen-induced stress effects on heritable variation in traits and reaction norms

Les pressions écologiques peuvent varier tant en nature qu'en intensité dans le temps et l'espace. C'est pourquoi, un phénotype unique ne peut pas forcément conférer la meilleure valeur sélective. La plasticité phénotypique peut être un moyen de s'accommoder de cette situation, en augmentant globalement la tolérance aux changements environnementaux. Comme pour tout trait de caractère, une variation génétique doit persister pour qu'évoluent les traits plastiques dans une population donnée. Cependant, les pressions extérieures peuvent affecter l'héritabilité, et la direction de ces changements peut dépendre du caractère en question, de l'espèce mais aussi du type de stress. Dans la présente thèse, nous avons cherché à élucider les effets des pressions pathogéniques sur les phénotypes et la génétique quantitative de plusieurs traits plastiques chez les embryons de deux salmonidés, la palée (Coregonus palaea), et la truite de rivière (Salmo trutta). Les salmonidés se prêtent à de telles études du fait de leur extraordinaire variabilité morphologique, comportementale et des traits d'histoire de vie. Par ailleurs, avec le déclin des salmonidés dans le monde, il est important de savoir combien la variabilité génétique persiste dans les normes de réaction afin d'aider à prédire leur capacité à répondre aux changements de leur milieu. Nous avons observé qu'une augmentation de la croissance des communautés microbiennes symbiotiques entraînait une mortalité accrue et une éclosion précoce chez la palée, et dévoilait la variance génétique additive pour ces deux caractères (Chapitres 1-2). Bien qu'aucune variation génétique n'ait été trouvée pour les normes de réaction, nous avons observé une variabilité de la plasticité d'éclosion. Néanmoins, on a trouvé que les temps d'éclosion étaient corrélés entre les environnements, ce qui pourrait limiter l'évolution de la norme de réaction. Le temps d'éclosion des embryons est lié à la taille des géniteurs mâles, ce qui indique des effets pléiotropiques. Dans le Chapitre 3, nous avons montré qu'une interaction triple entre la souche bactérienne {Pseudomonas fluorescens}, l'état de dévelopement de l'hôte ainsi que ses gènes ont une influence sur la mortalité, le temps d'éclosion et la taille des alevins de la palée. Nous avons démontré qu'une variation génétique subsistait généralement dans les normes de réaction des temps d'éclosion, mais rarement pour la taille des alevins, et jamais pour la mortalité. Dans le même temps, nous avons exhibé que des corrélations entre environnements dépendaient des caractères phénotypiques, mais contrairement au Chapitre 2, nous n'avons pas trouvé de preuve de corrélations transgénérationnelles. Le Chapitre 4 complète le chapitre précédent, en se plaçant du point de vue moléculaire, et décrit comment le traitement d'embryons avec P. fluorescens s'est traduit par une régulation négative d'expression du CMH-I indépendemment de la souche bactérienne. Nous avons non seulement trouvé une variation génétique des caractères phénotypiques moyens, mais aussi de la plasticité. Les deux derniers chapitres traitent de l'investigation, chez la truite de rivière, des différences spécifiques entre populations pour des normes de réaction induites par les pathogènes. Dans le Chapitre 5, nous avons illustré que le métissage entre des populations génétiquement distinctes n'affectait en rien la hauteur ou la forme des normes de réaction d'un trait précoce d'histoire de vie suite au traitement pathogénique. De surcroît, en dépit de l'éclosion tardive et de la réduction de la taille des alevins, le traitement n'a pas modifié la variation héritable des traits de caractère. D'autre part, dans le Chapitre 6, nous avons démontré que le traitement d'embryons avec des stimuli contenus dans l'eau de conspécifiques infectés a entraîné des réponses propre à chaque population en terme de temps d'éclosion ; néanmoins, nous avons observé peu de variabilité génétique des normes de réaction pour ce temps d'éclosion au sein des populations. - Ecological stressors can vary in type and intensity over space and time, and as such, a single phenotype may not confer the highest fitness. Phenotypic plasticity can act as a means to accommodate this situation, increasing overall tolerance to environmental change. As with any trait, for plastic traits to evolve in a population, genetic variation must persist. However, environmental stress can alter trait heritability, and the direction of this shift can be trait, species, and stressor-dependent. In this thesis, we sought to understand the effects of pathogen stressors on the phenotypes and genetic architecture of several plastic traits in the embryos of two salmonids, the whitefish (Coregonus palaea), and the brown trout (Salmo trutta). Salmonids lend themselves to such studies because their extraordinary variability in morphological, behavioral, and life-history traits. Also, with declines in salmonids worldwide, knowing how much genetic variability persists in reaction norms may help predict their ability to respond to environmental change. We found that increasing growth of symbiotic microbial communities increased mortality and induced hatching in whitefish, and released additive genetic variance for both traits (Chapters 1-2). While no genetic variation was found for survival reaction norms, we did find variability in hatching plasticity. Nevertheless, hatching time was correlated across environments, which could constrain evolution of the reaction norm. Hatching time in the induced environment was also correlated to sire size, indicating pleiotropic effects. In Chapter 3 we report that a three-way interaction between bacterial strain (Pseudomonas fluorescens), host developmental stage, and host genetics impacted mortality, hatching time, and hatchling size in whitefish. We also showed that genetic variation generally persisted in hatching age reaction norms, but rarely for hatchling length, and never for mortality. At the same time, we demonstrated that cross-environmental correlations were trait-dependent, and unlike Chapter 2, we found no evidence of cross-generational correlations. Chapter 4 expands on the previous chapter, moving to the molecular level, and describes how treatment of embryos with P. fluorescens resulted in strain-independent downregulation of MHC class I. Genetic variation was evident not only in trait means, but also in plasticity. In the last two chapters, we investigated population level differences in pathogen- induced reaction norms in brown trout. In Chapter 5, we found that interbreeding between genetically distinct populations did not affect the elevation or shapes of the reaction norms of early life-history traits after pathogen challenge. Moreover, despite delaying hatching and reducing larval length, treatment produced no discernable shifts in heritable variation in traits. On the other hand, in Chapter 6, we found that treatment of embryos with water-borne cues from infected conspecifics elicited population-specific responses in terms of hatching time; however, we found little evidence of genetic variability in hatching reaction norms within populations. We have made considerable progress in understanding how pathogen stressors affect various early life-history traits in salmonid embryos. We have demonstrated that the effect of a particular stressor on heritable variation in these traits can vary according to the trait and species under consideration, in addition to the developmental stage of the host. Moreover, we found evidence of genetic variability in some, but not all reaction norms in whitefish and brown trout.

On the tractability of the piecewise-linear approximation for general discrete-choice network revenue management

The choice network revenue management (RM) model incorporates customer purchase behavioras customers purchasing products with certain probabilities that are a function of the offeredassortment of products, and is the appropriate model for airline and hotel network revenuemanagement, dynamic sales of bundles, and dynamic assortment optimization. The underlyingstochastic dynamic program is intractable and even its certainty-equivalence approximation, inthe form of a linear program called Choice Deterministic Linear Program (CDLP) is difficultto solve in most cases. The separation problem for CDLP is NP-complete for MNL with justtwo segments when their consideration sets overlap; the affine approximation of the dynamicprogram is NP-complete for even a single-segment MNL. This is in contrast to the independentclass(perfect-segmentation) case where even the piecewise-linear approximation has been shownto be tractable. In this paper we investigate the piecewise-linear approximation for network RMunder a general discrete-choice model of demand. We show that the gap between the CDLP andthe piecewise-linear bounds is within a factor of at most 2. We then show that the piecewiselinearapproximation is polynomially-time solvable for a fixed consideration set size, bringing itinto the realm of tractability for small consideration sets; small consideration sets are a reasonablemodeling tradeoff in many practical applications. Our solution relies on showing that forany discrete-choice model the separation problem for the linear program of the piecewise-linearapproximation can be solved exactly by a Lagrangian relaxation. We give modeling extensionsand show by numerical experiments the improvements from using piecewise-linear approximationfunctions.

First detection of Waddlia chondrophila in Africa using SYBR Green real-time PCR on veterinary samples.

Waddlia chondrophila is a strict intracellular microorganism belonging to the order Chlamydiales that has been isolated twice from aborted bovine fetuses, once in USA and once in Germany. This bacterium is now considered as an abortigenic agent in cattle. However, no information is available regarding the presence of this bacterium in Africa. Given the low sensitivity of cell culture to recover such an obligate intracellular bacterium, molecular-based diagnostic approaches are warranted. This report describes the development of a quantitative SYBR Green real-time PCR assay targeting the recA gene of W. chondrophila. Analytical sensitivity was 10 copies of control plasmid DNA per reaction. No cross-amplification was observed when testing pathogens that can cause abortion in cattle. The PCR exhibited a good intra-run and inter-run reproducibility. This real-time PCR was then applied to 150 vaginal swabs taken from Tunisian cows that have aborted. Twelve samples revealed to be Waddlia positive, suggesting a possible role of this bacterium in this setting. This new real-time PCR assay represents a diagnostic tool that may be used to further study the prevalence of Waddlia infection.

Schizotypal perceptual aberrations of time: correlation between score, behavior and brain activity.

A fundamental trait of the human self is its continuum experience of space and time. Perceptual aberrations of this spatial and temporal continuity is a major characteristic of schizophrenia spectrum disturbances--including schizophrenia, schizotypal personality disorder and schizotypy. We have previously found the classical Perceptual Aberration Scale (PAS) scores, related to body and space, to be positively correlated with both behavior and temporo-parietal activation in healthy participants performing a task involving self-projection in space. However, not much is known about the relationship between temporal perceptual aberration, behavior and brain activity. To this aim, we composed a temporal Perceptual Aberration Scale (tPAS) similar to the traditional PAS. Testing on 170 participants suggested similar performance for PAS and tPAS. We then correlated tPAS and PAS scores to participants' performance and neural activity in a task of self-projection in time. tPAS scores correlated positively with reaction times across task conditions, as did PAS scores. Evoked potential mapping and electrical neuroimaging showed self-projection in time to recruit a network of brain regions at the left anterior temporal cortex, right temporo-parietal junction, and occipito-temporal cortex, and duration of activation in this network positively correlated with tPAS and PAS scores. These data demonstrate that schizotypal perceptual aberrations of both time and space, as reflected by tPAS and PAS scores, are positively correlated with performance and brain activation during self-projection in time in healthy individuals along the schizophrenia spectrum.

Evaluation of a polymerase chain reaction-restriction fragment length polymorphism assay for dermatophyte and nondermatophyte identification in onychomycosis.

BACKGROUND: Dermatophytes are the main cause of onychomycoses, but various nondermatophyte filamentous fungi are often isolated from abnormal nails. The correct identification of the aetiological agent of nail infections is necessary in order to recommend appropriate treatment. OBJECTIVE: To evaluate a rapid polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) assay based on 28S rDNA for fungal identification in nails on a large number of samples in comparison with cultures. METHODS: Infectious fungi were analysed using PCR-RFLP in 410 nail samples in which fungal elements were observed in situ by direct mycological examination (positive samples). The results were compared with those previously obtained by culture of fungi on Sabouraud agar from the same nail samples. RESULTS: PCR-RFLP identification of fungi in nails allowed validation of the results obtained in culture when Trichophyton spp. grew from infected samples. In addition, nondermatophyte filamentous fungi could be identified with certainty as the infectious agents in onychomycosis, and discriminated from dermatophytes as well as from transient contaminants. The specificity of the culture results relative to PCR-RFLP appeared to be 81%, 71%, 52% and 63% when Fusarium spp., Scopulariopsis brevicaulis, Aspergillus spp. and Candida spp., respectively, grew on Sabouraud agar. It was also possible to identify the infectious agent when direct nail mycological examination showed fungal elements, but negative results were obtained from fungal culture. CONCLUSIONS: Improved sensitivity for the detection of fungi in nails was obtained using the PCR-RFLP assay. Rapid and reliable molecular identification of the infectious fungus can be used routinely and presents several important advantages compared with culture in expediting the choice of appropriate antifungal therapy.

Changes in running mechanics and spring-mass behaviour during a 5-km time trial.

Research into the biomechanical manifestation of fatigue during exhaustive runs is increasingly popular but additional understanding of the adaptation of the spring-mass behaviour during the course of strenuous, self-paced exercises continues to be a challenge in order to develop optimized training and injury prevention programs. This study investigated continuous changes in running mechanics and spring-mass behaviour during a 5-km run. 12 competitive triathletes performed a 5-km running time trial (mean performance: ̴17 min 30 s) on a 200 m indoor track. Vertical and anterior-posterior ground reaction forces were measured every 200 m by a 5-m long force platform system, and used to determine spring-mass model characteristics. After a fast start, running velocity progressively decreased (- 11.6%; P<0.001) in the middle part of the race before an end spurt in the final 400-600 m. Stride length (- 7.4%; P<0.001) and frequency (- 4.1%; P=0.001) decreased over the 25 laps, while contact time (+ 8.9%; P<0.001) and total stride duration (+ 4.1%; P<0.001) progressively lengthened. Peak vertical forces (- 2.0%; P<0.01) and leg compression (- 4.3%; P<0.05), but not centre of mass vertical displacement (+ 3.2%; P>0.05), decreased with time. As a result, vertical stiffness decreased (- 6.0%; P<0.001) during the run, whereas leg stiffness changes were not significant (+ 1.3%; P>0.05). Spring-mass behaviour progressively changes during a 5-km time trial towards deteriorated vertical stiffness, which alters impact and force production characteristics.

Time-delayed fronts from biased random walks

We generalize a previous model of time-delayed reaction–diffusion fronts (Fort and Méndez 1999 Phys. Rev. Lett. 82 867) to allow for a bias in the microscopic random walk of particles or individuals. We also present a second model which takes the time order of events (diffusion and reproduction) into account. As an example, we apply them to the human invasion front across the USA in the 19th century. The corrections relative to the previous model are substantial. Our results are relevant to physical and biological systems with anisotropic fronts, including particle diffusion in disordered lattices, population invasions, the spread of epidemics, etc

Development of a duplex real-time PCR for the detection of Rickettsia spp. and typhus group rickettsia in clinical samples.

Molecular diagnosis using real-time polymerase chain reaction (PCR) may allow earlier diagnosis of rickettsiosis. We developed a duplex real-time PCR that amplifies (1) DNA of any rickettsial species and (2) DNA of both typhus group rickettsia, that is, Rickettsia prowazekii and Rickettsia typhi. Primers and probes were selected to amplify a segment of the 16S rRNA gene of Rickettsia spp. for the pan-rickettsial PCR and the citrate synthase gene (gltA) for the typhus group rickettsia PCR. Analytical sensitivity was 10 copies of control plasmid DNA per reaction. No cross-amplification was observed when testing human DNA and 22 pathogens or skin commensals. Real-time PCR was applied to 16 clinical samples. Rickettsial DNA was detected in the skin biopsies of three patients. In one patient with severe murine typhus, the typhus group PCR was positive in a skin biopsy from a petechial lesion and seroconversion was later documented. The two other patients with negative typhus group PCR suffered from Mediterranean and African spotted fever, respectively; in both cases, skin biopsy was performed on the eschar. Our duplex real-time PCR showed a good analytical sensitivity and specificity, allowing early diagnosis of rickettsiosis among three patients, and recognition of typhus in one of them.

A matter of time: delayed mate encounter postpones mating window initiation and reduces the strength of female choosiness

Reproductive success is determined by the presence and timing of encounter of mates. The latter depends on species-specific reproductive characteristics (e.g. initiation/duration of the mating window), season, and reproductive strategies (e.g. intensity of choosiness) that may potentially mitigate constraints imposed by mating windows. Despite their potentially crucial role for fitness and population dynamics, limited evidence exists about mating window initiation, duration and reproductive strategies. Here, we experimentally tested the mechanisms of initiation and the duration of the common lizard's Zootoca vivipara mating window, by manipulating the timing of mate encounter and analyzing its effect on (re-)mating probability. We furthermore tested treatment effects on female reproductive strategies, by measuring female choosiness. The timing of mate encounter and season did not significantly affect mating probability. However, a longer delay until mate encounter reduced female choosiness. Re-mating probability decreased with re-mating delay and was independent of mating delay. This indicates that mating window initiation depends on mate encounter, that its duration is fixed, and that plastic reproductive strategies exist. These findings contrast with previous beliefs and shows that mating windows per se may not necessarily constrain reproductive success, which is congruent with rapid range expansion and absence of positive density-effects on reproductive success (Allee effects). In summary, our results show that predicting the effect of mating windows on reproduction is complex and that experimental evidence is essential for evaluating their effect on reproduction and reproductive strategies, both being important determinants of population dynamics and the colonization of new habitats.

New diagnostic real-time PCR for specific detection of Parachlamydia acanthamoebae DNA in clinical samples

Given the low sensitivity of amoebal coculture, we developed a specific real-time PCR for the detection of Parachlamydia. The analytical sensitivity was high, and the inter- and intrarun variabilities were low. When the PCR was applied to nasopharyngeal aspirates, it was positive for six patients with bronchiolitis. Future studies should assess the role of Parachlamydia in bronchiolitis.

A case of acute retinal pigment epithelitis: spectral domain optical coherence tomography time course and physiopathologic hypothesis.

PURPOSE: To report the time course of retinal morphologic changes in a patient with acute retinal pigment epithelitis (ARPE) using spectral domain optical coherence tomography (SD-OCT). METHODS: A 30-year old man was referred for blurred vision of his right eye after five days that appeared suddenly 15 days after recovery from a flu-like syndrome. SD-OCT was performed immediately, followed by fluorescein and infracyanine angiography at eight days and then at three weeks. RESULTS: At presentation, a bubble of sub-macular deposit was observed on the right macula with central golden micronodules in a honeycomb pattern. SD-OCT showed an "anterior dislocation" of all the retinal layers up to the inner/outer segment (IS/OS) line and irregular deposits at the OS level together with thickening of the retinal pigment epithelial (RPE) layer. As visual acuity increased, eight days later, the OCT showed reduction of the sub-retinal deposits and an abnormal hyperflectivity of the sub-retinal and RPE layers was observed. The patient showed a positive serology for picornavirus. DISCUSSION: The acute SD-OCT sections of this patient with ARPE were compared with histological sections of a 35 day old Royal College of Surgeons rat. Similar findings could be observed, with preservation of the IS/OS line and accumulation of debris at the OS level, suggesting that ARPE symptoms could result from a transient phagocytic dysfunction of the RPE at the fovea, inducing reversible accumulation of undigested OS. Picornaviruses comprising enterovirus and coxsachievirus described as being associated with acute chorioretinitis. In this case, it was responsible for ARPE. CONCLUSION: We hypothesize that ARPE syndrome results from a transient dysfunction of RPE, which can occur as a post viral reaction.

Artificial evolution on network structures : how time and space influence dynamics

Abstract The main objective of this work is to show how the choice of the temporal dimension and of the spatial structure of the population influences an artificial evolutionary process. In the field of Artificial Evolution we can observe a common trend in synchronously evolv¬ing panmictic populations, i.e., populations in which any individual can be recombined with any other individual. Already in the '90s, the works of Spiessens and Manderick, Sarma and De Jong, and Gorges-Schleuter have pointed out that, if a population is struc¬tured according to a mono- or bi-dimensional regular lattice, the evolutionary process shows a different dynamic with respect to the panmictic case. In particular, Sarma and De Jong have studied the selection pressure (i.e., the diffusion of a best individual when the only selection operator is active) induced by a regular bi-dimensional structure of the population, proposing a logistic modeling of the selection pressure curves. This model supposes that the diffusion of a best individual in a population follows an exponential law. We show that such a model is inadequate to describe the process, since the growth speed must be quadratic or sub-quadratic in the case of a bi-dimensional regular lattice. New linear and sub-quadratic models are proposed for modeling the selection pressure curves in, respectively, mono- and bi-dimensional regu¬lar structures. These models are extended to describe the process when asynchronous evolutions are employed. Different dynamics of the populations imply different search strategies of the resulting algorithm, when the evolutionary process is used to solve optimisation problems. A benchmark of both discrete and continuous test problems is used to study the search characteristics of the different topologies and updates of the populations. In the last decade, the pioneering studies of Watts and Strogatz have shown that most real networks, both in the biological and sociological worlds as well as in man-made structures, have mathematical properties that set them apart from regular and random structures. In particular, they introduced the concepts of small-world graphs, and they showed that this new family of structures has interesting computing capabilities. Populations structured according to these new topologies are proposed, and their evolutionary dynamics are studied and modeled. We also propose asynchronous evolutions for these structures, and the resulting evolutionary behaviors are investigated. Many man-made networks have grown, and are still growing incrementally, and explanations have been proposed for their actual shape, such as Albert and Barabasi's preferential attachment growth rule. However, many actual networks seem to have undergone some kind of Darwinian variation and selection. Thus, how these networks might have come to be selected is an interesting yet unanswered question. In the last part of this work, we show how a simple evolutionary algorithm can enable the emrgence o these kinds of structures for two prototypical problems of the automata networks world, the majority classification and the synchronisation problems. Synopsis L'objectif principal de ce travail est de montrer l'influence du choix de la dimension temporelle et de la structure spatiale d'une population sur un processus évolutionnaire artificiel. Dans le domaine de l'Evolution Artificielle on peut observer une tendence à évoluer d'une façon synchrone des populations panmictiques, où chaque individu peut être récombiné avec tout autre individu dans la population. Déjà dans les année '90, Spiessens et Manderick, Sarma et De Jong, et Gorges-Schleuter ont observé que, si une population possède une structure régulière mono- ou bi-dimensionnelle, le processus évolutionnaire montre une dynamique différente de celle d'une population panmictique. En particulier, Sarma et De Jong ont étudié la pression de sélection (c-à-d la diffusion d'un individu optimal quand seul l'opérateur de sélection est actif) induite par une structure régulière bi-dimensionnelle de la population, proposant une modélisation logistique des courbes de pression de sélection. Ce modèle suppose que la diffusion d'un individu optimal suit une loi exponentielle. On montre que ce modèle est inadéquat pour décrire ce phénomène, étant donné que la vitesse de croissance doit obéir à une loi quadratique ou sous-quadratique dans le cas d'une structure régulière bi-dimensionnelle. De nouveaux modèles linéaires et sous-quadratique sont proposés pour des structures mono- et bi-dimensionnelles. Ces modèles sont étendus pour décrire des processus évolutionnaires asynchrones. Différentes dynamiques de la population impliquent strategies différentes de recherche de l'algorithme résultant lorsque le processus évolutionnaire est utilisé pour résoudre des problèmes d'optimisation. Un ensemble de problèmes discrets et continus est utilisé pour étudier les charactéristiques de recherche des différentes topologies et mises à jour des populations. Ces dernières années, les études de Watts et Strogatz ont montré que beaucoup de réseaux, aussi bien dans les mondes biologiques et sociologiques que dans les structures produites par l'homme, ont des propriétés mathématiques qui les séparent à la fois des structures régulières et des structures aléatoires. En particulier, ils ont introduit la notion de graphe sm,all-world et ont montré que cette nouvelle famille de structures possède des intéressantes propriétés dynamiques. Des populations ayant ces nouvelles topologies sont proposés, et leurs dynamiques évolutionnaires sont étudiées et modélisées. Pour des populations ayant ces structures, des méthodes d'évolution asynchrone sont proposées, et la dynamique résultante est étudiée. Beaucoup de réseaux produits par l'homme se sont formés d'une façon incrémentale, et des explications pour leur forme actuelle ont été proposées, comme le preferential attachment de Albert et Barabàsi. Toutefois, beaucoup de réseaux existants doivent être le produit d'un processus de variation et sélection darwiniennes. Ainsi, la façon dont ces structures ont pu être sélectionnées est une question intéressante restée sans réponse. Dans la dernière partie de ce travail, on montre comment un simple processus évolutif artificiel permet à ce type de topologies d'émerger dans le cas de deux problèmes prototypiques des réseaux d'automates, les tâches de densité et de synchronisation.

Selective catalytic reduction of nitrogen oxides with ammonia in forced unsteady state reactors - Case based reasoning and mathematical model simulation reasoning

The application of forced unsteady-state reactors in case of selective catalytic reduction of nitrogen oxides (NOx) with ammonia (NH3) is sustained by the fact that favorable temperature and composition distributions which cannot be achieved in any steady-state regime can be obtained by means of unsteady-state operations. In a normal way of operation the low exothermicity of the selective catalytic reduction (SCR) reaction (usually carried out in the range of 280-350°C) is not enough to maintain by itself the chemical reaction. A normal mode of operation usually requires supply of supplementary heat increasing in this way the overall process operation cost. Through forced unsteady-state operation, the main advantage that can be obtained when exothermic reactions take place is the possibility of trapping, beside the ammonia, the moving heat wave inside the catalytic bed. The unsteady state-operation enables the exploitation of the thermal storage capacity of the catalyticbed. The catalytic bed acts as a regenerative heat exchanger allowing auto-thermal behaviour when the adiabatic temperature rise is low. Finding the optimum reactor configuration, employing the most suitable operation model and identifying the reactor behavior are highly important steps in order to configure a proper device for industrial applications. The Reverse Flow Reactor (RFR) - a forced unsteady state reactor - corresponds to the above mentioned characteristics and may be employed as an efficient device for the treatment of dilute pollutant mixtures. As a main disadvantage, beside its advantages, the RFR presents the 'wash out' phenomena. This phenomenon represents emissions of unconverted reactants at every switch of the flow direction. As a consequence our attention was focused on finding an alternative reactor configuration for RFR which is not affected by the incontrollable emissions of unconverted reactants. In this respect the Reactor Network (RN) was investigated. Its configuration consists of several reactors connected in a closed sequence, simulating a moving bed by changing the reactants feeding position. In the RN the flow direction is maintained in the same way ensuring uniformcatalyst exploitation and in the same time the 'wash out' phenomena is annulated. The simulated moving bed (SMB) can operate in transient mode giving practically constant exit concentration and high conversion levels. The main advantage of the reactor network operation is emphasizedby the possibility to obtain auto-thermal behavior with nearly uniformcatalyst utilization. However, the reactor network presents only a small range of switching times which allow to reach and to maintain an ignited state. Even so a proper study of the complex behavior of the RN may give the necessary information to overcome all the difficulties that can appear in the RN operation. The unsteady-state reactors complexity arises from the fact that these reactor types are characterized by short contact times and complex interaction between heat and mass transportphenomena. Such complex interactions can give rise to a remarkable complex dynamic behavior characterized by a set of spatial-temporal patterns, chaotic changes in concentration and traveling waves of heat or chemical reactivity. The main efforts of the current research studies concern the improvement of contact modalities between reactants, the possibility of thermal wave storage inside the reactor and the improvement of the kinetic activity of the catalyst used. Paying attention to the above mentioned aspects is important when higher activity even at low feeding temperatures and low emissions of unconverted reactants are the main operation concerns. Also, the prediction of the reactor pseudo or steady-state performance (regarding the conversion, selectivity and thermal behavior) and the dynamicreactor response during exploitation are important aspects in finding the optimal control strategy for the forced unsteady state catalytic tubular reactors. The design of an adapted reactor requires knowledge about the influence of its operating conditions on the overall process performance and a precise evaluation of the operating parameters rage for which a sustained dynamic behavior is obtained. An apriori estimation of the system parameters result in diminution of the computational efforts. Usually the convergence of unsteady state reactor systems requires integration over hundreds of cycles depending on the initial guess of the parameter values. The investigation of various operation models and thermal transfer strategies give reliable means to obtain recuperative and regenerative devices which are capable to maintain an auto-thermal behavior in case of low exothermic reactions. In the present research work a gradual analysis of the SCR of NOx with ammonia process in forced unsteady-state reactors was realized. The investigation covers the presentationof the general problematic related to the effect of noxious emissions in the environment, the analysis of the suitable catalysts types for the process, the mathematical analysis approach for modeling and finding the system solutions and the experimental investigation of the device found to be more suitable for the present process. In order to gain information about the forced unsteady state reactor design, operation, important system parameters and their values, mathematical description, mathematicalmethod for solving systems of partial differential equations and other specific aspects, in a fast and easy way, and a case based reasoning (CBR) approach has been used. This approach, using the experience of past similarproblems and their adapted solutions, may provide a method for gaining informations and solutions for new problems related to the forced unsteady state reactors technology. As a consequence a CBR system was implemented and a corresponding tool was developed. Further on, grooving up the hypothesis of isothermal operation, the investigation by means of numerical simulation of the feasibility of the SCR of NOx with ammonia in the RFRand in the RN with variable feeding position was realized. The hypothesis of non-isothermal operation was taken into account because in our opinion ifa commercial catalyst is considered, is not possible to modify the chemical activity and its adsorptive capacity to improve the operation butis possible to change the operation regime. In order to identify the most suitable device for the unsteady state reduction of NOx with ammonia, considering the perspective of recuperative and regenerative devices, a comparative analysis of the above mentioned two devices performance was realized. The assumption of isothermal conditions in the beginningof the forced unsteadystate investigation allowed the simplification of the analysis enabling to focus on the impact of the conditions and mode of operation on the dynamic features caused by the trapping of one reactant in the reactor, without considering the impact of thermal effect on overall reactor performance. The non-isothermal system approach has been investigated in order to point out the important influence of the thermal effect on overall reactor performance, studying the possibility of RFR and RN utilization as recuperative and regenerative devices and the possibility of achieving a sustained auto-thermal behavior in case of lowexothermic reaction of SCR of NOx with ammonia and low temperature gasfeeding. Beside the influence of the thermal effect, the influence of the principal operating parameters, as switching time, inlet flow rate and initial catalyst temperature have been stressed. This analysis is important not only because it allows a comparison between the two devices and optimisation of the operation, but also the switching time is the main operating parameter. An appropriate choice of this parameter enables the fulfilment of the process constraints. The level of the conversions achieved, the more uniform temperature profiles, the uniformity ofcatalyst exploitation and the much simpler mode of operation imposed the RN as a much more suitable device for SCR of NOx with ammonia, in usual operation and also in the perspective of control strategy implementation. Theoretical simplified models have also been proposed in order to describe the forced unsteady state reactors performance and to estimate their internal temperature and concentration profiles. The general idea was to extend the study of catalytic reactor dynamics taking into account the perspectives that haven't been analyzed yet. The experimental investigation ofRN revealed a good agreement between the data obtained by model simulation and the ones obtained experimentally.

Discovery of a 29-gene panel in peripheral blood mononuclear cells for the detection of colorectal cancer and adenomas using high throughput real-time PCR.

Colorectal cancer (CRC) is the second leading cause of cancer-related death in developed countries. Early detection of CRC leads to decreased CRC mortality. A blood-based CRC screening test is highly desirable due to limited invasiveness and high acceptance rate among patients compared to currently used fecal occult blood testing and colonoscopy. Here we describe the discovery and validation of a 29-gene panel in peripheral blood mononuclear cells (PBMC) for the detection of CRC and adenomatous polyps (AP). Blood samples were prospectively collected from a multicenter, case-control clinical study. First, we profiled 93 samples with 667 candidate and 3 reference genes by high throughput real-time PCR (OpenArray system). After analysis, 160 genes were retained and tested again on 51 additional samples. Low expressed and unstable genes were discarded resulting in a final dataset of 144 samples profiled with 140 genes. To define which genes, alone or in combinations had the highest potential to discriminate AP and/or CRC from controls, data were analyzed by a combination of univariate and multivariate methods. A list of 29 potentially discriminant genes was compiled and evaluated for its predictive accuracy by penalized logistic regression and bootstrap. This method discriminated AP >1cm and CRC from controls with a sensitivity of 59% and 75%, respectively, with 91% specificity. The behavior of the 29-gene panel was validated with a LightCycler 480 real-time PCR platform, commonly adopted by clinical laboratories. In this work we identified a 29-gene panel expressed in PBMC that can be used for developing a novel minimally-invasive test for accurate detection of AP and CRC using a standard real-time PCR platform.