Strobe sequence design for haplotype assembly

Abstract Background Humans are diploid, carrying two copies of each chromosome, one from each parent. Separating the paternal and maternal chromosomes is an important component of genetic analyses such as determining genetic association, inferring evolutionary scenarios, computing recombination rates, and detecting cis-regulatory events. As the pair of chromosomes are mostly identical to each other, linking together of alleles at heterozygous sites is sufficient to phase, or separate the two chromosomes. In Haplotype Assembly, the linking is done by sequenced fragments that overlap two heterozygous sites. While there has been a lot of research on correcting errors to achieve accurate haplotypes via assembly, relatively little work has been done on designing sequencing experiments to get long haplotypes. Here, we describe the different design parameters that can be adjusted with next generation and upcoming sequencing technologies, and study the impact of design choice on the length of the haplotype. Results We show that a number of parameters influence haplotype length, with the most significant one being the advance length (distance between two fragments of a clone). Given technologies like strobe sequencing that allow for large variations in advance lengths, we design and implement a simulated annealing algorithm to sample a large space of distributions over advance-lengths. Extensive simulations on individual genomic sequences suggest that a non-trivial distribution over advance lengths results a 1-2 order of magnitude improvement in median haplotype length. Conclusions Our results suggest that haplotyping of large, biologically important genomic regions is feasible with current technologies

Epidemics on contact networks: a general stochastic approach

Dynamics on networks is considered from the perspective of Markov stochastic processes. We partially describe the state of the system through network motifs and infer any missing data using the available information. This versatile approach is especially well adapted for modelling spreading processes and/or population dynamics. In particular, the generality of our systematic framework and the fact that its assumptions are explicitly stated suggests that it could be used as a common ground for comparing existing epidemics models too complex for direct comparison, such as agent-based computer simulations. We provide many examples for the special cases of susceptible-infectious-susceptible (SIS) and susceptible-infectious-removed (SIR) dynamics (e.g., epidemics propagation) and we observe multiple situations where accurate results may be obtained at low computational cost. Our perspective reveals a subtle balance between the complex requirements of a realistic model and its basic assumptions.Comment: Main document: 16 pages, 7 figures. Electronic Supplementary Material (included): 6 pages, 1 tabl

The impact of contact tracing in clustered populations

The tracing of potentially infectious contacts has become an important part of the control strategy for many infectious diseases, from early cases of novel infections to endemic sexually transmitted infections. Here, we make use of mathematical models to consider the case of partner notification for sexually transmitted infection, however these models are sufficiently simple to allow more general conclusions to be drawn. We show that, when contact network structure is considered in addition to contact tracing, standard “mass action” models are generally inadequate. To consider the impact of mutual contacts (specifically clustering) we develop an improvement to existing pairwise network models, which we use to demonstrate that ceteris paribus, clustering improves the efficacy of contact tracing for a large region of parameter space. This result is sometimes reversed, however, for the case of highly effective contact tracing. We also develop stochastic simulations for comparison, using simple re-wiring methods that allow the generation of appropriate comparator networks. In this way we contribute to the general theory of network-based interventions against infectious disease

A rare case of arterial avulsion presenting with occult blood loss following total hip arthroplasty: a case report

INTRODUCTION: Iatrogenic arterial damage during total hip replacement is a rare but potentially life- or limb-threatening complication. To the best of our knowledge, this is the first reported case of an avulsion injury to a posterior branch of the profunda femoral artery during primary hip arthroplasty. CASE PRESENTATION: We describe the case of a 55-year-old Caucasian man who underwent a total hip replacement. The patient's hemoglobin levels dropped postoperatively, but there was no obvious bleeding, hemodynamic instability, pulsatile mass, or limb ischemia. The patient's hemoglobin levels continued to drop despite nine units of transfused blood. Three days after surgery, the patient underwent an angiography that showed an avulsion injury to a posterior branch of the profunda femoral artery. The avulsion was ligated and the hematoma was evacuated. CONCLUSION: Vascular damage may present in many ways including obvious bleeding, haemodynamic instability, a pulsatile mass, limb ischemia, and occult blood loss. Any of these signs in isolation or in combination could represent a vascular injury and an urgent angiogram should be considered

Remote sensing imagery segmentation: A hybrid approach

In remote sensing imagery, segmentation techniques fail to encounter multiple regions of interest due to challenges such as dense features, low illumination, uncertainties, and noise. Consequently, exploiting vast and redundant information makes segmentation a difficult task. Existing multilevel thresholding techniques achieve low segmentation accuracy with high temporal difficulty due to the absence of spatial information. To mitigate this issue, this paper presents a new Rényi’s entropy and modified cuckoo search-based robust automatic multi-thresholding algorithm for remote sensing image analysis. In the proposed method, the modified cuckoo search algorithm is combined with Rényi’s entropy thresholding criteria to determine optimal thresholds. In the modified cuckoo search algorithm, the Lévy flight step size was modified to improve the convergence rate. An experimental analysis was conducted to validate the proposed method, both qualitatively and quantitatively against existing metaheuristic-based thresholding methods. To do this, the performance of the proposed method was intensively examined on high-dimensional remote sensing imageries. Moreover, numerical parameter analysis is presented to compare the segmented results against the gray-level co-occurrence matrix, Otsu energy curve, minimum cross entropy, and Rényi’s entropy-based thresholding. Experiments demonstrated that the proposed approach is effective and successful in attaining accurate segmentation with low time complexity

Methods for adjusting population structure and familial relatedness in association test for collective effect of multiple rare variants on quantitative traits

Because of the low frequency of rare genetic variants in observed data, the statistical power of detecting their associations with target traits is usually low. The collapsing test of collective effect of multiple rare variants is an important and useful strategy to increase the power; in addition, family data may be enriched with causal rare variants and therefore provide extra power. However, when family data are used, both population structure and familial relatedness need to be adjusted for the possible inflation of false positives. Using a unified mixed linear model and family data, we compared six methods to detect the association between multiple rare variants and quantitative traits. Through the analysis of 200 replications of the quantitative trait Q2 from the Genetic Analysis Workshop 17 data set simulated for 697 subjects from 8 extended families, and based on quantile-quantile plots under the null and receiver operating characteristic curves, we compared the false-positive rate and power of these methods. We observed that adjusting for pedigree-based kinship gives the best control for false-positive rate, whereas adjusting for marker-based identity by state slightly outperforms in terms of power. An adjustment based on a principal components analysis slightly improves the false-positive rate and power. Taking into account type-1 error, power, and computational efficiency, we find that adjusting for pedigree-based kinship seems to be a good choice for the collective test of association between multiple rare variants and quantitative traits using family data

A measure of individual role in collective dynamics

Identifying key players in collective dynamics remains a challenge in several research fields, from the efficient dissemination of ideas to drug target discovery in biomedical problems. The difficulty lies at several levels: how to single out the role of individual elements in such intermingled systems, or which is the best way to quantify their importance. Centrality measures describe a node's importance by its position in a network. The key issue obviated is that the contribution of a node to the collective behavior is not uniquely determined by the structure of the system but it is a result of the interplay between dynamics and network structure. We show that dynamical influence measures explicitly how strongly a node's dynamical state affects collective behavior. For critical spreading, dynamical influence targets nodes according to their spreading capabilities. For diffusive processes it quantifies how efficiently real systems may be controlled by manipulating a single node.Comment: accepted for publication in Scientific Report

Comparison of contact patterns relevant for transmission of respiratory pathogens in Thailand and the Netherlands using respondent-driven sampling

Understanding infection dynamics of respiratory diseases requires the identification and quantification of behavioural, social and environmental factors that permit the transmission of these infections between humans. Little empirical information is available about contact patterns within real-world social networks, let alone on differences in these contact networks between populations that differ considerably on a socio-cultural level. Here we compared contact network data that were collected in the Netherlands and Thailand using a similar online respondent-driven method. By asking participants to recruit contact persons we studied network links relevant for the transmission of respiratory infections. We studied correlations between recruiter and recruited contacts to investigate mixing patterns in the observed social network components. In both countries, mixing patterns were assortative by demographic variables and random by total numbers of contacts. However, in Thailand participants reported overall more contacts which resulted in higher effective contact rates. Our findings provide new insights on numbers of contacts and mixing patterns in two different populations. These data could be used to improve parameterisation of mathematical models used to design control strategies. Although the spread of infections through populations depends on more factors, found similarities suggest that spread may be similar in the Netherlands and Thailand

Symptomatic snapping knee from biceps femoris tendon subluxation: an unusual case of lateral pain in a marathon runner

Snapping biceps femoris syndrome is an uncommon cause of lateral knee pain and may be difficult to diagnose, resulting in unsuccessful surgical intervention. In this report, we present an unusual case of a 37-year-old male marathon runner with unilateral snapping knee secondary to dislocation of the long head of the biceps femoris over the fibular head during knee flexion. The pain was great enough to interfere with his ability to practice sport. Possible causes of symptomatic snapping knee include multiple intra-articular or extra-articular pathology. Biceps femoris snapping over the fibular head is a rare condition. Reported causes include an anomalous insertion of the tendon into the tibia, trauma, and fibular-head abnormality. However, none of those conditions accounted for his symptoms. Failing conservative treatment, the patient underwent surgery for partial resection of the fibular head, with subsequent sudden resolution of symptoms and return to sport. Accurate knowledge and management of this rare condition is mandatory to avoid inappropriate therapy and unnecessary surgical procedures

Dynamics of multi-stage infections on networks

This paper investigates the dynamics of infectious diseases with a nonexponentially distributed infectious period. This is achieved by considering a multistage infection model on networks. Using pairwise approximation with a standard closure, a number of important characteristics of disease dynamics are derived analytically, including the final size of an epidemic and a threshold for epidemic outbreaks, and it is shown how these quantities depend on disease characteristics, as well as the number of disease stages. Stochastic simulations of dynamics on networks are performed and compared to output of pairwise models for several realistic examples of infectious diseases to illustrate the role played by the number of stages in the disease dynamics. These results show that a higher number of disease stages results in faster epidemic outbreaks with a higher peak prevalence and a larger final size of the epidemic. The agreement between the pairwise and simulation models is excellent in the cases we consider