RNA-seq Experiments Research Articles

Investigating the role of the ROS/CncC signaling pathway in the response to xenobiotics in Spodoptera frugiperda using Sf9 cells

Spodoptera frugiperda (fall armyworm, FAW) is an invasive polyphagous lepidopteran pest that has developed sophisticated resistance mechanisms involving detoxification enzymes to eliminate toxic compounds it encounters in its diet including insecticides. Although its inventory of detoxification enzymes is known, the mechanisms that enable an adapted response depending on the toxic compound remain largely unexplored. Sf9 cells were used to investigate the role of the transcription factors, Cap n′ collar isoform C (CncC) and musculoaponeurotic fibrosarcoma (Maf) in the regulation of the detoxification response. We overexpressed CncC, Maf or both genes, and knocked out (KO) CncC or its repressor Kelch-like ECH associated protein 1 (Keap1). Joint overexpression of CncC and Maf is required to confer increased tolerance to indole 3-carbinol (I3C), a plant secondary metabolite, and to methoprene, an insecticide. Both molecules induce reactive oxygen species (ROS) pulses in the different cell lines. The use of an antioxidant reversed ROS pulses and restored the tolerance to I3C and methoprene. The activity of detoxification enzymes varied according to the cell line. Suppression of Keap1 significantly increased the activity of cytochrome P450s, carboxylesterases and glutathione S-transferases. RNAseq experiments showed that CncC mainly regulates the expression of detoxification genes but is also at the crossroads of several signaling pathways (reproduction and immunity) maintaining homeostasis. We present new data in Sf9 cell lines suggesting that the CncC:Maf pathway plays a central role in FAW response to natural and synthetic xenobiotics. This knowledge helps to better understand detoxification gene expression and may help to design next-generation pest insect control measures.

Bacterial-fungal interactions promote parallel evolution of global transcriptional regulators in a widespread Staphylococcus species.

Experimental studies of microbial evolution have largely focused on monocultures of model organisms, but most microbes live in communities where interactions with other species may impact rates and modes of evolution. Using the cheese rind model microbial community, we determined how species interactions shape the evolution of the widespread food- and animal-associated bacterium Staphylococcus xylosus. We evolved S. xylosus for 450 generations alone or in co-culture with one of three microbes: the yeast Debaryomyces hansenii, the bacterium Brevibacterium aurantiacum, and the mold Penicillium solitum. We used the frequency of colony morphology mutants (pigment and colony texture phenotypes) and whole-genome sequencing of isolates to quantify phenotypic and genomic evolution. The yeast D. hansenii strongly promoted diversification of S. xylosus. By the end of the experiment, all populations co-cultured with the yeast were dominated by pigment and colony morphology mutant phenotypes. Populations of S. xylosus grown alone, with B. aurantiacum, or with P. solitum did not evolve novel phenotypic diversity. Whole-genome sequencing of individual mutant isolates across all four treatments identified numerous unique mutations in the operons for the SigB, Agr, and WalRK global regulators, but only in the D. hansenii treatment. Phenotyping and RNA-seq experiments highlighted altered pigment and biofilm production, spreading, stress tolerance, and metabolism ofS. xylosus mutants.Fitness experiments revealed antagonistic pleiotropy, where beneficial mutations that evolved in the presence of the yeast had strong negative fitness effects in other biotic environments. This work demonstrates that bacterial-fungal interactions can have long-term evolutionary consequences within multispecies microbiomes by facilitating the evolution of strain diversity.

AgAnimalGenomes: browsers for viewing and manually annotating farm animal genomes.

Current genome sequencing technologies have made it possible to generate highly contiguous genome assemblies for non-model animal species. Despite advances in genome assembly methods, there is still room for improvement in the delineation of specific gene features in the genomes. Here we present genome visualization and annotation tools to support seven livestock species (bovine, chicken, goat, horse, pig, sheep, and water buffalo), available in a new resource called AgAnimalGenomes. In addition to supporting the manual refinement of gene models, these browsers provide visualization tracks for hundreds of RNAseq experiments, as well as data generated by the Functional Annotation of Animal Genomes (FAANG) Consortium. For species with predicted gene sets from both Ensembl and RefSeq, the browsers provide special tracks showing the thousands of protein-coding genes that disagree across the two gene sources, serving as a valuable resource to alert researchers to gene model issues that may affect data interpretation. We describe the data and search methods available in the new genome browsers and how to use the provided tools to edit and create new gene models.

Crosstalk between KIF1C and PRKAR1A in left atrial myxoma

Cardiac myxoma (CM) is the most common benign cardiac tumor, and most CMs are left atrial myxomas (LAMs). Six variations of KIF1C, c.899 A > T, c.772 T > G, c.352 A > T, c.2895 C > T, c.3049 G > A, and c.*442_*443dup in left atrial myxoma tissues are identified by whole-exome sequencing (WES) and Sanger sequencing. RNA-seq and function experiments show the reduction of the expression of KIF1C and PRKAR1A caused by rare variations of KIF1C. KIF1C is observed to be located in the nucleus, bind to the promoter region of PRKAR1A, and regulate its transcription. Reduction of KIF1C decreases PRKAR1A expression and activates the PKA, which causes an increase in ERK1/2 phosphorylation and SRC-mediated STAT3 activation, a reduction of CDH1, TP53, CDKN1A, and BAX, and eventually promotes tumor formation both in vitro and in vivo. The results suggest that inhibition of KIF1C promotes the pathogenesis of LAM through positive feedback formed by the crosstalk between KIF1C and PRKAR1A.

Design, execution, and interpretation of plant RNA-seq analyses.

Genomics has transformed our understanding of the genetic architecture of traits and the genetic variation present in plants. Here, we present a review of how RNA-seq can be performed to tackle research challenges addressed by plant sciences. We discuss the importance of experimental design in RNA-seq, including considerations for sampling and replication, to avoid pitfalls and wasted resources. Approaches for processing RNA-seq data include quality control and counting features, and we describe common approaches and variations. Though differential gene expression analysis is the most common analysis of RNA-seq data, we review multiple methods for assessing gene expression, including detecting allele-specific gene expression and building co-expression networks. With the production of more RNA-seq data, strategies for integrating these data into genetic mapping pipelines is of increased interest. Finally, special considerations for RNA-seq analysis and interpretation in plants are needed, due to the high genome complexity common across plants. By incorporating informed decisions throughout an RNA-seq experiment, we can increase the knowledge gained.

Foreign RNA spike-ins enable accurate allele-specific expression analysis at scale.

Analysis of allele-specific expression is strongly affected by the technical noise present in RNA-seq experiments. Previously, we showed that technical replicates can be used for precise estimates of this noise, and we provided a tool for correction of technical noise in allele-specific expression analysis. This approach is very accurate but costly due to the need for two or more replicates of each library. Here, we develop a spike-in approach which is highly accurate at only a small fraction of the cost. We show that a distinct RNA added as a spike-in before library preparation reflects technical noise of the whole library and can be used in large batches of samples. We experimentally demonstrate the effectiveness of this approach using combinations of RNA from species distinguishable by alignment, namely, mouse, human, and Caenorhabditis elegans. Our new approach, controlFreq, enables highly accurate and computationally efficient analysis of allele-specific expression in (and between) arbitrarily large studies at an overall cost increase of ∼5%. Analysis pipeline for this approach is available at GitHub as R package controlFreq (github.com/gimelbrantlab/controlFreq).

A contamination focused approach for optimizing the single-cell RNA-seq experiment

Droplet-based single-cell RNA-seq (scRNA-seq) data are plagued by ambient contaminations caused by nucleic acid material released by dead and dying cells. This material is mixed into the buffer and is co-encapsulated with cells, leading to a lower signal-to-noise ratio. Although there exist computational methods to remove ambient contaminations post-hoc, the reliability of algorithms in generating high-quality data from low-quality sources remains uncertain. Here, we assess data quality before data filtering by a set of quantitative, contamination-based metrics that assess data quality more effectively than standard metrics. Through a series of controlled experiments, we report improvements that can minimize ambient contamination outside of tissue dissociation, via cell fixation, improved cell loading, microfluidic dilution, and nuclei versus cell preparation; many of these parameters are inaccessible on commercial platforms. We provide end-users with insights on factors that can guide their decision-making regarding optimizations that minimize ambient contamination, and metrics to assess data quality.

Oligo-dT anchored cDNA-SRAP and cDNA-SCoT aided identification of transcripts differentially expressed during the early stages of recovery of resurrection plant Haberlea rhodopensis Friv. from freezing-induced desiccation

Resurrection plants are a unique group of flora showing an extraordinary ability to survive extreme water shortages. Among them, the Bulgarian endemic plant Haberlea rhodopensis Friv. has been used as a model to study the mechanisms underlying desiccation tolerance induced by water and low temperature stress. Although extensive transcriptome data, including two RNAseq experiments, are available for drought induced desiccation and recovery, the gene expression during recovery after freezing-induced desiccation (RAF) remains largely unexplored. In this study we performed a differential screening using a new approach, based on oligo-dT anchored cDNA-SRAP (copy DNA sequence-related amplified polymorphism) and a modified oligo-dT anchored cDNA-SCoT (copy DNA start codon targeted) assays to provide new information on genes involved in the tolerance to low temperature stress and the recovery of H. rhodopensis. The putative functions of 19 identified transcript derived fragments were established by searching for homology with known sequences in public databases and the expression of 10 of them during the early stages of recovery was confirmed by quantitative reverse transcription polymerase chain reaction. Some of the identified genes have known functions in the response to abiotic stress in other plant species including resurrection plants, which confirms the effectiveness of the approach used in the present study. The results of this study will contribute to obtaining new information on the defence mechanisms in the recovery process of H. rhodopensis, especially at the early stages of recovery, which is the most vulnerable period for plants.

NUSAP1-LDHA-Glycolysis-Lactate feedforward loop promotes Warburg effect and metastasis in pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is characterized by hypoxia and hypovascular tumor microenvironment. Nucleolar and spindle associated protein 1 (NUSAP1) is a microtubule-associated protein that is known to be involved in cancer biology. Our study aimed to investigate the role of NUSAP1 in glycolytic metabolism and metastasis in PDAC. Expression and prognostic value of NUSAP1 in PDAC and common gastrointestinal tumors was evaluated. The function of NUSAP1 in PDAC progression was clarified by single-cell RNA-seq and further experiments in vitro, xenograft mouse model, spontaneous PDAC mice model and human tissue microarray. The downstream genes and signaling pathways regulated by NUSAP1 were explored by RNA-Seq. And the regulation of NUSAP1 on Lactate dehydrogenase A (LDHA)-mediated glycolysis and its underlying mechanism was further clarified by CHIP-seq. NUSAP1 was an independent unfavorable predictor of PDAC prognosis that playing a critical role in metastasis of PDAC by regulating LDHA-mediated glycolysis. Mechanically, NUSAP1 could bind to c-Myc and HIF-1α that forming a transcription regulatory complex localized to LDHA promoter region and enhanced its expression. Intriguingly, lactate upregulated NUSAP1 expression by inhibiting NUSAP1 protein degradation through lysine lactylated (Kla) modification, thus forming a NUSAP1-LDHA-glycolysis-lactate feedforward loop. The NUSAP1-LDHA-glycolysis-lactate feedforward loop is one of the underlying mechanisms to explain the metastasis and glycolytic metabolic potential in PDAC, which also provides a novel insights to understand the Warburg effect in cancer. Targeting NUSAP1 would be an attractive paradigm for PDAC treatment.

The effect of background noise and its removal on the analysis of single-cell expression data

BackgroundIn droplet-based single-cell and single-nucleus RNA-seq experiments, not all reads associated with one cell barcode originate from the encapsulated cell. Such background noise is attributed to spillage from cell-free ambient RNA or barcode swapping events.ResultsHere, we characterize this background noise exemplified by three scRNA-seq and two snRNA-seq replicates of mouse kidneys. For each experiment, cells from two mouse subspecies are pooled, allowing to identify cross-genotype contaminating molecules and thus profile background noise. Background noise is highly variable across replicates and cells, making up on average 3–35% of the total counts (UMIs) per cell and we find that noise levels are directly proportional to the specificity and detectability of marker genes. In search of the source of background noise, we find multiple lines of evidence that the majority of background molecules originates from ambient RNA. Finally, we use our genotype-based estimates to evaluate the performance of three methods (CellBender, DecontX, SoupX) that are designed to quantify and remove background noise. We find that CellBender provides the most precise estimates of background noise levels and also yields the highest improvement for marker gene detection. By contrast, clustering and classification of cells are fairly robust towards background noise and only small improvements can be achieved by background removal that may come at the cost of distortions in fine structure.ConclusionsOur findings help to better understand the extent, sources and impact of background noise in single-cell experiments and provide guidance on how to deal with it.

Unraveling the timeline of gene expression: A pseudotemporal trajectory analysis of single-cell RNA sequencing data

Background: Single-cell RNA sequencing (scRNA-seq) technologies have rapidly developed in recent years. The droplet-based single cell platforms enable the profiling of gene expression in tens of thousands of cells per sample. The goal of a typical scRNA-seq analysis is to identify different cell subpopulations and their respective marker genes. Additionally, trajectory analysis can be used to infer the developmental or differentiation trajectories of cells. Methods: This article demonstrates a comprehensive workflow for performing trajectory inference and time course analysis on a multi-sample single-cell RNA-seq experiment of the mouse mammary gland. The workflow uses open-source R software packages and covers all steps of the analysis pipeline, including quality control, doublet prediction, normalization, integration, dimension reduction, cell clustering, trajectory inference, and pseudo-bulk time course analysis. Sample integration and cell clustering follows the Seurat pipeline while the trajectory inference is conducted using the monocle3 package. The pseudo-bulk time course analysis uses the quasi-likelihood framework of edgeR. Results: Cells are ordered and positioned along a pseudotime trajectory that represented a biological process of cell differentiation and development. The study successfully identified genes that were significantly associated with pseudotime in the mouse mammary gland. Conclusions: The demonstrated workflow provides a valuable resource for researchers conducting scRNA-seq analysis using open-source software packages. The study successfully demonstrated the usefulness of trajectory analysis for understanding the developmental or differentiation trajectories of cells. This analysis can be applied to various biological processes such as cell development or disease progression, and can help identify potential biomarkers or therapeutic targets.

#3769 MECHANOSENSITIVE ION CHANNEL PIEZO1 MEDIATES BLOOD PRESSURE HOMEOSTASIS BY REGULATING RENIN EXPRESSION IN MICE

Abstract Background and Aims The renin-angiotensin system (RAS) plays a crucial role in the maintenance of fluid-electrolyte and blood pressure homeostasis. Renin, the step-limiting protease of the RAS, is produced, stored, and secreted by the renal juxtaglomerular (JG) cells, which are closely regulated by the blood pressure. Studies have shown that local blood flow-associated shear stress play an important role in renin expression and secretion. To date, however, the molecular identity through which JG cells respond to the blood pressure remains unclear. Method In this study, we determined the expression and subcellular localization of Piezo1 in JG cells by qPCR, Western blot, Immunofluorescence (IF), and multiplex immunohistochemistry. Then, a Fluo-4 AM-based calcium-imaging system was used to detect the dynamic changes of intracellular calcium in response to Piezo1 specific agonist Yoda1 and mechanical stress (MS)-induced by perfusion. Next, CRISPR/Cas9 technology was used to generate Piezo1 knockout JG cells to further confirm its contribution. Then, we used RNA-seq experiment to investigate the downstream signaling of Piezo1 in JG cells. Finally, we generated adeno-associated virus (AAV)-mediated kidney-specific Piezo1 knockdown mice to investigate the in vivo effect of Piezo1. Results We found that the mechanosensitive ion channel Piezo1 was co-expressed with renin in JG cells in mouse kidney slides as well as As4.1 cells, a commonly used JG cell line. Activation of Piezo1 by its agonist Yoda1 induced an intracellular calcium increase and reduced the expression of renin in these cells, while knockout of Piezo1 in JG cells abolished the effect of Yoda1. Meanwhile, fluid induced mechanical stress also led to an intracellular calcium increase and renin downregulation in wildtype but not Piezo1 knockout JG cells. Mechanistically, we demonstrated that activation of Piezo1 upregulated the Ptgs2 expression via the calcineurin-NFAT pathway and increased the production of Ptgs2 downstream molecule COX-2 and PGE2 in JG cells, which inhibited cAMP production and reduced renin expression in JG cells. In mice, we demonstrated that activation of Piezo1 significantly downregulated the blood pressure in wildtype but not kidney-specific Piezo1 knockdown mice. Conclusion In summary, these results revealed that activation of Piezo1 could lead to downregulation of the renin expression in vitro and in vivo, subsequently reduction of blood pressure, highlighting its therapeutic potential as a drug target of the renin-angiotensin system.

RNA Sequencing of Pooled Samples Effectively Identifies Differentially Expressed Genes.

Analysis of gene expression changes across the genome provides a powerful, unbiased tool for gaining insight into molecular mechanisms. We have effectively used RNA sequencing to identify differentially expressed genes in long-lived genetic mutants in C. elegans to advance our understanding of the genetic pathways that control longevity. Although RNA sequencing costs have come down, cost remains a barrier to examining multiple strains and time points with a sufficient number of biological replicates. To circumvent this, we have examined the efficacy of identifying differentially expressed genes by sequencing a pooled RNA sample from long-lived isp-1 mitochondrial mutant worms. We found that sequencing a pooled RNA sample could effectively identify genes that were found to be significantly upregulated in the two individually sequenced RNA-seq experiments. Finally, we compared the genes significantly upregulated in the two individually sequenced RNA-seq experiments to two previous microarray experiments to come up with a high-confidence list of modulated genes in long-lived isp-1 mutant worms. Overall, this work demonstrates that RNA sequencing of pooled RNA samples can be used to identify differentially expressed genes.

A multi-omic resource of wheat seed tissues for nutrient deposition and improvement for human health

As a globally important staple crop, wheat seeds provide us with nutrients and proteins. The trend of healthy dietary has become popular recently, emphasizing the consumption of whole-grain wheat products and the dietary benefits. However, the dynamic changes in nutritional profiles of different wheat seed regions (i.e., the embryo, endosperm and outer layers) during developmental stages and the molecular regulation have not been well studied. Here, we provide this multi-omic resource of wheat seeds and describe the generation, technical assessment and preliminary analyses. This resource includes a time-series RNA-seq dataset of the embryo, endosperm and outer layers of wheat seeds and their corresponding metabolomic dataset, covering the middle and late stages of seed development. Our RNA-seq experiments profile the expression of 63,708 genes, while the metabolomic data includes the abundance of 984 metabolites. We believe that this was the first reported transcriptome and metabolome dataset of wheat seeds that helps understand the molecular regulation of the deposition of beneficial nutrients and hence improvements for nutritional and processing quality traits.

Deficiency of Wdr60 and Wdr34 cause distinct neural tube malformation phenotypes in early embryos.

Cilia are specialized organelles that extend from plasma membrane, functioning as antennas for signal transduction and are involved in embryonic morphogenesis. Dysfunction of cilia lead to many developmental defects, including neural tube defects (NTDs). Heterodimer WDR60-WDR34 (WD repeat domain 60 and 34) are intermediate chains of motor protein dynein-2, which play important roles in ciliary retrograde transport. It has been reported that disruption of Wdr34 in mouse model results in NTDs and defects of Sonic Hedgehog (SHH) signaling. However, no Wdr60 deficiency mouse model has been reported yet. In this study, piggyBac (PB) transposon is used to interfere Wdr60 and Wdr34 expression respectively to establish Wdr60 PB/PB and Wdr34 PB/PB mouse models. We found that the expression of Wdr60 or Wdr34 is significantly decreased in the homozygote mice. Wdr60 homozygote mice die around E13.5 to E14.5, while Wdr34 homozygote mice die around E10.5 to E11.5. WDR60 is highly expressed in the head region at E10.5 and Wdr60 PB/PB embryos have head malformation. RNAseq and qRT-PCR experiments revealed that Sonic Hedgehog signaling is also downregulated in Wdr60 PB/PB head tissue, demonstrating that WDR60 is also required for promoting SHH signaling. Further experiments on mouse embryos also revealed that the expression levels of planar cell polarity (PCP) components such as CELSR1 and downstream signal molecule c-Jun were downregulated in WDR34 homozygotes compared to wildtype littermates. Coincidently, we observed much higher ratio of open cranial and caudal neural tube in Wdr34 PB/PB mice. CO-IP experiment showed that WDR60 and WDR34 both interact with IFT88, but only WDR34 interacts with IFT140. Taken together, WDR60 and WDR34 play overlapped and distinct functions in modulating neural tube development.

Phased grapevine genome sequence of an Rpv12 carrier for biotechnological exploration of resistance to Plasmopara viticola.

The downy mildew disease caused by the oomycete Plasmopara viticola is a serious threat for grapevine and can cause enormous yield losses in viticulture. The quantitative trait locus Rpv12, mediating resistance against P. viticola, was originally found in Asian Vitis amurensis. This locus and its genes were analyzed here in detail. A haplotype-separated genome sequence of the diploid Rpv12-carrier Gf.99-03 was created and annotated. The defense response against P. viticola was investigated in an infection time-course RNA-seq experiment, revealing approximately 600 upregulated Vitis genes during host-pathogen interaction. The Rpv12 regions of the resistance and the sensitivity encoding Gf.99-03 haplotype were structurally and functionally compared with each other. Two different clusters of resistance-related genes were identified within the Rpv12 locus. One cluster carries a set of four differentially expressed genes with three ACCELERATED CELL DEATH 6-like genes. The other cluster carries a set of six resistance gene analogs related to qualitative pathogen resistance. The Rpv12 locus and its candidate genes for P. viticola resistance provide a precious genetic resource for P. viticola resistance breeding. Newly developed co-segregating simple sequence repeat markers in close proximity to the R-genes enable its improved applicability in marker-assisted grapevine breeding.

Combined analyses of transcriptome and metabolome reveal the mechanism of exogenous strigolactone regulating the response of elephant grass to drought stress.

Elephant grass is widely used in feed production and ecological restoration because of its huge biomass and low occurrence of diseases and insect pets. However, drought seriously affects growth and development of this grass. Strigolactone (SL), a small molecular phytohormone, reportedly participates in improving resilience to cope with arid environment. But the mechanism of SL regulating elephant grass to response to drought stress remains unknown and needs further investigation. We conducted RNA-seq experiments and identified 84,296 genes including 765 and 2325 upregulated differential expression genes (DEGs) and 622 and 1826 downregulated DEGs, compared drought rehydration with spraying SL in roots and leaves, respectively. Combined with targeted phytohormones metabolite analysis, five hormones including 6-BA, ABA, MeSA, NAA, and JA had significant changes under re-watering and spraying SL stages. Moreover, a total of 17 co-expression modules were identified, of which eight modules had the most significant correlation with all physiological indicators with weighted gene co-expression network analysis. The venn analysis revealed the common genes between Kyoto Encyclopedia of Genes and Genomes enriched functional DEGs and the top 30 hub genes of higher weights in eight modules, respectively. Finally, 44 DEGs had been identified as key genes which played a major role in SL response to drought stress. After verification of its expression level by qPCR, six key genes in elephant grass including PpPEPCK, PpRuBPC, PpPGK, PpGAPDH, PpFBA, and PpSBPase genes regulated photosynthetic capacity under the SL treatment to respond to drought stress. Meanwhile, PpACAT, PpMFP2, PpAGT2, PpIVD, PpMCCA, and PpMCCB regulated root development and phytohormone crosstalk to respond to water deficit conditions. Our research led to a more comprehensive understanding about exogenous SL that plays a role in elephant grass response to drought stress and revealed insights into the SL regulating molecular mechanism in plants to adapt to the arid environment.

Splicing Characterization and Isoform Switch Events in Human Keratinocytes Carrying Oncogenes from High-Risk HPV-16 and Low-Risk HPV-84.

Infection of epithelial cells with high-risk HPV (HR-HPV) types, followed by expression of virus oncogenic proteins (E5, E6, and E7), leads to genomic imbalance, suppression of tumor inhibitors, and induction of oncogenes. Low-risk HPV (LR-HPV) may slow the rate at which cervical cancer spreads to an invasive stage since co-infection with LR-HPV is linked to a decreased risk of future invasive cancer than infection with HR-HPV alone. We then propose that cancer-progressing changes may be distinguished through identifying the functional differences between LR-HPV and HR-HPV. Lentiviral strategies were followed to establish HaCaT cells with constitutive expression of HPV oncogenes. RNAseq experiments were designed to analyze the transcriptome modulations caused by each of the E5, E6, and E7 oncogenes of HPV-16 and HPV-84 in HaCaT cells. We identified enhanced RNA degradation, spliceosome, and RNA polymerase pathways related to mRNA processing. ATTS (alternative transcription termination site) was discovered to be more prevalent in cells with HPV-16E5 than HPV-84E5. In HPV-16E6-infected cells, ATTS gain was significantly higher than ATTS loss. Cells with HPV-16E7 had more isoforms with intron retention (IR) than those with HPV-84E7. We identified switches in ADAM10, CLSPN, and RNPS1 that led to greater expression of the coding isoforms in HR-HPV. The results of this work highlight differences between LR-HPV and HR-HPV in mRNA processing. Moreover, crucial cervical cancer-related switch events were detected.

RNAcontacts: A Pipeline for Predicting Contacts from RNA Proximity Ligation Assays.

High-throughput RNA proximity ligation assays are molecular methods that are used to simultaneously analyze the spatial proximity of many RNAs in living cells. Their principle is based on cross-linking, fragmentation, and subsequent religation of RNAs, followed by high-throughput sequencing. The generated fragments have two different types of splits, one resulting from pre-mRNA splicing and the other formed by the ligation of spatially close RNA strands. Here, we present RNAcontacts, a universal pipeline for detecting RNA-RNA contacts in high-throughput RNA proximity ligation assays. RNAcontacts circumvents the inherent problem of mapping sequences with two distinct types of splits using a two-pass alignment, in which splice junctions are inferred from a control RNA-seq experiment on the first pass and then provided to the aligner as bona fide introns on the second pass. Compared to previously developed methods, our approach allows for a more sensitive detection of RNA contacts and has a higher specificity with respect to splice junctions that are present in the biological sample. RNAcontacts automatically extracts contacts, clusters their ligation points, computes the read support, and generates tracks for visualizing through the UCSC Genome Browser. The pipeline is implemented in Snakemake, a reproducible and scalable workflow management system for rapid and uniform processing of multiple datasets. RNAcontacts is a generic pipeline for the detection of RNA contacts that can be used with any proximity ligation method as long as one of the interacting partners is RNA. RNAcontacts is available via the GitHub repository https://github.com/smargasyuk/ RNAcontacts/.

An old dog with ‘new’ tricks: the MDR1 transporter in T cell-mediated immunity

Abstract Multidrug Resistance Protein 1 (MDR1) is an orphan ATP-binding cassette transporter, originally discovered for removing cytostatic drugs from resistant tumor cells. However, MDR1 is also functionally expressed in many normal cell types. Work from our group has shown that MDR1 is constitutively expressed in cytotoxic T cells (CTLs), where it promotes CTL-mediated immunity to infection. Adoptive transfer of Lymphocytic choriomeningitis virus (LCMV)-specific MDR1 (Abcb1a/b)-deficient naïve CTLs into wild-type recipients, followed by LCMV infection, has revealed that MDR1-deficient CTLs fail to accumulate to wild-type levels in response to either acute (Armstrong) or chronic (Clone-13) LCMV infection. This defective accumulation is observed as early as 5-days post-infection, suggesting that endogenous MDR1 transport activity may preferentially regulate early T cell receptor (TCR)-elicited CTL functions. New single cell (sc)RNA-seq experiments, coupled with RNA velocity and lineage trajectory analyses, suggest that MDR1 expression is induced rapidly during TCR activation, coinciding with the induction of other genes involved in mitochondrial and redox metabolism. In vitro-activated CTLs lacking MDR1 also experience oxidative stress, mitochondrial dysfunction and display mitochondrial localization of MDR1. Together, these data support a new model of MDR1 function, in which MDR1-dependent transport of endogenous substrates suppresses oxidative stress and supports redox homeostasis during periods of metabolic activation. In turn, these concepts provide a novel framework for improving the efficacy and durability of vaccines, and for elucidating metabolic vulnerabilities of multidrug-resistant cancers.