Random-set methods identify distinct aspects of the enrichment signal in gene-set analysis

A prespecified set of genes may be enriched, to varying degrees, for genes that have altered expression levels relative to two or more states of a cell. Knowing the enrichment of gene sets defined by functional categories, such as gene ontology (GO) annotations, is valuable for analyzing the biological signals in microarray expression data. A common approach to measuring enrichment is by cross-classifying genes according to membership in a functional category and membership on a selected list of significantly altered genes. A small Fisher's exact test $p$-value, for example, in this $2\times2$ table is indicative of enrichment. Other category analysis methods retain the quantitative gene-level scores and measure significance by referring a category-level statistic to a permutation distribution associated with the original differential expression problem. We describe a class of random-set scoring methods that measure distinct components of the enrichment signal. The class includes Fisher's test based on selected genes and also tests that average gene-level evidence across the category. Averaging and selection methods are compared empirically using Affymetrix data on expression in nasopharyngeal cancer tissue, and theoretically using a location model of differential expression. We find that each method has a domain of superiority in the state space of enrichment problems, and that both methods have benefits in practice. Our analysis also addresses two problems related to multiple-category inference, namely, that equally enriched categories are not detected with equal probability if they are of different sizes, and also that there is dependence among category statistics owing to shared genes. Random-set enrichment calculations do not require Monte Carlo for implementation. They are made available in the R package allez.Comment: Published at http://dx.doi.org/10.1214/07-AOAS104 in the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org

Equine Arteritis Virus Subgenomic RNA Transcription: UV Inactivation and Translation Inhibition Studies

AbstractThe expression of the genetic information of equine arteritis virus (EAV), an arterivirus, involves the synthesis of six subgenomic (sg) mRNAs. These are 5′ and 3′ coterminal since they are composed of a leader and a body sequence, which are identical to the 5′ and 3′ ends of the genome, respectively. Previously, it has been suggested thatcis-splicing of a genome-length precursor RNA is involved in their synthesis. This was reevaluated in a comparative analysis of the sg RNA synthesis of EAV, the coronavirus mouse hepatitis virus (MHV), and the alphavirus Sindbis virus. UV transcription mapping showed that the majority of the EAV sg RNAs made at later stages of infection is not derived from a genome-length precursor. However, complete independence of sg RNA synthesis from that of genomic RNA was never observed during the course of infection. The possibility that this resulted from UV irradiation-induced effects on the synthesis of the viral replicase was investigated by inhibiting translation using cycloheximide. For EAV, ongoing protein synthesis was found to be more important for the synthesis of sg RNA than for that of genomic RNA. In general, MHV transcription was extremely sensitive to translation inhibition, whereas EAV genomic RNA synthesis became independent ofde novoprotein synthesis late in infection

Statistical Use of Argonaute Expression and RISC Assembly in microRNA Target Identification

MicroRNAs (miRNAs) posttranscriptionally regulate targeted messenger RNAs (mRNAs) by inducing cleavage or otherwise repressing their translation. We address the problem of detecting m/miRNA targeting relationships in homo sapiens from microarray data by developing statistical models that are motivated by the biological mechanisms used by miRNAs. The focus of our modeling is the construction, activity, and mediation of RNA-induced silencing complexes (RISCs) competent for targeted mRNA cleavage. We demonstrate that regression models accommodating RISC abundance and controlling for other mediating factors fit the expression profiles of known target pairs substantially better than models based on m/miRNA expressions alone, and lead to verifications of computational target pair predictions that are more sensitive than those based on marginal expression levels. Because our models are fully independent of exogenous results from sequence-based computational methods, they are appropriate for use as either a primary or secondary source of information regarding m/miRNA target pair relationships, especially in conjunction with high-throughput expression studies.</p

Intersection of the Multivesicular Body Pathway and Lipid Homeostasis in RNA Replication by a Positive-Strand RNA Virus

Like many positive-strand RNA viruses, brome mosaic virus (BMV) RNA replication occurs in membrane-invaginated vesicular compartments. BMV RNA replication compartments show parallels with membrane-enveloped, budding retrovirus virions, whose release depends on the cellular multivesicular body (MVB) sorting pathway. BMV RNA replication compartments are not released from their parent membranes, but might depend on MVB functions for membrane invagination. Prior results show that BMV RNA replication is severely inhibited by deletion of the crucial MVB gene DOA4 or BRO1. We report here that involvement of DOA4 and BRO1 in BMV RNA replication is not dependent on the MVB pathway's membrane-shaping functions but rather is due to their roles in recycling ubiquitin from MVB cargos. We show that deleting DOA4 or BRO1 inhibits the ubiquitination- and proteasome-dependent activation of homologous transcription factors Mga2p and Spt23p, which regulate many lipid metabolism genes, including the fatty acid desaturase gene OLE1, which is essential for BMV RNA replication. However, Mga2p processing and BMV RNA replication are restored by supplementing free ubiquitin, which is depleted in doa4Δ and bro1Δ cells. The results identify Mga2p and Spt23p processing and lipid regulation as sensitive targets of ubiquitin depletion and correctly predict multiple effects of modulating additional host genes RFU1, UBP6, and UFD3. Our results also show that BMV RNA replication depends on additional Mga2p-regulated genes likely involved in lipid metabolism beyond OLE1. Among other points, these findings show the potential for blocking viral RNA replication by modulating lipid synthesis at multiple levels

Identification of an alternative triglyceride biosynthesis pathway

Triacylglycerols (TAGs) are the main source of stored energy in the body, providing an important substrate pool for mitochondrial beta-oxidation. Imbalances in the amount of TAGs are associated with obesity, cardiac disease and various other pathologies 1,2. In humans, TAGs are synthesized from excess, coenzyme A-conjugated fatty acids by diacylglycerol O-acyltransferases (DGAT1 and DGAT2) 3. In other organisms, this activity is complemented by additional enzymes 4, but whether such alternative pathways exist in humans remains unknown. Here we disrupt the DGAT pathway in haploid human cells and use iterative genetics to reveal an unrelated TAG-synthesizing system composed of a protein we called DIESL (also known as TMEM68, an acyltransferase of previously unknown function) and its regulator TMX1. Mechanistically, TMX1 binds to and controls DIESL at the endoplasmic reticulum, and loss of TMX1 leads to the unconstrained formation of DIESL-dependent lipid droplets. DIESL is an autonomous TAG synthase, and expression of human DIESL in Escherichia coli endows this organism with the ability to synthesize TAG. Although both DIESL and the DGATs function as diacylglycerol acyltransferases, they contribute to the cellular TAG pool under specific conditions. Functionally, DIESL synthesizes TAG at the expense of membrane phospholipids and maintains mitochondrial function during periods of extracellular lipid starvation. In mice, DIESL deficiency impedes rapid postnatal growth and affects energy homeostasis during changes in nutrient availability. We have therefore identified an alternative TAG biosynthetic pathway driven by DIESL under potent control by TMX1. </p

Statistical Use of Argonaute Expression and RISC Assembly in microRNA Target Identification

MicroRNAs (miRNAs) posttranscriptionally regulate targeted messenger RNAs (mRNAs) by inducing cleavage or otherwise repressing their translation. We address the problem of detecting m/miRNA targeting relationships in homo sapiens from microarray data by developing statistical models that are motivated by the biological mechanisms used by miRNAs. The focus of our modeling is the construction, activity, and mediation of RNA-induced silencing complexes (RISCs) competent for targeted mRNA cleavage. We demonstrate that regression models accommodating RISC abundance and controlling for other mediating factors fit the expression profiles of known target pairs substantially better than models based on m/miRNA expressions alone, and lead to verifications of computational target pair predictions that are more sensitive than those based on marginal expression levels. Because our models are fully independent of exogenous results from sequence-based computational methods, they are appropriate for use as either a primary or secondary source of information regarding m/miRNA target pair relationships, especially in conjunction with high-throughput expression studies

Psychosocial impact of undergoing prostate cancer screening for men with BRCA1 or BRCA2 mutations.

OBJECTIVES: To report the baseline results of a longitudinal psychosocial study that forms part of the IMPACT study, a multi-national investigation of targeted prostate cancer (PCa) screening among men with a known pathogenic germline mutation in the BRCA1 or BRCA2 genes. PARTICPANTS AND METHODS: Men enrolled in the IMPACT study were invited to complete a questionnaire at collaborating sites prior to each annual screening visit. The questionnaire included sociodemographic characteristics and the following measures: the Hospital Anxiety and Depression Scale (HADS), Impact of Event Scale (IES), 36-item short-form health survey (SF-36), Memorial Anxiety Scale for Prostate Cancer, Cancer Worry Scale-Revised, risk perception and knowledge. The results of the baseline questionnaire are presented. RESULTS: A total of 432 men completed questionnaires: 98 and 160 had mutations in BRCA1 and BRCA2 genes, respectively, and 174 were controls (familial mutation negative). Participants' perception of PCa risk was influenced by genetic status. Knowledge levels were high and unrelated to genetic status. Mean scores for the HADS and SF-36 were within reported general population norms and mean IES scores were within normal range. IES mean intrusion and avoidance scores were significantly higher in BRCA1/BRCA2 carriers than in controls and were higher in men with increased PCa risk perception. At the multivariate level, risk perception contributed more significantly to variance in IES scores than genetic status. CONCLUSION: This is the first study to report the psychosocial profile of men with BRCA1/BRCA2 mutations undergoing PCa screening. No clinically concerning levels of general or cancer-specific distress or poor quality of life were detected in the cohort as a whole. A small subset of participants reported higher levels of distress, suggesting the need for healthcare professionals offering PCa screening to identify these risk factors and offer additional information and support to men seeking PCa screening

Large expert-curated database for benchmarking document similarity detection in biomedical literature search

Document recommendation systems for locating relevant literature have mostly relied on methods developed a decade ago. This is largely due to the lack of a large offline gold-standard benchmark of relevant documents that cover a variety of research fields such that newly developed literature search techniques can be compared, improved and translated into practice. To overcome this bottleneck, we have established the RElevant LIterature SearcH consortium consisting of more than 1500 scientists from 84 countries, who have collectively annotated the relevance of over 180 000 PubMed-listed articles with regard to their respective seed (input) article/s. The majority of annotations were contributed by highly experienced, original authors of the seed articles. The collected data cover 76% of all unique PubMed Medical Subject Headings descriptors. No systematic biases were observed across different experience levels, research fields or time spent on annotations. More importantly, annotations of the same document pairs contributed by different scientists were highly concordant. We further show that the three representative baseline methods used to generate recommended articles for evaluation (Okapi Best Matching 25, Term Frequency-Inverse Document Frequency and PubMed Related Articles) had similar overall performances. Additionally, we found that these methods each tend to produce distinct collections of recommended articles, suggesting that a hybrid method may be required to completely capture all relevant articles. The established database server located at https://relishdb.ict.griffith.edu.au is freely available for the downloading of annotation data and the blind testing of new methods. We expect that this benchmark will be useful for stimulating the development of new powerful techniques for title and title/abstract-based search engines for relevant articles in biomedical research.Peer reviewe

Multifunctional Protein A Is the Only Viral Protein Required for Nodavirus RNA Replication Crown Formation

Positive-strand RNA virus RNA genome replication occurs in membrane-associated RNA replication complexes (RCs). Nodavirus RCs are outer mitochondrial membrane invaginations whose necked openings to the cytosol are &ldquo;crowned&rdquo; by a 12-fold symmetrical proteinaceous ring that functions as the main engine of RNA replication. Similar protein crowns recently visualized at the openings of alphavirus and coronavirus RCs highlight their broad conservation and functional importance. Using cryo-EM tomography, we earlier showed that the major nodavirus crown constituent is viral protein A, whose polymerase, RNA capping, membrane interaction and multimerization domains drive RC formation and function. Other viral proteins are strong candidates for unassigned EM density in the crown. RNA-binding RNAi inhibitor protein B2 co-immunoprecipitates with protein A and could form crown subdomains that protect nascent viral RNA and dsRNA templates. Capsid protein may interact with the crown since nodavirus virion assembly has spatial and other links to RNA replication. Using cryoelectron tomography and complementary approaches, we show that, even when formed in mammalian cells, nodavirus RC crowns generated without B2 and capsid proteins are functional and structurally indistinguishable from mature crowns in infected Drosophila cells expressing all viral proteins. Thus, the only nodaviral factors essential to form functional RCs and crowns are RNA replication protein A and an RNA template. We also resolve apparent conflicts in prior results on B2 localization in infected cells, revealing at least two distinguishable pools of B2. The results have significant implications for crown structure, assembly, function and control as an antiviral target