Search PubMed⌕ Search

Biomedical subjects

Ron Unger

Publications and source records attributed to Ron Unger.

13 recordsLinked to original sources

Genome-wide association, polygenic risk scores, and machine learning for chronic post-surgical pain risk stratification: A UK biobank study.

Chronic post-surgical pain is a prevalent and debilitating complication following surgery, representing a clinical challenge. Despite the established heritability of pain phenotypes, large-scale genetic studies remain limited. This study aimed to identify genetic variants associated with chronic post-surgical pain, develop polygenic risk scores, and integrate these with clinical features for risk prediction. UK Biobank data from 47,836 participants (2490 cases and 45,346 controls) were split into training (80%; n = 38,268) and validation (20%; n = 9568) sets prior to analysis. A genome-wide association study was conducted on the training set only, across 19 million variants, and polygenic risk scores were constructed and integrated with clinical features in a logistic regression framework. Two close, rare, imputed signals crossed the genome-wide significance threshold but lacked local linkage-disequilibrium support, while 220 variants crossed the suggestive threshold. In the held-out validation set, cases had higher mean polygenic risk scores than controls (0.138 vs. -0.021; Cohen's d = 0.16, p < 0.001). A logistic regression model integrating clinical features and polygenic risk scores achieved an area under the curve of 0.639 (95% CI: 0.583-0.693), higher than models using either feature set alone. The polygenic risk score for chronic post-surgical pain was among the most important predictors. Risk stratification revealed the top quartile had 3.84-fold higher odds of chronic post-surgical pain than the bottom quartile (95% CI: 2.00-7.37). These findings suggest a possible modest genetic contribution to chronic post-surgical pain. Polygenic risk scores may complement clinical factors in surgical risk stratification. PERSPECTIVE: Chronic post-surgical pain may have a modest genetic contribution. This UK Biobank study identified over 220 variants at suggestive significance and constructed a polygenic risk score that was significantly elevated in cases. A combined clinical-genomic model achieved a 3.84-fold difference in odds across predicted-risk quartiles.

Chronic post-surgical pain↗

Genome-wide analysis of C/D and H/ACA-like small nucleolar RNAs in Leishmania major indicates conservation among trypanosomatids in the repertoire and in their rRNA targets.

Small nucleolar RNAs (snoRNAs) are a large group of noncoding RNAs that exist in eukaryotes and archaea and guide modifications such as 2'-O-ribose methylations and pseudouridylation on rRNAs and snRNAs. Recently, we described a genome-wide screening approach with Trypanosoma brucei that revealed over 90 guide RNAs. In this study, we extended this approach to analyze the repertoire of the closely related human pathogen Leishmania major. We describe 23 clusters that encode 62 C/Ds that can potentially guide 79 methylations and 37 H/ACA-like RNAs that can potentially guide 30 pseudouridylation reactions. Like T. brucei, Leishmania also contains many modifications and guide RNAs relative to its genome size. This study describes 10 H/ACAs and 14 C/Ds that were not found in T. brucei. Mapping of 2'-O-methylations in rRNA regions rich in modifications suggests the existence of trypanosomatid-specific modifications conserved in T. brucei and Leishmania. Structural features of C/D snoRNAs, such as copy number, conservation of boxes, K turns, and intragenic and extragenic base pairing, were examined to elucidate the great variation in snoRNA abundance. This study highlights the power of comparative genomics for determining conserved features of noncoding RNAs.

Animals↗

Towards computing with proteins.

Can proteins be used as computational devices to address difficult computational problems? In recent years there has been much interest in biological computing, that is, building a general purpose computer from biological molecules. Most of the current efforts are based on DNA because of its ability to self-hybridize. The exquisite selectivity and specificity of complex protein-based networks motivated us to suggest that similar principles can be used to devise biological systems that will be able to directly implement any logical circuit as a parallel asynchronous computation. Such devices, powered by ATP molecules, would be able to perform, for medical applications, digital computation with natural interface to biological input conditions. We discuss how to design protein molecules that would serve as the basic computational element by functioning as a NAND logical gate, utilizing DNA tags for recognition, and phosphorylation and exonuclease reactions for information processing. A solution of these elements could carry out effective computation. Finally, the model and its robustness to errors were tested in a computer simulation.

Adenosine Triphosphate↗

Assessing societal impacts when planning restoration of large alluvial rivers: a case study of the Sacramento River project, California.

Studies have shown that ecological restoration projects are more likely to gain public support if they simultaneously increase important human services that natural resources provide to people. River restoration projects have the potential to influence many of the societal functions (e.g., flood control, water quality) that rivers provide, yet most projects fail to consider this in a comprehensive manner. Most river restoration projects also fail to take into account opportunities for revitalization of large-scale river processes, focusing instead on opportunities presented at individual parcels. In an effort to avoid these pitfalls while planning restoration of the Sacramento River, we conducted a set of coordinated studies to evaluate societal impacts of alternative restoration actions over a large geographic area. Our studies were designed to identify restoration actions that offer benefits to both society and the ecosystem and to meet the information needs of agency planning teams focusing on the area. We worked with local partners and public stakeholders to design and implement studies that assessed the effects of alternative restoration actions on flooding and erosion patterns, socioeconomics, cultural resources, and public access and recreation. We found that by explicitly and scientifically melding societal and ecosystem perspectives, it was possible to identify restoration actions that simultaneously improve both ecosystem health and the services (e.g., flood protection and recreation) that the Sacramento River and its floodplain provide to people. Further, we found that by directly engaging with local stakeholders to formulate, implement, and interpret the studies, we were able to develop a high level of trust that ultimately translated into widespread support for the project.

Animals↗

A genome-wide analysis of C/D and H/ACA-like small nucleolar RNAs in Trypanosoma brucei reveals a trypanosome-specific pattern of rRNA modification.

Small nucleolar RNAs (snoRNAs) constitute newly discovered noncoding small RNAs, most of which function in guiding modifications such as 2'-O-ribose methylation and pseudouridylation on rRNAs and snRNAs. To investigate the genome organization of Trypanosoma brucei snoRNAs and the pattern of rRNA modifications, we used a whole-genome approach to identify the repertoire of these guide RNAs. Twenty-one clusters encoding for 57 C/D snoRNAs and 34 H/ACA-like RNAs, which have the potential to direct 84 methylations and 32 pseudouridines, respectively, were identified. The number of 2'-O-methyls (Nms) identified on rRNA represent 80% of the expected modifications. The modifications guided by these RNAs suggest that trypanosomes contain many modifications and guide RNAs relative to their genome size. Interestingly, approximately 40% of the Nms are species-specific modifications that do not exist in yeast, humans, or plants, and 40% of the species-specific predicted modifications are located in unique positions outside the highly conserved domains. Although most of the guide RNAs were found in reiterated clusters, a few single-copy genes were identified. The large repertoire of modifications and guide RNAs in trypanosomes suggests that these modifications possibly play a central role in these parasites.

Animals↗

Small nucleolar RNAs that guide modification in trypanosomatids: repertoire, targets, genome organisation, and unique functions.

Small nucleolar RNAs constitute a family of newly discovered non-coding small RNAs, most of which function in guiding RNA modifications. Two prevalent types of modifications are 2'-O-methylation and pseudouridylation. The modification is directed by the formation of a canonical small nucleolar RNA-target duplex. Initially, RNA-guided modification was shown to take place on rRNA, but recent studies suggest that small nuclear RNA, mRNA, tRNA, and the trypanosome spliced leader RNA also undergo guided modifications. Trypanosomes contain more modifications and potentially more small nucleolar RNAs than yeast, and the increased number of modifications may help to preserve ribosome function under adverse environmental conditions during the cycling between the insect and mammalian host. The genome organisation in clusters carrying the two types of small nucleolar RNAs, C/D and H/ACA-like RNAs, resembles that in plants. However, the trypanosomatid H/ACA RNAs are similar to those found in Archaea and are composed of a single hairpin that may represent the primordial H/ACA RNA. In this review we summarise this new field of trypanosome small nucleolar RNAs, emphasising the open questions regarding the number of small nucleolar RNAs, the repertoire, genome organisation, and the unique function of guided modifications in these protozoan parasites.

Animals↗

Identification and functional characterization of lsm proteins in Trypanosoma brucei.

RNA interference of Sm proteins in Trypanosoma brucei demonstrated that the stability of the small nuclear RNAs (U1, U2, U4, U5) and the spliced leader RNA, but not U6 RNA, were affected upon Sm depletion (Mandelboim, M., Barth, S., Biton, M., Liang, X. H., and Michaeli, S. (2003) J. Biol. Chem. 278, 51469-51478), suggesting that Lsm proteins that bind and stabilize U6 RNA in other eukaryotes should exist in trypanosomes. In this study, we identified seven Lsm proteins (Lsm2p to Lsm8p) and examined the function of Lsm3p and Lsm8p by RNA interference silencing. Both Lsm proteins were found to be essential for U6 stability and mRNA decay. Silencing was lethal, and cis- and trans-splicing were inhibited. Importantly, silencing also affected the level of U4.U6 and the U4.U6/U5 tri-small nuclear ribonucleoprotein complexes. The presence of Lsm proteins in trypanosomes that diverged early in the eukaryotic lineage suggests that these proteins are highly conserved in both structure and function among eukaryotes. Interestingly, however, Lsm1p that is specific to the mRNA decay complex was not identified in the genome data base of any kinetoplastidae, and the Lsm8p that in other eukaryotes exclusively functions in U6 stability was found to function in trypanosomes also in mRNA decay. These data therefore suggest that in trypanosomes only a single Lsm complex may exist.

Amino Acid Sequence↗

Scaling law in sizes of protein sequence families: from super-families to orphan genes.

It has been observed that the size of protein sequence families is unevenly distributed, with few super families with a large number of members and many "orphan" proteins that do not belong to any family. Here it is shown that the distribution of sizes of protein families in different databases and classifications (Protomap, Prodom, Cog) follows a power-law behavior with similar scaling exponents, which is characteristic of self-organizing systems. Since large databases are used in this study, a more detailed analysis of the data than in previous studies was possible. Hence, it is shown that the size distribution is governed by two exponents, different for the super families and the orphan proteins. A simple model of protein evolution is proposed, in which proteins are dynamically generated and clustered into families. The model yields a scaling behavior very similar to the distribution observed in the actual sequence databases, including the two distinct regimes for the large and small families, and thus suggests that the existence of "super families" of proteins and "orphan" proteins are two manifestations of the same evolutionary process.

Algorithms↗

A rapid method for detection of putative RNAi target genes in genomic data.

RNAi, inhibition of gene expression by double stranded RNA molecules, has rapidly become a powerful laboratory technique to study gene function. The effectiveness of the procedure raised the question of whether this laboratory technique may actually mimic a natural cellular control mechanism that works on similar principles. Indeed recent evidence is accumulating to suggest that RNAi is a natural control mechanism that might even serve as a primitive immune response against RNA viruses and retroposons. Three different interference scenarios seem to be utilized by various RNAi mechanisms. One of the mechanisms involves degradation of mRNA molecules. Here we suggest a method to systematically scan entire genomes simultaneously for RNAi elements and the presence of cellular genes that are degraded by these RNAi elements via exact short base-pair matching. The method is based on scanning the genomes using a suffix tree data structure that was specifically modified to identify sets of combinations of repeated and inverted repeated sequences of 20 bp or more. Initial scan suggest that a large number, about 7% of C.elegans and 3% of C.briggsae genes, have the potential to be subject to natural RNAi control. Two methods are proposed to further analyze these genes to select the cases that are more likely to be actual cases of RNAi control. One method involves looking for ESTs that can provide direct evidence that RNAi control element are indeed expressed. The other method looks for synteny between C.elegans and C.briggsae assuming that genes that might be under RNAi control in both organisms are more likely to be biological significant. Taken together, supportive evidence was found for about 70 genes to be under RNAi control. Among these genes are: transposase, hormone receptors, homeobox proteins, defensin, actins, and several types of collagens. While our method is not capable of detecting all cases of natural RNAi control, it points to a large number of potential cases that can be further verified by experimental work.

Animals↗

RNA interference of signal peptide-binding protein SRP54 elicits deleterious effects and protein sorting defects in trypanosomes.

Trypanosomes are protozoan parasites that have a major impact on health. This family diverged very early from the eukaryotic lineage and possesses unique RNA processing mechanisms such as trans-splicing and RNA editing. The trypanosome signal recognition particle (SRP) has a unique composition compared with all known SRP complexes, because it contains two RNA molecules, the 7SL RNA and a tRNA-like molecule. RNA interference was utilized to elucidate the essentiality of the SRP pathway and its role in protein translocation in Trypanosoma brucei. The production of double stranded RNA specific for the signal peptide-binding protein SRP54 induced the degradation of the mRNA and a loss of the SRP54 protein. SRP54 depletion elicited inhibition in growth and cytokinesis, suggesting that the SRP pathway is essential. The translocation of four signal peptide-containing proteins was examined. Surprisingly, the proteins were translocated to the endoplasmic reticulum and properly processed. However, the surface EP procyclin, the lysosomal protein p67, and the flagellar pocket protein CRAM were mislocalized and accumulated in megavesicles, most likely because of a secondary effect on protein sorting. The translocation of these proteins to the endoplasmic reticulum under SRP54 depletion suggests that an alternative pathway for protein translocation exists in trypanosomes.

Amino Acid Sequence↗

Swaps in protein sequences.

An important question in protein evolution is to what extent proteins may have undergone swaps (switches of domain or fragment order) during evolution. Such events might have occurred in several forms: Swaps of short fragments, swaps of structural and functional motifs, or recombination of domains in multidomain proteins. This question is important for the theoretical understanding of the evolution of proteins, and has practical implications for using swaps as a design tool in protein engineering. In order to analyze the question systematically, we conducted a large scale survey of possible swaps and permutations among all pairs of protein from the Swissport database. A swap is defined as a specific kind of sequence mutation between two proteins in which two fragments that appear in both sequences have different relative order in the two sequences. For example, aXbYc and dYeXf are defined as a swap, where X and Y represent sequence fragments that switched their order. Identifying such swaps is difficult using standard sequence comparison packages. One of the main problems in the analysis stems from the fact that many sequences contain repeats, which may be identified as false-positive swaps. We have used two different approaches to detect pairs of proteins with swaps. The first approach is based on the predefined list of domains in Pfam. We identified all the proteins that share at least two domains and analyzed their relative order, looking for pairs in which the order of these domains was switched. We designed an algorithm to distinguish between real swaps and duplications. In the second approach, we used Blast to detect pairs of proteins that share several fragments. Then, we used an automatic procedure to select pairs that are likely to contain swaps. Those pairs were analyzed visually, using a graphical tool, to eliminate duplications. Combining these approaches, about 140 different cases of swaps in the Swissprot database were found (after eliminating multiple pairs within the same family). Some of the cases have been described in the literature, but many are novel examples. Although each new example identified may be interesting to analyze, our main conclusion is that cases of swaps are rare in protein evolution. This observation is at odds with the common view that proteins are very modular to the point that modules (e.g., domains) can be shuffled between proteins with minimal constraints. Our study suggests that sequential constraints, i.e., the relative order between domains, are highly conserved.

Algorithms↗

PESI--an intelligent system for prediction of enzyme-substrate interactions based on experimental constraints.

We present a system for predicting protein-protein modifications, and demonstrate its usefulness in the field of signal transduction research. Signal transduction is one of the most important areas of investigation in biological research. One of the major mechanisms frequently employed by cells to regulate signal transduction processes involves protein phosphorylation by various kinases. As many as 1,000 protein kinases and 500 protein phosphatases in the human genome are thought to be involved in phosphorylation processes which regulate all aspects of cell function. The complexity of such interactions stems from the enormous number of factors and interactions, which makes the identification of putative substrates for any given enzyme by straightforward experimentation increasingly difficult. We present here a data mining algorithm, based on the similarity between the modifier proteins and between the modified proteins, and on experimental constraints. The application presented here (PESI) focuses on substrate phosphorylation by various enzymes. This algorithm reduces the number of substrate candidates for experimental study by about two orders of magnitude. Moreover, this algorithm has already yielded predictions for previously unknown substrates of the enzymes PKCdelta and PKCeta, which we have confirmed experimentally.

Algorithms↗

Designing an A* algorithm for calculating edit distance between rooted-unordered trees.

Tree structures are useful for describing and analyzing biological objects and processes. Consequently, there is a need to design metrics and algorithms to compare trees. A natural comparison metric is the "Tree Edit Distance," the number of simple edit (insert/delete) operations needed to transform one tree into the other. Rooted-ordered trees, where the order between the siblings is significant, can be compared in polynomial time. Rooted-unordered trees are used to describe processes or objects where the topology, rather than the order or the identity of each node, is important. For example, in immunology, rooted-unordered trees describe the process of immunoglobulin (antibody) gene diversification in the germinal center over time. Comparing such trees has been proven to be a difficult computational problem that belongs to the set of NP-Complete problems. Comparing two trees can be viewed as a search problem in graphs. A* is a search algorithm that explores the search space in an efficient order. Using a good lower bound estimation of the degree of difference between the two trees, A* can reduce search time dramatically. We have designed and implemented a variant of the A* search algorithm suitable for calculating tree edit distance. We show here that A* is able to perform an edit distance measurement in reasonable time for trees with dozens of nodes.

Algorithms↗