A validated model of passive skeletal muscle to predict force and intramuscular pressure

The passive properties of skeletal muscle are often overlooked in muscle studies, yet they play a key role in tissue function in vivo. Studies analyzing and modeling muscle passive properties, while not uncommon, have never investigated the role of fluid content within the tissue. Additionally, intramuscular pressure (IMP) has been shown to correlate with muscle force in vivo and could be used to predict muscle force in the clinic. In this study, a novel model of skeletal muscle was developed and validated to predict both muscle stress and IMP under passive conditions for the New Zealand White Rabbit tibialis anterior. This model is the first to include fluid content within the tissue and uses whole muscle geometry. A nonlinear optimization scheme was highly effective at fitting model stress output to experimental stress data (normalized mean square error or NMSE fit value of 0.993) and validation showed very good agreement to experimental data (NMSE fit values of 0.955 and 0.860 for IMP and stress, respectively). While future work to include muscle activation would broaden the physiological application of this model, the passive implementation could be used to guide surgeries where passive muscle is stretched.     

Reliably Detecting Clinically Important Variants Requires Both Combined Variant Calls and Optimized Filtering Strategies

A diversity of tools is available for identification of variants from genome sequence data. Given the current complexity of incorporating external software into a genome analysis infrastructure, a tendency exists to rely on the results from a single tool alone. The quality of the output variant calls is highly variable however, depending on factors such as sequence library quality as well as the choice of short-read aligner, variant caller, and variant caller filtering strategy. Here we present a two-part study first using the high quality ‘genome in a bottle’ reference set to demonstrate the significant impact the choice of aligner, variant caller, and variant caller filtering strategy has on overall variant call quality and further how certain variant callers outperform others with increased sample contamination, an important consideration when analyzing sequenced cancer samples. This analysis confirms previous work showing that combining variant calls of multiple tools results in the best quality resultant variant set, for either specificity or sensitivity, depending on whether the intersection or union, of all variant calls is used respectively. Second, we analyze a melanoma cell line derived from a control lymphocyte sample to determine whether software choices affect the detection of clinically important melanoma risk-factor variants finding that only one of the three such variants is unanimously detected under all conditions. Finally, we describe a cogent strategy for implementing a clinical variant detection pipeline; a strategy that requires careful software selection, variant caller filtering optimizing, and combined variant calls in order to effectively minimize false negative variants. While implementing such features represents an increase in complexity and computation the results offer indisputable improvements in data quality.     

The interaction of decorin core protein fragments with type I collagen.

To further define the molecular interaction between decorin and type I collagen we generated a 20 kD fragment containing the N-terminal half of the core protein by Endoproteinase Arg C digestion and a 40 kD fragment including all leucine-rich repeats in the central part of decorin core by cleavage with 2-nitro-5-thiocyanobenzoate. The fragments did not influence collagen fibril formation, even at high concentration, and radioactive fragments showed little binding to collagen fibrils. Our observations suggest that neither the N-terminal half nor the central leucine-rich repeats of the decorin core protein can, by itself, interact fully with fibrillar collagen.     

How reliable are the Powell–Wetherall plot method and the maximum-length approach? Implications for length-based studies of growth and mortality

Length-based methods are the cornerstone of many population studies and stock assessments. This study tested two widely used methods: the Powell–Wetherall (P–W) plot and the L_max approach (i.e., estimating L_∞ directly from L_max). In most simulations, P–W estimates of the ratio total mortality/growth (Z/K ratio) were biased beyond acceptable limits (bias?>?30%). Bias in Z/K showed a complex behavior, without possible corrections. Estimates of asymptotic length (L_∞) were less biased than Z/K, but were very sensitive to intra-cohort variability in growth and to changes in the occurrence of large individuals in the sample. Exclusion of the largest size classes during the regression procedure or weighing by abundance does not solve these issues. Perfect linearization of the data and extremely narrow confidence intervals for Z/K will lead users to erroneous overconfidence in outputs. Clearly, the P–W method is not suitable for the assessment of Z/K ratios of natural populations. Estimation of L_∞ may be tentatively possible under very specific conditions, with necessary external verifications. Also, this study demonstrates that there is no way to estimate L_∞ directly from L_max, since there is no particular relationship to expect a priori between L_∞ and L_max. Errors in estimating L_∞ directly affect the estimate of the growth constant K and all other subsequent calculations in population studies, stock assessments and ecosystem models. New approaches are urgently needed for length-based studies of body growth (e.g., unconstrained curve fit with subsequent bootstrapping), that consider the inherent uncertainty regarding the underlying data and processes.     

Capturing pair-wise epistatic effects associated with three agronomic traits in barley

Genetic association mapping has been widely applied to determine genetic markers favorably associated with a trait of interest and provide information for marker-assisted selection. Many association mapping studies commonly focus on main effects due to intolerable computing intensity. This study aims to select several sets of DNA markers with potential epistasis to maximize genetic variations of some key agronomic traits in barley. By doing so, we integrated a MDR (multifactor dimensionality reduction) method with a forward variable selection approach. This integrated approach was used to determine single nucleotide polymorphism pairs with epistasis effects associated with three agronomic traits: heading date, plant height, and grain yield in barley from the barley Coordinated Agricultural Project. Our results showed that four, seven, and five SNP pairs accounted for 51.06, 45.66 and 40.42% for heading date, plant height, and grain yield, respectively with epistasis being considered, while corresponding contributions to these three traits were 45.32, 31.39, 31.31%, respectively without epistasis being included. The results suggested that epistasis model was more effective than non-epistasis model in this study and can be more preferred for other applications.     

Exotic grass invasion impacts fitness of an endangered prairie butterfly, <Emphasis Type="Italic">Icaricia icarioides fenderi</Emphasis>

Fender’s blue butterfly is an endangered species restricted to fragmented, grassland remnants that are becoming increasingly dominated by tall, invasive grasses in western Oregon, USA. I performed a removal experiment to assess the impacts of structural degradation accompanying the invasion of Arrhenatherum elatius, tall oat grass, on butterfly fitness and fitness related behaviors. Clipping of A. elatius to native grass sward height resulted in 2.5–5 times as many eggs laid per leaf of host plant. Both male and female butterflies basked more frequently in areas removed of A. elatius inflorescences and upon encountering the treatment edge butterflies had a high rate of return into a large area removed of the grass inflorescences. Although butterfly behavior appeared to be affected by the change in sward height on the treatment edge, there was no evidence for the edge causing a disproportionate egg load. Invasion and dominance by A. elatius appeared to diminish host plant apparency which may result in overloading of eggs on conspicuous host plants, increased incidence of emigration, and a decrease in the likelihood of colonization because female butterflies appeared indifferent to larval resources beneath A. elatius inflorescences. Dominance of natural shortgrass prairies by tall stature grasses like A. elatius may be an insidious form of habitat degradation for grassland Lepidoptera worldwide, but it may go largely unnoticed because larval and adult resources can persist under the unnaturally tall grass canopy.     

Multi-model recursive identification for nonlinear systems with non-uniformly sampling

A recursive least squares based on Multi-model is proposed for non-uniformly sampled-data nonlinear (NUSDN) systems. The corresponding state space model of an NUSDN system is derived using lifting technique. Taking advantage of the Fuzzy c-Mean Clustering algorithm, NUSDN is divided into several local models. The basic idea is that the NUSDN system is viewed as a model switching system under a given rule. Once the local models are identified, the global model is determined. A pH neutralization process validate the performance of the proposed algorithm.