Regulation of microRNA-1288 in colorectal cancer : altered expression and its clinicopathological significance

We aim to examine the miR-1288 expression in cancer cell lines and a large cohort of patients with colorectal cancer. Two colon cancer cell lines (SW480 and SW48) and one normal colonic epithelial cell line (FHC) were recruited. The miRNA expressions of miR-1288 were tested on these cell lines by using quantitative real-time polymerase chain reaction (qRT-PCR). An exogenous miR-1288 (mimic) was used to detect cell proliferation and cell cycle changes in SW480 using MTT calorimetric assay and flow cytometry, respectively. In addition, tissues from 122 patients with surgical resection of colorectum (82 adenocarcinomas, 20 adenomas, and 20 non-neoplastic tissues) were tested for miR-1288 expression by qRT-PCR. The colon cancer cell lines showed reduced expression of miR-1288 compared to normal colonic epithelial cell line. Over expression of miR-1288 in SW480 cell line showed increased cell proliferation and increased G2-M phase cells. In tissues, reduced miR-1288 expression was noted in majority of colorectal adenocarcinoma compared to colorectal adenoma and non-neoplastic tissues. Reduced or absent expression of miR-1288 was noted in 76% (n = 62/82) of the cancers. The expression levels of miR-1288 were higher in distal colorectal adenocarcinomas (P = 0.013) and in cancers of lower T staging (P = 0.033). To conclude, alternation of miR-1288 expression is important in the progression of colorectal cancer. The differential regulation of miR-1288 was found to be related to cancer location and pathological staging in colorectal cancers.

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.

Altered JS-2 expression in colorectal cancers and its clinical pathological relevance

JS-2 is a novel gene located at 5p15.2 and originally detected in primary oesophageal cancer. There is no study on the role of JS-2 in colorectal cancer. The aim of this study is to determine the gene copy number and expression of JS-2 in a large cohort of patients with colorectal tumours and correlate these to the clinicopathological features of the cancer patients. We evaluated the DNA copy number and mRNA expression of JS-2 in 176 colorectal tissues (116 adenocarcinomas, 30 adenomas and 30 non-neoplastic tissues) using real-time polymerase chain reaction. JS-2 expression was also evaluated in two colorectal cancer cell lines and a benign colorectal cell line. JS-2 amplification was noted in 35% of the colorectal adenocarcinomas. Significant differences in relative expression levels for JS-2 mRNA between different colorectal tissues were noted (p = 0.05). Distal colorectal adenocarcinoma had significantly higher copy number than proximal adenocarcinoma (p = 0.005). The relative expression level of JS-2 was different between colonic and rectal adenocarcinoma (p = 0.007). Mucinous adenocarcinoma showed higher JS-2 expression than non-mucinous adenocarcinoma (p = 0.02). Early T-stage cancers appear to have higher JS-2 copy number and lower expression of JS-2 mRNA than later stage cancers (p = 0.001 and 0.03 respectively). Colorectal cancer cell lines showed lower expression of JS-2 than the benign colorectal cell line. JS-2 copy number change and expression were shown for the first time to be altered in the carcinogenesis of colorectal cancer. In addition, genetic alteration of JS-2 was found to be related to location, pathological subtypes and staging of colorectal cancer.

The Saccharomyces cerevisiae RNA-binding protein Rbp29 functions in cytoplasmic mRNA metabolism

Here we report that the Saccharomyces cerevisiae RBP29 (SGN1, YIR001C) gene encodes a 29-kDa cytoplasmic protein that binds to mRNA in vivo. Rbp29p can be co-immunoprecipitated with the poly(A) tail-binding protein Pab1p from crude yeast extracts in a dosageand RNA-dependent manner. In addition, recombinant Rbp29p binds preferentially to poly(A) with nanomolar binding affinity in vitro. Although RBP29 is not essential for cell viability, its deletion exacerbates the slow growth phenotype of yeast strains harboring mutations in the eIF4G genes TIF4631 and TIF4632. Furthermore, overexpression of RBP29 suppresses the temperaturesensitive growth phenotype of specific tif4631, tif4632, and pab1 alleles. These data suggest that Rbp29p is an mRNA-binding protein that plays a role in modulating the expression of cytoplasmic mRNA.

Yhh1p/Cft1p directly links poly(A) site recognition and RNA polymerase II transcription termination

RNA polymerase II (pol II) transcription termination requires co‐transcriptional recognition of a functional polyadenylation signal, but the molecular mechanisms that transduce this signal to pol II remain unclear. We show that Yhh1p/Cft1p, the yeast homologue of the mammalian AAUAAA interacting protein CPSF 160, is an RNA‐binding protein and provide evidence that it participates in poly(A) site recognition. Interestingly, RNA binding is mediated by a central domain composed of predicted β‐propeller‐forming repeats, which occurs in proteins of diverse cellular functions. We also found that Yhh1p/Cft1p bound specifically to the phosphorylated C‐terminal domain (CTD) of pol II in vitro and in a two‐hybrid test in vivo. Furthermore, transcriptional run‐on analysis demonstrated that yhh1 mutants were defective in transcription termination, suggesting that Yhh1p/Cft1p functions in the coupling of transcription and 3′‐end formation. We propose that direct interactions of Yhh1p/Cft1p with both the RNA transcript and the CTD are required to communicate poly(A) site recognition to elongating pol II to initiate transcription termination.

miRPlant : an integrated tool for identification of plant miRNA from RNA sequencing data

Background Small RNA sequencing is commonly used to identify novel miRNAs and to determine their expression levels in plants. There are several miRNA identification tools for animals such as miRDeep, miRDeep2 and miRDeep*. miRDeep-P was developed to identify plant miRNA using miRDeep’s probabilistic model of miRNA biogenesis, but it depends on several third party tools and lacks a user-friendly interface. The objective of our miRPlant program is to predict novel plant miRNA, while providing a user-friendly interface with improved accuracy of prediction. Result We have developed a user-friendly plant miRNA prediction tool called miRPlant. We show using 16 plant miRNA datasets from four different plant species that miRPlant has at least a 10% improvement in accuracy compared to miRDeep-P, which is the most popular plant miRNA prediction tool. Furthermore, miRPlant uses a Graphical User Interface for data input and output, and identified miRNA are shown with all RNAseq reads in a hairpin diagram. Conclusions We have developed miRPlant which extends miRDeep* to various plant species by adopting suitable strategies to identify hairpin excision regions and hairpin structure filtering for plants. miRPlant does not require any third party tools such as mapping or RNA secondary structure prediction tools. miRPlant is also the first plant miRNA prediction tool that dynamically plots miRNA hairpin structure with small reads for identified novel miRNAs. This feature will enable biologists to visualize novel pre-miRNA structure and the location of small RNA reads relative to the hairpin. Moreover, miRPlant can be easily used by biologists with limited bioinformatics skills.

Identification of direct transcriptional targets of V600EBRAF/MEK signalling in melanoma

Oncogenic mutations in BRAF are common in melanoma and drive constitutive activation of the MEK/ERK pathway. To elucidate the transcriptional events downstream of V600EBRAF/MEK signalling we performed gene expression profiling of A375 melanoma cells treated with potent and selective inhibitors of V600EBRAF and MEK (PLX4720 and PD184352 respectively). Using a stringent Bayesian approach, we identified 69 transcripts that appear to be direct transcriptional targets of this pathway and whose expression changed after 6 h of pathway inhibition. We also identified several additional genes whose expression changed after 24 h of pathway inhibition and which are likely to be indirect transcriptional targets of the pathway. Several of these were confirmed by demonstrating their expression to be similarly regulated when BRAF was depleted using RNA interference, and by using qRT-PCR in other BRAF mutated melanoma lines. Many of these genes are transcription factors and feedback inhibitors of the ERK pathway and are also regulated by MEK signalling in NRAS mutant cells. This study provides a basis for understanding the molecular processes that are regulated by V600EBRAF/MEK signalling in melanoma cells.

Menstrual effluent induces epithelial-mesenchymal transitions in mesothelial cells

Background Menstrual effluent affects mesothelial cell (MC) morphology. We evaluated whether these changes were consistent with epithelial-mesenchymal transitions (EMT). Methods Monolayer cultures of MC were incubated overnight in conditioned media, prepared from cells isolated form menstrual effluent, with or without kinase and ATP inhibitors. Changes in cell morphology were monitored using time-lapse video microscopy and immunohistochemistry. Effects on the expression of EMT-associated molecules were evaluated using real-time RT-PCR and/or Western blot analysis. Results Incubation in conditioned media disrupted cell-cell contacts, and increased MC motility. The changes were reversible. During the changes the distribution of cytokeratins, fibrillar actin and α-tubulin changed. Sodium azide, an inhibitor of ATP production, and Genistein, a general tyrosine kinase inhibitor, antagonized these effects. Wortmannin, a phosphatidylinositol 3-kinase inhibitor, and SU6656, an Src tyrosine kinase inhibitor, only partially antagonized the effect. The expression of Snail and vimentin was markedly up-regulated, whereas the expression of E-cadherin was decreased and cytokeratins were altered. Conclusions In MC, menstrual effluent initiates a reversible, energy-dependent transition process from an epithelial to a mesenchymal phenotype. Involvement of the (Src) tyrosine kinase signalling pathway and the changes in the expression of cytokeratins, Snail, vimentin and E-cadherin demonstrate that the morphological changes are EMT.

High-yeld RNA-extraction method for saliva

BACKGROUND: The use of salivary diagnostics is increasing because of its noninvasiveness, ease of sampling, and the relatively low risk of contracting infectious organisms. Saliva has been used as a biological fluid to identify and validate RNA targets in head and neck cancer patients. The goal of this study was to develop a robust, easy, and cost-effective method for isolating high yields of total RNA from saliva for downstream expression studies. METHODS: Oral whole saliva (200 mu L) was collected from healthy controls (n = 6) and from patients with head and neck cancer (n = 8). The method developed in-house used QIAzol lysis reagent (Qiagen) to extract RNA from saliva (both cell-free supernatants and cell pellets), followed by isopropyl alcohol precipitation, cDNA synthesis, and real-time PCR analyses for the genes encoding beta-actin ("housekeeping" gene) and histatin (a salivary gland-specific gene). RESULTS: The in-house QIAzol lysis reagent produced a high yield of total RNA (0.89 -7.1 mu g) from saliva (cell-free saliva and cell pellet) after DNase treatment. The ratio of the absorbance measured at 260 nm to that at 280 nm ranged from 1.6 to 1.9. The commercial kit produced a 10-fold lower RNA yield. Using our method with the QIAzol lysis reagent, we were also able to isolate RNA from archived saliva samples that had been stored without RNase inhibitors at -80 degrees C for >2 years. CONCLUSIONS: Our in-house QIAzol method is robust, is simple, provides RNA at high yields, and can be implemented to allow saliva transcriptomic studies to be translated into a clinical setting.

A transcriptome resource for the koala (Phascolarctos cinereus) : insights into koala retrovirus transcription and sequence diversity

Background The koala, Phascolarctos cinereus, is a biologically unique and evolutionarily distinct Australian arboreal marsupial. The goal of this study was to sequence the transcriptome from several tissues of two geographically separate koalas, and to create the first comprehensive catalog of annotated transcripts for this species, enabling detailed analysis of the unique attributes of this threatened native marsupial, including infection by the koala retrovirus. Results RNA-Seq data was generated from a range of tissues from one male and one female koala and assembled de novo into transcripts using Velvet-Oases. Transcript abundance in each tissue was estimated. Transcripts were searched for likely protein-coding regions and a non-redundant set of 117,563 putative protein sequences was produced. In similarity searches there were 84,907 (72%) sequences that aligned to at least one sequence in the NCBI nr protein database. The best alignments were to sequences from other marsupials. After applying a reciprocal best hit requirement of koala sequences to those from tammar wallaby, Tasmanian devil and the gray short-tailed opossum, we estimate that our transcriptome dataset represents approximately 15,000 koala genes. The marsupial alignment information was used to look for potential gene duplications and we report evidence for copy number expansion of the alpha amylase gene, and of an aldehyde reductase gene. Koala retrovirus (KoRV) transcripts were detected in the transcriptomes. These were analysed in detail and the structure of the spliced envelope gene transcript was determined. There was appreciable sequence diversity within KoRV, with 233 sites in the KoRV genome showing small insertions/deletions or single nucleotide polymorphisms. Both koalas had sequences from the KoRV-A subtype, but the male koala transcriptome has, in addition, sequences more closely related to the KoRV-B subtype. This is the first report of a KoRV-B-like sequence in a wild population. Conclusions This transcriptomic dataset is a useful resource for molecular genetic studies of the koala, for evolutionary genetic studies of marsupials, for validation and annotation of the koala genome sequence, and for investigation of koala retrovirus. Annotated transcripts can be browsed and queried at http://koalagenome.org

Sexual selection in true fruit flies (Diptera: Tephritidae) : transcriptome and experimental evidences for phytochemicals increasing male competitive ability

In male tephritid fruit flies of the genus Bactrocera, feeding on secondary plant compounds (sensu lato male lures = methyl eugenol, raspberry ketone and zingerone) increases male mating success. Ingested male lures alter the male pheromonal blend, normally making it more attractive to females and this is considered the primary mechanism for the enhanced mating success. However, the male lures raspberry ketone and zingerone are known, across a diverse range of other organisms, to be involved in increasing energy metabolism. If this also occurs in Bactrocera, then this may represent an additional benefit to males as courtship is metabolically expensive and lure feeding may increase a fly's short-term energy. We tested this hypothesis by performing comparative RNA-seq analysis between zingerone-fed and unfed males of Bactrocera tryoni. We also carried out behavioural assays with zingerone- and cuelure-fed males to test whether they became more active. RNA-seq analysis revealed, in zingerone-fed flies, up-regulation of 3183 genes with homologues transcripts to those known to regulate intermale aggression, pheromone synthesis, mating and accessory gland proteins, along with significant enrichment of several energy metabolic pathways and gene ontology terms. Behavioural assays show significant increases in locomotor activity, weight reduction and successful mating after mounting; all direct/indirect measures of increased activity. These results suggest that feeding on lures leads to complex physiological changes, which result in more competitive males. These results do not negate the pheromone effect, but do strongly suggest that the phytochemical-induced sexual selection is governed by both female preference and male competitive mechanisms.

Sub-genomic RNA of defective interfering (D.I.) dengue viral particles is replicated in the same manner as full length genomes

The predicted secondary structure of sub-genomic RNA in dengue virus defective interfering (D.I.) particles from patients, or generated in vitro, resembled that of the 3′ and 5′ regions of wild type dengue virus (DENV) genomes. While these structures in the sub-genomic RNA were found to be essential for its replication, their nucleotide sequences were not, so long as any new sequences maintained wild type RNA secondary structure. These observations suggested that these sub-genomic fragments of RNA from dengue viruses were replicated in the same manner as the full length genomes of their wild type, “helper”, viruses and that they probably represent the smallest fragments of DENV RNA that can be replicated during a natural infection. While D.I. particles containing sub-genomic RNA are completely parasitic, the relationship between wild type and D.I. DENV may be symbiotic, with the D.I. particles enhancing the transmission of infectious DENV.

An evaluation of potential candidate genes involved in salinity tolerance in striped catfish (Pangasianodon hypopthalmus) using an RNA-SEQ approach

The project investigated the molecular response of Tra catfish (Pangasianodon hypophthalmus) to elevated salinity conditions. We employed Next generation sequencing platform to evaluate differential gene expression profiles of key genes under two salinity conditions. Results of the current project can form the basis for further studies to confirm the functional roles of specific genes that influence salinity tolerance in the target species and more broadly in other freshwater teleost fishes. Ultimately, the approach can contribute to developing superior culture stocks of the target species.

Expression and characterization of digestive enzyme genes from hepatopancreatic transcripts from redclaw crayfish (Cherax quadricarinatus)

Redclaw crayfish, Cherax quadricarinatus, possess a number of biological and commercial attributes that make them ideal for commercial aquaculture. While some studies have investigated digestive enzyme activity and nutritional requirements of this species, little information exists about the expression of digestive enzyme genes and their role in regulating digestive capacity. The current study therefore sequenced and annotated a RNASeq library constructed from a redclaw hepatopancreas to identify genes involved in digestive enzyme production. We observed that most of the transcripts that were annotated as digestive enzyme genes are associated with carbohydrate metabolism, thus confirming that redclaw have an innate capacity to digest a range of carbohydrate substrates. While endoglucanases were the most abundant group of digestive enzymes found, a number of novel transcripts were also detected. Here, we provide the first report for the presence and expression of endo-b-mannanase in freshwater crayfish. This novel gene showed significant alignment with a GH5 family protein from marine Limnoriids, wood borers that do not possess symbiotic microbes in their gut system. Overall, the data generated here provide an important resource to better understand the suite of digestive enzymes in redclaw that are very useful to fully utilize the species’ digestive capacity and will assist development of specific artificial feeds.

RNASeqBrowser: A genome browser for simultaneous visualization of raw strand specific RNAseq reads and UCSC genome browser custom tracks

Background Strand specific RNAseq data is now more common in RNAseq projects. Visualizing RNAseq data has become an important matter in Analysis of sequencing data. The most widely used visualization tool is the UCSC genome browser that introduced the custom track concept that enabled researchers to simultaneously visualize gene expression at a particular locus from multiple experiments. Our objective of the software tool is to provide friendly interface for visualization of RNAseq datasets. Results This paper introduces a visualization tool (RNASeqBrowser) that incorporates and extends the functionality of the UCSC genome browser. For example, RNASeqBrowser simultaneously displays read coverage, SNPs, InDels and raw read tracks with other BED and wiggle tracks -- all being dynamically built from the BAM file. Paired reads are also connected in the browser to enable easier identification of novel exon/intron borders and chimaeric transcripts. Strand specific RNAseq data is also supported by RNASeqBrowser that displays reads above (positive strand transcript) or below (negative strand transcripts) a central line. Finally, RNASeqBrowser was designed for ease of use for users with few bioinformatic skills, and incorporates the features of many genome browsers into one platform. Conclusions The features of RNASeqBrowser: (1) RNASeqBrowser integrates UCSC genome browser and NGS visualization tools such as IGV. It extends the functionality of the UCSC genome browser by adding several new types of tracks to show NGS data such as individual raw reads, SNPs and InDels. (2) RNASeqBrowser can dynamically generate RNA secondary structure. It is useful for identifying non-coding RNA such as miRNA. (3) Overlaying NGS wiggle data is helpful in displaying differential expression and is simple to implement in RNASeqBrowser. (4) NGS data accumulates a lot of raw reads. Thus, RNASeqBrowser collapses exact duplicate reads to reduce visualization space. Normal PC’s can show many windows of NGS individual raw reads without much delay. (5) Multiple popup windows of individual raw reads provide users with more viewing space. This avoids existing approaches (such as IGV) which squeeze all raw reads into one window. This will be helpful for visualizing multiple datasets simultaneously. RNASeqBrowser and its manual are freely available at http://www.australianprostatecentre.org/research/software/rnaseqbrowser webcite or http://sourceforge.net/projects/rnaseqbrowser/ webcite