A multi-scaled approach for simulating chemical reaction systems

In this paper we give an overview of some very recent work, as well as presenting a new approach, on the stochastic simulation of multi-scaled systems involving chemical reactions. In many biological systems (such as genetic regulation and cellular dynamics) there is a mix between small numbers of key regulatory proteins, and medium and large numbers of molecules. In addition, it is important to be able to follow the trajectories of individual molecules by taking proper account of the randomness inherent in such a system. We describe different types of simulation techniques (including the stochastic simulation algorithm, Poisson Runge-Kutta methods and the balanced Euler method) for treating simulations in the three different reaction regimes: slow, medium and fast. We then review some recent techniques on the treatment of coupled slow and fast reactions for stochastic chemical kinetics and present a new approach which couples the three regimes mentioned above. We then apply this approach to a biologically inspired problem involving the expression and activity of LacZ and LacY proteins in E coli, and conclude with a discussion on the significance of this work. (C) 2004 Elsevier Ltd. All rights reserved.

Identification of Mendel's white flower character

Background The genetic regulation of flower color has been widely studied, notably as a character used by Mendel and his predecessors in the study of inheritance in pea. Methodology/Principal Findings We used the genome sequence of model legumes, together with their known synteny to the pea genome to identify candidate genes for the A and A2 loci in pea. We then used a combination of genetic mapping, fast neutron mutant analysis, allelic diversity, transcript quantification and transient expression complementation studies to confirm the identity of the candidates. Conclusions/Significance We have identified the pea genes A and A2. A is the factor determining anthocyanin pigmentation in pea that was used by Gregor Mendel 150 years ago in his study of inheritance. The A gene encodes a bHLH transcription factor. The white flowered mutant allele most likely used by Mendel is a simple G to A transition in a splice donor site that leads to a mis-spliced mRNA with a premature stop codon, and we have identified a second rare mutant allele. The A2 gene encodes a WD40 protein that is part of an evolutionarily conserved regulatory complex.

DNA methylation at the novel CpG sites in the promoter of MED15/PCQAP gene as a biomarker for head and neck cancers

Head and neck cancers (HNCs) represent a significant and ever-growing burden to the modern society, mainly due to the lack of early diagnostic methods. A significant number of HNCs is often associated with drinking, smoking, chewing beetle nut, and human papilloma virus (HPV) infections. We have analyzed DNA methylation patterns in tumor and normal tissue samples collected from head and neck squamous cell carcinoma (HNSCC) patients who were smokers. We have identified novel methylation sites in the promoter of the mediator complex subunit 15 (MED15/PCQAP) gene (encoing a co-factor important for regulation of transcription initiation for promoters of many genes), hypermethylated specifically in tumor cells. Two clusters of CpG dinucleotides methylated in tumors, but not in normal tissue from the same patients, were identified. These CpG methylation events in saliva samples were further validated in a separate cohort of HNSCC patients (who developed cancer due to smoking or HPV infections) and healthy controls using methylation-specific PCR (MSP). We used saliva as a biological medium because of its non-invasive nature, close proximity to the tumors, easiness and it is an economically viable option for large-scale screening studies. The methylation levels for the two identified CpG clusters were significantly different between the saliva samples collected from healthy controls and HNSCC individuals (Welch's t-test returning P, 0.05 and Mann-Whitney test P, 0.01 for both). The developed MSP assays also provided a good discriminative ability with AUC values of 0.70 (P, 0.01) and 0.63 (P, 0.05). The identified novel CpG methylation sites may serve as potential non-invasive biomarkers for detecting HNSCC. © the authors.

A genome-wide association study identifies five loci influencing facial morphology in Europeans

Inter-individual variation in facial shape is one of the most noticeable phenotypes in humans, and it is clearly under genetic regulation; however, almost nothing is known about the genetic basis of normal human facial morphology. We therefore conducted a genome-wide association study for facial shape phenotypes in multiple discovery and replication cohorts, considering almost ten thousand individuals of European descent from several countries. Phenotyping of facial shape features was based on landmark data obtained from three-dimensional head magnetic resonance images (MRIs) and two-dimensional portrait images. We identified five independent genetic loci associated with different facial phenotypes, suggesting the involvement of five candidate genes-PRDM16, PAX3, TP63, C5orf50, and COL17A1-in the determination of the human face. Three of them have been implicated previously in vertebrate craniofacial development and disease, and the remaining two genes potentially represent novel players in the molecular networks governing facial development. Our finding at PAX3 influencing the position of the nasion replicates a recent GWAS of facial features. In addition to the reported GWA findings, we established links between common DNA variants previously associated with NSCL/P at 2p21, 8q24, 13q31, and 17q22 and normal facial-shape variations based on a candidate gene approach. Overall our study implies that DNA variants in genes essential for craniofacial development contribute with relatively small effect size to the spectrum of normal variation in human facial morphology. This observation has important consequences for future studies aiming to identify more genes involved in the human facial morphology, as well as for potential applications of DNA prediction of facial shape such as in future forensic applications.

Familial isolated hyperparathyroidism is linked to a 1.7 Mb region on chromosome 2p13.3-14

BACKGROUND: Familial isolated hyperparathyroidism (FIHP) is an autosomal dominantly inherited form of primary hyperparathyroidism. Although comprising only about 1% of cases of primary hyperparathyroidism, identification and functional analysis of a causative gene for FIHP is likely to advance our understanding of parathyroid physiology and pathophysiology. METHODS: A genome-wide screen of DNA from seven pedigrees with FIHP was undertaken in order to identify a region of genetic linkage with the disorder. RESULTS: Multipoint linkage analysis identified a region of suggestive linkage (LOD score 2.68) on chromosome 2. Fine mapping with the addition of three other families revealed significant linkage adjacent to D2S2368 (maximum multipoint LOD score 3.43). Recombination events defined a 1.7 Mb region of linkage between D2S2368 and D2S358 in nine pedigrees. Sequencing of the two most likely candidate genes in this region, however, did not identify a gene for FIHP. CONCLUSIONS: We conclude that a causative gene for FIHP lies within this interval on chromosome 2. This is a major step towards eventual precise identification of a gene for FIHP, likely to be a key component in the genetic regulation of calcium homeostasis.

Differential Hox Expression in Murine Embryonic Stem Cell Models of Normal and Malignant Hematopoiesis.

The Hox family are master transcriptional regulators of developmental processes, including hematopoiesis. The Hox regulators, caudal homeobox factors (Cdx1-4), and Meis1, along with several individual Hox proteins, are implicated in stem cell expansion during embryonic development, with gene dosage playing a significant role in the overall function of the integrated Hox network. To investigate the role of this network in normal and aberrant, early hematopoiesis, we employed an in vitro embryonic stem cell differentiation system, which recapitulates mouse developmental hematopoiesis. Expression profiles of Hox, Pbx1, and Meis1 genes were quantified at distinct stages during the hematopoietic differentiation process and compared with the effects of expressing the leukemic oncogene Tel/PDGFRß. During normal differentiation the Hoxa cluster, Pbx1 and Meis1 predominated, with a marked reduction in the majority of Hox genes (27/39) and Meis1 occurring during hematopoietic commitment. Only the posterior Hoxa cluster genes (a9, a10, a11, and a13) maintained or increased expression at the hematopoietic colony stage. Cdx4, Meis1, and a subset of Hox genes, including a7 and a9, were differentially expressed after short-term oncogenic (Tel/PDGFRß) induction. Whereas Hoxa4-10, b1, b2, b4, and b9 were upregulated during oncogenic driven myelomonocytic differentiation. Heterodimers between Hoxa7/Hoxa9, Meis1, and Pbx have previously been implicated in regulating target genes involved in hematopoietic stem cell (HSC) expansion and leukemic progression. These results provide direct evidence that transcriptional flux through the Hox network occurs at very early stages during hematopoietic differentiation and validates embryonic stem cell models for gaining insights into the genetic regulation of normal and malignant hematopoiesis.

Variation in Rice Cadmium Related to Human Exposure.

Cereal grains are the dominant source of cadmium in the human diet, with rice being to the fore. Here we explore the effect of geographic, genetic, and processing (milling) factors on rice grain cadmium and rice consumption rates that lead to dietary variance in cadmium intake. From a survey of 12 countries on four continents, cadmium levels in rice grain were the highest in Bangladesh and Sri Lanka, with both these countries also having high per capita rice intakes. For Bangladesh and Sri Lanka, there was high weekly intake of cadmium from rice, leading to intakes deemed unsafe by international and national regulators. While genetic variance, and to a lesser extent milling, provide strategies for reducing cadmium in rice, caution has to be used, as there is environmental regulation as well as genetic regulation of cadmium accumulation within rice grains. For countries that import rice, grain cadmium can be controlled by where that rice is sourced, but for countries with subsistence rice economies that have high levels of cadmium in rice grain, agronomic and breeding strategies are required to lower grain cadmium.

Silicon, the silver bullet for mitigating biotic and abiotic stress, and improving grain quality, in rice?

Adequate silicon fertilization greatly boosts rice yield and mitigates biotic and abiotic stress, and improves grain quality through lowering the content of cadmium and inorganic arsenic. This review on silicon dynamics in rice considers recent advances in our understanding of the role of silicon in rice, and the challenges of maintaining adequate silicon fertility within rice paddy systems. Silicon is increasingly considered as an element required for optimal plant performance, particularly in rice. Plants can survive with very low silicon under laboratory/glasshouse conditions, but this is highly artificial and, thus, silicon can be considered as essential for proper plant function in its environment. Silicon is incorporated into structural components of rice cell walls were it increases cell and tissue rigidity in the plant. Structural silicon provides physical protection to plants against microbial infection and insect attack as well as reducing the quality of the tissue to the predating organisms. The abiotic benefits are due to silicon's effect on overall organ strength. This helps protect against lodging, drought stress, high temperature (through efficient maintenance of transpiration), and photosynthesis by protecting against high UV. Furthermore, silicon also protects the plant from saline stress and against a range of toxic metal stresses (arsenic, cadmium, chromium, copper, nickel and zinc). Added to this, silicon application decreases grain concentrations of various human carcinogens, in particular arsenic, antimony and cadmium. As rice is efficient at stripping bioavailable silicon from the soil, recycling of silicon rich rice straw biomass or addition of inorganic silicon fertilizer, primarily obtained from iron and steel slag, needs careful management. Silicon in the soil may be lost if the silicon-cycle, traditionally achieved via composting of rice straw and returning it to the land, is being broken. As composting of rice straw and incorporation of composted or non-composted straw back to land are resource intensive activities, these activities are declining due to population shifts from the countryside to cities. Processes that accelerate rice straw composting, therefore, need to be identified to aid more efficient use of this resource. In addition, rice genetics may help address declining available silicon in paddy soils: for example by selecting for characteristics during breeding that lead to an increased ability of roots to access recalcitrant silicon sources from soil and/or via selection for traits that aid the maintenance of a high silicon status in shoots. Recent advances in understanding the genetic regulation of silicon uptake and transport by rice plants will aid these goals.

"C. elegans: modelo biológico do presente e do futuro"

Nos meados da década de 60, o Nobel da Medicina Sidney Brenner propôs a utilização do nemátode bacteriófago, do solo, Caenorhabditis elegans, também conhecido como C. elegans, para estudos de genética e desenvolvimento, dado possuir um conjunto de características que o tornavam o ideal para modelo biológico, nomeadamente a sua pequena dimensão (ca. de 1 mm), facilidade de observação, facilidade de cultivo e manutenção em laboratório, curto ciclo de vida com uma capacidade reprodutiva notável, presença de hermafroditas e machos (5%). O seu artigo seminal de 1974, “The genetics of C. elegans” abriu o caminho para a investigação nesta área. Hoje em dia, revistas como a Nature, Science, Genes and Development, etc… publicam frequentemente os resultados da investigação com recurso a este modelo. É de notar, no entanto, que os nemátodes têm sido utilizados há muito tempo como modelo de estudo, tais como van Beneden, em finais do século XIX, que observando células de Ascaris equorum, descobriu o fenómeno da meiose. O fenómeno da fertilização foi igualmente descoberto num nemátode. Desde a década de setenta até aos nossos dias, o C. elegans tem sido intensamente utilizado para estudos de anatomia interna e sua correlação com linhagens celulares e desenvolvimento. Assim, Sulston e Horvitz elucidaram a origem e desenvolvimento das 959 células somáticas que, de uma forma constante, se produzem nesta espécie (eutelia). Um dos primeiros sistemas a ser estudado foi o sistema nervoso, sendo este nemátode o primeiro animal de que se conhece perfeitamente a identidade de todos os neurónios, a sua linhagem e o circuito nervoso global. Para além de modelo de biologia do desenvolvimento, o C. elegans tem sido alvo de estudo do fenómeno do envelhecimento celular, tendo sido possível identificar os respectivos genes. Kennyon, nos anos 90, e mais recentemente Arantes e Oliveira, têm demonstrado ser possível “prolongar” a vida deste animal de 21 para mais de 180 dias, o que corresponde a 675 anos em vida humana, através de manipulações diversas, tais como agentes mutagénicos, ablação com raio laser, etc.. Também o mecanismo de morte celular programada ou “apoptose” tem sido estudado neste modelo. Em 2002, o Comité Nobel entendeu atribuir o prémio Nobel da Medicina a Brenner, Sulston e Horvitz “for their discoveries concerning genetic regulation of organ development and programmed cell death”. É interessante notar pela leitura de “The common thread”, de John Sulston, a importância decisiva da sequenciação do genoma de C. elegans, em 1998, para o avanço na sequenciação do genoma humano efectuada em 2000, e de cuja mega-equipa Sulston participou de forma decisiva. Finalmente, e como modelo pedagógico, o nemátode C. elegans constitui a escolha ideal para diversas disciplinas dos cursos de Biologia, tais como Biologia celular, Histologia, Biologia do Desenvolvimento, Genética, Etologia, etc….

Étude fonctionnelle et génétique d'une population de lymphocytes T CD4-CD8- impliquée dans la résistance au diabète auto-immun chez la souris

Le diabète auto-immun résulte de la destruction des cellules bêta pancréatiques sécrétrices d’insuline par les lymphocytes T du système immunitaire. Il s’ensuit une déficience hormonale qui peut être comblée par des injections quotidiennes d’insuline d’origine exogène, toutefois il demeure à ce jour impossible de guérir les patients atteints de la maladie. De façon générale, un système immunitaire sain reconnaît une multitude d’antigènes différents et assure ainsi notre défense à l’égard de différents pathogènes ou encore de cellules tumorales. Il arrive cependant que, pour des raisons génétiques et/ou environnementales, les lymphocytes T puissent s’activer de façon aberrante suite à la reconnaissance d’antigènes provenant du soi. C’est ce bris de tolérance qui mène au développement de pathologies auto-immunes telles que le diabète auto-immun. Afin de limiter l’auto-immunité, des mécanismes de sélection stricts permettent d’éliminer la majorité des lymphocytes T présentant une forte affinité envers des antigènes du soi lors de leur développement dans le thymus. Certains de ces lymphocytes réussissent toutefois à échapper à l’apoptose et migrent en périphérie afin d’y circuler en quête d’un antigène spécifiquement reconnu. Il est alors primordial que des mécanismes périphériques assurent le maintien de la tolérance immunitaire en faisant obstacle à l’activation et à la prolifération des lymphocytes T auto-réactifs. L’une des avenues afin d’inhiber le développement de réponses immunitaires aberrantes est la génération de lymphocytes T régulateurs. Ces cellules, d’origine thymique ou périphérique, peuvent arborer différents phénotypes et agissent via de multiples mécanismes afin d’inactiver et/ou éliminer les cellules impliquées dans l’apparition de pathologies auto-immunes. L’utilisation de modèles murins transgéniques a permis la mise en évidence d’une population peu caractérisée de lymphocytes T au potentiel régulateur. En effet, la proportion de ces cellules T n’exprimant pas les corécepteurs CD4 et CD8 (double négatives, DN) a été inversement corrélée à la prédisposition à l’auto-immunité chez ces ii souris. L’objectif principal de cette thèse est de démontrer la fonction immuno-régulatrice des lymphocytes T DN, tout en investiguant les facteurs génétiques responsables du maintien de cette population cellulaire. Nous avons observé que les lymphocytes T DN exercent une activité cytotoxique à l’égard des lymphocytes B de façon spécifique à l’antigène, via la libération de granules cytolytiques contenant du granzyme B et de la perforine. Par ailleurs, nous avons établi qu’un unique transfert adoptif de ces cellules est suffisant afin d’inhiber le développement du diabète auto-immun chez des hôtes transgéniques prédisposés à la maladie. Le recours à des souris déficientes pour l’expression du gène CD47 a permis de constater que la voie de signalisation CD47-Sirp est essentielle dans le maintien de la proportion des lymphocytes T DN. De plus, le locus murin de prédisposition au diabète auto-immun Idd13, qui contient le gène Sirp, a été identifié pour son rôle dans la régulation de la proportion de ces cellules. Finalement, une analyse génétique a révélé que d’autres intervalles génétiques sont impliqués dans le contrôle de la population des lymphocytes T DN. Parmi ceux-ci, un locus situé en région proximale du chromosome 12 a été validé grâce à la création de souris congéniques. Grâce aux résultats présentés dans cette thèse, notre compréhension de la biologie ainsi que de la régulation des lymphocytes T DN est approfondie. Ces connaissances constituent un pas important vers la création de thérapies cellulaires novatrices permettant de prévenir et de guérir diverses pathologies auto-immunes.

La cartographie des sites de régulation génétique à partir de données de débalancement allélique

En 1975, Wilson et King ont proposé que l'évolution opère non seulement via des changements affectant la structure des protéines, mais aussi via des mutations qui modifient la régulation génétique. L'étude des éléments régulateurs de l'expression génétique a un rôle important dans la compréhension de l'expression de différentes maladies et de la réponse thérapeutique. Nous avons développé un algorithme bio- informatique qui nous permet rapidement de trouver des sites de régulation génétique à travers tout le génome et pour une grande quantité de gènes. Notre approche consiste à trouver des sites polymorphes (SNPs) qui sont en déséquilibre de liaison avec le débalancement allélique (AI) afin de cartographier la région régulatrice et le site responsable. Notre méthode est avantageuse par rapport à d'autres méthodes, car elle n'a pas besoin des données « phasées». De plus, les données de débalancement allélique ne sont pas affectées par des facteurs externes étant donné qu'ils sont mesurés dans la même cellule. Nous avons démontré que notre approche est fiable et qu'elle peut détecter des sites loin du gène. De plus, il peut être appliqué à des données de génotypage sans avoir besoin de les « phaser » .

Molecular genetics of cork formation

La peridermis és una estructura complexa que protegeix els òrgans vegetals madurs (secundaris) i les zones que han sofert ferides de la pèrdua d'aigua i dels patògens. Aquesta funció barrera és deguda al fel·lema o súber, un teixit format per cèl·lules suberificades. Tant el fel·lema com la suberina són crucials per la vida de les plantes terrestres, però pràcticament no es coneix res dels processos moleculars que regulen la seva formació, probablement degut a la manca de models adequats. En aquesta tesi s'han identificat i caracteritzat gens induïts al fel·lema mitjançant la combinació de dues plantes models. L'escorça d'alzina surera (Quercus suber) s'ha utilitzat per aïllar gens candidats de la formació del fel·lema i per investigar el comportament d'alguns d'aquests gens durant l'estació de creixement, mentre que la pela de la patata (Solanum tuberosum) s'ha utilitzat en estudis de genètica reversa per demostrar la funció d'alguns gens reguladors al fel·lema.

Physiological changes in Campylobacter jejuni on entry into stationary phase

Campylobacter jejuni NCTC 11168 does not exhibit the general increase in cellular stress resistance on entry into stationary phase that is seen in most other bacteria. This is consistent with the lack of global stationary phase regulatory elements in this organism. deduced from an analysis of its genome sequence. We now show that C. jejuni NCTC 11168 does undergo certain changes in stationary phase, of a pattern not previously described. As cells entered stationary phase there was a change in membrane fatty acid composition, principally a decrease in the proportion of unsaturated fatty acids and an increase in the content of cyclopropane and short-chain fatty acids. These changes in membrane composition were accompanied by an increase in the resilience of the cell membrane towards loss of integrity caused by pressure and an increase in cellular pressure resistance. By contrast. there were no major changes in resistance to acid or heat treatment. A similar pattern of changes in stress resistance on entry, into stationary phase was seen in C. jejuni NCTC 11351, the type strain. These changes appear to represent a restricted Physiological response to the conditions existing in stationary phase cultures, in an organism having limited capacity for genetic regulation and adaptation to environment. © 2004 Elsevier B.V. All rights reserved.

Variation in rice cadmium related to human exposure.

Cereal grains are the dominant source of cadmium in the human diet, with rice being to the fore. Here we explore the effect of geographic, genetic, and processing (milling) factors on rice grain cadmium and rice consumption rates that lead to dietary variance in cadmium intake. From a survey of 12 countries on four continents, cadmium levels in rice grain were the highest in Bangladesh and Sri Lanka, with both these countries also having high per capita rice intakes. For Bangladesh and Sri Lanka, there was high weekly intake of cadmium from rice, leading to intakes deemed unsafe by international and national regulators. While genetic variance, and to a lesser extent milling, provide strategies for reducing cadmium in rice, caution has to be used, as there is environmental regulation as well as genetic regulation of cadmium accumulation within rice grains. For countries that import rice, grain cadmium can be controlled by where that rice is sourced, but for countries with subsistence rice economies that have high levels of cadmium in rice grain, agronomic and breeding strategies are required to lower grain cadmium.

A mathematical model of the sterol regulatory element binding protein 2 cholesterol biosynthesis pathway

Cholesterol is one of the key constituents for maintaining the cellular membrane and thus the integrity of the cell itself. In contrast high levels of cholesterol in the blood are known to be a major risk factor in the development of cardiovascular disease. We formulate a deterministic nonlinear ordinary differential equation model of the sterol regulatory element binding protein 2 (SREBP-2) cholesterol genetic regulatory pathway in an hepatocyte. The mathematical model includes a description of genetic transcription by SREBP-2 which is subsequently translated to mRNA leading to the formation of 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGCR), a main precursor of cholesterol synthesis. Cholesterol synthesis subsequently leads to the regulation of SREBP-2 via a negative feedback formulation. Parameterised with data from the literature, the model is used to understand how SREBP-2 transcription and regulation affects cellular cholesterol concentration. Model stability analysis shows that the only positive steady-state of the system exhibits purely oscillatory, damped oscillatory or monotic behaviour under certain parameter conditions. In light of our findings we postulate how cholesterol homestasis is maintained within the cell and the advantages of our model formulation are discussed with respect to other models of genetic regulation within the literature.