Mathematical Models of Tumor-Immune System Dynamics

"This collection of papers offers a broad synopsis of state-of-the-art mathematical methods used in modeling the interaction between tumors and the immune system. These papers were presented at the four-day workshop on Mathematical Models of Tumor-Immune System Dynamics held in Sydney, Australia from January 7th to January 10th, 2013. The workshop brought together applied mathematicians, biologists, and clinicians actively working in the field of cancer immunology to share their current research and to increase awareness of the innovative mathematical tools that are applicable to the growing field of cancer immunology. Recent progress in cancer immunology and advances in immunotherapy suggest that the immune system plays a fundamental role in host defense against tumors and could be utilized to prevent or cure cancer. Although theoretical and experimental studies of tumor-immune system dynamics have a long history, there are still many unanswered questions about the mechanisms that govern the interaction between the immune system and a growing tumor. The multidimensional nature of these complex interactions requires a cross-disciplinary approach to capture more realistic dynamics of the essential biology. The papers presented in this volume explore these issues and the results will be of interest to graduate students and researchers in a variety of fields within mathematical and biological sciences."--Publisher website

Longitudinal Assessment of Neuropathy in Diabetes using novel ophthalmic MARKers (LANDMark): Baseline findings

Purpose:Over the past decade, corneal nerve morphology and corneal sensation threshold have been explored as potential surrogate markers for the evaluation of diabetic neuropathy. We present the baseline findings of a Longitudinal Assessment of Neuropathy in Diabetes using novel ophthalmic Markers (LANDMark). Methods:The LANDMark Study is a 5-year, two-site, natural history (observational) study of individuals with Type 1 diabetes stratified into those with (T1W) and without (T1WO) neuropathy according to the Toronto criteria, and control subjects. All study participants undergo detailed annual assessment of neuropathy including corneal nerve parameters measured using corneal confocal microscopy and corneal sensitivity measured using non-contact corneal esthesiometry. Results:396 eligible individuals (208 in Brisbane and 188 in Manchester) were assessed: 76 T1W, 166 T1WO and 154 controls. Corneal sensation threshold (mbars) was significantly higher in T1W (1.0 ± 1.1) than T1WO (0.7 ± 0.7) and controls (0.6 ± 0.4) (P=0.002); post-hoc analysis (PHA) revealed no difference between T1WO and controls (Tukey HSD, P=0.502). Corneal nerve fiber length (mm/mm2) (CNFL) was lower in T1W (13.8 ± 6.4) than T1WO (19.1 ± 5.8) and controls (23.2 ± 6.3) (P<0.001); PHA revealed CNFL to be lower in T1W than T1WO, and lower in both of these groups than controls (P<0.001). Corneal nerve branch density (branches/mm2) (CNBD) was significantly lower in T1W (40 ± 32) than T1WO (62 ± 37) and controls (83 ± 46) (P<0.001); PHA showed CNBD was lower in T1W than T1WO, and lower in both groups than controls (P<0.001). Alcohol and cigarette consumption did not differ between groups, although age, BMI, BP, waist circumference, HbA1c, albumin-creatinine ratio, and cholesterol were slightly greater in T1W than T1WO (p<0.05). Some site differences were observed. Conclusions:The LANDMark baseline findings confirm that corneal sensitivity and corneal nerve morphometry can detect differences in neuropathy status in individuals with Type 1 diabetes and healthy controls. Corneal nerve morphology is significantly abnormal even in diabetic patients ‘without neuropathy’ compared to control participants. Results of the longitudinal trial will assess the capability of these tests for monitoring change in these parameters over time as potential surrogate markers for neuropathy.

Exploring Retinal Markers of Diabetic Neuropathy

Purpose : To investigate the application of retinal nerve fibre layer (RNFL) thickness as a marker for severity of diabetic peripheral neuropathy (DPN) in people with Type 2 diabetes. Methods : This was a cross-sectional study whereby 61 participants (mean age 61 [41-75 years], mean duration of diabetes 14 [1-40 years], 70% male) with Type 2 diabetes and DPN underwent optical coherence tomography (OCT) scans. Global and 4 quadrant (TSNI) RNFL thicknesses were measured at 3.45mm around the optic nerve head of one eye. Neuropathy disability score (NDS) was used to assess the severity of DPN on a 0 to 10 scale. Participants were divided into three age-matched groups representing mild (NDS=3-5), moderate (NDS=6-8) and severe (NDS=9-10) neuropathy. Two regression models were fitted for statistical analysis: 1) NDS scores as co-variate for global and quadrant RNFL thicknesses, 2) NDS groups as a factor for global RNFL thickness only. Results : Mean (SD) RNFL thickness (µm) was 103(9) for mild neuropathy (n=34), 101(10) for moderate neuropathy (n=16) and 95(13) in the group with severe neuropathy (n=11). Global RNFL thickness and NDS scores were statistically significantly related (b=-1.20, p=0.048). When neuropathy was assessed across groups, a trend of thinner mean RNFL thickness was observed with increasing severity of neuropathy; however, this result was not statistically significant (F=2.86, p=0.065). TSNI quadrant analysis showed that mean RNFL thickness reduction in the inferior quadrant was 2.55 µm per 1 unit increase in NDS score (p=0.005). However, the regression coefficients were not statistically significant for RNFL thickness in the superior (b=-1.0, p=0.271), temporal (b=-0.90, p=0.238) and nasal (b=-0.99, p=0.205) quadrants. Conclusions : RNFL thickness was reduced with increasing severity of DPN and the effect was most evident in the inferior quadrant. Measuring RNFL thickness using OCT may prove to be a useful, non-invasive technique for identifying severity of DPN and may also provide additional insight into common mechanisms for peripheral neuropathy and RNFL damage.

Comparative analysis of the Metarhizium anisopliae secretome in response to exposure to the greyback cane grub and grub cuticles

Metarhizium anisopliae is a well-characterized biocontrol agent of a wide range of insects including cane grubs. In this study, a two-dimensional (2D) electrophoresis was used to display secreted proteins of M. anisopliae strain FI-1045 growing on the whole greyback cane grubs and their isolated cuticles. Hydrolytic enzymes secreted by M. anisopliae play a key role in insect cuticle-degradation and initiation of the infection process. We have identified all the 101 protein spots displayed by cross-species identification (CSI) from the fungal kingdom. Among the identified proteins were 64-kDa serine carboxypeptidase, 1,3 beta-exoglucanase, Dynamin GTPase, THZ kinase, calcineurin like phosphoesterase, and phosphatidylinositol kinase secreted by M. ansiopliae (FI-1045) in response to exposure to the greyback cane grubs and their isolated cuticles. These proteins have not been previously identified from the culture supernatant of M. anisopliae during infection. To our knowledge, this the first proteomic map established to study the extracellular proteins secreted by M. ansiopliae (FI-1045) during infection of greyback cane grubs and its cuticles.

Adaptations to a subterranean environment and longevity revealed by the analysis of mole rat genomes

Subterranean mammals spend their lives in dark, unventilated environments that are rich in carbon dioxide and ammonia and low in oxygen. Many of these animals are also long-lived and exhibit reduced aging-associated diseases, such as neurodegenerative disorders and cancer. We sequenced the genome of the Damaraland mole rat (DMR, Fukomys damarensis) and improved the genome assembly of the naked mole rat (NMR, Heterocephalus glaber). Comparative genome analyses, along with the transcriptomes of related subterranean rodents, revealed candidate molecular adaptations for subterranean life and longevity, including a divergent insulin peptide, expression of oxygen-carrying globins in the brain, prevention of high CO2-induced pain perception, and enhanced ammonia detoxification. Juxtaposition of the genomes of DMR and other more conventional animals with the genome of NMR revealed several truly exceptional NMR features: unusual thermogenesis, an aberrant melatonin system, pain insensitivity, and unique processing of 28S rRNA. Together, these genomes and transcriptomes extend our understanding of subterranean adaptations, stress resistance, and longevity.

A comparative analysis of relevant crop carbon footprint calculators, with reference to cotton production in Australia

An increasing concern over the sustainability credentials of food and fiber crops require that farmers and their supply chain partners have access to appropriate and industry-friendly tools to be able to measure and improve the outcomes. This article focuses on one of the sustainability indicators, namely, greenhouse gas (GHG) emissions, and nine internationally accredited carbon footprint calculators were identified and compared on an outcomes basis against the same cropping data from a case study cotton farm. The purpose of this article is to identify the most “appropriate” methodology to be applied by cotton suppliers in this regard. From the analysis of the results, we subsequently propose a new integrated model as the basis for an internationally accredited carbon footprint tool for cotton and show how the model can be applied to evaluate the emission outcomes of different farming practices.

Innate transcriptional networks activated in bladder in response to uropathogenic Escherichia coli drive diverse biological pathways and rapid synthesis of IL-10 for defense against bacterial urinary tract infection

Early transcriptional activation events that occur in bladder immediately following bacterial urinary tract infection (UTI) are not well defined. In this study, we describe the whole bladder transcriptome of uropathogenic Escherichia coli (UPEC) cystitis in mice using genome-wide expression profiling to define the transcriptome of innate immune activation stemming from UPEC colonization of the bladder. Bladder RNA from female C57BL/6 mice, analyzed using 1.0 ST-Affymetrix microarrays, revealed extensive activation of diverse sets of innate immune response genes, including those that encode multiple IL-family members, receptors, metabolic regulators, MAPK activators, and lymphocyte signaling molecules. These were among 1564 genes differentially regulated at 2 h postinfection, highlighting a rapid and broad innate immune response to bladder colonization. Integrative systems-level analyses using InnateDB (http://www.innatedb.com) bioinformatics and ingenuity pathway analysis identified multiple distinct biological pathways in the bladder transcriptome with extensive involvement of lymphocyte signaling, cell cycle alterations, cytoskeletal, and metabolic changes. A key regulator of IL activity identified in the transcriptome was IL-10, which was analyzed functionally to reveal marked exacerbation of cystitis in IL-10–deficient mice. Studies of clinical UTI revealed significantly elevated urinary IL-10 in patients with UPEC cystitis, indicating a role for IL-10 in the innate response to human UTI. The whole bladder transcriptome presented in this work provides new insight into the diversity of innate factors that determine UTI on a genome-wide scale and will be valuable for further data mining. Identification of protective roles for other elements in the transcriptome will provide critical new insight into the complex cascade of events that underpin UTI.

Association of CYP1B1 codon 432 mutant allele in head and neck squamous cell cancer is reflected by somatic mutations of p53 in tumor tissue

Tobacco use is causally associated with head and neck squamous cell cancer (HNSCC). Here, we present the results of a case-control study that investigated the effects that the genetic variants of the cytochrome (CYP)1A1, CYP1B1, glutathione-S-transferase (GST)M1, GSTT1, and GSTP1 genes have on modifying the risk of smoking-related HNSCC. Allelisms of the CYP1A1, GSTT1, GSTM1, and GSTT1 genes alone were not associated with an increased risk. CYP1B1 codon 432 polymorphism was found to be a putative susceptibility factor in smoking-related HNSCC. The frequency of CYP1B1 polymorphism was significantly higher (P < 0.001) in the group of smoking cases when compared with smoking controls. Additionally, an odds ratio (OR) of 4.53 (2.62-7.98) was discovered when investigating smoking and nonsmoking cases for the susceptible genotype CYP1B1*2/*2, when compared with the presence of the genotype wild type. In combination with polymorphic variants of the GST genes, a synergistic-effect OR was observed. The calculated OR for the combined genotype CYP1B1*2/*2 and GSTM1*2/*2 was 12.8 (4.09-49.7). The calculated OR for the combined genotype was 13.4 (2.92-97.7) for CYP1B1*2/*2 and GSTT1*2/*2, and 24.1 (9.36-70.5) for the combination of CYP1B1*2/*2 and GSTT1-expressors. The impact of the polymorphic variants of the CYP1B1 gene on HNSCC risk is reflected by the strong association with the frequency of somatic mutations of the p53 gene. Smokers with susceptible genotype CYP1B1*2/*2 were 20 times more likely to show evidence of p53 mutations than were those with CYP1B1 wild type. Combined genotype analysis of CYP1B1 and GSTM1 or GSTT1 revealed interactive effects on the occurrence of p53 gene mutations. The results of the present study indicate that polymorphic variants of CYP1B1 relate significantly to the individual susceptibility of smokers to HNSCC.

Cytochrome p450 1b1, a new keystone in gene-environment interactions related to human head and neck cancer?

Alcohol consumption and tobacco smoking are major causes of head and neck cancers, and regional differences point to the importance of research into gene-environment interactions. Much interest has been focused on polymorphisms of CYP1A1 and of GSTM1 and GSTT1, but a number of studies have not demonstrated significant effects. This has mostly been ascribed to small sample sizes. In general, the impact of polymorphisms of metabolic enzymes appears inconsistent, with some reports of weak-to-moderate associations, and with others of no elevation of risks. The classical cytochrome P450 isoenzyme considered for metabolic activation of polycyclic aromatic hydrocarbons (PAH) is CYP1A1. A new member of the CYP1 family, CYP1B1, was cloned in 1994, currently representing the only member of the CYP1B subfamily. A number of single nucleotide polymorphisms of the CYP1B1 gene have been reported. The amino acid substitutions Val432Leu (CYP1B1*3) and Asn453Ser (CYP1B1*4), located in the heme binding domain of CYP1B1, appear as likely candidates to be linked with biological effects. CYP1B1 activates a wide range of PAH, aromatic and heterocyclic amines. Very recently, the CYP1B1 codon 432 polymorphism (CYP1B1*3) has been identified as a susceptibility factor in smoking-related head-and-neck squamous cell cancer. The impact of this polymorphic variant of CYP1B1 on cancer risk was also reflected by an association with the frequency of somatic mutations of the p53 gene. Combined genotype analysis of CYP1B1 and the glutathione transferases GSTM1 or GSTT1 has pointed to interactive effects. This provides new molecular evidence that tobacco smoke-specific compounds relevant to head and neck carcinogenesis are metabolically activated through CYP1B1 and is consistent with a major pathogenetic relevance of PAH as ingredients of tobacco smoke.

Markers of genetic susceptibility in human environmental hygiene and toxicology : the role of selected CYP, NAT and GST genes

Inherited genetic traits co-determine the susceptibility of an individual to a toxic chemical. Special emphasis has been put on individual responses to environmental and industrial carcinogens, but other chronic diseases are of increasing interest. Polymorphisms of relevant xenobiotic metabolising enzymes may be used as toxicological susceptibility markers. A growing number of genes encoding enzymes involved in biotransformation of toxicants and in cellular defence against toxicant-induced damage to the cells has been identified and cloned, leading to increased knowledge of allelic variants of genes and genetic defects that may result in a differential susceptibility toward environmental toxicants. "Low penetrating" polymorphisms in metabolism genes tend to be much more common in the population than allelic variants of "high penetrating" cancer genes, and are therefore of considerable importance from a public health point of view. Positive associations between cancer and CYP1A1 alleles, in particular the *2C I462V allele, were found for tissues following the aerodigestive tract. Again, in most cases, the effect of the variant CYP1A1 allele becomes apparent or clearer in connection with the GSTM1 null allele. The CYP1B1 codon 432 polymorphism (CYP1B1*3) has been identified as a susceptibility factor in smoking-related head-and-neck squameous cell cancer. The impact of this polymorphic variant of CYP1B1 on cancer risk was also reflected by an association with the frequency of somatic mutations of the p53 gene. Combined genotype analysis of CYP1B1 and the glutathione transferases GSTM1 or GSTT1 has also pointed to interactive effects. Of particular interest for the industrial and environmental field is the isozyme CYP2E1. Several genotypes of this isozyme have been characterised which seem to be associated with different levels of expression of enzyme activity. The acetylator status for NAT2 can be determined by genotyping or by phenotyping. In the pathogenesis of human bladder cancer due to occupational exposure to "classical" aromatic amines (benzidine, 4-aminodiphenyl, 1-naphthylamine) acetylation by NAT2 is regarded as a detoxication step. Interestingly, the underlying European findings of a higher susceptibility of slow acetylators towards aromatic amines are in contrast to findings in Chinese workers occupationally exposed to aromatic amines which points to different mechanisms of susceptibility between European and Chinese populations. Regarding human bladder cancer, the hypothesis has been put forward that genetic polymorphism of GSTM1 might be linked with the occurrence of this tumour type. This supports the hypothesis that exposure to PAH might causally be involved in urothelial cancers. The human polymorphic GST catalysing conjugation of halomethanes, dihalomethanes, ethylene oxide and a number of other industrial compounds could be characterised as a class theta enzyme (GSTT1) by means of molecular biology. "Conjugator" and "non-conjugator" phenotypes are coincident with the presence and absence of the GSTT1 gene. There are wide variations in the frequencies of GSTT1 deletion (GSTT1 *0/0) among different ethnicities. Human phenotyping is facilitated by the GST activity towards methyl bromide or ethylene oxide in erythrocytes which is representative of the metabolic GSTT1 competence of the entire organism. Inter-individual variations in xenobiotic metabolism capacities may be due to polymorphisms of the genes coding for the enzymes themselves or of the genes coding for the receptors or transcription factors which regulate the expression of the enzymes. Also, polymorphisms in several regions of genes may cause altered ligand affinity, transactivation activity or expression levels of the receptor subsequently influencing the expression of the downstream target genes. Studies of individual susceptibility to toxicants and gene-environment interaction are now emerging as an important component of molecular epidemiology.

Biological monitoring in workers in a nitrobenzene reduction plant : Haemoglobin versus serum albumin adducts

The high priority of monitoring workers exposed to nitrobenzene is a consequence of clear findings of experimental carcinogenicity of nitrobenzene and the associated evaluations by the International Agency for Research on Cancer. Eighty male employees of a nitrobenzene reduction plant, with potential skin contact with nitrobenzene and aniline, participated in a current medical surveillance programme. Blood samples were routinely taken and analysed for aniline, 4-aminodiphenyl (4-ADP) and benzidine adducts of haemoglobin (Hb) and human serum albumin (HSA). Also, levels of methaemoglobin (Met-Hb) and of carbon monoxide haemoglobin (CO-Hb) were monitored. Effects of smoking were straightforward. Using the rank sum test of Wilcoxon, we found that very clear-cut and statistically significant smoking effects (about 3-fold increases) were apparent on CO-Hb (P = 0.00085) and on the Hb adduct of 4-ADP (P = 0.0006). The mean aniline-Hb adduct level in smokers was 1.5 times higher than in non-smokers; the significance (P = 0.05375) was close to the 5% level. The strongest correlation was evident between the Hb and HSA adducts of aniline (rs = 0.846). Less pronounced correlations (but with P values < 0.02) appeared between aniline-Hb and 4-ADP-Hb adducts (rs = 0.388), between 4-ADP and 4-ADP-HSA adducts (rs = 0.373), and between 4-ADP-Hb and aniline-HSA adducts (rs = 0.275). In view of the proposal for additional use of the aniline-HSA adduct for biological monitoring, particularly in cases of acute overexposures or poisonings, the strong correlation of the Hb and HSA conjugates is noteworthy; the ratio aniline-HSA:aniline-Hb was 1:42 for the entire cohort.

DNA methylation at the novel CpG sites in the promoter of MED15/PCQAP gene as a biomarker for head and neck cancers

Head and neck cancers (HNCs) represent a significant and ever-growing burden to the modern society, mainly due to the lack of early diagnostic methods. A significant number of HNCs is often associated with drinking, smoking, chewing beetle nut, and human papilloma virus (HPV) infections. We have analyzed DNA methylation patterns in tumor and normal tissue samples collected from head and neck squamous cell carcinoma (HNSCC) patients who were smokers. We have identified novel methylation sites in the promoter of the mediator complex subunit 15 (MED15/PCQAP) gene (encoing a co-factor important for regulation of transcription initiation for promoters of many genes), hypermethylated specifically in tumor cells. Two clusters of CpG dinucleotides methylated in tumors, but not in normal tissue from the same patients, were identified. These CpG methylation events in saliva samples were further validated in a separate cohort of HNSCC patients (who developed cancer due to smoking or HPV infections) and healthy controls using methylation-specific PCR (MSP). We used saliva as a biological medium because of its non-invasive nature, close proximity to the tumors, easiness and it is an economically viable option for large-scale screening studies. The methylation levels for the two identified CpG clusters were significantly different between the saliva samples collected from healthy controls and HNSCC individuals (Welch's t-test returning P, 0.05 and Mann-Whitney test P, 0.01 for both). The developed MSP assays also provided a good discriminative ability with AUC values of 0.70 (P, 0.01) and 0.63 (P, 0.05). The identified novel CpG methylation sites may serve as potential non-invasive biomarkers for detecting HNSCC. © the authors.

Morphometric analysis of telencephalic structure in a variety of neognath and paleognath bird species reveals regional differences associated with specific behavioral traits

Birds exhibit a huge array of behavior, ecology and physiology, and occupy nearly every environment on earth, ranging from the desert outback of Australia to the tropical rain forests of Panama. Some birds have adopted a fully nocturnal lifestyle, such as the barn owl and kiwi, while others, such as the albatross, spend nearly their entire life flying over the ocean. Each species has evolved unique adaptations over millions of years to function in their respective niche. In order to increase processing power or network efficiency, many of these adaptations require enlargements and/or specializations of the brain as a whole or of specific brain regions. In this study, we examine the relative size and morphology of 9 telencephalic regions in a number of Paleognath and Neognath birds and relate the findings to differences in behavior and sensory ecology. We pay particular attention to those species that have undergone a relative enlargement of the telencephalon to determine whether this relative increase in telencephalic size is homogeneous across different brain regions or whether particular regions have become differentially enlarged. The analysis indicates that changes in the relative size of telencephalic regions are not homogeneous, with every species showing hypertrophy or hypotrophy of at least one of them. The three-dimensional structure of these regions in different species was also variable, in particular that of the mesopallium in kiwi. The findings from this study provide further evidence that the changes in relative brain size in birds reflect a process of mosaic evolution.

Acoustic identification of 12 species of echolocating bat by discriminant function analysis and artificial neural networks

We recorded echolocation calls from 14 sympatric species of bat in Britain. Once digitised, one temporal and four spectral features were measured from each call. The frequency-time course of each call was approximated by fitting eight mathematical functions, and the goodness of fit, represented by the mean-squared error, was calculated. Measurements were taken using an automated process that extracted a single call from background noise and measured all variables without intervention. Two species of Rhinolophus were easily identified from call duration and spectral measurements. For the remaining 12 species, discriminant function analysis and multilayer back-propagation perceptrons were used to classify calls to species level. Analyses were carried out with and without the inclusion of curve-fitting data to evaluate its usefulness in distinguishing among species. Discriminant function analysis achieved an overall correct classification rate of 79% with curve-fitting data included, while an artificial neural network achieved 87%. The removal of curve-fitting data improved the performance of the discriminant function analysis by 2 %, while the performance of a perceptron decreased by 2 %. However, an increase in correct identification rates when curve-fitting information was included was not found for all species. The use of a hierarchical classification system, whereby calls were first classified to genus level and then to species level, had little effect on correct classification rates by discriminant function analysis but did improve rates achieved by perceptrons. This is the first published study to use artificial neural networks to classify the echolocation calls of bats to species level. Our findings are discussed in terms of recent advances in recording and analysis technologies, and are related to factors causing convergence and divergence of echolocation call design in bats.