PILER-CR: Fast and accurate identification of CRISPR repeats

BACKGROUND: Sequencing of prokaryotic genomes has recently revealed the presence of CRISPR elements: short, highly conserved repeats separated by unique sequences of similar length. The distinctive sequence signature of CRISPR repeats can be found using general-purpose repeat- or pattern-finding software tools. However, the output of such tools is not always ideal for studying these repeats, and significant effort is sometimes needed to build additional tools and perform manual analysis of the output. RESULTS: We present PILER-CR, a program specifically designed for the identification and analysis of CRISPR repeats. The program executes rapidly, completing a 5 Mb genome in around 5 seconds on a current desktop computer. We validate the algorithm by manual curation and by comparison with published surveys of these repeats, finding that PILER-CR has both high sensitivity and high specificity. We also present a catalogue of putative CRISPR repeats identified in a comprehensive analysis of 346 prokaryotic genomes. CONCLUSION: PILER-CR is a useful tool for rapid identification and classification of CRISPR repeats. The software is donated to the public domain. Source code and a Linux binary are freely available at

Effect of α+-thalassaemia on episodes of fever due to malaria and other causes: a community-based cohort study in Tanzania

It is controversial to what degree α(+)-thalassaemia protects against episodes of uncomplicated malaria and febrile disease due to infections other than Plasmodium. In Tanzania, in children aged 6-60 months and height-for-age z-score < -1.5 SD (n = 612), rates of fevers due to malaria and other causes were compared between those with heterozygous or homozygotes α(+)-thalassaemia and those with a normal genotype, using Cox regression models that accounted for multiple events per child. The overall incidence of malaria was 3.0/child-year (1, 572/526 child-years); no differences were found in malaria rates between genotypes (hazard ratios, 95% CI: 0.93, 0.82-1.06 and 0.91, 0.73-1.14 for heterozygotes and homozygotes respectively, adjusted for baseline factors that were predictive for outcome). However, this association strongly depended on age: among children aged 6-17 months, those with α(+)-thalassaemia experienced episodes more frequently than those with a normal genotype (1.30, 1.02-1.65 and 1.15, 0.80-1.65 for heterozygotes and homozygotes respectively), whereas among their peers aged 18-60 months, α(+)-thalassaemia protected against malaria (0.80, 0.68-0.95 and 0.78, 0.60-1.03; p-value for interaction 0.001 and 0.10 for hetero- and homozygotes respectively). No effect was observed on non-malarial febrile episodes. In this population, the association between α(+)-thalassaemia and malaria depends on age. Our data suggest that protection by α(+)-thalassaemia is conferred by more efficient acquisition of malaria-specific immunity

Add-On Cannabidiol Treatment for Drug-Resistant Seizures in Tuberous Sclerosis Complex A Placebo-Controlled Randomized Clinical Trial

IMPORTANCE Efficacy of cannabidiol has been demonstrated in seizures associated with Lennox-Gastaut and Dravet syndromes but appears not yet to have been established in conditions with primarily focal seizures, such as tuberous sclerosis complex (TSC). OBJECTIVE To evaluate efficacy and safety of 25-mg/kg/day and 50-mg/kg/day cannabidiol dosages vs placebo against seizures associated with TSC. DESIGN, SETTING, AND PARTICIPANTS This double-blind, placebo-controlled randomized clinical trial (GWPCARE6) enrolled patients between April 6, 2016, and October 4, 2018; follow-up was completed on February 15, 2019. The trial was conducted at 46 sites in Australia, Poland, Spain, the Netherlands, United Kingdom, and United States. Eligible patients (aged 1-65 years) were those with a clinical diagnosis of TSC and medicationresistant epilepsy who had had at least 8 TSC-associated seizures during the 4-week baseline period, with at least 1 seizure occurring in at least 3 of the 4 weeks, and were currently taking at least 1 antiepileptic medication. INTERVENTIONS Patients received oral cannabidiol at 25mg/kg/day (CBD25) or 50 mg/kg/day (CBD50) or a matched placebo for 16 weeks. MAIN OUTCOMES AND MEASURES The prespecified primary outcomewas the change from baseline in number of TSC-associated seizures for cannabidiol vs placebo during the treatment period. RESULTS Of 255 patients screened for eligibility, 31 were excluded and 224 were randomized. Of the 224 included patients (median [range] age, 11.4 [1.1-56.8] years; 93 female patients [41.5%]), 75 were randomized to CBD25, 73 to CBD50, and 76 to placebo, with 201 completing treatment. The percentage reduction from baseline in the type of seizures considered the primary end point was 48.6%(95%CI, 40.4%-55.8%) for the CBD25 group, 47.5%(95%CI, 39.0%-54.8%) for the CBD50 group, and 26.5%(95%CI, 14.9%-36.5%) for the placebo group; the percentage reduction from placebo was 30.1% (95%CI, 13.9%-43.3%; P < .001) for the CBD25 group and 28.5%(95%CI, 11.9%-42.0%; nominal P = .002) for the CBD50 group. The most common adverse events were diarrhea (placebo group, 19 [25%]; CBD25 group, 23 [31%]; CBD50 group, 41 [56%]) and somnolence (placebo group, 7 [9%]; CBD25 group, 10 [13%]; CBD50 group, 19 [26%]), which occurred more frequently with cannabidiol than placebo. Eight patients in CBD25 group, 10 in CBD50 group, and 2 in the placebo group discontinued treatment because of adverse events. Twenty-eight patients taking cannabidiol (18.9%) had elevated liver transaminase levels vs none taking placebo. CONCLUSIONS AND RELEVANCE Cannabidiol significantly reduced TSC-associated seizures compared with placebo. The 25-mg/kg/day dosage had a better safety profile than the 50-mg/kg/day dosage. TRIAL REGISTRATION ClinicalTrials.gov Identifier: NCT0254476

Hot new directions for quasi-Monte Carlo research in step with applications

This article provides an overview of some interfaces between the theory of quasi-Monte Carlo (QMC) methods and applications. We summarize three QMC theoretical settings: first order QMC methods in the unit cube $[0,1]^s$ and in $\mathbb{R}^s$, and higher order QMC methods in the unit cube. One important feature is that their error bounds can be independent of the dimension $s$ under appropriate conditions on the function spaces. Another important feature is that good parameters for these QMC methods can be obtained by fast efficient algorithms even when $s$ is large. We outline three different applications and explain how they can tap into the different QMC theory. We also discuss three cost saving strategies that can be combined with QMC in these applications. Many of these recent QMC theory and methods are developed not in isolation, but in close connection with applications

The effects of weather and climate change on dengue

There is much uncertainty about the future impact of climate change on vector-borne diseases. Such uncertainty reflects the difficulties in modelling the complex interactions between disease, climatic and socioeconomic determinants. We used a comprehensive panel dataset from Mexico covering 23 years of province-specific dengue reports across nine climatic regions to estimate the impact of weather on dengue, accounting for the effects of non-climatic factors

Communication in production animal medicine: modelling a complex interaction with the example of dairy herd health medicine

<p>Abstract</p> <p>Background</p> <p>The importance of communication skills in veterinary medicine is increasingly recognised. Appropriate communication skills towards the client are of utmost importance in both companion animal practice and production animal field and consultancy work. The need for building a relationship with the client, alongside developing a structure for the consultation is widely recognised and applies to both types of veterinary practice.</p> <p>Results</p> <p>Veterinary advisory practice in production animal medicine is, however, characterised by a more complex communication on different levels. While the person-orientated communication is a permanent process between veterinarian and client with a rather personal perspective and defines the roles of interaction, the problem-orientated communication deals with emerging difficulties; the objective is to solve an acute health problem. The solution - orientated communication is a form of communication in which both veterinarian and client address longstanding situations or problems with the objective to improve herd health and subsequently productivity performance. All three forms of communication overlap.</p> <p>Conclusions</p> <p>Based on this model, it appears useful for a veterinary practice to offer both a curative and an advisory service, but to keep these two separated when deemed appropriate. In veterinary education, the strategies and techniques necessary for solution orientated communication should be included in the teaching of communication skills.</p

Individual differences in metabolomics: individualised responses and between-metabolite relationships

Many metabolomics studies aim to find ‘biomarkers’: sets of molecules that are consistently elevated or decreased upon experimental manipulation. Biological effects, however, often manifest themselves along a continuum of individual differences between the biological replicates in the experiment. Such differences are overlooked or even diminished by methods in standard use for metabolomics, although they may contain a wealth of information on the experiment. Properly understanding individual differences is crucial for generating knowledge in fields like personalised medicine, evolution and ecology. We propose to use simultaneous component analysis with individual differences constraints (SCA-IND), a data analysis method from psychology that focuses on these differences. This method constructs axes along the natural biochemical differences between biological replicates, comparable to principal components. The model may shed light on changes in the individual differences between experimental groups, but also on whether these differences correspond to, e.g., responders and non-responders or to distinct chemotypes. Moreover, SCA-IND reveals the individuals that respond most to a manipulation and are best suited for further experimentation. The method is illustrated by the analysis of individual differences in the metabolic response of cabbage plants to herbivory. The model reveals individual differences in the response to shoot herbivory, where two ‘response chemotypes’ may be identified. In the response to root herbivory the model shows that individual plants differ strongly in response dynamics. Thereby SCA-IND provides a hitherto unavailable view on the chemical diversity of the induced plant response, that greatly increases understanding of the system

Structural basis for CRISPR RNA-guided DNA recognition by Cascade

The CRISPR (clustered regularly interspaced short palindromic repeats) immune system in prokaryotes uses small guide RNAs to neutralize invading viruses and plasmids. In Escherichia coli, immunity depends on a ribonucleoprotein complex called Cascade. Here we present the composition and low-resolution structure of Cascade and show how it recognizes double-stranded DNA (dsDNA) targets in a sequence-specific manner. Cascade is a 405-kDa complex comprising five functionally essential CRISPR-associated (Cas) proteins (CasA1B2C6D1E1) and a 61-nucleotide CRISPR RNA (crRNA) with 5′-hydroxyl and 2′,3′-cyclic phosphate termini. The crRNA guides Cascade to dsDNA target sequences by forming base pairs with the complementary DNA strand while displacing the noncomplementary strand to form an R-loop. Cascade recognizes target DNA without consuming ATP, which suggests that continuous invader DNA surveillance takes place without energy investment. The structure of Cascade shows an unusual seahorse shape that undergoes conformational changes when it binds target DNA.

Genotype calling in tetraploid species from bi-allelic marker data using mixture models

<p>Abstract</p> <p>Background</p> <p>Automated genotype calling in tetraploid species was until recently not possible, which hampered genetic analysis. Modern genotyping assays often produce two signals, one for each allele of a bi-allelic marker. While ample software is available to obtain genotypes (homozygous for either allele, or heterozygous) for diploid species from these signals, such software is not available for tetraploid species which may be scored as five alternative genotypes (aaaa, baaa, bbaa, bbba and bbbb; nulliplex to quadruplex).</p> <p>Results</p> <p>We present a novel algorithm, implemented in the R package fitTetra, to assign genotypes for bi-allelic markers to tetraploid samples from genotyping assays that produce intensity signals for both alleles. The algorithm is based on the fitting of several mixture models with five components, one for each of the five possible genotypes. The models have different numbers of parameters specifying the relation between the five component means, and some of them impose a constraint on the mixing proportions to conform to Hardy-Weinberg equilibrium (HWE) ratios. The software rejects markers that do not allow a reliable genotyping for the majority of the samples, and it assigns a missing score to samples that cannot be scored into one of the five possible genotypes with sufficient confidence.</p> <p>Conclusions</p> <p>We have validated the software with data of a collection of 224 potato varieties assayed with an Illumina GoldenGate™ 384 SNP array and shown that all SNPs with informative ratio distributions are fitted. Almost all fitted models appear to be correct based on visual inspection and comparison with diploid samples. When the collection of potato varieties is analyzed as if it were a population, almost all markers seem to be in Hardy-Weinberg equilibrium. The R package fitTetra is freely available under the GNU Public License from <url>http://www.plantbreeding.wur.nl/UK/software_fitTetra.html</url> and as Additional files with this article.</p

Gate-tunable black phosphorus spin valve with nanosecond spin lifetimes

Two-dimensional materials offer new opportunities for both fundamental science and technological applications, by exploiting the electron spin. While graphene is very promising for spin communication due to its extraordinary electron mobility, the lack of a band gap restricts its prospects for semiconducting spin devices such as spin diodes and bipolar spin transistors. The recent emergence of 2D semiconductors could help overcome this basic challenge. In this letter we report the first important step towards making 2D semiconductor spin devices. We have fabricated a spin valve based on ultra-thin (5 nm) semiconducting black phosphorus (bP), and established fundamental spin properties of this spin channel material which supports all electrical spin injection, transport, precession and detection up to room temperature (RT). Inserting a few layers of boron nitride between the ferromagnetic electrodes and bP alleviates the notorious conductivity mismatch problem and allows efficient electrical spin injection into an n-type bP. In the non-local spin valve geometry we measure Hanle spin precession and observe spin relaxation times as high as 4 ns, with spin relaxation lengths exceeding 6 um. Our experimental results are in a very good agreement with first-principles calculations and demonstrate that Elliott-Yafet spin relaxation mechanism is dominant. We also demonstrate that spin transport in ultra-thin bP depends strongly on the charge carrier concentration, and can be manipulated by the electric field effect