The challenge of modeling nuclear receptor regulatory networks in mammalian cells.

Nuclear receptors are a major component of signal transduction in animals. They mediate the regulatory activities of many hormones, nutrients and metabolites on the homeostasis and physiology of cells and tissues. It is of high interest to model the corresponding regulatory networks. While molecular and cell biology studies of individual promoters have provided important mechanistic insight, a more complex picture is emerging from genome-wide studies. The regulatory circuitry of nuclear receptor regulated gene expression networks, and their response to cellular signaling, appear highly dynamic, and involve long as well as short range chromatin interactions. We review how progress in understanding the kinetics and regulation of cofactor recruitment, and the development of new genomic methods, provide opportunities but also a major challenge for modeling nuclear receptor mediated regulatory networks.

A limited role for beta-selection during gamma delta T cell development.

T cells belong to two distinct lineages expressing either alpha beta or gamma delta TCR. During alpha beta T cell development, it is clearly established that productive rearrangement at the TCR beta locus in immature precursor cells leads to the expression of a pre-TCR complex. Signaling through the pre-TCR results in the selective proliferation and maturation of TCR beta+ cells, a process that is known as beta-selection. However, the potential role of beta-selection during gamma delta T cell development is controversial. Whereas PCR-RFLP and sequencing techniques have provided evidence for a bias toward in-frame VDJ beta rearrangements in gamma delta cells (consistent with beta-selection), gamma delta cells apparently develop normally in mice that are unable to assemble a pre-TCR complex due to a deficiency in TCR beta or pT alpha genes. In this report, we have directly addressed the physiologic significance of beta-selection during gamma delta cell development in normal mice by quantitating intracellular TCR beta protein in gamma delta cells and correlating its presence with cell cycle status. Our results indicate that beta-selection plays a significant (although limited) role in gamma delta cell development by selectively amplifying a minor subset of gamma delta precursor cells with productively rearranged TCR beta genes.

Towards robust network based complex systems : from evolutionary cellular automata to biological models

Abstract Sitting between your past and your future doesn't mean you are in the present. Dakota Skye Complex systems science is an interdisciplinary field grouping under the same umbrella dynamical phenomena from social, natural or mathematical sciences. The emergence of a higher order organization or behavior, transcending that expected of the linear addition of the parts, is a key factor shared by all these systems. Most complex systems can be modeled as networks that represent the interactions amongst the system's components. In addition to the actual nature of the part's interactions, the intrinsic topological structure of underlying network is believed to play a crucial role in the remarkable emergent behaviors exhibited by the systems. Moreover, the topology is also a key a factor to explain the extraordinary flexibility and resilience to perturbations when applied to transmission and diffusion phenomena. In this work, we study the effect of different network structures on the performance and on the fault tolerance of systems in two different contexts. In the first part, we study cellular automata, which are a simple paradigm for distributed computation. Cellular automata are made of basic Boolean computational units, the cells; relying on simple rules and information from- the surrounding cells to perform a global task. The limited visibility of the cells can be modeled as a network, where interactions amongst cells are governed by an underlying structure, usually a regular one. In order to increase the performance of cellular automata, we chose to change its topology. We applied computational principles inspired by Darwinian evolution, called evolutionary algorithms, to alter the system's topological structure starting from either a regular or a random one. The outcome is remarkable, as the resulting topologies find themselves sharing properties of both regular and random network, and display similitudes Watts-Strogtz's small-world network found in social systems. Moreover, the performance and tolerance to probabilistic faults of our small-world like cellular automata surpasses that of regular ones. In the second part, we use the context of biological genetic regulatory networks and, in particular, Kauffman's random Boolean networks model. In some ways, this model is close to cellular automata, although is not expected to perform any task. Instead, it simulates the time-evolution of genetic regulation within living organisms under strict conditions. The original model, though very attractive by it's simplicity, suffered from important shortcomings unveiled by the recent advances in genetics and biology. We propose to use these new discoveries to improve the original model. Firstly, we have used artificial topologies believed to be closer to that of gene regulatory networks. We have also studied actual biological organisms, and used parts of their genetic regulatory networks in our models. Secondly, we have addressed the improbable full synchronicity of the event taking place on. Boolean networks and proposed a more biologically plausible cascading scheme. Finally, we tackled the actual Boolean functions of the model, i.e. the specifics of how genes activate according to the activity of upstream genes, and presented a new update function that takes into account the actual promoting and repressing effects of one gene on another. Our improved models demonstrate the expected, biologically sound, behavior of previous GRN model, yet with superior resistance to perturbations. We believe they are one step closer to the biological reality.

Androgenetic alopecia: identification of four genetic risk loci and evidence for the contribution of WNT signaling to its etiology.

The pathogenesis of androgenetic alopecia (AGA, male-pattern baldness) is driven by androgens, and genetic predisposition is the major prerequisite. Candidate gene and genome-wide association studies have reported that single-nucleotide polymorphisms (SNPs) at eight different genomic loci are associated with AGA development. However, a significant fraction of the overall heritable risk still awaits identification. Furthermore, the understanding of the pathophysiology of AGA is incomplete, and each newly associated locus may provide novel insights into contributing biological pathways. The aim of this study was to identify unknown AGA risk loci by replicating SNPs at the 12 genomic loci that showed suggestive association (5 × 10(-8)<P<10(-5)) with AGA in a recent meta-analysis. We analyzed a replication set comprising 2,759 cases and 2,661 controls of European descent to confirm the association with AGA at these loci. Combined analysis of the replication and the meta-analysis data identified four genome-wide significant risk loci for AGA on chromosomes 2q35, 3q25.1, 5q33.3, and 12p12.1. The strongest association signal was obtained for rs7349332 (P=3.55 × 10(-15)) on chr2q35, which is located intronically in WNT10A. Expression studies in human hair follicle tissue suggest that WNT10A has a functional role in AGA etiology. Thus, our study provides genetic evidence supporting an involvement of WNT signaling in AGA development.

A coiled coil trigger site is essential for rapid binding of synaptobrevin to the SNARE acceptor complex.

Exocytosis from synaptic vesicles is driven by stepwise formation of a tight alpha-helical complex between the fusing membranes. The complex is composed of the three SNAREs: synaptobrevin 2, SNAP-25, and syntaxin 1a. An important step in complex formation is fast binding of vesicular synaptobrevin to the preformed syntaxin 1.SNAP-25 dimer. Exactly how this step relates to neurotransmitter release is not well understood. Here, we combined different approaches to gain insights into this reaction. Using computational methods, we identified a stretch in synaptobrevin 2 that may function as a coiled coil "trigger site." This site is also present in many synaptobrevin homologs functioning in other trafficking steps. Point mutations in this stretch inhibited binding to the syntaxin 1.SNAP-25 dimer and slowed fusion of liposomes. Moreover, the point mutations severely inhibited secretion from chromaffin cells. Altogether, this demonstrates that the trigger site in synaptobrevin is crucial for productive SNARE zippering.

Low density lipoprotein signaling and pancreatic beta cells protection

SummaryLow-density lipoproteins (LDLs) have an important physiological role in organism transporting cholesterol and other fatty substances to target tissues. However, elevated LDL levels in the blood are associated with the formation of arterial plaques and consequently atherosclerosis. It is therefore important to characterize the intracellular pathways induced upon LDL stimulation as they might be involved in the pathological properties of these lipoproteins. It has been previously found that LDL stimulation of mouse embryonic fibroblasts activates p38 mitogen activated protein kinases (MAPKs). This leads to cell spreading and increase in the wound healing capabilities of the cells. These two responses might occur within atherosclerotic plaques.The aim of this project is to reveal the missing links between LDL particle and activation of p38 MAPK kinase. As previously shown in our lab activation of p38 MAPK kinase by the LDL particles occur independently of classical LDL receptor (LDLR). In this study we have shown that scavenger receptor type Β class I (SR-BI) is responsible for the signal transduction from the LDLs to the p38 MAPK. We have also shown that Mitogen activated kinase kinases (MKKs) that can directly activate ρ 38 MAPK in these conditions are MKK3 and MKK6 but not MKK4. We have also tested some of the intermediate components of the pathway like Ras and PI3 kinase but found that they do not play a role.The data obtained in this study showed a part of molecular mechanism responsible for p38 MAPK activation and subsequent wound healing and can contribute to our knowledge on function of the fibroblasts in the development of the atherosclerotic plaques.Diabetes Mellitus is a condition caused by disordered metabolism of blood glucose level. It is one of the most commonly spread disease in the western world, with the incidence reaching 8% of population in United States. Two most common types of diabetes are type 1 and 2 that differs slightly in the mechanism of the development. However in the basis of both types lies the cell death of pancreatic beta cells. The aim of this work is to improve beta cells survival in different pathophysiological settings. This could be extrapolated to the conditions in which Diabetes develops in humans. We decided to use RasGAP- derived fragment Ν with its strong antiapoptotic effect in beta cells. In our lab we have demonstrated that in the mild stress conditions RasGAP can be cleaved by caspases at the position 455 producing two fragments, fragment Ν and fragment C. Fragment Ν exerts

n-3 PUFAs modulate T-cell activation via protein kinase C-alpha and -epsilon and the NF-kappaB signaling pathway.

We elucidated the mechanisms of action of two n-3 PUFAs, eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA), in Jurkat T-cells. Both DHA and EPA were principally incorporated into phospholipids in the following order: phosphatidylcholine < phosphatidylethanolamine < phosphatidylinositol/phosphatidylserine. Furthermore, two isoforms of phospholipase A(2) (i.e., calcium-dependent and calcium-independent) were implicated in the release of DHA and EPA, respectively, during activation of these cells. The two fatty acids inhibited the phorbol 12-myristate 13-acetate (PMA)-induced plasma membrane translocation of protein kinase C (PKC)-alpha and -epsilon. The two n-3 PUFAs also inhibited the nuclear translocation of nuclear factor kappaB (NF-kappaB) and the transcription of the interleukin-2 (IL-2) gene in PMA-activated Jurkat T-cells. Together, these results demonstrate that DHA and EPA, being released by two isoforms of phospholipase A(2), modulate IL-2 gene expression by exerting their action on two PKC isoforms and NF-kappaB in Jurkat T-cells.

Pivotal role for a tail subunit of the RNA polymerase II mediator complex CgMed2 in azole tolerance and adherence in Candida glabrata.

Antifungal therapy failure can be associated with increased resistance to the employed antifungal agents. Candida glabrata, the second most common cause of invasive candidiasis, is intrinsically less susceptible to the azole class of antifungals and accounts for 15% of all Candida bloodstream infections. Here, we show that C. glabrata MED2 (CgMED2), which codes for a tail subunit of the RNA polymerase II Mediator complex, is required for resistance to azole antifungal drugs in C. glabrata. An inability to transcriptionally activate genes encoding a zinc finger transcriptional factor, CgPdr1, and multidrug efflux pump, CgCdr1, primarily contributes to the elevated susceptibility of the Cgmed2Δ mutant toward azole antifungals. We also report for the first time that the Cgmed2Δ mutant exhibits sensitivity to caspofungin, a constitutively activated protein kinase C-mediated cell wall integrity pathway, and elevated adherence to epithelial cells. The increased adherence of the Cgmed2Δ mutant was attributed to the elevated expression of the EPA1 and EPA7 genes. Further, our data demonstrate that CgMED2 is required for intracellular proliferation in human macrophages and modulates survival in a murine model of disseminated candidiasis. Lastly, we show an essential requirement for CgMed2, along with the Mediator middle subunit CgNut1 and the Mediator cyclin-dependent kinase/cyclin subunit CgSrb8, for the high-level fluconazole resistance conferred by the hyperactive allele of CgPdr1. Together, our findings underscore a pivotal role for CgMed2 in basal tolerance and acquired resistance to azole antifungals.

Error and attack tolerance of layered complex networks.

Many complex systems may be described by not one but a number of complex networks mapped on each other in a multi-layer structure. Because of the interactions and dependencies between these layers, the state of a single layer does not necessarily reflect well the state of the entire system. In this paper we study the robustness of five examples of two-layer complex systems: three real-life data sets in the fields of communication (the Internet), transportation (the European railway system), and biology (the human brain), and two models based on random graphs. In order to cover the whole range of features specific to these systems, we focus on two extreme policies of system's response to failures, no rerouting and full rerouting. Our main finding is that multi-layer systems are much more vulnerable to errors and intentional attacks than they appear from a single layer perspective.

Evolutionary games on complex networks : the emergence and maintenance of cooperation

Abstract The object of game theory lies in the analysis of situations where different social actors have conflicting requirements and where their individual decisions will all influence the global outcome. In this framework, several games have been invented to capture the essence of various dilemmas encountered in many common important socio-economic situations. Even though these games often succeed in helping us understand human or animal behavior in interactive settings, some experiments have shown that people tend to cooperate with each other in situations for which classical game theory strongly recommends them to do the exact opposite. Several mechanisms have been invoked to try to explain the emergence of this unexpected cooperative attitude. Among them, repeated interaction, reputation, and belonging to a recognizable group have often been mentioned. However, the work of Nowak and May (1992) showed that the simple fact of arranging the players according to a spatial structure and only allowing them to interact with their immediate neighbors is sufficient to sustain a certain amount of cooperation even when the game is played anonymously and without repetition. Nowak and May's study and much of the following work was based on regular structures such as two-dimensional grids. Axelrod et al. (2002) showed that by randomizing the choice of neighbors, i.e. by actually giving up a strictly local geographical structure, cooperation can still emerge, provided that the interaction patterns remain stable in time. This is a first step towards a social network structure. However, following pioneering work by sociologists in the sixties such as that of Milgram (1967), in the last few years it has become apparent that many social and biological interaction networks, and even some technological networks, have particular, and partly unexpected, properties that set them apart from regular or random graphs. Among other things, they usually display broad degree distributions, and show small-world topological structure. Roughly speaking, a small-world graph is a network where any individual is relatively close, in terms of social ties, to any other individual, a property also found in random graphs but not in regular lattices. However, in contrast with random graphs, small-world networks also have a certain amount of local structure, as measured, for instance, by a quantity called the clustering coefficient. In the same vein, many real conflicting situations in economy and sociology are not well described neither by a fixed geographical position of the individuals in a regular lattice, nor by a random graph. Furthermore, it is a known fact that network structure can highly influence dynamical phenomena such as the way diseases spread across a population and ideas or information get transmitted. Therefore, in the last decade, research attention has naturally shifted from random and regular graphs towards better models of social interaction structures. The primary goal of this work is to discover whether or not the underlying graph structure of real social networks could give explanations as to why one finds higher levels of cooperation in populations of human beings or animals than what is prescribed by classical game theory. To meet this objective, I start by thoroughly studying a real scientific coauthorship network and showing how it differs from biological or technological networks using divers statistical measurements. Furthermore, I extract and describe its community structure taking into account the intensity of a collaboration. Finally, I investigate the temporal evolution of the network, from its inception to its state at the time of the study in 2006, suggesting also an effective view of it as opposed to a historical one. Thereafter, I combine evolutionary game theory with several network models along with the studied coauthorship network in order to highlight which specific network properties foster cooperation and shed some light on the various mechanisms responsible for the maintenance of this same cooperation. I point out the fact that, to resist defection, cooperators take advantage, whenever possible, of the degree-heterogeneity of social networks and their underlying community structure. Finally, I show that cooperation level and stability depend not only on the game played, but also on the evolutionary dynamic rules used and the individual payoff calculations. Synopsis Le but de la théorie des jeux réside dans l'analyse de situations dans lesquelles différents acteurs sociaux, avec des objectifs souvent conflictuels, doivent individuellement prendre des décisions qui influenceront toutes le résultat global. Dans ce cadre, plusieurs jeux ont été inventés afin de saisir l'essence de divers dilemmes rencontrés dans d'importantes situations socio-économiques. Bien que ces jeux nous permettent souvent de comprendre le comportement d'êtres humains ou d'animaux en interactions, des expériences ont montré que les individus ont parfois tendance à coopérer dans des situations pour lesquelles la théorie classique des jeux prescrit de faire le contraire. Plusieurs mécanismes ont été invoqués pour tenter d'expliquer l'émergence de ce comportement coopératif inattendu. Parmi ceux-ci, la répétition des interactions, la réputation ou encore l'appartenance à des groupes reconnaissables ont souvent été mentionnés. Toutefois, les travaux de Nowak et May (1992) ont montré que le simple fait de disposer les joueurs selon une structure spatiale en leur permettant d'interagir uniquement avec leurs voisins directs est suffisant pour maintenir un certain niveau de coopération même si le jeu est joué de manière anonyme et sans répétitions. L'étude de Nowak et May, ainsi qu'un nombre substantiel de travaux qui ont suivi, étaient basés sur des structures régulières telles que des grilles à deux dimensions. Axelrod et al. (2002) ont montré qu'en randomisant le choix des voisins, i.e. en abandonnant une localisation géographique stricte, la coopération peut malgré tout émerger, pour autant que les schémas d'interactions restent stables au cours du temps. Ceci est un premier pas en direction d'une structure de réseau social. Toutefois, suite aux travaux précurseurs de sociologues des années soixante, tels que ceux de Milgram (1967), il est devenu clair ces dernières années qu'une grande partie des réseaux d'interactions sociaux et biologiques, et même quelques réseaux technologiques, possèdent des propriétés particulières, et partiellement inattendues, qui les distinguent de graphes réguliers ou aléatoires. Entre autres, ils affichent en général une distribution du degré relativement large ainsi qu'une structure de "petit-monde". Grossièrement parlant, un graphe "petit-monde" est un réseau où tout individu se trouve relativement près de tout autre individu en termes de distance sociale, une propriété également présente dans les graphes aléatoires mais absente des grilles régulières. Par contre, les réseaux "petit-monde" ont, contrairement aux graphes aléatoires, une certaine structure de localité, mesurée par exemple par une quantité appelée le "coefficient de clustering". Dans le même esprit, plusieurs situations réelles de conflit en économie et sociologie ne sont pas bien décrites ni par des positions géographiquement fixes des individus en grilles régulières, ni par des graphes aléatoires. De plus, il est bien connu que la structure même d'un réseau peut passablement influencer des phénomènes dynamiques tels que la manière qu'a une maladie de se répandre à travers une population, ou encore la façon dont des idées ou une information s'y propagent. Ainsi, durant cette dernière décennie, l'attention de la recherche s'est tout naturellement déplacée des graphes aléatoires et réguliers vers de meilleurs modèles de structure d'interactions sociales. L'objectif principal de ce travail est de découvrir si la structure sous-jacente de graphe de vrais réseaux sociaux peut fournir des explications quant aux raisons pour lesquelles on trouve, chez certains groupes d'êtres humains ou d'animaux, des niveaux de coopération supérieurs à ce qui est prescrit par la théorie classique des jeux. Dans l'optique d'atteindre ce but, je commence par étudier un véritable réseau de collaborations scientifiques et, en utilisant diverses mesures statistiques, je mets en évidence la manière dont il diffère de réseaux biologiques ou technologiques. De plus, j'extrais et je décris sa structure de communautés en tenant compte de l'intensité d'une collaboration. Finalement, j'examine l'évolution temporelle du réseau depuis son origine jusqu'à son état en 2006, date à laquelle l'étude a été effectuée, en suggérant également une vue effective du réseau par opposition à une vue historique. Par la suite, je combine la théorie évolutionnaire des jeux avec des réseaux comprenant plusieurs modèles et le réseau de collaboration susmentionné, afin de déterminer les propriétés structurelles utiles à la promotion de la coopération et les mécanismes responsables du maintien de celle-ci. Je mets en évidence le fait que, pour ne pas succomber à la défection, les coopérateurs exploitent dans la mesure du possible l'hétérogénéité des réseaux sociaux en termes de degré ainsi que la structure de communautés sous-jacente de ces mêmes réseaux. Finalement, je montre que le niveau de coopération et sa stabilité dépendent non seulement du jeu joué, mais aussi des règles de la dynamique évolutionnaire utilisées et du calcul du bénéfice d'un individu.

Volatile transport during the crystallization of anatectic melts: Oxygen, boron and hydrogen stable isotope study on the metamorphic complex of Naxos, Greece

Crystallization of anatectic melts in high-temperature metamorphic terrains releases volatile-rich magmas that can be transported into adjacent lithologies. This study addresses the variations in the oxygen, boron and hydrogen isotopic composition of aplite-pegmatite dikes that formed during the crystallization of anatectic melts in regional high-temperature metamorphism on the island of Naxos, Greece, and propagated upward into the overlying sequences of metamorphic schist. The transport distance of these dikes was increased through a significant horizontal component of travel that was imposed by contemporaneous low-angle extensional shearing. Laser fluorination oxygen isotope analyses of quartz, tourmaline, garnet, and biotite mineral separates from the aplite-pegmatite dikes show a progressive rise in delta(18)O values with increasing distance from the core. Oxygen isotope fractionations among quartz, tourmaline, and garnet show temperature variations from > 700degreesC down to similar to400degreesC. This range is considered to reflect isotopic fractionation beginning with crystallization at high temperatures in water-undersaturated conditions and then evolving through lower temperature crystallization and retrograde sub-solidus exchange. Two processes are examined for the cause of the progressive increase in delta(18)O values: (1) heterogeneous delta(18)O sources and (2) fluid-rock exchange between the aplite/pegmatite magmas and their host rock. Although the former process cannot be ruled out, there is as yet no evidence in the exposed sequences on Naxos for the presence of a suitable high delta(18)O magma source. In contrast, a tendency for the delta(18)O of quartz in the aplite/pegmatite dikes to approach that of the quartz in the metamorphic rock suggests that fluid-rock exchange with the host rock may potentially be an important process. Advection of fluid into the magma is examined based on Darcian fluid flow into an initially water-undersaturated buoyantly propagating aplitic dike magma. It is shown that such advective flow could only account for part of the O-18-enrichment, unless it were amplified by repeated injection of magma pulses, fluid recycling, and deformation-assisted post-crystallization exchange. The mechanism is, however, adequate to account for hydrogen isotope equilibration between dike and host rock. In contrast, variations in the delta(11)B values of tourmalines suggest that B-11/B-10 fractionation during crystallization and/or magma degassing was the major control of boron geochemistry rather than fluid-rock interaction and that the boron isotopic system was decoupled from that of oxygen. Copyright (C) 2003 Elsevier Ltd.

Pharmacological inhibition of fibroblast growth factor (FGF) receptor signaling ameliorates FGF23-mediated hypophosphatemic rickets.

Fibroblast growth factor 23 (FGF23) is a circulating factor secreted by osteocytes that is essential for phosphate homeostasis. In kidney proximal tubular cells FGF23 inhibits phosphate reabsorption and leads to decreased synthesis and enhanced catabolism of 1,25-dihydroxyvitamin D3 (1,25[OH]2 D3 ). Excess levels of FGF23 cause renal phosphate wasting and suppression of circulating 1,25(OH)2 D3 levels and are associated with several hereditary hypophosphatemic disorders with skeletal abnormalities, including X-linked hypophosphatemic rickets (XLH) and autosomal recessive hypophosphatemic rickets (ARHR). Currently, therapeutic approaches to these diseases are limited to treatment with activated vitamin D analogues and phosphate supplementation, often merely resulting in partial correction of the skeletal aberrations. In this study, we evaluate the use of FGFR inhibitors for the treatment of FGF23-mediated hypophosphatemic disorders using NVP-BGJ398, a novel selective, pan-specific FGFR inhibitor currently in Phase I clinical trials for cancer therapy. In two different hypophosphatemic mouse models, Hyp and Dmp1-null mice, resembling the human diseases XLH and ARHR, we find that pharmacological inhibition of FGFRs efficiently abrogates aberrant FGF23 signaling and normalizes the hypophosphatemic and hypocalcemic conditions of these mice. Correspondingly, long-term FGFR inhibition in Hyp mice leads to enhanced bone growth, increased mineralization, and reorganization of the disturbed growth plate structure. We therefore propose NVP-BGJ398 treatment as a novel approach for the therapy of FGF23-mediated hypophosphatemic diseases.

Wound repair and regeneration.

The repair of wounds is one of the most complex biological processes that occur during human life. After an injury, multiple biological pathways immediately become activated and are synchronized to respond. In human adults, the wound repair process commonly leads to a non-functioning mass of fibrotic tissue known as a scar. By contrast, early in gestation, injured fetal tissues can be completely recreated, without fibrosis, in a process resembling regeneration. Some organisms, however, retain the ability to regenerate tissue throughout adult life. Knowledge gained from studying such organisms might help to unlock latent regenerative pathways in humans, which would change medical practice as much as the introduction of antibiotics did in the twentieth century.

Nd-, O-, and H-isotopic evidence for complex, closed-system fluid evolution of the peralkaline Ilimaussaq intrusion, South Greenland

Relatively homogeneous oxygen isotope compositions of amphibole, clinopyroxene, and olivine separates (+5.2 to +5.7parts per thousand relative to VSMOW) and neodymium isotope compositions (epsilon(Nd(T)) = -0.9 to -1.8 for primary magmatic minerals and epsilon(Nd(T)) = -0.1 and -0.5 for mineral separates from late-stage pegmatites and hydrothermal veins) from the alkaline to agpaitic llimaussaq intrusion, South Greenland, indicate a closed system evolution of this igneous complex and support a mantle derivation of the magma. In contrast to the homogeneous oxygen and neodymium isotopic data, deltaD values for hand-picked amphibole separates vary between -92 and -232parts per thousand and are among the most deuterium-depleted values known from igneous amphiboles. The calculated fluid phase coexisting with these amphiboles has a homogeneous oxygen isotopic composition within the normal range of magmatic waters, but extremely heterogeneous and low D/H ratios, implying a decoupling of the oxygen- and hydrogen isotope systems. Of the several possibilities that can account for such unusually low deltaD values in amphiboles (e.g., late-stage hydrothermal exchange with meteoric water, extensive magmatic degassing, contamination with organic matter, and/or effects of Fe-content and pressure on amphibole-water fractionation) the most likely explanation for the range in deltaD values is that the amphiboles have been influenced by secondary interaction and reequilibration with D-depleted fluids obtained through late-magmatic oxidation of internally generated CH(4) and/or H(2). This interpretation is consistent with the known occurrence of abundant magmatic CH(4) in the Ilimaussaq rocks and with previous studies on the isotopic compositions of the rocks and fluids. Copyright (C) 2004 Elsevier Ltd.