ProbFAST: Probabilistic Functional Analysis System Tool

Background: The post-genomic era has brought new challenges regarding the understanding of the organization and function of the human genome. Many of these challenges are centered on the meaning of differential gene regulation under distinct biological conditions and can be performed by analyzing the Multiple Differential Expression (MDE) of genes associated with normal and abnormal biological processes. Currently MDE analyses are limited to usual methods of differential expression initially designed for paired analysis. Results: We proposed a web platform named ProbFAST for MDE analysis which uses Bayesian inference to identify key genes that are intuitively prioritized by means of probabilities. A simulated study revealed that our method gives a better performance when compared to other approaches and when applied to public expression data, we demonstrated its flexibility to obtain relevant genes biologically associated with normal and abnormal biological processes. Conclusions: ProbFAST is a free accessible web-based application that enables MDE analysis on a global scale. It offers an efficient methodological approach for MDE analysis of a set of genes that are turned on and off related to functional information during the evolution of a tumor or tissue differentiation. ProbFAST server can be accessed at http://gdm.fmrp.usp.br/probfast.

New Source of Muscle-Derived Stem Cells with Potential for Alveolar Bone Reconstruction in Cleft Lip and/or Palate Patients

Cleft lip and palate (CLP), one of the most frequent congenital malformations, affects the alveolar bone in the great majority of the cases, and the reconstruction of this defect still represents a challenge in the rehabilitation of these patients. One of the current most promising strategy to achieve this goal is the use of bone marrow stem cells (BMSC); however, isolation of BMSC or iliac bone, which is still the mostly used graft in the surgical repair of these patients, confers site morbidity to the donor. Therefore, in order to identify a new alternative source of stem cells with osteogenic potential without conferring morbidity to the donor, we have used orbicular oris muscle (OOM) fragments, which are regularly discarded during surgery repair (cheiloplasty) of CLP patients. We obtained cells from OOM fragments of four unrelated CLP patients (CLPMDSC) using previously described preplating technique. These cells, through flow cytometry analysis, were mainly positively marked for five mesenchymal stem cell antigens (CD29, CD90, CD105, SH3, and SH4), while negative for hematopoietic cell markers, CD14, CD34, CD45, and CD117, and for endothelial cell marker, CD31. After induction under appropriate cell culture conditions, these cells were capable to undergo chondrogenic, adipogenic, osteogenic, and skeletal muscle cell differentiation, as evidenced by immunohistochemistry. We also demonstrated that these cells together with a collagen membrane lead to bone tissue reconstruction in a critical-size cranial defects previously induced in non-immunocompromised rats. The presence of human DNA in the new bone was confirmed by PCR with human-specific primers and immunohistochemistry with human nuclei antibodies. In conclusion, we showed that cells from OOM have phenotypic and behavior characteristics similar to other adult stem cells, both in vitro and in vivo. Our findings suggest that these cells represent a promising source of stem cells for alveolar bone grafting treatment, particularly in young CLP patients.

Genetic parameters for pre-weaning traits in Braunvieh cattle

The objective of this study was to estimate genetic parameters for pre-weaning traits of Braunvieh cattle raised under tropical conditions in Brazil. The weight and weight gain parameters were birth weight (BW, N = 9955), weight at 120 days of age (W120, N = 5901), weaning weight at 205 days (WW, N = 6970), weight gain from birth to 205 days (GAIN205, N = 6013), weight gain from birth to 120 days (GAIN120, N = 5135), and weight gain from 120 to 205 days (GAIN85, N = 4482). Variance components were estimated using the animal model with the MTDFREML software. The relationship matrix included 35,188 animals; phenotypic measures were available for 18,688. Direct and maternal heritability increased from birth to weaning, with estimates of 0.23 +/- 0.037, 0.25 +/- 0.050, 0.41 +/- 0.059 for direct heritability for BW, W120 and WW, respectively, 0.08 +/- 0.012, 0.15 +/- 0.032, 0.22 +/- 0.036 for maternal genetic effects, and 0.18, 0.14 and 0.16 for total heritability estimates. For pre-weaning gains, estimates of heritability were 0.36 +/- 0.059, 0.30 +/- 0.059, 0.12 +/- 0.035 for direct genetic effects of the traits GAIN205, GAIN120 and GAIN85, respectively, 0.23 +/- 0.038, 0.17 +/- 0.037, 0.03 +/- 0.029 for estimates of maternal heritability, and 0.12, 0.13, 0.16 for total heritability, respectively. Genetic correlations between weights were greater between measures taken at shorter intervals. This information can be used to optimize the design of programs for genetic improvement of Braunvieh cattle raised under tropical conditions.

Identification of three distinguishable phenotypes in golden retriever muscular dystrophy

Duchenne muscular dystrophy (DMD) is a human disease characterized by progressive and irreversible skeletal muscle degeneration caused by mutations in genes coding for important muscle proteins. Unfortunately, there is no efficient treatment for this disease; it causes progressive loss of motor and muscular ability until death. The canine model (golden retriever muscular dystrophy) is similar to DMD, showing similar clinical signs. Fifteen dogs were followed from birth and closely observed for clinical signs. Dogs had their disease status confirmed by polymerase chain reaction analysis and genotyping. Clinical observations of musculoskeletal, morphological, gastrointestinal, respiratory, cardiovascular, and renal features allowed us to identify three distinguishable phenotypes in dystrophic dogs: mild (grade I), moderate (grade II) and severe (grade III). These three groups showed no difference in dystrophic alterations of muscle morphology and creatine kinase levels. This information will be useful for therapeutic trials, because DMD also shows significant, inter- and intra-familiar clinical variability. Additionally, being aware of phenotypic differences in this animal model is essential for correct interpretation and understanding of results obtained in pre-clinical trials.

Genetic analysis of average annual productivity of Nellore breeding cows (COWPROD)

With the aim of estimating the coefficient of heritability of average annual productivity of Nellore cows (COWPROD), a data set from 24,855 animals with known pedigree was analyzed. COWPROD is defined as the amount (in kilograms) of weaned calves produced yearly by one cow during her remaining time in herd ignoring a fixed period of 365 days. COWPROD was calculated regarding three standards: a) based on the post-weaning weight from the calves ignoring any kind of adjustment (COWPROD_NAJ), b) adjusted weight for the fixed effects (COWPROD_AJFIX) and c) adjusted weight for the fixed effects and for the genetic merit of the sire (COWPROD_AJFIN). The obtained heritabilities were 0.15, 0.15 and 0.16 for COWPROD_NAJ, COWPROD_AJFIX and COWPROD_AJFIN, respectively. A complete set composed of 105,158 COWPROD records on 130,740 animals in pedigree was also analyzed for predicting the genetic merit of all animals in the data set and for the calculation of the genetic, phenotypic and residual trends. Ranking correlation was high for the adjusted and non-adjusted data, yet, for some of the animals, the difference among the genetic values was large. This would be an indication that it would be better to work always with the adjusted weaning weights. The genetic trend was positive, but was of small magnitude (0.26% of the trait average) and the residual trend was negative as a consequence of the large intensification of the production system, which has been occurring in the last years in the farms studied. The phenotypic trend was also negative and intermediate between the genetic and the residual ones.

Additive genetic relationship of longevity with fertility and production traits in Nellore cattle based on bivariate models

Survival or longevity is an economically important trait in beef cattle. The main inconvenience for its inclusion in selection criteria is delayed recording of phenotypic data and the high computational demand for including survival in proportional hazard models. Thus, identification of a longevity-correlated trait that could be recorded early in life would be very useful for selection purposes. We estimated the genetic relationship of survival with productive and reproductive traits in Nellore cattle, including weaning weight (WW), post-weaning growth (PWG), muscularity (MUSC), scrotal circumference at 18 months (SC18), and heifer pregnancy (HP). Survival was measured in discrete time intervals and modeled through a sequential threshold model. Five independent bivariate Bayesian analyses were performed, accounting for cow survival and the five productive and reproductive traits. Posterior mean estimates for heritability (standard deviation in parentheses) were 0.55 (0.01) for WW, 0.25 (0.01) for PWG, 0.23 (0.01) for MUSC, and 0.48 (0.01) for SC18. The posterior mean estimates (95% confidence interval in parentheses) for the genetic correlation with survival were 0.16 (0.13-0.19), 0.30 (0.25-0.34), 0.31 (0.25-0.36), 0.07 (0.02-0.12), and 0.82 (0.78-0.86) for WW, PWG, MUSC, SC18, and HP, respectively. Based on the high genetic correlation and heritability (0.54) posterior mean estimates for HP, the expected progeny difference for HP can be used to select bulls for longevity, as well as for post-weaning gain and muscle score.

Genetic parameters for postweaning traits in Braunvieh cattle

Genetic parameters for traits related to postweaning growth in Braunvieh cattle, reared under tropical and sub-tropical conditions in Brazil, were studied. Weight traits were weight at 365 days of age (W365, N = 4055), at 450 days (W450, N = 3453), and at 550 days (W550, N = 1946), while weight gains were gain from weaning to 365 days of age (WGW365, N = 3060), from weaning to 450 days (WGW450, N = 2764), from weaning to 550 days (WGW550, N = 1531), from 365 to 550 days of age (WG365550, N = 1528), from 365 to 450 days (WG365450, N = 2401), and from 450 to 550 days (WG450550, N = 1563). A full animal model was used for estimating the variance components, using the MTDFREML software. The dataset contained 18,688 animals with phenotypic measures and 35,188 animals in the relationship matrix. Heritability estimates for postweaning weights decreased with age. For W365, W450 and W550, respectively, the direct heritability estimates were 0.29 +/- 0.061, 0.25 +/- 0.057, 0.16 +/- 0.060, maternal heritability was 0.20 +/- 0.035, 0.18 +/- 0.035, 0.13 +/- 0.052, and total heritability was 0.30, 0.35, 0.26. In this breed, maternal influence was found to be important up to 550 days of age. The greater genetic correlations between weights were observed for weights measured at shorter intervals. A large environmental effect was observed for weight gain between weaning and 550 days; this effect was greater for the gains between 365 and 550 days.

Comparison of different models to estimate genetic parameters for carcass traits in a commercial broiler line

The effect of genetic and non-genetic factors for carcass, breast meat and leg weights, and yields of a commercial broiler line were investigated using the restricted maximum likelihood method, considering four different animal models, including or excluding maternal genetic effect with covariance between direct and maternal genetic effects, and maternal permanent environmental effect. The likelihood ratio test was used to determine the most adequate model for each trait. For carcass, breast, and leg weight, and for carcass and breast yield, maternal genetic and permanent environmental effects as well as the covariance between direct and maternal genetic effects were significant. The estimates of direct and maternal heritability were 0.17 and 0.04 for carcass weight, 0.26 and 0.06 for breast weight, 0.22 and 0.02 for leg weight, 0.32 and 0.02 for carcass yield, and 0.52 and 0.04 for breast yield, respectively. For leg yield, maternal permanent environmental effect was important, in addition to direct genetic effects. For that trait, direct heritability and maternal permanent environmental variance as a proportion of the phenotypic variance were 0.43 and 0.02, respectively. The results indicate that ignoring maternal effects in the models, even though they were of small magnitude (0.02 to 0.06), tended to overestimate direct genetic variance and heritability for all traits.

Phylogeography of the common vampire bat (Desmodus rotundus): Marked population structure, Neotropical Pleistocene vicariance and incongruence between nuclear and mtDNA markers

Background: The common vampire bat Desmodus rotundus is an excellent model organism for studying ecological vicariance in the Neotropics due to its broad geographic range and its preference for forested areas as roosting sites. With the objective of testing for Pleistocene ecological vicariance, we sequenced a mitocondrial DNA (mtDNA) marker and two nuclear markers (RAG2 and DRB) to try to understand how Pleistocene glaciations affected the distribution of intraspecific lineages in this bat. Results: Five reciprocally monophyletic clades were evident in the mitochondrial gene tree, and in most cases with high bootstrap support: Central America (CA), Amazon and Cerrado (AMC), Pantanal (PAN), Northern Atlantic Forest (NAF) and Southern Atlantic Forest (SAF). The Atlantic forest clades formed a monophyletic clade with high bootstrap support, creating an east/west division for this species in South America. On the one hand, all coalescent and non-coalescent estimates point to a Pleistocene time of divergence between the clades. On the other hand, the nuclear markers showed extensive sharing of haplotypes between distant localities, a result compatible with male-biased gene flow. In order to test if the disparity between the mitochondrial and nuclear markers was due to the difference in mutation rate and effective size, we performed a coalescent simulation to examine the feasibility that, given the time of separation between the observed lineages, even with a gene flow rate close to zero, there would not be reciprocal monophyly for a neutral nuclear marker. We used the observed values of theta and an estimated mutation rate for the nuclear marker gene to perform 1000 iterations of the simulation. The results of this simulation were inconclusive: the number of iterations with and without reciprocal monophyly of one or more clades are similar. Conclusions: We therefore conclude that the pattern exhibited by the common vampire bat, with marked geographical structure for a mitochondrial marker and no phylogeographic structure for nuclear markers is compatible with a historical scenario of complete isolation of refuge-like populations during the Pleistocene. The results on demographic history on this species is compatible with the Carnaval-Moritz model of Pleistocene vicariance, with demographic expansions in the southern Atlantic forest.

The Liver Plays a Major Role in Clearance and Destruction of Blood Trypomastigotes in Trypanosoma cruzi Chronically Infected Mice

Intravenous challenge with Trypanosoma cruzi can be used to investigate the process and consequences of blood parasite clearance in experimental Chagas disease. One hour after intravenous challenge of chronically infected mice with 5610 6 trypomastigotes, the liver constituted a major site of parasite accumulation, as revealed by PCR. Intact parasites and/or parasite remnants were visualized at this time point scattered in the liver parenchyma. Moreover, at this time, many of liver-cleared parasites were viable, as estimated by the frequency of positive cultures, which considerably diminished after 48 h. Following clearance, the number of infiltrating cells in the hepatic tissue notably increased: initially (at 24 h) as diffuse infiltrates affecting the whole parenchyma, and at 48 h, in the form of large focal infiltrates in both the parenchyma and perivascular spaces. Phenotypic characterization of liver-infiltrating cells 24 h after challenge revealed an increase in Mac1(+), CD8(+) and CD4(+) cells, followed by natural killer (NK) cells. As evidence that liver-infiltrating CD4(+) and CD8(+) cells were activated, increased frequencies of CD69(+) CD8(+), CD69(+) CD4(+) and CD25(+) CD122(+) CD4(+) cells were observed at 24 and 48 h after challenge, and of CD25(-)CD122(+)CD4(+) cells at 48 h. The major role of CD4(+) cells in liver protection was suggested by data showing a very high frequency of interferon (IFN)-gamma-producing CD4(+) cells 24 h after challenge. In contrast, liver CD8(+) cells produced little IFN-gamma, even though they showed an enhanced potential for secreting this cytokine, as revealed by in vitro T cell receptor (TCR) stimulation. Confirming the effectiveness of the liver immune response in blood parasite control during the chronic phase of infection, no live parasites were detected in this organ 7 days after challenge.

Genome-wide SNP genotyping highlights the role of natural selection in Plasmodium falciparum population divergence

Background: The malaria parasite Plasmodium falciparum exhibits abundant genetic diversity, and this diversity is key to its success as a pathogen. Previous efforts to study genetic diversity in P. falciparum have begun to elucidate the demographic history of the species, as well as patterns of population structure and patterns of linkage disequilibrium within its genome. Such studies will be greatly enhanced by new genomic tools and recent large-scale efforts to map genomic variation. To that end, we have developed a high throughput single nucleotide polymorphism (SNP) genotyping platform for P. falciparum. Results: Using an Affymetrix 3,000 SNP assay array, we found roughly half the assays (1,638) yielded high quality, 100% accurate genotyping calls for both major and minor SNP alleles. Genotype data from 76 global isolates confirm significant genetic differentiation among continental populations and varying levels of SNP diversity and linkage disequilibrium according to geographic location and local epidemiological factors. We further discovered that nonsynonymous and silent (synonymous or noncoding) SNPs differ with respect to within-population diversity, interpopulation differentiation, and the degree to which allele frequencies are correlated between populations. Conclusions: The distinct population profile of nonsynonymous variants indicates that natural selection has a significant influence on genomic diversity in P. falciparum, and that many of these changes may reflect functional variants deserving of follow-up study. Our analysis demonstrates the potential for new high-throughput genotyping technologies to enhance studies of population structure, natural selection, and ultimately enable genome-wide association studies in P. falciparum to find genes underlying key phenotypic traits.

Molecular determinants of improved cathepsin B inhibition by new cystatins obtained by DNA shuffling

Background: Cystatins are inhibitors of cysteine proteases. The majority are only weak inhibitors of human cathepsin B, which has been associated with cancer, Alzheimer's disease and arthritis. Results: Starting from the sequences of oryzacystatin-1 and canecystatin-1, a shuffling library was designed and a hybrid clone obtained, which presented higher inhibitory activity towards cathepsin B. This clone presented two unanticipated point mutations as well as an N-terminal deletion. Reversing each point mutation independently or both simultaneously abolishes the inhibitory activity towards cathepsin B. Homology modeling together with experimental studies of the reverse mutants revealed the likely molecular determinants of the improved inhibitory activity to be related to decreased protein stability. Conclusion: A combination of experimental approaches including gene shuffling, enzyme assays and reverse mutation allied to molecular modeling has shed light upon the unexpected inhibitory properties of certain cystatin mutants against Cathepsin B. We conclude that mutations disrupting the hydrophobic core of phytocystatins increase the flexibility of the N-terminus, leading to an increase in inhibitory activity. Such mutations need not affect the inhibitory site directly but may be observed distant from it and manifest their effects via an uncoupling of its three components as a result of increased protein flexibility.

Diversity of free-living marine nematodes (Enoplida) from Baja California assessed by integrative taxonomy

We used morphological and molecular approaches to evaluate the diversity of free-living marine nematodes (order Enoplida) at four coastal sites in the Gulf of California and three on the Pacific coast of Baja California, Mexico. We identified 22 morphological species belonging to six families, of which Thoracostomopsidae and Oncholaimidae were the most diverse. The genus Mesacanthion (Thoracostomopsidae) was the most widespread and diverse. Five allopatric species, genetically and morphologically differentiated, were found in two localities in the Gulf of California (M. sp1 and M. sp2) and three in the Pacific coast (M. sp3, M. sp4 and M. sp5). Overall, we produced 19 and 20 sequences for the 18S and 28S genes, respectively. Neither gene displayed intraspecific polymorphisms, which allowed us to establish that some morphological variation was likely either ontogenetic or due to phenotypic plasticity. Although 18S and 28S phylogenies were topologically congruent (incongruence length difference test, P > 0.05), divergences between species were much higher in the 28S gene. Moreover, this gene possessed a stronger phylogenetic signal to resolve relationships involving Rhabdodemania and Bathylaimus. On the other hand, the close relationship of Pareurystomina (Enchilidiidae) with oncholaimids warrants further study. The 28S sequences (D2D3 domain) may be better suited for DNA barcoding of marine nematodes than those from the 18S rDNA, particularly for differentiating closely related or cryptic species. Finally, our results underline the relevance of adopting an integrative approach encompassing morphological and molecular analyses to improve the assessment of marine nematode diversity and advance their taxonomy.

Sequential injection analysis implementing multiple standard additions for As speciation by liquid chromatography and atomic fluorescence spectrometry (SIA-HPLC-AFS)

An analytical procedure for multiple standard additions of arsenic species using sequential injection analysis (SIA) is proposed for their quantification in seafood extracts. SIA presented flexibility for generating multiple specie standards at the ng mL(-1) concentration level by adding different volumes of As(III), As(V), monomethylarsonic (MMA) and dimethylarsinic (DMA) to the sample. The mixed sample plus standard solutions were delivered from SIA to fill the HPLC injection loop. Subsequently, As species were separated by HPLC and analyzed by atomic fluorescence spectrometry (AFS). The proposed system comprised two independently controlled modules, with the HPLC loop acting as the intermediary device. The analytical frequency was enhanced by combining the actions of both modules. While the added sample was flowing through the chromatographic column towards the detection system, the SIA program started performing the standard additions to another sample. The proposed method was applied to spoiled seafood extracts. Detection limits based on 3 sigma for As(III), As(V), MMA and DMA were 0.023, 0.39, 0.45 and 1.0 ng mL(-1), respectively. (C) 2011 Elsevier B.V. All rights reserved.

Variation in the Distribution of Trace Elements in Renal Cell Carcinoma

The development of cancer is a complex, multistage process during which a normal cell undergoes genetic changes that result in phenotypic alterations and in the acquisition of the ability to invade other sites. Inductively coupled plasma optical emission spectroscopy was used to estimate the contents of Al, Ca, Cd, Cr, Cu, Fe, K, Mg, Mn, Na, P, Pb, and Zn in healthy kidney and renal cell carcinoma (RCC), and significant differences were found for all elements. Along with the progression of the malignant disease, a progressive decrease of Cd and K was observed. In fact, for Cd, the concentration in stage T4 was 263.9 times lower than in stage T1, and for K, the concentration in stage T4 was 1.73 times lower than in stage T1. Progressive accumulation was detected for P, Pb, and Zn in stage T4. For P, the concentration in stage T4 was 11.1 times higher than in stage T1; for Pb, the concentration in stage T4 was 232.7 times higher than in T1; and for Zn, the concentration in T4 was 8.452 times higher than in T1. This study highlights the marked differences in the concentrations of selected trace metals in different malignant tumor stages. These findings indicate that some trace metals may play important roles in the pathogenesis of RCC.