Algorithms for sequence analysis via mutagenesis

Despite many successes of conventional DNA sequencing methods, some DNAs remain difficult or impossible to sequence. Unsequenceable regions occur in the genomes of many biologically important organisms, including the human genome. Such regions range in length from tens to millions of bases, and may contain valuable information such as the sequences of important genes. The authors have recently developed a technique that renders a wide range of problematic DNAs amenable to sequencing. The technique is known as sequence analysis via mutagenesis (SAM). This paper presents a number of algorithms for analysing and interpreting data generated by this technique.

Sequence analysis of rabbit haemorrhagic disease virus (RHDV) in Australia: alterations after its release

Liver samples from rabbits killed by RHDV, collected from five States in Australia in 1996 and 1997 were analysed by RT-PCR. A 398 bp fragment of the capsid protein (VP60) gene was amplified by PCR and directly sequenced. The alignment of the nucleotide and amino acid sequences and their comparison with the original strain of the virus released in Australia indicated genetic changes after two years have been small with 98.2% to 100% identity. The constructed phylogenetic tree suggests slight differences in nucleotide substitutions in various States but there is no clear evidence of clustering of sequences according to their geographic origin. In practical terms, sequencing of viral RNA provides a means of testing the efficacy of further releases and subsequent spread of the virus if such a strategy is employed as a means of enhancing RHD as a biological control of the wild rabbit in Australia.

Sequence analysis and expression of a virus-like particle protein,VLP2, from the parasitic wasp Venturia canescens

Endoparasitoid wasps produce maternal protein secretions, which are transported into the body of insect hosts at oviposition to regulate host physiology for successful development of their offspring. Venturia canescens calyx fluid contains so-called virus-like particles (VLPs) that are essential for immune evasion of the developing parasitoid inside the host. VLPs consist of four major proteins. In this paper, we describe the isolation and molecular cloning of a gene (vlp2) that is a constituent of VLPs and discuss its possible role in VLP structure and function.

Molecular characterization of the microbial species that colonize human ileal and colonic mucosa by using 16S rDNA sequence analysis

Aim: The aim of this study was to characterize the bacterial community adhering to the mucosa of the terminal ileum, and proximal and distal colon of the human digestive tract. Methods and Results: Pinch samples of the terminal ileum, proximal and distal colon were taken from a healthy 35-year-old, and a 68-year-old subject with mild diverticulosis. The 16S rDNA genes were amplified using a low number of PCR cycles, cloned, and sequenced. In total, 361 sequences were obtained comprising 70 operational taxonomic units (OTU), with a calculated coverage of 82.6%. Twenty-three per cent of OTU were common to the terminal ileum, proximal colon and distal colon, but 14% OTU were only found in the terminal ileum, and 43% were only associated with the proximal or distal colon. The most frequently represented clones were from the Clostridium group XIVa (24.7%), and the Bacteroidetes (Cytophaga-Flavobacteria-Bacteroides ) cluster (27.7%). Conclusion: Comparison of 16S rDNA clone libraries of the hindgut across mammalian species confirms that the distribution of phylogenetic groups is similar irrespective of the host species. Lesser site-related differences within groups or clusters of organisms, are probable. Significance and Impact: This study provides further evidence of the distribution of the bacteria on the mucosal surfaces of the human hindgut. Data contribute to the benchmarking of the microbial composition of the human digestive tract.

The applicability of recurrent neural networks for biological sequence analysis

Selection of machine learning techniques requires a certain sensitivity to the requirements of the problem. In particular, the problem can be made more tractable by deliberately using algorithms that are biased toward solutions of the requisite kind. In this paper, we argue that recurrent neural networks have a natural bias toward a problem domain of which biological sequence analysis tasks are a subset. We use experiments with synthetic data to illustrate this bias. We then demonstrate that this bias can be exploitable using a data set of protein sequences containing several classes of subcellular localization targeting peptides. The results show that, compared with feed forward, recurrent neural networks will generally perform better on sequence analysis tasks. Furthermore, as the patterns within the sequence become more ambiguous, the choice of specific recurrent architecture becomes more critical.

Cloning and sequence analysis of pituitary prolactin cDNA from the northern brown bandicoot (Isoodon macrourus)

The nuclectide sequence for pituitary prolactin cDNA from the marsupial bandicoot (Isoodon macrourus) was determined by reverse transcription-polymerase chain reaction and 5'/3' rapid amplification of cDNA ends. The deduced amino acid sequence showed high sequence identity with brushtail possum prolactin (95%) and all of the expected structural features of a quadruped prolactin. A prolactin gene tree was constructed and rates of evolution calculated for bandicoot, possum, opossum and several mammalian and non-mammalian prolactins. Bootstrap analysis provided strong support for marsupials as a sister group with eutherian mammals and weak support for opossum and bandicoot as an independent grouping from the brushtail possum. The rates of molecular evolution for marsupial prolactins were comparable to the slow rate seen in the majority of quadruped prolactins that have been sequenced. (c) 2005 Elsevier Inc. All rights reserved.

Seizure detection algorithm for neonates based on wave-sequence analysis

Objective: The description and evaluation of the performance of a new real-time seizure detection algorithm in the newborn infant. Methods: The algorithm includes parallel fragmentation of EEG signal into waves; wave-feature extraction and averaging; elementary, preliminary and final detection. The algorithm detects EEG waves with heightened regularity, using wave intervals, amplitudes and shapes. The performance of the algorithm was assessed with the use of event-based and liberal and conservative time-based approaches and compared with the performance of Gotman's and Liu's algorithms. Results: The algorithm was assessed on multi-channel EEG records of 55 neonates including 17 with seizures. The algorithm showed sensitivities ranging 83-95% with positive predictive values (PPV) 48-77%. There were 2.0 false positive detections per hour. In comparison, Gotman's algorithm (with 30 s gap-closing procedure) displayed sensitivities of 45-88% and PPV 29-56%; with 7.4 false positives per hour and Liu's algorithm displayed sensitivities of 96-99%, and PPV 10-25%; with 15.7 false positives per hour. Conclusions: The wave-sequence analysis based algorithm displayed higher sensitivity, higher PPV and a substantially lower level of false positives than two previously published algorithms. Significance: The proposed algorithm provides a basis for major improvements in neonatal seizure detection and monitoring. Published by Elsevier Ireland Ltd. on behalf of International Federation of Clinical Neurophysiology.

Analysis of the TGF beta functional pathway in epithelial ovarian carcinoma

Epithelial ovarian carcinoma is often diagnosed at an advanced stage of disease and is the leading cause of death from gynaecological neoplasia. The genetic changes that occur during the development of this carcinoma are poorly understood. It has been proposed that IGFIIR, TGF beta1 and TGF beta RII act as a functional unit in the TGF beta growth inhibitory pathway, and that somatic loss-of-function mutations in any one of these genes could lead to disruption of the pathway and subsequent loss of cell cycle control. We have examined these 3 genes in 25 epithelial ovarian carcinomas using single-stranded conformational polymorphism analysis and DNA sequence analysis. A total of 3 somatic missense mutations were found in the TGF beta RII gene, but none in IGFRII or TGF beta1. An association was found between TGF beta RII mutations and histology, with 2 out of 3 clear cell carcinomas having TGF beta RII mutations. This data supports other evidence from mutational analysis of the PTEN and beta -catenin genes that there are distinct developmental pathways responsible for the progression of different epithelial ovarian cancer histologic subtypes. (C) 2001 Cancer Research Campaign.

Molecular analysis of dimethyl sulphide dehydrogenase from Rhodovulum sulfidophilum: its place in the dimethyl sulphoxide reductase family of microbial molybdopterin-containing enzymes

Dimethyl sulphide dehydrogenase catalyses the oxidation of dimethyl sulphide to dimethyl sulphoxide (DMSO) during photoautotrophic growth of Rhodovulum sulfidophilum . Dimethyl sulphide dehydrogenase was shown to contain bis (molybdopterin guanine dinucleotide)Mo, the form of the pterin molybdenum cofactor unique to enzymes of the DMSO reductase family. Sequence analysis of the ddh gene cluster showed that the ddhA gene encodes a polypeptide with highest sequence similarity to the molybdop-terin-containing subunits of selenate reductase, ethylbenzene dehydrogenase. These polypeptides form a distinct clade within the DMSO reductase family. Further sequence analysis of the ddh gene cluster identified three genes, ddhB , ddhD and ddhC . DdhB showed sequence homology to NarH, suggesting that it contains multiple iron-sulphur clusters. Analysis of the N-terminal signal sequence of DdhA suggests that it is secreted via the Tat secretory system in complex with DdhB, whereas DdhC is probably secreted via a Sec-dependent mechanism. Analysis of a ddhA mutant showed that dimethyl sulphide dehydrogenase was essential for photolithotrophic growth of Rv. sulfidophilum on dimethyl sulphide but not for chemo-trophic growth on the same substrate. Mutational analysis showed that cytochrome c (2) mediated photosynthetic electron transfer from dimethyl sulphide dehydrogenase to the photochemical reaction centre, although this cytochrome was not essential for photoheterotrophic growth of the bacterium.