ELECTROPHORETIC AND MORPHOMETRIC CHARACTERS IN POPULATION DIFERENTIATION OF THE PEARL OYSTER, PINCTADA RADIATA (LEACH), FROM AROUND BAHRAIN

Variation at the leucine aminopeptidase (Lap), glucose phosphateisomerase (Gpi) and tetrazolium oxidase (To) loci was investigatedin samples of three populations, Al-Mayana (MAY), Shigita (SH)and Mina Salman (MS), of Pinctada radiata from pearl oysterbeds around Bahrain. The To locus was monomor-phic. SignificantLap and Gpi heterozygote deficiencies were evident and it issuggested that these were generated by selection. The MS population,to the East of Bahrain, differed significantly in Gpi allelefrequencies from both Northern populations (MAY, SH) and Nei'sgenetic identity indicates a close relationship between theNorthern populations. Measurements of shell morphometrics were used both as ratiosof one dimension to another, and as regressions of one dimensionon another to examine relatedness between populations. Boththese mor-phometric approaches gave different results from eachother and also differed from the electrophoretic data. It isconcluded that estimates of relatedness in pearl oysters basedon electrophoretic data will be more reliable than those basedon shell shape. (Received 20 November 1990; accepted 12 April 1991)     

Gut evacuation rates and ingestion rates of Eudiaptomus graciloides measured by means of the gut fluorescence method

Journal of Plankton Research, 8, 973–983, 1986 FIg. 2. Time-dependent changes in the gut content (percentageof initial ng pigment) of E. gro.ciloides at different temperaturesunder simultaneous feeding. Fig. 4. The relationship between instantaneous evacuation rateand temperature of E. graciloides. The regresston equation forfeeding animals: y = 0.0044 e(0.141 ) (r² = 0.90). For comparisonthe results of non-feeding animals are indicated with open circles.     

Unequal group variances in microarray data analyses

Motivation: In searching for differentially expressed (DE) genesin microarray data, we often observe a fraction of the genesto have unequal variability between groups. This is not an issuein large samples, where a valid test exists that uses individualvariances separately. The problem arises in the small-samplesetting, where the approximately valid Welch test lacks sensitivity,while the more sensitive moderated t-test assumes equal variance. Methods: We introduce a moderated Welch test (MWT) that allowsunequal variance between groups. It is based on (i) weightingof pooled and unpooled standard errors and (ii) improved estimationof the gene-level variance that exploits the information fromacross the genes. Results: When a non-trivial proportion of genes has unequalvariability, false discovery rate (FDR) estimates based on thestandard t and moderated t-tests are often too optimistic, whilethe standard Welch test has low sensitivity. The MWT is shownto (i) perform better than the standard t, the standard Welchand the moderated t-tests when the variances are unequal betweengroups and (ii) perform similarly to the moderated t, and betterthan the standard t and Welch tests when the group variancesare equal. These results mean that MWT is more reliable thanother existing tests over wider range of data conditions. Availability: R package to perform MWT is available at http://www.meb.ki.se/~yudpaw Contact: yudi.pawitan{at}ki.se Supplementary information: Supplementary data are availableat Bioinformatics online. Associate Editor: Martin Bishop     

Model-based deconvolution of genome-wide DNA binding

Motivation: Chromatin immunoprecipitation followed by hybridizationto a genomic tiling microarray (ChIP-chip) is a routinely usedprotocol for localizing the genomic targets of DNA-binding proteins.The resolution to which binding sites in this assay can be identifiedis commonly considered to be limited by two factors: (1) theresolution at which the genomic targets are tiled in the microarrayand (2) the large and variable lengths of the immunoprecipitatedDNA fragments. Results: We have developed a generative model of binding sitesin ChIP-chip data and an approach, MeDiChI, for efficientlyand robustly learning that model from diverse data sets. Wehave evaluated MeDiChI's performance using simulated data, aswell as on several diverse ChIP-chip data sets collected onwidely different tiling array platforms for two different organisms(Saccharomyces cerevisiae and Halobacterium salinarium NRC-1).We find that MeDiChI accurately predicts binding locations toa resolution greater than that of the probe spacing, even foroverlapping peaks, and can increase the effective resolutionof tiling array data by a factor of 5x or better. Moreover,the method's performance on simulated data provides insightsinto effectively optimizing the experimental design for increasedbinding site localization accuracy and efficacy. Availability: MeDiChI is available as an open-source R package,including all data, from http://baliga.systemsbiology.net/medichi. Contact: dreiss{at}systemsbiology.org Supplementary information: Supplementary data are availableat Bioinformatics online. Associate Editor: Martin Bishop     

Robust accurate identification of peptides (RAId): deciphering MS2 data using a structured library search with de novo based statistics

Motivation: The key to MS -based proteomics is peptide sequencing.The major challenge in peptide sequencing, whether library searchor de novo, is to better infer statistical significance andbetter attain noise reduction. Since the noise in a spectrumdepends on experimental conditions, the instrument used andmany other factors, it cannot be predicted even if the peptidesequence is known. The characteristics of the noise can onlybe uncovered once a spectrum is given. We wish to overcome suchissues. Results: We designed RAId to identify peptides from their associatedtandem mass spectrometry data. RAId performs a novel de novosequencing followed by a search in a peptide library that wecreated. Through de novo sequencing, we establish the spectrum-specificbackground score statistics for the library search. When thedatabase search fails to return significant hits, the top-rankingde novo sequences become potential candidates for new peptidesthat are not yet in the database. The use of spectrum-specificbackground statistics seems to enable RAId to perform well evenwhen the spectral quality is marginal. Other important featuresof RAId include its potential in de novo sequencing alone andthe ease of incorporating post-translational modifications. Availability: Programs implementing the methods described areavailable from the authors on request. Contact: yyu{at}ncbi.nlm.nih.gov Supplementary information: ftp://ftp.ncbi.nih.gov/pub/yyu/Proteomics/MSMS/RAId/MSMS_bioinfo_supp.pdf     

CORRIGENDUM

EAMES, F. E., 1968. New name for a Pakistan Eocene Turbonilla.Proc. malac. Soc. Lond. 38, 167. Line 4 : For J. Conch., Lond. p. 95, read Ann. Mag. Nat. Hist.(6) (6), 95.     

Contribution of Crenarchaeota and Euryarchaeota to the prokaryotic plankton in the coastal northwestern Black Sea

Journal of Plankton Research, 29,     

Observation of multiple folding pathways of beta-hairpin trpzip2 from independent continuous folding trajectories

Motivation: After 10-year investigations, the folding mechanismsof β-hairpins are still under debate. Experiments stronglysupport zip-out pathway, while most simulations prefer the hydrophobiccollapse model (including middle-out and zip-in pathways). Inthis article, we show that all pathways can occur during thefolding of β-hairpins but with different probabilities.The zip-out pathway is the most probable one. This is in agreementwith the experimental results. We came to our conclusions by38 100-ns room-temperature all-atom molecular dynamics simulationsof the β-hairpin trpzip2. Our results may help to clarifythe inconsistencies in the current pictures of β-hairpinfolding mechanisms. Contact: yxiao{at}mail.hust.edu.cn Supplementary information: Supplementary data are availableat Bioinformatics online. Associate Editor: Anna Tramontano     

'ANOVA-simultaneous component analysis (ASCA): a new tool for analyzing designed metabolomics data'

Smilde et al. Bioinformatics (2005), 21(13); 3043–3048 The above paper by Smilde et al. inappropriately quotes results     

PETALIFERA HABEI, A NEW SPECIES FROM JAPAN ADDENDUM

Since Petalifera habei was described in the Proceedings, 34,1, 12–18, April 1960, I have received from Dr. K. Babaa short account of the same species in Publ. Seto. mar. biol.Lab. 7, 3, 337–338. December, 1959, under the title "Thegenus Petalifera and a new species, P. ramosa, from Japan".I was unaware that Dr. Baba intended to describe it. As hispaper antedates mine, his name, Petalifera ramosa, must replaceP. Habei.     

A new algorithm for cluster analysis of genomic methylation: the Helicobacter pylori case

Motivation: The genomic methylation analysis is useful to typebacteria that have a high number of expressed type II methyltransferases.Methyltransferases are usually committed to Restriction andModification (R-M) systems, in which the restriction endonucleaseimposes high pressure on the expression of the cognate methyltransferasethat hinder R-M system loss. Conventional cluster methods donot reflect this tendency. An algorithm was developed for dendrogramconstruction reflecting the propensity for conservation of R-MType II systems. Results: The new algorithm was applied to 52 Helicobacter pyloristrains from different geographical regions and compared withconventional clustering methods. The algorithm works by firstgrouping strains that share a common minimum set of R-M systemsand gradually adds strains according to the number of the R-Msystems acquired. Dendrograms revealed a cluster of Africanstrains, which suggest that R-M systems are present in H.pylorigenome since its human host migrates from Africa. Availability: The software files are available at http://www.ff.ul.pt/paginas/jvitor/Bioinformatics/MCRM_algorithm.zip Contact: filipavale{at}fe.ucp.pt Supplementary information: Supplementary data are availableat Bioinformatics online. Associate Editor: Martin Bishop     

Effect of Irradiance on Chlorophyll Estimation with the Minolta SPAD-502 Leaf Chlorophyll Meter

Leaf chlorophyll content may be used as an indirect indicatorof crop nitrogen status. Chlorophyll meter values (SPAD values)taken with the Minolta SPAD-502 chlorophyll meter in the shadeplantOxalis acetosellaL. and in winter wheat (Triticum aestivumL.)varied by 15 and 8%, respectively, with variation in irradiance.The lowest SPAD-values were measured at high irradiance. Duringa natural night-day-night cycle SPAD values for winter wheatwere lowest in the middle of the day, highest at low irradianceat dusk and dawn and intermediate in darkness before dawn andafter dusk. The results indicate that irradiance during measurementshould be considered when using the Minolta SPAD-502 chlorophyllmeter for the estimation of crop N-status.Copyright 1998 Annalsof Botany Company Chlorophyll meter, nitrogen, irradiance,Oxalis acetosellaL.,Triticum aestivumL., winter wheat.     

Interaction Between Cryopreservation, Rewarming Rate and Seed Humidification on the Germination of Two Spanish Endemic Species

The effects of humidification, storage in liquid nitrogen (1or 30 d) and rewarming rate on seed germination were studiedin two Spanish endemics. Humidification resulted in higher germinationpercentages only in the species with hard covers, especiallyin slowly rewarmed seeds. In an experiment lasting 21 weeks,seeds stored in liquid nitrogen were removed for 10 min eachweek to mimic the withdrawal of samples from a seed bank; thishad no effect on germination.Copyright 1998 Annals of BotanyCompany Centaurea hyssopifolia,Limonium dichotomum, cryopreservation, cypsela, endemics, germination, humidification, seeds.     

AIMIE: a web-based environment for detection and interpretation of significant sequence motifs in prokaryotic genomes

Motivation: Genomes contain biologically significant informationthat extends beyond that encoded in genes. Some of this informationrelates to various short dispersed repeats distributed throughoutthe genome. The goal of this work was to combine tools for detectionof statistically significant dispersed repeats in DNA sequenceswith tools to aid development of hypotheses regarding theirpossible physiological functions in an easy-to-use web-basedenvironment. Results: Ab Initio Motif Identification Environment (AIMIE)was designed to facilitate investigations of dispersed sequencemotifs in prokaryotic genomes. We used AIMIE to analyze theEscherichia coli and Haemophilus influenzae genomes in orderto demonstrate the utility of the new environment. AIMIE detectedrepeated extragenic palindrome (REP) elements, CRISPR repeats,uptake signal sequences, intergenic dyad sequences and severalother over-represented sequence motifs. Distributional patternsof these motifs were analyzed using the tools included in AIMIE. Availability: AIMIE and the related software can be accessedat our web site http://www.cmbl.uga.edu/software.html. Contact: mrazek{at}uga.edu Associate Editor: Alex Bateman     

Mining SARS-CoV protease cleavage data using non-orthogonal decision trees: a novel method for decisive template selection

Motivation: Although the outbreak of the severe acute respiratorysyndrome (SARS) is currently over, it is expected that it willreturn to attack human beings. A critical challenge to scientistsfrom various disciplines worldwide is to study the specificityof cleavage activity of SARS-related coronavirus (SARS-CoV)and use the knowledge obtained from the study for effectiveinhibitor design to fight the disease. The most commonly usedinductive programming methods for knowledge discovery from dataassume that the elements of input patterns are orthogonal toeach other. Suppose a sub-sequence is denoted as P₂P₁P_1'P_2',the conventional inductive programming method may result ina rule like ‘if P₁ = Q, then the sub-sequence is cleaved,otherwise non-cleaved’. If the site P₁ is not orthogonalto the others (for instance, P₂, P_1' and P_2'), the predictionpower of these kind of rules may be limited. Therefore thisstudy is aimed at developing a novel method for constructingnon-orthogonal decision trees for mining protease data. Result: Eighteen sequences of coronavirus polyprotein were downloadedfrom NCBI (http://www.ncbi.nlm.nih.gov). Among these sequences,252 cleavage sites were experimentally determined. These sequenceswere scanned using a sliding window with size k to generateabout 50 000 k-mer sub-sequences (for short, k-mers). The valueof k varies from 4 to 12 with a gap of two. The bio-basis functionproposed by Thomson et al. is used to transform the k-mers toa high-dimensional numerical space on which an inductive programmingmethod is applied for the purpose of deriving a decision treefor decision-making. The process of this transform is referredto as a bio-mapping. The constructed decision trees select about10 out of 50 000 k-mers. This small set of selected k-mers isregarded as a set of decisive templates. By doing so, non-orthogonaldecision trees are constructed using the selected templatesand the prediction accuracy is significantly improved. Availability: The program for bio-mapping can be obtained byrequest to the author. Contact: z.r.yang{at}exeter.ac.uk     

GOTreePlus: an interactive gene ontology browser

Summary: We developed an interactive gene ontology (GO) browsernamed GOTreePlus that superimposes annotation information overGO structures. It can facilitate the identification of importantGO terms through interactive visualization of them in the GOstructure. The interactive pie chart summarizing an annotationdistribution for a selected GO term provides users with a succinctcontext-sensitive overview of their experimental results. Wetested our GOTreePlus using a proteome profiling dataset obtainedon differentiation of retinal pigment epithelial cells where399 proteins were quantified. Availability: http://bioinformatics.cnmcresearch.org/GOTreePlus/ Contact: jseo{at}cnmcresearch.org Associate Editor: John Quackenbush     

A fuzzy guided genetic algorithm for operon prediction

Motivation: The operon structure of the prokaryotic genome isa critical input for the reconstruction of regulatory networksat the whole genome level. As experimental methods for the detectionof operons are difficult and time-consuming, efforts are beingput into developing computational methods that can use availablebiological information to predict operons. Method: A genetic algorithm is developed to evolve a startingpopulation of putative operon maps of the genome into progressivelybetter predictions. Fuzzy scoring functions based on multiplecriteria are used for assessing the ‘fitness’ ofthe newly evolved operon maps and guiding their evolution. Results: The algorithm organizes the whole genome into operons.The fuzzy guided genetic algorithm-based approach makes it possibleto use diverse biological information like genome sequence data,functional annotations and conservation across multiple genomes,to guide the organization process. This approach does not requireany prior training with experimental operons. The predictionsfrom this algorithm for Escherchia coli K12 and Bacillus subtilisare evaluated against experimentally discovered operons forthese organisms. The accuracy of the method is evaluated usingan ROC (receiver operating characteristic) analysis. The areaunder the ROC curve is around 0.9, which indicates excellentaccuracy. Contact: roschen_csir{at}rediffmail.com     

3D-Garden: a system for modelling protein-protein complexes based on conformational refinement of ensembles generated with the marching cubes algorithm

Motivation: Reliable structural modelling of protein–proteincomplexes has widespread application, from drug design to advancingour knowledge of protein interactions and function. This workaddresses three important issues in protein–protein docking:implementing backbone flexibility, incorporating prior indicationsfrom experiment and bioinformatics, and providing public accessvia a server. 3D-Garden (Global And Restrained Docking ExplorationNexus), our benchmarked and server-ready flexible docking system,allows sophisticated programming of surface patches by the uservia a facet representation of the interactors’ molecularsurfaces (generated with the marching cubes algorithm). Flexibilityis implemented as a weighted exhaustive conformer search foreach clashing pair of molecular branches in a set of 5000 modelsfiltered from around 340 000 initially. Results: In a non-global assessment, carried out strictly accordingto the protocols for number of models considered and model qualityof the Critical Assessment of Protein Interactions (CAPRI) experiment,over the widely-used Benchmark 2.0 of 84 complexes, 3D-Gardenidentifies a set of ten models containing an acceptable or bettermodel in 29/45 test cases, including one with large conformationalchange. In 19/45 cases an acceptable or better model is rankedfirst or second out of 340 000 candidates. Availability: http://www.sbg.bio.ic.ac.uk/3dgarden (server) Contact: v.lesk{at}ic.ac.uk Supplementary information: Supplementary data are availableat Bioinformatics online. Associate Editor: Burkhard Rost