A crystallographic and modelling study of a human telomeric RNA (TERRA) quadruplex

DNA telomeric repeats in mammalian cells are transcribed to guanine-rich RNA sequences, which adopt parallel-stranded G-quadruplexes with a propeller-like fold. The successful crystallization and structure analysis of a bimolecular human telomeric RNA G-quadruplex, folded into the same crystalline environment as an equivalent DNA oligonucleotide sequence, is reported here. The structural basis of the increased stability of RNA telomeric quadruplexes over DNA ones and their preference for parallel topologies is described here. Our findings suggest that the 2'-OH hydroxyl groups in the RNA quadruplex play a significant role in redefining hydration structure in the grooves and the hydrogen bonding networks. The preference for specific nucleotides to populate the C3'-endo sugar pucker domain is accommodated by alterations in the phosphate backbone, which leads to greater stability through enhanced hydrogen bonding networks. Molecular dynamics simulations on the DNA and RNA quadruplexes are consistent with these findings. The computations, based on the native crystal structure, provide an explanation for RNA G-quadruplex ligand binding selectivity for a group of naphthalene diimide ligands as compared to the DNA G-quadruplex.

Detecting sequence dependent transcriptional pauses from RNA and protein number time series

Background: Evidence suggests that in prokaryotes sequence-dependent transcriptional pauses a?ect the dynamics of transcription and translation, as well as of small genetic circuits. So far, a few pause-prone sequences have been identi?ed from in vitro measurements of transcription elongation kinetics.

Results: Using a stochastic model of gene expression at the nucleotide and codon levels with realistic parameter values, we investigate three di?erent but related questions and present statistical methods for their analysis. First, we show that information from in vivo RNA and protein temporal numbers is su?cient to discriminate between models with and without a pause site in their coding sequence. Second, we demonstrate that it is possible to separate a large variety of models from each other with pauses of various durations and locations in the template by means of a hierarchical clustering and a random forest classi?er. Third, we introduce an approximate likelihood function that allows to estimate the location of a pause site.

Conclusions: This method can aid in detecting unknown pause-prone sequences from temporal measurements of RNA and protein numbers at a genome-wide scale and thus elucidate possible roles that these sequences play in the dynamics of genetic networks and phenotype.

Stable isotope probing analysis of the influence of liming on root exudate utilization by soil microorganisms

Rhizosphere microorganisms play an important role in soil carbon flow, through turnover of root exudates, but there is little information on which organisms are actively involved or on the influence of environmental conditions on active communities. In this study, a (CO2)-C-13 pulse labelling field experiment was performed in an upland grassland soil, followed by RNA-stable isotope probing (SIP) analysis, to determine the effect of liming on the structure of the rhizosphere microbial community metabolizing root exudates. The lower limit of detection for SIP was determined in soil samples inoculated with a range of concentrations of C-13-labelled Pseudomonas fluorescens and was found to lie between 10(5) and 10(6) cells per gram of soil. The technique was capable of detecting microbial communities actively assimilating root exudates derived from recent photo-assimilate in the field. Denaturing gradient gel electrophoresis (DGGE) profiles of bacteria, archaea and fungi derived from fractions obtained from caesium trifluoroacetate (CsTFA) density gradient ultracentrifugation indicated that active communities in limed soils were more complex than those in unlimed soils and were more active in utilization of recently exuded C-13 compounds. In limed soils, the majority of the community detected by standard RNA-DGGE analysis appeared to be utilizing root exudates. In unlimed soils, DGGE profiles from C-12 and C-13 RNA fractions differed, suggesting that a proportion of the active community was utilizing other sources of organic carbon. These differences may reflect differences in the amount of root exudation under the different conditions.

Rice-arsenate interactions in hydroponics:whole genome transcriptional analysis

Rice (Oryza sativa) varieties that are arsenate-tolerant (Bala) and -sensitive (Azucena) were used to conduct a transcriptome analysis of the response of rice seedlings to sodium arsenate (AsV) in hydroponic solution. RNA extracted from the roots of three replicate experiments of plants grown for 1 week in phosphate-free nutrient with or without 13.3 muM AsV was used to challenge the Affymetrix (52K) GeneChip Rice Genome array. A total of 576 probe sets were significantly up-regulated at least 2-fold in both varieties, whereas 622 were down-regulated. Ontological classification is presented. As expected, a large number of transcription factors, stress proteins, and transporters demonstrated differential expression. Striking is the lack of response of classic oxidative stress-responsive genes or phytochelatin synthases/synthatases. However, the large number of responses from genes involved in glutathione synthesis, metabolism, and transport suggests that glutathione conjugation and arsenate methylation may be important biochemical responses to arsenate challenge. In this report, no attempt is made to dissect differences in the response of the tolerant and sensitive variety, but analysis in a companion article will link gene expression to the known tolerance loci available in the BalaxAzucena mapping population.

Resolving candidate genes of mouse skeletal muscle QTL via RNA-Seq and expression network analyses

BACKGROUND:

We have recently identified a number of Quantitative Trait Loci (QTL) contributing to the 2-fold muscle weight difference between the LG/J and SM/J mouse strains and refined their confidence intervals. To facilitate nomination of the candidate genes responsible for these differences we examined the transcriptome of the tibialis anterior (TA) muscle of each strain by RNA-Seq.

13,726 genes were expressed in mouse skeletal muscle. Intersection of a set of 1061 differentially expressed transcripts with a mouse muscle Bayesian Network identified a coherent set of differentially expressed genes that we term the LG/J and SM/J Regulatory Network (LSRN). The integration of the QTL, transcriptome and the network analyses identified eight key drivers of the LSRN (Kdr, Plbd1, Mgp, Fah, Prss23, 2310014F06Rik, Grtp1, Stk10) residing within five QTL regions, which were either polymorphic or differentially expressed between the two strains and are strong candidates for quantitative trait genes (QTGs) underlying muscle mass. The insight gained from network analysis including the ability to make testable predictions is illustrated by annotating the LSRN with knowledge-based signatures and showing that the SM/J state of the network corresponds to a more oxidative state. We validated this prediction by NADH tetrazolium reductase staining in the TA muscle revealing higher oxidative potential of the SM/J compared to the LG/J strain (p<0.03).

Thus, integration of fine resolution QTL mapping, RNA-Seq transcriptome information and mouse muscle Bayesian Network analysis provides a novel and unbiased strategy for nomination of muscle QTGs.

Preparation of planar retinal specimens:verification by histology, mRNA profiling, and proteome analysis

Elucidation of the transcriptome and proteome of the normal retina will be difficult since it is comprised of at least 55 different cell types. However the characteristic layered cellular anatomy of the retina makes it amenable to planar sectioning, enabling the generation of enriched retinal cell populations. The aim of this study was to validate a reproducible method for preparing enriched retinal layers from porcine retina.

Automated tumour annotation and analysis for molecular pathology using TissueMark (R)

Objective: Molecular pathology relies on identifying anomalies using PCR or analysis of DNA/RNA. This is important in solid tumours where molecular stratification of patients define targeted treatment. These molecular biomarkers rely on examination of tumour, annotation for possible macro dissection/tumour cell enrichment and the estimation of % tumour. Manually marking up tumour is error prone. Method: We have developed a method for automated tumour mark-up and % cell calculations using image analysis called TissueMark® based on texture analysis for lung, colorectal and breast (cases=245, 100, 100 respectively). Pathologists marked slides for tumour and reviewed the automated analysis. A subset of slides was manually counted for tumour cells to provide a benchmark for automated image analysis. Results: There was a strong concordance between pathological and automated mark-up (100 % acceptance rate for macro-dissection). We also showed a strong concordance between manually/automatic drawn boundaries (median exclusion/inclusion error of 91.70 %/89 %). EGFR mutation analysis was precisely the same for manual and automated annotation-based macrodissection. The annotation accuracy rates in breast and colorectal cancer were 83 and 80 % respectively. Finally, region-based estimations of tumour percentage using image analysis showed significant correlation with actual cell counts. Conclusion: Image analysis can be used for macro-dissection to (i) annotate tissue for tumour and (ii) estimate the % tumour cells and represents an approach to standardising/improving molecular diagnostics.

Transcriptome analysis of a parasitic clade V nematode:comparative analysis of potential molecular anthelmintic targets in Cylicostephanus goldi

Clade V nematodes comprise several parasitic species that include the cyathostomins, primary helminth pathogens of horses. Next generation transcriptome datasets are available for eight parasitic clade V nematodes, although no equine parasites are included in this group. Here, we report next generation transcriptome sequencing analysis for the common cyathostomin species, Cylicostephanus goldi. A cDNA library was generated from RNA extracted from 17 C. goldi male and female adult parasites. Following sequencing using a 454 GS FLX pyrosequencer, a total of 475,215 sequencing reads were generated, which were assembled into 26,910 contigs. Using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases, 27% of the transcriptome was annotated. Further in-depth analysis was carried out by comparing the C. goldi dataset with the next generation transcriptomes and genomes of other clade V nematodes, with the Oesophagostomum dentatum transcriptome and the Haemonchus contortus genome showing the highest levels of sequence identity with the cyathostomin dataset (45%). The C. goldi transcriptome was mined for genes associated with anthelmintic mode of action and/or resistance. Sequences encoding proteins previously associated with the three major anthelmintic classes used in horses were identified, with the exception of the P-glycoprotein group. Targeted resequencing of the glutamate gated chloride channel α4 subunit (glc-3), one of the primary targets of the macrocyclic lactone anthelmintics, was performed for several cyathostomin species. We believe this study reports the first transcriptome dataset for an equine helminth parasite, providing the opportunity for in-depth analysis of these important parasites at the molecular level. Sequences encoding enzymes involved in key processes and genes associated with levamisole/pyrantel and macrocyclic lactone resistance, in particular the glutamate gated chloride channels, were identified. This novel data will inform cyathostomin biology and anthelmintic resistance studies in future.

Prevalence of germ-line mutations in p16, p19ARF, and CDK4 in familial melanoma:analysis of a clinic-based population

Five to ten percent of individuals with melanoma have another affected family member, suggesting familial predisposition. Germ-line mutations in the cyclin-dependent kinase (CDK) inhibitor p16 have been reported in a subset of melanoma pedigrees, but their prevalence is unknown in more common cases of familial melanoma that do not involve large families with multiple affected members. We screened for germ-line mutations in p16 and in two other candidate melanoma genes, p19ARF and CDK4, in 33 consecutive patients treated for melanoma; these patients had at least one affected first or second degree relative (28 independent families). Five independent, definitive p16 mutations were detected (18%, 95% confidence interval: 6%, 37%), including one nonsense, one disease-associated missense, and three small deletions. No mutations were detected in CDK4. Disease-associated mutations in p19ARF, whose transcript is derived in part from an alternative codon reading frame of p16, were only detected in patients who also had mutations inactivating p16. We conclude that germ-line p16 mutations are present in a significant fraction of individuals who have melanoma and a positive family history.

PTEN mRNA detection by chromogenic, RNA in situ technologies: a reliable alternative to PTEN immunohistochemistry

Immunohistochemical staining for phosphatase and tensin homolog (PTEN) does not have either an acceptable standard protocol or concordance of scoring between pathologists. Evaluation of PTEN mRNA with a unique and verified sequence probe may offer a realistic alternative providing a robust and reproducible protocol. In this study, we have evaluated an in situ hybridization (ISH) protocol for PTEN mRNA using RNAScope technology and compared it with a standard protocol for PTEN immunohistochemistry (IHC). PTEN mRNA expression by ISH was consistently more sensitive than PTEN IHC, with 56% of samples on a mixed-tumor tissue microarray (TMA) showing high expression by ISH compared with 42% by IHC. On a prostate TMA, 49% of cases showed high expression by ISH compared with 43% by IHC. Variations in PTEN mRNA expression within malignant epithelium were quantifiable using image analysis on the prostate TMAs. Within tumors, clear overexpression of PTEN mRNA on malignant epithelium compared with benign epithelium was frequently observed and quantified. The use of SpotStudio software in the mixed-tumor TMA allowed for clear demonstration of varying levels of PTEN mRNA between tumor samples by the mRNA methodology. This was evident by the quantifiable differences between distinct oropharyngeal tumors (up to 3-fold increase in average number of spots per cell between 2 cases). mRNA detection of PTEN or other biomarkers, for which optimal or standardized immunohistochemical techniques are not available, represents a means by which heterogeneity of expression within focal regions of tumor can be explored with more confidence.

Meta-analysis of prostate cancer gene expression data identifies a novel discriminatory signature enriched for glycosylating enzymes

BACKGROUND: Tumorigenesis is characterised by changes in transcriptional control. Extensive transcript expression data have been acquired over the last decade and used to classify prostate cancers. Prostate cancer is, however, a heterogeneous multifocal cancer and this poses challenges in identifying robust transcript biomarkers.

METHODS: In this study, we have undertaken a meta-analysis of publicly available transcriptomic data spanning datasets and technologies from the last decade and encompassing laser capture microdissected and macrodissected sample sets.

RESULTS: We identified a 33 gene signature that can discriminate between benign tissue controls and localised prostate cancers irrespective of detection platform or dissection status. These genes were significantly overexpressed in localised prostate cancer versus benign tissue in at least three datasets within the Oncomine Compendium of Expression Array Data. In addition, they were also overexpressed in a recent exon-array dataset as well a prostate cancer RNA-seq dataset generated as part of the The Cancer Genomics Atlas (TCGA) initiative. Biologically, glycosylation was the single enriched process associated with this 33 gene signature, encompassing four glycosylating enzymes. We went on to evaluate the performance of this signature against three individual markers of prostate cancer, v-ets avian erythroblastosis virus E26 oncogene homolog (ERG) expression, prostate specific antigen (PSA) expression and androgen receptor (AR) expression in an additional independent dataset. Our signature had greater discriminatory power than these markers both for localised cancer and metastatic disease relative to benign tissue, or in the case of metastasis, also localised prostate cancer.

CONCLUSION: In conclusion, robust transcript biomarkers are present within datasets assembled over many years and cohorts and our study provides both examples and a strategy for refining and comparing datasets to obtain additional markers as more data are generated.

Comprehensive molecular pathology analysis of small bowel adenocarcinoma reveals novel targets with potential for clinical utility

Small bowel accounts for only 0.5% of cancer cases in the US but incidence rates have been rising at 2.4% per year over the past decade. One-third of these are adenocarcinomas but little is known about their molecular pathology and no molecular markers are available for clinical use. Using a retrospective 28 patient matched normal-tumor cohort, next-generation sequencing, gene expression arrays and CpG methylation arrays were used for molecular profiling. Next-generation sequencing identified novel mutations in IDH1, CDH1, KIT, FGFR2, FLT3, NPM1, PTEN, MET, AKT1, RET, NOTCH1 and ERBB4. Array data revealed 17% of CpGs and 5% of RNA transcripts assayed to be differentially methylated and expressed respectively (p < 0.01). Merging gene expression and DNA methylation data revealed CHN2 as consistently hypermethylated and downregulated in this disease (Spearman -0.71, p < 0.001). Mutations in TP53 which were found in more than half of the cohort (15/28) and Kazald1 hypomethylation were both were indicative of poor survival (p = 0.03, HR = 3.2 and p = 0.01, HR = 4.9 respectively). By integrating high-throughput mutational, gene expression and DNA methylation data, this study reveals for the first time the distinct molecular profile of small bowel adenocarcinoma and highlights potential clinically exploitable markers.

Distribution, occurrence and biotoxin composition of the main shellfish toxin producing microalgae within European waters: A comparison of methods of analysis

Harmful algal blooms (HABs) are a natural global phenomena emerging in severity and extent. Incidents have many economic, ecological and human health impacts. Monitoring and providing early warning of toxic HABs are critical for protecting public health. Current monitoring programmes include measuring the number of toxic phytoplankton cells in the water and biotoxin levels in shellfish tissue. As these efforts are demanding and labour intensive, methods which improve the efficiency are essential. This study compares the utilisation of a multitoxin surface plasmon resonance (multitoxin SPR) biosensor with enzyme-linked immunosorbent assay (ELISA) and analytical methods such as high performance liquid chromatography with fluorescence detection (HPLC-FLD) and liquid chromatography-tandem mass spectrometry (LC-MS/MS) for toxic HAB monitoring efforts in Europe. Seawater samples (n = 256) from European waters, collected 2009-2011, were analysed for biotoxins: saxitoxin and analogues, okadaic acid and dinophysistoxins 1/2 (DTX1/DTX2) and domoic acid responsible for paralytic shellfish poisoning (PSP), diarrheic shellfish poisoning (DSP) and amnesic shellfish poisoning (ASP), respectively. Biotoxins were detected mainly in samples from Spain and Ireland. France and Norway appeared to have the lowest number of toxic samples. Both the multitoxin SPR biosensor and the RNA microarray were more sensitive at detecting toxic HABs than standard light microscopy phytoplankton monitoring. Correlations between each of the detection methods were performed with the overall agreement, based on statistical 2 × 2 comparison tables, between each testing platform ranging between 32% and 74% for all three toxin families illustrating that one individual testing method may not be an ideal solution. An efficient early warning monitoring system for the detection of toxic HABs could therefore be achieved by combining both the multitoxin SPR biosensor and RNA microarray.

Phylogenetic insights of Avena genus based on genomic analysis of 45S rDNA molecular organization

Tese de mestrado. Biologia (Biologia Molecular e Genética). Universidade de Lisboa, Faculdade de Ciências, 2014

RT-Bst: an integrated approach for reverse transcription and enrichment of cDNA from viral RNA

The synthesis of cDNA from RNA is challenging due to the inefficiency of reverse transcription (RT). In order to address this, a method was developed known as RT-Bst for sequential RT of RNA and Bst DNA polymerase amplification for enrichment of cDNA in a single tube reaction. Using genomic RNA from bacteriophage MS2, the yield of cDNA produced by RT alone and RT-Bst were compared by analysis of PCR-amplified products. Using random primers a superior performance was observed when amplifying MS2 RNA following RT-Bst compared to RT alone, indicating that greater quantities of cDNA were present after RT-Bst. RT-Bst was also compared with RT alone for their relative ability to produce sufficient cDNA to amplify 8 target regions spanning the respiratory syncytial virus (RSV) genome. Six out of 8 targets were amplified consistently by PCR subsequent to RT-Bst amplification whereas only 3 out of 8 targets could be amplified after RT alone. RSV sequences were selectively amplified using RSV specific primers from a mixed template containing an excess of MS2 RNA in a RT-Bst reaction without amplifying MS2 sequences. This suggests that RT-Bst can be used to amplify RNA sequences non-specifically using random primers and specifically using sequence specific primers and enhances the yield of cDNA when compared to RT alone.