Complementarities between mathematical and computational modeling: the case of the repeated Prisoners' Dilemma

We study the properties of the well known Replicator Dynamics when applied to a finitely repeated version of the Prisoners' Dilemma game. We characterize the behavior of such dynamics under strongly simplifying assumptions (i.e. only 3 strategies are available) and show that the basin of attraction of defection shrinks as the number of repetitions increases. After discussing the difficulties involved in trying to relax the 'strongly simplifying assumptions' above, we approach the same model by means of simulations based on genetic algorithms. The resulting simulations describe a behavior of the system very close to the one predicted by the replicator dynamics without imposing any of the assumptions of the mathematical model. Our main conclusion is that mathematical and computational models are good complements for research in social sciences. Indeed, while computational models are extremely useful to extend the scope of the analysis to complex scenarios hard to analyze mathematically, formal models can be useful to verify and to explain the outcomes of computational models.

Comparative Meiotic Studies in Triatoma sordida (Stål) and T. guasayana Wygodzinsky & Abalos (Reduviidae, Heteroptera)

Triatoma sordida and T. guasayana are competent Trypanosoma cruzi vectors, with overlapping distribution areas in Argentina. Both species are morphologically similar, and their immature stages are hard to discriminate. Cytogenetic studies in the genus Triatoma reveal scarce karyotypic variations, being 2n= 20 + XY the most frequent diploid number in males. In the present work the meiotic behaviour of different Argentinian populations of T. sordida and T. guasayana has been analyzed; the meiotic karyotype of both species has also been compared. The species differ in total chromosome area and in the relative area of the sex chromosomes. These meiotic karyotypic differences constitute an additional tool for the taxonomic characterization of T. sordida and T. guasayana. The analysis of an interpopulation hybrid of T. sordida (Brazil x Argentina) reveals a regular meiotic behaviour, despite the presence of heteromorphic bivalents. Our observations support the hypothesis that karyotype variations through the gain or loss of heterochromatin can not be considered as a primary mechanism of reproductive isolation in Triatoma.

The stabilizing effects of the acquired immunity on the schistosomiasis transmission modeling - the sensitivity Analysis

A mathematical model is proposed to analyze the effects of acquired immunity on the transmission of schistosomiasis in the human host. From this model the prevalence curve dependent on four parameters can be obtained. These parameters were estimated fitting the data by the maximum likelihood method. The model showed a good retrieving capacity of real data from two endemic areas of schistosomiasis: Touros, Brazil (Schistosoma mansoni) and Misungwi, Tanzania (S. haematobium). Also, the average worm burden per person and the dispersion of parasite per person in the community can be obtained from the model. In this paper, the stabilizing effects of the acquired immunity assumption in the model are assessed in terms of the epidemiological variables as follows. Regarded to the prevalence curve, we calculate the confidence interval, and related to the average worm burden and the worm dispersion in the community, the sensitivity analysis (the range of the variation) of both variables with respect to their parameters is performed.

A comparative study on IgG-ELISA, IgM-IFT and Kato-Katz methods for epidemiological purposes in a low endemic area for schistosomiasis

The high sensitivity and the possibility of automation of the enzyme-linked-immunosorbent-assay (ELISA) has indicated this technique as one of the most useful serological test for epidemiological studies. In the present study, an ELISA for detection of IgG antibodies against adult worm antigens (IgG-ELISA) was investigated for epidemiological purposes, in a rural area of the municipality of Itariri (São Paulo, Brazil). Blood on filter paper (1,180 samples) from about 650 schoolchildren were submitted to ELISA and the data compared to the results of the parasitological method of Kato-Katz and also to the IgM-IFT (immunofluorescence test for IgM antibodies to gut associated antigens). The prevalence rates respectively of 8.5%, 43.0%, and 56.2% by the Kato-Katz, IgG-ELISA, and IgM-IFT methods suggest the poor sensitivity of the parasitological method for detection of Schistosoma mansoni eggs in individuals with low worm burden, situation commonly observed in low endemic areas. These results can partially explain the poor degree of agreement between the IgG-ELISA and the Kato-Katz, as suggested by the Kappa index of 0.170. Otherwise, the Kappa index of 0.675 showed substantial agreement between the two serological tests. Some discrepancy of results between the two serological techniques must be better investigated.

Comparative Analysis by Polymerase Chain Reaction Amplified Minicircles of Kinetoplast DNA of a Stable Strain of Trypanosoma cruzi from São Felipe, Bahia, its Clones and Subclones: Possibility of Predominance of a Principal Clone in this Area

Molecular characterization of one stable strain of Trypanosoma cruzi, the 21 SF, representative of the pattern of strains isolated from the endemic area of São Felipe, State of Bahia, Brazil, maintained for 15 years in laboratory by serial passages in mice and classified as biodeme Type II and zymodeme 2 has been investigated. The kinetoplast DNA (kDNA) of parental strain, 5 clones and 14 subclones were analyzed. Schizodeme was established by comparative study of the fragments obtained from digestion of the 330-bp fragments amplified by polymerase chain reaction (PCR) from the variable regions of the minicicles, and digested by restriction endonucleases Rsa I and Hinf I. Our results show a high percentual of similarity between the restriction fragment lenght polymorphism (RFLP) for the parental strain and its clones and among these individual clones and their subclones at a level of 80 to 100%.This homology indicates a predominance of the same "principal clone" in the 21SF strain and confirms the homogeneity previously observed at biological and isozymic analysis. These results suggest the possibility that the T. cruzi strains with similar biological and isoenzymic patterns, circulating in this endemic area, are representative of one dominant clone. The presence of "principal clones" could be responsible for a predominant tropism of the parasites for specific organs and tissues and this could contribute to the pattern of clinico-pathological manifestations of Chagas's disease in one geographical area.

Comparative angiotensin II receptor blockade in healthy volunteers: the importance of dosing.

OBJECTIVES: We have reported previously that 80 mg valsartan and 50 mg losartan provide less receptor blockade than 150 mg irbesartan in normotensive subjects. In this study we investigated the importance of drug dosing in mediating these differences by comparing the AT(1)-receptor blockade induced by 3 doses of valsartan with that obtained with 3 other antagonists at given doses. METHODS: Valsartan (80, 160, and 320 mg), 50 mg losartan, 150 mg irbesartan, and 8 mg candesartan were administered to 24 healthy subjects in a randomized, open-label, 3-period crossover study. All doses were given once daily for 8 days. The angiotensin II receptor blockade was assessed with two techniques, the reactive rise in plasma renin activity and an in vitro radioreceptor binding assay that quantified the displacement of angiotensin II by the blocking agents. Measurements were obtained before and 4 and 24 hours after drug intake on days 1 and 8. RESULTS: At 4 and 24 hours, valsartan induced a dose-dependent "blockade" of AT(1) receptors. Compared with other antagonists, 80 mg valsartan and 50 mg losartan had a comparable profile. The 160-mg and 320-mg doses of valsartan blocked AT(1) receptors at 4 hours by 80%, which was similar to the effect of 150 mg irbesartan. At trough, however, the valsartan-induced blockade was slightly less than that obtained with irbesartan. With use of plasma renin activity as a marker of receptor blockade, on day 8, 160 mg valsartan was equivalent to 150 mg irbesartan and 8 mg candesartan. CONCLUSIONS: These results show that the differences in angiotensin II receptor blockade observed with the various AT(1) antagonists are explained mainly by differences in dosing. When 160-mg or 320-mg doses were investigated, the effects of valsartan hardly differed from those obtained with recommended doses of irbesartan and candesartan.

Policy Credibility and Delegation to Independent Regulatory Agencies: A Comparative Empirical Analysis

Independent regulatory agencies are one of the main institutional features of the 'rising regulatory state' in Western Europe. Governments are increasingly willing to abandon their regulatory competencies and to delegate them to specialized institutions that are at least partially beyond their control. This article examines the empirical consistency of one particular explanation of this phenomenon, namely the credibility hypothesis, claiming that governments delegate powers so as to enhance the credibility of their policies. Three observable implications are derived from the general hypothesis, linking credibility and delegation to veto players, complexity and interdependence. An independence index is developed to measure agency independence, which is then used in a multivariate analysis where the impact of credibility concerns on delegation is tested. The analysis relies on an original data set comprising independence scores for thirty-three regulators. Results show that the credibility hypothesis can explain a good deal of the variation in delegation. The economic nature of regulation is a strong determinant of agency independence, but is mediated by national institutions in the form of veto players.

Synchronous versus asynchronous modeling of gene regulatory networks.

MOTIVATION: In silico modeling of gene regulatory networks has gained some momentum recently due to increased interest in analyzing the dynamics of biological systems. This has been further facilitated by the increasing availability of experimental data on gene-gene, protein-protein and gene-protein interactions. The two dynamical properties that are often experimentally testable are perturbations and stable steady states. Although a lot of work has been done on the identification of steady states, not much work has been reported on in silico modeling of cellular differentiation processes. RESULTS: In this manuscript, we provide algorithms based on reduced ordered binary decision diagrams (ROBDDs) for Boolean modeling of gene regulatory networks. Algorithms for synchronous and asynchronous transition models have been proposed and their corresponding computational properties have been analyzed. These algorithms allow users to compute cyclic attractors of large networks that are currently not feasible using existing software. Hereby we provide a framework to analyze the effect of multiple gene perturbation protocols, and their effect on cell differentiation processes. These algorithms were validated on the T-helper model showing the correct steady state identification and Th1-Th2 cellular differentiation process. AVAILABILITY: The software binaries for Windows and Linux platforms can be downloaded from http://si2.epfl.ch/~garg/genysis.html.

Repair of non-circumferential cervical trachea defects by three different latissimus dorsi flaps. A comparative studyin an experimental setting.

BACKGROUND: Large intrathoracic airway defects may be closed using a pedicled latissimus dorsi (LD) flap, with rewarding results. This study addresses the question of whether this holds true for extrathoracic non-circumferential tracheal defects. METHODS: A cervical segment of the trachea of 4 x 1 cm was resected in 9 white male pigs. The defect was stented with a silicone stent for 3 months and closed either by an LD flap alone (group a, n = 3), an LD flap with an attached rib segment covered by pleura (group b, n = 3), or an LD flap reinforced by a perforated polylactide (MacroPore) plate (group c, n = 3). The trachea was assessed by rigid endoscopy at 3 and 4 months and histologically at 4 months postoperatively. RESULTS: The degree of stenosis at the level of the reconstruction at 4 months was 25, 50 and 75% in group a, 15, 50 and 60% in group b, and 20, 95 and 95% in group c, respectively. The percentage of the defect covered by columnar epithelium was 100% in all animals of group a, 60, 100 and 100% in group b, and 10, 0 and 0% in group c. Resorption of the rib was seen in all animals of group b and obstructive inflammatory polyps were found in 2 animals of group c. CONCLUSION: Pedicled LD flaps provided less satisfactory results for closure of large non-circumferential extrathoracic airway defects than observed after intrathoracic reconstruction. A pedicled rib segment added to the LD flap did not improve the results obtained from LD flap repair alone, and an embedded MacroPore prosthesis may result in severe airway stenosis due to plate migration and intense inflammatory reaction protruding into the tracheal lumen.

Comparative genetic analyses point to HCP5 as susceptibility locus for HCV-associated hepatocellular carcinoma.

BACKGROUND & AIMS: Recently, genetic variations in MICA (lead single nucleotide polymorphism [SNP] rs2596542) were identified by a genome-wide association study (GWAS) to be associated with hepatitis C virus (HCV)-related hepatocellular carcinoma (HCC) in Japanese patients. In the present study, we sought to determine whether this SNP is predictive of HCC development in the Caucasian population as well. METHODS: An extended region around rs2596542 was genotyped in 1924 HCV-infected patients from the Swiss Hepatitis C Cohort Study (SCCS). Pair-wise correlation between key SNPs was calculated both in the Japanese and European populations (HapMap3: CEU and JPT). RESULTS: To our surprise, the minor allele A of rs2596542 in proximity of MICA appeared to have a protective impact on HCC development in Caucasians, which represents an inverse association as compared to the one observed in the Japanese population. Detailed fine-mapping analyses revealed a new SNP in HCP5 (rs2244546) upstream of MICA as strong predictor of HCV-related HCC in the SCCS (univariable p=0.027; multivariable p=0.0002, odds ratio=3.96, 95% confidence interval=1.90-8.27). This newly identified SNP had a similarly directed effect on HCC in both Caucasian and Japanese populations, suggesting that rs2244546 may better tag a putative true variant than the originally identified SNPs. CONCLUSIONS: Our data confirms the MICA/HCP5 region as susceptibility locus for HCV-related HCC and identifies rs2244546 in HCP5 as a novel tagging SNP. In addition, our data exemplify the need for conducting meta-analyses of cohorts of different ethnicities in order to fine map GWAS signals.

Preliminary regional rockfall hazard mapping using LIDAR-based slope frequency distribution and CONEFALL modeling

A factor limiting preliminary rockfall hazard mapping at regional scale is often the lack of knowledge of potential source areas. Nowadays, high resolution topographic data (LiDAR) can account for realistic landscape details even at large scale. With such fine-scale morphological variability, quantitative geomorphometric analyses become a relevant approach for delineating potential rockfall instabilities. Using digital elevation model (DEM)-based ?slope families? concept over areas of similar lithology and cliffs and screes zones available from the 1:25,000 topographic map, a susceptibility rockfall hazard map was drawn up in the canton of Vaud, Switzerland, in order to provide a relevant hazard overview. Slope surfaces over morphometrically-defined thresholds angles were considered as rockfall source zones. 3D modelling (CONEFALL) was then applied on each of the estimated source zones in order to assess the maximum runout length. Comparison with known events and other rockfall hazard assessments are in good agreement, showing that it is possible to assess rockfall activities over large areas from DEM-based parameters and topographical elements.

Functional analysis and structural modeling of human APOBEC3G reveal the role of evolutionarily conserved elements in the inhibition of human immunodeficiency virus type 1 infection and Alu transposition.

Retroelements are important evolutionary forces but can be deleterious if left uncontrolled. Members of the human APOBEC3 family of cytidine deaminases can inhibit a wide range of endogenous, as well as exogenous, retroelements. These enzymes are structurally organized in one or two domains comprising a zinc-coordinating motif. APOBEC3G contains two such domains, only the C terminal of which is endowed with editing activity, while its N-terminal counterpart binds RNA, promotes homo-oligomerization, and is necessary for packaging into human immunodeficiency virus type 1 (HIV-1) virions. Here, we performed a large-scale mutagenesis-based analysis of the APOBEC3G N terminus, testing mutants for (i) inhibition of vif-defective HIV-1 infection and Alu retrotransposition, (ii) RNA binding, and (iii) oligomerization. Furthermore, in the absence of structural information on this domain, we used homology modeling to examine the positions of functionally important residues and of residues found to be under positive selection by phylogenetic analyses of primate APOBEC3G genes. Our results reveal the importance of a predicted RNA binding dimerization interface both for packaging into HIV-1 virions and inhibition of both HIV-1 infection and Alu transposition. We further found that the HIV-1-blocking activity of APOBEC3G N-terminal mutants defective for packaging can be almost entirely rescued if their virion incorporation is forced by fusion with Vpr, indicating that the corresponding region of APOBEC3G plays little role in other aspects of its action against this pathogen. Interestingly, residues forming the APOBEC3G dimer interface are highly conserved, contrasting with the rapid evolution of two neighboring surface-exposed amino acid patches, one targeted by the Vif protein of primate lentiviruses and the other of yet-undefined function.

A comparative study on the different staining methods and number of specimens for the detection of acid fast bacilli

The presence of acid fast bacilli in multiple specimens was investigated comparatively with Ziehl-Neelsen (ZN) and fluorescence microscopy (FM) staining in order to determine sensitivity in detecting tuberculosis (TB). A total of 465 specimens obtained from 295 patients were analysed at Harran University Medical School Hospital between March 1998 and March 2000. The culture was employed as the reference method. Sixty-eight patients (23.1%) were diagnosed as having TB by culture. The ZN and FM staining sensitivities were 67.6% (46/68) and 85.2% (58/68) respectively. Two hundred and one patients (68.1%) submitted one specimen to the laboratory. TB positivity was detected in 42 (20.9%) of these patients by culture. The sensitivities of ZN and FM stains were found to be 61% and 83% in these patients. However, in 18 patients (6.1%) who submitted two specimens to the laboratory, the TB was positive in six of them (33.3%) and ZN and FM sensitivities were 66% and 83% respectively. When three specimens or more were collected from the patients (76 patients, 25.8%), TB positivity was determined in 20 of them (26.3%) and the sensitivities were 80% and 92% in the ZN- and FM-stained smears, respectively. Our data indicate that in the diagnosis of TB, FM has greater sensitivity than ZN. In particular, in the case of a single specimen, the diagnostic value of FM is quite significant. It is, therefore, possible to conclude that both ZN and FM staining can be used for the diagnosis of TB when there are more than two specimens. However, if only one or two specimens are available, FM staining is preferable.