Approximate Bayesian Computation

Approximate Bayesian computation (ABC) constitutes a class of computational methods rooted in Bayesian statistics. In all model-based statistical inference, the likelihood function is of central importance, since it expresses the probability of the observed data under a particular statistical model, and thus quantifies the support data lend to particular values of parameters and to choices among different models. For simple models, an analytical formula for the likelihood function can typically be derived. However, for more complex models, an analytical formula might be elusive or the likelihood function might be computationally very costly to evaluate. ABC methods bypass the evaluation of the likelihood function. In this way, ABC methods widen the realm of models for which statistical inference can be considered. ABC methods are mathematically well-founded, but they inevitably make assumptions and approximations whose impact needs to be carefully assessed. Furthermore, the wider application domain of ABC exacerbates the challenges of parameter estimation and model selection. ABC has rapidly gained popularity over the last years and in particular for the analysis of complex problems arising in biological sciences (e.g., in population genetics, ecology, epidemiology, and systems biology).

This is a “Topic Page” article for PLOS Computational Biology.

     

Changing Mortality Profile among HIV-Infected Patients in Rio de Janeiro,Brazil: Shifting from AIDS to Non-AIDS Related Conditions in the HAART Era

Introduction

We describe temporal trends in the mortality rates and factors associated with AIDS and non-AIDS related mortality at the Evandro Chagas Clinical Research Institute (IPEC), Oswaldo Cruz Foundation (FIOCRUZ).

Methods

Adult patients enrolling from 1986 through 2009 with a minimum follow up of 60 days were included. Vital status was exhaustively checked using patients’ medical charts, through active contact with individuals and family members and by linkage with the Rio de Janeiro Mortality database using a previously validated algorithm. The CoDe protocol was used to establish the cause of death. Extended Cox proportional hazards models were used for multivariate modeling.

Results

A total of 3530 individuals met the inclusion criteria, out of which 868 (24.6%) deceased; median follow up per patient was 3.9 years (interquartile range 1.7–9.2 years). The dramatic decrease in the overall mortality rates was driven by AIDS-related causes that decreased from 9.19 deaths/100PYs n 1986–1991 to 1.35/100PYs in 2007–2009. Non-AIDS related mortality rates remained stable overtime, at around 1 death/100PYs. Immunodeficiency significantly increased the hazard of both AIDS-related and non-AIDS-related causes of death, while HAART use was strongly associated with a lower hazard of death from either cause.

Conclusions

Our results confirm the remarkable decrease in AIDS-related mortality as the HIV epidemic evolved and alerts to the conditions not traditionally related to HIV/AIDS which are now becoming more frequent, needing careful monitoring.     

Topological Progression in Proliferating Epithelia Is Driven by a Unique Variation in Polygon Distribution

Morphogenesis is consequence of lots of small coordinated variations that occur during development. In proliferating stages, tissue growth is coupled to changes in shape and organization. A number of studies have analyzed the topological properties of proliferating epithelia using the Drosophila wing disc as a model. These works are based in the existence of a fixed distribution of these epithelial cells according to their number of sides. Cell division, cell rearrangements or a combination of both mechanisms have been proposed to be responsible for this polygonal assembling. Here, we have used different system biology methods to compare images from two close proliferative stages that present high morphological similarity. This approach enables us to search for traces of epithelial organization. First, we show that geometrical and network characteristics of individual cells are mainly dependent on their number of sides. Second, we find a significant divergence between the distribution of polygons in epithelia from mid-third instar larva versus early prepupa. We show that this alteration propagates into changes in epithelial organization. Remarkably, only the variation in polygon distribution driven by morphogenesis leads to progression in epithelial organization. In addition, we identify the relevant features that characterize these rearrangements. Our results reveal signs of epithelial homogenization during the growing phase, before the planar cell polarity pathway leads to the hexagonal packing of the epithelium during pupal stages.     

Seascape genetics of the stalked kelp Pterygophora californica and comparative population genetics in the Santa Barbara Channel

We conducted a population genetic analysis of the stalked kelp, Pterygophora californica, in the Santa Barbara Channel, California, USA. The results were compared with previous work on the genetic differentiation of giant kelp, Macrocystis pyrifera, in the same region. These two sympatric kelps not only share many life history and dispersal characteristics but also differ in that dislodged P. californica does not produce floating rafts with buoyant fertile sporophytes, commonly observed for M. pyrifera. We used a comparative population genetic approach with these two species to test the hypothesis that the ability to produce floating rafts increases the genetic connectivity among kelp patches in the Santa Barbara Channel. We quantified the association of habitat continuity and oceanographic distance with the genetic differentiation observed in stalked kelp, like previously conducted for giant kelp. We compared both overall (across all patches) and pairwise (between patches) genetic differentiation. We found that oceanographic transit time, habitat continuity, and geographic distance were all associated with genetic connectivity in P. californica, supporting similar previous findings for M. pyrifera. Controlling for differences in heterozygosity between kelp species using Jost's D_EST, we showed that global differentiation and pairwise differentiation were similar among patches between the two kelp species, indicating that they have similar dispersal capabilities despite their differences in rafting ability. These results suggest that rafting sporophytes do not play a significant role in effective dispersal of M. pyrifera at ecologically relevant spatial and temporal scales.     

Evidence for Interrupted Biomagnification of Cadmium in Billfish Food Chain Based on Stable Carbon and Nitrogen Isotopes from Southwestern of Gulf of California

Biological Trace Element Research - We report cadmium (Cd) concentrations in muscle, liver, and blood of striped marlin (Kajikia audax) and blue marlin (Makaira nigricans), caught in the southwest...     

Characterization of C2C12 cells in simulated microgravity: Possible use for myoblast regeneration

Muscle loss is a major problem for many in lifetime. Muscle and bone degeneration has also been observed in individuals exposed to microgravity and in unloading conditions. C2C12 myoblst cells are able to form myotubes, and myofibers and these cells have been employed for muscle regeneration purposes and in myogenic regeneration and transplantation studies. We exposed C2C12 cells in an random position machine to simulate microgravity and study the energy and the biochemical challenges associated with this treatment. Simulated microgravity exposed C2C12 cells maintain positive proliferation indices and delay the differentiation process for several days. On the other hand this treatment significantly alters many of the biochemical and the metabolic characteristics of the cell cultures including calcium homeostasis. Recent data have shown that these perturbations are due to the inhibition of the ryanodine receptors on the membranes of intracellular calcium stores. We were able to reverse this perturbations treating cells with thapsigargin which prevents the segregation of intracellular calcium ions in the mitochondria and in the sarco/endoplasmic reticula. Calcium homeostasis appear a key target of microgravity exposure. In conclusion, in this study we reported some of the effects induced by the exposure of C2C12 cell cultures to simulated microgravity. The promising information obtained is of fundamental importance in the hope to employ this protocol in the field of regenerative medicine     

Trophic niche partitioning of five skate species of genus Bathyraja in northern and central Patagonia,Argentina

Overexploitation of marine communities can lead to modifications in the structure of the food web and can force organisms like elasmobranchs to change their feeding habits. To evaluate the impact that fisheries have on food webs and on the interactions between species, it is necessary to describe and quantify the diet of the species involved and follow it through time. This study compares the diet of five skate species using the data obtained from the by-catch of the Argentine hake (Merluccius hubbsi) fishery in north and central Patagonia, Argentina. Diet composition was assessed by analysing the digestive tract contents and trophic overlapping between species of the genus Bathyraja: Bathyraja albomaculata, Bathyraja brachyurops, Bathyraja macloviana, Bathyraja magellanica and Bathyraja multispinis. A total of 184 stomachs were analysed. The diets of B. albomaculata and B. macloviana mainly comprised annelids, whereas that of B. brachyurops primarily comprised fish, including hake heads discarded by the fishery. The diets of B. magellanica and B. multispinis were largely based on crustaceans. Despite the morphological similarities and their shared preference for benthic habitats, no complete diet overlaps were found between the different species. These results suggest that these skate species have undergone a process of diet specialisation. This is a common feeding strategy that occurs to successfully eliminate competition when resources are limited, which corresponds to the conditions found in an environment being affected by the pressures of overfishing.