GSVA: gene set variation analysis for microarray and RNA-seq data

Gene set enrichment (GSE) analysis is a popular framework for condensing information from gene expression profiles into a pathway or signature summary. The strengths of this approach over single gene analysis include noise and dimension reduction, as well as greater biological interpretability. As molecular profiling experiments move beyond simple case-control studies, robust and flexible GSE methodologies are needed that can model pathway activity within highly heterogeneous data sets. To address this challenge, we introduce Gene Set Variation Analysis (GSVA), a GSE method that estimates variation of pathway activity over a sample population in an unsupervised manner. We demonstrate the robustness of GSVA in a comparison with current state of the art sample-wise enrichment methods. Further, we provide examples of its utility in differential pathway activity and survival analysis. Lastly, we show how GSVA works analogously with data from both microarray and RNA-seq experiments. GSVA provides increased power to detect subtle pathway activity changes over a sample population in comparison to corresponding methods. While GSE methods are generally regarded as end points of a bioinformatic analysis, GSVA constitutes a starting point to build pathway-centric models of biology. Moreover, GSVA contributes to the current need of GSE methods for RNA-seq data. GSVA is an open source software package for R which forms part of the Bioconductor project and can be downloaded at http://www.bioconductor.org.

Enhanced Botrytis cinerea resistance of Arabidopsis plants grown in compost may be explained by increased expression of defense-related genes, as revealed by microarray analysis

Composts are the products obtained after the aerobic degradation of different types of organic matter waste and can be used as substrates or substrate/soil amendments for plant cultivation. There is a small but increasing number of reports that suggest that foliar diseases may be reduced when using compost, rather than standard substrates, as growing medium. The purpose of this study was to examine the gene expression alteration produced by the compost to gain knowledge of the mechanisms involved in compost-induced systemic resistance. A compost from olive marc and olive tree leaves was able to induce resistance against Botrytis cinerea in Arabidopsis, unlike the standard substrate, perlite. Microarray analyses revealed that 178 genes were differently expressed, with a fold change cut-off of 1, of which 155 were up-regulated and 23 were down-regulated in compost-grown, as against perlite-grown plants. A functional enrichment study of up-regulated genes revealed that 38 Gene Ontology terms were significantly enriched. Response to stress, biotic stimulus, other organism, bacterium, fungus, chemical and abiotic stimulus, SA and ABA stimulus, oxidative stress, water, temperature and cold were significantly enriched, as were immune and defense responses, systemic acquired resistance, secondary metabolic process and oxireductase activity. Interestingly, PR1 expression, which was equally enhanced by growing the plants in compost and by B. cinerea inoculation, was further boosted in compost-grown pathogen-inoculated plants. Compost triggered a plant response that shares similarities with both systemic acquired resistance and ABA-dependent/independent abiotic stress responses.

Characterization of HbpR binding by site-directed mutagenesis of its DNA-binding site and by deletion of the effector domain.

In the presence of 2-hydroxybiphenyl, the enhancer binding protein, HbpR, activates the sigma54-dependent P(hbpC) promoter and controls the initial steps of 2-hydroxybiphenyl degradation in Pseudomonas azelaica. In the activation process, an oligomeric HbpR complex of unknown subunit composition binds to an operator region containing two imperfect palindromic sequences. Here, the HbpR-DNA binding interactions were investigated by site-directed mutagenesis of the operator region and by DNA-binding assays using purified HbpR. Mutations that disrupted the twofold symmetry in the palindromes did not affect the binding affinity of HbpR, but various mutations along a 60 bp region, and also outside the direct palindromic sequences, decreased the binding affinity. Footprints of HbpR on mutant operator fragments showed that a partial loss of binding contacts occurs, suggesting that the binding of one HbpR 'protomer' in the oligomeric complex is impaired whilst leaving the other contacts intact. An HbpR variant, devoid of its N-terminal sensing A-domain, was unable to activate transcription from the hbpC promoter while maintaining protection of the operator DNA in footprints. Wild-type HbpR was unable to activate transcription from the hbpC promoter when delta A-HbpR was expressed in the same cell, suggesting the formation of (repressing) hetero-oligomers. This model implies that HbpR can self-associate on its operator DNA without effector recognition or ATP binding. Furthermore, our findings suggest that the N-terminal sensing domain of HbpR is needed to activate the central ATPase domain rather than to repress a constitutively active C domain, as is the case for the related regulatory protein XylR.

Saccharomyces cerevisiae Dmc1 and Rad51 Proteins Preferentially Function with Tid1 and Rad54 Proteins, Respectively, to Promote DNA Strand Invasion during Genetic Recombination.

The Saccharomyces cerevisiae Dmc1 and Tid1 proteins are required for the pairing of homologous chromosomes during meiotic recombination. This pairing is the precursor to the formation of crossovers between homologs, an event that is necessary for the accurate segregation of chromosomes. Failure to form crossovers can have serious consequences and may lead to chromosomal imbalance. Dmc1, a meiosis-specific paralog of Rad51, mediates the pairing of homologous chromosomes. Tid1, a Rad54 paralog, although not meiosis-specific, interacts with Dmc1 and promotes crossover formation between homologs. In this study, we show that purified Dmc1 and Tid1 interact physically and functionally. Dmc1 forms stable nucleoprotein filaments that can mediate DNA strand invasion. Tid1 stimulates Dmc1-mediated formation of joint molecules. Under conditions optimal for Dmc1 reactions, Rad51 is specifically stimulated by Rad54, establishing that Dmc1-Tid1 and Rad51-Rad54 function as specific pairs. Physical interaction studies show that specificity in function is not dictated by direct interactions between the proteins. Our data are consistent with the hypothesis that Rad51-Rad54 function together to promote intersister DNA strand exchange, whereas Dmc1-Tid1 tilt the bias toward interhomolog DNA strand exchange.

Phylogenetics of Olea (Oleaceae) based on plastid and nuclear ribosomal DNA sequences: tertiary climatic shifts and lineage differentiation times.

BACKGROUND AND AIMS: The genus Olea (Oleaceae) includes approx. 40 taxa of evergreen shrubs and trees classified in three subgenera, Olea, Paniculatae and Tetrapilus, the first of which has two sections (Olea and Ligustroides). Olive trees (the O. europaea complex) have been the subject of intensive research, whereas little is known about the phylogenetic relationships among the other species. To clarify the biogeographical history of this group, a molecular analysis of Olea and related genera of Oleaceae is thus necessary. METHODS: A phylogeny was built of Olea and related genera based on sequences of the nuclear ribosomal internal transcribed spacer-1 and four plastid regions. Lineage divergence and the evolution of abaxial peltate scales, the latter character linked to drought adaptation, were dated using a Bayesian method. KEY RESULTS: Olea is polyphyletic, with O. ambrensis and subgenus Tetrapilus not sharing a most recent common ancestor with the main Olea clade. Partial incongruence between nuclear and plastid phylogenetic reconstructions suggests a reticulation process in the evolution of subgenus Olea. Estimates of divergence times for major groups of Olea during the Tertiary were obtained. CONCLUSIONS: This study indicates the necessity of revising current taxonomic boundaries in Olea. The results also suggest that main lines of evolution were promoted by major Tertiary climatic shifts: (1) the split between subgenera Olea and Paniculatae appears to have taken place at the Miocene-Oligocene boundary; (2) the separation of sections Ligustroides and Olea may have occurred during the Early Miocene following the Mi-1 glaciation; and (3) the diversification within these sections (and the origin of dense abaxial indumentum in section Olea) was concomitant with the aridification of Africa in the Late Miocene.

Formation of knots in partially replicated DNA molecules.

Bacterial plasmids with two origins of replication in convergent orientation are frequently knotted in vivo. The knots formed are localised within the newly replicated DNA regions. Here, we analyse DNA knots tied within replication bubbles of such plasmids, and observe that the knots formed show predominantly positive signs of crossings. We propose that helical winding of replication bubbles in vivo leads to topoisomerase-mediated formation of knots on partially replicated DNA molecules.

Electrophoretic mobility of supercoiled, catenated and knotted DNA molecules.

We systematically varied conditions of two-dimensional (2D) agarose gel electrophoresis to optimize separation of DNA topoisomers that differ either by the extent of knotting, the extent of catenation or the extent of supercoiling. To this aim we compared electrophoretic behavior of three different families of DNA topoisomers: (i) supercoiled DNA molecules, where supercoiling covered the range extending from covalently closed relaxed up to naturally supercoiled DNA molecules; (ii) postreplicative catenanes with catenation number increasing from 1 to ∼15, where both catenated rings were nicked; (iii) knotted but nicked DNA molecules with a naturally arising spectrum of knots. For better comparison, we studied topoisomer families where each member had the same total molecular mass. For knotted and supercoiled molecules, we analyzed dimeric plasmids whereas catenanes were composed of monomeric forms of the same plasmid. We observed that catenated, knotted and supercoiled families of topoisomers showed different reactions to changes of agarose concentration and voltage during electrophoresis. These differences permitted us to optimize conditions for their separation and shed light on physical characteristics of these different types of DNA topoisomers during electrophoresis.

A genomic island present along the bacterial chromosome of the Parachlamydiaceae UWE25, an obligate amoebal endosymbiont, encodes a potentially functional F-like conjugative DNA transfer system.

BACKGROUND: The genome of Protochlamydia amoebophila UWE25, a Parachlamydia-related endosymbiont of free-living amoebae, was recently published, providing the opportunity to search for genomic islands (GIs). RESULTS: On the residual cumulative G+C content curve, a G+C-rich 19-kb region was observed. This sequence is part of a 100-kb chromosome region, containing 100 highly co-oriented ORFs, flanked by two 17-bp direct repeats. Two identical gly-tRNA genes in tandem are present at the proximal end of this genetic element. Several mobility genes encoding transposases and bacteriophage-related proteins are located within this chromosome region. Thus, this region largely fulfills the criteria of GIs. The G+C content analysis shows that several modules compose this GI. Surprisingly, one of them encodes all genes essential for F-like conjugative DNA transfer (traF, traG, traH, traN, traU, traW, and trbC), involved in sex pilus retraction and mating pair stabilization, strongly suggesting that, similarly to the other F-like operons, the parachlamydial tra unit is devoted to DNA transfer. A close relatedness of this tra unit to F-like tra operons involved in conjugative transfer is confirmed by phylogenetic analyses performed on concatenated genes and gene order conservation. These analyses and that of gly-tRNA distribution in 140 GIs suggest a proteobacterial origin of the parachlamydial tra unit. CONCLUSIONS: A GI of the UWE25 chromosome encodes a potentially functional F-like DNA conjugative system. This is the first hint of a putative conjugative system in chlamydiae. Conjugation most probably occurs within free-living amoebae, that may contain hundreds of Parachlamydia bacteria tightly packed in vacuoles. Such a conjugative system might be involved in DNA transfer between internalized bacteria. Since this system is absent from the sequenced genomes of Chlamydiaceae, we hypothesize that it was acquired after the divergence between Parachlamydiaceae and Chlamydiaceae, when the Parachlamydia-related symbiont was an intracellular bacteria. It suggests that this heterologous DNA was acquired from a phylogenetically-distant bacteria sharing an amoebal vacuole. Since Parachlamydiaceae are emerging agents of pneumonia, this GI might be involved in pathogenicity. In future, conjugative systems might be developed as genetic tools for Chlamydiales.

The DNA-binding domain (DBD) is essential for NR2E3 dimerization and interaction with Crx

Purpose:NR2E3 (PNR) is an orphan nuclear receptor essential for proper photoreceptor determination and differentiation. In humans, mutations in NR2E3 have been associated with the recessively inherited enhanced short wavelength sensitive (S-) cone syndrome (ESCS) and, more recently, with autosomal dominant retinitis pigmentosa (adRP). NR2E3 acts in concert with the transcription factors Crx and Nrl to repress cone-specific genes and activate rod-specific genes. NR2E3 and Crx have been shown to physically interact by their DNA-binding domain (DBD), which may also be implicated in the dimerization process of the nuclear receptor. However, neither NR2E3 homodimerization nor NR2E3/Crx complex formation has been investigated in detail. Methods:In this present work, we analyzed the dimerization of the NR2E3 protein and its interaction with Crx by bioluminescence resonance energy transfer (BRET2) which utilizes Renilla luciferase (hRluc) protein and its substrate DeepBlueC as an energy donor and a mutant green fluorescent protein (GFP2) as the acceptor. We investigated, on whole intact cells, the role of NR2E3 DBD-mutations in dimerization and association with Crx. Results:We clearly showed that NR2E3 formed homodimers in HEK-293T cells. Moreover, all causative NR2E3 mutations present in the DBD of the protein showed an alteration in dimerization, except for the R76Q and the R104W mutants. Interestingly, the adRP-linked G56R mutant was the only DBD-NR2E3 mutant that showed a correct interaction with Crx. Finally, we observed a decrease in rhodospin gene transactivation for all DBD-NR2E3 mutants tested and no potentiation for the adRP-linked G56R mutant. In addition, the p.G56R mutant enhanced the transrepression of M-opsin promoter, while all other DBD-NR2E3 mutants did not repress M-opsin transactivation. Conclusions:A defect, either in the dimer formation or in the interaction of NR2E3 with Crx, leads to abnormal transcriptional activity on rhodopsin and M-opsin promoter and to an atypical retinal development; while the titration of Crx by p.G56R-NR2E3 leads to low levels of rhodopsin and M-opsin expression and may be responsible for the strong adRP phenotype.

Satisfaction survey with DNA cards method to collect genetic samples for pharmacogenetics studies

Background: Pharmacogenetic studies are essential in understanding the interindividual variability of drug responses. DNA sample collection for genotyping is a critical step in genetic studies. A method using dried blood samples from finger-puncture, collected on DNA-cards, has been described as an alternative to the usual venepuncture technique. The purpose of this study is to evaluate the implementation of the DNA cards method in a multicentre clinical trial, and to assess the degree of investigators' satisfaction and the acceptance of the patients perceived by the investigators.Methods: Blood samples were collected on DNA-cards. The quality and quantity of DNA recovered were analyzed. Investigators were questioned regarding their general interest, previous experience, safety issues, preferences and perceived patient satisfaction. Results: 151 patients' blood samples were collected. Genotyping of GST polymorphisms was achieved in all samples (100%). 28 investigators completed the survey. Investigators perceived patient satisfaction as very good (60.7%) or good (39.3%), without reluctance to finger puncture. Investigators preferred this method, which was considered safer and better than the usual methods. All investigators would recommend using it in future genetic studies. Conclusion: Within the clinical trial setting, the DNA-cards method was very well accepted by investigators and patients (in perception of investigators), and was preferred to conventional methods due to its ease of use and safety.

Identification of Methylated Genes Associated with Aggressive Clinicopathological Features in Mantle Cell Lymphoma

Background: Mantle cell lymphoma (MCL) is genetically characterized by the t(11;14)(q13;q32) translocation and a high number of secondary chromosomal alterations. The contribution of DNA methylation to MCL lymphomagenesis is not well known. We sought to identify epigenetically silenced genes in these tumours that might have clinical relevance. Methodology/Principal Findings: To identify potential methylated genes in MCL we initially investigated seven MCL cell lines treated with epigenetic drugs and gene expression microarray profiling. The methylation status of selected candidate genes was validated by a quantitative assay and subsequently analyzed in a series of primary MCL (n=38). After pharmacological reversion we identified 252 potentially methylated genes. The methylation analysis of a subset of these genes (n=25) in the MCL cell lines and normal B lymphocytes confirmed that 80% of them were methylated in the cell lines but not in normal lymphocytes. The subsequent analysis in primary MCL identified five genes (SOX9,HOXA9,AHR,NR2F2 ,and ROBO1) frequently methylated in these tumours. The gene methylation events tended to occur in the same primary neoplasms and correlated with higher proliferation, increased number of chromosomal abnormalities, and shorter survival of the patients. Conclusions: We have identified a set of genes whose methylation degree and gene expression levels correlate with aggressive clinicopathological features of MCL. Our findings also suggest that a subset of MCL might show a CpG island methylator phenotype (CIMP) that may influence the behaviour of the tumours.