How to compare arc-annotated sequences: The alignment hierarchy

International audienceWe describe a new unifying framework to express comparison of arc-annotated sequences, which we call alignment of arc-annotated sequences. We first prove that this framework encompasses main existing models, which allows us to deduce complexity results for several cases from the literature. We also show that this framework gives rise to new relevant problems that have not been studied yet. We provide a thorough analysis of these novel cases by proposing two polynomial time algorithms and an NP-completeness proof. This leads to an almost exhaustive study of alignment of arc-annotated sequences

Changes of bivalent chromatin coincide with increased expression of developmental genes in cancer

Bivalent (poised or paused) chromatin comprises activating and repressing histone modifications at the same location. This combination of epigenetic marks at promoter or enhancer regions keeps genes expressed at low levels but poised for rapid activation. Typically, DNA at bivalent promoters is only lowly methylated in normal cells, but frequently shows elevated methylation levels in cancer samples. Here, we developed a universal classifier built from chromatin data that can identify cancer samples solely from hypermethylation of bivalent chromatin. Tested on over 7,000 DNA methylation data sets from several cancer types, it reaches an AUC of 0.92. Although higher levels of DNA methylation are often associated with transcriptional silencing, counter-intuitive positive statistical dependencies between DNA methylation and expression levels have been recently reported for two cancer types. Here, we re-analyze combined expression and DNA methylation data sets, comprising over 5,000 samples, and demonstrate that the conjunction of hypermethylation of bivalent chromatin and up-regulation of the corresponding genes is a general phenomenon in cancer. This up-regulation affects many developmental genes and transcription factors, including dozens of homeobox genes and other genes implicated in cancer. Thus, we reason that the disturbance of bivalent chromatin may be intimately linked to tumorigenesis

RNA secondary structure prediction from multi-aligned sequences

It has been well accepted that the RNA secondary structures of most functional non-coding RNAs (ncRNAs) are closely related to their functions and are conserved during evolution. Hence, prediction of conserved secondary structures from evolutionarily related sequences is one important task in RNA bioinformatics; the methods are useful not only to further functional analyses of ncRNAs but also to improve the accuracy of secondary structure predictions and to find novel functional RNAs from the genome. In this review, I focus on common secondary structure prediction from a given aligned RNA sequence, in which one secondary structure whose length is equal to that of the input alignment is predicted. I systematically review and classify existing tools and algorithms for the problem, by utilizing the information employed in the tools and by adopting a unified viewpoint based on maximum expected gain (MEG) estimators. I believe that this classification will allow a deeper understanding of each tool and provide users with useful information for selecting tools for common secondary structure predictions.Comment: A preprint of an invited review manuscript that will be published in a chapter of the book `Methods in Molecular Biology'. Note that this version of the manuscript may differ from the published versio

Using and Joining a Franchised Private Sector Provider Network in Myanmar

BACKGROUND: Quality is central to understanding provider motivations to join and remain within a social franchising network. Quality also appears as a key issue from the client's perspective, and may influence why a client chooses to use a franchised provider over another type of provider. The dynamic relationships between providers of social franchising clinics and clients who use these services have not been thoroughly investigated in the context of Myanmar, which has an established social franchising network. This study examines client motivations to use a Sun Quality Health network provider and provider motivations to join and remain in the Sun Quality Health network. Taken together, these two aims provide an opportunity to explore the symbiotic relationship between client satisfaction and provider incentives to increase the utilization of reproductive health care services. METHODS AND FINDINGS: Results from a series of focus group discussions with clients of reproductive health services and franchised providers shows that women chose health services provided by franchised private sector general practitioners because of its perceived higher quality, associated with the availability of effective, affordable, drugs. A key finding of the study is associated with providers. Provider focus group discussions indicate that a principle determinate for joining and remaining in the Sun Quality Health Network was serving the poor

Transfer Matrices and Partition-Function Zeros for Antiferromagnetic Potts Models. V. Further Results for the Square-Lattice Chromatic Polynomial

We derive some new structural results for the transfer matrix of square-lattice Potts models with free and cylindrical boundary conditions. In particular, we obtain explicit closed-form expressions for the dominant (at large |q|) diagonal entry in the transfer matrix, for arbitrary widths m, as the solution of a special one-dimensional polymer model. We also obtain the large-q expansion of the bulk and surface (resp. corner) free energies for the zero-temperature antiferromagnet (= chromatic polynomial) through order q^{-47} (resp. q^{-46}). Finally, we compute chromatic roots for strips of widths 9 <= m <= 12 with free boundary conditions and locate roughly the limiting curves.Comment: 111 pages (LaTeX2e). Includes tex file, three sty files, and 19 Postscript figures. Also included are Mathematica files data_CYL.m and data_FREE.m. Many changes from version 1: new material on series expansions and their analysis, and several proofs of previously conjectured results. Final version to be published in J. Stat. Phy

T‐cell prolymphocytic leukemia is associated with deregulation of oncogenic microRNAs on transcriptional and epigenetic level

Deregulation of micro(mi)-RNAs is a common mechanism in tumorigenesis. We investigated the expression of 2083 miRNAs in T-cell prolymphocytic leukemia (T-PLL). Compared to physiologic CD4+ and CD8+ T-cell subsets, 111 miRNAs were differentially expressed in T-PLL. Of these, 33 belonged to miRNA gene clusters linked to cancer. Genomic variants affecting miRNAs were infrequent with the notable exception of copy number aberrations. Remarkably, we found strong upregulation of the miR-200c/-141 cluster in T-PLL to be associated with DNA hypomethylation and active promoter marks. Our findings suggest that copy number aberrations and epigenetic changes could contribute to miRNA deregulation in T-PLL

Assessing the Utility of Thermodynamic Features for microRNA Target Prediction under Relaxed Seed and No Conservation Requirements

BACKGROUND: Many computational microRNA target prediction tools are focused on several key features, including complementarity to 5'seed of miRNAs and evolutionary conservation. While these features allow for successful target identification, not all miRNA target sites are conserved and adhere to canonical seed complementarity. Several studies have propagated the use of energy features of mRNA:miRNA duplexes as an alternative feature. However, different independent evaluations reported conflicting results on the reliability of energy-based predictions. Here, we reassess the usefulness of energy features for mammalian target prediction, aiming to relax or eliminate the need for perfect seed matches and conservation requirement. METHODOLOGY/PRINCIPAL FINDINGS: We detect significant differences of energy features at experimentally supported human miRNA target sites and at genome-wide sites of AGO protein interaction. This trend is confirmed on datasets that assay the effect of miRNAs on mRNA and protein expression changes, and a simple linear regression model leads to significant correlation of predicted versus observed expression change. Compared to 6-mer seed matches as baseline, application of our energy-based model leads to ∼3-5-fold enrichment on highly down-regulated targets, and allows for prediction of strictly imperfect targets with enrichment above baseline. CONCLUSIONS/SIGNIFICANCE: In conclusion, our results indicate significant promise for energy-based miRNA target prediction that includes a broader range of targets without having to use conservation or impose stringent seed match rules

The Hominin Sites and Paleolakes Drilling Project:High-Resolution Paleoclimate Records from the East African Rift System and Their Implications for Understanding the Environmental Context of Hominin Evolution

The possibility of a causal relationship between Earth history processes and hominin evolution in Africa has been the subject of intensive paleoanthropological research for the last 25 years. One fundamental question is: can any geohistorical processes, in particular, climatic ones, be characterized with sufficient precision to enable temporal correlation with events in hominin evolution and provide support for a possible causal mechanism for evolutionary changes? Previous attempts to link paleoclimate and hominin evolution have centered on evidence from the outcrops where the hominin fossils are found, as understanding whether and how hominin populations responded to habitat change must be examined at the local basinal scale. However, these outcrop records typically provide incomplete, low-resolution climate and environmental histories, and surface weathering often precludes the application of highly sensitive, state-of-the-art paleoenvironmental methods. Continuous and well-preserved deep-sea drill core records have provided an alternative approach to reconstructing the context of hominin evolution, but have been collected at great distances from hominin sites and typically integrate information over vast spatial scales. The goal of the Hominin Sites and Paleolakes Drilling Project (HSPDP) is to analyze climate and other Earth system dynamics using detailed paleoenvironmental data acquired through scientific drilling of lacustrine depocenters at or near six key paleoanthropological sites in Kenya and Ethiopia. This review provides an overview of a unique collaboration of paleoanthropologists and earth scientists who have joined together to explicitly explore key hypotheses linking environmental history and mammalian (including hominin) evolution and potentially develop new testable hypotheses. With a focus on continuous, high-resolution proxies at timescales relevant to both biological and cultural evolution, the HSPDP aims to dramatically expand our understanding of the environmental history of eastern Africa during a significant portion of the Late Neogene and Quaternary, and to generate useful models of long-term environmental dynamics in the regionpublishersversionPeer reviewe

ViennaRNA Package 2.0

<p>Abstract</p> <p>Background</p> <p>Secondary structure forms an important intermediate level of description of nucleic acids that encapsulates the dominating part of the folding energy, is often well conserved in evolution, and is routinely used as a basis to explain experimental findings. Based on carefully measured thermodynamic parameters, exact dynamic programming algorithms can be used to compute ground states, base pairing probabilities, as well as thermodynamic properties.</p> <p>Results</p> <p>The <monospace>ViennaRNA</monospace> Package has been a widely used compilation of RNA secondary structure related computer programs for nearly two decades. Major changes in the structure of the standard energy model, the <it>Turner 2004 </it>parameters, the pervasive use of multi-core CPUs, and an increasing number of algorithmic variants prompted a major technical overhaul of both the underlying <monospace>RNAlib</monospace> and the interactive user programs. New features include an expanded repertoire of tools to assess RNA-RNA interactions and restricted ensembles of structures, additional output information such as <it>centroid </it>structures and <it>maximum expected accuracy </it>structures derived from base pairing probabilities, or <it>z</it>-<it>scores </it>for locally stable secondary structures, and support for input in <monospace>fasta</monospace> format. Updates were implemented without compromising the computational efficiency of the core algorithms and ensuring compatibility with earlier versions.</p> <p>Conclusions</p> <p>The <monospace>ViennaRNA Package 2.0</monospace>, supporting concurrent computations <monospace>via OpenMP</monospace>, can be downloaded from <url>http://www.tbi.univie.ac.at/RNA</url>.</p