Computation of the Transient in Max-Plus Linear Systems via SMT-Solving

This paper proposes a new approach, grounded in Satisfiability Modulo Theories (SMT), to study the transient of a Max-Plus Linear (MPL) system, that is the number of steps leading to its periodic regime. Differently from state-of-the-art techniques, our approach allows the analysis of periodic behaviors for subsets of initial states, as well as the characterization of sets of initial states exhibiting the same specific periodic behavior and transient. Our experiments show that the proposed technique dramatically outperforms state-of-the-art methods based on max-plus algebra computations for systems of large dimensions.Comment: The paper consists of 22 pages (including references and Appendix). It is accepted in FORMATS 2020 First revisio

Megabank found? Flanks record sea level

On Leg 101, the first international voyage for the Ocean Drilling Program, the deep-sea drilling ship JOIOES Resolution (SEDCO/BP 471) left Miami, Fla., on Jan. 31 to investigate the geology of the Bahamas. (Leg 100 tested the Resolution's readiness. See July Geotimes.) Before returning to Miami on March 14, the crew had drilled 19 holes al 11 sites and recovered 46.2% of the cored section (about 1.5 of 3.1 km cored). The scientific party wanted to test conflicting hypotheses about the development of the modern shallow water carbonate banks and intervening deep -water throughs in the Bahamas, and to study the growth patterns of carbonate slopes and their response to sea-level fluctuations. Those objectives (the 'deep ' and the 'shallow') were selected beause recent advances in interpreting the micropaleontology of shallow-water carbonate platforms, coupled with data from previous sedimentological investigations and regional and site-specific seismic surveys, now permit consistent stratigraphic comparisons in the Bahamas

The systemic lupus erythematosus IRF5 risk haplotype is associated with systemic sclerosis

Systemic sclerosis (SSc) is a fibrotic autoimmune disease in which the genetic component plays an important role. One of the strongest SSc association signals outside the human leukocyte antigen (HLA) region corresponds to interferon (IFN) regulatory factor 5 (IRF5), a major regulator of the type I IFN pathway. In this study we aimed to evaluate whether three different haplotypic blocks within this locus, which have been shown to alter the protein function influencing systemic lupus erythematosus (SLE) susceptibility, are involved in SSc susceptibility and clinical phenotypes. For that purpose, we genotyped one representative single-nucleotide polymorphism (SNP) of each block (rs10488631, rs2004640, and rs4728142) in a total of 3,361 SSc patients and 4,012 unaffected controls of Caucasian origin from Spain, Germany, The Netherlands, Italy and United Kingdom. A meta-analysis of the allele frequencies was performed to analyse the overall effect of these IRF5 genetic variants on SSc. Allelic combination and dependency tests were also carried out. The three SNPs showed strong associations with the global disease (rs4728142: P = 1.34×10<sup>−8</sup>, OR = 1.22, CI 95% = 1.14–1.30; rs2004640: P = 4.60×10<sup>−7</sup>, OR = 0.84, CI 95% = 0.78–0.90; rs10488631: P = 7.53×10<sup>−20</sup>, OR = 1.63, CI 95% = 1.47–1.81). However, the association of rs2004640 with SSc was not independent of rs4728142 (conditioned P = 0.598). The haplotype containing the risk alleles (rs4728142*A-rs2004640*T-rs10488631*C: P = 9.04×10<sup>−22</sup>, OR = 1.75, CI 95% = 1.56–1.97) better explained the observed association (likelihood P-value = 1.48×10<sup>−4</sup>), suggesting an additive effect of the three haplotypic blocks. No statistical significance was observed in the comparisons amongst SSc patients with and without the main clinical characteristics. Our data clearly indicate that the SLE risk haplotype also influences SSc predisposition, and that this association is not sub-phenotype-specific

A genome-wide association study follow-up suggests a possible role for PPARG in systemic sclerosis susceptibility

Introduction: A recent genome-wide association study (GWAS) comprising a French cohort of systemic sclerosis (SSc) reported several non-HLA single-nucleotide polymorphisms (SNPs) showing a nominal association in the discovery phase. We aimed to identify previously overlooked susceptibility variants by using a follow-up strategy.<p></p> Methods: Sixty-six non-HLA SNPs showing a P value <10-4 in the discovery phase of the French SSc GWAS were analyzed in the first step of this study, performing a meta-analysis that combined data from the two published SSc GWASs. A total of 2,921 SSc patients and 6,963 healthy controls were included in this first phase. Two SNPs, PPARG rs310746 and CHRNA9 rs6832151, were selected for genotyping in the replication cohort (1,068 SSc patients and 6,762 healthy controls) based on the results of the first step. Genotyping was performed by using TaqMan SNP genotyping assays. Results: We observed nominal associations for both PPARG rs310746 (PMH = 1.90 × 10-6, OR, 1.28) and CHRNA9 rs6832151 (PMH = 4.30 × 10-6, OR, 1.17) genetic variants with SSc in the first step of our study. In the replication phase, we observed a trend of association for PPARG rs310746 (P value = 0.066; OR, 1.17). The combined overall Mantel-Haenszel meta-analysis of all the cohorts included in the present study revealed that PPARG rs310746 remained associated with SSc with a nominal non-genome-wide significant P value (PMH = 5.00 × 10-7; OR, 1.25). No evidence of association was observed for CHRNA9 rs6832151 either in the replication phase or in the overall pooled analysis.<p></p> Conclusion: Our results suggest a role of PPARG gene in the development of SSc

Exceptional preservation of palaeozoic steroids in a diagenetic continuum

The occurrence of intact sterols has been restricted to immature Cretaceous (~125 Ma) sediments with one report from the Late Jurassic (~165 Ma). Here we report the oldest occurrence of intact sterols in a Crustacean fossil preserved for ca. 380 Ma within a Devonian concretion. The exceptional preservation of the biomass is attributed to microbially induced carbonate encapsulation, preventing full decomposition and transformation thus extending sterol occurrences in the geosphere by 250 Ma. A suite of diagenetic transformation products of sterols was also identified in the concretion, demonstrating the remarkable coexistence of biomolecules and geomolecules in the same sample. Most importantly the original biolipids were found to be the most abundant steroids in the sample. We attribute the coexistence of steroids in a diagenetic continuum-ranging from stenols to triaromatic steroids-to microbially mediated eogenetic processes

Measuring Global Credibility with Application to Local Sequence Alignment

Computational biology is replete with high-dimensional (high-D) discrete prediction and inference problems, including sequence alignment, RNA structure prediction, phylogenetic inference, motif finding, prediction of pathways, and model selection problems in statistical genetics. Even though prediction and inference in these settings are uncertain, little attention has been focused on the development of global measures of uncertainty. Regardless of the procedure employed to produce a prediction, when a procedure delivers a single answer, that answer is a point estimate selected from the solution ensemble, the set of all possible solutions. For high-D discrete space, these ensembles are immense, and thus there is considerable uncertainty. We recommend the use of Bayesian credibility limits to describe this uncertainty, where a (1−α)%, 0≤α≤1, credibility limit is the minimum Hamming distance radius of a hyper-sphere containing (1−α)% of the posterior distribution. Because sequence alignment is arguably the most extensively used procedure in computational biology, we employ it here to make these general concepts more concrete. The maximum similarity estimator (i.e., the alignment that maximizes the likelihood) and the centroid estimator (i.e., the alignment that minimizes the mean Hamming distance from the posterior weighted ensemble of alignments) are used to demonstrate the application of Bayesian credibility limits to alignment estimators. Application of Bayesian credibility limits to the alignment of 20 human/rodent orthologous sequence pairs and 125 orthologous sequence pairs from six Shewanella species shows that credibility limits of the alignments of promoter sequences of these species vary widely, and that centroid alignments dependably have tighter credibility limits than traditional maximum similarity alignments

Boolean dynamics revisited through feedback interconnections

Boolean models of physical or biological systems describe the global dynamics of the system and their attractors typically represent asymptotic behaviors. In the case of large networks composed of several modules, it may be difficult to identify all the attractors. To explore Boolean dynamics from a novel viewpoint, we will analyse the dynamics emerging from the composition of two known Boolean modules. The state transition graphs and attractors for each of the modules can be combined to construct a new asymptotic graph which will (1) provide a reliable method for attractor computation with partial information; (2) illustrate the differences in dynamical behavior induced by the updating strategy (asynchronous, synchronous, or mixed); and (3) show the inherited organization/structure of the original network’s state transition graph.publishe

The Insulator Protein SU(HW) Fine-Tunes Nuclear Lamina Interactions of the Drosophila Genome

Specific interactions of the genome with the nuclear lamina (NL) are thought to assist chromosome folding inside the nucleus and to contribute to the regulation of gene expression. High-resolution mapping has recently identified hundreds of large, sharply defined lamina-associated domains (LADs) in the human genome, and suggested that the insulator protein CTCF may help to demarcate these domains. Here, we report the detailed structure of LADs in Drosophila cells, and investigate the putative roles of five insulator proteins in LAD organization. We found that the Drosophila genome is also organized in discrete LADs, which are about five times smaller than human LADs but contain on average a similar number of genes. Systematic comparison to new and published insulator binding maps shows that only SU(HW) binds preferentially at LAD borders and at specific positions inside LADs, while GAF, CTCF, BEAF-32 and DWG are mostly absent from these regions. By knockdown and overexpression studies we demonstrate that SU(HW) weakens genome – NL interactions through a local antagonistic effect, but we did not obtain evidence that it is essential for border formation. Our results provide insights into the evolution of LAD organization and identify SU(HW) as a fine-tuner of genome – NL interactions