Use of principal components to aggregate rare variants in case-control and family-based association studies in the presence of multiple covariates

Rare variants may help to explain some of the missing heritability of complex diseases. Technological advances in next-generation sequencing give us the opportunity to test this hypothesis. We propose two new methods (one for case-control studies and one for family-based studies) that combine aggregated rare variants and common variants located within a region through principal components analysis and allow for covariate adjustment. We analyzed 200 replicates consisting of 209 case subjects and 488 control subjects and compared the results to weight-based and step-up aggregation methods. The principal components and collapsing method showed an association between the gene FLT1 and the quantitative trait Q1 (P<10−30) in a fraction of the computation time of the other methods. The proposed family-based test has inconclusive results. The two methods provide a fast way to analyze simultaneously rare and common variants at the gene level while adjusting for covariates. However, further evaluation of the statistical efficiency of this approach is warranted

A self-organized model for cell-differentiation based on variations of molecular decay rates

Systemic properties of living cells are the result of molecular dynamics governed by so-called genetic regulatory networks (GRN). These networks capture all possible features of cells and are responsible for the immense levels of adaptation characteristic to living systems. At any point in time only small subsets of these networks are active. Any active subset of the GRN leads to the expression of particular sets of molecules (expression modes). The subsets of active networks change over time, leading to the observed complex dynamics of expression patterns. Understanding of this dynamics becomes increasingly important in systems biology and medicine. While the importance of transcription rates and catalytic interactions has been widely recognized in modeling genetic regulatory systems, the understanding of the role of degradation of biochemical agents (mRNA, protein) in regulatory dynamics remains limited. Recent experimental data suggests that there exists a functional relation between mRNA and protein decay rates and expression modes. In this paper we propose a model for the dynamics of successions of sequences of active subnetworks of the GRN. The model is able to reproduce key characteristics of molecular dynamics, including homeostasis, multi-stability, periodic dynamics, alternating activity, differentiability, and self-organized critical dynamics. Moreover the model allows to naturally understand the mechanism behind the relation between decay rates and expression modes. The model explains recent experimental observations that decay-rates (or turnovers) vary between differentiated tissue-classes at a general systemic level and highlights the role of intracellular decay rate control mechanisms in cell differentiation.Comment: 16 pages, 5 figure

How consistent are the transcriptome changes associated with cold acclimation in two species of the Drosophila virilis group?

This work was financially support by a Marie Curie Initial Training Network grant, “Understanding the evolutionary origin of biological diversity” (ITN-2008–213780 SPECIATION), grants from the Academy of Finland to A.H. (project 132619) and M.K. (projects 268214 and 272927), a grant from NERC, UK to M.G.R. (grant NE/J020818/1), and NERC, UK PhD studentship to D.J.P. (NE/I528634/1).For many organisms the ability to cold acclimate with the onset of seasonal cold has major implications for their fitness. In insects, where this ability is widespread, the physiological changes associated with increased cold tolerance have been well studied. Despite this, little work has been done to trace changes in gene expression during cold acclimation that lead to an increase in cold tolerance. We used an RNA-Seq approach to investigate this in two species of the Drosophila virilis group. We found that the majority of genes that are differentially expressed during cold acclimation differ between the two species. Despite this, the biological processes associated with the differentially expressed genes were broadly similar in the two species. These included: metabolism, cell membrane composition, and circadian rhythms, which are largely consistent with previous work on cold acclimation/cold tolerance. In addition, we also found evidence of the involvement of the rhodopsin pathway in cold acclimation, a pathway that has been recently linked to thermotaxis. Interestingly, we found no evidence of differential expression of stress genes implying that long-term cold acclimation and short-term stress response may have a different physiological basis.PostprintPeer reviewe

A CLINICAL STUDY OF INHALANT ANAESTHESIA IN DOGS

A clinical trial was undertaken using three different inhalant anaesthetic agents and one intravenous anaesthetic agent in dogs undergoing routine desexing surgery. Healthy adult dogs undergoing either ovariohysterectomy or castration were assessed as to their demeanour, with the more excitable dogs being placed in groups receiving premedication with acepromazine and morphine. All dogs were then randomly assigned an anaesthetic agent for induction of general anaesthesia. The agents were the inhalants halothane, isoflurane and sevoflurane, and the intravenous agent propofol. Inhalant inductions were undertaken using a tight fitting mask attached to a standard anaesthetic machine with a rebreathing circuit, with the maximum dose of inhalant available from a standard vaporiser. Propofol inductions were undertaken via intravenous catheter. Dogs induced with propofol were randomly assigned one of the three inhalant agents for maintenance. Those induced by inhalant agent were maintained using the same agent. The surgical procedure was undertaken in standard fashion, as was recovery from anaesthesia. All dogs received the non-steroidal anti-inflammatory agent meloxicam. Data collection was divided into three stages: induction, maintenance, and recovery from anaesthesia. Variables measured at induction of anaesthesia were time to intubation, number of intubation attempts, tolerance of mask, quality of induction and quality of transfer to the maintenance stage. Standard variables for monitoring of anaesthesia were recorded throughout the maintenance of anaesthesia. Variables measured at recovery were time to righting, time to standing and quality of recovery. The mean time to intubation when using the newer inhalant sevoflurane (196.2 ± 14.8sec, mean ± SE) was not significantly different to that for halothane (221.4 ± 14.0sec) or isoflurane (172.4 ± 15.0sec). Time to intubation with isoflurane was significantly faster than with halothane. Mean time to intubation with propofol (85.4 ± 7.7sec) was significantly faster than that for any of the three inhalants. Choice of inhalant had no effect on quality of induction. The use of premedication significantly improved the quality of induction. The use of propofol for induction likewise significantly improved the quality of induction. Standard cardiorespiratory variables measured during the maintenance phase of anaesthesia remained within normal clinical ranges for all three inhalants, and were therefore not further analysed. Choice of inhalant agent had no significant effect on the time to righting or standing in recovery. The use of propofol for induction had no effect on these variables. Animals placed in groups receiving premedication had significantly longer times to righting and standing. The oesophageal temperature at the end of the procedure had a significant effect on times to righting and standing, with lower temperatures contributing to slower recoveries. Independent of procedure time, male dogs had shorter times to righting than female dogs

Predicting species dominance shifts across elevation gradients in mountain forests in Greece under a warmer and drier climate

The Mediterranean Basin is expected to face warmer and drier conditions in the future, following projected increases in temperature and declines in precipitation. The aim of this study is to explore how forests dominated by Abies borisii-regis, Abies cephalonica, Fagus sylvatica, Pinus nigra and Quercus frainetto will respond under such conditions. We combined an individual-based model (GREFOS), with a novel tree ring data set in order to constrain tree diameter growth and to account for inter- and intraspecific growth variability. We used wood density data to infer tree longevity, taking into account inter- and intraspecific variability. The model was applied at three 500-m-wide elevation gradients at Taygetos in Peloponnese, at Agrafa on Southern Pindos and at Valia Kalda on Northern Pindos in Greece. Simulations adequately represented species distribution and abundance across the elevation gradients under current climate. We subsequently used the model to estimate species and functional trait shifts under warmer and drier future conditions based on the IPCC A1B scenario. In all three sites, a retreat of less drought-tolerant species and an upward shift of more drought-tolerant species were simulated. These shifts were also associated with changes in two key functional traits, in particular maximum radial growth rate and wood density. Drought-tolerant species presented an increase in their average maximal growth and decrease in their average wood density, in contrast to less drought-tolerant species

High-Throughput NMR Assessment of the Tertiary Structure of Food Allergens

In vitro component-resolved diagnosis of food allergy requires purified allergens that have to meet high standards of quality. These include the authentication of their conformation, which is relevant for the recognition by specific IgE antibodies from allergic patients. Therefore, highly sensitive and reliable screening methods for the analysis of proteins/allergens are required to assess their structural integrity. In the present study one-dimensional 1H Nuclear Magnetic Resonance (1D 1H-NMR) analysis was adopted for the assessment of overall structural and dynamic properties and authentication of a set of relevant food allergens, including non-specific lipid transfer proteins from apple, peach and hazelnut, 7/8S seed storage globulins from hazelnut and peanut, 11S seed storage globulins from hazelnut and peanut, caseins from cows' and goats' milk and tropomyosin from shrimp.Two sets of 1D 1H-NMR experiments, using 700 MHz and 600 MHz instruments at 298 K were carried out to determine the presence and the extent of tertiary structure. Structural similarity among members of the individual allergen families was also assessed and changes under thermal stress investigated. The nuclear magnetic resonance (NMR) results were compared with structural information available either from the literature, Protein Data Bank entries, or derived from molecular models.1D (1)H-NMR analysis of food allergens allowed their classification into molecules with rigid, extended and ordered tertiary structures, molecules without a rigid tertiary structure and molecules which displayed both features. Differences in thermal stability were also detected. In summary, 1D (1)H-NMR gives insights into molecular fold of proteins and offers an independent method for assessing structural properties of proteins

Insect Repellents: Modulators of Mosquito Odorant Receptor Activity

Background: DEET, 2-undecanone (2-U), IR3535 and Picaridin are widely used as insect repellents to prevent interactions between humans and many arthropods including mosquitoes. Their molecular action has only recently been studied, yielding seemingly contradictory theories including odorant-dependent inhibitory and odorant-independent excitatory activities on insect olfactory sensory neurons (OSNs) and odorant receptor proteins (ORs). Methodology/Principal Findings: Here we characterize the action of these repellents on two Aedes aegypti ORs, AaOR2 and AaOR8, individually co-expressed with the common co-receptor AaOR7 in Xenopus oocytes; these ORs are respectively activated by the odors indole (AaOR2) and (R)-(2)-1-octen3-ol (AaOR8), odorants used to locate oviposition sites and host animals. In the absence of odorants, DEET activates AaOR2 but not AaOR8, while 2-U activates AaOR8 but not AaOR2; IR3535 and Picaridin do not activate these ORs. In the presence of odors, DEET strongly inhibits AaOR8 but not AaOR2, while 2-U strongly inhibits AaOR2 but not AaOR8; IR3535 and Picaridin strongly inhibit both ORs. Conclusions/Significance: These data demonstrate that repellents can act as olfactory agonists or antagonists, thus modulating OR activity, bringing concordance to conflicting models

Passive immunoprophylaxis and therapy with humanized monoclonal antibody specific for influenza A H5 hemagglutinin in mice

BACKGROUND: Highly pathogenic avian H5N1 influenza virus is a major public health concern. Given the lack of effective vaccine and recent evidence of antiviral drug resistance in some isolates, alternative strategies for containment of a possible future pandemic are needed. Humanized monoclonal antibodies (mAbs) that neutralize H5N1 virus could be used as prophylaxis and treatment to aid in the containment of such a pandemic. METHODS: Neutralizing mAbs against H5 hemagglutinin were humanized and introduced into C57BL/6 mice (1, 5, or 10 mg/kg bodyweight) one day prior to-, one day post- and three days post-lethal challenge with H5N1 A/Vietnam/1203/04 virus. Efficacy was determined by observation of weight loss as well as survival. RESULTS: Two mAbs neutralizing for antigenically variant H5N1 viruses, A/Vietnam/1203/04 and A/Hong Kong/213/03 were identified and humanized without loss of specificity. Both antibodies exhibited prophylactic efficacy in mice, however, VN04-2-huG1 performed better requiring only 1 mg/kg bodyweight for complete protection. When used to treat infection VN04-2-huG1 was also completely protective, even when introduced three days post infection, although higher dose of antibody was required. CONCLUSION: Prophylaxis and treatment using neutralizing humanized mAbs is efficacious against lethal challenge with A/Vietnam/1203/04, providing proof of principle for the use of passive antibody therapy as a containment option in the event of pandemic influenza

To what extent can behaviour change techniques be identified within an adaptable implementation package for primary care? A prospective directed content analysis

Interpreting evaluations of complex interventions can be difficult without sufficient description of key intervention content. We aimed to develop an implementation package for primary care which could be delivered using typically available resources and could be adapted to target determinants of behaviour for each of four quality indicators: diabetes control, blood pressure control, anticoagulation for atrial fibrillation and risky prescribing. We describe the development and prospective verification of behaviour change techniques (BCTs) embedded within the adaptable implementation packages