Measuring gametic disequilibrium from multilocus data

We describe a Bayesian approach to analyzing multilocus genotype or haplotype data to assess departures from gametic (linkage) equilibrium. Our approach employs a Markov chain Monte Carlo (MCMC) algorithm to approximate the posterior probability distributions of disequilibrium parameters. The distributions are computed exactly in some simple settings. Among other advantages, posterior distributions can be presented visually, which allows the uncertainties in parameter estimates to be readily assessed. In addition, background knowledge can be incorporated, where available, to improve the precision of inferences. The method is illustrated by application to previously published datasets; implications for multilocus forensic match probabilities and for simple association-based gene mapping are also discussed.

Model predictive control using dynamic integrated system optimisation and parameter estimation (DISOPE)

DISOPE is a technique for solving optimal control problems where there are differences in structure and parameter values between reality and the model employed in the computations. The model reality differences can also allow for deliberate simplification of model characteristics and performance indices in order to facilitate the solution of the optimal control problem. The technique was developed originally in continuous time and later extended to discrete time. The main property of the procedure is that by iterating on appropriately modified model based problems the correct optimal solution is achieved in spite of the model-reality differences. Algorithms have been developed in both continuous and discrete time for a general nonlinear optimal control problem with terminal weighting, bounded controls and terminal constraints. The aim of this paper is to show how the DISOPE technique can aid receding horizon optimal control computation in nonlinear model predictive control.

Hierarchical analysis of production efficiency in a coastal trawl fishery

We present, pedagogically, the Bayesian approach to composed error models under alternative, hierarchical characterizations; demonstrate, briefly, the Bayesian approach to model comparison using recent advances in Markov Chain Monte Carlo (MCMC) methods; and illustrate, empirically, the value of these techniques to natural resource economics and coastal fisheries management, in particular. The Bayesian approach to fisheries efficiency analysis is interesting for at least three reasons. First, it is a robust and highly flexible alternative to commonly applied, frequentist procedures, which dominate the literature. Second,the Bayesian approach is extremely simple to implement, requiring only a modest addition to most natural-resource economist tool-kits. Third, despite its attractions, applications of Bayesian methodology in coastal fisheries management are few.

Bayesian ranking and slection of fishing boat efficiencies

The steadily accumulating literature on technical efficiency in fisheries attests to the importance of efficiency as an indicator of fleet condition and as an object of management concern. In this paper, we extend previous work by presenting a Bayesian hierarchical approach that yields both efficiency estimates and, as a byproduct of the estimation algorithm, probabilistic rankings of the relative technical efficiencies of fishing boats. The estimation algorithm is based on recent advances in Markov Chain Monte Carlo (MCMC) methods— Gibbs sampling, in particular—which have not been widely used in fisheries economics. We apply the method to a sample of 10,865 boat trips in the US Pacific hake (or whiting) fishery during 1987–2003. We uncover systematic differences between efficiency rankings based on sample mean efficiency estimates and those that exploit the full posterior distributions of boat efficiencies to estimate the probability that a given boat has the highest true mean efficiency.

Isolation and characterization of an AINTEGUMENTA-like gene in different coconut (Cocos nucifera L.) varieties from Sri Lanka

Development of an efficient tissue culture protocol in coconut is hampered by numerous technical constraints. Thus a greater understanding of the fundamental aspects of embryogenesis is essential. The role of AINTEGUMENTA-like genes in embryogenesis has been elucidated not only in model plants but also in economically important crops. A coconut gene, CnANT, that encodes two APETALA2 (AP2) domains and a conserved linker region similar to those of the BABY BOOM transcription factor was cloned, characterized, and its tissue specific expression was examined. The full-length cDNA of 1,780 bp contains a 1,425-bp open reading frame that encodes a putative peptide of 474 amino acids. The genomic DNA sequence includes 2,317 bp and consists of nine exons interrupted by eight introns. The exon/intron organization of CnANT is similar to that of homologous genes in other plant species. Analysis of differential tissue expression by real-time polymerase chain reaction indicated that CnANT is expressed more highly in in vitro grown tissues than in other vegetative tissues. Sequence comparison of the genomic sequence of CnANT in different coconut varieties revealed one single nucleotide polymorphism and one indel in the first exon and first intron, respectively, which differentiate the Tall group of trees from Dwarfs. The indel sequence, which can be considered a simple sequence repeats marker, was successfully used to distinguish the Tall and Dwarf groups as well as to develop a marker system, which may be of value in the identification of parental varieties that are used in coconut breeding programs in Sri Lanka.

Microcontroller-based peak current mode control using digital slope compensation

Microcontroller-based peak current mode control of a buck converter is investigated. The new solution uses a discrete time controller with digital slope compensation. This is implemented using only a single-chip microcontroller to achieve desirable cycle-by-cycle peak current limiting. The digital controller is implemented as a two-pole, two-zero linear difference equation designed using a continuous time model of the buck converter and a discrete time transform. Subharmonic oscillations are removed with digital slope compensation using a discrete staircase ramp. A 16 W hardware implementation directly compares analog and digital control. Frequency response measurements are taken and it is shown that the crossover frequency and expected phase margin of the digital control system match that of its analog counterpart.

Environmentally-acquired bacteria influence microbial diversity and natural innate immune responses at gut surfaces

Background: Early microbial colonization of the gut reduces the incidence of infectious, inflammatory and autoimmune diseases. Recent population studies reveal that childhood hygiene is a significant risk factor for development of inflammatory bowel disease, thereby reinforcing the hygiene hypothesis and the potential importance of microbial colonization during early life. The extent to which early-life environment impacts on microbial diversity of the adult gut and subsequent immune processes has not been comprehensively investigated thus far. We addressed this important question using the pig as a model to evaluate the impact of early-life environment on microbe/host gut interactions during development. Results: Genetically-related piglets were housed in either indoor or outdoor environments or in experimental isolators. Analysis of over 3,000 16S rRNA sequences revealed major differences in mucosa-adherent microbial diversity in the ileum of adult pigs attributable to differences in earlylife environment. Pigs housed in a natural outdoor environment showed a dominance of Firmicutes, in particular Lactobacillus, whereas animals housed in a hygienic indoor environment had reduced Lactobacillus and higher numbers of potentially pathogenic phylotypes. Our analysis revealed a strong negative correlation between the abundance of Firmicutes and pathogenic bacterial populations in the gut. These differences were exaggerated in animals housed in experimental isolators. Affymetrix microarray technology and Real-time Polymerase Chain Reaction revealed significant gut-specific gene responses also related to early-life environment. Significantly, indoorhoused pigs displayed increased expression of Type 1 interferon genes, Major Histocompatibility Complex class I and several chemokines. Gene Ontology and pathway analysis further confirmed these results.

Linguistic evidence supports date for Homeric epics

The Homeric epics are among the greatest masterpieces of literature, but when they were produced is not known with certainty. Here we apply evolutionary-linguistic phylogenetic statistical methods to differences in Homeric, Modern Greek and ancient Hittite vocabulary items to estimate a date of approximately 710–760 BCE for these great works. Our analysis compared a common set of vocabulary items among the three pairs of languages, recording for each item whether the words in the two languages were cognate – derived from a shared ancestral word – or not. We then used a likelihood-based Markov chain Monte Carlo procedure to estimate the most probable times in years separating these languages given the percentage of words they shared, combined with knowledge of the rates at which different words change. Our date for the epics is in close agreement with historians' and classicists' beliefs derived from historical and archaeological sources.

Pulsating equilibria: stability through migration

This paper is to present a model of spatial equilibrium using a nonlinear generalization of Markov-chain type model, and to show the dynamic stability of a unique equilibrium. Even at an equilibrium, people continue to migrate among regions as well as among agent-types, and yet their overall distribution remain unchanged. The model is also adapted to suggest a theory of traffic distribution in a city.

A generalised Bayesian instrumental variable approach under student t-distributed errors with application

Bayesian analysis is given of an instrumental variable model that allows for heteroscedasticity in both the structural equation and the instrument equation. Specifically, the approach for dealing with heteroscedastic errors in Geweke (1993) is extended to the Bayesian instrumental variable estimator outlined in Rossi et al. (2005). Heteroscedasticity is treated by modelling the variance for each error using a hierarchical prior that is Gamma distributed. The computation is carried out by using a Markov chain Monte Carlo sampling algorithm with an augmented draw for the heteroscedastic case. An example using real data illustrates the approach and shows that ignoring heteroscedasticity in the instrument equation when it exists may lead to biased estimates.

How large is congressional dependence in agriculture. Bayesian inference about ‘scale’ and ‘scope’ in measuring a spatial externality

The political economy literature on agriculture emphasizes influence over political outcomes via lobbying conduits in general, political action committee contributions in particular and the pervasive view that political preferences with respect to agricultural issues are inherently geographic. In this context, ‘interdependence’ in Congressional vote behaviour manifests itself in two dimensions. One dimension is the intensity by which neighboring vote propensities influence one another and the second is the geographic extent of voter influence. We estimate these facets of dependence using data on a Congressional vote on the 2001 Farm Bill using routine Markov chain Monte Carlo procedures and Bayesian model averaging, in particular. In so doing, we develop a novel procedure to examine both the reliability and the consequences of different model representations for measuring both the ‘scale’ and the ‘scope’ of spatial (geographic) co-relations in voting behaviour.

Effects of phthalic acid esters on the liver and thyroid

The effects, over periods from 3 days to 9 months of administration, of diets containing di-2-ethylhexyl phthalate are very similar to those observed in rats administered diets containing hypolipidemic drugs such as clofibrate. Changes occur in a characteristic order commencing with alterations in the distribution of lipid within the liver, quickly followed by proliferation of hepatic peroxisomes and induction of the specialized P-450 isoenzyme(s) catalyzing omega oxidation of fatty acids. There follows a phase of mild liver damage indicated by induction of glucose-6-phosphatase activity and a loss of glycogen, eventually leading to the formation of enlarged lysosomes through autophagy and the accumulation of lipofuscin. Associated changes are found in the kidney and thyroid. The renal changes are limited to the proximal convoluted tubules and are generally similar to changes found in the liver. The effects on the thyroid are more marked. Although the levels of thyroxine in plasma fail to about half normal values, serum triiodothyronine remains close to normal values while the appearance of the thyroid varies, very marked hyperactivity being noted 7 days after commencement of treatment, this is less marked at 14 days, but even after 9 months treatment there is clear cut evidence for hyperactivity with colloid changes which indicate this has persisted for some time. Straight chain analogs of di-2-ethylhexyl phthalate, di-n-hexyl phthalate and di-n-oxtyl phthalate differ entirely in their short-term effects on the liver and kidney but have similar effects on the thyroid. The short-term in vivo hepatic effects of the three phthalate esters can be reproduced in hepatocytes in tissue culture. All three phthalate esters, as well as clofibrate, have early marked effects on the metabolism of fatty acids in isolated hepatocytes. The nature of these changes is such as to increase storage of lipid in the liver. A hypothesis is presented to explain the progress from these initial metabolic effects to the final formation of liver tumors.

Smart-phone based electrocardiogram wavelet decomposition and neural network classification

This paper discusses ECG classification after parametrizing the ECG waveforms in the wavelet domain. The aim of the work is to develop an accurate classification algorithm that can be used to diagnose cardiac beat abnormalities detected using a mobile platform such as smart-phones. Continuous time recurrent neural network classifiers are considered for this task. Records from the European ST-T Database are decomposed in the wavelet domain using discrete wavelet transform (DWT) filter banks and the resulting DWT coefficients are filtered and used as inputs for training the neural network classifier. Advantages of the proposed methodology are the reduced memory requirement for the signals which is of relevance to mobile applications as well as an improvement in the ability of the neural network in its generalization ability due to the more parsimonious representation of the signal to its inputs.

The endophytic fungus Piriformospora indica protects wheat from Fusarium crown rot disease in simulated UK autumn conditions

The root endophytic fungus Piriformospora indica (Sebacinacea) forms mutualistic symbioses with a broad range of host plants, increasing their biomass production and resistance to fungal pathogens. We evaluated the effect of P. indica on Fusarium crown rot disease of wheat, under in vitro and glasshouse conditions. Interaction of P. indica and Fusarium isolates under axenic culture conditions indicated no direct antagonistic activity of P. indica against Fusarium isolates. Seedlings of wheat were inoculated with P. indica and pathogenic Fusarium culmorum or F. graminearum and grown in sterilised soil-free medium or in a non-sterilised mix of soil and sand. Fusarium alone reduced emergence and led to visible browning and reduced root growth. Roots of seedlings in pots inoculated with both Fusarium isolates and P. indica were free of visible symptoms; seed emergence and root biomass were equivalent to the uninoculated. DNA was quantified by real-time polymerase chain reaction (qPCR). The ratio of Fusarium DNA to wheat DNA rose rapidly in the plants inoculated with Fusarium alone; isolates and species were not significantly different. P. indica inoculation reduced the ratio of Fusarium to host DNA in the root systems. The reduction increased with time. The ratio of P. indica to wheat DNA initially rose but then declined in root systems without Fusarium. With Fusarium, the ratio rose throughout the experiment. The absolute amount of Fusarium DNA in root systems increased in the absence of P. indica but was static in plants co-inoculated with P. indica.

Expression of Circadian Clock and Melatonin Receptors within Cultured Rat Cardiomyocytes

Melatonin, the pineal gland hormone, provides entrainment of many circadian rhythms to the ambient light/dark cycle. Recently, cardiovascular studies have demostrated melatonin interactions with many physiological processes and diseases, such as hypertension and cardiopathologies. Although membrane melatonin receptors (MT1, MT2) and the transcriptional factor ROR alpha have been reported to be expressed in the heart, there is no evidence of the cell-type expressing receptors as well as the possible role of melatonin on the expression of the circadian clock of cardiomyocytes, which play an important role in cardiac metabolism and function. Therefore, the aim of this study was to evaluate the mRNA and protein expressions of MT1, MT2, and ROR alpha and to determine whether melatonin directly influences expression of circadian clocks within cultured rat cardiomyocytes. Adult rat cardiomyocyte cultures were created, and the cells were stimulated with 1 nM melatonin or vehicle. Gene expressions were assayed by real-time polymerase chain reaction (PCR). The mRNA and protein expressions of membrane melatonin receptors and RORa were established within adult rat cardiomyocytes. Two hours of melatonin stimulation did not alter the expression pattern of the analyzed genes. However, given at the proper time, melatonin kept Rev-erb alpha expression chronically high, specifically 12 h after melatonin treatment, avoiding the rhythmic decline of Rev-erb alpha mRNA. The blockage of MT1 and MT2 by luzindole did not alter the observed melatonin-induced expression of Rev-erb alpha mRNA, suggesting the nonparticipation of MT1 and MT2 on the melatonin effect within cardiomyocytes. It is possible to speculate that melatonin, in adult rat cardiomyocytes, may play an important role in the light signal transduction to peripheral organs, such as the heart, modulating its intrinsic rhythmicity. (Author correspondence: cipolla@icb.usp.br)