A metric for evaluating biological information in gene sets and its application to identify co-expressed gene clusters in PBMC.

Recent technological advances have made the gathering of comprehensive gene expression datasets a commodity. This has shifted the limiting step of transcriptomic studies from the accumulation of data to their analyses and interpretation. The main problem in analyzing transcriptomics data is that the number of independent samples is typically much lower (<100) than the number of genes whose expression is quantified (typically >14,000). To address this, it would be desirable to reduce the gathered data's dimensionality without losing information. Clustering genes into discrete modules is one of the most commonly used tools to accomplish this task. While there are multiple clustering approaches, there is a lack of informative metrics available to evaluate the resultant clusters' biological quality. Here we present a metric that incorporates known ground truth gene sets to quantify gene clusters' biological quality derived from standard clustering techniques. The GECO (Ground truth Evaluation of Clustering Outcomes) metric demonstrates that quantitative and repeatable scoring of gene clusters is not only possible but computationally lightweight and robust. Unlike current methods, it allows direct comparison between gene clusters generated by different clustering techniques. It also reveals that current cluster analysis techniques often underestimate the number of clusters that should be formed from a dataset, which leads to fewer clusters of lower quality. As a test case, we applied GECO combined with k-means clustering to derive an optimal set of co-expressed gene modules derived from PBMC, which we show to be superior to previously generated modules generated on whole-blood. Overall, GECO provides a rational metric to test and compare different clustering approaches to analyze high-dimensional transcriptomic data.

HLA-DR Marks Recently Divided Antigen-Specific Effector CD4 T Cells in Active Tuberculosis Patients.

Upon Ag encounter, T cells can rapidly divide and form an effector population, which plays an important role in fighting acute infections. In humans, little is known about the molecular markers that distinguish such effector cells from other T cell populations. To address this, we investigated the molecular profile of T cells present in individuals with active tuberculosis (ATB), where we expect Ag encounter and expansion of effector cells to occur at higher frequency in contrast to Mycobacterium tuberculosis-sensitized healthy IGRA+ individuals. We found that the frequency of HLA-DR+ cells was increased in circulating CD4 T cells of ATB patients, and was dominantly expressed in M. tuberculosis Ag-specific CD4 T cells. We tested and confirmed that HLA-DR is a marker of recently divided CD4 T cells upon M. tuberculosis Ag exposure using an in vitro model examining the response of resting memory T cells from healthy IGRA+ to Ags. Thus, HLA-DR marks a CD4 T cell population that can be directly detected ex vivo in human peripheral blood, whose frequency is increased during ATB disease and contains recently divided Ag-specific effector T cells. These findings will facilitate the monitoring and study of disease-specific effector T cell responses in the context of ATB and other infections.

A system-view of Bordetella pertussis booster vaccine responses in adults primed with whole-cell versus acellular vaccine in infancy.

The increased incidence of whooping cough worldwide suggests that current vaccination against Bordetella pertussis infection has limitations in quality and duration of protection. The resurgence of infection has been linked to the introduction of acellular vaccines (aP), which have an improved safety profile compared with the previously used whole-cell (wP) vaccines. To determine immunological differences between aP and wP priming in infancy, we performed a systems approach of the immune response to booster vaccination. Transcriptomic, proteomic, cytometric, and serologic profiling revealed multiple shared immune responses with different kinetics across cohorts, including an increase of blood monocyte frequencies and strong antigen-specific IgG responses. Additionally, we found a prominent subset of aP-primed individuals (30%) with a strong differential signature, including higher levels of expression for CCL3, NFKBIA, and ICAM1. Contrary to the wP individuals, this subset displayed increased PT-specific IgE responses after boost and higher antigen-specific IgG4 and IgG3 antibodies against FHA and FIM2/3 at baseline and after boost. Overall, the results show that, while broad immune response patterns to Tdap boost overlap between aP- and wP-primed individuals, a subset of aP-primed individuals present a divergent response. These findings provide candidate targets to study the causes and correlates of waning immunity after aP vaccination.

Quantitative and Qualitative Perturbations of CD8+ MAITs in Healthy Mycobacterium tuberculosis-Infected Individuals.

CD8 T cells are considered important contributors to the immune response against Mycobacterium tuberculosis, yet limited information is currently known regarding their specific immune signature and phenotype. In this study, we applied a cell population transcriptomics strategy to define immune signatures of human latent tuberculosis infection (LTBI) in memory CD8 T cells. We found a 41-gene signature that discriminates between memory CD8 T cells from healthy LTBI subjects and uninfected controls. The gene signature was dominated by genes associated with mucosal-associated invariant T cells (MAITs) and reflected the lower frequency of MAITs observed in individuals with LTBI. There was no evidence for a conventional CD8 T cell-specific signature between the two cohorts. We, therefore, investigated MAITs in more detail based on Vα7.2 and CD161 expression and staining with an MHC-related protein 1 (MR1) tetramer. This revealed two distinct populations of CD8+Vα7.2+CD161+ MAITs: MR1 tetramer+ and MR1 tetramer-, which both had distinct gene expression compared with memory CD8 T cells. Transcriptomic analysis of LTBI versus noninfected individuals did not reveal significant differences for MR1 tetramer+ MAITs. However, gene expression of MR1 tetramer- MAITs showed large interindividual diversity and a tuberculosis-specific signature. This was further strengthened by a more diverse TCR-α and -ß repertoire of MR1 tetramer- cells as compared with MR1 tetramer+ Thus, circulating memory CD8 T cells in subjects with latent tuberculosis have a reduced number of conventional MR1 tetramer+ MAITs as well as a difference in phenotype in the rare population of MR1 tetramer- MAITs compared with uninfected controls.

The Challenge of Distinguishing Cell-Cell Complexes from Singlet Cells in Non-Imaging Flow Cytometry and Single-Cell Sorting.

Our recent work has highlighted that care needs to be taken when interpreting single cell data originating from flow cytometry acquisition or cell sorting: We found that doublets of T cells bound to other immune cells are often present in the live singlet gate of human peripheral blood samples acquired by flow cytometry. This hidden "contamination" generates atypical gene signatures of mixed cell lineage in what is assumed to be single cells, which can lead to data misinterpretation, such as the description of novel immune cell types. Here, based on the example of T cell-monocyte complexes, we identify experimental and data analysis strategies to help distinguishing between singlets and cell-cell complexes in non-imaging flow cytometry and single-cell sorting. We found robust molecular signatures in both T cell-monocyte and T cell-B cell complexes that can distinguish them from singlets at both protein and mRNA levels. Imaging flow cytometry with appropriate gating strategy (matching the one used in cell sorting) and direct microscopy imaging after cell sorting were the two methods of choice to detect the presence of cell-cell complexes in suspicious dual-expressing cells. We finally applied this knowledge to highlight the likely presence of T cell-B cell complexes in a recently published dataset describing a novel cell population with mixed T cell and B cell lineage properties. © 2020 The Authors. Cytometry Part A published by Wiley Periodicals, Inc. on behalf of International Society for Advancement of Cytometry.

Identifying inaccuracies in gene expression estimates from unstranded RNA-seq data.

RNA-seq methods are widely utilized for transcriptomic profiling of biological samples. However, there are known caveats of this technology which can skew the gene expression estimates. Specifically, if the library preparation protocol does not retain RNA strand information then some genes can be erroneously quantitated. Although strand-specific protocols have been established, a significant portion of RNA-seq data is generated in non-strand-specific manner. We used a comprehensive stranded RNA-seq dataset of 15 blood cell types to identify genes for which expression would be erroneously estimated if strand information was not available. We found that about 10% of all genes and 2.5% of protein coding genes have a two-fold or higher difference in estimated expression when strand information of the reads was ignored. We used parameters of read alignments of these genes to construct a machine learning model that can identify which genes in an unstranded dataset might have incorrect expression estimates and which ones do not. We also show that differential expression analysis of genes with biased expression estimates in unstranded read data can be recovered by limiting the reads considered to those which span exonic boundaries. The resulting approach is implemented as a package available at https://github.com/mikpom/uslcount .

Circulating T cell-monocyte complexes are markers of immune perturbations.

Our results highlight for the first time that a significant proportion of cell doublets in flow cytometry, previously believed to be the result of technical artifacts and thus ignored in data acquisition and analysis, are the result of biological interaction between immune cells. In particular, we show that cell:cell doublets pairing a T cell and a monocyte can be directly isolated from human blood, and high resolution microscopy shows polarized distribution of LFA1/ICAM1 in many doublets, suggesting in vivo formation. Intriguingly, T cell-monocyte complex frequency and phenotype fluctuate with the onset of immune perturbations such as infection or immunization, reflecting expected polarization of immune responses. Overall these data suggest that cell doublets reflecting T cell-monocyte in vivo immune interactions can be detected in human blood and that the common approach in flow cytometry to avoid studying cell:cell complexes should be re-visited.

Host Transcriptomics as a Tool to Identify Diagnostic and Mechanistic Immune Signatures of Tuberculosis.

Tuberculosis (TB) is a major infectious disease worldwide, and is associated with several challenges for control and eradication. First, more accurate diagnostic tools that better represent the spectrum of infection states are required; in particular, identify the latent TB infected individuals with high risk of developing active TB. Second, we need to better understand, from a mechanistic point of view, why the immune system is unsuccessful in some cases for control and elimination of the pathogen. Host transcriptomics is a powerful approach to identify both diagnostic and mechanistic immune signatures of diseases. We have recently reported that optimal study design for these two purposes should be guided by different sets of criteria. Here, based on already published transcriptomics signatures of tuberculosis, we further develop these guidelines and identify additional factors to consider for obtaining diagnostic vs. mechanistic signatures in terms of cohorts, samples, data generation and analysis. Diagnostic studies should aim to identify small disease signatures with high discriminatory power across all affected populations, and against similar pathologies to TB. Specific focus should be made on improving the diagnosis of infected individuals at risk of developing active disease. Conversely, mechanistic studies should focus on tissues biopsies, immune relevant cell subsets, state of the art transcriptomic techniques and bioinformatics tools to understand the biological meaning of identified gene signatures that could facilitate therapeutic interventions. Finally, investigators should ensure their data are made publicly available along with complete annotations to facilitate metadata and cross-study analyses.

GOnet: a tool for interactive Gene Ontology analysis.

BACKGROUND: Biological interpretation of gene/protein lists resulting from -omics experiments can be a complex task. A common approach consists of reviewing Gene Ontology (GO) annotations for entries in such lists and searching for enrichment patterns. Unfortunately, there is a gap between machine-readable output of GO software and its human-interpretable form. This gap can be bridged by allowing users to simultaneously visualize and interact with term-term and gene-term relationships. RESULTS: We created the open-source GOnet web-application (available at http://tools.dice-database.org/GOnet/ ), which takes a list of gene or protein entries from human or mouse data and performs GO term annotation analysis (mapping of provided entries to GO subsets) or GO term enrichment analysis (scanning for GO categories overrepresented in the input list). The application is capable of producing parsable data formats and importantly, interactive visualizations of the GO analysis results. The interactive results allow exploration of genes and GO terms as a graph that depicts the natural hierarchy of the terms and retains relationships between terms and genes/proteins. As a result, GOnet provides insight into the functional interconnection of the submitted entries. CONCLUSIONS: The application can be used for GO analysis of any biological data sources resulting in gene/protein lists. It can be helpful for experimentalists as well as computational biologists working on biological interpretation of -omics data resulting in such lists.

Whole transcriptome profiling of adult and infective stages of the trematode Opisthorchis felineus.

Opisthorchis felineus, the trematode belonging to the family Opisthorchiidae, is a causative agent of the infection called opisthorchiasis or liver fluke infection. Being a close relative of Opisthorchis viverrini and Clonorchis sinensis (oriental liver flukes) it is encountered in northern Eurasia, especially in Russia, Kazakhstan, Belarus, Ukraine, and Baltic countries. Whole genome data for oriental liver flukes revealed their adaptations for life in the bile duct but our knowledge of O. felineus is scarce. To address this knowledge gap and uncover evolutionary aspect of the adaptations on the transcriptomic level, we used RNA-sequencing approach to investigate two stages of the parasite residing in different hosts. Bioinformatic analysis revealed specific features affecting various biochemical pathways and gene networks. Namely, we observed the loss of genes involved in polyamine synthesis, methionine salvage and peroxisome biogenesis. Some of the gene families, like MD-2 lipid binding proteins, calmodulins and cathepsins on the contrary have expanded compared to free living eukaryotes. We identified significant differences between the stages in homeodomain-containing genes, G-protein coupled receptors, and neuroactive signaling systems. Granulin-like growth factors specific for O. felineus were also identified. In this work, we provide the first whole transcriptome investigation of this parasite. We also hope that these results will create a background for further molecular research of helminth infections and opisthorchiasis in particular.

Comparison of the three-dimensional organization of sperm and fibroblast genomes using the Hi-C approach.

BACKGROUND: The three-dimensional organization of the genome is tightly connected to its biological function. The Hi-C approach was recently introduced as a method that can be used to identify higher-order chromatin interactions genome-wide. The aim of this study was to determine genome-wide chromatin interaction frequencies using the Hi-C approach in mouse sperm cells and embryonic fibroblasts. RESULTS: The obtained data demonstrate that the three-dimensional genome organizations of sperm and fibroblast cells show a high degree of similarity both with each other and with the previously described mouse embryonic stem cells. Both A- and B-compartments and topologically associated domains are present in spermatozoa and fibroblasts. Nevertheless, sperm cells and fibroblasts exhibit statistically significant differences between each other in the contact probabilities of defined loci. Tight packaging of the sperm genome results in an enrichment of long-range contacts compared with the fibroblasts. However, only 30% of the differences in the number of contacts are based on differences in the densities of their genome packages; the main source of the differences is the gain or loss of contacts that are specific for defined genome regions. We find that the dependence of the contact probability on genomic distance for sperm is close to the dependence predicted for the fractal globular folding of chromatin. CONCLUSIONS: Overall, we can conclude that the three-dimensional structure of the genome is passed through generations without being dramatically changed in sperm cells.

Adult Opisthorchis felineus major protein fractions deduced from transcripts: comparison with liver flukes Opisthorchis viverrini and Clonorchis sinensis.

The epidemiologically important liver flukes Opisthorchis felineus, Opisthorchis viverrini, and Clonorchis sinensis are of interest to health professionals, epidemiologists, pharmacologists, and molecular biologists. Recently the transcriptomes of the latter two species were intensively investigated. However our knowledge on molecular biology of O. felineus is scarce. We report the first results of the O. felineus transcriptome analysis. We isolated and annotated a total of 2560 expressed sequence tag (EST) sequences from adult O. felineus (deposited within the database of expressed sequence tags (dbEST), under accession numbers GenBank: JK624271-JK626790, JK006511-JK006547, JK649790-JK649792). Clustering and analysis resulted in the detection of 267 contigs. Of the protein sequences deduced from these, 82% had homologs in the NCBI (nr) protein database and 63% contained conserved domains, allowing the functions to be interpreted using the Gene Ontology terms. Comprehensive analysis of Opisthorchiidae- and Trematoda-specific substitutions within amino acid sequences deduced for the proteins myoglobin, vitelline precursor protein, cathepsin F, and 28kDa glutathione transferase was carried out. The gene set of the 32 ribosomal proteins for the three Opisthorchiidae species with the addition of available Schistosoma and Fasciola orthologs was created and is provided in the supplementary. The orthologous gene set created was used for inferring phylogeny within the Trematoda with special attention to interrelations within the Opisthorchiidae. The phylogenetic analysis revealed a closer relationship between C. sinensis and O. viverrini and some divergence of O. felineus from either O. viverrini or C. sinensis.