A proof of the DBRF-MEGN method, an algorithm for deducing minimum equivalent gene networks

Background

We previously developed the DBRF-MEGN (difference-based regulation finding-minimum equivalent gene network) method, which deduces the most parsimonious signed directed graphs (SDGs) consistent with expression profiles of single-gene deletion mutants. However, until the present study, we have not presented the details of the method's algorithm or a proof of the algorithm.

Results

We describe in detail the algorithm of the DBRF-MEGN method and prove that the algorithm deduces all of the exact solutions of the most parsimonious SDGs consistent with expression profiles of gene deletion mutants.

Conclusions

The DBRF-MEGN method provides all of the exact solutions of the most parsimonious SDGs consistent with expression profiles of gene deletion mutants.     

A special issue on DNA damage response and genome stability

A novel, cancer-fighting function was recently discovered for Smad ubiquitination regulatory factor 2 (Smurf2).     

Genome-scale identification of Caenorhabditis elegans regulatory elements by tiling-array mapping of DNase I hypersensitive sites

Background

High throughput techniques have generated a huge set of biological data, which are deposited in various databases. Efficient exploitation of these databases is often hampered by a lack of appropriate tools, which allow easy and reliable identification of genes that miss functional characterization but are correlated with specific biological conditions (e.g. organotypic expression).

Results

We have developed a simple algorithm (DGSA = Database-dependent Gene Selection and Analysis) to identify genes with unknown functions involved in organ development concentrating on the heart. Using our approach, we identified a large number of yet uncharacterized genes, which are expressed during heart development. An initial functional characterization of genes by loss-of-function analysis employing morpholino injections into zebrafish embryos disclosed severe developmental defects indicating a decisive function of selected genes for developmental processes.

Conclusion

We conclude that DGSA is a versatile tool for database mining allowing efficient selection of uncharacterized genes for functional analysis.     

Multiscale Hy3S: Hybrid stochastic simulation for supercomputers

Background

DNA copy number alterations are one of the main characteristics of the cancer cell karyotype and can contribute to the complex phenotype of these cells. These alterations can lead to gains in cellular oncogenes as well as losses in tumor suppressor genes and can span small intervals as well as involve entire chromosomes. The ability to accurately detect these changes is central to understanding how they impact the biology of the cell.

Results

We describe a novel algorithm called CARAT (Copy Number Analysis with Regression And Tree) that uses probe intensity information to infer copy number in an allele-specific manner from high density DNA oligonuceotide arrays designed to genotype over 100, 000 SNPs. Total and allele-specific copy number estimations using CARAT are independently evaluated for a subset of SNPs using quantitative PCR and allelic TaqMan reactions with several human breast cancer cell lines. The sensitivity and specificity of the algorithm are characterized using DNA samples containing differing numbers of X chromosomes as well as a test set of normal individuals. Results from the algorithm show a high degree of agreement with results from independent verification methods.

Conclusion

Overall, CARAT automatically detects regions with copy number variations and assigns a significance score to each alteration as well as generating allele-specific output. When coupled with SNP genotype calls from the same array, CARAT provides additional detail into the structure of genome wide alterations that can contribute to allelic imbalance.     

Absence of carious lesions at margins of glass-ionomer cement and amalgam restorations: An update of systematic review evidence

Background

A recent analysis of protein sequences deposited in the NCBI RefSeq database indicates that ~8.5 million protein sequences are encoded in prokaryotic and eukaryotic genomes, where ~30% are explicitly annotated as "hypothetical" or "uncharacterized" protein. Our Comparison of Protein Active-Site Structures (CPASS v.2) database and software compares the sequence and structural characteristics of experimentally determined ligand binding sites to infer a functional relationship in the absence of global sequence or structure similarity. CPASS is an important component of our Functional Annotation Screening Technology by NMR (FAST-NMR) protocol and has been successfully applied to aid the annotation of a number of proteins of unknown function.

Findings

We report a major upgrade to our CPASS software and database that significantly improves its broad utility. CPASS v.2 is designed with a layered architecture to increase flexibility and portability that also enables job distribution over the Open Science Grid (OSG) to increase speed. Similarly, the CPASS interface was enhanced to provide more user flexibility in submitting a CPASS query. CPASS v.2 now allows for both automatic and manual definition of ligand-binding sites and permits pair-wise, one versus all, one versus list, or list versus list comparisons. Solvent accessible surface area, ligand root-mean square difference, and Cβ distances have been incorporated into the CPASS similarity function to improve the quality of the results. The CPASS database has also been updated.

Conclusions

CPASS v.2 is more than an order of magnitude faster than the original implementation, and allows for multiple simultaneous job submissions. Similarly, the CPASS database of ligand-defined binding sites has increased in size by ~ 38%, dramatically increasing the likelihood of a positive search result. The modification to the CPASS similarity function is effective in reducing CPASS similarity scores for false positives by ~30%, while leaving true positives unaffected. Importantly, receiver operating characteristics (ROC) curves demonstrate the high correlation between CPASS similarity scores and an accurate functional assignment. As indicated by distribution curves, scores ≥ 30% infer a functional similarity. Software URL: http://cpass.unl.edu.     

Towards a noninvasive anatomical and functional diagnostic work-up of patients with suspected coronary artery disease

Combining multidetector computed tomography and cardiovascular magnetic resonance imaging provides the clinician a strategy to comprehensively evaluate coronary morphology and function noninvasively. In the MARCC trial (Magnetic Resonance and CT in suspected CAD) a new noninvasive diagnostic work-up for patients with suspected coronary artery disease will be developed, involving the sequential use of both imaging techniques. (Neth Heart J 2010;18:270-3.)     

Osteopathic manipulative treatment and its relationship to autonomic nervous system activity as demonstrated by heart rate variability: a repeated measures study

Background

Osteopathic manipulative treatment (OMT) and ultrasound physical therapy (UPT) are commonly used for chronic low back pain. Although there is evidence from a systematic review and meta-analysis that OMT generally reduces low back pain, there are no large clinical trials that specifically assess OMT efficacy in chronic low back pain. Similarly, there is a lack of evidence involving UPT for chronic low back pain.

Methods

The OSTEOPAThic Health outcomes In Chronic low back pain (OSTEOPATHIC) Trial is a Phase III randomized controlled trial that seeks to study 488 subjects between August 2006 and June 2010. It uses a 2 × 2 factorial design to independently assess the efficacy of OMT and UPT for chronic low back pain. The primary outcome is a visual analogue scale score for pain. Secondary outcomes include back-specific functioning, generic health, work disability, and satisfaction with back care.

Conclusion

This randomized controlled trial will potentially be the largest involving OMT. It will provide long awaited data on the efficacy of OMT and UPT for chronic low back pain.

Trial registration

http://www.clinicaltrials.gov, NCT00315120     

Somatic embryogenesis of tissue cultures of Papaver somniferum and Papaver orientale and its relationship to alkaloid and lipid metabolism

Transfer from complete to 2,4-D free Gamborg's B5-medium efficiently induced somatic embryogenesis in Papaver tissue cultures (P. somniferum and P. orientale). Embryogenesis was preceded by a strong temporary accumulation of triacylglycerols. In both tissue cultures large amounts of sanguinarine type alkaloids were present, which disappeared during regeneration in the P. orientale cultures but persisted in the P. somniferum cultures. In the P. somniferum cultures protopine and morphine type alkaloids (morphine, codeine, thebaine) appeared about 45 days after exchanging the medium. Thebaine was the main alkaloid in the P. somniferum embryoids accumulating up to 0.2 % of dry weight.     

Genome skimming identifies polymorphism in tern populations and species

Background

Large-scale protein structure alignment, an indispensable tool to structural bioinformatics, poses a tremendous challenge on computational resources. To ensure structure alignment accuracy and efficiency, efforts have been made to parallelize traditional alignment algorithms in grid environments. However, these solutions are costly and of limited accessibility. Others trade alignment quality for speedup by using high-level characteristics of structure fragments for structure comparisons.

Findings

We present ppsAlign, a parallel protein structure Alignment framework designed and optimized to exploit the parallelism of Graphics Processing Units (GPUs). As a general-purpose GPU platform, ppsAlign could take many concurrent methods, such as TM-align and Fr-TM-align, into the parallelized algorithm design. We evaluated ppsAlign on an NVIDIA Tesla C2050 GPU card, and compared it with existing software solutions running on an AMD dual-core CPU. We observed a 36-fold speedup over TM-align, a 65-fold speedup over Fr-TM-align, and a 40-fold speedup over MAMMOTH.

Conclusions

ppsAlign is a high-performance protein structure alignment tool designed to tackle the computational complexity issues from protein structural data. The solution presented in this paper allows large-scale structure comparisons to be performed using massive parallel computing power of GPU.     

Inactivation of aphidicolin by plant cells

Plant cells are endowed with an aphidicolin inactivating activity. Data on cultured cells show that the rate of inactivation depends on the cell type, Daucus carota cells being the most effective among the other tested materials (Oryza sativa and Nicotiana plumbaginifolia). Also germinating seedling of Haplopappus gracilis and of Citrullus vulgaris inactivate aphidicolin. Inactivation, which may lead to unexpected results when a prolonged incubation with the drug is required, as in the case of the induction of synchrony of the cell cycle by aphidicolin, can be controlled by appropriately choosing the experimental conditions.     

Indole alkaloids from cell suspension cultures of Tabernaemontana divaricata and Tabernanthe iboga

From cell suspension cultures of Tabernaemontana divaricata and Tabernanthe iboga grown under standard conditions, six monoterpenoid indole alkaloids have been isolated and identified. T. divaricata synthesized apparicine, catharanthine, coronaridine, conoflorine, tubotaiwine and vinervine, whereas T. iboga produced tubotaiwine and conoflorine. Both cultures are a reasonable source for conoflorine, which is expected to be a good candidate for studying the mechanism of Aspidosperma type alkaloid formation at the cell-free level.     

Large-scale preparation of active caspase-3 in E. coli by designing its thrombin-activatable precursors

Background

Combination of CHD (chromo-helicase-DNA binding protein)-specific polymerase chain reaction (PCR) with electrophoresis (PCR/electrophoresis) is the most common avian molecular sexing technique but it is lab-intensive and gel-required. Gender determination often fails when the difference in length between the PCR products of CHD-Z and CHD-W genes is too short to be resolved.

Results

Here, we are the first to introduce a PCR-melting curve analysis (PCR/MCA) to identify the gender of birds by genomic DNA, which is gel-free, quick, and inexpensive. Spilornis cheela hoya (S. c. hoya) and Pycnonotus sinensis (P. sinensis) were used to illustrate this novel molecular sexing technique. The difference in the length of CHD genes in S. c. hoya and P. sinensis is 13-, and 52-bp, respectively. Using Griffiths' P2/P8 primers, molecular sexing failed both in PCR/electrophoresis of S. c. hoya and in PCR/MCA of S. c. hoya and P. sinensis. In contrast, we redesigned sex-specific primers to yield 185- and 112-bp PCR products for the CHD-Z and CHD-W genes of S. c. hoya, respectively, using PCR/MCA. Using this specific primer set, at least 13 samples of S. c. hoya were examined simultaneously and the Tm peaks of CHD-Z and CHD-W PCR products were distinguished.

Conclusion

In this study, we introduced a high-throughput avian molecular sexing technique and successfully applied it to two species. This new method holds a great potential for use in high throughput sexing of other avian species, as well.     

Isolation and characterization of DNA sequences from Triticum aestivum which function as origins of replication in Saccharomyces cerevisiae

Recombinant YIp5 plasmids with the DNA from Triticum aestivum are capable of autonomous replication in Saccharomyces cerevisiae. The URA transformants are unstable without selection pressure, and transformation of yeast cells with these plasmids occurs at high frequency. The cloned sequences were characterized and analyzed to state their belonging to Triticum tribe.     

The COG database: an updated version includes eukaryotes

Background

The availability of multiple, essentially complete genome sequences of prokaryotes and eukaryotes spurred both the demand and the opportunity for the construction of an evolutionary classification of genes from these genomes. Such a classification system based on orthologous relationships between genes appears to be a natural framework for comparative genomics and should facilitate both functional annotation of genomes and large-scale evolutionary studies.

Results

We describe here a major update of the previously developed system for delineation of Clusters of Orthologous Groups of proteins (COGs) from the sequenced genomes of prokaryotes and unicellular eukaryotes and the construction of clusters of predicted orthologs for 7 eukaryotic genomes, which we named KOGs after eukaryotic orthologous groups. The COG collection currently consists of 138,458 proteins, which form 4873 COGs and comprise 75% of the 185,505 (predicted) proteins encoded in 66 genomes of unicellular organisms. The eukaryotic orthologous groups (KOGs) include proteins from 7 eukaryotic genomes: three animals (the nematode Caenorhabditis elegans, the fruit fly Drosophila melanogaster and Homo sapiens), one plant, Arabidopsis thaliana, two fungi (Saccharomyces cerevisiae and Schizosaccharomyces pombe), and the intracellular microsporidian parasite Encephalitozoon cuniculi. The current KOG set consists of 4852 clusters of orthologs, which include 59,838 proteins, or ~54% of the analyzed eukaryotic 110,655 gene products. Compared to the coverage of the prokaryotic genomes with COGs, a considerably smaller fraction of eukaryotic genes could be included into the KOGs; addition of new eukaryotic genomes is expected to result in substantial increase in the coverage of eukaryotic genomes with KOGs. Examination of the phyletic patterns of KOGs reveals a conserved core represented in all analyzed species and consisting of ~20% of the KOG set. This conserved portion of the KOG set is much greater than the ubiquitous portion of the COG set (~1% of the COGs). In part, this difference is probably due to the small number of included eukaryotic genomes, but it could also reflect the relative compactness of eukaryotes as a clade and the greater evolutionary stability of eukaryotic genomes.

Conclusion

The updated collection of orthologous protein sets for prokaryotes and eukaryotes is expected to be a useful platform for functional annotation of newly sequenced genomes, including those of complex eukaryotes, and genome-wide evolutionary studies.     

6-OHDA generated ROS induces DNA damage and p53- and PUMA-dependent cell death

Background

Parkinson's disease (PD) is characterized by the selective loss of dopaminergic neurons in the substantia nigra (SN), resulting in tremor, rigidity, and bradykinesia. Although the etiology is unknown, insight into the disease process comes from the dopamine (DA) derivative, 6-hydroxydopamine (6-OHDA), which produces PD-like symptoms. Studies show that 6-OHDA activates stress pathways, such as the unfolded protein response (UPR), triggers mitochondrial release of cytochrome-c, and activates caspases, such as caspase-3. Because the BH3-only protein, Puma (p53-upregulated mediator of apoptosis), is activated in response to UPR, it is thought to be a link between cell stress and apoptosis.

Results

To test the hypothesis that Puma serves such a role in 6-OHDA-mediated cell death, we compared the response of dopaminergic neurons from wild-type and Puma-null mice to 6-OHDA. Results indicate that Puma is required for 6-OHDA-induced cell death, in primary dissociated midbrain cultures as well as in vivo. In these cultures, 6-OHDA-induced DNA damage and p53 were required for 6-OHDA-induced cell death. In contrast, while 6-OHDA led to upregulation of UPR markers, loss of ATF3 did not protect against 6-OHDA.

Conclusions

Together, our results indicate that 6-OHDA-induced upregulation of Puma and cell death are independent of UPR. Instead, p53 and DNA damage repair pathways mediate 6-OHDA-induced toxicity.