Bayesian adaptive sequence alignment algorithms

The selection of a scoring matrix and gap penalty parameters continues to be an important problem in sequence alignment. We describe here an algorithm, the 'Bayes block aligner, which bypasses this requirement. Instead of requiring a fixed set of parameter settings, this algorithm returns the Bayesian posterior probability for the number of gaps and for the scoring matrices in any series of interest. Furthermore, instead of returning the single best alignment for the chosen parameter settings, this algorithm returns the posterior distribution of all alignments considering the full range of gapping and scoring matrices selected, weighing each in proportion to its probability based on the data. We compared the Bayes aligner with the popular Smith-Waterman algorithm with parameter settings from the literature which had been optimized for the identification of structural neighbors, and found that the Bayes aligner correctly identified more structural neighbors. In a detailed examination of the alignment of a pair of kinase and a pair of GTPase sequences, we illustrate the algorithm's potential to identify subsequences that are conserved to different degrees. In addition, this example shows that the Bayes aligner returns an alignment-free assessment of the distance between a pair of sequences.      

Effects of plant height and row spacing on kenaf forage potential with multiple harvests

Kenaf (Hibiscus cannabinus L.) forage potential can be enhanced through its regrowth capacity and higher production in narrow rows. A field experiment was conducted in Matamoros, Coahuila, Mexico, during 2 growing seasons (2004 and 2005) to study the effects of plant height and row spacing on kenaf forage potential with multiple harvests. This study evaluated the effects of (1) 2 plant heights at cutting (1.0-1.2 m and 1.8-2.0 m) and (2) 4 inter row spacings (0.19, 0.38, 0.57 and 0.76 m) using a 2 x 4 factorial arrangement of treatments in a completely randomized block design with 4 replications. Dry matter (DM) and crude protein (CP) yields, DM partitioning, neutral detergent fiber (NDF) and CP concentrations were determined. Heights at cutting × row spacing interactions were not significant for the monitored variables (p>0.05). Kenaf response to treatments was only relevant for main effects (p≤0.05). Row spacing and plant height affected DM and CP yields (p≤0.05), whereas only plant height affected chemical composition and DM partitioning (p≤0.05). Dry matter (17.0%-26.0%), and CP (12.4%-15.6%) yields were higher (p≤0.05) when plant heights had reached 1.8 to 2.0 m. Row spacing reduction from 0.76 m to 0.38 and 0.19 m increased DM yield (20.4-33.4%) and CP yield (24.2-38.5%) (p≤0.05). Kenaf forage potential increases when planted in narrow rows and harvested 2 or 3 times during the growing season.     

Vegetal cover estimated by digital photographs related to biomass in a grassland site in northern Mexico

A regression model was used to determine the relationship between aerial herbaceous biomass and vegetation coverage estimated by digital images. Four samplings (n=36 each date) of vegetation cover and herbaceous biomass were performed during the growing season in 2011 in a grassland dominated by Bouteloua gracilis in La Cieneguilla, Municipality of Villa Hidalgo, Durango. Average production of dry biomass was 37.36 ± 9.66 g/m², and mean vegetation cover 30.02%. Dry biomass data were tested for normality using the test of Kolmogorov Smirnov, finding a lack of fit. The data were subjected to a logarithmic transformation and the model Ln(y) = 1.637926 + 0.08501X - 0.000586X² with an adjusted R² = 0.89 was found. In order to validate this model, another five samplings were carried out in 2013 at the same site during summer and autumn, using the same sampling size for each date as in 2011. Data collected in 2013 were analyzed with the model Ln (y) = β₀ + β₁X + β₂X². A comparison of regression coefficients was carried out between the 2011 and 2013 models with t (180+144-9-11-2=302, p<0.05) = 1.967. The results indicated that it is possible to use the 2011 regression model to estimate herbaceous aerial biomass from vegetation cover measurements with aerial photographs in La Cieneguilla site during summer and fall.     

A novel scheme to assess factors involved in the reproducibility of DNA-microarray data

Background  

In research laboratories using DNA-microarrays, usually a number of researchers perform experiments, each generating possible sources of error. There is a need for a quick and robust method to assess data quality and sources of errors in DNA-microarray experiments. To this end, a novel and cost-effective validation scheme was devised, implemented, and employed.     

Quantitative proteomic analysis of amniocytes reveals potentially dysregulated molecular networks in Down syndrome

Background

Down syndrome (DS), caused by an extra copy of chromosome 21, affects 1 in 750 live births and is characterized by cognitive impairment and a constellation of congenital defects. Currently, little is known about the molecular pathogenesis and no direct genotype-phenotype relationship has yet been confirmed. Since DS amniocytes are expected to have a distinct biological behaviour compared to normal amniocytes, we hypothesize that relative quantification of proteins produced from trisomy and euploid (chromosomally normal) amniocytes will reveal dysregulated molecular pathways.

Results

Chromosomally normal- and Trisomy 21-amniocytes were quantitatively analyzed by using Stable Isotope Labeling of Amino acids in Cell culture and tandem mass spectrometry. A total of 4919 unique proteins were identified from the supernatant and cell lysate proteome. More specifically, 4548 unique proteins were identified from the lysate, and 91% of these proteins were quantified based on MS/MS spectra ratios of peptides containing isotope-labeled amino acids. A total of 904 proteins showed significant differential expression and were involved in 25 molecular pathways, each containing a minimum of 16 proteins. Sixty of these proteins consistently showed aberrant expression from trisomy 21 affected amniocytes, indicating their potential role in DS pathogenesis. Nine proteins were analyzed with a multiplex selected reaction monitoring assay in an independent set of Trisomy 21-amniocyte samples and two of them (SOD1 and NES) showed a consistent differential expression.

Conclusions

The most extensive proteome of amniocytes and amniotic fluid has been generated and differentially expressed proteins from amniocytes with Trisomy 21 revealed molecular pathways that seem to be most significantly affected by the presence of an extra copy of chromosome 21.     

Chromosome phylogeny of the subfamily Pitheciinae (Platyrrhini,Primates) by classic cytogenetics and chromosome painting

Background  

The New World monkey (Platyrrhini) subfamily Pitheciinae is represented by the genera Pithecia, Chiropotes and Cacajao. In this work we studied the karyotypes of Pithecia irrorata (2n = 48) and Cacajao calvus rubicundus (2n = 45 in males and 2n = 46 in females) by G- and C-banding, NOR staining and chromosome painting using human and Saguinus oedipus whole chromosome probes. The karyotypes of both species were compared with each other and with Chiropotes utahicki (2n = 54) from the literature.     

Molecular modeling and analysis of human and plant endo-β-N-acetyl- glucosaminidases for mutations effects on function

Endo- β-N-acetylgucosaminidases (ENGases) are the enzymes that catalyze both hydrolysis and transglycosylation reactions. It is of interest to study ENGases because of their ability to synthesize glycopeptides. Homology models of Human, Arabidopsis thaliana and Sorghum ENGases were developed and their active sites marked based on information available from Arthrobacter protophormiae (PDB ID: 3FHQ) ENGase. Further, these models were docked with the natural substrate GlcNAc-Asn and the inhibitor Man₃GlcNAc-thiazoline. The catalytic triad of Asn, Glu and Tyr (N171, E173 and Y205 of bacteria) were found to be conserved across the phyla. The crucial Y299F mutation showing 3 times higher transglycosylation activity than in wild type Endo-A is known. The hydrolytic activity remained unchanged in bacteria, while the transglycosylation activity increased. This Y to F change is found to be naturally evolved and should be attributing higher transglycosylation rates in human and Arabidopsis thaliana ENGases. Ligand interactions Ligplots revealed the interaction of amino acids with hydrophobic side chains and polar uncharged side chain amino acids. Thus, structure based molecular model-ligand interactions provide insights into the catalytic mechanism of ENGases and assist in the rational engineering of ENGases.