Computerized optimization of elastic booster autopilots. Volume 1: Technical manual

The philosophy and the mathematical basis of the nonlinear programming algorithm underlying the development of the COEBRA program were given. A User's Manual was given in a separate document. The purpose of this work was to convert the COEBRA program from the CDC 6400/6500 digital computer system to the UNIVAC 1108 at the George C. Marshall Space Flight Center and to provide an instruction manual on the use of the program

Lateral gene transfer, rearrangement, reconciliation

Background. Models of ancestral gene order reconstruction have progressively integrated different evolutionary patterns and processes such as unequal gene content, gene duplications, and implicitly sequence evolution via reconciled gene trees. These models have so far ignored lateral gene transfer, even though in unicellular organisms it can have an important confounding effect, and can be a rich source of information on the function of genes through the detection of transfers of clusters of genes. Result. We report an algorithm together with its implementation, DeCoLT, that reconstructs ancestral genome organization based on reconciled gene trees which summarize information on sequence evolution, gene origination, duplication, loss, and lateral transfer. DeCoLT optimizes in polynomial time on the number of rearrangements, computed as the number of gains and breakages of adjacencies between pairs of genes. We apply DeCoLT to 1099 gene families from 36 cyanobacteria genomes. Conclusion. DeCoLT is able to reconstruct adjacencies in 35 ancestral bacterial genomes with a thousand gene families in a few hours, and detects clusters of co-transferred genes. DeCoLT may also be used with any relationship between genes instead of adjacencies, to reconstruct ancestral interactions, functions or complexes. Availability. http://pbil.univ-lyon1.fr/software/DeCoLT

Congruent evolution of genetic and environmental robustness in microRNA

Genetic robustness, the preservation of an optimal phenotype in the face of mutations, is critical to the understanding of evolution as phenotypically expressed genetic variation is the fuel of natural selection. The origin of genetic robustness, whether it evolves directly by natural selection or it is a correlated byproduct of other phenotypic traits, is, however, unresolved. Examining microRNA (miRNA) genes of several eukaryotic species, Borenstein and Ruppin (Borenstein et al. 2006, PNAS 103: 6593), showed that the structure of miRNA precursor stem-loops exhibits significantly increased mutational robustness in comparison with a sample of random RNA sequences with the same stem-loop structure. The observed robustness was found to be uncorrelated with traditional measures of environmental robustness -- implying that miRNA sequences show evidence of the direct evolution of genetic robustness. These findings are surprising as theoretical results indicate that the direct evolution of robustness requires high mutation rates and/or large effective population sizes only found among RNA viruses, not multicellular eukaryotes. Introducing a novel measure of environmental robustness based on the equilibrium thermodynamic ensemble of secondary structures of the miRNA precursor sequences we demonstrate that the biophysics of RNA folding, induces a high level of correlation between genetic (mutational) and environmental (thermodynamic) robustness, as expected from the theory of plastogenetic congruence introduced by Ancel and Fontana (Ancel et al. 2000, J. Exp. Zool. 288: 242). In light of theoretical considerations we believe that this correlation strongly suggests that genetic robustness observed in miRNA sequences is the byproduct of selection for environmental robustness.Comment: Accepted for publication in Mol. Biol. Evol. Supplemental Information available as a separate pdf file from http://angel.elte.hu:/~ssolo/miRNA_supp_mat.pd

Unifying Parsimonious Tree Reconciliation

Evolution is a process that is influenced by various environmental factors, e.g. the interactions between different species, genes, and biogeographical properties. Hence, it is interesting to study the combined evolutionary history of multiple species, their genes, and the environment they live in. A common approach to address this research problem is to describe each individual evolution as a phylogenetic tree and construct a tree reconciliation which is parsimonious with respect to a given event model. Unfortunately, most of the previous approaches are designed only either for host-parasite systems, for gene tree/species tree reconciliation, or biogeography. Hence, a method is desirable, which addresses the general problem of mapping phylogenetic trees and covering all varieties of coevolving systems, including e.g., predator-prey and symbiotic relationships. To overcome this gap, we introduce a generalized cophylogenetic event model considering the combinatorial complete set of local coevolutionary events. We give a dynamic programming based heuristic for solving the maximum parsimony reconciliation problem in time O(n^2), for two phylogenies each with at most n leaves. Furthermore, we present an exact branch-and-bound algorithm which uses the results from the dynamic programming heuristic for discarding partial reconciliations. The approach has been implemented as a Java application which is freely available from http://pacosy.informatik.uni-leipzig.de/coresym.Comment: Peer-reviewed and presented as part of the 13th Workshop on Algorithms in Bioinformatics (WABI2013

Some gating potentiators, including VX-770, diminish ΔF508-CFTR functional expression.

Cystic fibrosis (CF) is caused by mutations in the CF transmembrane regulator (CFTR) that result in reduced anion conductance at the apical membrane of secretory epithelia. Treatment of CF patients carrying the G551D gating mutation with the potentiator VX-770 (ivacaftor) largely restores channel activity and has shown substantial clinical benefit. However, most CF patients carry the ΔF508 mutation, which impairs CFTR folding, processing, function, and stability. Studies in homozygous ΔF508 CF patients indicated little clinical benefit of monotherapy with the investigational corrector VX-809 (lumacaftor) or VX-770, whereas combination clinical trials show limited but significant improvements in lung function. We show that VX-770, as well as most other potentiators, reduces the correction efficacy of VX-809 and another investigational corrector, VX-661. To mimic the administration of VX-770 alone or in combination with VX-809, we examined its long-term effect in immortalized and primary human respiratory epithelia. VX-770 diminished the folding efficiency and the metabolic stability of ΔF508-CFTR at the endoplasmic reticulum (ER) and post-ER compartments, respectively, causing reduced cell surface ΔF508-CFTR density and function. VX-770-induced destabilization of ΔF508-CFTR was influenced by second-site suppressor mutations of the folding defect and was prevented by stabilization of the nucleotide-binding domain 1 (NBD1)-NBD2 interface. The reduced correction efficiency of ΔF508-CFTR, as well as of two other processing mutations in the presence of VX-770, suggests the need for further optimization of potentiators to maximize the clinical benefit of corrector-potentiator combination therapy in CF

Optimal tumor sampling for immunostaining of biomarkers in breast carcinoma

IntroductionBiomarkers, such as Estrogen Receptor, are used to determine therapy and prognosis in breast carcinoma. Immunostaining assays of biomarker expression have a high rate of inaccuracy; for example, estimates are as high as 20% for Estrogen Receptor. Biomarkers have been shown to be heterogeneously expressed in breast tumors and this heterogeneity may contribute to the inaccuracy of immunostaining assays. Currently, no evidence-based standards exist for the amount of tumor that must be sampled in order to correct for biomarker heterogeneity. The aim of this study was to determine the optimal number of 20X fields that are necessary to estimate a representative measurement of expression in a whole tissue section for selected biomarkers: ER, HER-2, AKT, ERK, S6K1, GAPDH, Cytokeratin, and MAP-Tau.MethodsTwo collections of whole tissue sections of breast carcinoma were immunostained for biomarkers. Expression was quantified using the Automated Quantitative Analysis (AQUA) method of quantitative immunofluorescence. Simulated sampling of various numbers of fields (ranging from one to thirty five) was performed for each marker. The optimal number was selected for each marker via resampling techniques and minimization of prediction error over an independent test set.ResultsThe optimal number of 20X fields varied by biomarker, ranging between three to fourteen fields. More heterogeneous markers, such as MAP-Tau protein, required a larger sample of 20X fields to produce representative measurement.ConclusionsThe optimal number of 20X fields that must be sampled to produce a representative measurement of biomarker expression varies by marker with more heterogeneous markers requiring a larger number. The clinical implication of these findings is that breast biopsies consisting of a small number of fields may be inadequate to represent whole tumor biomarker expression for many markers. Additionally, for biomarkers newly introduced into clinical use, especially if therapeutic response is dictated by level of expression, the optimal size of tissue sample must be determined on a marker-by-marker basis

FISH and immunohistochemical status of the hepatocyte growth factor receptor (c-Met) in 184 invasive breast tumors

Grants PI05/0961 and PI06/1513 from Ministerio de Sanidad y Consumo ISCIII and RTICC 06/0020/1