Confluence reduction for Markov automata

Markov automata are a novel formalism for specifying systems exhibiting nondeterminism, probabilistic choices and Markovian rates. Recently, the process algebra MAPA was introduced to efficiently model such systems. As always, the state space explosion threatens the analysability of the models generated by such specifications. We therefore introduce confluence reduction for Markov automata, a powerful reduction technique to keep these models small. We define the notion of confluence directly on Markov automata, and discuss how to syntactically detect confluence on the MAPA language as well. That way, Markov automata generated by MAPA specifications can be reduced on-the-fly while preserving divergence-sensitive branching bisimulation. Three case studies demonstrate the significance of our approach, with reductions in analysis time up to an order of magnitude

Use of partial least squares regression to impute SNP genotypes in Italian Cattle breeds

Background The objective of the present study was to test the ability of the partial least squares regression technique to impute genotypes from low density single nucleotide polymorphisms (SNP) panels i.e. 3K or 7K to a high density panel with 50K SNP. No pedigree information was used. Methods Data consisted of 2093 Holstein, 749 Brown Swiss and 479 Simmental bulls genotyped with the Illumina 50K Beadchip. First, a single-breed approach was applied by using only data from Holstein animals. Then, to enlarge the training population, data from the three breeds were combined and a multi-breed analysis was performed. Accuracies of genotypes imputed using the partial least squares regression method were compared with those obtained by using the Beagle software. The impact of genotype imputation on breeding value prediction was evaluated for milk yield, fat content and protein content. Results In the single-breed approach, the accuracy of imputation using partial least squares regression was around 90 and 94% for the 3K and 7K platforms, respectively; corresponding accuracies obtained with Beagle were around 85% and 90%. Moreover, computing time required by the partial least squares regression method was on average around 10 times lower than computing time required by Beagle. Using the partial least squares regression method in the multi-breed resulted in lower imputation accuracies than using single-breed data. The impact of the SNP-genotype imputation on the accuracy of direct genomic breeding values was small. The correlation between estimates of genetic merit obtained by using imputed versus actual genotypes was around 0.96 for the 7K chip. Conclusions Results of the present work suggested that the partial least squares regression imputation method could be useful to impute SNP genotypes when pedigree information is not available

Properties of Foreshocks and Aftershocks of the Non-Conservative SOC Olami-Feder-Christensen Model: Triggered or Critical Earthquakes?

Following Hergarten and Neugebauer [2002] who discovered aftershock and foreshock sequences in the Olami-Feder-Christensen (OFC) discrete block-spring earthquake model, we investigate to what degree the simple toppling mechanism of this model is sufficient to account for the properties of earthquake clustering in time and space. Our main finding is that synthetic catalogs generated by the OFC model share practically all properties of real seismicity at a qualitative level, with however significant quantitative differences. We find that OFC catalogs can be in large part described by the concept of triggered seismicity but the properties of foreshocks depend on the mainshock magnitude, in qualitative agreement with the critical earthquake model and in disagreement with simple models of triggered seismicity such as the Epidemic Type Aftershock Sequence (ETAS) model [Ogata, 1988]. Many other features of OFC catalogs can be reproduced with the ETAS model with a weaker clustering than real seismicity, i.e. for a very small average number of triggered earthquakes of first generation per mother-earthquake.Comment: revtex, 19 pages, 8 eps figure

Uncertainty quantification for kinetic models in socio-economic and life sciences

Kinetic equations play a major rule in modeling large systems of interacting particles. Recently the legacy of classical kinetic theory found novel applications in socio-economic and life sciences, where processes characterized by large groups of agents exhibit spontaneous emergence of social structures. Well-known examples are the formation of clusters in opinion dynamics, the appearance of inequalities in wealth distributions, flocking and milling behaviors in swarming models, synchronization phenomena in biological systems and lane formation in pedestrian traffic. The construction of kinetic models describing the above processes, however, has to face the difficulty of the lack of fundamental principles since physical forces are replaced by empirical social forces. These empirical forces are typically constructed with the aim to reproduce qualitatively the observed system behaviors, like the emergence of social structures, and are at best known in terms of statistical information of the modeling parameters. For this reason the presence of random inputs characterizing the parameters uncertainty should be considered as an essential feature in the modeling process. In this survey we introduce several examples of such kinetic models, that are mathematically described by nonlinear Vlasov and Fokker--Planck equations, and present different numerical approaches for uncertainty quantification which preserve the main features of the kinetic solution.Comment: To appear in "Uncertainty Quantification for Hyperbolic and Kinetic Equations

BOVITA: a first overview on genome-wide genetic diversity of Italian autochthonous cattle breeds

Analysis of genomic data is increasingly becoming part of the livestock industry and is an invaluable resource for effective management of breeding programs in small populations. The recent availability of genome-wide SNP panels allows providing background information concerning genome structure in domestic animals, opening new perspectives to livestock genetics. BOVITA was established to join local efforts and resources for the genomic characterization of Italian local cattle breeds. Despite the growing diffusion of some cosmopolite specialized breeds, several autochthonous breeds are still bred in Italy. The main aim of the BOVITA is to investigate the genomic structure of Italian local cattle breeds, to provide information on their genetic status that will be useful for the management of the genetic variability, as a contribution to biodiversity conservation and prioritization actions. A total of about 800 animals (20-32 per breed) belonging to thirty Italian cattle breeds (Agerolese, Bar\ue0-Pustertaler, Burlina, Cabannina, Calvana, Chianina, Cinisara, Garfagnina, Italian Brown, Italian Holstein, Italian Simmental, Marchigiana, Maremmana, Modenese, Modicana, Mucca Pisana, Pezzata Rossa d\u2019Oropa, Piedmontese, Pinzgau, Podolica, Pontremolese, Pustertaler, Reggiana, Rendena, Romagnola, Rossa Siciliana, Sarda, Sardo-Bruna, Sardo-Modicana and Ottonese-Varzese) and two cosmopolitan breeds (Charolaise and Limousine) genotyped with the Illumina BovineSNP50 v2 BeadChip array were collected for the analysis. The genotypes of several breeds were detected in the frame of the project, whereas for some breeds these data are derived by previous studies. The dataset will be analyzed to: study several aspects of population genetic diversity, multi-dimensional scaling plot, population structure, linkage disequilibrium, and runs of homozygosity. In addition, comparative analysis of conserved haplotypes will be conducted to identify genomic segments under selection pressure. Such information also provides important insights into the mechanisms of evolution and is useful for the annotation of significant functional genomics regions. Data analysis will also be useful to select SNPs suitable for parentage test and breed genetic traceability. The analysis of the data will pinpoint the genetic distinctiveness of Italian breeds. Moreover, the obtained results contribute to a better characterization of history and genetic structure of Italian cattle breeds

Theoretical and technological building blocks for an innovation accelerator

The scientific system that we use today was devised centuries ago and is inadequate for our current ICT-based society: the peer review system encourages conservatism, journal publications are monolithic and slow, data is often not available to other scientists, and the independent validation of results is limited. Building on the Innovation Accelerator paper by Helbing and Balietti (2011) this paper takes the initial global vision and reviews the theoretical and technological building blocks that can be used for implementing an innovation (in first place: science) accelerator platform driven by re-imagining the science system. The envisioned platform would rest on four pillars: (i) Redesign the incentive scheme to reduce behavior such as conservatism, herding and hyping; (ii) Advance scientific publications by breaking up the monolithic paper unit and introducing other building blocks such as data, tools, experiment workflows, resources; (iii) Use machine readable semantics for publications, debate structures, provenance etc. in order to include the computer as a partner in the scientific process, and (iv) Build an online platform for collaboration, including a network of trust and reputation among the different types of stakeholders in the scientific system: scientists, educators, funding agencies, policy makers, students and industrial innovators among others. Any such improvements to the scientific system must support the entire scientific process (unlike current tools that chop up the scientific process into disconnected pieces), must facilitate and encourage collaboration and interdisciplinarity (again unlike current tools), must facilitate the inclusion of intelligent computing in the scientific process, must facilitate not only the core scientific process, but also accommodate other stakeholders such science policy makers, industrial innovators, and the general public

Whole mitochondrial genomes unveil the impact of domestication on goat matrilineal variability

Background: The current extensive use of the domestic goat (Capra hircus) is the result of its medium size and high adaptability as multiple breeds. The extent to which its genetic variability was influenced by early domestication practices is largely unknown. A common standard by which to analyze maternally-inherited variability of livestock species is through complete sequencing of the entire mitogenome (mitochondrial DNA, mtDNA). Results: We present the first extensive survey of goat mitogenomic variability based on 84 complete sequences selected from an initial collection of 758 samples that represent 60 different breeds of C. hircus, as well as its wild sister species, bezoar (Capra aegagrus) from Iran. Our phylogenetic analyses dated the most recent common ancestor of C. hircus to ~460,000 years (ka) ago and identified five distinctive domestic haplogroups (A, B1, C1a, D1 and G). More than 90 % of goats examined were in haplogroup A. These domestic lineages are predominantly nested within C. aegagrus branches, diverged concomitantly at the interface between the Epipaleolithic and early Neolithic periods, and underwent a dramatic expansion starting from ~12–10 ka ago. Conclusions: Domestic goat mitogenomes descended from a small number of founding haplotypes that underwent domestication after surviving the last glacial maximum in the Near Eastern refuges. All modern haplotypes A probably descended from a single (or at most a few closely related) female C. aegagrus. Zooarchaelogical data indicate that domestication first occurred in Southeastern Anatolia. Goats accompanying the first Neolithic migration waves into the Mediterranean were already characterized by two ancestral A and C variants. The ancient separation of the C branch (~130 ka ago) suggests a genetically distinct population that could have been involved in a second event of domestication. The novel diagnostic mutational motifs defined here, which distinguish wild and domestic haplogroups, could be used to understand phylogenetic relationships among modern breeds and ancient remains and to evaluate whether selection differentially affected mitochondrial genome variants during the development of economically important breeds

Mitogenomes from Egyptian Cattle Breeds: New Clues on the Origin of Haplogroup Q and the Early Spread of Bos taurus from the Near East

Background Genetic studies support the scenario that Bos taurus domestication occurred in the Near East during the Neolithic transition about 10 thousand years (ky) ago, with the likely exception of a minor secondary event in Italy. However, despite the proven effectiveness of whole mitochondrial genome data in providing valuable information concerning the origin of taurine cattle, until now no population surveys have been carried out at the level of mitogenomes in local breeds from the Near East or surrounding areas. Egypt is in close geographic and cultural proximity to the Near East, in particular the Nile Delta region, and was one of the first neighboring areas to adopt the Neolithic package. Thus, a survey of mitogenome variation of autochthonous taurine breeds from the Nile Delta region might provide new insights on the early spread of cattle rearing outside the Near East. Methodology Using Illumina high-throughput sequencing we characterized the mitogenomes from two cattle breeds, Menofi (N = 17) and Domiaty (N = 14), from the Nile Delta region. Phylogenetic and Bayesian analyses were subsequently performed. Conclusions Phylogenetic analyses of the 31 mitogenomes confirmed the prevalence of haplogroup T1, similar to most African cattle breeds, but showed also high frequencies for haplogroups T2, T3 and Q1, and an extremely high haplotype diversity, while Bayesian skyline plots pointed to a main episode of population growth ~12.5 ky ago. Comparisons of Nile Delta mitogenomes with those from other geographic areas revealed that (i) most Egyptian mtDNAs are probably direct local derivatives from the founder domestic herds which first arrived from the Near East and the extent of gene flow from and towards the Nile Delta region was limited after the initial founding event(s); (ii) haplogroup Q1 was among these founders, thus proving that it underwent domestication in the Near East together with the founders of the T clades