Computing and visually analyzing mutual information in molecular co-evolution

<p>Abstract</p> <p>Background</p> <p>Selective pressure in molecular evolution leads to uneven distributions of amino acids and nucleotides. In fact one observes correlations among such constituents due to a large number of biophysical mechanisms (folding properties, electrostatics, ...). To quantify these correlations the mutual information -after proper normalization - has proven most effective. The challenge is to navigate the large amount of data, which in a study for a typical protein cannot simply be plotted.</p> <p>Results</p> <p>To visually analyze mutual information we developed a matrix visualization tool that allows different views on the mutual information matrix: filtering, sorting, and weighting are among them. The user can interactively navigate a huge matrix in real-time and search e.g., for patterns and unusual high or low values. A computation of the mutual information matrix for a sequence alignment in FASTA-format is possible. The respective stand-alone program computes in addition proper normalizations for a null model of neutral evolution and maps the mutual information to <it>Z</it>-scores with respect to the null model.</p> <p>Conclusions</p> <p>The new tool allows to compute and visually analyze sequence data for possible co-evolutionary signals. The tool has already been successfully employed in evolutionary studies on HIV1 protease and acetylcholinesterase. The functionality of the tool was defined by users using the tool in real-world research. The software can also be used for visual analysis of other matrix-like data, such as information obtained by DNA microarray experiments. The package is platform-independently implemented in <monospace>Java</monospace> and free for academic use under a GPL license.</p

Dispelling urban myths about default uncertainty factors in chemical risk assessment - Sufficient protection against mixture effects?

© 2013 Martin et al.; licensee BioMed Central LtdThis article has been made available through the Brunel Open Access Publishing Fund.Assessing the detrimental health effects of chemicals requires the extrapolation of experimental data in animals to human populations. This is achieved by applying a default uncertainty factor of 100 to doses not found to be associated with observable effects in laboratory animals. It is commonly assumed that the toxicokinetic and toxicodynamic sub-components of this default uncertainty factor represent worst-case scenarios and that the multiplication of those components yields conservative estimates of safe levels for humans. It is sometimes claimed that this conservatism also offers adequate protection from mixture effects. By analysing the evolution of uncertainty factors from a historical perspective, we expose that the default factor and its sub-components are intended to represent adequate rather than worst-case scenarios. The intention of using assessment factors for mixture effects was abandoned thirty years ago. It is also often ignored that the conservatism (or otherwise) of uncertainty factors can only be considered in relation to a defined level of protection. A protection equivalent to an effect magnitude of 0.001-0.0001% over background incidence is generally considered acceptable. However, it is impossible to say whether this level of protection is in fact realised with the tolerable doses that are derived by employing uncertainty factors. Accordingly, it is difficult to assess whether uncertainty factors overestimate or underestimate the sensitivity differences in human populations. It is also often not appreciated that the outcome of probabilistic approaches to the multiplication of sub-factors is dependent on the choice of probability distributions. Therefore, the idea that default uncertainty factors are overly conservative worst-case scenarios which can account both for the lack of statistical power in animal experiments and protect against potential mixture effects is ill-founded. We contend that precautionary regulation should provide an incentive to generate better data and recommend adopting a pragmatic, but scientifically better founded approach to mixture risk assessment. © 2013 Martin et al.; licensee BioMed Central Ltd.Oak Foundatio

Mathematical modelling of lymphatic filariasis elimination programmes in India: Required duration of mass drug administration and post-treatment level of infection indicators

Background: India has made great progress towards the elimination of lymphatic filariasis. By 2015, most endemic districts had completed at least five annual rounds of mass drug administration (MDA). The next challenge is to determine when MDA can be stopped. We performed a simulation study with the individual-based model LYMFASIM to help clarify this. Methods: We used a model-variant for Indian settings. We considered different hypotheses on detectability of antigenaemia (Ag) in relation to underlying adult worm burden, choosing the most likely hypothesis by comparing the model predicted association between community-level microfilaraemia (Mf) and antigenaemia (Ag) prevalence levels to observed data (collated from literature). Next, we estimated how long MDA must be continued in order to achieve elimination in different transmission settings and what Mf and Ag prevalence may still remain 1 year after the last required MDA round. The robustness of key-outcomes was assessed in a sensitivity analysis. Results: Our model matched observed data qualitatively well when we assumed an Ag detection rate of 50 % for single worm infections, which increases with the number of adult worms (modelled by relating detection to the presence of female worms). The required duration of annual MDA increased with higher baseline endemicity and lower coverage (varying between 2 and 12 rounds), while the remaining residual infection 1 year after the last required treatment declined with transmission intensity. For low and high transmission settings, the median residual infection levels were 1.0 % and 0.4 % (Mf prevalence in the 5+ population), and 3.5 % and 2.0 % (Ag prevalence in 6-7 year-old children). Conclusion: To achieve elimination in high transmission settings, MDA must be continued longer and infection levels must be reduced to lower levels than in low-endemic communities. Although our simulations were for Indian settings, qualitatively similar patterns are also expected in other areas. This should be taken into account in decision algorithms to define whether MDA can be interrupted. Transmission assessment surveys should ideally be targeted to communities with the highest pre-control transmission levels, to minimize the risk of programme failure

Applicability of a short/rapid 13C-urea breath test for Helicobacter pylori: retrospective multicenter chart review study

<p>Abstract</p> <p>Background</p> <p>Carbon labeled urea breath tests usually entail a two point sampling with a 20 to 30-minute gap. Our aim was to evaluate the duration of time needed for diagnosing <it>Helicobacter pylori </it>by the BreathID^®System.</p> <p>Methods</p> <p>This is a retrospective multicenter chart review study. Test location, date, delta over baseline, and duration of the entire test were recorded. Consecutively ¹³C urea breath tests results were extracted from the files over a nine year period.</p> <p>Results</p> <p>Of the 12,791 tests results, 35.1% were positively diagnosed and only 0.1% were inconclusive. A statistically significant difference in prevalence among the countries was found: Germany showing the lowest, 13.3%, and Israel the highest, 44.1%. Significant differences were found in time to diagnosis: a positive diagnosis had the shortest and an inconclusive result had the longest. Overall test duration averaged 15.1 minutes in Germany versus approximately 13 minutes in other countries. Diagnosis was achieved after approximately 9 minutes in Israel, Italy and Switzerland, but after 10 on average in the others. The mean delta over baseline value for a negative diagnosis was 1.03 ± 0.86, (range, 0.9 - 5), versus 20.2 ± 18.9, (range, 5.1 - 159.4) for a positive one.</p> <p>Conclusions</p> <p>The BreathID^®System used in diagnosing <it>Helicobacter pylori </it>can safely shorten test duration on average of 10-13 minutes without any loss of sensitivity or specificity and with no test lasting more than 21 minutes.</p

Discovery of nucleotide polymorphisms in the Musa gene pool by Ecotilling

Musa (banana and plantain) is an important genus for the global export market and in local markets where it provides staple food for approximately 400 million people. Hybridization and polyploidization of several (sub)species, combined with vegetative propagation and human selection have produced a complex genetic history. We describe the application of the Ecotilling method for the discovery and characterization of nucleotide polymorphisms in diploid and polyploid accessions of Musa. We discovered over 800 novel alleles in 80 accessions. Sequencing and band evaluation shows Ecotilling to be a robust and accurate platform for the discovery of polymorphisms in homologous and homeologous gene targets. In the process of validating the method, we identified two single nucleotide polymorphisms that may be deleterious for the function of a gene putatively important for phototropism. Evaluation of heterozygous polymorphism and haplotype blocks revealed a high level of nucleotide diversity in Musa accessions. We further applied a strategy for the simultaneous discovery of heterozygous and homozygous polymorphisms in diploid accessions to rapidly evaluate nucleotide diversity in accessions of the same genome type. This strategy can be used to develop hypotheses for inheritance patterns of nucleotide polymorphisms within and between genome types. We conclude that Ecotilling is suitable for diversity studies in Musa, that it can be considered for functional genomics studies and as tool in selecting germplasm for traditional and mutation breeding approaches

Histopathology of Growth Anomaly Affecting the Coral, Montipora capitata: Implications on Biological Functions and Population Viability

Growth anomalies (GAs) affect the coral, Montipora capitata, at Wai'ōpae, southeast Hawai'i Island. Our histopathological analysis of this disease revealed that the GA tissue undergoes changes which compromise anatomical machinery for biological functions such as defense, feeding, digestion, and reproduction. GA tissue exhibited significant reductions in density of ova (66.1–93.7%), symbiotic dinoflagellates (38.8–67.5%), mesenterial filaments (11.2–29.0%), and nematocytes (28.8–46.0%). Hyperplasia of the basal body wall but no abnormal levels of necrosis and algal or fungal invasion was found in GA tissue. Skeletal density along the basal body wall was significantly reduced in GAs compared to healthy or unaffected sections. The reductions in density of the above histological features in GA tissue were collated with disease severity data to quantify the impact of this disease at the colony and population level. Resulting calculations showed this disease reduces the fecundity of M. capitata colonies at Wai'ōpae by 0.7–49.6%, depending on GA severity, and the overall population fecundity by 2.41±0.29%. In sum, GA in this M. capitata population reduces the coral's critical biological functions and increases susceptibility to erosion, clearly defining itself as a disease and an ecological threat

Transient Receptor Potential Ion Channels Control Thermoregulatory Behaviour in Reptiles

Biological functions are governed by thermodynamics, and animals regulate their body temperature to optimise cellular performance and to avoid harmful extremes. The capacity to sense environmental and internal temperatures is a prerequisite for the evolution of thermoregulation. However, the mechanisms that enable ectothermic vertebrates to sense heat remain unknown. The recently discovered thermal characteristics of transient receptor potential ion channels (TRP) render these proteins suitable to act as temperature sensors. Here we test the hypothesis that TRPs are present in reptiles and function to control thermoregulatory behaviour. We show that the hot-sensing TRPV1 is expressed in a crocodile (Crocodylus porosus), an agamid (Amphibolurus muricatus) and a scincid (Pseudemoia entrecasteauxii) lizard, as well as in the quail and zebrafinch (Coturnix chinensis and Poephila guttata). The TRPV1 genes from all reptiles form a unique clade that is delineated from the mammalian and the ancestral Xenopus sequences by an insertion of two amino acids. TRPV1 and the cool-sensing TRPM8 are expressed in liver, muscle (transversospinalis complex), and heart tissues of the crocodile, and have the potential to act as internal thermometer and as external temperatures sensors. Inhibition of TRPV1 and TRPM8 in C. porosus abolishes the typically reptilian shuttling behaviour between cooling and heating environments, and leads to significantly altered body temperature patterns. Our results provide the proximate mechanism of thermal selection in terrestrial ectotherms, which heralds a fundamental change in interpretation, because TRPs provide the mechanism for a tissue-specific input into the animals' thermoregulatory response

Ricin Toxicokinetics and Its Sensitive Detection in Mouse Sera or Feces Using Immuno-PCR

Ricin (also called RCA-II or RCA(60)), one of the most potent toxins and documented bioweapons, is derived from castor beans of Ricinus communis. Several in vitro methods have been designed for ricin detection in complex food matrices in the event of intentional contamination. Recently, a novel Immuno-PCR (IPCR) assay was developed with a limit of detection of 10 fg/ml in a buffer matrix and about 10-1000-fold greater sensitivity than other methods in various food matrices.In order to devise a better diagnostic test for ricin, the IPCR assay was adapted for the detection of ricin in biological samples collected from mice after intoxication. The limit of detection in both mouse sera and feces was as low as 1 pg/ml. Using the mouse intravenous (iv) model for ricin intoxication, a biphasic half-life of ricin, with a rapid t(1/2)α of 4 min and a slower t(1/2)β of 86 min were observed. The molecular biodistribution time for ricin following oral ingestion was estimated using an antibody neutralization assay. Ricin was detected in the blood stream starting at approximately 6-7 h post- oral intoxication. Whole animal histopathological analysis was performed on mice treated orally or systemically with ricin. Severe lesions were observed in the pancreas, spleen and intestinal mesenteric lymph nodes, but no severe pathology in other major organs was observed.The determination of in vivo toxicokinetics and pathological effects of ricin following systemic and oral intoxication provide a better understanding of the etiology of intoxication and will help in the future design of more effective diagnostic and therapeutic methods

Transposon Excision from an Atypical Site: A Mechanism of Evolution of Novel Transposable Elements

The role of transposable elements in sculpting the genome is well appreciated but remains poorly understood. Some organisms, such as humans, do not have active transposons; however, transposable elements were presumably active in their ancestral genomes. Of specific interest is whether the DNA surrounding the sites of transposon excision become recombinogenic, thus bringing about homologous recombination. Previous studies in maize and Drosophila have provided conflicting evidence on whether transposon excision is correlated with homologous recombination. Here we take advantage of an atypical Dissociation (Ds) element, a maize transposon that can be mobilized by the Ac transposase gene in Arabidopsis thaliana, to address questions on the mechanism of Ds excision. This atypical Ds element contains an adjacent 598 base pairs (bp) inverted repeat; the element was allowed to excise by the introduction of an unlinked Ac transposase source through mating. Footprints at the excision site suggest a micro-homology mediated non-homologous end joining reminiscent of V(D)J recombination involving the formation of intra-helix 3′ to 5′ trans-esterification as an intermediate, a mechanism consistent with previous observations in maize, Antirrhinum and in certain insects. The proposed mechanism suggests that the broken chromosome at the excision site should not allow recombinational interaction with the homologous chromosome, and that the linked inverted repeat should also be mobilizable. To test the first prediction, we measured recombination of flanking chromosomal arms selected for the excision of Ds. In congruence with the model, Ds excision did not influence crossover recombination. Furthermore, evidence for correlated movement of the adjacent inverted repeat sequence is presented; its origin and movement suggest a novel mechanism for the evolution of repeated elements. Taken together these results suggest that the movement of transposable elements themselves may not directly influence linkage. Possibility remains, however, for novel repeated DNA sequences produced as a consequence of transposon movement to influence crossover in subsequent generations