Zircon ages in granulite facies rocks: decoupling from geochemistry above 850 °C?

Granulite facies rocks frequently show a large spread in their zircon ages, the interpretation of which raises questions: Has the isotopic system been disturbed? By what process(es) and conditions did the alteration occur? Can the dates be regarded as real ages, reflecting several growth episodes? Furthermore, under some circumstances of (ultra-)high-temperature metamorphism, decoupling of zircon U–Pb dates from their trace element geochemistry has been reported. Understanding these processes is crucial to help interpret such dates in the context of the P–T history. Our study presents evidence for decoupling in zircon from the highest grade metapelites (> 850 °C) taken along a continuous high-temperature metamorphic field gradient in the Ivrea Zone (NW Italy). These rocks represent a well-characterised segment of Permian lower continental crust with a protracted high-temperature history. Cathodoluminescence images reveal that zircons in the mid-amphibolite facies preserve mainly detrital cores with narrow overgrowths. In the upper amphibolite and granulite facies, preserved detrital cores decrease and metamorphic zircon increases in quantity. Across all samples we document a sequence of four rim generations based on textures. U–Pb dates, Th/U ratios and Ti-in-zircon concentrations show an essentially continuous evolution with increasing metamorphic grade, except in the samples from the granulite facies, which display significant scatter in age and chemistry. We associate the observed decoupling of zircon systematics in high-grade non-metamict zircon with disturbance processes related to differences in behaviour of non-formula elements (i.e. Pb, Th, U, Ti) at high-temperature conditions, notably differences in compatibility within the crystal structure

Topological Structure of the Space of Phenotypes: The Case of RNA Neutral Networks

The evolution and adaptation of molecular populations is constrained by the diversity accessible through mutational processes. RNA is a paradigmatic example of biopolymer where genotype (sequence) and phenotype (approximated by the secondary structure fold) are identified in a single molecule. The extreme redundancy of the genotype-phenotype map leads to large ensembles of RNA sequences that fold into the same secondary structure and can be connected through single-point mutations. These ensembles define neutral networks of phenotypes in sequence space. Here we analyze the topological properties of neutral networks formed by 12-nucleotides RNA sequences, obtained through the exhaustive folding of sequence space. A total of 412 sequences fragments into 645 subnetworks that correspond to 57 different secondary structures. The topological analysis reveals that each subnetwork is far from being random: it has a degree distribution with a well-defined average and a small dispersion, a high clustering coefficient, and an average shortest path between nodes close to its minimum possible value, i.e. the Hamming distance between sequences. RNA neutral networks are assortative due to the correlation in the composition of neighboring sequences, a feature that together with the symmetries inherent to the folding process explains the existence of communities. Several topological relationships can be analytically derived attending to structural restrictions and generic properties of the folding process. The average degree of these phenotypic networks grows logarithmically with their size, such that abundant phenotypes have the additional advantage of being more robust to mutations. This property prevents fragmentation of neutral networks and thus enhances the navigability of sequence space. In summary, RNA neutral networks show unique topological properties, unknown to other networks previously described

Efficient Algorithms for Probing the RNA Mutation Landscape

The diversity and importance of the role played by RNAs in the regulation and development of the cell are now well-known and well-documented. This broad range of functions is achieved through specific structures that have been (presumably) optimized through evolution. State-of-the-art methods, such as McCaskill's algorithm, use a statistical mechanics framework based on the computation of the partition function over the canonical ensemble of all possible secondary structures on a given sequence. Although secondary structure predictions from thermodynamics-based algorithms are not as accurate as methods employing comparative genomics, the former methods are the only available tools to investigate novel RNAs, such as the many RNAs of unknown function recently reported by the ENCODE consortium. In this paper, we generalize the McCaskill partition function algorithm to sum over the grand canonical ensemble of all secondary structures of all mutants of the given sequence. Specifically, our new program, RNAmutants, simultaneously computes for each integer k the minimum free energy structure MFE(k) and the partition function Z(k) over all secondary structures of all k-point mutants, even allowing the user to specify certain positions required not to mutate and certain positions required to base-pair or remain unpaired. This technically important extension allows us to study the resilience of an RNA molecule to pointwise mutations. By computing the mutation profile of a sequence, a novel graphical representation of the mutational tendency of nucleotide positions, we analyze the deleterious nature of mutating specific nucleotide positions or groups of positions. We have successfully applied RNAmutants to investigate deleterious mutations (mutations that radically modify the secondary structure) in the Hepatitis C virus cis-acting replication element and to evaluate the evolutionary pressure applied on different regions of the HIV trans-activation response element. In particular, we show qualitative agreement between published Hepatitis C and HIV experimental mutagenesis studies and our analysis of deleterious mutations using RNAmutants. Our work also predicts other deleterious mutations, which could be verified experimentally. Finally, we provide evidence that the 3′ UTR of the GB RNA virus C has been optimized to preserve evolutionarily conserved stem regions from a deleterious effect of pointwise mutations. We hope that there will be long-term potential applications of RNAmutants in de novo RNA design and drug design against RNA viruses. This work also suggests potential applications for large-scale exploration of the RNA sequence-structure network. Binary distributions are available at http://RNAmutants.csail.mit.edu/

Integrative inference of gene-regulatory networks in Escherichia coli using information theoretic concepts and sequence analysis

<p>Abstract</p> <p>Background</p> <p>Although <it>Escherichia coli </it>is one of the best studied model organisms, a comprehensive understanding of its gene regulation is not yet achieved. There exist many approaches to reconstruct regulatory interaction networks from gene expression experiments. Mutual information based approaches are most useful for large-scale network inference.</p> <p>Results</p> <p>We used a three-step approach in which we combined gene regulatory network inference based on directed information (DTI) and sequence analysis. DTI values were calculated on a set of gene expression profiles from 19 time course experiments extracted from the Many Microbes Microarray Database. Focusing on influences between pairs of genes in which one partner encodes a transcription factor (TF) we derived a network which contains 878 TF - gene interactions of which 166 are known according to RegulonDB. Afterward, we selected a subset of 109 interactions that could be confirmed by the presence of a phylogenetically conserved binding site of the respective regulator. By this second step, the fraction of known interactions increased from 19% to 60%. In the last step, we checked the 44 of the 109 interactions not yet included in RegulonDB for functional relationships between the regulator and the target and, thus, obtained ten TF - target gene interactions. Five of them concern the regulator LexA and have already been reported in the literature. The remaining five influences describe regulations by Fis (with two novel targets), PhdR, PhoP, and KdgR. For the validation of our approach, one of them, the regulation of lipoate synthase (LipA) by the pyruvate-sensing pyruvate dehydrogenate repressor (PdhR), was experimentally checked and confirmed.</p> <p>Conclusions</p> <p>We predicted a set of five novel TF - target gene interactions in <it>E. coli</it>. One of them, the regulation of <it>lipA </it>by the transcriptional regulator PdhR was validated experimentally. Furthermore, we developed DTInfer, a new R-package for the inference of gene-regulatory networks from microarrays using directed information.</p

Quantitative global studies of reactomes and metabolomes using a vectorial representation of reactions and chemical compounds

<p>Abstract</p> <p>Background</p> <p>Global studies of the protein repertories of organisms are providing important information on the characteristics of the protein space. Many of these studies entail classification of the protein repertory on the basis of structure and/or sequence similarities. The situation is different for metabolism. Because there is no good way of measuring similarities between chemical reactions, there is a barrier to the development of global classifications of "metabolic space" and subsequent studies comparable to those done for protein sequences and structures.</p> <p>Results</p> <p>In this work, we propose a vectorial representation of chemical reactions, which allows them to be compared and classified. In this representation, chemical compounds, reactions and pathways may be represented in the same vectorial space. We show that the representation of chemical compounds reflects their physicochemical properties and can be used for predictive purposes. We use the vectorial representations of reactions to perform a global classification of the reactome of the model organism <it>E. coli</it>.</p> <p>Conclusions</p> <p>We show that this unsupervised clustering results in groups of enzymes more coherent in biological terms than equivalent groupings obtained from the EC hierarchy. This hierarchical clustering produces an optimal set of 21 groups which we analyzed for their biological meaning.</p

Implementing Routine HIV Testing: The Role of State Law

In September 2006, the Centers for Disease Control and Prevention (CDC) recommended routine HIV testing for all Americans aged 13–64, which would eliminate requirements for written consent and pretest counseling as previously required. However, this approach may conflict with state requirements concerning pretest counseling and informed consent for HIV testing. Our survey of state HIV testing laws demonstrates that the majority of states have HIV testing requirements that are inconsistent with the CDC's recommendations. Moreover, states that have recently amended their laws have not eased the requirements for pretest counseling and informed consent. The reasons for the persistence of these legal requirements must be understood to effect policy changes to increase HIV testing

Cross-Talk between Signaling Pathways Can Generate Robust Oscillations in Calcium and cAMP

BACKGROUND:To control and manipulate cellular signaling, we need to understand cellular strategies for information transfer, integration, and decision-making. A key feature of signal transduction is the generation of only a few intracellular messengers by many extracellular stimuli. METHODOLOGY/PRINCIPAL FINDINGS:Here we model molecular cross-talk between two classic second messengers, cyclic AMP (cAMP) and calcium, and show that the dynamical complexity of the response of both messengers increases substantially through their interaction. In our model of a non-excitable cell, both cAMP and calcium concentrations can oscillate. If mutually inhibitory, cross-talk between the two second messengers can increase the range of agonist concentrations for which oscillations occur. If mutually activating, cross-talk decreases the oscillation range, but can generate 'bursting' oscillations of calcium and may enable better filtering of noise. CONCLUSION:We postulate that this increased dynamical complexity allows the cell to encode more information, particularly if both second messengers encode signals. In their native environments, it is unlikely that cells are exposed to one stimulus at a time, and cross-talk may help generate sufficiently complex responses to allow the cell to discriminate between different combinations and concentrations of extracellular agonists

Enhanced pharmacological efficacy of sumatriptan due to modification of its physicochemical properties by inclusion in selected cyclodextrins

The study focused on the pharmacological action of sumatriptan, in particular its antiallodynic and antihyperalgesic properties, as an effect of cyclodextrinic inclusion of sumatriptan, resulting in changes of its physicochemical qualities such as dissolution and permeability through artificial biological membranes, which had previously been examined in vitro in a gastro-intestinal model. The inclusion of sumatriptan into β-cyclodextrin and 2-hydroxylpropylo-β-cyclodextrin by kneading was confirmed with the use of spectral (fourier-transform infrared spectroscopy (FT-IR); solid state nuclear magnetic resonance spectroscopy with magic angle spinning condition, 1H and 13C MAS NMR) and thermal (differential scanning calorimetry (DSC)) methods. A precise indication of the domains of sumatriptan responsible for its interaction with cyclodextrin cavities was possible due to a theoretical approach to the analysis of experimental spectra. A high-performance liquid chromatography with a diode-array detector method (HPLC-DAD) was employed to determine changes in the concentration of sumatriptan during dissolution and permeability experiments. The inclusion of sumatriptan in complex with cyclodextrins was found to significantly modify its dissolution profiles by increasing the concentration of sumatriptan in complexed form in an acceptor solution compared to in its free form. Following complexation, sumatriptan manifested an enhanced ability to permeate through artificial biological membranes in a gastro-intestinal model for both cyclodextrins at all pH values. As a consequence of the greater permeability of sumatriptan and its increased dissolution from the complexes, an improved pharmacological response was observed when cyclodextrin complexes were applied

A Genetic Screen for Dihydropyridine (DHP)-Resistant Worms Reveals New Residues Required for DHP-Blockage of Mammalian Calcium Channels

Dihydropyridines (DHPs) are L-type calcium channel (Cav1) blockers prescribed to treat several diseases including hypertension. Cav1 channels normally exist in three states: a resting closed state, an open state that is triggered by membrane depolarization, followed by a non-conducting inactivated state that is triggered by the influx of calcium ions, and a rapid change in voltage. DHP binding is thought to alter the conformation of the channel, possibly by engaging a mechanism similar to voltage dependent inactivation, and locking a calcium ion in the pore, thereby blocking channel conductance. As a Cav1 channel crystal structure is lacking, the current model of DHP action has largely been achieved by investigating the role of candidate Cav1 residues in mediating DHP-sensitivity. To better understand DHP-block and identify additional Cav1 residues important for DHP-sensitivity, we screened 440,000 randomly mutated Caenorhabditis elegans genomes for worms resistant to DHP-induced growth defects. We identified 30 missense mutations in the worm Cav1 pore-forming (α1) subunit, including eleven in conserved residues known to be necessary for DHP-binding. The remaining polymorphisms are in eight conserved residues not previously associated with DHP-sensitivity. Intriguingly, all of the worm mutants that we analyzed phenotypically exhibited increased channel activity. We also created orthologous mutations in the rat α1C subunit and examined the DHP-block of current through the mutant channels in culture. Six of the seven mutant channels examined either decreased the DHP-sensitivity of the channel and/or exhibited significant residual current at DHP concentrations sufficient to block wild-type channels. Our results further support the idea that DHP-block is intimately associated with voltage dependent inactivation and underscores the utility of C. elegans as a screening tool to identify residues important for DHP interaction with mammalian Cav1 channels

Next generation transcriptomes for next generation genomes using est2assembly

<p>Abstract</p> <p>Background</p> <p>The decreasing costs of capillary-based Sanger sequencing and next generation technologies, such as 454 pyrosequencing, have prompted an explosion of transcriptome projects in non-model species, where even shallow sequencing of transcriptomes can now be used to examine a range of research questions. This rapid growth in data has outstripped the ability of researchers working on non-model species to analyze and mine transcriptome data efficiently.</p> <p>Results</p> <p>Here we present a semi-automated platform '<it>est2assembly</it>' that processes raw sequence data from Sanger or 454 sequencing into a hybrid <it>de-novo </it>assembly, annotates it and produces GMOD compatible output, including a SeqFeature database suitable for GBrowse. Users are able to parameterize assembler variables, judge assembly quality and determine the optimal assembly for their specific needs. We used <it>est2assembly </it>to process <it>Drosophila </it>and <it>Bicyclus </it>public Sanger EST data and then compared them to published 454 data as well as eight new insect transcriptome collections.</p> <p>Conclusions</p> <p>Analysis of such a wide variety of data allows us to understand how these new technologies can assist EST project design. We determine that assembler parameterization is as essential as standardized methods to judge the output of ESTs projects. Further, even shallow sequencing using 454 produces sufficient data to be of wide use to the community. <it>est2assembly </it>is an important tool to assist manual curation for gene models, an important resource in their own right but especially for species which are due to acquire a genome project using Next Generation Sequencing.</p