Global parameter identification of stochastic reaction networks from single trajectories

We consider the problem of inferring the unknown parameters of a stochastic biochemical network model from a single measured time-course of the concentration of some of the involved species. Such measurements are available, e.g., from live-cell fluorescence microscopy in image-based systems biology. In addition, fluctuation time-courses from, e.g., fluorescence correlation spectroscopy provide additional information about the system dynamics that can be used to more robustly infer parameters than when considering only mean concentrations. Estimating model parameters from a single experimental trajectory enables single-cell measurements and quantification of cell--cell variability. We propose a novel combination of an adaptive Monte Carlo sampler, called Gaussian Adaptation, and efficient exact stochastic simulation algorithms that allows parameter identification from single stochastic trajectories. We benchmark the proposed method on a linear and a non-linear reaction network at steady state and during transient phases. In addition, we demonstrate that the present method also provides an ellipsoidal volume estimate of the viable part of parameter space and is able to estimate the physical volume of the compartment in which the observed reactions take place.Comment: Article in print as a book chapter in Springer's "Advances in Systems Biology

Study protocol of the LARK (TROG 17.03) clinical trial: a phase II trial investigating the dosimetric impact of Liver Ablative Radiotherapy using Kilovoltage intrafraction monitoring

BACKGROUND: Stereotactic Ablative Body Radiotherapy (SABR) is a non-invasive treatment which allows delivery of an ablative radiation dose with high accuracy and precision. SABR is an established treatment for both primary and secondary liver malignancies, and technological advances have improved its efficacy and safety. Respiratory motion management to reduce tumour motion and image guidance to achieve targeting accuracy are crucial elements of liver SABR. This phase II multi-institutional TROG 17.03 study, Liver Ablative Radiotherapy using Kilovoltage intrafraction monitoring (LARK), aims to investigate and assess the dosimetric impact of the KIM real-time image guidance technology. KIM utilises standard linear accelerator equipment and therefore has the potential to be a widely available real-time image guidance technology for liver SABR. METHODS: Forty-six patients with either hepatocellular carcinoma or oligometastatic disease to the liver suitable for and treated with SABR using Kilovoltage Intrafraction Monitoring (KIM) guidance will be included in the study. The dosimetric impact will be assessed by quantifying accumulated patient dose distribution with or without the KIM intervention. The patient treatment outcomes of local control, toxicity and quality of life will be measured. DISCUSSION: Liver SABR is a highly effective treatment, but precise dose delivery is challenging due to organ motion. Currently, there is a lack of widely available options for performing real-time tumour localisation to assist with accurate delivery of liver SABR. This study will provide an assessment of the impact of KIM as a potential solution for real-time image guidance in liver SABR. TRIAL REGISTRATION: This trial was registered on December 7th 2016 on ClinicalTrials.gov under the trial-ID NCT02984566

A unified framework for multi-locus association analysis of both common and rare variants

<p>Abstract</p> <p>Background</p> <p>Common, complex diseases are hypothesized to result from a combination of common and rare genetic variants. We developed a unified framework for the joint association testing of both types of variants. Within the framework, we developed a union-intersection test suitable for genome-wide analysis of single nucleotide polymorphisms (SNPs), candidate gene data, as well as medical sequencing data. The union-intersection test is a composite test of association of genotype frequencies and differential correlation among markers.</p> <p>Results</p> <p>We demonstrated by computer simulation that the false positive error rate was controlled at the expected level. We also demonstrated scenarios in which the multi-locus test was more powerful than traditional single marker analysis. To illustrate use of the union-intersection test with real data, we analyzed a publically available data set of 319,813 autosomal SNPs genotyped for 938 cases of Parkinson disease and 863 neurologically normal controls for which no genome-wide significant results were found by traditional single marker analysis. We also analyzed an independent follow-up sample of 183 cases and 248 controls for replication.</p> <p>Conclusions</p> <p>We identified a single risk haplotype with a directionally consistent effect in both samples in the gene <it>GAK</it>, which is involved in clathrin-mediated membrane trafficking. We also found suggestive evidence that directionally inconsistent marginal effects from single marker analysis appeared to result from risk being driven by different haplotypes in the two samples for the genes <it>SYN3 </it>and <it>NGLY1</it>, which are involved in neurotransmitter release and proteasomal degradation, respectively. These results illustrate the utility of our unified framework for genome-wide association analysis of common, complex diseases.</p

Environmental variables, habitat discontinuity and life history shaping the genetic structure of Pomatoschistus marmoratus

Coastal lagoons are semi-isolated ecosystems exposed to wide fluctuations of environmental conditions and showing habitat fragmentation. These features may play an important role in separating species into different populations, even at small spatial scales. In this study, we evaluate the concordance between mitochondrial (previous published data) and nuclear data analyzing the genetic variability of Pomatoschistus marmoratus in five localities, inside and outside the Mar Menor coastal lagoon (SE Spain) using eight microsatellites. High genetic diversity and similar levels of allele richness were observed across all loci and localities, although significant genic and genotypic differentiation was found between populations inside and outside the lagoon. In contrast to the FST values obtained from previous mitochondrial DNA analyses (control region), the microsatellite data exhibited significant differentiation among samples inside the Mar Menor and between lagoonal and marine samples. This pattern was corroborated using Cavalli-Sforza genetic distances. The habitat fragmentation inside the coastal lagoon and among lagoon and marine localities could be acting as a barrier to gene flow and contributing to the observed genetic structure. Our results from generalized additive models point a significant link between extreme lagoonal environmental conditions (mainly maximum salinity) and P. marmoratus genetic composition. Thereby, these environmental features could be also acting on genetic structure of coastal lagoon populations of P. marmoratus favoring their genetic divergence. The mating strategy of P. marmoratus could be also influencing our results obtained from mitochondrial and nuclear DNA. Therefore, a special consideration must be done in the selection of the DNA markers depending on the reproductive strategy of the species

An increased fraction of circulating miR-363 and miR-16 is particle bound in patients with chronic lymphocytic leukaemia as compared to normal subjects.

In vitro culture studies have shown that miR-363 is enriched in extracellular vesicles from chronic lymphocytic leukaemia cells. We wondered whether miR-363 was detectable in plasma, which is an essential precondition for further studies to assess its usefulness as a biomarker. Using samples from two clinical trials: one enrolling patients with advanced disease and the other asymptomatic patients with early stage disease, we determined plasma miR-363 levels and secondly investigated the distribution of this miRNA between plasma and particle bound fractions in patients and normal subjects.Advanced disease (n = 95) was associated with higher levels of miR-363 than early stage disease (n = 45) or normal subjects (n = 11) but there was no association with markers of prognosis. The distribution of specific miRNA between particle bound and plasma protein fractions was investigated using size exclusion chromatography on plasma from patients (n = 4) and normal subjects (n = 3). ~ 20% of total miR-16 and miR-363 is particle bound in patients while there was no detectable particle bound material in normal subjects. Our work demonstrates that miR-363 levels are raised in chronic lymphocytic leukaemia patients and raises the possibility that distribution of circulating miRNA between plasma fractions differs in health and disease

Potential climatic transitions with profound impact on Europe

We discuss potential transitions of six climatic subsystems with large-scale impact on Europe, sometimes denoted as tipping elements. These are the ice sheets on Greenland and West Antarctica, the Atlantic thermohaline circulation, Arctic sea ice, Alpine glaciers and northern hemisphere stratospheric ozone. Each system is represented by co-authors actively publishing in the corresponding field. For each subsystem we summarize the mechanism of a potential transition in a warmer climate along with its impact on Europe and assess the likelihood for such a transition based on published scientific literature. As a summary, the ‘tipping’ potential for each system is provided as a function of global mean temperature increase which required some subjective interpretation of scientific facts by the authors and should be considered as a snapshot of our current understanding. <br/

The influence of tumor size and environment on gene expression in commonly used human tumor lines

BACKGROUND: The expression profiles of solid tumor models in rodents have been only minimally studied despite their extensive use to develop anticancer agents. We have applied RNA expression profiling using Affymetrix U95A GeneChips to address fundamental biological questions about human tumor lines. METHODS: To determine whether gene expression changed significantly as a tumor increased in size, we analyzed samples from two human colon carcinoma lines (Colo205 and HCT-116) at three different sizes (200 mg, 500 mg and 1000 mg). To investigate whether gene expression was influenced by the strain of mouse, tumor samples isolated from C.B-17 SCID and Nu/Nu mice were also compared. Finally, the gene expression differences between tissue culture and in vivo samples were investigated by comparing profiles from lines grown in both environments. RESULTS: Multidimensional scaling and analysis of variance demonstrated that the tumor lines were dramatically different from each other and that gene expression remained constant as the tumors increased in size. Statistical analysis revealed that 63 genes were differentially expressed due to the strain of mouse the tumor was grown in but the function of the encoded proteins did not link to any distinct biological pathways. Hierarchical clustering of tissue culture and xenograft samples demonstrated that for each individual tumor line, the in vivo and in vitro profiles were more similar to each other than any other profile. We identified 36 genes with a pattern of high expression in xenograft samples that encoded proteins involved in extracellular matrix, cell surface receptors and transcription factors. An additional 17 genes were identified with a pattern of high expression in tissue culture samples and encoded proteins involved in cell division, cell cycle and RNA production. CONCLUSIONS: The environment a tumor line is grown in can have a significant effect on gene expression but tumor size has little or no effect for subcutaneously grown solid tumors. Furthermore, an individual tumor line has an RNA expression pattern that clearly defines it from other lines even when grown in different environments. This could be used as a quality control tool for preclinical oncology studies

Duplication of 7q34 is specific to juvenile pilocytic astrocytomas and a hallmark of cerebellar and optic pathway tumours

BACKGROUND: Juvenile pilocytic astrocytomas (JPA), a subgroup of low-grade astrocytomas (LGA), are common, heterogeneous and poorly understood subset of brain tumours in children. Chromosomal 7q34 duplication leading to fusion genes formed between KIAA1549 and BRAF and subsequent constitutive activation of BRAF was recently identified in a proportion of LGA, and may be involved in their pathogenesis. Our aim was to investigate additional chromosomal unbalances in LGA and whether incidence of 7q34 duplication is associated with tumour type or location. METHODS AND RESULTS: Using Illumina-Human-Hap300-Duo and 610-Quad high-resolution-SNP-based arrays and quantitative PCR on genes of interest, we investigated 84 paediatric LGA. We demonstrate that 7q34 duplication is specific to sporadic JPA (35 of 53-66%) and does not occur in other LGA subtypes (0 of 27) or NFI-associated-JPA (0 of 4). We also establish that it is site specific as it occurs in the majority of cerebellar JPA (24 of 30-80%) followed by brainstem, hypothalamic/optic pathway JPA (10 of 16-62.5%) and is rare in hemispheric JPA (1 of 7-14%). The MAP-kinase pathway, assessed through ERK phosphorylation, was active in all tumours regardless of 7q34 duplication. Gain of function studies performed on hTERT-immortalised astrocytes show that overexpression of wild-type BRAF does not increase cell proliferation or baseline MAPK signalling even if it sensitises cells to EGFR stimulation. CONCLUSIONS AND INTERPRETATION: Our results suggest that variants of JPA might arise from a unique site-restricted progenitor cell where 7q34 duplication, a hallmark of this tumour-type in association to MAPK-kinase pathway activation, potentially plays a site-specific role in their pathogenesis. Importantly, gain of function abnormalities in components of MAP-Kinase signalling are potentially present in all JPA making this tumour amenable to therapeutic targeting of this pathway. British Journal of Cancer (2009) 101, 722-733. doi: 10.1038/sj.bjc.6605179 www.bjcancer.com Published online 14 July 2009 (C) 2009 Cancer Research U

Retracing the history and planning the future of the red squirrel (Sciurus vulgaris) in Ireland using non-invasive genetics

The Eurasian red squirrel’s (Sciurus vulgaris) history in Ireland is largely unknown, but the original population is thought to have been driven to extinction by humans in the 17th Century, and multiple records exist for its subsequent reintroduction in the 19th 4 Century. However, it is currently unknown how these reintroductions affect the red squirrel population today, or may do so in the future. In this study, we report on the development of a DNA toolkit for the non-invasive genetic study of the red squirrel. Non-invasively collected red squirrel samples were combined with other samples collected throughout Ireland and previously published mitochondrial DNA (mtDNA) data from Ireland, Great Britain and continental Europe to give an insight into population genetics and historical introductions of the red squirrel in Ireland. Our findings demonstrate that the Irish red squirrel population is on a national scale quite genetically diverse, but at a local level contains relatively low levels of genetic diversity and evidence of genetic structure. This is likely an artefact of the introduction of a small number of genetically similar animals to specific sites. A lack of continuous woodland cover in Ireland has prevented further mixing with animals of different origins that may have been introduced even to neighbouring sites. Consequently, some of these genetically isolated populations are or may in the future be at risk of extinction. The Irish red squirrel population contains mtDNA haplotypes of both a British and Continental European origin, the former of which are now extinct or simply not recorded in contemporary Great Britain. The Irish population is therefore important in terms of red squirrel conservation not only in Ireland, but also for Great Britain, and should be appropriately managed

Allelic spectrum of the natural variation in CRP

With the recent completion of the International HapMap Project, many tools are in hand for genetic association studies seeking to test the common variant/common disease hypothesis. In contrast, very few tools and resources are in place for genotype–phenotype studies hypothesizing that rare variation has a large impact on the phenotype of interest. To create these tools for rare variant/common disease studies, much interest is being generated towards investing in re-sequencing either large sample sizes of random chromosomes or smaller sample sizes of patients with extreme phenotypes. As a case study for rare variant discovery in random chromosomes, we have re-sequenced ~1,000 chromosomes representing diverse populations for the gene C-reactive protein (CRP). CRP is an important gene in the fields of cardiovascular and inflammation genetics, and its size (~2 kb) makes it particularly amenable medical or deep re-sequencing. With these data, we explore several issues related to the present-day candidate gene association study including the benefits of complete SNP discovery, the effects of tagSNP selection across diverse populations, and completeness of dbSNP for CRP. Also, we show that while deep re-sequencing uncovers potentially medically relevant coding SNPs, these SNPs are fleetingly rare when genotyped in a population-based survey of 7,000 Americans (NHANES III). Collectively, these data suggest that several different types re-sequencing and genotyping approaches may be required to fully understand the complete spectrum of alleles that impact human phenotypes. ELECTRONIC SUPPLEMENTARY MATERIAL: Supplementary material is available for this article at http://dx.doi.org/10.1007/s00439-006-0160-y and is accessible for authorized users