Multiclass cancer classification and biomarker discovery using GA-based algorithms

MOTIVATION: The development of microarray-based high-throughput gene profiling has led to the hope that this technology could provide an efficient and accurate means of diagnosing and classifying tumors, as well as predicting prognoses and effective treatments. However, the large amount of data generated by microarrays requires effective reduction of discriminant gene features into reliable sets of tumor biomarkers for such multiclass tumor discrimination. The availability of reliable sets of biomarkers, especially serum biomarkers, should have a major impact on our understanding and treatment of cancer. RESULTS: We have combined genetic algorithm (GA) and all paired (AP) support vector machine (SVM) methods for multiclass cancer categorization. Predictive features can be automatically determined through iterative GA/SVM, leading to very compact sets of non-redundant cancer-relevant genes with the best classification performance reported to date. Interestingly, these different classifier sets harbor only modest overlapping gene features but have similar levels of accuracy in leave-one-out cross-validations (LOOCV). Further characterization of these optimal tumor discriminant features, including the use of nearest shrunken centroids (NSC), analysis of annotations and literature text mining, reveals previously unappreciated tumor subclasses and a series of genes that could be used as cancer biomarkers. With this approach, we believe that microarray-based multiclass molecular analysis can be an effective tool for cancer biomarker discovery and subsequent molecular cancer diagnosis.     

PII signal transduction proteins: sensors of alpha-ketoglutarate that regulate nitrogen metabolism

PII proteins are small homotrimeric signal transduction proteins that regulate the activities of metabolic enzymes and permeases, and control the activities of signal transduction enzymes. The protein family shows high conservation, with examples in eukaryota (plants and eukaryotic algae), archaea, and bacteria. This distribution indicates that PII is one of the most ancient signalling proteins known.     

Epidemiological evaluations of the efficacy of slow-released praziquantel-medicated bars for dogs in the prevention and control of cystic echinococcosis in man and animals

To assess the epidemiological efficacy of type SRP III slow-released praziquantel-medicated bars for dogs in the prevention and control of cystic echinococcosis in man and livestock, praziquantel-medicated bars were implanted subcutaneously in over 90% of dogs in villages in north Xinjiang, China, where cystic echinococcosis is highly endemic. After implantation, infection rate of Echinococcus granulosus in dogs, specific antibodies in children and prevalence of echinococcosis in one-year-old lambs were observed for 3 years. Coproantigen of E. granulosus was positive in 41.2% of the dogs at the start of experiment. In the second and third year after medicated-bar implantation, coproantigen was undetectable in any dogs examined, while 3.0% of dogs were positive at the end of the fourth year. The antibody positive rate in 7-year-old pupils, that was 41.2% before the experiment, declined gradually and it was 5.4% in the fourth year, while children in the non-intervention control area showed 30.6% seropositivity. The prevalence of hydatid disease in children aged 7–16 years also declined significantly. The prevalence of hydatidosis in lambs one year of age was 44.8% in the first year, dropping to 10.7% in the fourth year, while in the non-intervention control area the level of infection was 46.4%. These results demonstrated not only that the slow released praziquantel-medicated bars efficiently blocked reinfection in dogs at least for 2 years, but also the measure was effective in preventing transmission of cystic echinococcosis to both man and livestock.     

Characterization of a protein-based adhesive elastomer secreted by the Australian frog Notaden bennetti

When provoked, Notaden bennetti frogs secrete an exudate which rapidly forms a tacky elastic solid ("frog glue"). This protein-based material acts as a promiscuous pressure-sensitive adhesive that functions even in wet conditions. We conducted macroscopic tests in air to assess the tensile strength of moist glue (up to 78 +/- 8 kPa) and the shear strength of dry glue (1.7 +/- 0.3 MPa). We also performed nanomechanical measurements in water to determine the adhesion (1.9-7.2 nN or greater), resilience (43-56%), and elastic modulus (170-1035 kPa) of solid glue collected in different ways. Dry glue contains little carbohydrate and consists mainly of protein. The protein complement is rich in Gly (15.8 mol %), Pro (8.8 mol %), and Glu/Gln (14.1 mol %); it also contains some 4-hydroxyproline (4.6 mol %) but no 5-hydroxylysine or 3,4-dihydroxyphenylalanine (L-Dopa). Denaturing gel electrophoresis of the glue reveals a characteristic pattern of proteins spanning 13-400 kDa. The largest protein (Nb-1R, apparent molecular mass 350-500 kDa) is also the most abundant, and this protein appears to be the key structural component. The solid glue can be dissolved in dilute acids; raising the ionic strength causes the glue components to self-assemble spontaneously into a solid which resembles the starting material. We describe scattering studies on dissolved and solid glue and provide microscopy images of glue surfaces and sections, revealing a porous interior that is consistent with the high water content (85-90 wt %) of moist glue. In addition to compositional similarities with other biological adhesives and well-known elastomeric proteins, the circular dichroism spectrum of dissolved glue is almost identical to that for soluble elastin and electron and scanning probe microscopy images invite comparison with silk fibroins. Covalent cross-linking does not seem to be necessary for the glue to set.     

Further investigation on potassium-induced conformation transition of Nephila spidroin film with two-dimensional infrared correlation spectroscopy

We used two-dimensional (2D) correlation infrared spectroscopy to study further the potassium-induced conformation transition in Nephila spidroin films. It provided increased resolution and important new information on the sequence of events in the conformation transition process, showing that beta-sheet formed from the helical component before they formed from random coil. It also showed more evidence that formation of the 1691 cm(-1) (turn/bend) peak did not proceed with the same kinetics as the 1620 cm(-1) (antiparallel beta-sheet component) one, so we attribute the 1691 cm(-1) peak to turns which formed with different kinetics as the antiparallel beta-sheets. We present a single coherent and detailed hypothesis for the assembly and secondary structural transition of silk proteins in vivo and in vitro based on our findings and on evidence from other laboratories.     

Na,K-ATPase from mice lacking the gamma subunit (FXYD2) exhibits altered Na+ affinity and decreased thermal stability

The gamma subunit of the Na,K-ATPase, a 7-kDa single-span membrane protein, is a member of the FXYD gene family. Several FXYD proteins have been shown to bind to Na,K-ATPase and modulate its properties, and each FXYD protein appears to alter enzyme kinetics differently. Different results have sometimes been obtained with different experimental systems, however. To test for effects of gamma in a native tissue environment, mice lacking a functional gamma subunit gene (Fxyd2) were generated. These mice were viable and without observable pathology. Prior work in the mouse embryo showed that gamma is expressed at the blastocyst stage. However, there was no delay in blastocele formation, and the expected Mendelian ratios of offspring were obtained even with Fxyd2-/- dams. In adult Fxyd2-/- mouse kidney, splice variants of gamma that have different nephron segment-specific expression patterns were absent. Purified gamma-deficient renal Na,K-ATPase displayed higher apparent affinity for Na+ without significant change in apparent affinity for K+. Affinity for ATP, which was expected to be decreased, was instead slightly increased. The results suggest that regulation of Na+ sensitivity is a major functional role for this protein, whereas regulation of ATP affinity may be context-specific. Most importantly, this implies that gamma and other FXYD proteins have their effects by local and not global conformation change. Na,K-ATPase lacking the gamma subunit had increased thermal lability. Combined with other evidence that gamma participates in an early step of thermal denaturation, this indicates that FXYD proteins may play an important structural role in the enzyme complex.     

Crystal structure of human vacuolar protein sorting protein 29 reveals a phosphodiesterase/nuclease-like fold and two protein-protein interaction sites

Vacuolar protein sorting protein 29 (Vps29p), which is involved in retrograde trafficking from prevacuolar endosomes to the trans-Golgi network, performs its biological functions by participating in the formation of a "retromer complex." In human cells, this complex comprises four conserved proteins: hVps35p, hVps29p, hVps26p, and sorting nexin 1 protein (SNX1). Here, we report the crystal structure of hVps29p at 2.1 Angstroms resolution, the first three-dimensional structure of the retromer subunits. This novel structure adopts a four-layered alpha-beta-beta-alpha sandwich fold. hVps29p contains a metal-binding site that is very similar to the active sites of some proteins of the phosphodiesterase/nuclease protein family, indicating that hVps29p may carry out chemically similar functions. Structure and sequence conservation analysis suggests that hVps29p contains two protein-protein interaction sites. One site, which potentially serves as the interface between hVps29p and hVps35p, comprises 5 conserved hydrophobic and 8 hydrophilic residues. The other site is relatively more hydrophilic and may serve as a binding interface with hVps26p, SNX1, or other target proteins.     

Genetic mapping and comparative analysis of seven mutants related to seed fiber development in cotton

Mapping of genes that play major roles in cotton fiber development is an important step toward their cloning and manipulation, and provides a test of their relationships (if any) to agriculturally-important QTLs. Seven previously identified fiber mutants, four dominant (Li ₁, Li ₂, N ₁ and Fbl) and three recessive (n ₂, sma-4(h _a), and sma-4(fz)), were genetically mapped in six F₂ populations comprising 124 or more plants each. For those mutants previously assigned to chromosomes by using aneuploids or by linkage to other morphological markers, all map locations were concordant except n ₂, which mapped to the homoeolog of the chromosome previously reported. Three mutations with primary effects on fuzz fibers (N ₁, Fbl, n ₂) mapped near the likelihood peaks for QTLs that affected lint fiber productivity in the same populations, perhaps suggesting pleiotropic effects on both fiber types. However, only Li ₁ mapped within the likelihood interval for 191 previously detected lint fiber QTLs discovered in non-mutant crosses, suggesting that these mutations may occur in genes that played early roles in cotton fiber evolution, and for which new allelic variants are quickly eliminated from improved germplasm. A close positional association between sma-4(h _a ), two leaf and stem-borne trichome mutants (t ₁ , t ₂), and a gene previously implicated in fiber development, sucrose synthase, raises questions about the possibility that these genes may be functionally related. Increasing knowledge of the correspondence of the cotton and Arabidopsis genomes provides several avenues by which genetic dissection of cotton fiber development may be accelerated.     

Molecular analysis of the O-antigen gene cluster of Escherichia coli O86:B7 and characterization of the chain length determinant gene (wzz)

Escherichia coli O86:B7 has long been used as a model bacterial strain to study the generation of natural blood group antibody in humans, and it has been shown to possess high human blood B activity. The O-antigen structure of O86:B7 was solved recently in our laboratory. Comparison with the published structure of O86:H2 showed that both O86 subtypes shared the same O unit, yet each of the O antigens is polymerized from a different terminal sugar in a different glycosidic linkage. To determine the genetic basis for the O-antigen differences between the two O86 strains, we report the complete sequence of O86:B7 O-antigen gene cluster between galF and hisI, each gene was identified based on homology to other genes in the GenBank databases. Comparison of the two O86 O-antigen gene clusters revealed that the encoding regions between galF and gnd are identical, including wzy genes. However, deletion of the two wzy genes revealed that wzy in O86:B7 is responsible for the polymerization of the O antigen, while the deletion of wzy in O86:H2 has no effect on O-antigen biosynthesis. Therefore, we proposed that there must be another functional wzy gene outside the O86:H2 O-antigen gene cluster. Wzz proteins determine the degree of polymerization of the O antigen. When separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis, the lipopolysaccharide (LPS) of O86:B7 exhibited a modal distribution of LPS bands with relatively short O units attached to lipid A-core, which differs from the LPS pattern of O86:H2. We proved that the wzz genes are responsible for the different LPS patterns found in the two O86 subtypes, and we also showed that the very short type of LPS is responsible for the serum sensitivity of the O86:B7 strain.     

Digital biology: an emerging and promising discipline

This article examines the role of computation and quantitative methods in modern biomedical research to identify emerging scientific, technical, policy and organizational trends. It identifies common concerns and practices in the emerging community of computationally-oriented bio-scientists by reviewing a national symposium, Digital Biology: the Emerging Paradigm, held at the National Institutes of Health in Bethesda, Maryland, November 6th and 7th 2003. This meeting showed how biomedical computing promises scientific breakthroughs that will yield significant health benefits. Three key areas that define the emerging discipline of digital biology are: scientific data integration, multi-scale modeling and networked science. Each area faces unique technical challenges and information policy issues that must be addressed as the field matures. Here we summarize the emergent challenges and offer suggestions to academia, industry and government on how best to expand the role of computation in their scientific activities.