Mechanical Engineering - Research Publications

Permanent URI for this collection

Search Results

Now showing 1 - 10 of 21
  • Item
    Thumbnail Image
    Discovering the pharmacodynamics of conolidine and cannabidiol using a cultured neuronal network based workflow
    Mendis, GDC ; Berecki, G ; Morrisroe, E ; Pachernegg, S ; Li, M ; Varney, M ; Osborne, PB ; Reid, CA ; Halgamuge, S ; Petrou, S (NATURE PUBLISHING GROUP, 2019-01-15)
    Determining the mechanism of action (MOA) of novel or naturally occurring compounds mostly relies on assays tailored for individual target proteins. Here we explore an alternative approach based on pattern matching response profiles obtained using cultured neuronal networks. Conolidine and cannabidiol are plant-derivatives with known antinociceptive activity but unknown MOA. Application of conolidine/cannabidiol to cultured neuronal networks altered network firing in a highly reproducible manner and created similar impact on network properties suggesting engagement with a common biological target. We used principal component analysis (PCA) and multi-dimensional scaling (MDS) to compare network activity profiles of conolidine/cannabidiol to a series of well-studied compounds with known MOA. Network activity profiles evoked by conolidine and cannabidiol closely matched that of ω-conotoxin CVIE, a potent and selective Cav2.2 calcium channel blocker with proposed antinociceptive action suggesting that they too would block this channel. To verify this, Cav2.2 channels were heterologously expressed, recorded with whole-cell patch clamp and conolidine/cannabidiol was applied. Remarkably, conolidine and cannabidiol both inhibited Cav2.2, providing a glimpse into the MOA that could underlie their antinociceptive action. These data highlight the utility of cultured neuronal network-based workflows to efficiently identify MOA of drugs in a highly scalable assay.
  • Item
    Thumbnail Image
    Unsupervised discovery of microbial population structure within metagenomes using nucleotide base composition
    Saeed, I ; Tang, S-L ; Halgamuge, SK (OXFORD UNIV PRESS, 2012-03)
    An approach to infer the unknown microbial population structure within a metagenome is to cluster nucleotide sequences based on common patterns in base composition, otherwise referred to as binning. When functional roles are assigned to the identified populations, a deeper understanding of microbial communities can be attained, more so than gene-centric approaches that explore overall functionality. In this study, we propose an unsupervised, model-based binning method with two clustering tiers, which uses a novel transformation of the oligonucleotide frequency-derived error gradient and GC content to generate coarse groups at the first tier of clustering; and tetranucleotide frequency to refine these groups at the secondary clustering tier. The proposed method has a demonstrated improvement over PhyloPythia, S-GSOM, TACOA and TaxSOM on all three benchmarks that were used for evaluation in this study. The proposed method is then applied to a pyrosequenced metagenomic library of mud volcano sediment sampled in southwestern Taiwan, with the inferred population structure validated against complementary sequencing of 16S ribosomal RNA marker genes. Finally, the proposed method was further validated against four publicly available metagenomes, including a highly complex Antarctic whale-fall bone sample, which was previously assumed to be too complex for binning prior to functional analysis.
  • Item
    Thumbnail Image
    CoNVEX: copy number variation estimation in exome sequencing data using HMM
    Amarasinghe, KC ; Li, J ; Halgamuge, SK (BioMed Central, 2013)
    BACKGROUND: One of the main types of genetic variations in cancer is Copy Number Variations (CNV). Whole exome sequencing (WES) is a popular alternative to whole genome sequencing (WGS) to study disease specific genomic variations. However, finding CNV in Cancer samples using WES data has not been fully explored. RESULTS: We present a new method, called CoNVEX, to estimate copy number variation in whole exome sequencing data. It uses ratio of tumour and matched normal average read depths at each exonic region, to predict the copy gain or loss. The useful signal produced by WES data will be hindered by the intrinsic noise present in the data itself. This limits its capacity to be used as a highly reliable CNV detection source. Here, we propose a method that consists of discrete wavelet transform (DWT) to reduce noise. The identification of copy number gains/losses of each targeted region is performed by a Hidden Markov Model (HMM). CONCLUSION: HMM is frequently used to identify CNV in data produced by various technologies including Array Comparative Genomic Hybridization (aCGH) and WGS. Here, we propose an HMM to detect CNV in cancer exome data. We used modified data from 1000 Genomes project to evaluate the performance of the proposed method. Using these data we have shown that CoNVEX outperforms the existing methods significantly in terms of precision. Overall, CoNVEX achieved a sensitivity of more than 92% and a precision of more than 50%.
  • Item
    Thumbnail Image
    Bioinformatics Pipelines for Targeted Resequencing and Whole-Exome Sequencing of Human and Mouse Genomes: A Virtual Appliance Approach for Instant Deployment
    Li, J ; Doyle, MA ; Saeed, I ; Wong, SQ ; Mar, V ; Goode, DL ; Caramia, F ; Doig, K ; Ryland, GL ; Thompson, ER ; Hunter, SM ; Halgamuge, SK ; Ellul, J ; Dobrovic, A ; Campbell, IG ; Papenfuss, AT ; McArthur, GA ; Tothill, RW ; Calogero, RA (PUBLIC LIBRARY SCIENCE, 2014-04-21)
    Targeted resequencing by massively parallel sequencing has become an effective and affordable way to survey small to large portions of the genome for genetic variation. Despite the rapid development in open source software for analysis of such data, the practical implementation of these tools through construction of sequencing analysis pipelines still remains a challenging and laborious activity, and a major hurdle for many small research and clinical laboratories. We developed TREVA (Targeted REsequencing Virtual Appliance), making pre-built pipelines immediately available as a virtual appliance. Based on virtual machine technologies, TREVA is a solution for rapid and efficient deployment of complex bioinformatics pipelines to laboratories of all sizes, enabling reproducible results. The analyses that are supported in TREVA include: somatic and germline single-nucleotide and insertion/deletion variant calling, copy number analysis, and cohort-based analyses such as pathway and significantly mutated genes analyses. TREVA is flexible and easy to use, and can be customised by Linux-based extensions if required. TREVA can also be deployed on the cloud (cloud computing), enabling instant access without investment overheads for additional hardware. TREVA is available at http://bioinformatics.petermac.org/treva/.
  • Item
    Thumbnail Image
    Prokaryotic assemblages and metagenomes in pelagic zones of the South China Sea
    Tseng, C-H ; Chiang, P-W ; Lai, H-C ; Shiah, F-K ; Hsu, T-C ; Chen, Y-L ; Wen, L-S ; Tseng, C-M ; Shieh, W-Y ; Saeed, I ; Halgamuge, S ; Tang, S-L (BIOMED CENTRAL LTD, 2015-03-20)
    BACKGROUND: Prokaryotic microbes, the most abundant organisms in the ocean, are remarkably diverse. Despite numerous studies of marine prokaryotes, the zonation of their communities in pelagic zones has been poorly delineated. By exploiting the persistent stratification of the South China Sea (SCS), we performed a 2-year, large spatial scale (10, 100, 1000, and 3000 m) survey, which included a pilot study in 2006 and comprehensive sampling in 2007, to investigate the biological zonation of bacteria and archaea using 16S rRNA tag and shotgun metagenome sequencing. RESULTS: Alphaproteobacteria dominated the bacterial community in the surface SCS, where the abundance of Betaproteobacteria was seemingly associated with climatic activity. Gammaproteobacteria thrived in the deep SCS, where a noticeable amount of Cyanobacteria were also detected. Marine Groups II and III Euryarchaeota were predominant in the archaeal communities in the surface and deep SCS, respectively. Bacterial diversity was higher than archaeal diversity at all sampling depths in the SCS, and peaked at mid-depths, agreeing with the diversity pattern found in global water columns. Metagenomic analysis not only showed differential %GC values and genome sizes between the surface and deep SCS, but also demonstrated depth-dependent metabolic potentials, such as cobalamin biosynthesis at 10 m, osmoregulation at 100 m, signal transduction at 1000 m, and plasmid and phage replication at 3000 m. When compared with other oceans, urease at 10 m and both exonuclease and permease at 3000 m were more abundant in the SCS. Finally, enriched genes associated with nutrient assimilation in the sea surface and transposase in the deep-sea metagenomes exemplified the functional zonation in global oceans. CONCLUSIONS: Prokaryotic communities in the SCS stratified with depth, with maximal bacterial diversity at mid-depth, in accordance with global water columns. The SCS had functional zonation among depths and endemically enriched metabolic potentials at the study site, in contrast to other oceans.
  • Item
    Thumbnail Image
    Accurate reconstruction of viral quasispecies spectra through improved estimation of strain richness
    Jayasundara, D ; Saeed, I ; Chang, BC ; Tang, S-L ; Halgamuge, SK (BMC, 2015-12-09)
    BACKGROUND: Estimating the number of different species (richness) in a mixed microbial population has been a main focus in metagenomic research. Existing methods of species richness estimation ride on the assumption that the reads in each assembled contig correspond to only one of the microbial genomes in the population. This assumption and the underlying probabilistic formulations of existing methods are not useful for quasispecies populations where the strains are highly genetically related. RESULTS: On benchmark data sets, our estimation method provided accurate richness estimates (< 0.2 median estimation error) and improved the precision of ViQuaS by 2%-13% and F-score by 1%-9% without compromising the recall rates. We also demonstrate that our estimation method can be used to improve the precision and F-score of ShoRAH by 0%-7% and 0%-5% respectively. CONCLUSIONS: The proposed probabilistic estimation method can be used to estimate the richness of viral populations with a quasispecies behavior and to improve the accuracy of the quasispecies spectra reconstructed by the existing methods ViQuaS and ShoRAH in the presence of a moderate level of technical sequencing errors. AVAILABILITY: http://sourceforge.net/projects/viquas/.
  • Item
    Thumbnail Image
    Assessing Species Diversity Using Metavirome Data: Methods and Challenges
    Herath, D ; Jayasundara, D ; Ackland, D ; Saeed, I ; Tang, S-L ; Halgamuge, S (ELSEVIER SCIENCE BV, 2017)
    Assessing biodiversity is an important step in the study of microbial ecology associated with a given environment. Multiple indices have been used to quantify species diversity, which is a key biodiversity measure. Measuring species diversity of viruses in different environments remains a challenge relative to measuring the diversity of other microbial communities. Metagenomics has played an important role in elucidating viral diversity by conducting metavirome studies; however, metavirome data are of high complexity requiring robust data preprocessing and analysis methods. In this review, existing bioinformatics methods for measuring species diversity using metavirome data are categorised broadly as either sequence similarity-dependent methods or sequence similarity-independent methods. The former includes a comparison of DNA fragments or assemblies generated in the experiment against reference databases for quantifying species diversity, whereas estimates from the latter are independent of the knowledge of existing sequence data. Current methods and tools are discussed in detail, including their applications and limitations. Drawbacks of the state-of-the-art method are demonstrated through results from a simulation. In addition, alternative approaches are proposed to overcome the challenges in estimating species diversity measures using metavirome data.
  • Item
    Thumbnail Image
    The node-weighted Steiner tree approach to identify elements of cancer-related signaling pathways
    Sun, Y ; Ma, C ; Halgamuge, S (BMC, 2017-12-28)
    BACKGROUND: Cancer constitutes a momentous health burden in our society. Critical information on cancer may be hidden in its signaling pathways. However, even though a large amount of money has been spent on cancer research, some critical information on cancer-related signaling pathways still remains elusive. Hence, new works towards a complete understanding of cancer-related signaling pathways will greatly benefit the prevention, diagnosis, and treatment of cancer. RESULTS: We propose the node-weighted Steiner tree approach to identify important elements of cancer-related signaling pathways at the level of proteins. This new approach has advantages over previous approaches since it is fast in processing large protein-protein interaction networks. We apply this new approach to identify important elements of two well-known cancer-related signaling pathways: PI3K/Akt and MAPK. First, we generate a node-weighted protein-protein interaction network using protein and signaling pathway data. Second, we modify and use two preprocessing techniques and a state-of-the-art Steiner tree algorithm to identify a subnetwork in the generated network. Third, we propose two new metrics to select important elements from this subnetwork. On a commonly used personal computer, this new approach takes less than 2 s to identify the important elements of PI3K/Akt and MAPK signaling pathways in a large node-weighted protein-protein interaction network with 16,843 vertices and 1,736,922 edges. We further analyze and demonstrate the significance of these identified elements to cancer signal transduction by exploring previously reported experimental evidences. CONCLUSIONS: Our node-weighted Steiner tree approach is shown to be both fast and effective to identify important elements of cancer-related signaling pathways. Furthermore, it may provide new perspectives into the identification of signaling pathways for other human diseases.
  • Item
    Thumbnail Image
    Estimating neuronal conductance model parameters using dynamic action potential clamp
    Deerasooriya, Y ; Berecki, G ; Kaplan, D ; Forster, IC ; Halgamuge, S ; Petrou, S (Elsevier, 2019-09-01)
    Parameterization of neuronal membrane conductance models relies on data acquired from current clamp (CC) or voltage clamp (VC) recordings. Although the CC approach provides key information on a neuron’s firing properties, it is often difficult to disentangle the influence of multiple conductances that contribute to the excitation properties of a real neuron. Isolation of a single conductance using pharmacological agents or heterologous expression simplifies analysis but requires extensive VC evaluation to explore the complete state behavior of the channel of interest.
  • Item
    Thumbnail Image
    Inferring copy number and genotype in tumour exome data
    Amarasinghe, KC ; Li, J ; Hunter, SM ; Ryland, GL ; Cowin, PA ; Campbell, IG ; Halgamuge, SK (BMC, 2014-08-28)
    BACKGROUND: Using whole exome sequencing to predict aberrations in tumours is a cost effective alternative to whole genome sequencing, however is predominantly used for variant detection and infrequently utilised for detection of somatic copy number variation. RESULTS: We propose a new method to infer copy number and genotypes using whole exome data from paired tumour/normal samples. Our algorithm uses two Hidden Markov Models to predict copy number and genotypes and computationally resolves polyploidy/aneuploidy, normal cell contamination and signal baseline shift. Our method makes explicit detection on chromosome arm level events, which are commonly found in tumour samples. The methods are combined into a package named ADTEx (Aberration Detection in Tumour Exome). We applied our algorithm to a cohort of 17 in-house generated and 18 TCGA paired ovarian cancer/normal exomes and evaluated the performance by comparing against the copy number variations and genotypes predicted using Affymetrix SNP 6.0 data of the same samples. Further, we carried out a comparison study to show that ADTEx outperformed its competitors in terms of precision and F-measure. CONCLUSIONS: Our proposed method, ADTEx, uses both depth of coverage ratios and B allele frequencies calculated from whole exome sequencing data, to predict copy number variations along with their genotypes. ADTEx is implemented as a user friendly software package using Python and R statistical language. Source code and sample data are freely available under GNU license (GPLv3) at http://adtex.sourceforge.net/.