Salmonella Typhimurium (S. Tm) is a leading cause of diarrhea. The disease is triggered by pathogen invasion into the gut epithelium. Invasion is attributed to the SPI-1 type 3 secretion system (T1). T1 injects effector proteins into epithelial cells and thereby elicits rearrangements of the host cellular actin cytoskeleton and pathogen invasion. The T1 effector proteins SopE, SopB, SopE2 and SipA are contributing to this. However, the host cell factors contributing to invasion are still not completely understood. To address this question comprehensively, we used Hela tissue culture cells, a genome-wide siRNA library, a modified gentamicin protection assay and S. TmSipA, a sopBsopE2sopE mutant which strongly relies on the T1 effector protein SipA to invade host cells. We found that S. TmSipA invasion does not elicit membrane ruffles, nor promote the entry of non-invasive bacteria "in trans". However, SipA-mediated infection involved the SPIRE family of actin nucleators, besides well-established host cell factors (WRC, ARP2/3, RhoGTPases, COPI). Stage-specific follow-up assays and knockout fibroblasts indicated that SPIRE1 and SPIRE2 are involved in different steps of the S. Tm infection process. Whereas SPIRE1 interferes with bacterial binding, SPIRE2 influences intracellular replication of S. Tm. Hence, these two proteins might fulfill non-redundant functions in the pathogen-host interaction. The lack of co-localization hints to a short, direct interaction between S. Tm and SPIRE proteins or to an indirect effect.

Cancer has long been understood as a somatic evolutionary process, but many details of tumor progression remain elusive. Here, we present BitPhylogenyBitPhylogeny, a probabilistic framework to reconstruct intra-tumor evolutionary pathways. Using a full Bayesian approach, we jointly estimate the number and composition of clones in the sample as well as the most likely tree connecting them. We validate our approach in the controlled setting of a simulation study and compare it against several competing methods. In two case studies, we demonstrate how BitPhylogeny BitPhylogeny reconstructs tumor phylogenies from methylation patterns in colon cancer and from single-cell exomes in myeloproliferative neoplasm.

Reconstructing the evolution of tumors is a key aspect towards the identification of appropriate cancer therapies. The task is challenging because tumors evolve as heterogeneous cell populations. Single-cell sequencing holds the promise of resolving the heterogeneity of tumors; however, it has its own challenges including elevated error rates, allelic drop-out, and uneven coverage. Here, we develop a new approach to mutation detection in individual tumor cells by leveraging the evolutionary relationship among cells. Our method, called SCIΦ, jointly calls mutations in individual cells and estimates the tumor phylogeny among these cells. Employing a Markov Chain Monte Carlo scheme enables us to reliably call mutations in each single cell even in experiments with high drop-out rates and missing data. We show that SCIΦ outperforms existing methods on simulated data and applied it to different real-world datasets, namely a whole exome breast cancer as well as a panel acute lymphoblastic leukemia dataset.

Annotation and interpretation of DNA aberrations identified through next-generation sequencing is becoming an increasingly important task. Even more so in the context of data analysis pipelines for medical applications, where genomic aberrations are associated with phenotypic and clinical features. Here we describe a workflow to identify potential gene targets in aberrated genes or pathways and their corresponding drugs. To this end, we provide the R/Bioconductor package rDGIdb, an R wrapper to query the drug-gene interaction database (DGIdb). DGIdb accumulates drug-gene interaction data from 15 different resources and allows filtering on different levels. The rDGIdb package makes these resources and tools available to R users. Moreover, rDGIdb queries can be automated through incorporation of the rDGIdb package into NGS sequencing pipelines.

Personalized, precision, P4, or stratified medicine is understood as a medical approach in which patients are stratified based on their disease subtype, risk, prognosis, or treatment response using specialized diagnostic tests. The key idea is to base medical decisions on individual patient characteristics, including molecular and behavioral biomarkers, rather than on population averages. Personalized medicine is deeply connected to and dependent on data science, specifically machine learning (often named Artificial Intelligence in the mainstream media). While during recent years there has been a lot of enthusiasm about the potential of 'big data' and machine learning-based solutions, there exist only few examples that impact current clinical practice. The lack of impact on clinical practice can largely be attributed to insufficient performance of predictive models, difficulties to interpret complex model predictions, and lack of validation via prospective clinical trials that demonstrate a clear benefit compared to the standard of care. In this paper, we review the potential of state-of-the-art data science approaches for personalized medicine, discuss open challenges, and highlight directions that may help to overcome them in the future.

Recent advancements of sequencing technology have opened up unprecedented opportunities in many application areas. Virus samples can now be sequenced efficiently with very deep coverage to infer the genetic diversity of the underlying virus populations. Several sequencing platforms with different underlying technologies and performance characteristics are available for viral diversity studies. Here, we investigate how the differences between two common platforms provided by 454/Roche and Illumina affect viral diversity estimation and the reconstruction of viral haplotypes. Using a mixture of ten HIV clones sequenced with both platforms and additional simulation experiments, we assessed the trade-off between sequencing coverage, read length, and error rate. For fixed costs, short Illumina reads can be generated at higher coverage and allow for detecting variants at lower frequencies. They can also be sufficient to assess the diversity of the sample if sequences are dissimilar enough, but, in general, assembly of full-length haplotypes is feasible only with the longer 454/Roche reads. The quantitative comparison highlights the advantages and disadvantages of both platforms and provides guidance for the design of viral diversity studies.

Many viruses, including the clinically relevant RNA viruses HIV (human immunodeficiency virus) and HCV (hepatitis C virus), exist in large populations and display high genetic heterogeneity within and between infected hosts. Assessing intra-patient viral genetic diversity is essential for understanding the evolutionary dynamics of viruses, for designing effective vaccines, and for the success of antiviral therapy. Next-generation sequencing (NGS) technologies allow the rapid and cost-effective acquisition of thousands to millions of short DNA sequences from a single sample. However, this approach entails several challenges in experimental design and computational data analysis. Here, we review the entire process of inferring viral diversity from sample collection to computing measures of genetic diversity. We discuss sample preparation, including reverse transcription and amplification, and the effect of experimental conditions on diversity estimates due to in vitro base substitutions, insertions, deletions, and recombination. The use of different NGS platforms and their sequencing error profiles are compared in the context of various applications of diversity estimation, ranging from the detection of single nucleotide variants (SNVs) to the reconstruction of whole-genome haplotypes. We describe the statistical and computational challenges arising from these technical artifacts, and we review existing approaches, including available software, for their solution. Finally, we discuss open problems, and highlight successful biomedical applications and potential future clinical use of NGS to estimate viral diversity.

HIV latency is a major obstacle to curing infection. Current strategies to eradicate HIV aim at increasing transcription of the latent provirus. In the present study we observed that latently infected CD4+ T cells from HIV-infected individuals failed to produce viral particles upon ex vivo exposure to SAHA (vorinostat), despite effective inhibition of histone deacetylases. To identify steps that were not susceptible to the action of SAHA or other latency reverting agents, we used a primary CD4+ T cell model, joint host and viral RNA sequencing, and a viral-encoded reporter. This model served to investigate the characteristics of latently infected cells, the dynamics of HIV latency, and the process of reactivation induced by various stimuli. During latency, we observed persistence of viral transcripts but only limited viral translation. Similarly, the reactivating agents SAHA and disulfiram successfully increased viral transcription, but failed to effectively enhance viral translation, mirroring the ex vivo data. This study highlights the importance of post-transcriptional blocks as one mechanism leading to HIV latency that needs to be relieved in order to purge the viral reservoir.

Next-generation sequencing (NGS) is a valuable tool for the detection and quantification of HIV-1 variants in vivo. However, these technologies require detailed characterization and control of artificially induced errors to be applicable for accurate haplotype reconstruction. To investigate the occurrence of substitutions, insertions, and deletions at the individual steps of RT-PCR and NGS, 454 pyrosequencing was performed on amplified and non-amplified HIV-1 genomes. Artificial recombination was explored by mixing five different HIV-1 clonal strains (5-virus-mix) and applying different RT-PCR conditions followed by 454 pyrosequencing. Error rates ranged from 0.04-0.66% and were similar in amplified and non-amplified samples. Discrepancies were observed between forward and reverse reads, indicating that most errors were introduced during the pyrosequencing step. Using the 5-virus-mix, non-optimized, standard RT-PCR conditions introduced artificial recombinants in a fraction of at least 30% of the reads that subsequently led to an underestimation of true haplotype frequencies. We minimized the fraction of recombinants down to 0.9-2.6% by optimized, artifact-reducing RT-PCR conditions. This approach enabled correct haplotype reconstruction and frequency estimations consistent with reference data obtained by single genome amplification. RT-PCR conditions are crucial for correct frequency estimation and analysis of haplotypes in heterogeneous virus populations. We developed an RT-PCR procedure to generate NGS data useful for reliable haplotype reconstruction and quantification.

Intra-tumor heterogeneity poses substantial challenges for cancer treatment. A tumor's composition can be deduced by reconstructing its mutational history. Central to current approaches is the infinite sites assumption that every genomic position can only mutate once over the lifetime of a tumor. The validity of this assumption has never been quantitatively assessed. We developed a rigorous statistical framework to test the infinite sites assumption with single-cell sequencing data. Our framework accounts for the high noise and contamination present in such data. We found strong evidence for the same genomic position being mutationally affected multiple times in individual tumors for 11 of 12 single-cell sequencing data sets from a variety of human cancers. Seven cases involved the loss of earlier mutations, five of which occurred at sites unaffected by large-scale genomic deletions. Four cases exhibited a parallel mutation, potentially indicating convergent evolution at the base pair level. Our results refute the general validity of the infinite sites assumption and indicate that more complex models are needed to adequately quantify intra-tumor heterogeneity for more effective cancer treatment.

HIV-1 infects CD4+ T cells and completes its replication cycle in approximately 24 hours. We employed repeated measurements in a standardized cell system and rigorous mathematical modeling to characterize the emergence of the viral replication intermediates and their impact on the cellular transcriptional response with high temporal resolution. We observed 7,991 (73%) of the 10,958 expressed genes to be modulated in concordance with key steps of viral replication. Fifty-two percent of the overall variability in the host transcriptome was explained by linear regression on the viral life cycle. This profound perturbation of cellular physiology was investigated in the light of several regulatory mechanisms, including transcription factors, miRNAs, host-pathogen interaction, and proviral integration. Key features were validated in primary CD4+ T cells, and with viral constructs using alternative entry strategies. We propose a model of early massive cellular shutdown and progressive upregulation of the cellular machinery to complete the viral life cycle.

Next-generation sequencing technologies can be used to analyse genetically heterogeneous samples at unprecedented detail. The high coverage achievable with these methods enables the detection of many low-frequency variants. However, sequencing errors complicate the analysis of mixed populations and result in inflated estimates of genetic diversity. We developed a probabilistic Bayesian approach to minimize the effect of errors on the detection of minority variants. We applied it to pyrosequencing data obtained from a 1.5-kb-fragment of the HIV-1 gag/pol gene in two control and two clinical samples. The effect of PCR amplification was analysed. Error correction resulted in a two- and five-fold decrease of the pyrosequencing base substitution rate, from 0.05% to 0.03% and from 0.25% to 0.05% in the non-PCR and PCR-amplified samples, respectively. We were able to detect viral clones as rare as 0.1% with perfect sequence reconstruction. Probabilistic haplotype inference outperforms the counting-based calling method in both precision and recall. Genetic diversity observed within and between two clinical samples resulted in various patterns of phenotypic drug resistance and suggests a close epidemiological link. We conclude that pyrosequencing can be used to investigate genetically diverse samples with high accuracy if technical errors are properly treated.

Cells respond to stress by remodeling their transcriptome through transcription and degradation. Xrn1p-dependent degradation in P-bodies is the most prevalent decay pathway, yet, P-bodies may facilitate not only decay, but also act as a storage compartment. However, which and how mRNAs are selected into different degradation pathways and what determines the fate of any given mRNA in P-bodies remain largely unknown. We devised a new method to identify both common and stress-specific mRNA subsets associated with P-bodies. mRNAs targeted for degradation to P-bodies, decayed with different kinetics. Moreover, the localization of a specific set of mRNAs to P-bodies under glucose deprivation was obligatory to prevent decay. Depending on its client mRNA, the RNA-binding protein Puf5p either promoted or inhibited decay. Furthermore, the Puf5p-dependent storage of a subset of mRNAs in P-bodies under glucose starvation may be beneficial with respect to chronological lifespan.

MiRNAs and other small noncoding RNAs (sncRNAs) are key players in post-transcriptional gene regulation. HIV-1 derived small noncoding RNAs (sncRNAs) have been described in HIV-1 infected cells, but their biological functions still remain to be elucidated. Here, we approached the question whether viral sncRNAs may play a role in the RNA interference (RNAi) pathway or whether viral mRNAs are targeted by cellular miRNAs in human monocyte derived macrophages (MDM).

Clonal diversity is a consequence of cancer cell evolution driven by Darwinian selection. Precise characterization of clonal architecture is essential to understand the evolutionary history of tumor development and its association with treatment resistance. Here, using a single-cell DNA sequencing, we report the clonal architecture and mutational histories of 123 acute myeloid leukemia (AML) patients. The single-cell data reveals cell-level mutation co-occurrence and enables reconstruction of mutational histories characterized by linear and branching patterns of clonal evolution, with the latter including convergent evolution. Through xenotransplantion, we show leukemia initiating capabilities of individual subclones evolving in parallel. Also, by simultaneous single-cell DNA and cell surface protein analysis, we illustrate both genetic and phenotypic evolution in AML. Lastly, single-cell analysis of longitudinal samples reveals underlying evolutionary process of therapeutic resistance. Together, these data unravel clonal diversity and evolution patterns of AML, and highlight their clinical relevance in the era of precision medicine.

The HIV-1 reservoir is the major hurdle to a cure. We here evaluate viral and host characteristics associated with reservoir size and long-term dynamics in 1,057 individuals on suppressive antiretroviral therapy for a median of 5.4 years. At the population level, the reservoir decreases with diminishing differences over time, but increases in 26.6% of individuals. Viral blips and low-level viremia are significantly associated with slower reservoir decay. Initiation of ART within the first year of infection, pretreatment viral load, and ethnicity affect reservoir size, but less so long-term dynamics. Viral blips and low-level viremia are thus relevant for reservoir and cure studies.

Germinal centers (GCs) are specialized compartments within the secondary lymphoid organs where B cells proliferate, differentiate, and mutate their antibody genes in response to the presence of foreign antigens. Through the GC lifespan, interclonal competition between B cells leads to increased affinity of the B cell receptors for antigens accompanied by a loss of clonal diversity, although the mechanisms underlying clonal dynamics are not completely understood. We present here a multi-scale quantitative model of the GC reaction that integrates an intracellular component, accounting for the genetic events that shape B cell differentiation, and an extracellular stochastic component, which accounts for the random cellular interactions within the GC. In addition, B cell receptors are represented as sequences of nucleotides that mature and diversify through somatic hypermutations. We exploit extensive experimental characterizations of the GC dynamics to parameterize our model, and visualize affinity maturation by means of evolutionary phylogenetic trees. Our explicit modeling of B cell maturation enables us to characterise the evolutionary processes and competition at the heart of the GC dynamics, and explains the emergence of clonal dominance as a result of initially small stochastic advantages in the affinity to antigen. Interestingly, a subset of the GC undergoes massive expansion of higher-affinity B cell variants (clonal bursts), leading to a loss of clonal diversity at a significantly faster rate than in GCs that do not exhibit clonal dominance. Our work contributes towards an in silico vaccine design, and has implications for the better understanding of the mechanisms underlying autoimmune disease and GC-derived lymphomas.

In melanoma, a switch from a proliferative melanocytic to an invasive mesenchymal phenotype is based on dramatic transcriptional reprogramming which involves complex interactions between a variety of signaling pathways and their downstream transcriptional regulators. TGFβ/SMAD, Hippo/YAP/TAZ, and Wnt/β-catenin signaling pathways are major inducers of transcriptional reprogramming and converge at several levels. Here, we report that TGFβ/SMAD, YAP/TAZ, and β-catenin are all required for a proliferative-to-invasive phenotype switch. Loss and gain of function experimentation, global gene expression analysis, and computational nested effects models revealed the hierarchy between these signaling pathways and identified shared target genes. SMAD-mediated transcription at the top of the hierarchy leads to the activation of YAP/TAZ and of β-catenin, with YAP/TAZ governing an essential subprogram of TGFβ-induced phenotype switching. Wnt/β-catenin signaling is situated further downstream and exerts a dual role: it promotes the proliferative, differentiated melanoma cell phenotype and it is essential but not sufficient for SMAD or YAP/TAZ-induced phenotype switching. The results identify epistatic interactions among the signaling pathways underlying melanoma phenotype switching and highlight the priorities in targets for melanoma therapy.

In human beings, there are five reported variants of concern of severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2). However, in contrast to human beings, descriptions of infections of animals with specific variants are still rare. The aim of this study is to systematically investigate SARS-CoV-2 infections in companion animals in close contact with SARS-CoV-2-positive owners ("COVID-19 households") with a focus on the Delta variant. Samples, obtained from companion animals and their owners were analyzed using a real-time reverse transcriptase-polymerase chain reaction (RT-qPCR) and next-generation sequencing (NGS). Animals were also tested for antibodies and neutralizing activity against SARS-CoV-2. Eleven cats and three dogs in nine COVID-19-positive households were RT-qPCR and/or serologically positive for the SARS-CoV-2 Delta variant. For seven animals, the genetic sequence could be determined. The animals were infected by one of the pangolin lineages B.1.617.2, AY.4, AY.43 and AY.129 and between zero and three single-nucleotide polymorphisms (SNPs) were detected between the viral genomes of animals and their owners, indicating within-household transmission between animal and owner and in multi-pet households also between the animals. NGS data identified SNPs that occur at a higher frequency in the viral sequences of companion animals than in viral sequences of humans, as well as SNPs, which were exclusively found in the animals investigated in the current study and not in their owners. In conclusion, our study is the first to describe the SARS-CoV-2 Delta variant transmission to animals in Switzerland and provides the first-ever description of Delta-variant pangolin lineages AY.129 and AY.4 in animals. Our results reinforce the need of a One Health approach in the monitoring of SARS-CoV-2 in animals.

Forward genetic screens with genome-wide CRISPR libraries are powerful tools for resolving cellular circuits and signaling pathways. Applying this technology to organoids, however, has been hampered by technical limitations. Here we report improved accuracy and robustness for pooled-library CRISPR screens by capturing sgRNA integrations in single organoids, substantially reducing required cell numbers for genome-scale screening. We applied our approach to wild-type and APC mutant human intestinal organoids to identify genes involved in resistance to TGF-β-mediated growth restriction, a key process during colorectal cancer progression, and validated hits including multiple subunits of the tumor-suppressive SWI/SNF chromatin remodeling complex. Mutations within these genes require concurrent inactivation of APC to promote TGF-β resistance and attenuate TGF-β target gene transcription. Our approach can be applied to a variety of assays and organoid types to facilitate biological discovery in primary 3D tissue models.