Learning genetic epistasis using Bayesian network scoring criteria

<p>Abstract</p> <p>Background</p> <p>Gene-gene epistatic interactions likely play an important role in the genetic basis of many common diseases. Recently, machine-learning and data mining methods have been developed for learning epistatic relationships from data. A well-known combinatorial method that has been successfully applied for detecting epistasis is <it>Multifactor Dimensionality Reduction </it>(MDR). Jiang et al. created a combinatorial epistasis learning method called <it>BNMBL </it>to learn Bayesian network (BN) epistatic models. They compared BNMBL to MDR using simulated data sets. Each of these data sets was generated from a model that associates two SNPs with a disease and includes 18 unrelated SNPs. For each data set, BNMBL and MDR were used to score all 2-SNP models, and BNMBL learned significantly more correct models. In real data sets, we ordinarily do not know the number of SNPs that influence phenotype. BNMBL may not perform as well if we also scored models containing more than two SNPs. Furthermore, a number of other BN scoring criteria have been developed. They may detect epistatic interactions even better than BNMBL.</p> <p>Although BNs are a promising tool for learning epistatic relationships from data, we cannot confidently use them in this domain until we determine which scoring criteria work best or even well when we try learning the correct model without knowledge of the number of SNPs in that model.</p> <p>Results</p> <p>We evaluated the performance of 22 BN scoring criteria using 28,000 simulated data sets and a real Alzheimer's GWAS data set. Our results were surprising in that the Bayesian scoring criterion with large values of a hyperparameter called α performed best. This score performed better than other BN scoring criteria and MDR at <it>recall </it>using simulated data sets, at detecting the hardest-to-detect models using simulated data sets, and at substantiating previous results using the real Alzheimer's data set.</p> <p>Conclusions</p> <p>We conclude that representing epistatic interactions using BN models and scoring them using a BN scoring criterion holds promise for identifying epistatic genetic variants in data. In particular, the Bayesian scoring criterion with large values of a hyperparameter α appears more promising than a number of alternatives.</p

Observations on the agonistic behaviour of Lepadogaster lepadogaster purpurea (Pisces: Gobiesocidae)

Agonistic behaviour of captive Lepadogaster lepadogaster purpurea was described, forming the first record of such behaviour in the family Gobiesocidae. Both male and female engaged in similar behaviours. The agonistic interactions were related mainly to the control of shelter sites and larger fish tended to chase away and take over the shelter of smaller fish

The Cluster Variation Method for Efficient Linkage Analysis on Extended Pedigrees

BACKGROUND: Computing exact multipoint LOD scores for extended pedigrees rapidly becomes infeasible as the number of markers and untyped individuals increase. When markers are excluded from the computation, significant power may be lost. Therefore accurate approximate methods which take into account all markers are desirable. METHODS: We present a novel method for efficient estimation of LOD scores on extended pedigrees. Our approach is based on the Cluster Variation Method, which deterministically estimates likelihoods by performing exact computations on tractable subsets of variables (clusters) of a Bayesian network. First a distribution over inheritances on the marker loci is approximated with the Cluster Variation Method. Then this distribution is used to estimate the LOD score for each location of the trait locus. RESULTS: First we demonstrate that significant power may be lost if markers are ignored in the multi-point analysis. On a set of pedigrees where exact computation is possible we compare the estimates of the LOD scores obtained with our method to the exact LOD scores. Secondly, we compare our method to a state of the art MCMC sampler. When both methods are given equal computation time, our method is more efficient. Finally, we show that CVM scales to large problem instances. CONCLUSION: We conclude that the Cluster Variation Method is as accurate as MCMC and generally is more efficient. Our method is a promising alternative to approaches based on MCMC sampling

Night Shift: Expansion of Temporal Niche Use Following Reductions in Predator Density

Predation shapes many fundamental aspects of ecology. Uncertainty remains, however, about whether predators can influence patterns of temporal niche construction at ecologically relevant timescales. Partitioning of time is an important mechanism by which prey avoid interactions with predators. However, the traits that control a prey organism's capacity to operate during a particular portion of the diel cycle are diverse and complex. Thus, diel prey niches are often assumed to be relatively unlikely to respond to changes in predation risk at short timescales. Here we present evidence to the contrary. We report results that suggest that the anthropogenic depletion of daytime active predators (species that are either diurnal or cathemeral) in a coral reef ecosystem is associated with rapid temporal niche expansions in a multi-species assemblage of nocturnal prey fishes. Diurnal comparisons of nocturnal prey fish abundance in predator rich and predator depleted reefs at two atolls revealed that nocturnal fish were approximately six (biomass) and eight (density) times more common during the day on predator depleted reefs. Amongst these, the prey species that likely were the most specialized for nocturnal living, and thus the most vulnerable to predation (i.e. those with greatest eye size to body length ratio), showed the strongest diurnal increases at sites where daytime active predators were rare. While we were unable to determine whether these observed increases in diurnal abundance by nocturnal prey were the result of a numerical or behavioral response, either effect could be ecologically significant. These results raise the possibility that predation may play an important role in regulating the partitioning of time by prey and that anthropogenic depletions of predators may be capable of causing rapid changes to key properties of temporal community architecture

Combining SPR with atomic-force microscopy enables single-molecule insights into activation and suppression of the complement cascade

This work was supported by Leverhulme Trust Grant RPG-2015-109.Activation and suppression of the complement system compete on every serum-exposed surface, host or foreign. Potentially harmful outcomes of this competition depend on surface molecules through mechanisms that remain incompletely understood. Combining surface plasmon resonance (SPR) with atomic force microscopy (AFM), here we studied two complement system proteins at the single-molecule level: C3b, the proteolytically activated form of C3, and factor H (FH), the surface-sensing C3b-binding complement regulator. We used SPR to monitor complement initiation occurring through a positive-feedback loop wherein surface-deposited C3b participates in convertases that cleave C3, thereby depositing more C3b. Over multiple cycles of flowing factor B, factor D, and C3 over the SPR chip, we amplified C3b from ∼20 to ∼220 molecules·μm−2. AFM revealed C3b clusters of up to 20 molecules and solitary C3b molecules deposited up to 200 nm away from the clusters. A force of 0.17 ± 0.02 nanonewtons was needed to pull a single FH molecule, anchored to the AFM probe, from its complex with surface-attached C3b. The extent to which FH molecules stretched before detachment varied widely among complexes. Performing force-distance measurements with FH(D1119G), a variant lacking one of the C3b-binding sites and causing atypical hemolytic uremic syndrome, we found that it detached more uniformly and easily. In further SPR experiments, KD values between FH and C3b on a custom-made chip surface were 5-fold tighter than on commercial chips and similar to those on erythrocytes. These results suggest that the chemistry at the surface on which FH acts drives conformational adjustments that are functionally critical.Publisher PDFPeer reviewe

Pathways to ischemic neuronal cell death: are sex differences relevant?

We have known for some time that the epidemiology of human stroke is sexually dimorphic until late in life, well beyond the years of reproductive senescence and menopause. Now, a new concept is emerging: the mechanisms and outcome of cerebral ischemic injury are influenced strongly by biological sex as well as the availability of sex steroids to the brain. The principal mammalian estrogen (17 β estradiol or E2) is neuroprotective in many types of brain injury and has been the major focus of investigation over the past several decades. However, it is becoming increasingly clear that although hormones are a major contributor to sex-specific outcomes, they do not fully account for sex-specific responses to cerebral ischemia. The purpose of this review is to highlight recent studies in cell culture and animal models that suggest that genetic sex determines experimental stroke outcome and that divergent cell death pathways are activated after an ischemic insult. These sex differences need to be identified if we are to develop efficacious neuroprotective agents for use in stroke patients

Neuroactive substances specifically modulate rhythmic body contractions in the nerveless metazoon Tethya wilhelma (Demospongiae, Porifera)

BACKGROUND: Sponges (Porifera) are nerve- and muscleless metazoa, but display coordinated motor reactions. Therefore, they represent a valuable phylum to investigate coordination systems, which evolved in a hypothetical Urmetazoon prior to the central nervous system (CNS) of later metazoa. We have chosen the contractile and locomotive species Tethya wilhelma (Demospongiae, Hadromerida) as a model system for our research, using quantitative analysis based on digital time lapse imaging. In order to evaluate candidate coordination pathways, we extracorporeally tested a number of chemical messengers, agonists and antagonists known from chemical signalling pathways in animals with CNS. RESULTS: Sponge body contraction of T. wilhelma was induced by caffeine, glycine, serotonine, nitric oxide (NO) and extracellular cyclic adenosine monophosphate (cAMP). The induction by glycine and cAMP followed patterns varying from other substances. Induction by cAMP was delayed, while glycine lead to a bi-phasic contraction response. The frequency of the endogenous contraction rhythm of T. wilhelma was significantly decreased by adrenaline and NO, with the same tendency for cAMP and acetylcholine. In contrast, caffeine and glycine increased the contraction frequency. The endogenous rhythm appeared irregular during application of caffeine, adrenaline, NO and cAMP. Caffeine, glycine and NO attenuated the contraction amplitude. All effects on the endogenous rhythm were neutralised by the washout of the substances from the experimental reactor system. CONCLUSION: Our study demonstrates that a number of chemical messengers, agonists and antagonists induce contraction and/or modulate the endogenous contraction rhythm and amplitude of our nerveless model metazoon T. wilhelma. We conclude that a relatively complex system of chemical messengers regulates the contraction behaviour through auto- and paracrine signalling, which is presented in a hypothetical model. We assume that adrenergic, adenosynergic and glycinergic pathways, as well as pathways based on NO and extracellular cAMP are candidates for the regulation and timing of the endogenous contraction rhythm within pacemaker cells, while GABA, glutamate and serotonine are candidates for the direct coordination of the contractile cells

Monitoring an Alien Invasion: DNA Barcoding and the Identification of Lionfish and Their Prey on Coral Reefs of the Mexican Caribbean

BACKGROUND: In the Mexican Caribbean, the exotic lionfish Pterois volitans has become a species of great concern because of their predatory habits and rapid expansion onto the Mesoamerican coral reef, the second largest continuous reef system in the world. This is the first report of DNA identification of stomach contents of lionfish using the barcode of life reference database (BOLD). METHODOLOGY/PRINCIPAL FINDINGS: We confirm with barcoding that only Pterois volitans is apparently present in the Mexican Caribbean. We analyzed the stomach contents of 157 specimens of P. volitans from various locations in the region. Based on DNA matches in the Barcode of Life Database (BOLD) and GenBank, we identified fishes from five orders, 14 families, 22 genera and 34 species in the stomach contents. The families with the most species represented were Gobiidae and Apogonidae. Some prey taxa are commercially important species. Seven species were new records for the Mexican Caribbean: Apogon mosavi, Coryphopterus venezuelae, C. thrix, C. tortugae, Lythrypnus minimus, Starksia langi and S. ocellata. DNA matches, as well as the presence of intact lionfish in the stomach contents, indicate some degree of cannibalism, a behavior confirmed in this species by the first time. We obtained 45 distinct crustacean prey sequences, from which only 20 taxa could be identified from the BOLD and GenBank databases. The matches were primarily to Decapoda but only a single taxon could be identified to the species level, Euphausia americana. CONCLUSIONS/SIGNIFICANCE: This technique proved to be an efficient and useful method, especially since prey species could be identified from partially-digested remains. The primary limitation is the lack of comprehensive coverage of potential prey species in the region in the BOLD and GenBank databases, especially among invertebrates