Deep sequencing-based analysis of gene expression in bovine mammary epithelial cells after Staphylococcus aureus, Escherichia coli, and Klebsiella pneumoniae infection.

The goal of this study was to characterize the transcriptome of primary bovine mammalian epithelial cells (pBMECs) and to identify candidate genes for response and resistance to Staphylococcus aureus (strain S108), Escherichia coli (strain E23), and Klebsiella pneumoniae (strain K96) infection. Using Solexa sequencing, approximately 4.9 million total sequence tags were obtained from each of the three infected libraries and the control library. Gene Ontology (GO) analysis of the S108-infected pBMECs showed differentially expressed genes (DEGs) were significantly involved in metabolic processes. In E23-infected pBMECs, DEGs were predominantly associated with cell death and programmed cell death GO terms, while in K96-infected pBMECs, DEGs were primarily involved in metabolic processes and in utero embryonic development. Analysis of the cluster of orthologous groups of proteins showed that the S108-infected, E23-infected and K96-infected pBMECs were significantly involved in "Translation, ribosomal structure and biogenesis", "General function prediction only" and "Replication, recombination and repair". The transcriptome sequences were also annotated for KEGG orthology, and it was found that DEGs in S108-infected pBMECs were significantly involved in oxidative phosphorylation and Parkinson's disease. The clustered pathway terms of the DEGs of the E23-infected pBMECs were found to involve the NOD-like receptor signaling pathway and oxidative phosphorylation, while those of the K96-infected pBMECs were primarily involved in oxidative phosphorylation and apoptosis. Our results have identified a number of immune-related genes that showed changes in gene expression after bacterial infection, and provided insight into the interactions between pBMECs and the bacteria.

Modelling and simulation of porcine liver tissue indentation using finite element method and uniaxial stress-strain data.

We hypothesize that both compression and elongation stress-strain data should be considered for modeling and simulation of soft tissue indentation. Uniaxial stress-strain data were obtained from in vitro loading experiments of porcine liver tissue. An axisymmetric finite element model was used to simulate liver tissue indentation with tissue material represented by hyperelastic models. The material parameters were derived from uniaxial stress-strain data of compressions, elongations, and combined compression and elongation of porcine liver samples. in vitro indentation tests were used to validate the finite element simulation. Stress-strain data from the simulation with material parameters derived from the combined compression and elongation data match the experimental data best. This is due to its better ability in modeling 3D deformation since the behavior of biological soft tissue under indentation is affected by both its compressive and tensile characteristics. The combined logarithmic and polynomial model is somewhat better than the 5-constant Mooney-Rivlin model as the constitutive model for this indentation simulation.

Liver tissue characterization from uniaxial stress-strain data using probabilistic and inverse finite element methods.

Biological soft tissue is highly inhomogeneous with scattered stress-strain curves. Assuming that the instantaneous strain at a specific stress varies according to a normal distribution, a nondeterministic approach is proposed to model the scattered stress-strain relationship of the tissue samples under compression. Material parameters of the liver tissue modeled using Mooney-Rivlin hyperelastic constitutive equation were represented by a statistical function with normal distribution. Mean and standard deviation of the material parameters were determined using inverse finite element method and inverse mean-value first-order second-moment (IMVFOSM) method respectively. This method was verified using computer simulation based on direct Monte-Carlo (MC) method. The simulated cumulative distribution function (CDF) corresponded well with that of the experimental stress-strain data. The resultant nondeterministic material parameters were able to model the stress-strain curves from other separately conducted liver tissue compression tests. Stress-strain data from these new tests could be predicted using the nondeterministic material parameters.

Accurate two-dimensional cardiac strain calculation using adaptive windowed Fourier transform and Gabor wavelet transform.

PURPOSE: Cardiac strain calculated from tagged magnetic resonance (MR) images provides clinicians information about abnormalities of heart-wall motion in patients. It is important to develop an accurate method to determine the cardiac strain efficiently. An adaptive windowed harmonic phase (AWHARP) method is proposed for cardiac strain calculation. MATERIALS AND METHODS: AWHARP is based on adaptive windowed Fourier transform (AWFT) and 2D Gabor wavelet transform (2D-GWT). The AWFT provides a spatially varying representation of the signal spectra, which allows the harmonic phase (HARP) image to be extracted with high accuracy. Instantaneous spatial frequencies are calculated using 2D-GWT, and the widths of the adaptive windows are then determined according to the instantaneous spatial frequencies for multi-resolution analysis of phase extraction. The proposed method was studied using simulated images and patients' MR images. Both single tagged images (SPAMM) and subtracted tagged images (CSPAMM) were generated using our simulation method, and their results calculated using AWHARP and HARP methods were compared. Normal and pathological tagged MR images were also processed to evaluate the performance of our method. RESULTS: Our experimental results show that the accuracies of phase and strain images calculated using the AWHARP method are higher than that calculated using the HARP method especially for large tag line deformation. The improvement in accuracies can be up to 3.2 strain (E1) and 17.3 calculation from MR images reveals that the cardiac strain in the end-systolic state is significantly reduced for patients with hypertrophic cardiomyopathy (HCM) compared to that of healthy subjects. CONCLUSION: The proposed AWHARP is an accurate and efficient method for cardiac strain estimation from MR images. This new algorithm can help clinicians to detect left ventricle dysfunctions and myocardial diseases with accurate cardiac strain analysis.

Comparison of five DNA extraction methods for molecular analysis of Jerusalem artichoke (Helianthus tuberosus).

DNA extraction is an essential step for molecular analysis of an organism, but it is difficult to acquire a sufficient amount of pure DNA from plant tissue with high levels of phenolic compounds, carbohydrates, proteins, and secondary metabolites. Jerusalem artichoke (Helianthus tuberosus) has high levels of such substances. We compared five commonly used methods of extracting genomic DNA in tests made with leaves and seed of four Jerusalem artichoke genotypes: 1) modified method of Tai and Tanksley, 2) method of Doyle and Doyle, 3) method of Porebski, 4) modified method of Storchová, and 5) Plant DNA Kit of Omega Bio-tek. The quality and quantity of extracted DNAs were assessed by photometric assay, electrophoresis on 1% agarose gel and a PCR-based technique. The modified method of Tai and Tanksley was found to be superior for both young leaves and seed. The quality of the extracted DNA was confirmed by sequence-related amplified polymorphism. This information will be useful for molecular analyses of Jerusalem artichoke and other related Helianthus species.

Effects of imperfections on localized bulging in inflated membrane tubes.

The problem of localized bulging in inflated membrane tubes shares the same features with a variety of other localization problems such as formation of kink bands in fibre-reinforced composites and layered structures. This type of localization is known to be very sensitive to imperfections, but the precise nature of such sensitivity has not so far been quantified. In this paper, we study effects of localized wall thinning/thickening on the onset of localized bulging in inflated membrane tubes as a prototypical example. It is shown that localized wall thinning may reduce the critical pressure or circumferential stretch by an amount of the order of the square root of maximum wall thickness reduction. As a typical example, a 10 per cent maximum wall thinning may reduce the critical circumferential stretch by 19 per cent. This square root law complements the well-known Koiter's two-thirds power law for subcritical periodic bifurcations. The relevance of our results to mathematical modelling of aneurysm formation in human arteries is also discussed.

Genetic diversity of worldwide Jerusalem artichoke (Helianthus tuberosus) germplasm as revealed by RAPD markers.

Jerusalem artichoke (Helianthus tuberosus) is a wild relative of the cultivated sunflower (H. annuus); it is an old tuber crop that has recently received renewed interest. We used RAPD markers to characterize 147 Jerusalem artichoke accessions from nine countries. Thirty RAPD primers were screened; 13 of them detected 357 reproducible RAPD bands, of which 337 were polymorphic. Various diversity analyses revealed several different patterns of RAPD variation. More than 93% of the RAPD variation was found within accessions of a country. Weak genetic differentiation was observed between wild and cultivated accessions. Six groups were detected in this germplasm set. Four ancestral groups were found for the Canadian germplasm. The most genetically distinct accessions were identified. These findings provide useful diversity information for understanding the Jerusalem artichoke gene pool, for conserving Jerusalem artichoke germplasm, and for choosing germplasm for genetic improvement.

Motion tracking and strain map computation for quasi-static magnetic resonance elastography.

This paper presents a new imaging method for quasi-static magnetic resonance elastography (MRE). Tagged magnetic resonance (MR) imaging of human lower leg was acquired with probe indentation using a MR-compatible actuation system. Indentation force was recorded for soft tissue elasticity reconstruction. Motion tracking and strain map of human lower leg are calculated using a harmonic phase (HARP)-based method. Simulated tagged MR images were constructed and analyzed to validate the HARP-based method. Our results show that the proposed imaging method can be used to generate accurate motion distribution and strain maps of the targeted soft tissue.

Point estimation and graphical inference of marginal dominance for two viability loci controlling inbreeding depression.

A deterministic analysis is conducted to examine marginal dominance for two linked viability loci influencing inbreeding depression and its graphical inferences. Four estimators of marginal dominance are derived, assuming a biallelic marker locus completely linked to one of the viability loci, and the biases in expected estimates due to the other deleterious locus are discussed. Three conditions under which apparent partial dominance or underdominance could occur are found, i.e. when two multiplicative, partially recessive loci are linked in coupling phase and when two synergistic, highly overdominant loci are linked in coupling or repulsion phases. Expected frequencies of the three marker genotypes in selfed progeny are derived, considering two linkage phases, two types of marker locus position with respect to the viability loci, and the multiplicative and synergistic fitness models. Segregation ratios are generated for the marker locus linked to either two overdominant or partially recessive loci and plotted in gene action graphs to examine the robustness of the graphical inferences of gene action due to the presence of an additional linked viability locus. Under a multiplicative fitness model, the presence of an additional partially recessive or overdominant locus in the vicinity of the marker locus does not greatly affect the graphical inferences of the relative role of partially recessive or overdominant genes in expression of inbreeding depression. A marker linked to two synergistic, highly overdominant loci can behave as though linked to a partially recessive, partially dominant or underdominant locus, even with relatively weak synergism.

Marker-based inferences about epistasis for genes influencing inbreeding depression.

We describe a multilocus, marker-based regression method for inferring interactions between genes controlling inbreeding depression in self-fertile organisms. It is based upon selfing a parent heterozygous for several unlinked codominant markers, then analyzing the fitness of progeny marker genotypes. If loci causing inbreeding depression are linked to marker loci, then viability selection is manifested by distorted segregation of markers, and fecundity selection by dependence of the fecundity character upon the marker genotype. To characterize this selection, fitness is regressed on the proportion of loci homozygous for markers linked to deleterious alleles, and epistasis is detected by nonlinearity of the regression. Alternatively, fitness can be regressed on the proportion of heterozygous loci. Other modes of selection can be incorporated with a bivariate regression involving both homozygote and heterozygote marker genotypes. The advantage of this marker-based approach is that "purging" is minimized and specific chromosomal segments are identified; its disadvantage lies in low statistical power when linkage is not strong and/or the linkage phase between marker and selected loci is uncertain. Using this method in the wildflower Mimulus guttatus, we found predominant multiplicative gene interaction determining fecundity and some negative synergistic (nonmultiplicative) interaction for viability.

Marker-based inferences about fecundity genes contributing to inbreeding depression in Mimulus guttatus.

Eight unlinked isozyme loci were used as genetic markers to characterize fecundity genes contributing to inbreeding depression in two selfed progeny arrays of Mimulus guttatus. Five fecundity traits were measured. Six of eight marked chromosomal segments were significantly associated with the expression of these traits. The number of genes detected for five traits in two progeny arrays varied, with an average of 2.8 genes per trait. Individual segments explained 1.44-9.29%, and together accounted for 3.85-11.32%, of phenotypic variation. Of 20 significant associations, 10 could be interpreted as exhibiting partial dominance, 7 overdominance, 3 partial recessivity, and 0 underdominance. Significant pairwise epistasis was rare. The results of this study suggest that inbreeding depression is caused by many deleterious genes of relatively small, partially dominant effects.

Evidence for the partial dominance of viability genes contributing to inbreeding depression in Mimulus guttatus.

The relative importance of different modes of gene expression of viability genes contributing to inbreeding depression was investigated in the wild plant, Mimulus guttatus. Viability genes were identified by self-fertilizing 31 outbred plants, each heterozygous for three to nine unlinked allozyme markers, and analyzing segregation ratios of selfed progeny at maturity for deviations from 1:2:1 ratios. In this study, 24 linkages of viability genes to marker loci were detected. To infer the nature of gene action for these viability genes, a "model-free" graphical method was developed that examines the "space" of segregation ratios allowed by each of seven selection models (i.e., overdominance, complete recessivity, partial recessivity, additivity, partial dominance, complete dominance and underdominance). Using this method, we found that, of 24 linkages detected, 18 were consistent with either partial dominance, complete dominance or underdominance. Six were consistent with either partial recessivity, complete recessivity or overdominance. This finding indicates that, in these chromosomal segments identified by allozyme markers, partial dominance plays the predominant role in inbreeding depression. This is inconsistent with either the dominance or overdominance hypotheses proposed to account for inbreeding depression.

On estimating the linkage of marker genes to viability genes controlling inbreeding depression.

Statistical properties and extensions of Hedrick and Muona's method for mapping viability alleles causing inbreeding depression are discussed in this paper. Their method uses the segregation ratios among selfed progeny of a marker-locus heterozygote to estimate the viability reduction, "s", of an allele and its recombination fraction, "c", with the marker. Explicit estimators are derived for c and s, including expressions for their variances. The degree of estimation bias is examined for cases when (1) the viability allele is partially recessive and (2) the marker locus is linked to two viability loci. If linkage or viability reduction is moderate, very large sample sizes are required to obtain reliable estimates of c and s, in part because these estimates show a statistical correlation close to unity. Power is further reduced because alleles causing viability reduction often occur at low frequency at specific loci in a population. To increase power, we present a statistical model for the joint analysis of several selfed progeny arrays selected at random from a population. Assuming a fixed total number of progeny, we determine the optimal number of progeny arrays versus number of progeny per array under this model. We also examine the increase of information provided by a second, flanking marker. Two flanking markers provide vastly superior estimation properties, reducing sample sizes by approximately 95% from those required by a single marker.

Pollen pool heterogeneity in jack pine (Pinus banksiana Lamb.): a problem for estimating outcrossing rates?

Pollen pool heterogeneity, which violates an assumption of the mixed-mating model, is a major potential problem in measuring plant mating systems. In this study, isozyme markers were used to examine pollen pool heterogeneity in two natural populations of jack pine, Pinus banksiana Lamb., in northwestern Ontario, Canada. Population multilocus estimates of outcrossing rate ranged from 0.83 to 0.95 and differed significantly between populations. Single-tree multilocus outcrossing rates were found to be homogeneous among trees in both populations. Computer simulation studies indicated that a consanguineous pollen pool (pollen gametes related to the mother tree) was capable of biasing population outcrossing estimates downward. Random pollen pool heterogeneity (uncorrelated with maternal genotypes) did not appear to affect population outcrossing estimates in the simulations. Heterogeneity G-tests and Spearman rank tests showed that pollen pool heterogeneity existed in the two natural populations examined; however, it did not have a major effect on population outcrossing estimates, since the consanguineous pollen pool detected was probably a relatively minor component of the outcross pollen pool in both populations. In addition, heterogeneity G-tests were found to be not sensitive in detecting pollen pool heterogeneity caused by consanguineous pollen pool.