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.      

Enkephalin-encoding herpes simplex virus-1 decreases inflammation and hotplate sensitivity in a chronic pancreatitis model

Background

A chronic pancreatitis model was developed in young male Lewis rats fed a high-fat and alcohol liquid diet beginning at three weeks. The model was used to assess time course and efficacy of a replication defective herpes simplex virus type 1 vector construct delivering human cDNA encoding preproenkephalin (HSV-ENK).

Results

Most surprising was the relative lack of inflammation and tissue disruption after HSV-ENK treatment compared to the histopathology consistent with pancreatitis (inflammatory cell infiltration, edema, acinar cell hypertrophy, fibrosis) present as a result of the high-fat and alcohol diet in controls. The HSV-ENK vector delivered to the pancreatic surface at week 3 reversed pancreatitis-associated hotplate hypersensitive responses for 4–6 weeks, while control virus encoding β-galactosidase cDNA (HSV-β-gal) had no effect. Increased Fos expression seen bilaterally in pain processing regions in control animals with pancreatitis was absent in HSV-ENK-treated animals. Increased met-enkephalin staining was evident in pancreas and lower thoracic spinal cord laminae I–II in the HSV-ENK-treated rats.

Conclusion

Thus, clear evidence is provided that site specific HSV-mediated transgene delivery of human cDNA encoding preproenkephalin ameliorates pancreatic inflammation and significantly reduces hypersensitive hotplate responses for an extended time consistent with HSV mediated overexpression, without tolerance or evidence of other opiate related side effects.     

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.     

M5 muscarinic receptors are needed for slow activation of dopamine neurons and for rewarding brain stimulation

Mesopontine cholinergic neurons (Ch5 and Ch6 cell groups) activate the cerebral cortex via thalamic projections, and activate locomotion and reward via dopamine neurons in the substantia nigra and ventral tegmental area (VTA). Nicotinic receptors in VTA activate dopamine neurons quickly, and are needed for the stimulant and rewarding effects of nicotine in rats. Muscarinic receptors in VTA activate dopamine neurons slowly, and are needed for the rewarding effects of hypothalamic stimulation, but do not increase locomotion. Antisense oligonucleotides targetting M5 mRNA, when infused into the VTA, inhibited M5 receptor binding and rewarding hypothalamic stimulation. Mutant mice with truncated M5 muscarinic receptor genes drank more water than wild-type controls. Spontaneous locomotion and locomotor responses to amphetamine and scopolamine were unchanged. Electrical stimulation near Ch6 induced dopamine release in the nucleus accumbens in two phases, an early phase (0-2 min after stimulation) dependent on nicotinic and gluatamatergic receptors in VTA, and a late phase (8-50 min after stimulation) dependent on muscarinic receptors in VTA. The late phase was lost in M5 mutant mice, while the early phase was unchanged. M5 muscarinic receptors bind slowly to muscarinic ligands, and appear to mediate slow secretions.     

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.