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.

Flippase (FLP) recombinase-mediated marker recycling in the dermatophyte Arthroderma vanbreuseghemii.

Biological processes can be elucidated by investigating complex networks of relevant factors and genes. However, this is not possible in species for which dominant selectable markers for genetic studies are unavailable. To overcome the limitation in selectable markers for the dermatophyte Arthroderma vanbreuseghemii (anamorph: Trichophyton mentagrophytes), we adapted the flippase (FLP) recombinase-recombination target (FRT) site-specific recombination system from the yeast Saccharomyces cerevisiae as a selectable marker recycling system for this fungus. Taking into account practical applicability, we designed FLP/FRT modules carrying two FRT sequences as well as the flp gene adapted to the pathogenic yeast Candida albicans (caflp) or a synthetic codon-optimized flp (avflp) gene with neomycin resistance (nptII) cassette for one-step marker excision. Both flp genes were under control of the Trichophyton rubrum copper-repressible promoter (PCTR4). Molecular analyses of resultant transformants showed that only the avflp-harbouring module was functional in A. vanbreuseghemii. Applying this system, we successfully produced the Ku80 recessive mutant strain devoid of any selectable markers. This strain was subsequently used as the recipient for sequential multiple disruptions of secreted metalloprotease (fungalysin) (MEP) or serine protease (SUB) genes, producing mutant strains with double MEP or triple SUB gene deletions. These results confirmed the feasibility of this system for broad-scale genetic manipulation of dermatophytes, advancing our understanding of functions and networks of individual genes in these fungi.

Influence of the ankle arthropathies on the 3D kinematics of the lower limbs during gait: analysis with an ambulatory system.

Introduction: Ankle arthropathy is associated with a decreased motion of the ankle-hindfoot during ambulation. Ankle arthrodesis was shown to result in degeneration of the neighbour joints of the foot. Inversely, total ankle arthroplasty conceptually preserves the adjacent joints because of the residual mobility of the ankle but this has not been demonstrated yet in vivo. It has also been reported that degenerative ankle diseases, and even arthrodesis, do not result in alteration of the knee and hip joints. We present the preliminary results of a new approach of this problem based on ambulatory gait analysis. Patients and Methods: Motion analysis of the lower limbs was performed using a Physilog® (BioAGM, CH) system consisting of three-dimensional (3D) accelerometer and gyroscope, coupled to a magnetic system (Liberty©, Polhemus, USA). Both systems have been validated. Three groups of two patients were included into this pilot study and compared to healthy subjects (controls) during level walking: patients with ankle osteoarthritis (group 1), patients treated by ankle arthrodesis (group 2), patients treated by total ankle prosthesis (group 3). Results: Motion patterns of all analyzed joints over more than 20 gait cycles in each subject were highly repeatable. Motion amplitude of the ankle-hindfoot in control patients was similar to recently reported results. Ankle arthrodesis limited the motion of the ankle-hindfoot in the sagittal and horizontal planes. The prosthetic ankle allowed a more physiologic movement in the sagittal plane only. Ankle arthritis and its treatments did not influence the range of motion of the knee and hip joint during stance phase, excepted for a slight decrease of the hip flexion in groups 1 and 2. Conclusion: The reliability of the system was shown by the repeatability of the consecutive measurements. The results of this preliminary study were similar to those obtained through laboratory gait analysis. However, our system has the advantage to allow ambulatory analysis of 3D kinematics of the lower limbs outside of a gait laboratory and in real life conditions. To our knowledge this is a new concept in the analysis of ankle arthropathy and its treatments. Therefore, there is a potential to address specific questions like the difficult comparison of the benefits of ankle arthroplasty versus arthrodesis. The encouraging results of this pilot study offer the perspective to analyze the consequences of ankle arthropathy and its treatments on the biomechanics of the lower limbs ambulatory, in vivo and in daily life conditions.

A modular magnetic anastomotic device for minimally invasive digestive anastomosis: proof of concept and preliminary data in the pig model.

BACKGROUND: The aim of our study was to assess the feasibility of minimally invasive digestive anastomosis using a modular flexible magnetic anastomotic device made up of a set of two flexible chains of magnetic elements. The assembly possesses a non-deployed linear configuration which allows it to be introduced through a dedicated small-sized applicator into the bowel where it takes the deployed form. A centering suture allows the mating between the two parts to be controlled in order to include the viscerotomy between the two magnetic rings and the connected viscera. METHODS AND PROCEDURES: Eight pigs were involved in a 2-week survival experimental study. In five colorectal anastomoses, the proximal device was inserted by a percutaneous endoscopic technique, and the colon was divided below the magnet. The distal magnet was delivered transanally to connect with the proximal magnet. In three jejunojejunostomies, the first magnetic chain was injected in its linear configuration through a small enterotomy. Once delivered, the device self-assembled into a ring shape. A second magnet was injected more distally through the same port. The centering sutures were tied together extracorporeally and, using a knot pusher, magnets were connected. Ex vivo strain testing to determine the compression force delivered by the magnetic device, burst pressure of the anastomosis, and histology were performed. RESULTS: Mean operative time including endoscopy was 69.2 ± 21.9 min, and average time to full patency was 5 days for colorectal anastomosis. Operative times for jejunojejunostomies were 125, 80, and 35 min, respectively. The postoperative period was uneventful. Burst pressure of all anastomoses was ≥ 110 mmHg. Mean strain force to detach the devices was 6.1 ± 0.98 and 12.88 ± 1.34 N in colorectal and jejunojejunal connections, respectively. Pathology showed a mild-to-moderate inflammation score. CONCLUSIONS: The modular magnetic system showed enormous potential to create minimally invasive digestive anastomoses, and may represent an alternative to stapled anastomoses, being easy to deliver, effective, and low cost.

Synchronizability of EEG-Based Functional Networks in Early Alzheimer's Disease.

Recently graph theory and complex networks have been widely used as a mean to model functionality of the brain. Among different neuroimaging techniques available for constructing the brain functional networks, electroencephalography (EEG) with its high temporal resolution is a useful instrument of the analysis of functional interdependencies between different brain regions. Alzheimer's disease (AD) is a neurodegenerative disease, which leads to substantial cognitive decline, and eventually, dementia in aged people. To achieve a deeper insight into the behavior of functional cerebral networks in AD, here we study their synchronizability in 17 newly diagnosed AD patients compared to 17 healthy control subjects at no-task, eyes-closed condition. The cross-correlation of artifact-free EEGs was used to construct brain functional networks. The extracted networks were then tested for their synchronization properties by calculating the eigenratio of the Laplacian matrix of the connection graph, i.e., the largest eigenvalue divided by the second smallest one. In AD patients, we found an increase in the eigenratio, i.e., a decrease in the synchronizability of brain networks across delta, alpha, beta, and gamma EEG frequencies within the wide range of network costs. The finding indicates the destruction of functional brain networks in early AD.

Protein-protein interactions between adenovirus DNA polymerase and nuclear factor I mediate formation of the DNA replication preinitiation complex.

The in vitro adenovirus (Ad) DNA replication system provides an assay to study the interaction of viral and host replication proteins with the DNA template in the formation of the preinitiation complex. This initiation system requires in addition to the origin DNA sequences 1) Ad DNA polymerase (Pol), 2) Ad preterminal protein (pTP), the covalent acceptor for protein-primed DNA replication, and 3) nuclear factor I (NFI), a host cell protein identical to the CCAAT box-binding transcription factor. The interactions of these proteins were studied by coimmunoprecipitation and Ad origin DNA binding assays. The Ad Pol can bind to origin sequences only in the presence of another protein which can be either pTP or NFI. While NFI alone can bind to its origin recognition sequence, pTP does not specifically recognize DNA unless Ad Pol is present. Thus, protein-protein interactions are necessary for the targetting of either Ad Pol or pTP to the preinitiation complex. DNA footprinting demonstrated that the Ad DNA site recognized by the pTP.Pol complex was within the first 18 bases at the end of the template which constitutes the minimal origin of replication. Mutagenesis studies have defined the Ad Pol interaction site on NFI between amino acids 68-150, which overlaps the DNA binding and replication activation domain of this factor. A putative zinc finger on the Ad Pol has been mutated to a product that fails to bind the Ad origin sequences but still interacts with pTP. These results indicate that both protein-protein and protein-DNA interactions mediate specific recognition of the replication origin by Ad DNA polymerase.

Local optima networks of hard combinatorial landscapes

Combinatorial optimization involves finding an optimal solution in a finite set of options; many everyday life problems are of this kind. However, the number of options grows exponentially with the size of the problem, such that an exhaustive search for the best solution is practically infeasible beyond a certain problem size. When efficient algorithms are not available, a practical approach to obtain an approximate solution to the problem at hand, is to start with an educated guess and gradually refine it until we have a good-enough solution. Roughly speaking, this is how local search heuristics work. These stochastic algorithms navigate the problem search space by iteratively turning the current solution into new candidate solutions, guiding the search towards better solutions. The search performance, therefore, depends on structural aspects of the search space, which in turn depend on the move operator being used to modify solutions. A common way to characterize the search space of a problem is through the study of its fitness landscape, a mathematical object comprising the space of all possible solutions, their value with respect to the optimization objective, and a relationship of neighborhood defined by the move operator. The landscape metaphor is used to explain the search dynamics as a sort of potential function. The concept is indeed similar to that of potential energy surfaces in physical chemistry. Borrowing ideas from that field, we propose to extend to combinatorial landscapes the notion of the inherent network formed by energy minima in energy landscapes. In our case, energy minima are the local optima of the combinatorial problem, and we explore several definitions for the network edges. At first, we perform an exhaustive sampling of local optima basins of attraction, and define weighted transitions between basins by accounting for all the possible ways of crossing the basins frontier via one random move. Then, we reduce the computational burden by only counting the chances of escaping a given basin via random kick moves that start at the local optimum. Finally, we approximate network edges from the search trajectory of simple search heuristics, mining the frequency and inter-arrival time with which the heuristic visits local optima. Through these methodologies, we build a weighted directed graph that provides a synthetic view of the whole landscape, and that we can characterize using the tools of complex networks science. We argue that the network characterization can advance our understanding of the structural and dynamical properties of hard combinatorial landscapes. We apply our approach to prototypical problems such as the Quadratic Assignment Problem, the NK model of rugged landscapes, and the Permutation Flow-shop Scheduling Problem. We show that some network metrics can differentiate problem classes, correlate with problem non-linearity, and predict problem hardness as measured from the performances of trajectory-based local search heuristics.

E2F activation of S phase promoters via association with HCF-1 and the MLL family of histone H3K4 methyltransferases.

E2F transcriptional regulators control human-cell proliferation by repressing and activating the transcription of genes required for cell-cycle progression, particularly the S phase. E2F proteins repress transcription in association with retinoblastoma pocket proteins, but less is known about how they activate transcription. Here, we show that the human G1 phase regulator HCF-1 associates with both activator (E2F1 and E2F3a) and repressor (E2F4) E2F proteins, properties that are conserved in insect cells. Human HCF-1-E2F interactions are versatile: their associations and binding to E2F-responsive promoters are cell-cycle selective, and HCF-1 displays coactivator properties when bound to the E2F1 activator and corepressor properties when bound to the E2F4 repressor. During the G1-to-S phase transition, HCF-1 recruits the mixed-lineage leukemia (MLL) and Set-1 histone H3 lysine 4 methyltransferases to E2F-responsive promoters and induces histone methylation and transcriptional activation. These results suggest that HCF-1 induces cell-cycle-specific transcriptional activation by E2F proteins to promote cell proliferation.

Plant defense against herbivory: progress in identifying synergism, redundancy, and antagonism between resistance traits.

Plants respond to herbivore attack through a complex and variable system of defense, involving different physical barriers, toxic chemicals, and recruitment of natural enemies. To fully understand the relative role of each type of defense, their synergisms, redundancies, or antagonisms between traits, a variety of methods of enquiry, commonly used in plant physiology and ecology, have been employed. By overexpressing or silencing genes of interest, it is possible to understand the specific role of a particular defensive molecule or mode of action. We argue, however, that these types of experiments alone are not enough to holistically understand the physiological as well as ecological role of plant defenses. We thus advocate for the use of a combination of methods, including genetic modification, quantitative genetics, and phylogenetically controlled comparative studies.

Arbuscular mycorrhizal fungi mediate below-ground plant-herbivore interactions: a phylogenetic study

Ecological interactions are complex networks, but have typically been studied in a pairwise fashion. Examining how third-party species can modify the outcome of pairwise interactions may allow us to better predict their outcomes in realistic systems. For instance, arbuscular mycorrhizal fungi (AMF) can affect plant interactions with other organisms, including below-ground herbivores, but the mechanisms underlying these effects remain unclear. Here, we use a comparative, phylogenetically controlled approach to test the relative importance of mycorrhizal colonization and plant chemical defences (cardenolides) in predicting plant survival and the abundance of a generalist below-ground herbivore across 14 species of milkweeds (Asclepias spp.). Plants were inoculated with a mixture of four generalist AMF species or left uninoculated. After 1month, larvae of Bradysia sp. (Diptera: Sciaridae), a generalist below-ground herbivore, colonized plant roots. We performed phylogenetically controlled analyses to assess the influence of AMF colonization and toxic cardenolides on plant growth, mortality and infestation by fungus gnats. Overall, plants inoculated with AMF exhibited greater survival than did uninoculated plants. Additionally, surviving inoculated plants had lower numbers of larvae in their roots and fewer non-AM fungi than surviving uninoculated plants. In phylogenetic controlled regressions, gnat density in roots was better predicted by the extent of root colonized by AMF than by root cardenolide concentration. Taken as a whole, AMF modify the effect of below-ground herbivores on plants in a species-specific manner, independent of changes in chemical defence. This study adds to the growing body of literature demonstrating that mycorrhizal fungi may improve plant fitness by conferring protection against antagonists, rather than growth benefits. In addition, we advocate using comparative analyses to disentangle the roles of shared history and ecology in shaping trait expression and to better predict the outcomes of complex multitrophic interactions.

A high-resolution anatomical atlas of the transcriptome in the mouse embryo.

Ascertaining when and where genes are expressed is of crucial importance to understanding or predicting the physiological role of genes and proteins and how they interact to form the complex networks that underlie organ development and function. It is, therefore, crucial to determine on a genome-wide level, the spatio-temporal gene expression profiles at cellular resolution. This information is provided by colorimetric RNA in situ hybridization that can elucidate expression of genes in their native context and does so at cellular resolution. We generated what is to our knowledge the first genome-wide transcriptome atlas by RNA in situ hybridization of an entire mammalian organism, the developing mouse at embryonic day 14.5. This digital transcriptome atlas, the Eurexpress atlas (http://www.eurexpress.org), consists of a searchable database of annotated images that can be interactively viewed. We generated anatomy-based expression profiles for over 18,000 coding genes and over 400 microRNAs. We identified 1,002 tissue-specific genes that are a source of novel tissue-specific markers for 37 different anatomical structures. The quality and the resolution of the data revealed novel molecular domains for several developing structures, such as the telencephalon, a novel organization for the hypothalamus, and insight on the Wnt network involved in renal epithelial differentiation during kidney development. The digital transcriptome atlas is a powerful resource to determine co-expression of genes, to identify cell populations and lineages, and to identify functional associations between genes relevant to development and disease.

Bucket-handle tear of the triangular fibrocartilage complex : case report of a complex peripheral injury with separation of the distal radioulnar ligaments from the articular disc.

Palmer previously proposed a classification system of triangular fibrocartilage complex (TFCC) injuries that proved to be useful in directing clinical management. However, dorsal peripheral tears (variants of class 1C) were not described and have rarely been reported in the literature since. We herewith present a rare case of bucket-handle tear of the TFCC. To our knowledge, this is the first case demonstrating partial separation of both the palmar and dorsal distal radioulnar ligaments (DRULs) from the articular disc. The particular wrist magnetic resonance (MR) arthrographic findings of this unusual complex peripheral TFCC tear (a variant of both class 1B and 1C) were nicely appreciated upon sagittal reformatted images.

Arrested kinetic Li isotope fractionation at the margin of the llimaussaq complex, South Greenland: Evidence for open-system processes during final cooling of peralkaline igneous rocks

Li contents [Li] and isotopic composition (delta Li-7) of mafic minerals (mainly amphibole and clinopyroxene) from the alkaline to peralkaline Ilimaussaq plutonic complex, South Greenland, track the behavior of Li and its isotopes during magmatic differentiation and final cooling of an alkaline igneous system. [Li] in amphibole increase from < 10 ppm in Caamphiboles of the least differentiated unit to >3000 ppm in Na-amphiboles of the highly evolved units. In contrast, [Li] in clinopyroxene are comparatively low (<85 ppm) and do not vary systematically with differentiation. The distribution of Li between amphibole and pyroxene is controlled by the major element composition of the minerals (Ca-rich and Na-rich, respectively) and changes in oxygen fugacity (due to Li incorporation via coupled substitution with ferric iron) during magmatic differentiation. delta(7) Li values of all minerals span a wide range from + 17 to - 8 parts per thousand, with the different intrusive units of the complex having distinct Li isotopic systematics. Amphiboles, which dominate the Li budget of whole-rocks from the inner part of the complex, have constant delta Li-7 of + 1.8 +/- 2.2 parts per thousand (2 sigma, n = 15). This value reflects a homogeneous melt reservoir and is consistent with their mantle derivation, in agreement with published O and Nd isotopic data. Clinopyroxenes of these samples are consistently lighter, with Delta Li-7(amph-cpx). as large as 8 parts per thousand and are thus not in Li isotope equilibrium. These low values probably reflect late-stage diffusion of Li into clinopyroxene during final cooling of the rocks, thus enriching the clinopyroxene in 6 Li. At the margin of the complex delta(7) Li in the syenites increases systematically, from +2 to high values of + 14 parts per thousand. This, coupled with the observed Li isotope systematics of the granitic country rocks, reflects post-magmatic open-system processes occurring during final cooling of the intrusion. Although the shape and magnitude of the Li isotope and elemental profiles through syenite and country rock are suggestive of diffusion-driven isotope fractionation, they cannot be modeled by one-dimensional diffusive transport and point to circulation of a fluid having a high 67 Li value (possibly seawater) along the chilled contact. In all, this study demonstrates that Li isotopes can be used to identify complex fluid- and diffusion-governed processes taking place during the final cooling of such rocks. (c) 2007 Elsevier B.V All rights reserved.