Smart histogram analysis applied to the skull-stripping problem in T1-weighted MRI

In this paper we address the "skull-stripping" problem in 3D MR images. We propose a new method that employs an efficient and unique histogram analysis. A fundamental component of this analysis is an algorithm for partitioning a histogram based on the position of the maximum deviation from a Gaussian fit. In our experiments we use a comprehensive image database, including both synthetic and real MRI. and compare our method with other two well-known methods, namely BSE and BET. For all datasets we achieved superior results. Our method is also highly independent of parameter tuning and very robust across considerable variations of noise ratio.

Analysis of plant LTR-retrotransposons at the fine-scale family level reveals individual molecular patterns

Background: Sugarcane is an important crop worldwide for sugar production and increasingly, as a renewable energy source. Modern cultivars have polyploid, large complex genomes, with highly unequal contributions from ancestral genomes. Long Terminal Repeat retrotransposons (LTR-RTs) are the single largest components of most plant genomes and can substantially impact the genome in many ways. It is therefore crucial to understand their contribution to the genome and transcriptome, however a detailed study of LTR-RTs in sugarcane has not been previously carried out. Results: Sixty complete LTR-RT elements were classified into 35 families within four Copia and three Gypsy lineages. Structurally, within lineages elements were similar, between lineages there were large size differences. FISH analysis resulted in the expected pattern of Gypsy/heterochromatin, Copia/euchromatin, but in two lineages there was localized clustering on some chromosomes. Analysis of related ESTs and RT-PCR showed transcriptional variation between tissues and families. Four distinct patterns were observed in sRNA mapping, the most unusual of which was that of Ale1, with very large numbers of 24nt sRNAs in the coding region. The results presented support the conclusion that distinct small RNA-regulated pathways in sugarcane target the lineages of LTR-RT elements. Conclusions: Individual LTR-RT sugarcane families have distinct structures, and transcriptional and regulatory signatures. Our results indicate that in sugarcane individual LTR-RT families have distinct behaviors and can potentially impact the genome in diverse ways. For instance, these transposable elements may affect nearby genes by generating a diverse set of small RNA's that trigger gene silencing mechanisms. There is also some evidence that ancestral genomes contribute significantly different element numbers from particular LTR-RT lineages to the modern sugarcane cultivar genome.

Automatic design of decision-tree induction algorithms tailored to flexible-receptor docking data

Background: This paper addresses the prediction of the free energy of binding of a drug candidate with enzyme InhA associated with Mycobacterium tuberculosis. This problem is found within rational drug design, where interactions between drug candidates and target proteins are verified through molecular docking simulations. In this application, it is important not only to correctly predict the free energy of binding, but also to provide a comprehensible model that could be validated by a domain specialist. Decision-tree induction algorithms have been successfully used in drug-design related applications, specially considering that decision trees are simple to understand, interpret, and validate. There are several decision-tree induction algorithms available for general-use, but each one has a bias that makes it more suitable for a particular data distribution. In this article, we propose and investigate the automatic design of decision-tree induction algorithms tailored to particular drug-enzyme binding data sets. We investigate the performance of our new method for evaluating binding conformations of different drug candidates to InhA, and we analyze our findings with respect to decision tree accuracy, comprehensibility, and biological relevance. Results: The empirical analysis indicates that our method is capable of automatically generating decision-tree induction algorithms that significantly outperform the traditional C4.5 algorithm with respect to both accuracy and comprehensibility. In addition, we provide the biological interpretation of the rules generated by our approach, reinforcing the importance of comprehensible predictive models in this particular bioinformatics application. Conclusions: We conclude that automatically designing a decision-tree algorithm tailored to molecular docking data is a promising alternative for the prediction of the free energy from the binding of a drug candidate with a flexible-receptor.

Different morphology, stage and treatment affect immune cell infiltration and long-term outcome in patients with non-small-cell lung carcinoma

Aims: Development of effective immune-based therapies for patients with non-small-cell lung carcinoma (NSCLC) depends on an accurate characterization of complex interactions that occur between immune cells and the tumour environment. Methods and results: Innate and adaptive immune responses were evaluated in relation to prognosis in 65 patients with surgically excised NSCLC. Immunohistochemistry and morphometry were used to determine the abundance and distribution of immune cells. We found low numbers of immune cells and levels of cytokines in the tumour environment when compared with surrounding parenchyma. Smoking was associated inversely with the adaptive immune response and directly with innate immunity. We observed a prominent adaptive immune response in squamous cell carcinomas (SCC) but greater innate immune responses in adenocarcinomas and large cell carcinomas. Cox model analysis showed a low risk of death for smoking <41 packs/year, N-0 tambour stage, squamous carcinoma, CD4(+) > 16.81% and macrophages/monocytes >4.5%. Collectively, the data indicate that in NSCLC there is not a substantive local immune cell infiltrate within the tumour. Conclusion: Although immune cell infiltration is limited in NSCLC it appears to have an impact on prognosis and this may be of relevance for new immunotherapeutic approaches.

Automatic aspect discrimination in data clustering

The attributes describing a data set may often be arranged in meaningful subsets, each of which corresponds to a different aspect of the data. An unsupervised algorithm (SCAD) that simultaneously performs fuzzy clustering and aspects weighting was proposed in the literature. However, SCAD may fail and halt given certain conditions. To fix this problem, its steps are modified and then reordered to reduce the number of parameters required to be set by the user. In this paper we prove that each step of the resulting algorithm, named ASCAD, globally minimizes its cost-function with respect to the argument being optimized. The asymptotic analysis of ASCAD leads to a time complexity which is the same as that of fuzzy c-means. A hard version of the algorithm and a novel validity criterion that considers aspect weights in order to estimate the number of clusters are also described. The proposed method is assessed over several artificial and real data sets.

Comparative Genome Analysis of Trichophyton rubrum and Related Dermatophytes Reveals Candidate Genes Involved in Infection

The major cause of athlete's foot is Trichophyton rubrum, a dermatophyte or fungal pathogen of human skin. To facilitate molecular analyses of the dermatophytes, we sequenced T. rubrum and four related species, Trichophyton tonsurans, Trichophyton equinum, Microsporum canis, and Microsporum gypseum. These species differ in host range, mating, and disease progression. The dermatophyte genomes are highly colinear yet contain gene family expansions not found in other human-associated fungi. Dermatophyte genomes are enriched for gene families containing the LysM domain, which binds chitin and potentially related carbohydrates. These LysM domains differ in sequence from those in other species in regions of the peptide that could affect substrate binding. The dermatophytes also encode novel sets of fungus-specific kinases with unknown specificity, including nonfunctional pseudokinases, which may inhibit phosphorylation by competing for kinase sites within substrates, acting as allosteric effectors, or acting as scaffolds for signaling. The dermatophytes are also enriched for a large number of enzymes that synthesize secondary metabolites, including dermatophyte-specific genes that could synthesize novel compounds. Finally, dermatophytes are enriched in several classes of proteases that are necessary for fungal growth and nutrient acquisition on keratinized tissues. Despite differences in mating ability, genes involved in mating and meiosis are conserved across species, suggesting the possibility of cryptic mating in species where it has not been previously detected. These genome analyses identify gene families that are important to our understanding of how dermatophytes cause chronic infections, how they interact with epithelial cells, and how they respond to the host immune response. IMPORTANCE Athlete's foot, jock itch, ringworm, and nail infections are common fungal infections, all caused by fungi known as dermatophytes (fungi that infect skin). This report presents the genome sequences of Trichophyton rubrum, the most frequent cause of athlete's foot, as well as four other common dermatophytes. Dermatophyte genomes are enriched for four gene classes that may contribute to the ability of these fungi to cause disease. These include (i) proteases secreted to degrade skin; (ii) kinases, including pseudokinases, that are involved in signaling necessary for adapting to skin; (iii) secondary metabolites, compounds that act as toxins or signals in the interactions between fungus and host; and (iv) a class of proteins (LysM) that appear to bind and mask cell wall components and carbohydrates, thus avoiding the host's immune response to the fungi. These genome sequences provide a strong foundation for future work in understanding how dermatophytes cause disease.

Effect of pulse repetition rate and number of pulses in the analysis of polypropylene and high density polyethylene by nanosecond infrared laser induced breakdown spectroscopy

Pulse repetition rates and the number of laser pulses are among the most important parameters that do affect the analysis of solid materials by laser induced breakdown spectroscopy, and the knowledge of their effects is of fundamental importance for suggesting analytical strategies when dealing with laser ablation processes of polymers. In this contribution, the influence of these parameters in the ablated mass and in the features of craters was evaluated in polypropylene and high density polyethylene plates containing pigment-based PbCrO4. Surface characterization and craters profile were carried out by perfilometry and scanning electron microscopy. Area, volume and profile of craters were obtained using Taylor Map software. A laser induced breakdown spectroscopy system consisted of a Q-Switched Nd:YAG laser (1064 nm, 5 ns) and an Echelle spectrometer equipped with ICCD detector were used. The evaluated operating conditions consisted of 10, 25 and 50 laser pulses at 1, 5 and 10 Hz, 250 mJ/pulse (85 J cm(-2)), 2 mu s delay time and 6 mu s integration time gate. Differences in the topographical features among craters of both polymers were observed. The decrease in the repetition rate resulted in irregular craters and formation of edges, especially in polypropylene sample. The differences in the topographical features and ablated masses were attributed to the influence of the degree of crystallinity, crystalline melting temperature and glass transition temperature in the ablation process of the high density polyethylene and polypropylene. It was also observed that the intensities of chromium and lead emission signals obtained at 10 Hz were two times higher than at 5 Hz by keeping the number of laser pulses constant. (C) 2011 Elsevier B. V. All rights reserved.

Combinatorial-topological framework for the analysis of global dynamics

We discuss an algorithmic framework based on efficient graph algorithms and algebraic-topological computational tools. The framework is aimed at automatic computation of a database of global dynamics of a given m-parameter semidynamical system with discrete time on a bounded subset of the n-dimensional phase space. We introduce the mathematical background, which is based upon Conley's topological approach to dynamics, describe the algorithms for the analysis of the dynamics using rectangular grids both in phase space and parameter space, and show two sample applications. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4767672]

The I405V and Taq1B polymorphisms of the CETP gene differentially affect sub-clinical carotid atherosclerosis

Background: Cholesteryl ester transfer protein (CETP) plays a major role in lipid metabolism, but studies on the association of CETP polymorphisms with risks of cardiovascular disease are inconsistent. This study investigated whether the CETP gene I405V and Taq1B polymorphisms modified subclinical atherosclerosis in an asymptomatic Brazilian population sample. Methods: The polymorphisms were analyzed using polymerase chain reaction in 207 adult volunteers. Serum lipid profiles, oxLDL Ab titers, C-reactive protein and tumor necrosis factor-a concentrations and CETP and phospholipid transfer protein (PLTP) activities were determined, and common carotid artery intima-media thickness (cIMT) was measured using ultrasonography. Results: No differences in cIMT were observed between the presence or absence of the minor B2 and V alleles in either polymorphism. However, inverse correlations between mean cIMT and CETP activity in the presence of these polymorphisms were observed, and positive correlations of these polymorphisms with PLTP activity and oxLDL Ab titers were identified. Moreover, logistic multivariate analysis revealed that the presence of the B2 allele was associated with a 5.1-fold (CI 95%, OR: 1.26 - 21.06) increased risk for cIMT, which was equal and above the 66th percentile and positively interacted with age. However, no associations with the V allele or CETP and PLTP activities were observed. Conclusions: None of the studied parameters, including CETP activity, explained the different relationships between these polymorphisms and cIMT, suggesting that other non-determined factors were affected by the genotypes and related to carotid atherosclerotic disease.

Germline mutation analysis of Tpit in Poodle dogs with ACTH-dependent hyperadrenocorticism

There is a high incidence of pituitary-dependent hyperadrenocorticism (PDH) in Poodle dogs, with family members being affected by the disease, suggesting a genetic involvement. Tpit is an obligate transcription factor for the expression of pro-opiomelanocortingene and for corticotroph terminal differentiation. The aim of the present study was to screen the Tpit gene for germline mutations in Poodles with PDH. Fifty Poodle dogs (33 female, 8.71 +/- 2.8 years) with PDH and 50 healthy Poodle dogs (32 females, 9.4241 2.8 years) were studied. Genomic DNA was isolated from peripheral blood, amplified by PCR and submitted to automatic sequence. No mutation in the coding region of Tpit was found, whereas the new single nucleotide polymorphism p.S343G, in heterozygous state, was found in the same frequency in both PDH and control groups. We concluded that Tpit gain-of-function mutations are not involved in the etiology of PDH in Poodle dogs.

A simple analysis of the influence of the solvent-induced electronic polarization on the N-15 magnetic shielding of pyridine in water

Electronic polarization induced by the interaction of a reference molecule with a liquid environment is expected to affect the magnetic shielding constants. Understanding this effect using realistic theoretical models is important for proper use of nuclear magnetic resonance in molecular characterization. In this work, we consider the pyridine molecule in water as a model system to briefly investigate this aspect. Thus, Monte Carlo simulations and quantum mechanics calculations based on the B3LYP/6-311++G (d,p) are used to analyze different aspects of the solvent effects on the N-15 magnetic shielding constant of pyridine in water. This includes in special the geometry relaxation and the electronic polarization of the solute by the solvent. The polarization effect is found to be very important, but, as expected for pyridine, the geometry relaxation contribution is essentially negligible. Using an average electrostatic model of the solvent, the magnetic shielding constant is calculated as -58.7 ppm, in good agreement with the experimental value of -56.3 ppm. The explicit inclusion of hydrogen-bonded water molecules embedded in the electrostatic field of the remaining solvent molecules gives the value of -61.8 ppm.

Pedigree analysis and inbreeding depression on growth traits in Brazilian Marchigiana and Bonsmara breeds

The study of population structure by pedigree analysis is useful to identify important circumstances that affect the genetic history of populations. The intensive use of a small number of superior individuals may reduce the genetic diversity of populations. This situation is very common for the beef cattle breeds. Therefore, the objectives of the present study were to analyze the pedigree and possible inbreeding depression on traits of economic interest in the Marchigiana and Bonsmara breeds and to test the inclusion of the individual inbreeding coefficient (F-i) or individual increases in inbreeding coefficient (Delta F-i) in the genetic evaluation model for the quantification of inbreeding depression. The complete pedigree file of the Marchigiana breed included 29,411 animals born between 1950 and 2003. For the Bonsmara breed, the pedigree file included 18,695 animals born between 1988 and 2006. Only animals with at least 2 equivalent generations of known pedigree were kept in the analyses of inbreeding effect on birth weight, weaning weight measured at about 205 d, and BW at 14 mo in the Marchigiana breed, and on birth weight, weaning weight, and scro-tal circumference measured at 12 mo in the Bonsmara breed. The degree of pedigree knowledge was greater for Marchigiana than for Bonsmara animals. The average generation interval was 7.02 and 3.19 for the Marchigiana and Bonsmara breed, respectively. The average inbreeding coefficient was 1.33% for Marchigiana and 0.26% for Bonsmara. The number of ancestors explaining 50% of the gene pool and effective population size computed via individual increase in coancestry were 13 and 97.79 for Marchigiana and 41 and 54.57 for Bonsmara, respectively. These estimates indicate reduction in genetic variability in both breeds. Inbreeding depression was observed for most of the growth traits. The model including Delta F-i can be considered more adequate to quantify inbreeding depression. The inclusion of F-i or Delta F-i in the genetic evaluation model may not result in better fit to the data. A genetic evaluation with simultaneous estimation of inbreeding depression can be performed in Marchigiana and Bonsmara breeds, providing additional information to producers and breeders.

Development of a high sensitive automatic setup for screening of microcystins in surface waters by employing a LED-based photometric detector

In this manuscript, an automatic setup for screening of microcystins in surface waters by employing photometric detection is described. Microcystins are toxins delivered by cyanobacteria within an aquatic environment, which have been considered strongly poisonous for humans. For that reason, the World Health Organization (WHO) has proposed a provisional guideline value for drinking water of 1 mu g L-1. In this work, we developed an automated equipment setup, which allows the screening of water for concentration of microcystins below 0.1 mu g V. The photometric method was based on the enzyme-linked immunosorbent assay (ELISA) and the analytical signal was monitored at 458 nm using a homemade LED-based photometer. The proposed system was employed for the detection of microcystins in rivers and lakes waters. Accuracy was assessed by processing samples using a reference method and applying the paired t-test between results. No significant difference at the 95% confidence level was observed. Other useful features including a linear response ranging from 0.05 up to 2.00 mu g L-1 (R-2 =0.999) and a detection limit of 0.03 mu g L-1 microcystins were achieved. (C) 2011 Elsevier B.V. All rights reserved.

Parametric Analysis of Effective Material Properties of Thickness-Shear Piezoelectric Macro-Fibre Composites

A previous study on the characterization of effective material properties of a d(15) thickness-shear piezoelectric Macro-Fibre Composite (MFC) made of seven layers (Kapton, Acrylic, Electrode, Piezoceramic Fibre and Epoxy Composite, Electrode, Acrylic, Kapton) using a finite element homogenization method has shown that the packaging reduces significantly the shear stiffness of the piezoceramic material and, thus, leads to significantly smaller effective electromechanical coupling coefficient k(15) and piezoelectric stress constant e(15) when compared to the piezoceramic fibre properties. Therefore, the main objective of this work is to perform a parametric analysis in which the effect of the variations of fibre volume fraction, Epoxy elastic modulus, electrode thickness and active layer thickness on the MFC effective material properties is evaluated. Results indicate that an effective d(15) MFC should use relatively thick fibres having relatively high shear modulus and relatively stiff epoxy filler. On the other hand, the electrode thickness does not affect significantly the MFC performance.

Analysis of nonlinear dynamics of vocal folds using high-speed video observation and biomechanical modeling

Primary voice production occurs in the larynx through vibrational movements carried out by vocal folds. However, many problems can affect this complex system resulting in voice disorders. In this context, time-frequency-shape analysis based on embedding phase space plots and nonlinear dynamics methods have been used to evaluate the vocal fold dynamics during phonation. For this purpose, the present work used high-speed video to record the vocal fold movements of three subjects and extract the glottal area time series using an image segmentation algorithm. This signal is used for an optimization method which combines genetic algorithms and a quasi-Newton method to optimize the parameters of a biomechanical model of vocal folds based on lumped elements (masses, springs and dampers). After optimization, this model is capable of simulating the dynamics of recorded vocal folds and their glottal pulse. Bifurcation diagrams and phase space analysis were used to evaluate the behavior of this deterministic system in different circumstances. The results showed that this methodology can be used to extract some physiological parameters of vocal folds and reproduce some complex behaviors of these structures contributing to the scientific and clinical evaluation of voice production. (C) 2010 Elsevier Inc. All rights reserved.