A transcriptomic analysis of striped catfish (Pangasianodon hypophthalmus) in response to salinity adaptation: De novo assembly, gene annotation and marker discovery

The striped catfish (Pangasianodon hypophthalmus) culture industry in the Mekong Delta in Vietnam has developed rapidly over the past decade. The culture industry now however, faces some significant challenges, especially related to climate change impacts notably from predicted extensive saltwater intrusion into many low topographical coastal provinces across the Mekong Delta. This problem highlights a need for development of culture stocks that can tolerate more saline culture environments as a response to expansion of saline water-intruded land. While a traditional artificial selection program can potentially address this need, understanding the genomic basis of salinity tolerance can assist development of more productive culture lines. The current study applied a transcriptomic approach using Ion PGM technology to generate expressed sequence tag (EST) resources from the intestine and swim bladder from striped catfish reared at a salinity level of 9 ppt which showed best growth performance. Total sequence data generated was 467.8 Mbp, consisting of 4,116,424 reads with an average length of 112 bp. De novo assembly was employed that generated 51,188 contigs, and allowed identification of 16,116 putative genes based on the GenBank non-redundant database. GO annotation, KEGG pathway mapping, and functional annotation of the EST sequences recovered with a wide diversity of biological functions and processes. In addition, more than 11,600 simple sequence repeats were also detected. This is the first comprehensive analysis of a striped catfish transcriptome, and provides a valuable genomic resource for future selective breeding programs and functional or evolutionary studies of genes that influence salinity tolerance in this important culture species.

How to sequence and annotate insect mitochondrial genomes for systematic and comparative genomics research

Over the past decade the mitochondrial (mt) genome has become the most widely used genomic resource available for systematic entomology. While the availability of other types of ‘–omics’ data – in particular transcriptomes – is increasing rapidly, mt genomes are still vastly cheaper to sequence and are far less demanding of high quality templates. Furthermore, almost all other ‘–omics’ approaches also sequence the mt genome, and so it can form a bridge between legacy and contemporary datasets. Mitochondrial genomes have now been sequenced for all insect orders, and in many instances representatives of each major lineage within orders (suborders, series or superfamilies depending on the group). They have also been applied to systematic questions at all taxonomic scales from resolving interordinal relationships (e.g. Cameron et al., 2009; Wan et al., 2012; Wang et al., 2012), through many intraordinal (e.g. Dowton et al., 2009; Timmermans et al., 2010; Zhao et al. 2013a) and family-level studies (e.g. Nelson et al., 2012; Zhao et al., 2013b) to population/biogeographic studies (e.g. Ma et al., 2012). Methodological issues around the use of mt genomes in insect phylogenetic analyses and the empirical results found to date have recently been reviewed by Cameron (2014), yet the technical aspects of sequencing and annotating mt genomes were not covered. Most papers which generate new mt genome report their methods in a simplified form which can be difficult to replicate without specific knowledge of the field. Published studies utilize a sufficiently wide range of approaches, usually without justification for the one chosen, that confusion about commonly used jargon such as ‘long PCR’ and ‘primer walking’ could be a serious barrier to entry. Furthermore, sequenced mt genomes have been annotated (gene locations defined) to wildly varying standards and improving data quality through consistent annotation procedures will benefit all downstream users of these datasets. The aims of this review are therefore to: 1. Describe in detail the various sequencing methods used on insect mt genomes; 2. Explore the strengths/weakness of different approaches; 3. Outline the procedures and software used for insect mt genome annotation, and; 4. Highlight quality control steps used for new annotations, and to improve the re-annotation of previously sequenced mt genomes used in systematic or comparative research.

Human single-stranded DNA binding protein 1 (hSSB1/NABP2) is required for the stability and repair of stalled replication forks

Aberrant DNA replication is a primary cause of mutations that are associated with pathological disorders including cancer. During DNA metabolism, the primary causes of replication fork stalling include secondary DNA structures, highly transcribed regions and damaged DNA. The restart of stalled replication forks is critical for the timely progression of the cell cycle and ultimately for the maintenance of genomic stability. Our previous work has implicated the single-stranded DNA binding protein, hSSB1/NABP2, in the repair of DNA double-strand breaks via homologous recombination. Here, we demonstrate that hSSB1 relocates to hydroxyurea (HU)-damaged replication forks where it is required for ATR and Chk1 activation and recruitment of Mre11 and Rad51. Consequently, hSSB1-depleted cells fail to repair and restart stalled replication forks. hSSB1 deficiency causes accumulation of DNA strand breaks and results in chromosome aberrations observed in mitosis, ultimately resulting in hSSB1 being required for survival to HU and camptothecin. Overall, our findings demonstrate the importance of hSSB1 in maintaining and repairing DNA replication forks and for overall genomic stability.

Chemotherapeutic compounds targeting the DNA double-strand break repair pathways : the good, the bad, and the promising

The repair of DNA double-strand breaks (DSBs) is a critical cellular mechanism that exists to ensure genomic stability. DNA DSBs are the most deleterious type of insult to a cell’s genetic material and can lead to genomic instability, apoptosis, or senescence. Incorrectly repaired DNA DSBs have the potential to produce chromosomal translocations and genomic instability, potentially leading to cancer. The prevalence of DNA DSBs in cancer due to unregulated growth and errors in repair opens up a potential therapeutic window in the treatment of cancers. The cellular response to DNA DSBs is comprised of two pathways to ensure DNA breaks are repaired: homologous recombination and non-homologous end joining. Identifying chemotherapeutic compounds targeting proteins involved in these DNA repair pathways has shown promise as a cancer therapy for patients, either as a monotherapy or in combination with genotoxic drugs. From the beginning, there have been a number of chemotherapeutic compounds that have yielded successful responses in the clinic, a number that have failed (CGK-733 and iniparib), and a number of promising targets for future studies identified. This review looks in detail at how the cell responds to these DNA DSBs and investigates the chemotherapeutic avenues that have been and are currently being explored to target this repair process.

A control theoretic paradigm for cell signaling networks : a simple complexity for a sensitive robustness

The fields of molecular biology and cell biology are being flooded with complex genomic and proteomic datasets of large dimensions. We now recognize that each molecule in the cell and tissue can no longer be viewed as an isolated entity. Instead, each molecule must be considered as one member of an interacting network. Consequently, there is an urgent need for mathematical models to understand the behavior of cell signaling networks in health and in disease.

Transforming the future of treatment for ovarian cancer

As for many other cancers, metastasis is the leading cause of death of patients with ovarian cancer. Vigorous basic and clinical research is being performed to initiate more efficacious treatment strategies to improve the poor outcome of women with this cancer. Current treatment for ovarian cancer includes advanced cyto-reductive surgery and traditional platinum and taxane combined chemotherapy. Clinical trials using novel cytotoxic reagents and tyrosine kinase inhibitors have also been progressing. In parallel, the application of robust unbiased high throughput research platforms using transcriptomic and proteomic approaches has identified that not only individual cell signalling pathways, but a network of molecular pathways, play an important role in the biology of ovarian cancer. Furthermore, intensive genomic and epigenetic analyses have also revealed single nucleotide polymorphisms associated with risk and/or aetiology of this cancer including patient response to treatment. Taken together, these approaches, that are advancing our understanding, will have an impact on the generation of new therapeutic approaches and strategies for improving the outcome and quality of life of patients with ovarian cancer in the near future.

Multimedia copyright protection

Internet and its widespread usage for multimedia document distribution put the copyright issue in a complete new setting. Multimedia documents, specifically those installed on a web page, are no longer passive as they typically include active applets. Copyright protection safeguards the intellectual property (IP) of multimedia documents, which are either sold or distributed free of charge. In this Chapter, the basic tools for copyright protection are discussed. First, general concepts and the vocabulary used in copyright protection of multimedia documents are discussed. Later, taxonomy of watermarking and fingerprinting techniques are studied. This part is concluded by a review of the literature dealing with IP security. The main part of the chapter discusses the generic watermarking scheme and illustrates it on three specific examples: collusion-free watermarking, spread spectrum watermarking, and software fingerprinting. Future trends and conclusions close the chapter.

Drosophila 3′ UTRs are more complex than protein-coding sequences

The 3′ UTRs of eukaryotic genes participate in a variety of post-transcriptional (and some transcriptional) regulatory interactions. Some of these interactions are well characterised, but an undetermined number remain to be discovered. While some regulatory sequences in 3′ UTRs may be conserved over long evolutionary time scales, others may have only ephemeral functional significance as regulatory profiles respond to changing selective pressures. Here we propose a sensitive segmentation methodology for investigating patterns of composition and conservation in 3′ UTRs based on comparison of closely related species. We describe encodings of pairwise and three-way alignments integrating information about conservation, GC content and transition/transversion ratios and apply the method to three closely related Drosophila species: D. melanogaster, D. simulans and D. yakuba. Incorporating multiple data types greatly increased the number of segment classes identified compared to similar methods based on conservation or GC content alone. We propose that the number of segments and number of types of segment identified by the method can be used as proxies for functional complexity. Our main finding is that the number of segments and segment classes identified in 3′ UTRs is greater than in the same length of protein-coding sequence, suggesting greater functional complexity in 3′ UTRs. There is thus a need for sustained and extensive efforts by bioinformaticians to delineate functional elements in this important genomic fraction. C code, data and results are available upon request.

Single nucleotide polymorphisms and mRNA expression of VEGF-A in papillary thyroid carcinoma: Potential markers for aggressive phenotypes

BACKGROUND AND OBJECTIVES Polymorphisms of the VEGF gene are known to affect the biological behaviour of cancers but have seldom been studied in thyroid cancer. The aim of the current study is to evaluate the prevalence and relevance of VEGF-A polymorphisms and mRNA expression in papillary thyroid carcinoma (PTC). MATERIALS AND METHODS Genomic DNA and total RNA were isolated from paraffin-embedded tissue from 91 PTC (51 conventional PTC and 40 follicular variant) and 78 control thyroid tissues. Three DNA polymorphisms (+936C > T, +405C > G and -141A > C) in the 3' and 5' untranslated region (3'-UTR, 5'-UTR) of VEGF-A were studied using PCR and RFLP. Also, the mRNA expression of VEGF-A in these tissues was studied by real-time PCR. RESULTS Distribution of polymorphisms in the 5'-UTR (VEGF-VEGF -141A > C and +405C > G) and 3'-UTR (VEGF +936C > T) were all significantly different in PTC and benign thyroid tissue (p = 0.0001, 0.001 and 0.028 respectively). The VEGF -141 C allele was more common in PTC with lymph node metastases (p = 0.026). VEGF + 405 Galleles andVEGF +936 CC genotype were more common in PTC of advanced pathological staging (p = 0.018 and 0.017 respectively). Also, increased expression of VEGF-A mRNA was noted in PTC compared to control (p = 0.009). Within the group of patients with conventional PTC, those with lymph nodal metastases had a higher level of VEGF-A mRNA expression than other patients (p = 0.0003). CONCLUSION These findings suggest that VEGF polymorphisms and mRNA expression may predict the aggressiveness behaviour of thyroid cancer.

Optimizing I/O cost and managing memory for composition vector method based on correlation matrix calculation in bioinformatics

The generation of a correlation matrix for set of genomic sequences is a common requirement in many bioinformatics problems such as phylogenetic analysis. Each sequence may be millions of bases long and there may be thousands of such sequences which we wish to compare, so not all sequences may fit into main memory at the same time. Each sequence needs to be compared with every other sequence, so we will generally need to page some sequences in and out more than once. In order to minimize execution time we need to minimize this I/O. This paper develops an approach for faster and scalable computing of large-size correlation matrices through the maximal exploitation of available memory and reducing the number of I/O operations. The approach is scalable in the sense that the same algorithms can be executed on different computing platforms with different amounts of memory and can be applied to different bioinformatics problems with different correlation matrix sizes. The significant performance improvement of the approach over previous work is demonstrated through benchmark examples.

Assembly and annotation of a non-model gastropod (Nerita melanotragus) transcriptome: a comparison of De novo assemblers

Background The sequencing, de novo assembly and annotation of transcriptome datasets generated with next generation sequencing (NGS) has enabled biologists to answer genomic questions in non-model species with unprecedented ease. Reliable and accurate de novo assembly and annotation of transcriptomes, however, is a critically important step for transcriptome assemblies generated from short read sequences. Typical benchmarks for assembly and annotation reliability have been performed with model species. To address the reliability and accuracy of de novo transcriptome assembly in non-model species, we generated an RNAseq dataset for an intertidal gastropod mollusc species, Nerita melanotragus, and compared the assembly produced by four different de novo transcriptome assemblers; Velvet, Oases, Geneious and Trinity, for a number of quality metrics and redundancy. Results Transcriptome sequencing on the Ion Torrent PGM™ produced 1,883,624 raw reads with a mean length of 133 base pairs (bp). Both the Trinity and Oases de novo assemblers produced the best assemblies based on all quality metrics including fewer contigs, increased N50 and average contig length and contigs of greater length. Overall the BLAST and annotation success of our assemblies was not high with only 15-19% of contigs assigned a putative function. Conclusions We believe that any improvement in annotation success of gastropod species will require more gastropod genome sequences, but in particular an increase in mollusc protein sequences in public databases. Overall, this paper demonstrates that reliable and accurate de novo transcriptome assemblies can be generated from short read sequencers with the right assembly algorithms. Keywords: Nerita melanotragus; De novo assembly; Transcriptome; Heat shock protein; Ion torrent

Muscle gene electrotransfer is increased by the antioxidant tempol in mice

Electropermeabilization (EP) is an effective method of gene transfer into different tissues. During EP, reactive oxygen species (ROS) are formed, which could affect transfection efficiency. The role of generated ROS and the role of antioxidants in electrotransfer in myoblasts in vitro and in Musculus tibialis cranialis in mice were, therefore, investigated. We demonstrate in the study that during EP of C2C12 myoblasts, ROS are generated on the surface of the cells, which do not induce long-term genomic DNA damage. Plasmid DNA for transfection (pEGFP-N1), which is present outside the cells during EP, neutralizes the generated ROS. The ROS generation is proportional to the amplitude of the electric pulses and can be scavenged by antioxidants, such as vitamin C or tempol. When antioxidants were used during gene electrotransfer, the transfection efficiency of C2C12 myoblasts was statistically significantly increased 1.6-fold with tempol. Also in vivo, the transfection efficiency of M. tibialis cranialis in mice was statistically significantly increased 1.4-fold by tempol. The study indicates that ROS are generated on cells during EP and can be scavenged by antioxidants. Specifically, tempol can be used to improve gene electrotransfer into the muscle and possibly also to other tissues.

Integrin-linked kinase as a target for ERG-mediated invasive properties in prostate cancer models

Approximately half of prostate cancers (PCa) carry TMPRSS2-ERG translocations; however, the clinical impact of this genomic alteration remains enigmatic. Expression of v-ets erythroblastosis virus E26 oncogene like (avian) gene (ERG) promotes prostatic epithelial dysplasia in transgenic mice and acquisition of epithelial-to-mesenchymal transition (EMT) characteristics in human prostatic epithelial cells (PrECs). To explore whether ERG-induced EMT in PrECs was associated with therapeutically targetable transformation characteristics, we established stable populations of BPH-1, PNT1B and RWPE-1 immortalized human PrEC lines that constitutively express flag-tagged ERG3 (fERG). All fERG-expressing populations exhibited characteristics of in vitro and in vivo transformation. Microarray analysis revealed >2000 commonly dysregulated genes in the fERG-PrEC lines. Functional analysis revealed evidence that fERG cells underwent EMT and acquired invasive characteristics. The fERG-induced EMT transcript signature was exemplified by suppressed expression of E-cadherin and keratins 5, 8, 14 and 18; elevated expression of N-cadherin, N-cadherin 2 and vimentin, and of the EMT transcriptional regulators Snail, Zeb1 and Zeb2, and lymphoid enhancer-binding factor-1 (LEF-1). In BPH-1 and RWPE-1-fERG cells, fERG expression is correlated with increased expression of integrin-linked kinase (ILK) and its downstream effectors Snail and LEF-1. Interfering RNA suppression of ERG decreased expression of ILK, Snail and LEF-1, whereas small interfering RNA suppression of ILK did not alter fERG expression. Interfering RNA suppression of ERG or ILK impaired fERG-PrEC Matrigel invasion. Treating fERG-BPH-1 cells with the small molecule ILK inhibitor, QLT-0267, resulted in dose-dependent suppression of Snail and LEF-1 expression, Matrigel invasion and reversion of anchorage-independent growth. These results suggest that ILK is a therapeutically targetable mediator of ERG-induced EMT and transformation in PCa.

The Anadara trapezia transcriptome : a resource for molluscan physiological genomics

In this study we undertook deep sequencing of the blood cockle, Anadara trapezia, transcriptome to generate genomic resources for future functional genomics analyses. Over 27 million high quality paired end reads were assembled into 75 024 contigs. Of these contigs, 29 013 (38.7%) received significant BLASTx hits and gene ontology (GO) terms were assigned to 13 718 of these sequences. This resourcewill facilitate physiological genomic studies to test the gene expression response of A. trapezia to various environmental stresses.

G9a histone methyltransferase contributes to imprinting in the mouse placenta

Whereas DNA methylation is essential for genomic imprinting, the importance of histone methylation in the allelic expression of imprinted genes is unclear. Imprinting control regions (ICRs), however, are marked by histone H3-K9 methylation on their DNA-methylated allele. In the placenta, the paternal silencing along the Kcnq1 domain on distal chromosome 7 also correlates with the presence of H3-K9 methylation, but imprinted repression at these genes is maintained independently of DNA methylation. To explore which histone methyltransferase (HMT) could mediate the allelic H3-K9 methylation on distal chromosome 7, and at ICRs, we generated mouse conceptuses deficient for the SET domain protein G9a. We found that in the embryo and placenta, the differential DNA methylation at ICRs and imprinted genes is maintained in the absence of G9a. Accordingly, in embryos, imprinted gene expression was unchanged at the domains analyzed, in spite of a global loss of H3-K9 dimethylation (H3K9me2). In contrast, the placenta-specific imprinting of genes on distal chromosome 7 is impaired in the absence of G9a, and this correlates with reduced levels of H3K9me2 and H3K9me3. These findings provide the first evidence for the involvement of an HMT and suggest that histone methylation contributes to imprinted gene repression in the trophoblast.