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.     

Simple one-medium formulation regeneration of fingerroot [<Emphasis Type="Italic">Boesenbergia rotunda</Emphasis> (L.) mansf. Kulturpfl.] via somatic embryogenesis

Summary Most published protocols necessitate different media formulations for multistep somatic embryogenesis. This study aims to establish a simple but effective formulation for the regeneration of plantlets of the pharmaceutically active Boesenbergia rotunda (L.) Mansf. Kulturpfl, formerly Boesenbergia/Kaempferia pandurata (Schult), to ensure a superior and consistent supply of materials for commercialization purposes. In this study, a single-medium formulation of Murashige and Skoog (MS) supplemented with 13.54μM 2,4-dichlorophenoxyacetic acid (2,4-D) was found to be the only medium out of eight formulations to promote the complete somatic embryogenesis process for the culture of B. rotunda (L.). Callus cultures were initiated from a total of 280 explants of rhizome meristem. The percentage of cultures forming embryogenic callus was 23.3 ±4.3% on this MS medium augmented by 13.54μM 2,4-D. The best plantlet regeneration rate was attained from the first subcultured callus with a mean of 6.6±0.1 plantlets per 1 cm diameter aggregate of callus. Somatic embryogenesis characteristic of monocots was evident from histological studies. The regenerated plantlets have been successfully established in soil.     

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.     

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.     

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.      

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.