Natural variation in the parameters of innate immune cells is preferentially driven by genetic factors.

The quantification and characterization of circulating immune cells provide key indicators of human health and disease. To identify the relative effects of environmental and genetic factors on variation in the parameters of innate and adaptive immune cells in homeostatic conditions, we combined standardized flow cytometry of blood leukocytes and genome-wide DNA genotyping of 1,000 healthy, unrelated people of Western European ancestry. We found that smoking, together with age, sex and latent infection with cytomegalovirus, were the main non-genetic factors that affected variation in parameters of human immune cells. Genome-wide association studies of 166 immunophenotypes identified 15 loci that showed enrichment for disease-associated variants. Finally, we demonstrated that the parameters of innate cells were more strongly controlled by genetic variation than were those of adaptive cells, which were driven by mainly environmental exposure. Our data establish a resource that will generate new hypotheses in immunology and highlight the role of innate immunity in susceptibility to common autoimmune diseases.

Publisher Correction: Natural variation in the parameters of innate immune cells is preferentially driven by genetic factors.

In the version of this article initially published, the name of one author was incorrect (James P. Santo). The correct name is James P. Di Santo. The error has been corrected in the HTML and PDF versions of the article.

Polygenic risk for skin autoimmunity impacts immune checkpoint blockade in bladder cancer.

PD-1 and PD-L1 act to restrict T cell responses in cancer and contribute to self-tolerance. Consistent with this role, PD-1 checkpoint inhibitors have been associated with immune-related adverse events (irAEs), immune toxicities thought to be autoimmune in origin. Analyses of dermatological irAEs have identified an association with improved overall survival (OS) following anti-PD-(L)1 therapy, but the factors that contribute to this relationship are poorly understood. We collected germline whole-genome sequencing data from IMvigor211, a recent phase 3 randomized controlled trial comparing atezolizumab (anti-PD-L1) monotherapy to chemotherapy in bladder cancer. We found that high vitiligo, high psoriasis, and low atopic dermatitis polygenic risk scores (PRSs) were associated with longer OS under anti-PD-L1 monotherapy as compared to chemotherapy, reflecting the Th17 polarization of these diseases. PRSs were not correlated with tumor mutation burden, PD-L1 immunohistochemistry, nor T-effector gene signatures. Shared genetic factors impact risk for dermatological autoimmunity and anti-PD-L1 monotherapy in bladder cancer.

SMSSVD: SubMatrix Selection Singular Value Decomposition.

Motivation: High throughput biomedical measurements normally capture multiple overlaid biologically relevant signals and often also signals representing different types of technical artefacts like e.g. batch effects. Signal identification and decomposition are accordingly main objectives in statistical biomedical modeling and data analysis. Existing methods, aimed at signal reconstruction and deconvolution, in general, are either supervised, contain parameters that need to be estimated or present other types of ad hoc features. We here introduce SubMatrix Selection Singular Value Decomposition (SMSSVD), a parameter-free unsupervised signal decomposition and dimension reduction method, designed to reduce noise, adaptively for each low-rank-signal in a given data matrix, and represent the signals in the data in a way that enable unbiased exploratory analysis and reconstruction of multiple overlaid signals, including identifying groups of variables that drive different signals. Results: The SMSSVD method produces a denoised signal decomposition from a given data matrix. It also guarantees orthogonality between signal components in a straightforward manner and it is designed to make automation possible. We illustrate SMSSVD by applying it to several real and synthetic datasets and compare its performance to golden standard methods like PCA (Principal Component Analysis) and SPC (Sparse Principal Components, using Lasso constraints). The SMSSVD is computationally efficient and despite being a parameter-free method, in general, outperforms existing statistical learning methods. Availability and implementation: A Julia implementation of SMSSVD is openly available on GitHub (https://github.com/rasmushenningsson/SubMatrixSelectionSVD.jl). Supplementary information: Supplementary data are available at Bioinformatics online.

Increasing Clinical Severity during a Dengue Virus Type 3 Cuban Epidemic: Deep Sequencing of Evolving Viral Populations.

UNLABELLED: During the dengue virus type 3 (DENV-3) epidemic that occurred in Havana in 2001 to 2002, severe disease was associated with the infection sequence DENV-1 followed by DENV-3 (DENV-1/DENV-3), while the sequence DENV-2/DENV-3 was associated with mild/asymptomatic infections. To determine the role of the virus in the increasing severity demonstrated during the epidemic, serum samples collected at different time points were studied. A total of 22 full-length sequences were obtained using a deep-sequencing approach. Bayesian phylogenetic analysis of consensus sequences revealed that two DENV-3 lineages were circulating in Havana at that time, both grouped within genotype III. The predominant lineage is closely related to Peruvian and Ecuadorian strains, while the minor lineage is related to Venezuelan strains. According to consensus sequences, relatively few nonsynonymous mutations were observed; only one was fixed during the epidemic at position 4380 in the NS2B gene. Intrahost genetic analysis indicated that a significant minor population was selected and became predominant toward the end of the epidemic. In conclusion, greater variability was detected during the epidemic's progression in terms of significant minority variants, particularly in the nonstructural genes. An increasing trend of genetic diversity toward the end of the epidemic was observed only for synonymous variant allele rates, with higher variability in secondary cases. Remarkably, significant intrahost genetic variation was demonstrated within the same patient during the course of secondary infection with DENV-1/DENV-3, including changes in the structural proteins premembrane (PrM) and envelope (E). Therefore, the dynamic of evolving viral populations in the context of heterotypic antibodies could be related to the increasing clinical severity observed during the epidemic. IMPORTANCE: Based on the evidence that DENV fitness is context dependent, our research has focused on the study of viral factors associated with intraepidemic increasing severity in a unique epidemiological setting. Here, we investigated the intrahost genetic diversity in acute human samples collected at different time points during the DENV-3 epidemic that occurred in Cuba in 2001 to 2002 using a deep-sequencing approach. We concluded that greater variability in significant minor populations occurred as the epidemic progressed, particularly in the nonstructural genes, with higher variability observed in secondary infection cases. Remarkably, for the first time significant intrahost genetic variation was demonstrated within the same patient during the course of secondary infection with DENV-1/DENV-3, including changes in structural proteins. These findings indicate that high-resolution approaches are needed to unravel molecular mechanisms involved in dengue pathogenesis.

Group Selection and Contribution of Minority Variants during Virus Adaptation Determines Virus Fitness and Phenotype.

Understanding how a pathogen colonizes and adapts to a new host environment is a primary aim in studying emerging infectious diseases. Adaptive mutations arise among the thousands of variants generated during RNA virus infection, and identifying these variants will shed light onto how changes in tropism and species jumps can occur. Here, we adapted Coxsackie virus B3 to a highly permissive and less permissive environment. Using deep sequencing and bioinformatics, we identified a multi-step adaptive process to adaptation involving residues in the receptor footprints that correlated with receptor availability and with increase in virus fitness in an environment-specific manner. We show that adaptation occurs by selection of a dominant mutation followed by group selection of minority variants that together, confer the fitness increase observed in the population, rather than selection of a single dominant genotype.

What is a "unimodal" cell population? Using statistical tests as criteria for unimodality in automated gating and quality control.

Many automated gating algorithms for flow cytometry data are based on the concept of unimodal cell populations. However, in this article, we show that criteria previously used to make decisions on unimodality cannot adequately distinguish unimodal from bimodal densities. We show that dip and bandwidth tests for unimodality, taken from the statistics literature, can do this with consistent and low error rates. These tests also have the possibility to adjust the significance level to handle the trade-off between failing to detect a second mode and seeing a second mode when there is none. The differences between the dip and bandwidth tests are elucidated using real data from the FlowCAP I challenge, also guidelines for flow cytometry data preprocessing are given. © 2017 International Society for Advancement of Cytometry.

BayesFlow: latent modeling of flow cytometry cell populations.

BACKGROUND: Flow cytometry is a widespread single-cell measurement technology with a multitude of clinical and research applications. Interpretation of flow cytometry data is hard; the instrumentation is delicate and can not render absolute measurements, hence samples can only be interpreted in relation to each other while at the same time comparisons are confounded by inter-sample variation. Despite this, most automated flow cytometry data analysis methods either treat samples individually or ignore the variation by for example pooling the data. A key requirement for models that include multiple samples is the ability to visualize and assess inferred variation, since what could be technical variation in one setting would be different phenotypes in another. RESULTS: We introduce BayesFlow, a pipeline for latent modeling of flow cytometry cell populations built upon a Bayesian hierarchical model. The model systematizes variation in location as well as shape. Expert knowledge can be incorporated through informative priors and the results can be supervised through compact and comprehensive visualizations. BayesFlow is applied to two synthetic and two real flow cytometry data sets. For the first real data set, taken from the FlowCAP I challenge, BayesFlow does not only give a gating which would place it among the top performers in FlowCAP I for this dataset, it also gives a more consistent treatment of different samples than either manual gating or other automated gating methods. The second real data set contains replicated flow cytometry measurements of samples from healthy individuals. BayesFlow gives here cell populations with clear expression patterns and small technical intra-donor variation as compared to biological inter-donor variation. CONCLUSIONS: Modeling latent relations between samples through BayesFlow enables a systematic analysis of inter-sample variation. As opposed to other joint gating methods, effort is put at ensuring that the obtained partition of the data corresponds to actual cell populations, and the result is therefore directly biologically interpretable. BayesFlow is freely available at GitHub.

The ABCs of viral hepatitis that define biomarker signatures of acute viral hepatitis.

UNLABELLED: Viral hepatitis is the leading cause of liver disease worldwide and can be caused by several agents, including hepatitis A (HAV), B (HBV), and C (HCV) virus. We employed multiplexed protein immune assays to identify biomarker signatures of viral hepatitis in order to define unique and common responses for three different acute viral infections of the liver. We performed multianalyte profiling, measuring the concentrations of 182 serum proteins obtained from acute HAV- (18), HBV- (18), and HCV-infected (28) individuals, recruited as part of a hospital-based surveillance program in Cairo, Egypt. Virus-specific biomarker signatures were identified and validation was performed using a unique patient population. A core signature of 46 plasma proteins was commonly modulated in all three infections, as compared to healthy controls. Principle component analysis (PCA) revealed a host response based upon 34 proteins, which could distinguish HCV patients from HAV- and HBV-infected individuals or healthy controls. When HAV and HBV groups were compared directly, 34 differentially expressed serum proteins allowed the separation of these two patient groups. A validation study was performed on an additional 111 patients, confirming the relevance of our initial findings, and defining the 17 analytes that reproducibly segregated the patient populations. CONCLUSIONS: This combined discovery and biomarker validation approach revealed a previously unrecognized virus-specific induction of host proteins. The identification of hepatitis virus specific signatures provides a foundation for functional studies and the identification of potential correlates of viral clearance.

A framework for list representation, enabling list stabilization through incorporation of gene exchangeabilities.

Analysis of multivariate data sets from, for example, microarray studies frequently results in lists of genes which are associated with some response of interest. The biological interpretation is often complicated by the statistical instability of the obtained gene lists, which may partly be due to the functional redundancy among genes, implying that multiple genes can play exchangeable roles in the cell. In this paper, we use the concept of exchangeability of random variables to model this functional redundancy and thereby account for the instability. We present a flexible framework to incorporate the exchangeability into the representation of lists. The proposed framework supports straightforward comparison between any 2 lists. It can also be used to generate new more stable gene rankings incorporating more information from the experimental data. Using 2 microarray data sets, we show that the proposed method provides more robust gene rankings than existing methods with respect to sampling variations, without compromising the biological significance of the rankings.

The correlation pattern of acquired copy number changes in 164 ETV6/RUNX1-positive childhood acute lymphoblastic leukemias.

The ETV6/RUNX1 fusion gene, present in 25% of B-lineage childhood acute lymphoblastic leukemia (ALL), is thought to represent an initiating event, which requires additional genetic changes for leukemia development. To identify additional genetic alterations, 24 ETV6/RUNX1-positive ALLs were analyzed using 500K single nucleotide polymorphism arrays. The results were combined with previously published data sets, allowing us to ascertain genomic copy number aberrations (CNAs) in 164 cases. In total, 45 recurrent CNAs were identified with an average number of 3.5 recurrent changes per case (range 0-13). Twenty-six percent of cases displayed a set of recurrent CNAs identical to that of other cases in the data set. The majority (74%), however, displayed a unique pattern of recurrent CNAs, indicating a large heterogeneity within this ALL subtype. As previously demonstrated, alterations targeting genes involved in B-cell development were common (present in 28% of cases). However, the combined analysis also identified alterations affecting nuclear hormone response (24%) to be a characteristic feature of ETV6/RUNX1-positive ALL. Studying the correlation pattern of the CNAs allowed us to highlight significant positive and negative correlations between specific aberrations. Furthermore, oncogenetic tree models identified ETV6, CDKN2A/B, PAX5, del(6q) and +16 as possible early events in the leukemogenic process.

The projection score--an evaluation criterion for variable subset selection in PCA visualization.

BACKGROUND: In many scientific domains, it is becoming increasingly common to collect high-dimensional data sets, often with an exploratory aim, to generate new and relevant hypotheses. The exploratory perspective often makes statistically guided visualization methods, such as Principal Component Analysis (PCA), the methods of choice. However, the clarity of the obtained visualizations, and thereby the potential to use them to formulate relevant hypotheses, may be confounded by the presence of the many non-informative variables. For microarray data, more easily interpretable visualizations are often obtained by filtering the variable set, for example by removing the variables with the smallest variances or by only including the variables most highly related to a specific response. The resulting visualization may depend heavily on the inclusion criterion, that is, effectively the number of retained variables. To our knowledge, there exists no objective method for determining the optimal inclusion criterion in the context of visualization. RESULTS: We present the projection score, which is a straightforward, intuitively appealing measure of the informativeness of a variable subset with respect to PCA visualization. This measure can be universally applied to find suitable inclusion criteria for any type of variable filtering. We apply the presented measure to find optimal variable subsets for different filtering methods in both microarray data sets and synthetic data sets. We note also that the projection score can be applied in general contexts, to compare the informativeness of any variable subsets with respect to visualization by PCA. CONCLUSIONS: We conclude that the projection score provides an easily interpretable and universally applicable measure of the informativeness of a variable subset with respect to visualization by PCA, that can be used to systematically find the most interpretable PCA visualization in practical exploratory analysis.

Transcriptional regulation of aquaporins in accessions of Arabidopsis in response to drought stress.

Aquaporins facilitate water transport over cellular membranes, and are therefore believed to play an important role in water homeostasis. In higher plants aquaporin-like proteins, also called major intrinsic proteins (MIPs), are divided into five subfamilies. We have previously shown that MIP transcription in Arabidopsis thaliana is generally downregulated in leaves upon drought stress, apart from two members of the plasma membrane intrinsic protein (PIP) subfamily, AtPIP1;4 and AtPIP2;5, which are upregulated. In order to assess whether this regulation is general or accession-specific we monitored the gene expression of all PIPs in five Arabidopsis accessions. The overall drought regulation of PIPs was well conserved for all five accessions tested, suggesting a general and fundamental physiological role of this drought response. In addition, significant differences among accessions were identified for transcripts of three PIP genes. Principal component analysis showed that most of the PIP transcriptional variation during drought stress could be explained by one variable linked to leaf water content. Promoter-GUS constructs of AtPIP1;4, AtPIP2;5 and also AtPIP2;6, which is unresponsive to drought stress, had distinct expression patterns concentrated in the base of the leaf petioles and parts of the flowers. The presence of drought stress response elements within the 1.6-kb promoter regions of AtPIP1;4 and AtPIP2;5 was demonstrated by comparing transcription of the promoter reporter construct and the endogenous gene upon drought stress. Analysis by ATTED-II and other web-based bioinformatical tools showed that several of the MIPs downregulated upon drought are strongly co-expressed, whereas AtPIP1;4, AtPIP2;5 and AtPIP2;6 are not co-expressed.

Immune checkpoint inhibitors increase T cell immunity during SARS-CoV-2 infection.

The COVID-19 pandemic has spread worldwide, yet the role of antiviral T cell immunity during infection and the contribution of immune checkpoints remain unclear. By prospectively following a cohort of 292 patients with melanoma, half of which treated with immune checkpoint inhibitors (ICIs), we identified 15 patients with acute or convalescent COVID-19 and investigated their transcriptomic, proteomic, and cellular profiles. We found that ICI treatment was not associated with severe COVID-19 and did not alter the induction of inflammatory and type I interferon responses. In-depth phenotyping demonstrated expansion of CD8 effector memory T cells, enhanced T cell activation, and impaired plasmablast induction in ICI-treated COVID-19 patients. The evaluation of specific adaptive immunity in convalescent patients showed higher spike (S), nucleoprotein (N), and membrane (M) antigen-specific T cell responses and similar induction of spike-specific antibody responses. Our findings provide evidence that ICI during COVID-19 enhanced T cell immunity without exacerbating inflammation.

Integrative analysis of gene expression and copy number alterations using canonical correlation analysis.

BACKGROUND: With the rapid development of new genetic measurement methods, several types of genetic alterations can be quantified in a high-throughput manner. While the initial focus has been on investigating each data set separately, there is an increasing interest in studying the correlation structure between two or more data sets. Multivariate methods based on Canonical Correlation Analysis (CCA) have been proposed for integrating paired genetic data sets. The high dimensionality of microarray data imposes computational difficulties, which have been addressed for instance by studying the covariance structure of the data, or by reducing the number of variables prior to applying the CCA. In this work, we propose a new method for analyzing high-dimensional paired genetic data sets, which mainly emphasizes the correlation structure and still permits efficient application to very large data sets. The method is implemented by translating a regularized CCA to its dual form, where the computational complexity depends mainly on the number of samples instead of the number of variables. The optimal regularization parameters are chosen by cross-validation. We apply the regularized dual CCA, as well as a classical CCA preceded by a dimension-reducing Principal Components Analysis (PCA), to a paired data set of gene expression changes and copy number alterations in leukemia. RESULTS: Using the correlation-maximizing methods, regularized dual CCA and PCA+CCA, we show that without pre-selection of known disease-relevant genes, and without using information about clinical class membership, an exploratory analysis singles out two patient groups, corresponding to well-known leukemia subtypes. Furthermore, the variables showing the highest relevance to the extracted features agree with previous biological knowledge concerning copy number alterations and gene expression changes in these subtypes. Finally, the correlation-maximizing methods are shown to yield results which are more biologically interpretable than those resulting from a covariance-maximizing method, and provide different insight compared to when each variable set is studied separately using PCA. CONCLUSIONS: We conclude that regularized dual CCA as well as PCA+CCA are useful methods for exploratory analysis of paired genetic data sets, and can be efficiently implemented also when the number of variables is very large.

Early changes in the hypothalamic region in prodromal Huntington disease revealed by MRI analysis.

Huntington disease (HD) is a fatal neurodegenerative disorder caused by an expanded CAG repeat. Its length can be used to estimate the time of clinical diagnosis, which is defined by overt motor symptoms. Non-motor symptoms begin before motor onset, and involve changes in hypothalamus-regulated functions such as sleep, emotion and metabolism. Therefore we hypothesized that hypothalamic changes occur already prior to the clinical diagnosis. We performed voxel-based morphometry and logistic regression analyses of cross-sectional MR images from 220 HD gene carriers and 75 controls in the Predict-HD study. We show that changes in the hypothalamic region are detectable before clinical diagnosis and that its grey matter contents alone are sufficient to distinguish HD gene carriers from control cases. In conclusion, our study shows, for the first time, that alterations in grey matter contents in the hypothalamic region occur at least a decade before clinical diagnosis in HD using MRI.

Machine learning can identify newly diagnosed patients with CLL at high risk of infection.

Infections have become the major cause of morbidity and mortality among patients with chronic lymphocytic leukemia (CLL) due to immune dysfunction and cytotoxic CLL treatment. Yet, predictive models for infection are missing. In this work, we develop the CLL Treatment-Infection Model (CLL-TIM) that identifies patients at risk of infection or CLL treatment within 2 years of diagnosis as validated on both internal and external cohorts. CLL-TIM is an ensemble algorithm composed of 28 machine learning algorithms based on data from 4,149 patients with CLL. The model is capable of dealing with heterogeneous data, including the high rates of missing data to be expected in the real-world setting, with a precision of 72% and a recall of 75%. To address concerns regarding the use of complex machine learning algorithms in the clinic, for each patient with CLL, CLL-TIM provides explainable predictions through uncertainty estimates and personalized risk factors.

Sex differences in IL-17 contribute to chronicity in male versus female urinary tract infection.

Sex-based differences influence incidence and outcome of infectious disease. Women have a significantly greater incidence of urinary tract infection (UTI) than men, yet, conversely, male UTI is more persistent with greater associated morbidity. Mechanisms underlying these sex-based differences are unknown, in part due to a lack of experimental models. We optimized a model to transurethrally infect male mice and directly compared UTI in both sexes. Although both sexes were initially equally colonized by uropathogenic E. coli, only male and testosterone-treated female mice remained chronically infected for up to 4 weeks. Female mice had more robust innate responses, including higher IL-17 expression, and increased γδ T cells and group 3 innate lymphoid cells in the bladder following infection. Accordingly, neutralizing IL-17 abolished resolution in female mice, identifying a cytokine pathway necessary for bacterial clearance. Our findings support the concept that sex-based responses to UTI contribute to impaired innate immunity in males and provide a rationale for non-antibiotic-based immune targeting to improve the response to UTI.

Human thymopoiesis is influenced by a common genetic variant within the TCRA-TCRD locus.

The thymus is the primary lymphoid organ where naïve T cells are generated; however, with the exception of age, the parameters that govern its function in healthy humans remain unknown. We characterized the variability of thymic function among 1000 age- and sex-stratified healthy adults of the Milieu Intérieur cohort, using quantification of T cell receptor excision circles (TRECs) in peripheral blood T cells as a surrogate marker of thymopoiesis. Age and sex were the only nonheritable factors identified that affect thymic function. TREC amounts decreased with age and were higher in women compared to men. In addition, a genome-wide association study revealed a common variant (rs2204985) within the T cell receptor TCRA-TCRD locus, between the DD2 and DD3 gene segments, which associated with TREC amounts. Strikingly, transplantation of human hematopoietic stem cells with the rs2204985 GG genotype into immunodeficient mice led to thymopoiesis with higher TRECs, increased thymocyte counts, and a higher TCR repertoire diversity. Our population immunology approach revealed a genetic locus that influences thymopoiesis in healthy adults, with potentially broad implications in precision medicine.