BioPAX support in CellDesigner

MOTIVATION: BioPAX is a standard language for representing and exchanging models of biological processes at the molecular and cellular levels. It is widely used by different pathway databases and genomics data analysis software. Currently, the primary source of BioPAX data is direct exports from the curated pathway databases. It is still uncommon for wet-lab biologists to share and exchange pathway knowledge using BioPAX. Instead, pathways are usually represented as informal diagrams in the literature. In order to encourage formal representation of pathways, we describe a software package that allows users to create pathway diagrams using CellDesigner, a user-friendly graphical pathway-editing tool and save the pathway data in BioPAX Level 3 format. AVAILABILITY: The plug-in is freely available and can be downloaded at ftp://ftp.pantherdb.org/CellDesigner/plugins/BioPAX/ CONTACT: huaiyumi@usc.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

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.     

Rice plant development: from zygote to spikelet

Rice is becoming a model plant in monocotyledons and a model cereal crop. For better understanding of the rice plant, it is essential to elucidate the developmental programs of the life cycle. To date, several attempts have been made in rice to categorize the developmental processes of some organs into substages. These studies are based exclusively on the morphological and anatomical viewpoints. Recent advancement in genetics and molecular biology has given us new aspects of developmental processes. In this review, we first describe the phasic development of the rice plant, and then describe in detail the developmental courses of major organs, leaf, root and spikelet, and specific organs/tissues. Also, for the facility of future studies, we propose a staging system for each organ.     

A novel cytochrome P450 is implicated in brassinosteroid biosynthesis via the characterization of a rice dwarf mutant, dwarf11, with reduced seed length

We have characterized a rice (Oryza sativa) dwarf mutant, dwarf11 (d11), that bears seeds of reduced length. To understand the mechanism by which seed length is regulated, the D11 gene was isolated by a map-based cloning method. The gene was found to encode a novel cytochrome P450 (CYP724B1), which showed homology to enzymes involved in brassinosteroid (BR) biosynthesis. The dwarf phenotype of d11 mutants was restored by the application of the brassinolide (BL). Compared with wild-type plants, the aberrant D11 mRNA accumulated at higher levels in d11 mutants and was dramatically reduced by treatment with BL, implying that the gene is feedback-regulated by BL. Precise determination of the defective step(s) in BR synthesis in d11 mutants proved intractable because of tissue specificity and the complex control of BR accumulation in plants. However, 6-deoxotyphasterol (6-DeoxoTY) and typhasterol (TY), but not any upstream intermediates before these compounds, effectively restored BR response in d11 mutants in a lamina joint bending assay. Multiple lines of evidence together suggest that the D11/CYP724B1 gene plays a role in BR synthesis and may be involved in the supply of 6-DeoxoTY and TY in the BR biosynthesis network in rice.     

Using process diagrams for the graphical representation of biological networks

With the increased interest in understanding biological networks, such as protein-protein interaction networks and gene regulatory networks, methods for representing and communicating such networks in both human- and machine-readable form have become increasingly important. Although there has been significant progress in machine-readable representation of networks, as exemplified by the Systems Biology Mark-up Language (SBML) (http://www.sbml.org) issues in human-readable representation have been largely ignored. This article discusses human-readable diagrammatic representations and proposes a set of notations that enhances the formality and richness of the information represented. The process diagram is a fully state transition-based diagram that can be translated into machine-readable forms such as SBML in a straightforward way. It is supported by CellDesigner, a diagrammatic network editing software (http://www.celldesigner.org/), and has been used to represent a variety of networks of various sizes (from only a few components to several hundred components).     

Evolutionary conservation of 5' upstream sequence of nine genes between human and great apes

Nucleotide sequences of nine 5' upstream gene regions for human, chimpanzee, gorilla, and orangutan were determined. We estimated nucleotide differences (d) for each region between human and great apes. The overall d was 0.027 (ranged from 0.004 to 0.052). Rates of nucleotide substitution were estimated by using d and divergence times of human, chimpanzee, gorilla, and orangutan. The overall rate of nucleotide substitution between human and other hominoids was estimated to be 0.52-0.85 x 10(-9). This rate in 5' upstream regions was lower than that of synonymous sites, suggesting that 5' upstream regions have evolved under some functional constraints. Because lower rates have been reported for coding sequences in primates compared to rodents, we also estimated the rate (1.17-1.76 x 10(-9)) of nucleotide substitutions for the corresponding 5' upstream regions in rodents (mouse/rat comparison). Thus the primate rate was lower than rodent rate also for the 5' upstream regions.     

<Emphasis Type="Italic">OsMADS22</Emphasis>, an <Emphasis Type="Italic">STMADS11</Emphasis>-like MADS-box gene of rice,is expressed in non-vegetative tissues and its ectopic expression induces spikelet meristem indeterminacy

We report the cDNA sequence and gene expression patterns of OsMADS22, a novel member of the STMADS11-like family of MADS-box genes, from rice. In contrast to previously reported STMADS11-like genes, whose expression is detected in vegetative tissues, OsMADS22 is mainly expressed during embryogenesis and flower development. In situ hybridization analysis revealed that OsMADS22 expression is localized in the L1 layer of embryos and in developing stamen primordia. Ectopic expression of OsMADS22 in transgenic rice plants resulted in aberrant floral morphogenesis, characterized by a disorganized palea, an elongated glume, and a two-floret spikelet. The results are discussed in terms of rice spikelet development and a novel non-vegetative role for a STMADS11-like gene.     

Expression cloning of a human cDNA restoring sphingomyelin synthesis and cell growth in sphingomyelin synthase-defective lymphoid cells

Sphingomyelin (SM) synthase has been assumed to be involved in both cell death and survival by regulating pro-apoptotic mediator ceramide and pro-survival mediator diacylglycerol. However, its precise functions are ambiguous due to the lack of molecular cloning of SM synthase gene(s). We isolated WR19L/Fas-SM(-) mouse lymphoid cells, which show a defect of SM at the plasma membrane due to the lack of SM synthase activity and resistance to cell death induced by an SM-directed cytolytic protein lysenin. WR19L/Fas-SM(-) cells were also highly susceptible to methyl-beta-cyclodextrin (MbetaCD) as compared with the WR19L/Fas-SM(+) cells, which are capable of SM synthesis. By expression cloning method using WR19L/Fas-SM(-) cells and MbetaCD-based selection, we have succeeded in cloning of a human cDNA responsible for SM synthase activity. The cDNA encodes a peptide of 413 amino acids named SMS1 (putative molecular mass, 48.6 kDa), which contains a sterile alpha motif domain near the N-terminal region and four predicted transmembrane domains. WR19L/Fas-SM(-) cells expressing SMS1 cDNA (WR19L/Fas-SMS1) restored the resistance against MbetaCD, the accumulation of SM at the plasma membrane, and SM synthesis by transferring phosphocholine from phosphatidylcholine to ceramide. Furthermore, WR19L/Fas-SMS1 cells, as well as WR19L/Fas-SM(-) cells supplemented with exogenous SM, restored cell growth ability in serum-free conditions, where the growth of WR19L/Fas-SM(-) cells was severely inhibited. The results suggest that SMS1 is responsible for SM synthase activity in mammalian cells and plays a critical role in cell growth of mouse lymphoid cells.     

Polymorphisms in the trace amine receptor 4 (TRAR4) gene on chromosome 6q23.2 are associated with susceptibility to schizophrenia

Several linkage studies across multiple population groups provide convergent support for a susceptibility locus for schizophrenia--and, more recently, for bipolar disorder--on chromosome 6q13-q26. We genotyped 192 European-ancestry and African American (AA) pedigrees with schizophrenia from samples that previously showed linkage evidence to 6q13-q26, focusing on the MOXD1-STX7-TRARs gene cluster at 6q23.2, which contains a number of prime candidate genes for schizophrenia. Thirty-one screening single-nucleotide polymorphisms (SNPs) were selected, providing a minimum coverage of at least 1 SNP/20 kb. The association observed with rs4305745 (P=.0014) within the TRAR4 (trace amine receptor 4) gene remained significant after correction for multiple testing. Evidence for association was proportionally stronger in the smaller AA sample. We performed database searches and sequenced genomic DNA in a 30-proband subsample to obtain a high-density map of 23 SNPs spanning 21.6 kb of this gene. Single-SNP analyses and also haplotype analyses revealed that rs4305745 and/or two other polymorphisms in perfect linkage disequilibrium (LD) with rs4305745 appear to be the most likely variants underlying the association of the TRAR4 region with schizophrenia. Comparative genomic analyses further revealed that rs4305745 and/or the associated polymorphisms in complete LD with rs4305745 could potentially affect gene expression. Moreover, RT-PCR studies of various human tissues, including brain, confirm that TRAR4 is preferentially expressed in those brain regions that have been implicated in the pathophysiology of schizophrenia. These data provide strong preliminary evidence that TRAR4 is a candidate gene for schizophrenia; replication is currently being attempted in additional clinical samples.