samExploreR: exploring reproducibility and robustness of RNA-seq results based on SAM files

MOTIVATION: Data from RNA-seq experiments provide us with many new possibilities to gain insights into biological and disease mechanisms of cellular functioning. However, the reproducibility and robustness of RNA-seq data analysis results is often unclear. This is in part attributed to the two counter acting goals of (a) a cost efficient and (b) an optimal experimental design leading to a compromise, e.g., in the sequencing depth of experiments.

RESULTS: We introduce an R package called samExploreR that allows the subsampling (m out of n bootstraping) of short-reads based on SAM files facilitating the investigation of sequencing depth related questions for the experimental design. Overall, this provides a systematic way for exploring the reproducibility and robustness of general RNA-seq studies. We exemplify the usage of samExploreR by studying the influence of the sequencing depth and the annotation on the identification of differentially expressed genes.

AVAILABILITY: Availability: samExploreR is available as an R package from Bioconductor (after acceptance of the paper, download link: http://www.bio-complexity.com/samExploreR_1.0.0.tar.gz).

Phylum-wide transcriptome analysis of oogenesis and early embryogenesis in selected nematode species

Oogenesis is a prerequisite for embryogenesis in Metazoa. During both biological processes important decisions must be made to form the embryo and hence ensure the next generation: (1) Maternal gene products (mRNAs, proteins and nutrients) must be supplied to the embryo. (2) Polarity must be established and axes must be specified. While incorporation of maternal gene products occurs during oogenesis, the time point of polarity establishment and axis specification varies among species, as it is accomplished either prior, during, or after fertilisation. But not only the time point when these events take place varies among species but also the underlying mechanisms by which they are triggered. For the nematode model Caenorhabditis elegans the underlying pathways and gene regulatory networks (GRNs) are well understood. It is known that there the sperm entry point initiates a primary polarity in the 1-celled egg and with it the establishment of the anteroposterior axis. However, studies of other nematodes demonstrated that polarity establishment can be independent of sperm entry (Goldstein et al., 1998; Lahl et al., 2006) and that cleavage patterns, symmetry formation and cell specification also differ from C. elegans. In contrast to the studied Chromadorea (more derived nematodes including C. elegans), embryos of some marine Enoplea (more basal representatives) even show no discernible early polarity and blastomeres can adopt variable cell fates (Voronov and Panchin 1998). The underlying pathways controlling the obviously variant embryonic processes in non-Caenorhabditis nematodes are essentially unknown. In this thesis I addressed this issue by performing a detailed unbiased comparative transcriptome analysis based on microarrays and RNA sequencing of selected developmental stages in a variety of nematodes from different phylogenetic branches with C. elegans as a reference system and a nematomorph as an outgroup representative. In addition, I made use of available genomic data to determine the presence or absence of genes for which no expression had been detected. In particular, I focussed on components of selected pathways or GRNs which are known to play essential roles during C. elegans development and/or other invertebrate or vertebrate model systems. Oogenesis must be regulated differently in non-Caenorhabditis nematodes, as crucial controlling components of Wnt and sex determination signaling are absent in these species. In this respect, I identified female-specific expression of potential polarity associated genes during gonad development and oogenesis in the Enoplean nematode Romanomermis culicivorax. I could show that known downstream components of the polarity complexes PAR-3/-6/PKC-3 and PAR-1/-2 are absent in non-Caenorhabditis species. Even PAR-2 as part of the polarity complex does not exist in these nematodes. Instead, transcriptomes of nematodes (including C. elegans), show expression of other polarity-associated complexes such as the Lgl (Lethal giant larvae) complex. This result could pose an alternative route for nematodes and nematomorphs to initiate polarity during early embryogenesis. I could show that crucial pathways of axis specification, such as Wnt and BMP are very different in C. elegans compared to other nematodes. In the former, Wnt signaling, for instance, is mediated by four paralogous beta-catenins, while other Chromadorea have fewer and Enoplea only one beta-catenin. The transcriptomes of R. culicivorax and the nematomorph show that regulators of BMP (e.g. Chordin), are specifically expressed during early embryogenesis only in Enoplea and the close outgroup of nematomorphs. In conclusion, my results demonstrate that the molecular machinery controlling oogenesis and embryogenesis in nematodes is unexpectedly variable and C. elegans cannot be taken as a general model for nematode development. Under this perspective, Enoplean nematodes show more similarities with outgroups than with C. elegans. It appears that certain pathway components were lost or gained during evolution and others adopted new functions. Based on my findings I can conjecture, which pathway components may be ancestral and which were newly acquired in the course of nematode evolution.

Transcriptome-Wide Mapping of Pea Seed Ageing Reveals a Pivotal Role for Genes Related to Oxidative Stress and Programmed Cell Death

Understanding of seed ageing, which leads to viability loss during storage, is vital for ex situ plant conservation and agriculture alike. Yet the potential for regulation at the transcriptional level has not been fully investigated. Here, we studied the relationship between seed viability, gene expression and glutathione redox status during artificial ageing of pea (Pisum sativum) seeds. Transcriptome-wide analysis using microarrays was complemented with qRT-PCR analysis of selected genes and a multilevel analysis of the antioxidant glutathione. Partial degradation of DNA and RNA occurred from the onset of artificial ageing at 60% RH and 50 degrees C, and transcriptome profiling showed that the expression of genes associated with programmed cell death, oxidative stress and protein ubiquitination were altered prior to any sign of viability loss. After 25 days of ageing viability started to decline in conjunction with progressively oxidising cellular conditions, as indicated by a shift of the glutathione redox state towards more positive values (>-190 mV). The unravelling of the molecular basis of seed ageing revealed that transcriptome reprogramming is a key component of the ageing process, which influences the progression of programmed cell death and decline in antioxidant capacity that ultimately lead to seed viability loss.

Messenger ribonucleic acid of cerebral nuclei.

1. RNA was isolated from crude nuclear preparations and from ribosomes derived from rat brain and liver. Nuclear RNA was obtained by lysis of the nuclei with sodium dodecyl sulphate, followed by denaturation and removal of DNA and protein with hot phenol. 2. Base composition analyses indicated that the cerebral nuclear RNA preparation contained a higher proportion of non-ribosomal RNA than the analogous hepatic preparation. 3. Sucrose-density-gradient analyses revealed a heterogeneous profile for each nuclear RNA preparation, with two major peaks possessing the sedimentation properties of ribosomal RNA (18s and 28s). 4. Template activities of both preparations were widely distributed through the sucrose density gradients. 5. The cerebral nuclear RNA preparation was more active than the hepatic nuclear RNA preparation in promoting amino acid incorporation in cell-free systems from Escherichia coli and rat brain. 6. Cerebral nuclear RNA stimulated amino acid incorporation in a cerebral ribosomal system even in the presence of an excess of purified E. coli transfer RNA. 7. It is concluded that a significant proportion of cerebral nuclear RNA has the characteristics of messenger RNA.

Exploring protein-RNA interactions with site-directed mutagenesis and phage display

The importance of RNA as a mediator of genetic information is widely appreciated. RNA molecules also participate in the regulation of various post-transcriptional activities, such as mRNA splicing, editing, RNA stability and transport. Their regulatory roles for these activities are highly dependent on finely tuned associations with cognate proteins. The RNA recognition motif (RRM) is an ancient RNA binding module that participates in hundreds of essential activities where specific RNA recognition is required. We have applied phage display and site-directed mutagenesis to dissect principles of RRM-controlled RNA recognition. The model systems we are investigating are U1A and CUG-BP1. In this dissertation, the molecular basis of the binding affinity of U1A-RNA beyond individual contacts was investigated. We have identified and evaluated the contributions of the local cooperativity formed by three neighboring residues (Asn15, Asn16 and Glu19) to the stability of the U1A-RNA complex. The localized cooperative network was mapped by double-mutant cycles and explored using phage display. We also showed that a cluster of these residues forms a “hot spot” on the surface of U1A; a single substitution at position 19 with Gln or His can alter the binding properties of U1A to recognize a non-cognate G4U RNA. Finally, we applied a deletion analysis of CUG-BP1 to define the contributions of individual RRMs and RRM combinations to the stability of the complex formed between CUG-BP1 and the GRE sequence. The preliminary results showed RRM3 of CUG-BP1 is a key domain for RNA binding. It possibly binds to the GRE sequence cooperatively with RRM2 of CUG-BP1. RRM1 of CUG-BP1 is not required for GRE recognition, but may be important for maintaining the stability of the full-length CUG-BP1.

Whole Genome Analysis of Gene Expression Reveals Coordinated Activation of Signaling and Metabolic Pathways during Pollen-Pistil Interactions in Arabidopsis

Plant reproduction depends on the concerted activation of many genes to ensure correct communication between pollen and pistil. Here, we queried the whole transcriptome of Arabidopsis (Arabidopsis thaliana) in order to identify genes with specific reproductive functions. We used the Affymetrix ATH1 whole genome array to profile wild-type unpollinated pistils and unfertilized ovules. By comparing the expression profile of pistils at 0.5, 3.5, and 8.0 h after pollination and applying a number of statistical and bioinformatics criteria, we found 1,373 genes differentially regulated during pollen-pistil interactions. Robust clustering analysis grouped these genes in 16 time-course clusters representing distinct patterns of regulation. Coregulation within each cluster suggests the presence of distinct genetic pathways, which might be under the control of specific transcriptional regulators. A total of 78% of the regulated genes were expressed initially in unpollinated pistil and/or ovules, 15% were initially detected in the pollen data sets as enriched or preferentially expressed, and 7% were induced upon pollination. Among those, we found a particular enrichment for unknown transcripts predicted to encode secreted proteins or representing signaling and cell wall-related proteins, which may function by remodeling the extracellular matrix or as extracellular signaling molecules. A strict regulatory control in various metabolic pathways suggests that fine-tuning of the biochemical and physiological cellular environment is crucial for reproductive success. Our study provides a unique and detailed temporal and spatial gene expression profile of in vivo pollen-pistil interactions, providing a framework to better understand the basis of the molecular mechanisms operating during the reproductive process in higher plants.

Transcriptome Analysis of Capsicum Chlorosis Virus-Induced Hypersensitive Resistance Response in Bell Capsicum

Background Capsicum chlorosis virus (CaCV) is an emerging pathogen of capsicum, tomato and peanut crops in Australia and South-East Asia. Commercial capsicum cultivars with CaCV resistance are not yet available, but CaCV resistance identified in Capsicum chinense is being introgressed into commercial Bell capsicum. However, our knowledge of the molecular mechanisms leading to the resistance response to CaCV infection is limited. Therefore, transcriptome and expression profiling data provide an important resource to better understand CaCV resistance mechanisms. Methodology/Principal Findings We assembled capsicum transcriptomes and analysed gene expression using Illumina HiSeq platform combined with a tag-based digital gene expression system. Total RNA extracted from CaCV/mock inoculated CaCV resistant (R) and susceptible (S) capsicum at the time point when R line showed a strong hypersensitive response to CaCV infection was used in transcriptome assembly. Gene expression profiles of R and S capsicum in CaCV- and buffer-inoculated conditions were compared. None of the genes were differentially expressed (DE) between R and S cultivars when mock-inoculated, while 2484 genes were DE when inoculated with CaCV. Functional classification revealed that the most highly up-regulated DE genes in R capsicum included pathogenesis-related genes, cell death-associated genes, genes associated with hormone-mediated signalling pathways and genes encoding enzymes involved in synthesis of defense-related secondary metabolites. We selected 15 genes to confirm DE expression levels by real-time quantitative PCR. Conclusion/Significance DE transcript profiling data provided comprehensive gene expression information to gain an understanding of the underlying CaCV resistance mechanisms. Further, we identified candidate CaCV resistance genes in the CaCV-resistant C. annuum x C. chinense breeding line. This knowledge will be useful in future for fine mapping of the CaCV resistance locus and potential genetic engineering of resistance into CaCV-susceptible crops.

Alternate Conformations Regulate Ribosomal Recoding in a Positive-sense RNA Virus

Positive-sense RNA viruses are important animal, plant, insect and bacteria pathogens and constitute the largest group of RNA viruses. Due to the relatively small size of their genomes, these viruses have evolved a variety of non-canonical translation mechanisms to optimize coding capacity expanding their proteome diversity. One such strategy is codon redefinition or recoding. First described in viruses, recoding is a programmed translation event in which codon alterations are context dependent. Recoding takes place in a subset of messenger RNA (mRNAs) with some products reflecting new, and some reflecting standard, meanings. The ratio between the two is both critical and highly regulated. While a variety of recoding mechanisms have been documented, (ribosome shunting, stop-carry on, termination-reinitiation, and translational bypassing), the two most extensively employed by RNA viruses are Programmed Ribosomal Frameshifting (PRF) and Programmed Ribosomal Readthrough (PRT). While both PRT and PRF subvert normal decoding for expression of C-terminal extension products, the former involves an alteration of reading frame, and the latter requires decoding of a non-sense codon. Both processes occur at a low but defined frequency, and both require Recoding Stimulatory Elements (RSE) for regulation and optimum functionality. These stimulatory signals can be embedded in the RNA in the form of sequence or secondary structure, or trans-acting factors outside the mRNA such as proteins or micro RNAs (miRNA). Despite 40+ years of study, the precise mechanisms by which viral RSE mediate ribosome recoding for the synthesis of their proteins, or how the ratio of these products is maintained, is poorly defined. This study reveals that in addition to a long distance RNA:RNA interaction, three alternate conformations and a phylogenetically conserved pseudoknot regulate PRT in the carmovirus Turnip crinkle virus (TCV).

Raman microscopy of formamide-intercalated kaolinites treated by controlled-rate thermal analysis technology

Raman spectroscopy of formamide-intercalated kaolinites treated using controlled-rate thermal analysis technology (CRTA), allowing the separation of adsorbed formamide from intercalated formamide in formamide-intercalated kaolinites, is reported. The Raman spectra of the CRTA-treated formamide-intercalated kaolinites are significantly different from those of the intercalated kaolinites, which display a combination of both intercalated and adsorbed formamide. An intense band is observed at 3629 cm-1, attributed to the inner surface hydroxyls hydrogen bonded to the formamide. Broad bands are observed at 3600 and 3639 cm-1, assigned to the inner surface hydroxyls, which are hydrogen bonded to the adsorbed water molecules. The hydroxyl-stretching band of the inner hydroxyl is observed at 3621 cm-1 in the Raman spectra of the CRTA-treated formamide-intercalated kaolinites. The results of thermal analysis show that the amount of intercalated formamide between the kaolinite layers is independent of the presence of water. Significant differences are observed in the CO stretching region between the adsorbed and intercalated formamide.

Detailed analysis of a conservative difference approximation for the time fractional diffusion equation

Diffusion equations that use time fractional derivatives are attractive because they describe a wealth of problems involving non-Markovian Random walks. The time fractional diffusion equation (TFDE) is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α ∈ (0, 1). Developing numerical methods for solving fractional partial differential equations is a new research field and the theoretical analysis of the numerical methods associated with them is not fully developed. In this paper an explicit conservative difference approximation (ECDA) for TFDE is proposed. We give a detailed analysis for this ECDA and generate discrete models of random walk suitable for simulating random variables whose spatial probability density evolves in time according to this fractional diffusion equation. The stability and convergence of the ECDA for TFDE in a bounded domain are discussed. Finally, some numerical examples are presented to show the application of the present technique.

Determination of preventive maintenance lead time using hybrid analysis

The time for conducting Preventive Maintenance (PM) on an asset is often determined using a predefined alarm limit based on trends of a hazard function. In this paper, the authors propose using both hazard and reliability functions to improve the accuracy of the prediction particularly when the failure characteristic of the asset whole life is modelled using different failure distributions for the different stages of the life of the asset. The proposed method is validated using simulations and case studies.