CNVassoc: Association analysis of CNV data using R

Background: Copy number variants (CNV) are a potentially important component of the genetic contribution to risk of common complex diseases. Analysis of the association between CNVs and disease requires that uncertainty in CNV copy-number calls, which can be substantial, be taken into account; failure to consider this uncertainty can lead to biased results. Therefore, there is a need to develop and use appropriate statistical tools. To address this issue, we have developed CNVassoc, an R package for carrying out association analysis of common copy number variants in population-based studies. This package includes functions for testing for association with different classes of response variables (e.g. class status, censored data, counts) under a series of study designs (case-control, cohort, etc) and inheritance models, adjusting for covariates. The package includes functions for inferring copy number (CNV genotype calling), but can also accept copy number data generated by other algorithms (e.g. CANARY, CGHcall, IMPUTE). Results: Here we present a new R package, CNVassoc, that can deal with different types of CNV arising from different platforms such as MLPA o aCGH. Through a real data example we illustrate that our method is able to incorporate uncertainty in the association process. We also show how our package can also be useful when analyzing imputed data when analyzing imputed SNPs. Through a simulation study we show that CNVassoc outperforms CNVtools in terms of computing time as well as in convergence failure rate. Conclusions: We provide a package that outperforms the existing ones in terms of modelling flexibility, power, convergence rate, ease of covariate adjustment, and requirements for sample size and signal quality. Therefore, we offer CNVassoc as a method for routine use in CNV association studiesThis work has been supported by the Spanish Ministry of Science and Innovation (MTM2008-02457 to JRG, BIO2009-12458 to RD-U and statistical genetics network MTM2010-09526-E (subprograma MTM) to JRG, IS, GL and RD-U). GL is supported by the Juan de la Cierva Program of the Spanish Ministry of Science and Innovation

Thermodynamic analysis of black hole solutions in gravitating nonlinear electrodynamics

We perform a general study of the thermodynamic properties of static electrically charged black hole solutions of nonlinear electrodynamics minimally coupled to gravitation in three space dimensions. The Lagrangian densities governing the dynamics of these models in flat space are defined as arbitrary functions of the gauge field invariants, constrained by some requirements for physical admissibility. The exhaustive classification of these theories in flat space, in terms of the behaviour of the Lagrangian densities in vacuum and on the boundary of their domain of definition, defines twelve families of admissible models. When these models are coupled to gravity, the flat space classification leads to a complete characterization of the associated sets of gravitating electrostatic spherically symmetric solutions by their central and asymptotic behaviours. We focus on nine of these families, which support asymptotically Schwarzschild-like black hole configurations, for which the thermodynamic analysis is possible and pertinent. In this way, the thermodynamic laws are extended to the sets of black hole solutions of these families, for which the generic behaviours of the relevant state variables are classified and thoroughly analyzed in terms of the aforementioned boundary properties of the Lagrangians. Moreover, we find universal scaling laws (which hold and are the same for all the black hole solutions of models belonging to any of the nine families) running the thermodynamic variables with the electric charge and the horizon radius. These scale transformations form a one-parameter multiplicative group, leading to universal "renormalization group"-like first-order differential equations. The beams of characteristics of these equations generate the full set of black hole states associated to any of these gravitating nonlinear electrodynamics...Comment: 51 single column pages, 19 postscript figures, 2 tables, GRG tex style; minor corrections added; final version appearing in General Relativity and Gravitatio

Comparison of flipped learning and traditional lecture method for teaching digestive system diseases in undergraduate medicine: A prospective non-randomized controlled trial

Introduction: This study examined the effects of a large-scale flipped learning (FL) approach in an undergraduate course of Digestive System Diseases. Methods: This prospective non-randomized trial recruited 404 students over three academic years. In 2016, the course was taught entirely in a Traditional Lecture (TL) style, in 2017 half of the course (Medical topics) was replaced by FL while the remaining half (Surgical topics) was taught by TL and in 2018, the whole course was taught entirely by FL. Academic performance, class attendance and student’s satisfaction surveys were compared between cohorts. Results: Test scores were higher in the FL module (Medical) than in the TL module (Surgical) in the 2017 cohort but were not different when both components were taught entirely by TL (2016) or by FL (2018). Also, FL increased the probability of reaching superior grades (scores >7.0) and improved class attendance and students’ satisfaction. Conclusion: The holistic FL model is more effective for teaching undergraduate clinical gastroenterology compared to traditional teaching methods and has a positive impact on classroom attendances

Superconformal Flavor Simplified

A simple explanation of the flavor hierarchies can arise if matter fields interact with a conformal sector and different generations have different anomalous dimensions under the CFT. However, in the original study by Nelson and Strassler many supersymmetric models of this type were considered to be 'incalculable' because the R-charges were not sufficiently constrained by the superpotential. We point out that nearly all such models are calculable with the use of a-maximization. Utilizing this, we construct the simplest vector-like flavor models and discuss their viability. A significant constraint on these models comes from requiring that the visible gauge couplings remain perturbative throughout the conformal window needed to generate the hierarchies. However, we find that there is a small class of simple flavor models that can evade this bound.Comment: 43 pages, 1 figure; V3: small corrections and clarifications, references adde

Diagnostic delay for giant cell arteritis – a systematic review and meta-analysis

Background Giant cell arteritis (GCA), if untreated, can lead to blindness and stroke. The study’s objectives were to (1) determine a new evidence-based benchmark of the extent of diagnostic delay for GCA and (2) examine the role of GCA-specific characteristics on diagnostic delay. Methods Medical literature databases were searched from inception to November 2015. Articles were included if reporting a time-period of diagnostic delay between onset of GCA symptoms and diagnosis. Two reviewers assessed the quality of the final articles and extracted data from these. Random-effects meta-analysis was used to pool the mean time-period (95% confidence interval (CI)) between GCA symptom onset and diagnosis, and the delay observed for GCA-specific characteristics. Heterogeneity was assessed by I 2 and by 95% prediction interval (PI). Results Of 4128 articles initially identified, 16 provided data for meta-analysis. Mean diagnostic delay was 9.0 weeks (95% CI, 6.5 to 11.4) between symptom onset and GCA diagnosis (I 2 = 96.0%; P < 0.001; 95% PI, 0 to 19.2 weeks). Patients with a cranial presentation of GCA received a diagnosis after 7.7 (95% CI, 2.7 to 12.8) weeks (I 2 = 98.4%; P < 0.001; 95% PI, 0 to 27.6 weeks) and those with non-cranial GCA after 17.6 (95% CI, 9.7 to 25.5) weeks (I 2 = 96.6%; P < 0.001; 95% PI, 0 to 46.1 weeks). Conclusions The mean delay from symptom onset to GCA diagnosis was 9 weeks, or longer when cranial symptoms were absent. Our research provides an evidence-based benchmark for diagnostic delay of GCA and supports the need for improved public awareness and fast-track diagnostic pathways

Network meta‐analysis of post‐exposure prophylaxis randomized clinical trials

Objectives: We performed a network meta‐analysis of PEP randomized clinical trials to evaluate the best regimen. / Methods: After MEDLINE/Pubmed search, studies were included if: (1) were randomized, (2) comparing at least 2 PEP three‐drug regimens and, (3) reported completion rates or discontinuation at 28 days. Five studies with 1105 PEP initiations were included and compared ritonavir‐boosted lopinavir (LPV/r) vs. atazanavir (ATV) (one study), cobicistat‐boosted elvitegravir (EVG/c) (one study), raltegravir (RAL) (one study) or maraviroc (MVC) (two studies). We estimated the probability of each treatment of being the best based on the evaluation of five outcomes: PEP non‐completion at day 28, PEP discontinuation due to adverse events, PEP switching due to any cause, lost to follow‐up and adverse events. / Results: Participants were mostly men who have sex with men (n = 832, 75%) with non‐occupational exposure to HIV (89.86%). Four‐hundred fifty‐four (41%) participants failed to complete their PEP course for any reason. The Odds Ratio (OR) for PEP non‐completion at day 28 in each antiretroviral compared to LPV/r was: ATV 0.95 (95% CI 0.58–1.56; EVG/c: OR 0.65 95% CI 0.30–1.37; RAL: OR 0.68 95% CI 0.41–1.13; and MVC: OR 0.69 95% CI 0.47–1.01. In addition, the rankogram showed that EVG/c had the highest probability of being the best treatment for the lowest rates in PEP non‐completion at day 28, switching, lost to follow‐up or adverse events and MVC for PEP discontinuations due to adverse events. / Conclusions: Our study shows the advantages of integrase inhibitors when used as PEP, particularly EVG as a Single‐Tablet Regimen

Bioturbation in a Declining Oxygen Environment, in situ Observations from Wormcam

Bioturbation, the displacement and mixing of sediment particles by fauna or flora, facilitates life supporting processes by increasing the quality of marine sediments. In the marine environment bioturbation is primarily mediated by infaunal organisms, which are susceptible to perturbations in their surrounding environment due to their sedentary life history traits. Of particular concern is hypoxia, dissolved oxygen (DO) concentrations ≤2.8 mg l−1, a prevalent and persistent problem that affects both pelagic and benthic fauna. A benthic observing system (Wormcam) consisting of a buoy, telemetering electronics, sediment profile camera, and water quality datasonde was developed and deployed in the Rappahannock River, VA, USA, in an area known to experience seasonal hypoxia from early spring to late fall. Wormcam transmitted a time series of in situ images and water quality data, to a website via wireless internet modem, for 5 months spanning normoxic and hypoxic periods. Hypoxia was found to significantly reduce bioturbation through reductions in burrow lengths, burrow production, and burrowing depth. Although infaunal activity was greatly reduced during hypoxic and near anoxic conditions, some individuals remained active. Low concentrations of DO in the water column limited bioturbation by infaunal burrowers and likely reduced redox cycling between aerobic and anaerobic states. This study emphasizes the importance of in situ observations for understanding how components of an ecosystem respond to hypoxia

Empirical Relationship between Intra-Purine and Intra-Pyrimidine Differences in Conserved Gene Sequences

DNA sequences seen in the normal character-based representation appear to have a formidable mixing of the four nucleotides without any apparent order. Nucleotide frequencies and distributions in the sequences have been studied extensively, since the simple rule given by Chargaff almost a century ago that equates the total number of purines to the pyrimidines in a duplex DNA sequence. While it is difficult to trace any relationship between the bases from studies in the character representation of a DNA sequence, graphical representations may provide a clue. These novel representations of DNA sequences have been useful in providing an overview of base distribution and composition of the sequences and providing insights into many hidden structures. We report here our observation based on a graphical representation that the intra-purine and intra-pyrimidine differences in sequences of conserved genes generally follow a quadratic distribution relationship and show that this may have arisen from mutations in the sequences over evolutionary time scales. From this hitherto undescribed relationship for the gene sequences considered in this report we hypothesize that such relationships may be characteristic of these sequences and therefore could become a barrier to large scale sequence alterations that override such characteristics, perhaps through some monitoring process inbuilt in the DNA sequences. Such relationship also raises the possibility of intron sequences playing an important role in maintaining the characteristics and could be indicative of possible intron-late phenomena

TRY plant trait database - enhanced coverage and open access

Plant traits-the morphological, anatomical, physiological, biochemical and phenological characteristics of plants-determine how plants respond to environmental factors, affect other trophic levels, and influence ecosystem properties and their benefits and detriments to people. Plant trait data thus represent the basis for a vast area of research spanning from evolutionary biology, community and functional ecology, to biodiversity conservation, ecosystem and landscape management, restoration, biogeography and earth system modelling. Since its foundation in 2007, the TRY database of plant traits has grown continuously. It now provides unprecedented data coverage under an open access data policy and is the main plant trait database used by the research community worldwide. Increasingly, the TRY database also supports new frontiers of trait-based plant research, including the identification of data gaps and the subsequent mobilization or measurement of new data. To support this development, in this article we evaluate the extent of the trait data compiled in TRY and analyse emerging patterns of data coverage and representativeness. Best species coverage is achieved for categorical traits-almost complete coverage for 'plant growth form'. However, most traits relevant for ecology and vegetation modelling are characterized by continuous intraspecific variation and trait-environmental relationships. These traits have to be measured on individual plants in their respective environment. Despite unprecedented data coverage, we observe a humbling lack of completeness and representativeness of these continuous traits in many aspects. We, therefore, conclude that reducing data gaps and biases in the TRY database remains a key challenge and requires a coordinated approach to data mobilization and trait measurements. This can only be achieved in collaboration with other initiatives