Alkali-calcic and alkaline post-orogenic (PO) granite magmatism: petrologic constraints and geodynamic settings

The end of an orogenic Wilson cycle corresponds to amalgamation of terranes into a Pangaea and is marked by widespread magmatism dominated by granitoids. The post-collision event starts with magmatic processes still influenced by subducted crustal materials. The dominantly calc-alkaline suites show a shift from normal to high-K to very high-K associations. Source regions are composed of depleted and later enriched orogenic subcontinental lithospheric mantle, affected by dehydration melting and generating more and more K- and LILE-rich magmas. In the vicinity of intra-crustal magma chambers, anatexis by incongruent melting of hydrous minerals may generate peraluminous granitoids bearing mafic enclaves. The post-collision event ends with emplacement of bimodal post-orogenic (PO) suites along transcurrent fault zones. Two suites are defined, (i) the alkali-calcic monzonite-monzogranite-syenogranite-alkali feldspar granite association characterised by [biotite + plagioclase] fractionation and moderate [LILE + HFSE] enrichments and (ii) the alkaline monzonite-syenite-alkali feldspar granite association characterised by [amphibole + alkali feldspar] fractionation and displaying two evolutionary trends, one peralkaline with sodic mafic mineralogy and higher enrichments in HFSE than in LILE, and the other aluminous biotite-bearing marked by HFSE depletion relative to LILE due to accessory mineral precipitation. Alkali-calcic and alkaline suites differ essentially in the amounts of water present within intra-crustal magma chambers, promoting crystallisation of various mineral assemblages. The ultimate enriched and not depleted mantle source is identical for the two PO suites. The more primitive LILE and HFSE-rich source rapidly replaces the older orogenic mantle source during lithosphere delamination and becomes progressively the thermal boundary layer of the new lithosphere. Present rock compositions are a mixture of major mantle contribution and various crustal components carried by F-rich aqueous fluids circulating within convective cells created around magma chambers. In favourable areas, PO suites pre-date a new orogenic Wilson cycle. (C) 1998 Elsevier Science B.V. All rights reserved.

MAP/ERK kinase kinase 1 (MEKK1) mediates transcriptional repression by interacting with polycystic kidney disease-1 (PKD1) promoter-bound p53 tumor suppressor protein.

Mitogen-activated protein kinase (MAPK) cascades regulate a wide variety of cellular processes that ultimately depend on changes in gene expression. We have found a novel mechanism whereby one of the key MAP3 kinases, Mekk1, regulates transcriptional activity through an interaction with p53. The tumor suppressor protein p53 down-regulates a number of genes, including the gene most frequently mutated in autosomal dominant polycystic kidney disease (PKD1). We have discovered that Mekk1 translocates to the nucleus and acts as a co-repressor with p53 to down-regulate PKD1 transcriptional activity. This repression does not require Mekk1 kinase activity, excluding the need for an Mekk1 phosphorylation cascade. However, this PKD1 repression can also be induced by the stress-pathway stimuli, including TNFα, suggesting that Mekk1 activation induces both JNK-dependent and JNK-independent pathways that target the PKD1 gene. An Mekk1-p53 interaction at the PKD1 promoter suggests a new mechanism by which abnormally elevated stress-pathway stimuli might directly down-regulate the PKD1 gene, possibly causing haploinsufficiency and cyst formation.

Comparative genomics of the vertebrate insulin/TOR signal transduction pathway: A network-level analysis of selective pressures

Complexity of biological function relies on large networks of interacting molecules. However, the evolutionary properties of these networks are not fully understood. It has been shown that selective pressures depend on the position of genes in the network. We have previously shown that in the Drosophila insulin/target of rapamycin (TOR) signal transduction pathway there is a correlation between the pathway position and the strength of purifying selection, with the downstream genes being most constrained. In this study, we investigated the evolutionary dynamics of this well-characterized pathway in vertebrates. More specifically, we determined the impact of natural selection on the evolution of 72 genes of this pathway. We found that in vertebrates there is a similar gradient of selective constraint in the insulin/TOR pathway to that found in Drosophila. This feature is neither the result of a polarity in the impact of positive selection nor of a series of factors affecting selective constraint levels (gene expression level and breadth, codon bias, protein length, and connectivity). We also found that pathway genes encoding physically interacting proteins tend to evolve under similar selective constraints. The results indicate that the architecture of the vertebrate insulin/TOR pathway constrains the molecular evolution of its components. Therefore, the polarity detected in Drosophila is neither specific nor incidental of this genus. Hence, although the underlying biological mechanisms remain unclear, these may be similar in both vertebrates and Drosophila.

Role of the JNK-interacting protein 1/islet brain 1 in cell degeneration in Alzheimer disease and diabetes.

Numerous epidemiological studies and some pharmacological clinical trials show the close connection between Alzheimer disease (AD) and type 2 diabetes (T2D) and thereby, shed more light into the existence of possible similar pathogenic mechanisms between these two diseases. Diabetes increases the risk of developing AD and sensitizers of insulin currently used as diabetes drugs can efficiently slow cognitive decline of the neurological disorder. Deposits of amyloid aggregate and hyperphosphorylation of tau, which are hallmarks of AD, have been also found in degenerating pancreatic islets beta-cells of patients with T2D. These events may have a causal role in the pathogenesis of the two diseases. Increased c-Jun NH(2)-terminal kinase (JNK) activity is found in neurofibrillary tangles (NFT) of AD and promotes programmed cell death of beta-cells exposed to a diabetic environment. The JNK-interacting protein 1 (JIP-1), also called islet brain 1 (IB1) because it is mostly expressed in the brain and islets, is a key regulator of the JNK pathway in neuronal and beta-cells. JNK, hyperphosphorylated tau and IB1/JIP-1 all co-localize with amyloids deposits in NFT and islets of AD and patients with T2D. This review discusses the role of the IB1/JIP-1 and the JNK pathway in the molecular pathogenesis of AD and T2D.

Integrated methodologies to study human transcriptional regulation from functional genomics data

Abstract : The human body is composed of a huge number of cells acting together in a concerted manner. The current understanding is that proteins perform most of the necessary activities in keeping a cell alive. The DNA, on the other hand, stores the information on how to produce the different proteins in the genome. Regulating gene transcription is the first important step that can thus affect the life of a cell, modify its functions and its responses to the environment. Regulation is a complex operation that involves specialized proteins, the transcription factors. Transcription factors (TFs) can bind to DNA and activate the processes leading to the expression of genes into new proteins. Errors in this process may lead to diseases. In particular, some transcription factors have been associated with a lethal pathological state, commonly known as cancer, associated with uncontrolled cellular proliferation, invasiveness of healthy tissues and abnormal responses to stimuli. Understanding cancer-related regulatory programs is a difficult task, often involving several TFs interacting together and influencing each other's activity. This Thesis presents new computational methodologies to study gene regulation. In addition we present applications of our methods to the understanding of cancer-related regulatory programs. The understanding of transcriptional regulation is a major challenge. We address this difficult question combining computational approaches with large collections of heterogeneous experimental data. In detail, we design signal processing tools to recover transcription factors binding sites on the DNA from genome-wide surveys like chromatin immunoprecipitation assays on tiling arrays (ChIP-chip). We then use the localization about the binding of TFs to explain expression levels of regulated genes. In this way we identify a regulatory synergy between two TFs, the oncogene C-MYC and SP1. C-MYC and SP1 bind preferentially at promoters and when SP1 binds next to C-NIYC on the DNA, the nearby gene is strongly expressed. The association between the two TFs at promoters is reflected by the binding sites conservation across mammals, by the permissive underlying chromatin states 'it represents an important control mechanism involved in cellular proliferation, thereby involved in cancer. Secondly, we identify the characteristics of TF estrogen receptor alpha (hERa) target genes and we study the influence of hERa in regulating transcription. hERa, upon hormone estrogen signaling, binds to DNA to regulate transcription of its targets in concert with its co-factors. To overcome the scarce experimental data about the binding sites of other TFs that may interact with hERa, we conduct in silico analysis of the sequences underlying the ChIP sites using the collection of position weight matrices (PWMs) of hERa partners, TFs FOXA1 and SP1. We combine ChIP-chip and ChIP-paired-end-diTags (ChIP-pet) data about hERa binding on DNA with the sequence information to explain gene expression levels in a large collection of cancer tissue samples and also on studies about the response of cells to estrogen. We confirm that hERa binding sites are distributed anywhere on the genome. However, we distinguish between binding sites near promoters and binding sites along the transcripts. The first group shows weak binding of hERa and high occurrence of SP1 motifs, in particular near estrogen responsive genes. The second group shows strong binding of hERa and significant correlation between the number of binding sites along a gene and the strength of gene induction in presence of estrogen. Some binding sites of the second group also show presence of FOXA1, but the role of this TF still needs to be investigated. Different mechanisms have been proposed to explain hERa-mediated induction of gene expression. Our work supports the model of hERa activating gene expression from distal binding sites by interacting with promoter bound TFs, like SP1. hERa has been associated with survival rates of breast cancer patients, though explanatory models are still incomplete: this result is important to better understand how hERa can control gene expression. Thirdly, we address the difficult question of regulatory network inference. We tackle this problem analyzing time-series of biological measurements such as quantification of mRNA levels or protein concentrations. Our approach uses the well-established penalized linear regression models where we impose sparseness on the connectivity of the regulatory network. We extend this method enforcing the coherence of the regulatory dependencies: a TF must coherently behave as an activator, or a repressor on all its targets. This requirement is implemented as constraints on the signs of the regressed coefficients in the penalized linear regression model. Our approach is better at reconstructing meaningful biological networks than previous methods based on penalized regression. The method is tested on the DREAM2 challenge of reconstructing a five-genes/TFs regulatory network obtaining the best performance in the "undirected signed excitatory" category. Thus, these bioinformatics methods, which are reliable, interpretable and fast enough to cover large biological dataset, have enabled us to better understand gene regulation in humans.

The hourglass and the early conservation models--co-existing patterns of developmental constraints in vertebrates.

Developmental constraints have been postulated to limit the space of feasible phenotypes and thus shape animal evolution. These constraints have been suggested to be the strongest during either early or mid-embryogenesis, which corresponds to the early conservation model or the hourglass model, respectively. Conflicting results have been reported, but in recent studies of animal transcriptomes the hourglass model has been favored. Studies usually report descriptive statistics calculated for all genes over all developmental time points. This introduces dependencies between the sets of compared genes and may lead to biased results. Here we overcome this problem using an alternative modular analysis. We used the Iterative Signature Algorithm to identify distinct modules of genes co-expressed specifically in consecutive stages of zebrafish development. We then performed a detailed comparison of several gene properties between modules, allowing for a less biased and more powerful analysis. Notably, our analysis corroborated the hourglass pattern at the regulatory level, with sequences of regulatory regions being most conserved for genes expressed in mid-development but not at the level of gene sequence, age, or expression, in contrast to some previous studies. The early conservation model was supported with gene duplication and birth that were the most rare for genes expressed in early development. Finally, for all gene properties, we observed the least conservation for genes expressed in late development or adult, consistent with both models. Overall, with the modular approach, we showed that different levels of molecular evolution follow different patterns of developmental constraints. Thus both models are valid, but with respect to different genomic features.

Islet-brain1/JNK-interacting protein-1 is required for early embryogenesis in mice.

Islet-brain1/JNK-interacting protein-1 (IB1/JIP-1) is a scaffold protein that organizes the JNK, MKK7, and MLK1 to allow signaling specificity. Targeted disruption of the gene MAPK8IP1 encoding IB1/JIP-1 in mice led to embryonic death prior to blastocyst implantation. In culture, no IB1/JIP-1(-/-) embryos were identified indicating that accelerated cell death occurred during the first cell cycles. IB1/JIP-1 expression was detected in unfertilized oocytes, in spermatozoa, and in different stages of embryo development. Thus, despite the maternal and paternal transmission of the IB1/JIP-1 protein, early transcription of the MAPK8IP1 gene is required for the survival of the fertilized oocytes.

Juxtaposition of melt impregnation and high-temperature shear zones in the upper mantle; field and petrological constraints from the Lanzo peridotite (Northern Italy)

Results of a field and microstructural study between the northern and the central bodies of the Lanzo plagioclase peridotite massif (NW Italy) indicate that the spatial distribution of deformation is asymmetric across kilometre-scale mantle shear zones. The southwestern part of the shear zone (footwall) shows a gradually increasing degree of deformation from porphyroclastic peridotites to mylonite, whereas the northeastern part (hanging wall) quickly grades into weakly deformed peridotites. Discordant gabbroic and basaltic dykes are asymmetrically distributed and far more abundant in the footwall of the shear zone. The porphyroclastic peridotite displays porphyroclastic zones and domains of igneous crystallization whereas mylonites are characterized by elongated porphyroclasts, embedded between fine-grained, polycrystalline bands of olivine, plagioclase, clinopyroxene, orthopyroxene, spinel, rare titanian pargasite, and domains of recrystallized olivine. Two types of melt impregnation textures have been found: (1) clinopyroxene porphyroclasts incongruently reacted with migrating melt to form orthopyroxene plagioclase; (2) olivine porphyroclasts are partially replaced by interstitial orthopyroxene. The meltrock reaction textures tend to disappear in the mylonites, indicating that deformation in the mylonite continued under subsolidus conditions. The pyroxene chemistry is correlated with grain size. High-Al pyroxene cores indicate high temperatures (11001030C), whereas low-Al neoblasts display lower final equilibration temperatures (860C). The spinel Cr-number [molar Cr/(Cr Al)] and TiO2 concentrations show extreme variability covering almost the entire range known from abyssal peridotites. The spinel compositions of porphyroclastic peridotites from the central body are more variable than spinel from mylonite, mylonite with ultra-mylonite bands, and porphyroclastic rocks of the northern body. The spinel compositions probably indicate disequilibrium and would favour rapid cooling, and a faster exhumation of the central peridotite body, relative to the northern one. Our results indicate that melt migration and high-temperature deformation are juxtaposed both in time and space. Meltrock reaction may have caused grain-size reduction, which in turn led to localization of deformation. It is likely that melt-lubricated, actively deforming peridotites acted as melt focusing zones, with permeabilities higher than the surrounding, less deformed peridotites. Later, under subsolidus conditions, pinning in polycrystalline bands in the mylonites inhibited substantial grain growth and led to permanent weak zones in the upper mantle peridotite, with a permeability that is lower than in the weakly deformed peridotites. Such an inversion in permeability might explain why actively deforming, fine-grained peridotite mylonite acted as a permeability barrier and why ascending mafic melts might terminate and crystallize as gabbros along actively deforming shear zones. Melt-lubricated mantle shear zones provide a mechanism for explaining the discontinuous distribution of gabbros in oceancontinent transition zones, oceanic core complexes and ultraslow-spreading ridges.

Constraints on the permeability structure of alluvial aquifers from the poro-elastic inversion of multi-frequency P-wave velocity logs

We have explored the possibility of obtaining first-order permeability estimates for saturated alluvial sediments based on the poro-elastic interpretation of the P-wave velocity dispersion inferred from sonic logs. Modern sonic logging tools designed for environmental and engineering applications allow one for P-wave velocity measurements at multiple emitter frequencies over a bandwidth covering 5 to 10 octaves. Methodological considerations indicate that, for saturated unconsolidated sediments in the silt to sand range and typical emitter frequencies ranging from approximately 1 to 30 kHz, the observable velocity dispersion should be sufficiently pronounced to allow one for reliable first-order estimations of the permeability structure. The corresponding predictions have been tested on and verified for a borehole penetrating a typical surficial alluvial aquifer. In addition to multifrequency sonic logs, a comprehensive suite of nuclear and electrical logs, an S-wave log, a litholog, and a limited number laboratory measurements of the permeability from retrieved core material were also available. This complementary information was found to be essential for parameterizing the poro-elastic inversion procedure and for assessing the uncertainty and internal consistency of corresponding permeability estimates. Our results indicate that the thus obtained permeability estimates are largely consistent with those expected based on the corresponding granulometric characteristics, as well as with the available evidence form laboratory measurements. These findings are also consistent with evidence from ocean acoustics, which indicate that, over a frequency range of several orders-of-magnitude, the classical theory of poro-elasticity is generally capable of explaining the observed P-wave velocity dispersion in medium- to fine-grained seabed sediments

The Upper Oligocene of Montgat (Catalan Coastal Ranges, Spain): new age constraints to the western Mediterranean Basin opening

The Oligocene deposits of Montgat are integrated in a small outcrop made up of Cenozoic and Mesozoic rocks located in the Garraf-Montnegre horst, close to the major Barcelona fault. The Oligocene of Montgat consists of detrital sediments of continental origin mainly deposited in alluvial fan environments; these deposits are folded and affected by thrusts and strike-slip faults. They can be divided in two lithostratigraphic units separated by a minor southwest-directed thrust: (i) the Turó de Montgat Unit composed of litharenites and lithorudites with high contents of quartz, feldspar, plutonic and limestone rock fragments; and (ii) the Pla de la Concòrdia Unit composed of calcilitharenites and calcilithorudites with high contents of dolosparite and dolomicrite rock fragments. The petrological composition of both units indicates that sediments were derived from the erosion of Triassic (Buntsandstein, Muschelkalk and Keuper facies), Jurassic and Lower Cretaceous rocks (Barremian to Aptian in age). Stratigraphic and petrological data suggest that these units correspond to two coalescent alluvial fans with a source area located northwestwards in the adjoining Collserola and Montnegre inner areas. Micromammal fossils (Archaeomys sp.) found in a mudstone layer of the Pla de la Concòrdia Unit assign a Chattian age (Late Oligocene) to the studied materials. Thus, the Montgat deposits are the youngest dated deposits affected by the contractional deformation that led to the development of the Catalan Intraplate Chain. Taking into account that the oldest syn-rift deposits in the Catalan Coastal Ranges are Aquitanian in age, this allows to precise that the change from a compressive to an extensional regime in this area took place during latest Oligocene-earliest Aquitanian times. This age indicates that the onset of crustal extension related to the opening of the western Mediterranean Basin started in southern France during latest Eocene-early Oligocene and propagated southwestward, affecting the Catalan Coastal Ranges and the northeastern part of the Valencia trough during the latest Chattian-earliest Aquitanian times.

Stochastic delay equations with non-negativity constraints driven by fractional Brownian motion

In this note we prove an existence and uniqueness result for the solution of multidimensional stochastic delay differential equations with normal reflection. The equations are driven by a fractional Brownian motion with Hurst parameter H > 1/2. The stochastic integral with respect to the fractional Brownian motion is a pathwise Riemann¿Stieltjes integral.