Fractal properties of relaxation clusters and phase transition in a stochastic sandpile automaton

We study numerically the spatial properties of relaxation clusters in a two dimensional sandpile automaton with dynamic rules depending stochastically on a parameter p, which models the effects of static friction. In the limiting cases p=1 and p=0 the model reduces to the critical height model and critical slope model, respectively. At p=p_c, a continuous phase transition occurs to the state characterized by a nonzero average slope. Our analysis reveals that the loss of finite average slope at the transition is accompanied by the loss of fractal properties of the relaxation clusters.Comment: 11 page

ViFi: accurate detection of viral integration and mRNA fusion reveals indiscriminate and unregulated transcription in proximal genomic regions in cervical cancer.

The integration of viral sequences into the host genome is an important driver of tumorigenesis in many viral mediated cancers, notably cervical cancer and hepatocellular carcinoma. We present ViFi, a computational method that combines phylogenetic methods with reference-based read mapping to detect viral integrations. In contrast with read-based reference mapping approaches, ViFi is faster, and shows high precision and sensitivity on both simulated and biological data, even when the integrated virus is a novel strain or highly mutated. We applied ViFi to matched genomic and mRNA data from 68 cervical cancer samples from TCGA and found high concordance between the two. Surprisingly, viral integration resulted in a dramatic transcriptional upregulation in all proximal elements, including LINEs and LTRs that are not normally transcribed. This upregulation is highly correlated with the presence of a viral gene fused with a downstream human element. Moreover, genomic rearrangements suggest the formation of apparent circular extrachromosomal (ecDNA) human-viral structures. Our results suggest the presence of apparent small circular fusion viral/human ecDNA, which correlates with indiscriminate and unregulated expression of proximal genomic elements, potentially contributing to the pathogenesis of HPV-associated cervical cancers. ViFi is available at https://github.com/namphuon/ViFi

Number and Size Distribution of Colorectal Adenomas under the Multistage Clonal Expansion Model of Cancer

Colorectal cancer (CRC) is believed to arise from mutant stem cells in colonic crypts that undergo a well-characterized progression involving benign adenoma, the precursor to invasive carcinoma. Although a number of (epi)genetic events have been identified as drivers of this process, little is known about the dynamics involved in the stage-wise progression from the first appearance of an adenoma to its ultimate conversion to malignant cancer. By the time adenomas become endoscopically detectable (i.e., are in the range of 1–2 mm in diameter), adenomas are already comprised of hundreds of thousands of cells and may have been in existence for several years if not decades. Thus, a large fraction of adenomas may actually remain undetected during endoscopic screening and, at least in principle, could give rise to cancer before they are detected. It is therefore of importance to establish what fraction of adenomas is detectable, both as a function of when the colon is screened for neoplasia and as a function of the achievable detection limit. To this end, we have derived mathematical expressions for the detectable adenoma number and size distributions based on a recently developed stochastic model of CRC. Our results and illustrations using these expressions suggest (1) that screening efficacy is critically dependent on the detection threshold and implicit knowledge of the relevant stem cell fraction in adenomas, (2) that a large fraction of non-extinct adenomas remains likely undetected assuming plausible detection thresholds and cell division rates, and (3), under a realistic description of adenoma initiation, growth and progression to CRC, the empirical prevalence of adenomas is likely inflated with lesions that are not on the pathway to cancer

Lyapunov exponents and transport in the Zhang model of Self-Organized Criticality

We discuss the role played by the Lyapunov exponents in the dynamics of Zhang's model of Self-Organized Criticality. We show that a large part of the spectrum (slowest modes) is associated with the energy transpor in the lattice. In particular, we give bounds on the first negative Lyapunov exponent in terms of the energy flux dissipated at the boundaries per unit of time. We then establish an explicit formula for the transport modes that appear as diffusion modes in a landscape where the metric is given by the density of active sites. We use a finite size scaling ansatz for the Lyapunov spectrum and relate the scaling exponent to the scaling of quantities like avalanche size, duration, density of active sites, etc ...Comment: 33 pages, 6 figures, 1 table (to appear

Mapping clustered mutations in cancer reveals APOBEC3 mutagenesis of ecDNA

Clustered somatic mutations are common in cancer genomes and previous analyses reveal several types of clustered single-base substitutions, which include doublet- and multi-base substitutions1–5, diffuse hypermutation termed omikli6, and longer strand-coordinated events termed kataegis3,7–9. Here we provide a comprehensive characterization of clustered substitutions and clustered small insertions and deletions (indels) across 2,583 whole-genome-sequenced cancers from 30 types of cancer10. Clustered mutations were highly enriched in driver genes and associated with differential gene expression and changes in overall survival. Several distinct mutational processes gave rise to clustered indels, including signatures that were enriched in tobacco smokers and homologous-recombination-deficient cancers. Doublet-base substitutions were caused by at least 12 mutational processes, whereas most multi-base substitutions were generated by either tobacco smoking or exposure to ultraviolet light. Omikli events, which have previously been attributed to APOBEC3 activity6, accounted for a large proportion of clustered substitutions; however, only 16.2% of omikli matched APOBEC3 patterns. Kataegis was generated by multiple mutational processes, and 76.1% of all kataegic events exhibited mutational patterns that are associated with the activation-induced deaminase (AID) and APOBEC3 family of deaminases. Co-occurrence of APOBEC3 kataegis and extrachromosomal DNA (ecDNA), termed kyklonas (Greek for cyclone), was found in 31% of samples with ecDNA. Multiple distinct kyklonic events were observed on most mutated ecDNA. ecDNA containing known cancer genes exhibited both positive selection and kyklonic hypermutation. Our results reveal the diversity of clustered mutational processes in human cancer and the role of APOBEC3 in recurrently mutating and fuelling the evolution of ecDNA

Targeted and Untargeted Approaches Unravel Novel Candidate Genes and Diagnostic SNPs for Quantitative Resistance of the Potato (Solanum tuberosum L.) to Phytophthora infestans Causing the Late Blight Disease

The oomycete Phytophthora infestans causes late blight of potato, which can completely destroy the crop. Therefore, for the past 160 years, late blight has been the most important potato disease worldwide. The identification of cultivars with high and durable field resistance to P. infestans is an objective of most potato breeding programs. This type of resistance is polygenic and therefore quantitative. Its evaluation requires multi-year and location trials. Furthermore, quantitative resistance to late blight correlates with late plant maturity, a negative agricultural trait. Knowledge of the molecular genetic basis of quantitative resistance to late blight not compromised by late maturity is very limited. It is however essential for developing diagnostic DNA markers that facilitate the efficient combination of superior resistance alleles in improved cultivars. We used association genetics in a population of 184 tetraploid potato cultivars in order to identify single nucleotide polymorphisms (SNPs) that are associated with maturity corrected resistance (MCR) to late blight. The population was genotyped for almost 9000 SNPs from three different sources. The first source was candidate genes specifically selected for their function in the jasmonate pathway. The second source was novel candidate genes selected based on comparative transcript profiling (RNA-Seq) of groups of genotypes with contrasting levels of quantitative resistance to P. infestans. The third source was the first generation 8.3k SolCAP SNP genotyping array available in potato for genome wide association studies (GWAS). Twenty seven SNPs from all three sources showed robust association with MCR. Some of those were located in genes that are strong candidates for directly controlling quantitative resistance, based on functional annotation. Most important were: a lipoxygenase (jasmonate pathway), a 3-hydroxy-3-methylglutaryl coenzyme A reductase (mevalonate pathway), a P450 protein (terpene biosynthesis), a transcription factor and a homolog of a major gene for resistance to P. infestans from the wild potato species Solanum venturii. The candidate gene approach and GWAS complemented each other as they identified different genes. The results of this study provide new insight in the molecular genetic basis of quantitative resistance in potato and a toolbox of diagnostic SNP markers for breeding applications.</p

Combinatorial CRISPR-Cas9 screens for de novo mapping of genetic interactions.

We developed a systematic approach to map human genetic networks by combinatorial CRISPR-Cas9 perturbations coupled to robust analysis of growth kinetics. We targeted all pairs of 73 cancer genes with dual guide RNAs in three cell lines, comprising 141,912 tests of interaction. Numerous therapeutically relevant interactions were identified, and these patterns replicated with combinatorial drugs at 75% precision. From these results, we anticipate that cellular context will be critical to synthetic-lethal therapies

Systems biological and mechanistic modelling of radiation-induced cancer

This paper summarises the five presentations at the First International Workshop on Systems Radiation Biology that were concerned with mechanistic models for carcinogenesis. The mathematical description of various hypotheses about the carcinogenic process, and its comparison with available data is an example of systems biology. It promises better understanding of effects at the whole body level based on properties of cells and signalling mechanisms between them. Of these five presentations, three dealt with multistage carcinogenesis within the framework of stochastic multistage clonal expansion models, another presented a deterministic multistage model incorporating chromosomal aberrations and neoplastic transformation, and the last presented a model of DNA double-strand break repair pathways for second breast cancers following radiation therapy

A Heuristic Solution of the Identifiability Problem of the Age-Period-Cohort Analysis of Cancer Occurrence: Lung Cancer Example

Background: The Age–Period–Cohort (APC) analysis is aimed at estimating the following effects on disease incidence: (i) the age of the subject at the time of disease diagnosis; (ii) the time period, when the disease occurred; and (iii) the date of birth of the subject. These effects can help in evaluating the biological events leading to the disease, in estimating the influence of distinct risk factors on disease occurrence, and in the development of new strategies for disease prevention and treatment. Methodology/Principal Findings: We developed a novel approach for estimating the APC effects on disease incidence rates in the frame of the Log-Linear Age-Period-Cohort (LLAPC) model. Since the APC effects are linearly interdependent and cannot be uniquely estimated, solving this identifiability problem requires setting four redundant parameters within a set of unknown parameters. By setting three parameters (one of the time-period and the birth-cohort effects and the corresponding age effect) to zero, we reduced this problem to the problem of determining one redundant parameter and, used as such, the effect of the time-period adjacent to the anchored time period. By varying this identification parameter, a family of estimates of the APC effects can be obtained. Using a heuristic assumption that the differences between the adjacent birth-cohort effects are small, we developed a numerical method for determining the optimal value of the identification parameter, by which a unique set of all APC effects is determined and the identifiability problem is solved