Better prediction of protein contact number using a support vector regression analysis of amino acid sequence

Background: Protein tertiary structure can be partly characterized via each amino acid's contact number measuring how residues are spatially arranged. The contact number of a residue in a folded protein is a measure of its exposure to the local environment, and is defined as the number of C-beta atoms in other residues within a sphere around the C-beta atom of the residue of interest. Contact number is partly conserved between protein folds and thus is useful for protein fold and structure prediction. In turn, each residue's contact number can be partially predicted from primary amino acid sequence, assisting tertiary fold analysis from sequence data. In this study, we provide a more accurate contact number prediction method from protein primary sequence. Results: We predict contact number from protein sequence using a novel support vector regression algorithm. Using protein local sequences with multiple sequence alignments (PSI-BLAST profiles), we demonstrate a correlation coefficient between predicted and observed contact numbers of 0.70, which outperforms previously achieved accuracies. Including additional information about sequence weight and amino acid composition further improves prediction accuracies significantly with the correlation coefficient reaching 0.73. If residues are classified as being either contacted or non-contacted, the prediction accuracies are all greater than 77%, regardless of the choice of classification thresholds. Conclusion: The successful application of support vector regression to the prediction of protein contact number reported here, together with previous applications of this approach to the prediction of protein accessible surface area and B-factor profile, suggests that a support vector regression approach may be very useful for determining the structure-function relation between primary sequence and higher order consecutive protein structural and functional properties.

Expressed sequence tag analysis of genes expressed during development of the tropical abalone Haliotis asinina

The tropical abalone. Haliotis asinina. is,in ideal species to investigate the molecular mechanisms that control development. growth, reproduction and shell formation in all cultured haliotids. Here we describe the analysis of 232 expressed sequence tags (EST) obtained front a developmental H. asinina cDNA library intended for future microarray studies. From this data set we identified 183 unique gene Clusters. Of these, 90 clusters showed significant homology with sequences lodged in GenBank, ranging in function from general housekeeping to signal transduction, gene regulation and cell-cell communication. Seventy-one clusters possessed completely novel ORFs greater than 50 codons in length, highlighting the paucity of sequence data from molluscs and other lophotrochozoans. This study of developmental gene expression in H. asinina provides the foundation for further detailed analyses of abalone growth, development and reproduction.

Genorne sequence of the thermostable Newcastle disease virus (strain I-2) reveals a possible phenotypic locus

The complete genome sequence of the Australian 1-2 heat-tolerant Newcastle disease virus (NDV) vaccine (master seed stocks) was determined and compared to the sequence of the parent virus from which it had been derived after exposure of the parent stock at 56 degrees C for 30 min. Nucleotide changes were observed at a number of positions with synonymous mutations being greater than those observed for non-synonymous mutations. Sequence data for the HN gene of a parental culture of V4 and two heat-tolerant variants of V4 were obtained. These were compared with the data for the 1-2 viruses and with published sequences for parental V4 and for a number of ND vaccine strains. Sequence analyses did not reveal the ARG 303 deletion in the HN protein, previously claimed to be responsible for the thermostable phenotype. No consistent changes were detected that would indicate involvement of the HN protein in heat resistance. The majority of alterations were observed in the L protein of the virus and it is proposed that these alterations were responsible for the heat-tolerant phenotype of the 1-2 NDV vaccine. (c) 2005 Elsevier B.V. All rights reserved.

Mechanisms of Thermal Adaptation Revealed From the Genomes of the Antarctic Archaea Methanogenium frigidum and Methanococcoides burtonii

We generated draft genome sequences for two cold-adapted Archaea, Methanogenium frigidum and Methanococcoides burtonii, to identify genotypic characteristics that distinguish them from Archaea with a higher optimal growth temperature (OGT). Comparative genomics revealed trends in amino acid and tRNA composition, and structural features of proteins. Proteins from the cold-adapted Archaea are characterized by a higher content of noncharged polar amino acids, particularly Gin and Thr and a lower content of hydrophobic amino acids, particularly Leu. Sequence data from nine methanogen genomes (OGT 15degrees-98degreesC) were used to generate IIII modeled protein structures. Analysis of the models from the cold-adapted Archaea showed a strong tendency in the solvent-accessible area for more Gin, Thr, and hydrophobic residues and fewer charged residues. A cold shock domain (CSD) protein (CspA homolog) was identified in M. frigidum, two hypothetical proteins with CSD-folds in M. burtonii, and a unique winged helix DNA-binding domain protein in M. burtonii. This suggests that these types of nucleic acid binding proteins have a critical role in cold-adapted Archaea. Structural analysis of tRNA sequences from the Archaea indicated that GC content is the major factor influencing tRNA stability in hyperthermophiles, but not in the psychrophiles, mesophiles or moderate thermophiles. Below an OGT of 60degreesC, the GC content in tRNA was largely unchanged, indicating that any requirement for flexibility of tRNA in psychrophiles is mediated by other means. This is the first time that comparisons have been performed with genome data from Archaea spanning the growth temperature extremes. from psychrophiles to hyperthermophiles

Sugarcane moth borers (Lepidoptera: Noctuidae and Pyraloidea): phylogenetics constructed using COII and 16S mitochondrial partial gene sequences

Sugarcane moth borers are a diverse group of species occurring in several genera, but predominately within the Noctuidae and Pyraloidea. They cause economic loss in sugarcane and other crops through damage to stems and stalks by larval boring. Partial sequence data from two mitochondrial genes, COII and 16S, were used to construct a molecular phylogeny based on 26 species from ten genera and six tribes. The Noctuidae were found to be monophyletic, providing molecular support for the taxonomy within this subfamily. However, the Pyraloidea are paraphyletic, with the noctuids splitting Galleriinae and Schoenobiinae from the Crambinae. This supports the separation of the Pyralidae and Crambinae, but does not support the concept of the incorporation of the Schoenobiinae in the Crambidae. Of the three crambine genera examined, Diatraea was monophyletic, Chilo paraphyletic, and Eoreuma was basal to the other two genera. Within the Noctuidae, Sesamia and Bathytricha were monophyletic, with Busseola basal to Bathytricha. Many species in this study (both noctuids and pyraloids) had different biotypes within collection localities and across their distribution; however the individual biotypes were not phylogenetically informative. These data highlight the need for taxonomic revisions at all taxon levels and provide a basis for the development of DNA-based diagnostics for rapidly identifying many species at any developmental stage. This ability is vital, as the species are an incursion threat to Australia and have the potential to cause significant losses to the sugar industry.

First record of Passalora calotropidis in Australia and its generic position

Passalora calotropidis has been found for the first time in Australia on rubber bush (Calotropis procera) in northern Queensland where it was associated with a damaging leaf spot disease. Analysis of sequence data of the ITS region indicated that P. calotropidis belonged to a group that consisted of species of Pseudocercospora. The generic position of P. calotropidis and its potential for biological control are discussed.

Evolution of duplicate control regions in the mitochondrial genomes of metazoa: A case study with Australasian Ixodes ticks

To investigate the evolution pattern and phylogenetic utility of duplicate control regions (CRs) in mitochondrial (mt) genomes, we sequenced the entire mt genomes of three Ixodes species and part of the mt genomes of another I I species. All the species from the Australasian lineage have duplicate CRs, whereas the other species have one CR. Sequence analyses indicate that the two CRs of the Australasian Ixodes ticks have evolved in concert in each species. In addition to the Australasian Ixodes ticks, species from seven other lineages of metazoa also have mt genomes with duplicate CRs. Accumulated mtDNA sequence data from these metazoans and two recent experiments on replication of mt genomes in human cell lines with duplicate CRs allowed us to re-examine four intriguing questions about the presence of duplicate CRs in the mt genomes of metazoa: (1) Why do some mt genomes, but not others, have duplicate CRs? (2) How did mt genomes with duplicate CRs evolve? (3) How could the nucleotide sequences of duplicate CRs remain identical or very similar over evolutionary time? (4) Are duplicate CRs phylogenetic markers? It appears that mt genomes with duplicate CRs have a selective advantage in replication over mt genomes with one CR. Tandem duplication followed by deletion of genes is the most plausible mechanism for the generation of mt genomes with duplicate CRs. Once duplicate CRs occur in an mt genome, they tend to evolve in concert, probably by gene conversion. However, there are lineages where gene conversion may not always occur, and, thus, the two CRs may evolve independently in these lineages. Duplicate CRs have much potential as phylogenetic markers at low taxonomic levels, such as within genera, within families, or among families, but not at high taxonomic levels, such as among orders.

The use and implications of ribosomal DNA sequencing for the discrimination of digenean species

In just over a decade, the use of molecular approaches for the recognition of parasites has become commonplace. For trematodes, the internal transcribed spacer region of ribosomal DNA (ITS rDNA) has become the default region of choice. Here, we review the findings of 63 studies that report ITS rDNA sequence data for about 155 digenean species from 19 families, and then review the levels of variation that have been reported and how the variation has been interpreted. Overall, complete ITS sequences (or ITS1 or ITS2 regions alone) usually distinguish trematode species clearly, including combinations for which morphology gives ambiguous results. Closely related species may have few base differences and in at least one convincing case the ITS2 sequences of two good species are identical. In some cases, the ITS1 region gives greater resolution than the ITS2 because of the presence of variable repeat units that are generally lacking in the ITS2. Intraspecific variation is usually low and frequently apparently absent. Information on geographical variation of digeneans is limited but at least some of the reported variation probably reflects the presence of multiple species. Despite the accepted dogma that concerted evolution makes the individual representative of the entire species, a significant number of studies have reported at least some intraspecific variation. The significance of such variation is difficult to assess a posteriori, but it seems likely that identification and sequencing errors account for some of it and failure to recognise separate species may also be significant. Some reported variation clearly requires further analysis. The use of a yardstick to determine when separate species should be recognised is flawed. Instead, we argue that consistent genetic differences that are associated with consistent morphological or biological traits should be considered the marker for separate species. We propose a generalised approach to the use of rDNA to distinguish trematode species.

Discovery of a new source of rifamycin antibiotics in marine sponge actinobacteria by phylogenetic prediction

Phylogenetic analysis of the ketosynthase (KS) gene sequences of marine sponge-derived Salinispora strains of actinobacteria indicated that the polyketide synthase (PKS) gene sequence most closely related to that of Salinispora was the rifamycin B synthase of Amycolatopsis mediterranei. This result was not expected from taxonomic species tree phylogenetics using 16S rRNA sequences. From the PKS sequence data generated from our sponge-derived Salinispora strains, we predicted that such strains might synthesize rifamycin-like compounds. Liquid chromatography-tandem mass spectrometry (LC/MS/MS) analysis was applied to one sponge-derived Salinispora strain to test the hypothesis of rifamycin synthesis. The analysis reported here demonstrates that this Salinispora isolate does produce compounds of the rifamycin class, including rifamycin B and rifamycin SV. A rifamycin-specific KS primer set was designed, and that primer set increased the number of rifamycin-positive strains detected by PCR screening relative to the number detectable using a conserved KS-specific set. Thus, the Salinispora group of actinobacteria represents a potential new source of rifamycins outside the genus Amycolatopsis and the first recorded source of rifamycins from marine bacteria.

Characterization of Kudoa monodactyli n. sp (Myxosporea : Multivalvulida) from the muscle of Monodactylus argenteus (Teleostei : Monodactylidae) from Moreton Bay, Queensland, Australia

Kudoa monodactyli n. sp. is described from the somatic musculature of Monodactylus argenteus from several localities in southern Queensland, Australia. This is the first record of a myxozoan parasite from the family Monodactylidae. The spores typically have five polar capsules, making this species similar to the four other five-valved Kudoa species (K. neurophila, K. muscularis, K. shulmani, K. cutanea) that have been described to date. However, morphometric measurements particularly of spore length and width make the species from M. argenteus distinct from the other species. Comparison of the small subunit ribosomal DNA sequence of this species with its congeners for which sequence data are available, provides further evidence of novelty. Kudoa monodactyli n. sp. displays 38 (of 1,554) nucleotide differences compared with rDNA sequence of Kudoa neurophila, which on phylogenetic analysis places these species in clades exclusive of each other. Phylogenetic analyses also provide evidence that the number of valves per spore in this genus is an imperfect indicator of relatedness.

Evolutionary dynamics of wAu-Like Wolbachia variants in neotropical Drosophila spp.

Wolbachia bacteria are common intracellular symbionts of arthropods and have been extensively studied in Drosophila. Most research focuses on two Old Word hosts, Drosophila melanogaster and Drosophila simulans, and does not take into account that some of the Wolbachia associations in these species may have evolved only after their fast global expansion and after the exposure to Wolbachia of previously isolated habitats. Here we looked at Wolbachia of Neotropical Drosophila species. Seventy-one lines of 16 Neotropical Drosophild species sampled in different regions and at different time points were analyzed. Wolbachia is absent in lines of Drosophild willistoni collected before the 1970s, but more recent samples are infected with a strain designated wWiL Wolbachia is absent in all other species of the willistoni group. Polymorphic wWil-related strains were detected in some saltans group species, with D. septentriosaltans being coinfected with at least four variants. Based on wsp and ftsZ sequence data, wWil of D. willistoni is identical to wAu, a strain isolated from D. simulans, but can be discriminated when using a polymorphic minisatellite marker. In contrast to wAu, which infects both germ line and somatic tissues of D. simulans, wWil is found exclusively in the primordial germ line cells of D. willistoni embryos. We report on a pool of closely related Wolbachia strains in Neotropical Drosophila species as a potential source for the wAu strain in D. simulans. Possible evolutionary scenarios reconstructing the infection history of wAu-like Wolbachia in Neotropical Drosophild species and the Old World species D. simulans are discussed.

Molecular characterization of Thelastomatoidea (Nematoda : Oxyurida) from cockroaches in Australia

A molecular approach was used to genetically characterize 5 species (Aoruroides queenslandensis. Blattophila sphaerolaima, Cordonicola gibsoni, Desmicola ornato and Leidynemella fusiformis) belonging to the superfamily. Thelastomatoidea fi (Nematoda: Oxyurida), a group of pinworms that parasitizes terrestrial arthropods. The D3 domain of the large subunit Of nuclear ribosomal RNA (LSU) was sequenced for individual specimens, and the analysis of the sequence data allowed the genetic relationships of the 5 species to be studied dagger. The sequence variation in the D3 domain within individual species (0-1-8%) was significantly less than the differences among species (4(.)3-12(.)4%). Phylogenetic analyses, Using maximum parsimony, maximum likelihood, and neighbour-joining, tree-building methods, established relationships among the 5 species of Thelastomatoidea and Oxyuris equi (a species of the order Oxyurida). The molecular approach employed provides the prospect for developing DNA tools for the specific identification of the Thelastomatoidea, irrespective of developmental stage and sex, as a basis for systematic, ecological and/or population genetic investigations of members within this superfamily.

Molecular systematics of the Philippine malaria vector Anopheles flavirostris

Allozyme and molecular sequence data from the malaria vector Anopheles flavirostris (Ludlow) (Diptera: Culicidae) were analysed from 34 sites throughout the Philippines, including the type locality, to test the hypothesis that this taxon is a single panmictic species. A finer-scaled allozyme study, of mainly Luzon samples, revealed no fixed genetic differences in sympatric sites and only low levels of variation. We obtained data from partial sequences for the internal transcribed spacer 2 (ITS2) (483 bp), the third domain (D3) (330 bp) of the 28S ribosomal DNA subunit and cytochrome c oxidase subunit I (COI) of mitochondrial DNA (261 bp). No sequence variation was observed for ITS2, only a one base pair difference was observed between Philippine and Indonesian D3 sequences and An. flavirostris sequences were unique, confirming their diagnostic value for this taxon. Sixteen COI haplotypes were identified, giving 25 parsimony informative sites. Neighbour-Joining, Maximum Parsimony, Maximum Likelihood and Bayesian phylogenetic analysis of COI sequences for An. flavirostris and outgroup taxa revealed strong branch support for the monophyly of An. flavirostris, thus confirming that Philippine populations of this taxon comprise a single separate species within the Minimus Subgroup of the Funestus Group. Variation in the behaviour of An. flavirostris is likely to be intraspecific rather than interspecific in origin. © 2006 The Royal Entomological Society.

Sequencing a new target genome: The Boophilus microplus (Acari : Ixodidae) genome project

The southern cattle tick, Boophilus microplus (Canestrini), causes annual economic losses in the hundreds of millions of dollars to cattle producers throughout the world, and ranks as the most economically important tick from a global perspective. Control failures attributable to the development of pesticide resistance have become commonplace, and novel control technologies are needed. The availability of the genome sequence will facilitate the development of these new technologies, and we are proposing sequencing to a 4-6X draft coverage. Many existing biological resources are available to facilitate a genome sequencing project, including several inbred laboratory tick strains, a database of approximate to 45,000 expressed sequence tags compiled into a B. microplus Gene Index, a bacterial artificial chromosome (BAC) library, an established B. microplus cell line, and genomic DNA suitable for library synthesis. Collaborative projects are underway to map BACs and cDNAs to specific chromosomes and to sequence selected BAC clones. When completed, the genome sequences from the cow, B. microphis, and the B. microphis-borne pathogens Babesia bovis and Anaplasma marginale will enhance studies of host-vector-pathogen systems. Genes involved in the regeneration of amputated tick limbs and transitions through developmental stages are largely unknown. Studies of these and other interesting biological questions will be advanced by tick genome sequence data. Comparative genomics offers the prospect of new insight into many, perhaps all, aspects of the biology of ticks and the pathogens they transmit to farm animals and people. The B. microplus genome sequence will fill a major gap in comparative genomics: a sequence from the Metastriata lineage of ticks. The purpose of the article is to synergize interest in and provide rationales for sequencing the genome of B. microplus and for publicizing currently available genomic resources for this tick.