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.      

Freezing in Conifer Xylem: II. PIT ASPIRATION AND BUBBLE FORMATION

A scanning electron microscope equipped with a freezing stagewas used to examine the effects of slow freezing on pit aspirationand bubble formation in living tree stems. The size (approximately 2.0 µm diameter) and the sphericalor ellipsoidal shape of the bubbles found in the centre of frozenlumens indicated freezing rates greater than 25 µm s^–1.Both unaspirated and aspirated bordered pits were found in thefrozen xylem. The technique used did not reveal enough pitsto determine whether unaspirated pits were more prevalent thanaspirated pits. These results are compared with hypotheses and results fromprevious work on freezing in conifer xylem. Key words: Freezing, conifer xylem, bordered pits, bubbles     

Expression of ATP6V1C1 during oral carcinogenesis

We investigated the gene and protein expressions of V-type ATPase protein subunit C1 (ATP6V1C1) in cases of oral squamous cell carcinoma (OSCC) and contralateral normal mucosa in smokers, nonsmokers and former smokers. Subjects were separated into five groups of 15: group 1, smokers with OSCC; group 2, normal contralateral mucosa of OSCC patients; group 3, chronic smokers; group 4, former smokers who had stopped smoking 1 year earlier; group 5, individuals who had never smoked. Exfoliative cytology specimens from oral mucosa of smokers, former smokers and nonsmokers showed normal gene and protein expression. We found significantly greater gene expression in the OSCC group than in the nonsmoker groups. No difference in gene expression was observed between normal contralateral mucosa and nonsmoker groups, smoker and nonsmoker groups or former smoker and nonsmoker groups. We observed intense immunostaining for ATP6V1C1 protein in all cases of OSCC and weak or no staining in smoker, former smoker and nonsmoker groups. Significantly greater expression of ATP6V1C1 protein was observed in the OSCC group compared to the other groups, which supports the role of ATP6V1C1 in effecting changes associated with oral cancer. Analysis of the mucosae of chronic smokers, former smokers and the normal contralateral mucosa of patients with OSCC showed unaltered ATP6V1C1 gene and protein expression. Early stages of carcinogenesis, represented by altered epithelium of chronic smokers, had neither gene nor protein alterations as seen in OSCC. Therefore, we infer that the changes in ATP6V1C1 occur during later stages of carcinogenesis. Our preliminary study provides a basis for future studies of using ATP6V1C1 levels for detecting early stage OSCC.     

Effects of age,replicative lifespan and growth rate of human nucleus pulposus cells on selecting age range for cell-based biological therapies for degenerative disc diseases

Autologous disc cell implantation, growth factors and gene therapy appear to be promising therapies for disc regeneration. Unfortunately, the replicative lifespan and growth kinetics of human nucleus pulposus (NP) cells related to host age are unclear. We investigated the potential relations among age, replicative lifespan and growth rate of NP cells, and determined the age range that is suitable for cell-based biological therapies for degenerative disc diseases. We used NP tissues classified by decade into five age groups: 30s, 40s, 50s, 60s and 70s. The mean cumulative population doubling level (PDL) and population doubling rate (PDR) of NP cells were assessed by decade. We also investigated correlations between cumulative PDL and age, and between PDR and age. The mean cumulative PDL and PDR decreased significantly in patients in their 60s. The mean cumulative PDL and PDR in the younger groups (30s, 40s and 50s) were significantly higher than those in the older groups (60s and 70s). There also were significant negative correlations between cumulative PDL and age, and between PDR and age. We found that the replicative lifespan and growth rate of human NP cells decreased with age. The replicative potential of NP cells decreased significantly in patients 60 years old and older. Young individuals less than 60 years old may be suitable candidates for NP cell-based biological therapies for treating degenerative disc diseases.     

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.     

中国金丝桃属金丝桃组(藤黄科)新分类群

描述了中国云南北部的金丝桃属Hypericum L.(藤黄科)两个新种,即楚雄金丝桃H. fosteri N. Robson和漾濞金丝桃H. wardianum N. Robson.二者均与短柱金丝桃H. hookerianum Wight & Arn.相近.楚雄金丝桃为一地方特有种,目前仅在模式地点发现;漾濞金丝桃在缅甸东北部还有分布.二者主要以蒴果伸长,非卵状近球形,萼片无中脉而有别于短柱金丝桃.此外,将宽萼金丝桃H. latisepalum (N. Robson) N. Robson从亚种提升为种;将变异较大的纤枝金丝桃H. lagarocladum N. Robson划分为两个亚种,即叶较宽、具有开张习性的原亚种纤枝金丝桃H. lagarocladum ssp. lagarocladum和叶较窄、更为直立的新亚种狭叶金丝桃H. lagarocladum ssp. angustifolium N. Robson.     

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.     

Generation of whole genome sequences of new Cryptosporidium hominis and Cryptosporidium parvum isolates directly from stool samples

Background

Whole genome sequencing (WGS) of Cryptosporidium spp. has previously relied on propagation of the parasite in animals to generate enough oocysts from which to extract DNA of sufficient quantity and purity for analysis. We have developed and validated a method for preparation of genomic Cryptosporidium DNA suitable for WGS directly from human stool samples and used it to generate 10 high-quality whole Cryptosporidium genome assemblies. Our method uses a combination of salt flotation, immunomagnetic separation (IMS), and surface sterilisation of oocysts prior to DNA extraction, with subsequent use of the transposome-based Nextera XT kit to generate libraries for sequencing on Illumina platforms. IMS was found to be superior to caesium chloride density centrifugation for purification of oocysts from small volume stool samples and for reducing levels of contaminant DNA.

Results

The IMS-based method was used initially to sequence whole genomes of Cryptosporidium hominis gp60 subtype IbA10G2 and Cryptosporidium parvum gp60 subtype IIaA19G1R2 from small amounts of stool left over from diagnostic testing of clinical cases of cryptosporidiosis. The C. parvum isolate was sequenced to a mean depth of 51.8X with reads covering 100 % of the bases of the C. parvum Iowa II reference genome (Bioproject PRJNA 15586), while the C. hominis isolate was sequenced to a mean depth of 34.7X with reads covering 98 % of the bases of the C. hominis TU502 v1 reference genome (Bioproject PRJNA 15585).The method was then applied to a further 17 stools, successfully generating another eight new whole genome sequences, of which two were C. hominis (gp60 subtypes IbA10G2 and IaA14R3) and six C. parvum (gp60 subtypes IIaA15G2R1 from three samples, and one each of IIaA17G1R1, IIaA18G2R1, and IIdA22G1), demonstrating the utility of this method to sequence Cryptosporidium genomes directly from clinical samples. This development is especially important as it reduces the requirement to propagate Cryptosporidium oocysts in animal models prior to genome sequencing.

Conclusion

This represents the first report of high-quality whole genome sequencing of Cryptosporidium isolates prepared directly from human stool samples.