Postmortem toxicology findings from the Camden Opioid Research Initiative.

The United States continues to be impacted by decades of an opioid misuse epidemic, worsened by the COVID-19 pandemic and by the growing prevalence of highly potent synthetic opioids (HPSO) such as fentanyl. In instances of a toxicity event, first-response administration of reversal medications such as naloxone can be insufficient to fully counteract the effects of HPSO, particularly when there is co-occurring substance use. In an effort to characterize and study this multi-faceted problem, the Camden Opioid Research Initiative (CORI) has been formed. The CORI study has collected and analyzed post-mortem toxicology data from 42 cases of decedents who expired from opioid-related toxicity in the South New Jersey region to characterize substance use profiles. Co-occurring substance use, whether by intent or through possible contamination of the illicit opioid supply, is pervasive among deaths due to opioid toxicity, and evidence of medication-assisted treatment is scarce. Nearly all (98%) of the toxicology cases show the presence of the HPSO, fentanyl, and very few (7%) results detected evidence of medication-assisted treatment for opioid use disorder, such as buprenorphine or methadone, at the time of death. The opioid toxicity reversal drug, naloxone, was detected in 19% of cases, but 100% of cases expressed one or more stimulants, and sedatives including xylazine were detected in 48% of cases. These results showing complex substance use profiles indicate that efforts at mitigating the opioid misuse epidemic must address the complications presented by co-occurring stimulant and other substance use, and reduce barriers to and stigmas of seeking effective medication-assisted treatments.

Patient Perceptions and Potential Utility of Pharmacogenetic Testing in Chronic Pain Management and Opioid Use Disorder in the Camden Opioid Research Initiative.

Pharmacogenetics (PGx) has the potential to improve opioid medication management. Here, we present patient perception data, pharmacogenetic data and medication management trends in patients with chronic pain (arm 1) and opioid use disorder (arm 2) treated at Cooper University Health Care in Camden City, NJ. Our results demonstrate that the majority of patients in both arms of the study (55% and 65%, respectively) are open to pharmacogenetic testing, and most (66% and 69%, respectively) believe that genetic testing has the potential to improve their medical care. Our results further support the potential for CYP2D6 PGx testing to inform chronic pain medication management for poor metabolizers (PMs) and ultrarapid metabolizers (UMs). Future efforts to implement PGx testing in chronic pain management, however, must address patient concerns about genetic test result access and genetic discrimination.

The Genomics of Opioid Addiction Longitudinal Study (GOALS): study design for a prospective evaluation of genetic and non-genetic factors for development of and recovery from opioid use disorder.

BACKGROUND: The opioid use disorder and overdose crisis in the United States affects public health as well as social and economic welfare. While several genetic and non-genetic risk factors for opioid use disorder have been identified, many of the genetic associations have not been independently replicated, and it is not well understood how these factors interact. This study is designed to evaluate relationships among these factors prospectively to develop future interventions to help prevent or treat opioid use disorder. METHODS: The Genomics of Opioid Addiction Longitudinal Study (GOALS) is a prospective observational study assessing the interplay of genetic and non-genetic by collecting comprehensive genetic and non-genetic information on 400 participants receiving medication for opioid use disorder. Participants will be assessed at four time points over 1 year. A saliva sample will be collected for large-scale genetic data analyses. Non-genetic assessments include validated surveys measuring addiction severity, depression, anxiety, and adverse childhood experiences, as well as treatment outcomes such as urine toxicology results, visit frequency, and number of pre and post-treatment overdoses extracted from electronic medical records. DISCUSSION: We will use these complex data to investigate the relative contributions of genetic and non-genetic risk factors to opioid use disorder and related treatment outcomes.

Correction to: DNA copy number evolution in Drosophila cell lines.

Following publication of the original article [1], the authors reported the following errors.

DNA copy number evolution in Drosophila cell lines.

BACKGROUND: Structural rearrangements of the genome resulting in genic imbalance due to copy number change are often deleterious at the organismal level, but are common in immortalized cell lines and tumors, where they may be an advantage to cells. In order to explore the biological consequences of copy number changes in the Drosophila genome, we resequenced the genomes of 19 tissue-culture cell lines and generated RNA-Seq profiles. RESULTS: Our work revealed dramatic duplications and deletions in all cell lines. We found three lines of evidence indicating that copy number changes were due to selection during tissue culture. First, we found that copy numbers correlated to maintain stoichiometric balance in protein complexes and biochemical pathways, consistent with the gene balance hypothesis. Second, while most copy number changes were cell line-specific, we identified some copy number changes shared by many of the independent cell lines. These included dramatic recurrence of increased copy number of the PDGF/VEGF receptor, which is also over-expressed in many cancer cells, and of bantam, an anti-apoptosis miRNA. Third, even when copy number changes seemed distinct between lines, there was strong evidence that they supported a common phenotypic outcome. For example, we found that proto-oncogenes were over-represented in one cell line (S2-DRSC), whereas tumor suppressor genes were under-represented in another (Kc167). CONCLUSION: Our study illustrates how genome structure changes may contribute to selection of cell lines in vitro. This has implications for other cell-level natural selection progressions, including tumorigenesis.

Identification of neoblast- and regeneration-specific miRNAs in the planarian Schmidtea mediterranea.

In recent years, the planarian Schmidtea mediterranea has emerged as a tractable model system to study stem cell biology and regeneration. MicroRNAs are small RNA species that control gene expression by modulating translational repression and mRNA stability and have been implicated in the regulation of various cellular processes. Though recent studies have identified several miRNAs in S. mediterranea, their expression in neoblast subpopulations and during regeneration has not been examined. Here, we identify several miRNAs whose expression is enriched in different neoblast subpopulations and in regenerating tissue at different time points in S. mediterranea. Some of these miRNAs were enriched within 3 h post-amputation and may, therefore, play a role in wound healing and/or neoblast migration. Our results also revealed miRNAs, such as sme-miR-2d-3p and the sme-miR-124 family, whose expression is enriched in the cephalic ganglia, are also expressed in the brain primordium during CNS regeneration. These results provide new insight into the potential biological functions of miRNAs in neoblasts and regeneration in planarians.

XACT, a long noncoding transcript coating the active X chromosome in human pluripotent cells.

X-chromosome inactivation (XCI) in mammals relies on XIST, a long noncoding transcript that coats and silences the X chromosome in cis. Here we report the discovery of a long noncoding RNA, XACT, that is expressed from and coats the active X chromosome specifically in human pluripotent cells. In the absence of XIST, XACT is expressed from both X chromosomes in humans but not in mice, suggesting a unique role for XACT in the control of human XCI initiation.

Transcriptome analysis reveals strain-specific and conserved stemness genes in Schmidtea mediterranea.

The planarian Schmidtea mediterranea is a powerful model organism for studying stem cell biology due to its extraordinary regenerative ability mediated by neoblasts, a population of adult somatic stem cells. Elucidation of the S. mediterranea transcriptome and the dynamics of transcript expression will increase our understanding of the gene regulatory programs that regulate stem cell function and differentiation. Here, we have used RNA-Seq to characterize the S. mediterranea transcriptome in sexual and asexual animals and in purified neoblast and differentiated cell populations. Our analysis identified many uncharacterized genes, transcripts, and alternatively spliced isoforms that are differentially expressed in a strain or cell type-specific manner. Transcriptome profiling of purified neoblasts and differentiated cells identified neoblast-enriched transcripts, many of which likely play important roles in regeneration and stem cell function. Strikingly, many of the neoblast-enriched genes are orthologs of genes whose expression is enriched in human embryonic stem cells, suggesting that a core set of genes that regulate stem cell function are conserved across metazoan species.

Small RNA pathways in Schmidtea mediterranea.

Planarians are bilaterally symmetrical fresh water organisms capable of regenerating body parts from small fragments following bodily injury. Planarians possess a specialized population of pluripotent cells called neoblasts, which are responsible for their unique regenerative ability. The study of planarian stem cell biology and regeneration has traditionally focused on the transcription factors and proteins that regulate signal transduction pathways. New evidence shows that small RNA molecules are important players in stem cell function and regeneration, yet little is known about the exact nature of their regulatory roles during the regenerative process. In this review, we discuss biogenesis of microRNAs and piwiRNAs and their functional role in key developmental pathways in vertebrates and invertebrates with an emphasis on recent studies on planarian small RNA pathways.

Essential features and rational design of CRISPR RNAs that function with the Cas RAMP module complex to cleave RNAs.

Small RNAs target invaders for silencing in the CRISPR-Cas pathways that protect bacteria and archaea from viruses and plasmids. The CRISPR RNAs (crRNAs) contain sequence elements acquired from invaders that guide CRISPR-associated (Cas) proteins back to the complementary invading DNA or RNA. Here, we have analyzed essential features of the crRNAs associated with the Cas RAMP module (Cmr) effector complex, which cleaves targeted RNAs. We show that Cmr crRNAs contain an 8 nucleotide 5' sequence tag (also found on crRNAs associated with other CRISPR-Cas pathways) that is critical for crRNA function and can be used to engineer crRNAs that direct cleavage of novel targets. We also present data that indicate that the Cmr complex cleaves an endogenous complementary RNA in Pyrococcus furiosus, providing direct in vivo evidence of RNA targeting by the CRISPR-Cas system. Our findings indicate that the CRISPR RNA-Cmr protein pathway may be exploited to cleave RNAs of interest.

Loss of DNA mismatch repair imparts a selective advantage in planarian adult stem cells.

Lynch syndrome (LS) leads to an increased risk of early-onset colorectal and other types of cancer and is caused by germline mutations in DNA mismatch repair (MMR) genes. Loss of MMR function results in a mutator phenotype that likely underlies its role in tumorigenesis. However, loss of MMR also results in the elimination of a DNA damage-induced checkpoint/apoptosis activation barrier that may allow damaged cells to grow unchecked. A fundamental question is whether loss of MMR provides pre-cancerous stem cells an immediate selective advantage in addition to establishing a mutator phenotype. To test this hypothesis in an in vivo system, we utilized the planarian Schmidtea mediterranea which contains a significant population of identifiable adult stem cells. We identified a planarian homolog of human MSH2, a MMR gene which is mutated in 38% of LS cases. The planarian Smed-msh2 is expressed in stem cells and some progeny. We depleted Smed-msh2 mRNA levels by RNA-interference and found a striking survival advantage in these animals treated with a cytotoxic DNA alkylating agent compared to control animals. We demonstrated that this tolerance to DNA damage is due to the survival of mitotically active, MMR-deficient stem cells. Our results suggest that loss of MMR provides an in vivo survival advantage to the stem cell population in the presence of DNA damage that may have implications for tumorigenesis.

Evolution of alternative and constitutive regions of mammalian 5'UTRs.

BACKGROUND: Alternative splicing (AS) in protein-coding sequences has emerged as an important mechanism of regulation and diversification of animal gene function. By contrast, the extent and roles of alternative events including AS and alternative transcription initiation (ATI) within the 5'-untranslated regions (5'UTRs) of mammalian genes are not well characterized. RESULTS: We evaluated the abundance, conservation and evolution of putative regulatory control elements, namely, upstream start codons (uAUGs) and open reading frames (uORFs), in the 5'UTRs of human and mouse genes impacted by alternative events. For genes with alternative 5'UTRs, the fraction of alternative sequences (those present in a subset of the transcripts) is much greater than that in the corresponding coding sequence, conceivably, because 5'UTRs are not bound by constraints on protein structure that limit AS in coding regions. Alternative regions of mammalian 5'UTRs evolve faster and are subject to a weaker purifying selection than constitutive portions. This relatively weak selection results in over-abundance of uAUGs and uORFs in the alternative regions of 5'UTRs compared to constitutive regions. Nevertheless, even in alternative regions, uORFs evolve under a stronger selection than the rest of the sequences, indicating that some of the uORFs are conserved regulatory elements; some of the non-conserved uORFs could be involved in species-specific regulation. CONCLUSION: The findings on the evolution and selection in alternative and constitutive regions presented here are consistent with the hypothesis that alternative events, namely, AS and ATI, in 5'UTRs of mammalian genes are likely to contribute to the regulation of translation.

Solving the problem of Trans-Genomic Query with alignment tables.

The trans-genomic query (TGQ) problem--enabling the free query of biological information, even across genomes--is a central challenge facing bioinformatics. Solutions to this problem can alter the nature of the field, moving it beyond the jungle of data integration and expanding the number and scope of questions that can be answered. An alignment table is a binary relationship on locations (sequence segments). An important special case of alignment tables are hit tables ? tables of pairs of highly similar segments produced by alignment tools like BLAST. However, alignment tables also include general binary relationships, and can represent any useful connection between sequence locations. They can be curated, and provide a high-quality queryable backbone of connections between biological information. Alignment tables thus can be a natural foundation for TGQ, as they permit a central part of the TGQ problem to be reduced to purely technical problems involving tables of locations.Key challenges in implementing alignment tables include efficient representation and indexing of sequence locations. We define a location datatype that can be incorporated naturally into common off-the-shelf database systems. We also describe an implementation of alignment tables in BLASTGRES, an extension of the open-source POSTGRESQL database system that provides indexing and operators on locations required for querying alignment tables. This paper also reviews several successful large-scale applications of alignment tables for Trans-Genomic Query. Tables with millions of alignments have been used in queries about alternative splicing, an area of genomic analysis concerning the way in which a single gene can yield multiple transcripts. Comparative genomics is a large potential application area for TGQ and alignment tables.

DNA BARCODING IN LAND PLANTS: DEVELOPING STANDARDS TO QUANTIFY AND MAXIMIZE SUCCESS.

The selection of a DNA barcode in plants has been impeded in part due to the relatively low rates of nucleotide substitution observed at the most accessible plastid markers. However, the absence of consensus also reflects a lack of standards for comparing potential barcode markers. While many publications have suggested a host of plant DNA barcodes, the studies cannot be readily compared with each other through any quantitative or statistical parameter, partly because they put forward no single compelling rationale relevant to the adoption of a DNA barcode in plants. Here, we argue that the efficacy of any particular plant DNA barcode selection should reflect the anticipated performance of the resulting barcode database in assignment of a query sequence to species. While legitimate scientific disagreement exists over the criteria relevant to "database performance", the notion gives a unifying rationale for prioritizing selection criteria. Accordingly, we suggest a measure of barcode efficacy based on the rationale of database performance, "the probability of correct identification" (PCI). Moreover, the definition of PCI is left flexible enough to handle most of the scientific disagreement over how to best evaluate DNA barcodes. Finally, we consider how different types of barcodes might require different methods of analysis and database design and indicate how the analysis might affect the selection of the most broadly effective barcode for land plants.

Widespread positive selection in synonymous sites of mammalian genes.

Evolution of protein sequences is largely governed by purifying selection, with a small fraction of proteins evolving under positive selection. The evolution at synonymous positions in protein-coding genes is not nearly as well understood, with the extent and types of selection remaining, largely, unclear. A statistical test to identify purifying and positive selection at synonymous sites in protein-coding genes was developed. The method compares the rate of evolution at synonymous sites (Ks) to that in intron sequences of the same gene after sampling the aligned intron sequences to mimic the statistical properties of coding sequences. We detected purifying selection at synonymous sites in approximately 28% of the 1,562 analyzed orthologous genes from mouse and rat, and positive selection in approximately 12% of the genes. Thus, the fraction of genes with readily detectable positive selection at synonymous sites is much greater than the fraction of genes with comparable positive selection at nonsynonymous sites, i.e., at the level of the protein sequence. Unlike other genes, the genes with positive selection at synonymous sites showed no correlation between Ks and the rate of evolution in nonsynonymous sites (Ka), indicating that evolution of synonymous sites under positive selection is decoupled from protein evolution. The genes with purifying selection at synonymous sites showed significant anticorrelation between Ks and expression level and breadth, indicating that highly expressed genes evolve slowly. The genes with positive selection at synonymous sites showed the opposite trend, i.e., highly expressed genes had, on average, higher Ks. For the genes with positive selection at synonymous sites, a significantly lower mRNA stability is predicted compared to the genes with negative selection. Thus, mRNA destabilization could be an important factor driving positive selection in nonsynonymous sites, probably, through regulation of expression at the level of mRNA degradation and, possibly, also translation rate. So, unexpectedly, we found that positive selection at synonymous sites of mammalian genes is substantially more common than positive selection at the level of protein sequences. Positive selection at synonymous sites might act through mRNA destabilization affecting mRNA levels and translation.

Assessing the impact of alternative splicing on domain interactions in the human proteome.

We have constructed a database of alternatively spliced protein forms (ASP), consisting of 13,384 protein isoform sequences of 4422 human genes (www.bioinformatics.ucla.edu/ASP). We identified fifty protein domain types that were selectively removed by alternative splicing at much higher frequencies than average (p-value < 0.01). These include many well-known protein-interaction domains (e.g., KRAB; ankyrin repeats; Kelch) including some that have been previously shown to be regulated functionally by alternative splicing (e.g., collagen domain). We present a number of novel examples (Kruppel transcription factors; Pbx2; Enc1) from the ASP database, illustrating how this pattern of alternative splicing changes the structure of a biological pathway, by redirecting protein interaction networks at key switch points. Our bioinformatics analysis indicates that a major impact of alternative splicing is removal of protein-protein interaction domains that mediate key linkages in protein interaction networks. ASP expands the available dataset of human alternatively spliced protein forms from 1989 human genes (SwissProt release 42) to 5413 (nonredundant set, ASP + SwissProt), a nearly 3-fold increase. ASP will enhance the existing pool of protein sequences that are searched by mass spectroscopy software during the identification of peptide fragments.

Evidence for a subpopulation of conserved alternative splicing events under selection pressure for protein reading frame preservation.

Recently there has been much interest in assessing the role of alternative splicing in evolution. We have sought to measure functional selection pressure on alternatively spliced single-exon skips, by calculating the fraction that are an exact multiple of 3 nt in length and therefore preserve protein reading-frame in both the exon-inclusion and exon-skip splice forms. The frame-preservation ratio (defined as the number of exons that are an exact multiple of three in length, divided by the number of exons that are not) was slightly above random for both constitutive exons and alternatively spliced exons as a whole in human and mouse. However, orthologous exons that were observed to be alternatively spliced in the expressed sequence tag data from two or more organisms showed a substantially increased bias to be frame-preserving. This effect held true only for exons within the protein coding region, and not the untranslated region. In five animal genomes (human, mouse, rat, zebrafish, Drosophila), we observed an association between these conserved alternative splicing events and increased selection pressure for frame-preservation. Surprisingly, this effect became stronger as a function of decreasing exon inclusion level: for alternatively spliced exons that were included in a majority of the gene's transcripts, the frame-preservation bias was no higher than that of constitutive exons, whereas for alternatively spliced exons that were included in only a minority of the gene's transcripts, the frame-preservation bias increased nearly 20-fold. These data indicate that a subpopulation of modern alternative splicing events was present in the common ancestors of these genomes, and was under functional selection pressure to preserve the protein reading frame.

The multiassembly problem: reconstructing multiple transcript isoforms from EST fragment mixtures.

Recent evidence of abundant transcript variation (e.g., alternative splicing, alternative initiation, alternative polyadenylation) in complex genomes indicates that cataloging the complete set of transcripts from an organism is an important project. One challenge is the fact that most high-throughput experimental methods for characterizing transcripts (such as EST sequencing) give highly detailed information about short fragments of transcripts or protein products, instead of a complete characterization of a full-length form. We analyze this "multiassembly problem"-reconstructing the most likely set of full-length isoform sequences from a mixture of EST fragment data-and present a graph-based algorithm for solving it. In a variety of tests, we demonstrate that this algorithm deals appropriately with coupling of distinct alternative splicing events, increasing fragmentation of the input data and different types of transcript variation (such as alternative splicing, initiation, polyadenylation, and intron retention). To test the method's performance on pure fragment (EST) data, we removed all mRNA sequences, and found it produced no errors in 40 cases tested. Using this algorithm, we have constructed an Alternatively Spliced Proteins database (ASP) from analysis of human expressed and genomic sequences, consisting of 13,384 protein isoforms of 4422 genes, yielding an average of 3.0 protein isoforms per gene.