Robust algorithms for solving stochastic partial differential equations

A robust semi-implicit central partial difference algorithm for the numerical solution of coupled stochastic parabolic partial differential equations (PDEs) is described. This can be used for calculating correlation functions of systems of interacting stochastic fields. Such field equations can arise in the description of Hamiltonian and open systems in the physics of nonlinear processes, and may include multiplicative noise sources. The algorithm can be used for studying the properties of nonlinear quantum or classical field theories. The general approach is outlined and applied to a specific example, namely the quantum statistical fluctuations of ultra-short optical pulses in chi((2)) parametric waveguides. This example uses a non-diagonal coherent state representation, and correctly predicts the sub-shot noise level spectral fluctuations observed in homodyne detection measurements. It is expected that the methods used wilt be applicable for higher-order correlation functions and other physical problems as well. A stochastic differencing technique for reducing sampling errors is also introduced. This involves solving nonlinear stochastic parabolic PDEs in combination with a reference process, which uses the Wigner representation in the example presented here. A computer implementation on MIMD parallel architectures is discussed. (C) 1997 Academic Press.

Genetic population structure of the catadromous perciform: Macquaria novemaculeata (Percichthyidae)

Genetic population structure in the catadromous Australian bass Macquaria novemaculeata was investigated using samples from four locations spanning 600 km along the eastern Australian coastline. Both allozymes and mtDNA control region sequences were examined. Population subdivision estimates based on allozymes revealed low levels of population structuring (G(st)=0.043, P<0.05). However, mtDNA indicated moderate levels of geographic population structure (G(st)=0.146, P<0.01). Phylogenetic analysis of mtDNA control region sequences (mean sequence divergence 1.9%) indicated little phylogeographic structuring. Results suggested that genotypic variation within each river population, while bring affected primarily by genetic drift, was also prevented from more significant divergence by homogenizing levels of gene flow-synonymous with a one-dimensional stepping-stone model of population structure. The catadromous life history of Macquaria novemaculeata was considered to br influential on the pattern of population structure displayed. Results were compared to the few population genetic studies involving catadromous fishes, indicating that catadromy alone is unlikely to be a good predictor of population structure. A more comprehensive suite of biological characteristics than simple life-history traits must be considered fully to allow reliable predictive models of population structure to be formulated. (C) 1997 The Fisheries Society of the British Isles.

Strategies for genetic model specification in the screening of genome-wide meta-analysis signals for further replication

Background Meta-analysis is increasingly being employed as a screening procedure in large-scale association studies to select promising variants for follow-up studies. However, standard methods for meta-analysis require the assumption of an underlying genetic model, which is typically unknown a priori. This drawback can introduce model misspecifications, causing power to be suboptimal, or the evaluation of multiple genetic models, which augments the number of false-positive associations, ultimately leading to waste of resources with fruitless replication studies. We used simulated meta-analyses of large genetic association studies to investigate naive strategies of genetic model specification to optimize screenings of genome-wide meta-analysis signals for further replication. Methods Different methods, meta-analytical models and strategies were compared in terms of power and type-I error. Simulations were carried out for a binary trait in a wide range of true genetic models, genome-wide thresholds, minor allele frequencies (MAFs), odds ratios and between-study heterogeneity (tau(2)). Results Among the investigated strategies, a simple Bonferroni-corrected approach that fits both multiplicative and recessive models was found to be optimal in most examined scenarios, reducing the likelihood of false discoveries and enhancing power in scenarios with small MAFs either in the presence or in absence of heterogeneity. Nonetheless, this strategy is sensitive to tau(2) whenever the susceptibility allele is common (MAF epsilon 30%), resulting in an increased number of false-positive associations compared with an analysis that considers only the multiplicative model. Conclusion Invoking a simple Bonferroni adjustment and testing for both multiplicative and recessive models is fast and an optimal strategy in large meta-analysis-based screenings. However, care must be taken when examined variants are common, where specification of a multiplicative model alone may be preferable.

Most of the patients presenting myocardial infarction would not be eligible for intensive lipid-lowering based on clinical algorithms or plasma C-reactive protein

Objective: The study we assessed how often patients who are manifesting a myocardial infarction (MI) would not be considered candidates for intensive lipid-lowering therapy based on the current guidelines. Methods: In 355 consecutive patients manifesting ST elevation MI (STEMI), admission plasma C-reactive protein (CRP) was measured and Framingham risk score (FRS), PROCAM risk score, Reynolds risk score, ASSIGN risk score, QRISK, and SCORE algorithms were applied. Cardiac computed tomography and carotid ultrasound were performed to assess the coronary artery calcium score (CAC), carotid intima-media thickness (cIMT) and the presence of carotid plaques. Results: Less than 50% of STEMI patients would be identified as having high risk before the event by any of these algorithms. With the exception of FRS (9%), all other algorithms would assign low risk to about half of the enrolled patients. Plasma CRP was <1.0 mg/L in 70% and >2 mg/L in 14% of the patients. The average cIMT was 0.8 +/- 0.2 mm and only in 24% of patients was >= 1.0 mm. Carotid plaques were found in 74% of patients. CAC > 100 was found in 66% of patients. Adding CAC >100 plus the presence of carotid plaque, a high-risk condition would be identified in 100% of the patients using any of the above mentioned algorithms. Conclusion: More than half of patients manifesting STEMI would not be considered as candidates for intensive preventive therapy by the current clinical algorithms. The addition of anatomical parameters such as CAC and the presence of carotid plaques can substantially reduce the CVD risk underestimation. (C) 2010 Elsevier Ireland Ltd. All rights reserved.

On the extent of genetic variation for transpiration efficiency in sorghum

A glasshouse study examined 49 diverse sorghum lines for variation in transpiration efficiency. Three of the 49 lines grown were Sorghum spp, native to Australia; one was the major weed Johnson grass (Sorghum halepense), and the remaining 45 lines were cultivars of Sorghum bicolor. All plants were grown under non-limiting water and nutrient conditions using a semi-automatic pot watering system designed to facilitate accurate measurement of water use. Plants were harvested 56-58 days after sowing and dry weights of plant parts were determined. Transpiration efficiency differed significantly among cultivars. The 3 Australian native sorghums had much lower transpiration efficiency than the other 46 cultivars, which ranged from 7.7 to 6.0 g/kg. For the 46 diverse cultivars, the ratio of range in transpiration efficiency to its l.s.d. was 2.0, which was similar to that found among more adapted cultivars in a previous study. This is a significant finding as it suggests that there is likely to be little pay-off from pursuing screening of unadapted material for increased variation in transpiration efficiency. It is necessary, however, also to examine absolute levels of transpiration efficiency to determine whether increased levels have been found. The cultivar with greatest transpiration efficiency in this study (IS9710) had a value 9% greater (P < 0.05) than the accepted standard for adapted sorghum cultivars. The potential impact of such an increase in transpiration efficiency warrants continued effort to capture it. Transpiration efficiency has been related theoretically and experimentally to the degree of carbon isotope discrimination in leaf tissue in sorghum, which thus offers a relatively simple selection index. In this study, the variation in transpiration efficiency was not related simply to carbon isotope discrimination. Significant associations of transpiration efficiency with ash content and indices of photosynthetic capacity were found. However, the associations were not strong. These results suggest that a simple screening technique could not be based on any of the measures or indices analysed in this study. A better understanding of the physiological basis of the observed genetic differences in transpiration efficiency may assist in developing reliable selection indices. It was concluded that the potential value of the improvement in transpiration efficiency over the accepted standard and the degree of genetic variation found warrant further study on this subject. It was suggested that screening for genetic variation under water-limiting conditions may provide useful insights and should be pursued.

Genetic basis of the formation and identity of type I and type II neurons in Drosophila embryos

The embryonic peripheral nervous system of Drosophila contains two main types of sensory neurons: type I neurons, which innervate external sense organs and chordotonal organs, and type II multidendritic neurons, Here, we analyse the origin of the difference between type I and type II in the case of the neurons that depend on the proneural genes of the achaete-scute complex (ASC), We show that, in Notch(-) embryos, the type I neurons are missing while type nr neurons are produced in excess, indicating that the type I/type II choice relies on Notch-mediated cell communication, In contrast, both type I and type II neurons are absent in numb(-) embryos and after ubiquitous expression of tramtrack, indicating that the activity of numb and the absence of tramtrack are required to produce both external sense organ and multidendritic neural fates, The analysis of string(-) embryos reveals that when the precursors are unable to divide they differentiate mostly into type II neurons, indicating that the type II is the default neuronal fate, We also report a new mutant phenotype where the ASC-dependent neurons are converted into-type II neurons, providing evidence for the existence of one or more genes required for maintaining the alternative (type I) fate, Our results suggest that the same mechanism of type I/type II specification may operate at a late step of the ASC-dependent lineages, when multidendritic neurons arise as siblings of the external sense organ neurons and, at an early step, when other multidendritic neurons precursors arise as siblings of external sense organ precursors.

A two-dimensional fuzzy random model of soil pore structure

A new conceptual model for soil pore-solid structure is formalized. Soil pore-solid structure is proposed to comprise spatially abutting elements each with a value which is its membership to the fuzzy set ''pore,'' termed porosity. These values have a range between zero (all solid) and unity (all pore). Images are used to represent structures in which the elements are pixels and the value of each is a porosity. Two-dimensional random fields are generated by allocating each pixel a porosity by independently sampling a statistical distribution. These random fields are reorganized into other pore-solid structural types by selecting parent points which have a specified local region of influence. Pixels of larger or smaller porosity are aggregated about the parent points and within the region of interest by controlled swapping of pixels in the image. This creates local regions of homogeneity within the random field. This is similar to the process known as simulated annealing. The resulting structures are characterized using one-and two-dimensional variograms and functions describing their connectivity. A variety of examples of structures created by the model is presented and compared. Extension to three dimensions presents no theoretical difficulties and is currently under development.

A high-resolution genetic map of the familial Mediterranean fever candidate region allows identification of haplotype-sharing among ethnic groups

Familial Mediterranean fever (FMF) is a recessive disorder of inflammation caused by mutations in a gene (designated MEFV) on chromosome 16p13.3, We have recently constructed a 1-Mb cosmid contig that includes the FMF critical region. Here we show genotype data for 12 markers from our physical map, including 5 newly identified microsatellites, in FMF families. Intrafamilial recombinations placed MEFV in the similar to 285 kb between D16S468/D16S3070 and D16S3376. We observed significant linkage disequilibrium in the North African Jewish population, and historical recombinants in the founder haplotype placed MEFV between D16S3082 and D16S3373 (similar to 200 kb). In smaller panels of Iraqi Jewish, Arab, and Armenian families, there were significant allelic associations only for D16S3370 and D16S2617 among the Armenians. A sizable minority of Iraqi Jewish and Armenian carrier chromosomes appeared to be derived from the North African Jewish ancestral haplotype. We observed a unique FMF haplotype common to Iraqi Jews, Arabs, and Armenians and two other haplotypes restricted to either the Iraqi Jewish or the Armenian population. These data support the view that a few major mutations account for a large percentage of the cases of FMF and suggest that same of these mutations arose before the affected Middle Eastern populations diverged from one another. (C) 1997 Academic Press.

Genetic heterogeneity in familial acute myelogenous leukemia: Evidence for a second locus at chromosome 16q21-23.2

The identification of genes responsible for the rare cases of familial leukemia may afford insight into the mechanism underlying the more common sporadic occurrences. Here we test a single family with 11 relevant meioses transmitting autosomal dominant acute myelogenous leukemia (AML) and myelodysplasia for linkage to three potential candidate loci. In a different family with inherited AML, linkage to chromosome 21q22.1-22.2 was recently reported; we exclude linkage to 21q22.1-22.2, demonstrating that familial AML is a heterogeneous disease. After reviewing familial leukemia and observing anticipation in the form of a declining age of onset with each generation, we had proposed 9p21-22 and 16q22 as additional candidate loci. Whereas linkage to 9p21-22 can be excluded, the finding of a maximum two-point LOD score of 2.82 with the microsatellite marker D16S522 at a recombination fraction theta = 0 provides evidence supporting linkage to 16q22. Haplotype analysis reveals a 23.5-cM (17.9-Mb) commonly inherited region among all affected family members extending from D16S451 to D1GS289, In order to extract maximum linkage information with missing individuals, incomplete informativeness with individual markers in this interval, and possible deviance from strict autosomal dominant inheritance, we performed nonparametric linkage analysis (NPL) and found a maximum NPL statistic corresponding to a P-value of .00098, close to the maximum conditional probability of linkage expected for a pedigree with this structure. Mutational analysis in this region specifically excludes expansion of the AT-rich minisatellite repeat FRA16B fragile site and the CAG trinucleotide repeat in the E2F-4 transcription factor. The ''repeat expansion detection'' method, capable of detecting dynamic mutation associated with anticipation, more generally excludes large CAG repeat expansion as a cause of leukemia in this family.

Molecular epidemiology and genetic diversity of hepatitis B virus genotype E in an isolated Afro-Colombian community

Hepatitis B virus (HBV) infection is a significant public health concern with 350 million chronic carriers worldwide. Eight HBV genotypes (A-H) have been described so far. Genotype E (HBV/E) is widely distributed in West Africa and has rarely been found in other continents, except for a few cases in individuals with an African background. In this study, we characterized HBV genotypes in Quibdo, Colombia, by partial S/P gene sequencing, and found, for the first time, HBV/E circulating in nine Afro-Colombian patients who had no recent contact with Africa. The presence of HBV/E in this community as a monophyletic group suggests that it was a result of a recent introduction by some Afro-descendent contact or, alternatively, that the virus came with slaves brought to Colombia. By using sequences with sampling dates, we estimated the substitution rate to be about 3.2x10(-4) substitutions per site per year, which resulted in a time to the most recent common ancestor (TMRCA) of 29 years. In parallel, we also estimated the TMRCA for HBV/E by using two previously estimated substitution rates (7.7x10(-4) and 1.5x10(-5) substitutions per site per year). The TMRCA was around 35 years under the higher rate and 1500 years under the slower rate. In sum, this work reports for the first time the presence of an exclusively African HBV genotype circulating in South America. We also discuss the time of the entry of this virus into America based on different substitution rates estimated for HBV.

Frequency of concurrent autoimmune disorders in patients with autoimmune hepatitis - Effect of age, gender, and genetic background

Background: Concurrent autoimmune disorders (CAIDs) have been shown to occur in 22% to 34% of the patients with autoimmune hepatitis (AIH). Their presence has been linked to female gender, older age, and to certain HLA antigens, namely HLA-A11. DRB1*04, and DRB4*01. Aims: To assess the frequency and nature of CAID in Brazilian patients with AIH types 1 (AIH-1) and 2 (AIH-2) and to investigate the influence of age, gender, and genetic background in their occurrence. Patients and Methods: The presence and nature of CAID was studied in 143 patients [117 females, median age 11 (1.3 to 69)] with AIH-1 (n = 125) and AIH-2 (n = 28). HLA typing and tumor necrosis factor a gene promoter and exon I cytotoxic T lymphocyte associated antigen 4 (CTLA-4) gene polymorphisms were determined by polymerase chain reaction-based techniques. Results: The frequency of CAID was similar in patients with AIH-1 (14%) and AIH-2 (18%), but their nature was shown to vary. Arthritis was seen in half of the patients (n = 8) with CAID and AIH-1 and in none of those with AIH-2. Subjects with AIH-1 and CAID were shown to be older [24 (1.3 to 6 1) vs. 11 (1.3 to 69) y P = 0.02] and to have more often circulating antinuclear antibody (76% vs. 40%, P = 0.008) and less frequently antiactin antibodies (33% vs. 75%, P = 0.008) when compared with their counterparts without CAID. No particular HLA-DR and DQ alleles, as well as tumor necrosis factor a and CTLA-4 genotypes, were associated with CAID. Conclusions: The nature, but not the frequency, of CAID was shown to vary in AIH-1 and AIH-2. In subjects with AIH-1, CAID was linked to older subjects and to the presence of antinuclear antibody. No predisposition to CAID was associated to HLA-DRB1*04 or DDB4*01 alleles. The observed lower frequency of CAID could be attributed to the lower age of disease onset in Brazilians and to differences in HLA-encoded susceptibility to AIH-1 observed in South America.

Model for Differential Nursing Diagnosis of Alterations in Urinary Elimination Based on Fuzzy Logic

Nursing diagnoses associated with alterations of urinary elimination require different interventions, Nurses, who are not specialists, require support to diagnose and manage patients with disturbances of urine elimination. The aim of this study was to present a model based on fuzzy logic for differential diagnosis of alterations in urinary elimination, considering nursing diagnosis approved by the North American Nursing Diagnosis Association, 2001-2002. Fuzzy relations and the maximum-minimum composition approach were used to develop the system. The model performance was evaluated with 195 cases from the database of a previous study, resulting in 79.0% of total concordance and 19.5% of partial concordance, when compared with the panel of experts. Total discordance was observed in only three cases (1.5%). The agreement between model and experts was excellent (kappa = 0.98, P < .0001) or substantial (kappa = 0.69, P < .0001) when considering the overestimative accordance (accordance was considered when at least one diagnosis was equal) and the underestimative discordance (discordance was considered when at least one diagnosis was different), respectively. The model herein presented showed good performance and a simple theoretical structure, therefore demanding few computational resources.

A Bayesian approach to fuzzy hypotheses testing for the estimation of optimal age for vaccination against measles

Fuzzy Bayesian tests were performed to evaluate whether the mother`s seroprevalence and children`s seroconversion to measles vaccine could be considered as ""high"" or ""low"". The results of the tests were aggregated into a fuzzy rule-based model structure, which would allow an expert to influence the model results. The linguistic model was developed considering four input variables. As the model output, we obtain the recommended age-specific vaccine coverage. The inputs of the fuzzy rules are fuzzy sets and the outputs are constant functions, performing the simplest Takagi-Sugeno-Kang model. This fuzzy approach is compared to a classical one, where the classical Bayes test was performed. Although the fuzzy and classical performances were similar, the fuzzy approach was more detailed and revealed important differences. In addition to taking into account subjective information in the form of fuzzy hypotheses it can be intuitively grasped by the decision maker. Finally, we show that the Bayesian test of fuzzy hypotheses is an interesting approach from the theoretical point of view, in the sense that it combines two complementary areas of investigation, normally seen as competitive. (C) 2007 IMACS. Published by Elsevier B.V. All rights reserved.

A fuzzy Reed-Frost model for epidemic spreading

In this paper, we present a fuzzy approach to the Reed-Frost model for epidemic spreading taking into account uncertainties in the diagnostic of the infection. The heterogeneities in the infected group is based on the clinical signals of the individuals (symptoms, laboratorial exams, medical findings, etc.), which are incorporated into the dynamic of the epidemic. The infectivity level is time-varying and the classification of the individuals is performed through fuzzy relations. Simulations considering a real problem with data of the viral epidemic in a children daycare are performed and the results are compared with a stochastic Reed-Frost generalization.