Illuminating Choices for Library Prep: A Comparison of Library Preparation Methods for Whole Genome Sequencing of Cryptococcus neoformans Using Illumina HiSeq.

The industry of next-generation sequencing is constantly evolving, with novel library preparation methods and new sequencing machines being released by the major sequencing technology companies annually. The Illumina TruSeq v2 library preparation method was the most widely used kit and the market leader; however, it has now been discontinued, and in 2013 was replaced by the TruSeq Nano and TruSeq PCR-free methods, leaving a gap in knowledge regarding which is the most appropriate library preparation method to use. Here, we used isolates from the pathogenic fungi Cryptococcus neoformans var. grubii and sequenced them using the existing TruSeq DNA v2 kit (Illumina), along with two new kits: the TruSeq Nano DNA kit (Illumina) and the NEBNext Ultra DNA kit (New England Biolabs) to provide a comparison. Compared to the original TruSeq DNA v2 kit, both newer kits gave equivalent or better sequencing data, with increased coverage. When comparing the two newer kits, we found little difference in cost and workflow, with the NEBNext Ultra both slightly cheaper and faster than the TruSeq Nano. However, the quality of data generated using the TruSeq Nano DNA kit was superior due to higher coverage at regions of low GC content, and more SNPs identified. Researchers should therefore evaluate their resources and the type of application (and hence data quality) being considered when ultimately deciding on which library prep method to use

How Many Subpopulations is Too Many? Exponential Lower Bounds for Inferring Population Histories

Reconstruction of population histories is a central problem in population genetics. Existing coalescent-based methods, like the seminal work of Li and Durbin (Nature, 2011), attempt to solve this problem using sequence data but have no rigorous guarantees. Determining the amount of data needed to correctly reconstruct population histories is a major challenge. Using a variety of tools from information theory, the theory of extremal polynomials, and approximation theory, we prove new sharp information-theoretic lower bounds on the problem of reconstructing population structure -- the history of multiple subpopulations that merge, split and change sizes over time. Our lower bounds are exponential in the number of subpopulations, even when reconstructing recent histories. We demonstrate the sharpness of our lower bounds by providing algorithms for distinguishing and learning population histories with matching dependence on the number of subpopulations. Along the way and of independent interest, we essentially determine the optimal number of samples needed to learn an exponential mixture distribution information-theoretically, proving the upper bound by analyzing natural (and efficient) algorithms for this problem.Comment: 38 pages, Appeared in RECOMB 201

Treatment of secondary hyperparathyroidism in haemodialysis patients: a randomised clinical trial comparing paricalcitol and alfacalcidol

<p>Abstract</p> <p>Background</p> <p>Secondary hyperparathyroidism is a common feature in patients with chronic kidney disease. Its serious clinical consequences include renal osteodystrophy, calcific uremic arteriolopathy, and vascular calcifications that increase morbidity and mortality.</p> <p>Reduced synthesis of active vitamin D contributes to secondary hyperparathyroidism. Therefore, this condition is managed with activated vitamin D. However, hypercalcemia and hyperphosphatemia limit the use of activated vitamin D.</p> <p>In Denmark alfacalcidol is the primary choice of vitamin D analog.</p> <p>A new vitamin D analog, paricalcitol, may be less prone to induce hypercalcemia and hyperphosphatemia.</p> <p>However, a randomised controlled clinical study comparing alfacalcidol and paricalcitol has never been performed.</p> <p>The primary objective of this study is to compare alfacalcidol and paricalcitol. We evaluate the suppression of the secondary hyperparathyroidism and the tendency towards hyperphosphatemia and hypercalcemia.</p> <p>Methods/Design</p> <p>This is an investigator-initiated cross-over study. Nine Danish haemodialysis units will recruit 117 patients with end stage renal failure on maintenance haemodialysis therapy.</p> <p>Patients are randomised into two treatment arms. After a wash out period of 6 weeks they receive increasing doses of alfacalcidol or paricalcitol for a period of 16 weeks and after a further wash out period of 6 weeks they receive the contrary treatment (paricalcitol or alfacalcidol) for 16 weeks.</p> <p>Discussion</p> <p>Hyperparathyroidism, hypercalcemia and hyperphosphatemia are associated with increased cardiovascular mortality in patients with chronic kidney disease.</p> <p>If there is any difference in the ability of these two vitamin D analogs to decrease the secondary hyperparathyroidism without causing hypercalcemia and hyperphosphatemia, there may also be a difference in the risk of cardiovascular mortality depending on which vitamin D analog that are used. This has potential major importance for this group of patients.</p> <p>Trial registration</p> <p>ClinicalTrials.gov NCT004695</p

Wide-Field Motion Integration in Fly VS Cells: Insights from an Inverse Approach

Fly lobula plate tangential cells are known to perform wide-field motion integration. It is assumed that the shape of these neurons, and in particular the shape of the subclass of VS cells, is responsible for this type of computation. We employed an inverse approach to investigate the morphology-function relationship underlying wide-field motion integration in VS cells. In the inverse approach detailed, model neurons are optimized to perform a predefined computation: here, wide-field motion integration. We embedded the model neurons to be optimized in a biologically plausible model of fly motion detection to provide realistic inputs, and subsequently optimized model neuron with and without active conductances (gNa, gK, gK(Na)) along their dendrites to perform this computation. We found that both passive and active optimized model neurons perform well as wide-field motion integrators. In addition, all optimized morphologies share the same blueprint as real VS cells. In addition, we also found a recurring blueprint for the distribution of gK and gNa in the active models. Moreover, we demonstrate how this morphology and distribution of conductances contribute to wide-field motion integration. As such, by using the inverse approach we can predict the still unknown distribution of gK and gNa and their role in motion integration in VS cells

Recombination rate and selection strength in HIV intra-patient evolution

The evolutionary dynamics of HIV during the chronic phase of infection is driven by the host immune response and by selective pressures exerted through drug treatment. To understand and model the evolution of HIV quantitatively, the parameters governing genetic diversification and the strength of selection need to be known. While mutation rates can be measured in single replication cycles, the relevant effective recombination rate depends on the probability of coinfection of a cell with more than one virus and can only be inferred from population data. However, most population genetic estimators for recombination rates assume absence of selection and are hence of limited applicability to HIV, since positive and purifying selection are important in HIV evolution. Here, we estimate the rate of recombination and the distribution of selection coefficients from time-resolved sequence data tracking the evolution of HIV within single patients. By examining temporal changes in the genetic composition of the population, we estimate the effective recombination to be r=1.4e-5 recombinations per site and generation. Furthermore, we provide evidence that selection coefficients of at least 15% of the observed non-synonymous polymorphisms exceed 0.8% per generation. These results provide a basis for a more detailed understanding of the evolution of HIV. A particularly interesting case is evolution in response to drug treatment, where recombination can facilitate the rapid acquisition of multiple resistance mutations. With the methods developed here, more precise and more detailed studies will be possible, as soon as data with higher time resolution and greater sample sizes is available.Comment: to appear in PLoS Computational Biolog

Dependence of paracentric inversion rate on tract length

BACKGROUND: We develop a Bayesian method based on MCMC for estimating the relative rates of pericentric and paracentric inversions from marker data from two species. The method also allows estimation of the distribution of inversion tract lengths. RESULTS: We apply the method to data from Drosophila melanogaster and D. yakuba. We find that pericentric inversions occur at a much lower rate compared to paracentric inversions. The average paracentric inversion tract length is approx. 4.8 Mb with small inversions being more frequent than large inversions. If the two breakpoints defining a paracentric inversion tract are uniformly and independently distributed over chromosome arms there will be more short tract-length inversions than long; we find an even greater preponderance of short tract lengths than this would predict. Thus there appears to be a correlation between the positions of breakpoints which favors shorter tract lengths. CONCLUSION: The method developed in this paper provides the first statistical estimator for estimating the distribution of inversion tract lengths from marker data. Application of this method for a number of data sets may help elucidate the relationship between the length of an inversion and the chance that it will get accepted

Delineation of Culicoides species by morphology and barcode exemplified by three new species of the subgenus Culicoides (Diptera: Ceratopogonidae) from Scandinavia

BACKGROUND: Culicoides biting midges (Diptera: Ceratopogonidae) cause biting nuisance to livestock and humans and are vectors of a range of pathogens of medical and veterinary importance. Despite their economic significance, the delineation and identification of species where only morphology is considered, as well as the evolutionary relationships between species within this genus remains problematic. In recent years molecular barcoding has assisted substantially in the identification of biting midges in the multiple entomological survey projects which were initiated in many European countries following the bluetongue outbreak in 2006–2009. These studies revealed potentially new species and “species-complexes” with large genetic and morphological variability. Here we use molecular barcoding, together with morphological analysis, to study subgenus Culicoides Latreille from Scandinavia with focus on three potentially new species. METHODS: Biting midges were collected at various sites in Denmark and Sweden. Culicoides specimens were described by variation of a fragment of their cytochrome c oxidase subunit 1 (COI) gene sequence and wing, palp and antennal characters. RESULTS: It is shown that three new species initially separated by DNA barcoding with mitochondrial COI can be distinguished by morphological characters. In this context a key to Scandinavian subgenus Culicoides using wing and maxillary palp characters is presented. The key is including the three new species Culicoides boyi, Culicoides selandicus and Culicoides kalix. CONCLUSION: Three new species of Culicoides biting midges were identified and could be identified by both molecular and morphological differences. Evaluation of differences between and within taxa of biting midges using COI barcode yielded a rough estimate of species delineation; interspecies differences across Culicoides subgenera approaches 20%, whereas intraspecies differences are below 4% and in most cases below 1%. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s13071-015-0750-4) contains supplementary material, which is available to authorized users

Extracellular Hsp72 concentration relates to a minimum endogenous criteria during acute exercise-heat exposure

Extracellular heat-shock protein 72 (eHsp72) concentration increases during exercise-heat stress when conditions elicit physiological strain. Differences in severity of environmental and exercise stimuli have elicited varied response to stress. The present study aimed to quantify the extent of increased eHsp72 with increased exogenous heat stress, and determine related endogenous markers of strain in an exercise-heat model. Ten males cycled for 90 min at 50% O2peak in three conditions (TEMP, 20°C/63% RH; HOT, 30.2°C/51%RH; VHOT, 40.0°C/37%RH). Plasma was analysed for eHsp72 pre, immediately post and 24-h post each trial utilising a commercially available ELISA. Increased eHsp72 concentration was observed post VHOT trial (+172.4%) (P<0.05), but not TEMP (-1.9%) or HOT (+25.7%) conditions. eHsp72 returned to baseline values within 24hrs in all conditions. Changes were observed in rectal temperature (Trec), rate of Trec increase, area under the curve for Trec of 38.5°C and 39.0°C, duration Trec ≥ 38.5°C and ≥ 39.0°C, and change in muscle temperature, between VHOT, and TEMP and HOT, but not between TEMP and HOT. Each condition also elicited significantly increasing physiological strain, described by sweat rate, heart rate, physiological strain index, rating of perceived exertion and thermal sensation. Stepwise multiple regression reported rate of Trec increase and change in Trec to be predictors of increased eHsp72 concentration. Data suggests eHsp72 concentration increases once systemic temperature and sympathetic activity exceeds a minimum endogenous criteria elicited during VHOT conditions and is likely to be modulated by large, rapid changes in core temperature