Rapid screening for phenotype-genotype associations by linear transformations of genomic evaluations

Background

Currently, association studies are analysed using statistical mixed models, with marker effects estimated by a linear transformation of genomic breeding values. The variances of marker effects are needed when performing the tests of association. However, approaches used to estimate the parameters rely on a prior variance or on a constant estimate of the additive variance. Alternatively, we propose a standardized test of association using the variance of each marker effect, which generally differ among each other. Random breeding values from a mixed model including fixed effects and a genomic covariance matrix are linearly transformed to estimate the marker effects.

Results

The standardized test was neither conservative nor liberal with respect to type I error rate (false-positives), compared to a similar test using Predictor Error Variance, a method that was too conservative. Furthermore, genomic predictions are solved efficiently by the procedure, and the p-values are virtually identical to those calculated from tests for one marker effect at a time. Moreover, the standardized test reduces computing time and memory requirements.The following steps are used to locate genome segments displaying strong association. The marker with the highest − log(p-value) in each chromosome is selected, and the segment is expanded one Mb upstream and one Mb downstream of the marker. A genomic matrix is calculated using the information from those markers only, which is used as the variance-covariance of the segment effects in a model that also includes fixed effects and random genomic breeding values. The likelihood ratio is then calculated to test for the effect in every chromosome against a reduced model with fixed effects and genomic breeding values. In a case study with pigs, a significant segment from chromosome 6 explained 11% of total genetic variance.

Conclusions

The standardized test of marker effects using their own variance helps in detecting specific genomic regions involved in the additive variance, and in reducing false positives. Moreover, genome scanning of candidate segments can be used in meta-analyses of genome-wide association studies, as it enables the detection of specific genome regions that affect an economically relevant trait when using multiple populations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-246) contains supplementary material, which is available to authorized users.     

Comparative genomics and transcriptomics in ants provide new insights into the evolution and function of odorant binding and chemosensory proteins

Background

The complex societies of ants and other social insects rely on sophisticated chemical communication. Two families of small soluble proteins, the odorant binding and chemosensory proteins (OBPs and CSPs), are believed to be important in insect chemosensation. To better understand the role of these proteins in ant olfaction, we examined their evolution and expression across the ants using phylogenetics and sex- and tissue-specific RNA-seq.

Results

We find that subsets of both OBPs and CSPs are expressed in the antennae, contradicting the previous hypothesis that CSPs have replaced OBPs in ant olfaction. Both protein families have several highly conserved clades with a single ortholog in all eusocial hymenopterans, as well as clades with more dynamic evolution and many taxon-specific radiations. The dynamically evolving OBPs and CSPs have been hypothesized to function in chemical communication. Intriguingly, we find that seven members of the conserved clades are expressed specifically in the antennae of the clonal raider ant Cerapachys biroi, whereas only one dynamically evolving CSP is antenna specific. The orthologs of the conserved, antenna-specific C. biroi genes are also expressed in antennae of the ants Camponotus floridanus and Harpegnathos saltator, indicating that antenna-specific expression of these OBPs and CSPs is conserved across ants. Most members of the dynamically evolving clades in both protein families are expressed primarily in non-chemosensory tissues and thus likely do not fulfill chemosensory functions.

Conclusions

Our results identify candidate OBPs and CSPs that are likely involved in conserved aspects of ant olfaction, and suggest that OBPs and CSPs may not rapidly evolve to recognize species-specific signals.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-718) contains supplementary material, which is available to authorized users.     

Defining adult asthma endotypes by clinical features and patterns of volatile organic compounds in exhaled air

Background

Several classifications of adult asthma patients using cluster analyses based on clinical and demographic information has resulted in clinical phenotypic clusters that do not address molecular mechanisms. Volatile organic compounds (VOC) in exhaled air are released during inflammation in response to oxidative stress as a result of activated leukocytes. VOC profiles in exhaled air could distinguish between asthma patients and healthy subjects. In this study, we aimed to classify new asthma endotypes by combining inflammatory mechanisms investigated by VOC profiles in exhaled air and clinical information of asthma patients.

Methods

Breath samples were analyzed for VOC profiles by gas chromatography–mass spectrometry from asthma patients (n = 195) and healthy controls (n = 40). A total of 945 determined compounds were subjected to discriminant analysis to find those that could discriminate healthy from asthmatic subjects. 2-step cluster analysis based on clinical information and VOCs in exhaled air were used to form asthma endotypes.

Results

We identified 16 VOCs, which could distinguish between healthy and asthma subjects with a sensitivity of 100% and a specificity of 91.1%. Cluster analysis based on VOCs in exhaled air and the clinical parameters FEV1, FEV1 change after 3 weeks of hospitalization, allergic sensitization, Junipers symptoms score and asthma medications resulted in the formation of 7 different asthma endotype clusters. We identified asthma clusters with different VOC profiles but similar clinical characteristics and endotypes with similar VOC profiles, but distinct clinical characteristics.

Conclusion

This study demonstrates that both, clinical presentation of asthma and inflammatory mechanisms in the airways should be considered for classification of asthma subtypes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12931-014-0136-8) contains supplementary material, which is available to authorized users.     

Conceptual-level workflow modeling of scientific experiments using NMR as a case study

Background  

Scientific workflows improve the process of scientific experiments by making computations explicit, underscoring data flow, and emphasizing the participation of humans in the process when intuition and human reasoning are required. Workflows for experiments also highlight transitions among experimental phases, allowing intermediate results to be verified and supporting the proper handling of semantic mismatches and different file formats among the various tools used in the scientific process. Thus, scientific workflows are important for the modeling and subsequent capture of bioinformatics-related data. While much research has been conducted on the implementation of scientific workflows, the initial process of actually designing and generating the workflow at the conceptual level has received little consideration.     

Acute effects of prolonged intermittent low-intensity isometric warm-up schemes on jump,sprint, and agility performance in collegiate soccer players

The aim of the present study was to compare the effects of different warm-up interventions on jump, sprint and agility performance in collegiate soccer players. Twenty-one healthy male college soccer players (age: 20.14 ± 1.65 years; body height: 179.9 ± 8.34 cm; body mass: 74.4 ± 13.0 kg; % body fat: 9.45 ± 4.8) participated in the study. Subjects underwent four different randomized warm-up protocols separated by at least 48 hours. The warm-up schemes were: 1. no conditioning contraction protocol (NCC); 2. dynamic stretching (DS); 3. prolonged intermittent low-intensity isometric exercise (ST); and, 4. ST with an additional external load equal to 30% of body weight (ST + 30% BW). All interventions were preceded by a general warm-up. Results from one-way repeated measures ANOVA demonstrated a significant difference in countermovement jump (CMJ) at F(3,60) = 10.2, ηρ² = 0.337, p < 0.01. Post hoc analysis revealed a significant difference in CMJ performance in DS when compared to NCC and ST + 30% BW. No significant difference in CMJ was observed between DS and ST. CMJ scores in NCC, ST, and ST + 30% BW were non-significant. There was a significant difference in speed; F(3, 60) = 6.61, ηρ² = 0.248, p < 0.01. Post hoc analysis revealed significantly better time in DS than NCC and ST. However, no difference in speed was observed between DS and ST + 30% BW. Similarly, speed was similar in NCC, ST and ST + 30% BW. A significant difference in agility performance was also observed; F(3, 60) = 24.1, ηρ²= 0.546, p < 0.01. Post hoc analysis revealed significantly greater performance gains in DS than NCC. No significant difference in agility was observed in DS, ST and ST + 30% BW. In conclusion, a prolonged intermittent low-intensity isometric protocol using bodyweight only showed similar benefits with dynamic stretching in countermovement jump performance. When the same isometric condition with additional load equal to 30% of bodyweight was applied, effects in speed and agility were similar to dynamic stretching.     

CFAP54 is required for proper ciliary motility and assembly of the central pair apparatus in mice

Motile cilia and flagella play critical roles in fluid clearance and cell motility, and dysfunction commonly results in the pediatric syndrome primary ciliary dyskinesia (PCD). CFAP221, also known as PCDP1, is required for ciliary and flagellar function in mice and Chlamydomonas reinhardtii, where it localizes to the C1d projection of the central microtubule apparatus and functions in a complex that regulates flagellar motility in a calcium-dependent manner. We demonstrate that the genes encoding the mouse homologues of the other C. reinhardtii C1d complex members are primarily expressed in motile ciliated tissues, suggesting a conserved function in mammalian motile cilia. The requirement for one of these C1d complex members, CFAP54, was identified in a mouse line with a gene-trapped allele. Homozygous mice have PCD characterized by hydrocephalus, male infertility, and mucus accumulation. The infertility results from defects in spermatogenesis. Motile cilia have a structural defect in the C1d projection, indicating that the C1d assembly mechanism requires CFAP54. This structural defect results in decreased ciliary beat frequency and perturbed cilia-driven flow. This study identifies a critical role for CFAP54 in proper assembly and function of mammalian cilia and flagella and establishes the gene-trapped allele as a new model of PCD.     

Kinesin Motor Enzymology: Chemistry,Structure, and Physics of Nanoscale Molecular Machines

Biophysical Reviews - Molecular motors are enzymes that convert chemical potential energy into controlled kinetic energy for mechanical work inside cells. Understanding the biophysics of these...     

Performance evaluation of cetacean species distribution models developed using generalized additive models and boosted regression trees

Species distribution models (SDMs) are important management tools for highly mobile marine species because they provide spatially and temporally explicit information on animal distribution. Two prevalent modeling frameworks used to develop SDMs for marine species are generalized additive models (GAMs) and boosted regression trees (BRTs), but comparative studies have rarely been conducted; most rely on presence‐only data; and few have explored how features such as species distribution characteristics affect model performance. Since the majority of marine species BRTs have been used to predict habitat suitability, we first compared BRTs to GAMs that used presence/absence as the response variable. We then compared results from these habitat suitability models to GAMs that predict species density (animals per km²) because density models built with a subset of the data used here have previously received extensive validation. We compared both the explanatory power (i.e., model goodness of fit) and predictive power (i.e., performance on a novel dataset) of the GAMs and BRTs for a taxonomically diverse suite of cetacean species using a robust set of systematic survey data (1991–2014) within the California Current Ecosystem. Both BRTs and GAMs were successful at describing overall distribution patterns throughout the study area for the majority of species considered, but when predicting on novel data, the density GAMs exhibited substantially greater predictive power than both the presence/absence GAMs and BRTs, likely due to both the different response variables and fitting algorithms. Our results provide an improved understanding of some of the strengths and limitations of models developed using these two methods. These results can be used by modelers developing SDMs and resource managers tasked with the spatial management of marine species to determine the best modeling technique for their question of interest.     

Adenosine activation of A(2B) receptor(s) is essential for stimulated epithelial ciliary motility and clearance

Mucociliary clearance, vital to lung clearance, is dependent on cilia beat frequency (CBF), coordination of cilia, and the maintenance of periciliary fluid. Adenosine, the metabolic breakdown product of ATP, is an important modulator of ciliary motility. However, the contributions of specific adenosine receptors to key airway ciliary motility processes are unclear. We hypothesized that adenosine modulates ciliary motility via activation of its cell surface receptors (A(1), A(2A), A(2B), or A(3)). To test this hypothesis, mouse tracheal rings (MTRs) excised from wild-type and adenosine receptor knockout mice (A(1), A(2A), A(2B), or A(3), respectively), and bovine ciliated bronchial epithelial cells (BBECs) were stimulated with known cilia activators, isoproterenol (ISO; 10 μM) and/or procaterol (10 μM), in the presence or absence of 5'-(N-ethylcarboxamido) adenosine (NECA), a nonselective adenosine receptor agonist [100 nM (A(1), A(2A), A(3)); 10 μM (A(2B))], and CBF was measured. Cells and MTRs were also stimulated with NECA (100 nM or 10 μM) in the presence and absence of adenosine deaminase inhibitor, erythro-9- (2-hydroxy-3-nonyl) adenine hydrochloride (10 μM). Both ISO and procaterol stimulated CBF in untreated cells and/or MTRs from both wild-type and adenosine knockout mice by ~3 Hz. Likewise, CBF significantly increased ~2-3 Hz in BBECs and wild-type MTRs stimulated with NECA. MTRs from A(1), A(2A), and A(3) knockout mice stimulated with NECA also demonstrated an increase in CBF. However, NECA failed to stimulate CBF in MTRs from A(2B) knockout mice. To confirm the mechanism by which adenosine modulates CBF, protein kinase activity assays were conducted. The data revealed that NECA-stimulated CBF is mediated by the activation of cAMP-dependent PKA. Collectively, these data indicate that purinergic stimulation of CBF requires A(2B) adenosine receptor activation, likely via a PKA-dependent pathway.     

Withdrawal of inhaled corticosteroids in individuals with COPD - a systematic review and comment on trial methodology

Inhaled corticosteroids (ICS) reduce COPD exacerbation frequency and slow decline in health related quality of life but have little effect on lung function, do not reduce mortality, and increase the risk of pneumonia. We systematically reviewed trials in which ICS have been withdrawn from patients with COPD, with the aim of determining the effect of withdrawal, understanding the differing results between trials, and making recommendations for improving methodology in future trials where medication is withdrawn. Trials were identified by two independent reviewers using MEDLINE, EMBASE and CINAHL, citations of identified studies were checked, and experts contacted to identify further studies. Data extraction was completed independently by two reviewers. The methodological quality of each trial was determined by assessing possible sources of systematic bias as recommended by the Cochrane collaboration. We included four trials; the quality of three was adequate. In all trials, outcomes were generally worse for patients who had had ICS withdrawn, but differences between outcomes for these patients and patients who continued with medication were mostly small and not statistically significant. Due to data paucity we performed only one meta-analysis; this indicated that patients who had had medication withdrawn were 1.11 (95% CI 0.84 to 1.46) times more likely to have an exacerbation in the following year, but the definition of exacerbations was not consistent between the three trials, and the impact of withdrawal was smaller in recent trials which were also trials conducted under conditions that reflected routine practice. There is no evidence from this review that withdrawing ICS in routine practice results in important deterioration in patient outcomes. Furthermore, the extent of increase in exacerbations depends on the way exacerbations are defined and managed and may depend on the use of other medication. In trials where medication is withdrawn, investigators should report other medication use, definitions of exacerbations and management of patients clearly. Intention to treat analyses should be used and interpreted appropriately.