Evaluation of ultra-low input RNA sequencing for the study of human T cell transcriptome

Comparative analysis of cell-cell communication at single-cell resolution

Inference of cell-cell communication from single-cell RNA sequencing data is a powerful technique to uncover intercellular communication pathways, yet existing methods perform this analysis at the level of the cell type or cluster, discarding single-cell-level information. Here we present Scriabin, a flexible and scalable framework for comparative analysis of cell-cell communication at single-cell resolution that is performed without cell aggregation or downsampling. We use multiple published atlas-scale datasets, genetic perturbation screens and direct experimental validation to show that Scriabin accurately recovers expected cell-cell communication edges and identifies communication networks that can be obscured by agglomerative methods. Additionally, we use spatial transcriptomic data to show that Scriabin can uncover spatial features of interaction from dissociated data alone. Finally, we demonstrate applications to longitudinal datasets to follow communication pathways operating between timepoints. Our approach represents a broadly applicable strategy to reveal the full structure of niche-phenotype relationships in health and disease.

Dissecting the Immune System through Gene Regulation

The immune system plays a dual role in human health, functioning both as a protector against pathogens and, at times, as a contributor to disease. This feature emphasizes the importance to uncover the underlying causes of its malfunctions, necessitating an in-depth analysis in both pathological and physiological conditions to better understand the immune system and immune disorders. Recent advances in scientific technology have enabled extensive investigations into gene regulation, a crucial mechanism governing cellular functionality. Studying gene regulatory mechanisms within the immune system is a promising avenue for enhancing our understanding of immune cells and the immune system as a whole. The gene regulatory mechanisms, revealed through various methodologies, and their implications in the field of immunology are discussed in this chapter.

Gene expression changes in therapeutic ultrasound-treated venous leg ulcers

Introduction

Low-frequency, low-intensity ultrasound has been previously shown to promote healing of chronic wounds in humans, but mechanisms behind these effects are poorly understood. The purpose of this study was to evaluate gene expression differences in debrided human venous ulcer tissue from patients treated with low-frequency (20 kHz), low-intensity (100 mW/cm2) ultrasound compared to a sham treatment in an effort to better understand the potential biological mechanisms.

Methods

Debrided venous ulcer tissue was collected from 32 subjects one week after sham treatment or low-frequency, low-intensity ultrasound treatment. Of these samples, 7 samples (3 ultrasound treated and 4 sham treated) yielded sufficient quality total RNA for analysis by ultra-high multiplexed PCR (Ampliseq) and expression of more than 24,000 genes was analyzed. 477 genes were found to be significantly differentially expressed between the ultrasound and sham groups using cut-off values of p < 0.05 and fold change of 2.

Results and Discussion

The top differentially expressed genes included those involved in regulation of cell metabolism, proliferation, and immune cell signaling. Gene set enrichment analysis identified 20 significantly enriched gene sets from upregulated genes and 4 significantly enriched gene sets from downregulated genes. Most of the enriched gene sets from upregulated genes were related to cell-cell signaling pathways. The most significantly enriched gene set from downregulated genes was the inflammatory response gene set. These findings show that therapeutic ultrasound influences cellular behavior in chronic wounds as early as 1 week after application. Considering the well-known role of chronic inflammation in impairing wound healing in chronic wounds, these results suggest that a downregulation of inflammatory genes is a possible biological mechanism of ultrasound-mediated venous chronic wound healing. Such increased understanding may ultimately lead to the enhancement of ultrasound devices to accelerate chronic wound healing and increase patient quality of life.

Meso-seq for in-depth transcriptomics in ultra-low amounts of FACS-purified neuronal nuclei

Profiling of gene expression in sparse populations of genetically defined neurons is essential for dissecting the molecular mechanisms that control the development and plasticity of neural circuits. However, current transcriptomic approaches are ill suited for detailed mechanistic studies in sparse neuronal populations, as they either are technically complex and relatively expensive (e.g., single-cell RNA sequencing [RNA-seq]) or require large amounts of input material (e.g., traditional bulk RNA-seq). Thus, we established Meso-seq, a meso-scale protocol for identifying more than 10,000 robustly expressed genes in as little as 50 FACS-sorted neuronal nuclei. We demonstrate that Meso-seq works well for multiple neuroscience applications, including transcriptomics in antibody-labeled cortical neurons in mice and non-human primates, analyses of experience-regulated gene programs, and RNA-seq from visual cortex neurons labeled ultra-sparsely with viruses. Given its simplicity, robustness, and relatively low costs, Meso-seq is well suited for molecular-mechanistic studies in ultra-sparse neuronal populations in the brain.

HIV-1 Drug Resistance Assay Using Ion Torrent Next Generation Sequencing and On-Instrument End-to-End Analysis Software

HIV-1 antiretroviral therapy management requires sequencing the protease, reverse transcriptase, and integrase portions of the HIV-1 pol gene. Most resistance testing is performed with Sanger sequencing, which has limited ability to detect minor variants. Next generation sequencing (NGS) platforms enable variant detection at frequencies as low as 1% allowing for earlier detection of resistance and modification of therapy. Implementation of NGS assays in the clinical laboratory is hindered by complicated assay design, cumbersome wet bench procedures, and the complexity of data analysis and bioinformatics. We developed a complete NGS protocol and companion analysis and reporting pipeline using AmpliSeq multiplex PCR, Ion Torrent S5 XL sequencing, and Stanford's HIVdb resistance algorithm. Implemented as a Torrent Suite software plugin, the pipeline runs automatically after sequencing. An optimum variant frequency threshold of 10% was determined by comparing Sanger sequences of archived samples from ViroSeq testing, resulting in a sensitivity of 98.2% and specificity of 99.0%. The majority (91%) of drug resistance mutations were detected by both Sanger and NGS, with 1.7% only by Sanger and 7.3% only by NGS. Variant calls were highly reproducible and there was no cross-reactivity to VZV, HBV, CMV, EBV, and HCV. The limit of detection was 500 copies/mL. The NGS assay performance was comparable to ViroSeq Sanger sequencing and has several advantages, including a publicly available end-to-end analysis and reporting plugin. The assay provides a straightforward path for implementation of NGS for HIV drug resistance testing in the laboratory setting without additional investment in bioinformatics infrastructure and resources.

Metabolic and functional impairment of CD8+ T cells from the lungs of influenza-infected obese mice

Obesity is an independent risk factor for morbidity and mortality in response to influenza infection. However, the underlying mechanisms by which obesity impairs immunity are unclear. Herein, we investigated the effects of diet-induced obesity on pulmonary CD8+ T cell metabolism, cytokine production, and transcriptome as a potential mechanism of impairment during influenza virus infection in mice. Male C57BL/6J lean and obese mice were infected with sub-lethal mouse-adapted A/PR/8/34 influenza virus, generating a pulmonary anti-viral and inflammatory response. Extracellular metabolic flux analyses revealed pulmonary CD8+ T cells from obese mice, compared with lean controls, had suppressed oxidative and glycolytic metabolism at day 10 post-infection. Flow cytometry showed the impairment in pulmonary CD8+ T cell metabolism with obesity was independent of changes in glucose or fatty acid uptake, but concomitant with decreased CD8+ GrB+ IFNγ+ populations. Notably, the percent of pulmonary effector CD8+ GrB+ IFNγ+ T cells at day 10 post-infection correlated positively with total CD8+ basal extracellular acidification rate and basal oxygen consumption rate. Finally, next-generation RNA sequencing revealed complex and unique transcriptional regulation of sorted effector pulmonary CD8+ CD44+ T cells from obese mice compared to lean mice following influenza infection. Collectively, the data suggest diet-induced obesity increases influenza virus pathogenesis, in part, through CD8+ T cell-mediated metabolic reprogramming and impaired effector CD8+ T cell function.

Sensitive and accurate analysis of gene expression signatures enabled by oligonucleotide-labelled cDNA

High-throughput RNA sequencing offers a comprehensive analysis of transcriptome complexity originated from regulatory events, such as differential gene expression, alternative polyadenylation and others, and allows the increase in diagnostic capacity and precision. For gene expression profiling applications that do not specifically require information on alternative splicing events, the mRNA 3' termini counting approach is a cost-effective alternative to whole transcriptome sequencing. Here, we report MTAS-seq (mRNA sequencing via terminator-assisted synthesis) - a novel RNA-seq library preparation method directed towards mRNA 3' termini. We demonstrate the specific enrichment for 3'-terminal regions by simple and quick single-tube protocol with built-in molecular barcoding to enable accurate estimation of transcript abundance. To achieve that, we synthesized oligonucleotide-modified dideoxynucleotides which enable the generation of cDNA libraries at the reverse transcription step. We validated the performance of MTAS-seq on well-characterized reference bulk RNA and further tested it with eukaryotic cell lysates.

Independent variations in genome-wide expression, alternative splicing, and DNA methylation in brain tissues among castes of the buff-tailed bumblebee, Bombus terrestris

Caste differentiation in social hymenopterans is an intriguing example of phenotypic plasticity. However, the co-ordination among gene regulatory factors to mediate caste differentiation remains inconclusive. In this study, we determined the role of gene regulation and related epigenetic processes in pre-imaginal caste differentiation in the primitively eusocial bumblebee Bombus terrestris. By combining RNA-Seq data from Illumina and PacBio and accurately quantifying methylation at whole-genomic base pair resolution, we found that queens, workers, and drones mainly differentiate in gene expression but not in alternative splicing and DNA methylation. Gynes are the most distinct with the lowest global level of whole-genomic methylation and with the largest number of caste-specific transcripts and alternative splicing events. By contrast, workers exhibit few uniquely expressed or alternatively spliced genes. Moreover, several genes involved in hormone and neurotransmitter metabolism are related to caste differentiation, whereas several neuropeptides are linked with sex differentiation. Despite little genome-wide association among differential gene expression, splicing, and differential DNA methylation, the overlapped gene ontology (GO) terms point to nutrition-related activity. Therefore, variations in gene regulation correlate with the behavioral differences among castes and highlight the specialization of toolkit genes in bumblebee gynes at the beginning of the adult stage.

Transcriptional Profiling of CD8+ CMV-Specific T Cell Functional Subsets Obtained Using a Modified Method for Isolating High-Quality RNA From Fixed and Permeabilized Cells

Previous studies suggest that the presence of antigen-specific polyfunctional T cells is correlated with improved pathogen clearance, disease control, and clinical outcomes; however, the molecular mechanisms responsible for the generation, function, and survival of polyfunctional T cells remain unknown. The study of polyfunctional T cells has been, in part, limited by the need for intracellular cytokine staining (ICS), necessitating fixation and cell membrane permeabilization that leads to unacceptable degradation of RNA. Adopting elements from prior research efforts, we developed and optimized a modified protocol for the isolation of high-quality RNA (i.e., RIN > 7) from primary human T cells following aldehyde-fixation, detergent-based permeabilization, intracellular cytokines staining, and sorting. Additionally, this method also demonstrated utility preserving RNA when staining for transcription factors. This modified protocol utilizes an optimized combination of an RNase inhibitor and high-salt buffer that is cost-effective while maintaining the ability to identify and resolve cell populations for sorting. Overall, this protocol resulted in minimal loss of RNA integrity, quality, and quantity during cytoplasmic staining of cytokines and subsequent flourescence-activated cell sorting. Using this technique, we obtained the transcriptional profiles of functional subsets (i.e., non-functional, monofunctional, bifunctional, polyfunctional) of CMV-specific CD8+T cells. Our analyses demonstrated that these functional subsets are molecularly distinct, and that polyfunctional T cells are uniquely enriched for transcripts involved in viral response, inflammation, cell survival, proliferation, and metabolism when compared to monofunctional cells. Polyfunctional T cells demonstrate reduced activation-induced cell death and increased proliferation after antigen re-challenge. Further in silico analysis of transcriptional data suggested a critical role for STAT5 transcriptional activity in polyfunctional cell activation. Pharmacologic inhibition of STAT5 was associated with a significant reduction in polyfunctional cell cytokine expression and proliferation, demonstrating the requirement of STAT5 activity not only for proliferation and cell survival, but also cytokine expression. Finally, we confirmed this association between CMV-specific CD8+ polyfunctionality with STAT5 signaling also exists in immunosuppressed transplant recipients using single cell transcriptomics, indicating that results from this study may translate to this vulnerable patient population. Collectively, these results shed light on the mechanisms governing polyfunctional T cell function and survival and may ultimately inform multiple areas of immunology, including but not limited to the development of new vaccines, CAR-T cell therapies, and adoptive T cell transfer.

Advances in Genomics for Drug Development

Drug development (target identification, advancing drug leads to candidates for preclinical and clinical studies) can be facilitated by genetic and genomic knowledge. Here, we review the contribution of population genomics to target identification, the value of bulk and single cell gene expression analysis for understanding the biological relevance of a drug target, and genome-wide CRISPR editing for the prioritization of drug targets. In genomics, we discuss the different scope of genome-wide association studies using genotyping arrays, versus exome and whole genome sequencing. In transcriptomics, we discuss the information from drug perturbation and the selection of biomarkers. For CRISPR screens, we discuss target discovery, mechanism of action and the concept of gene to drug mapping. Harnessing genetic support increases the probability of drug developability and approval.

Functional Genomics of Healthy and Pathological Fetal Membranes

Premature preterm rupture of membranes (PPROM), rupture of fetal membranes before 37 weeks of gestation, is the leading identifiable cause of spontaneous preterm births. Often there is no obvious cause that is identified in a patient who presents with PPROM. Identifying the upstream molecular events that lead to fetal membrane weakening presents potentially actionable mechanisms which could lead to the identification of at-risk patients and to the development of new therapeutic interventions. Functional genomic studies have transformed understanding of the role of gene regulation in diverse cells and tissues involved health and disease. Here, we review the results of those studies in the context of fetal membranes. We will highlight relevant results from major coordinated functional genomics efforts and from targeted studies focused on individual cell or tissue models. Studies comparing gene expression and DNA methylation between healthy and pathological fetal membranes have found differential regulation between labor and quiescent tissue as well as in preterm births, preeclampsia, and recurrent pregnancy loss. Whole genome and exome sequencing studies have identified common and rare fetal variants associated with preterm births. However, few fetal membrane tissue studies have modeled the response to stimuli relevant to pregnancy. Fetal membranes are readily adaptable to cell culture and relevant cellular phenotypes are readily observable. For these reasons, this is now an unrealized opportunity for genomic studies isolating the effect of cell signaling cascades and mapping the fetal membrane responses that lead to PPROM and other pregnancy complications.