Search

Abstract Deep single-cell multi-omic profiling offers a promising approach to understand and overcome drug resistance in relapsed or refractory (rr) acute myeloid leukemia (AML). Here, we combine single-cell ex vivo drug profiling (pharmacoscopy) with single-cell and bulk DNA, RNA, and protein analyses, alongside clinical data from 21 rrAML patients. Unsupervised data integration reveals reduced ex vivo response to the Bcl-2 inhibitor venetoclax (VEN) in patients treated with both a hypomethylating agent (HMA) and VEN, compared to those pre-exposed to chemotherapy or HMA alone. Integrative analysis identifies both known and unreported mechanisms of innate and treatment-related VEN resistance and suggests alternative treatments, like targeting increased proliferation with the PLK inhibitor volasertib. Additionally, high CD36 expression in VEN-resistant blasts associates with sensitivity to CD36-targeted antibody treatment ex vivo. This study demonstrates how single-cell multi-omic profiling can uncover drug resistance mechanisms and treatment vulnerabilities, providing a valuable resource for future AML research.

Abstract Understanding the complex background of cancer requires genotype-phenotype information in single-cell resolution. Here, we perform long-read single-cell RNA sequencing (scRNA-seq) on clinical samples from three ovarian cancer patients presenting with omental metastasis and increase the PacBio sequencing depth to 12,000 reads per cell. Our approach captures 152,000 isoforms, of which over 52,000 were not previously reported. Isoform-level analysis accounting for non-coding isoforms reveals 20% overestimation of protein-coding gene expression on average. We also detect cell type-specific isoform and poly-adenylation site usage in tumor and mesothelial cells, and find that mesothelial cells transition into cancer-associated fibroblasts in the metastasis, partly through the TGF-β/miR-29/Collagen axis. Furthermore, we identify gene fusions, including an experimentally validated IGF2BP2::TESPA1 fusion, which is misclassified as high TESPA1 expression in matched short-read data, and call mutations confirmed by targeted NGS cancer gene panel results. With these findings, we envision long-read scRNA-seq to become increasingly relevant in oncology and personalized medicine.

Abstract Multi-omics profiling by CITE-seq bridges the RNA-protein gap in single-cell analysis but has been largely applied to liquid biopsies. Applying CITE-seq to clinically relevant solid biopsies to characterize healthy tissue and the tumor microenvironment is an essential next step in single-cell translational studies. In this study, gating of cell populations based on their transcriptome signatures for use in cell type-specific ridge plots allowed identification of positive antibody signals and setting of manual thresholds. Next, we compare five skin dissociation protocols by taking into account dissociation efficiency, captured cell type heterogeneity and recovered surface proteome. To assess the effect of enzymatic digestion on transcriptome and epitope expression in immune cell populations, we analyze peripheral blood mononuclear cells (PBMCs) with and without dissociation. To further assess the RNA-protein gap, RNA-protein we perform codetection and correlation analyses on thresholded protein values. Finally, in a proof-of-concept study, using protein abundance analysis on selected surface markers in a cohort of healthy skin, primary, and metastatic melanoma we identify CD56 surface marker expression on metastatic melanoma cells, which was further confirmed by multiplex immunohistochemistry. This work provides practical guidelines for processing and analysis of clinically relevant solid tissue biopsies for biomarker discovery.

Geographical hotspots with high frequency of multi-drug resistant tuberculosis (MDR-TB) have been observed in several locations, such as the country of Georgia. Here, the authors analyse genomic sequences from tuberculosis bacteria isolated from Georgia to show that the transmission fitness of MDR-TB strains is heterogeneous, and highly drug-resistant and transmissible isolates contribute to the emergence and maintenance of MDR-TB hotspots.

Many different animal species are susceptible to SARS-CoV-2, including a few Canidae (domestic dog and raccoon dog). So far, only experimental evidence is available concerning SARS-CoV-2 infections in red foxes (Vulpes vulpes). This is the first report of SARS-CoV-2 RNA detection in a sample from a red fox. The RT-qPCR-positive fox was zoo-kept together with another fox and two bears in the Swiss Canton of Zurich. Combined material from a conjunctival and nasal swab collected for canine distemper virus diagnostics tested positive for SARS-CoV-2 RNA with Ct values of 36.9 (E gene assay) and 35.7 (RdRp gene assay). The sample was analysed for SARS-CoV-2 within a research project testing residual routine diagnostic samples from different animal species submitted between spring 2020 and December 2022 to improve knowledge on SARS-CoV-2 infections within different animal species and investigate their potential role in a One Health context. Within this project, 246 samples from 153 different animals from Swiss zoos and other wild animal species all tested SARS-CoV-2 RT-qPCR and/or serologically negative so far, except for the reported fox. The source of SARS-CoV-2 in the fox is unknown. The fox disappeared within the naturally structured enclosure, and the cadaver was not found. No further control measures were undertaken.

Fresh vs. slow-frozen tissues from various malignancies are compared for the establishment of 2D, 3D and ex vivo cultures, as well as for scRNAseq analysis, and found to be comparable and suitable for cancer research.

Abstract Background The continued spread of SARS-CoV-2 and emergence of new variants with higher transmission rates and/or partial resistance to vaccines has further highlighted the need for large-scale testing and genomic surveillance. However, current diagnostic testing (e.g., PCR) and genomic surveillance methods (e.g., whole genome sequencing) are performed separately, thus limiting the detection and tracing of SARS-CoV-2 and emerging variants. Results Here, we developed DeepSARS, a high-throughput platform for simultaneous diagnostic detection and genomic surveillance of SARS-CoV-2 by the integration of molecular barcoding, targeted deep sequencing, and computational phylogenetics. DeepSARS enables highly sensitive viral detection, while also capturing genomic diversity and viral evolution. We show that DeepSARS can be rapidly adapted for identification of emerging variants, such as alpha, beta, gamma, and delta strains, and profile mutational changes at the population level. Conclusions DeepSARS sets the foundation for quantitative diagnostics that capture viral evolution and diversity. Graphical abstract DeepSARS uses molecular barcodes (BCs) and multiplexed targeted deep sequencing (NGS) to enable simultaneous diagnostic detection and genomic surveillance of SARS-CoV-2. Image was created using Biorender.com .

Summary: The various stages of epithelial-mesenchymal transition (EMT) generate phenotypically heterogeneous populations of cells. Here, we detail a dual recombinase lineage tracing system using a transgenic mouse model of metastatic breast cancer to trace and characterize breast cancer cells at different EMT stages. We describe analytical steps to label cancer cells at an early partial or a late full EMT state, followed by tracking their behavior in tumor slice cultures. We then characterize their transcriptome by five-cell RNA sequencing.For complete details on the use and execution of this protocol, please refer to Luond et al. (2021). : Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Abstract Background Tumor-specific genomic aberrations are routinely determined by high-throughput genomic measurements. It remains unclear how complex genome alterations affect molecular networks through changing protein levels and consequently biochemical states of tumor tissues. Results Here, we investigate the propagation of genomic effects along the axis of gene expression during prostate cancer progression. We quantify genomic, transcriptomic, and proteomic alterations based on 105 prostate samples, consisting of benign prostatic hyperplasia regions and malignant tumors, from 39 prostate cancer patients. Our analysis reveals the convergent effects of distinct copy number alterations impacting on common downstream proteins, which are important for establishing the tumor phenotype. We devise a network-based approach that integrates perturbations across different molecular layers, which identifies a sub-network consisting of nine genes whose joint activity positively correlates with increasingly aggressive tumor phenotypes and is predictive of recurrence-free survival. Further, our data reveal a wide spectrum of intra-patient network effects, ranging from similar to very distinct alterations on different molecular layers. Conclusions This study uncovers molecular networks with considerable convergent alterations across tumor sites and patients. It also exposes a diversity of network effects: we could not identify a single sub-network that is perturbed in all high-grade tumor regions.

Current methods to generate sequence-function data at large scale are either technically complex or limited to specific applications. Here the authors introduce DNA-based phenotypic recording to overcome these limitations and enable deep learning for accurate prediction of function from sequence.

Abstract Background The genus Macrostomum consists of small free-living flatworms and contains Macrostomum lignano, which has been used in investigations of ageing, stem cell biology, bioadhesion, karyology, and sexual selection in hermaphrodites. Two types of mating behaviour occur within this genus. Some species, including M. lignano, mate via reciprocal copulation, where, in a single mating, both partners insert their male copulatory organ into the female storage organ and simultaneously donate and receive sperm. Other species mate via hypodermic insemination, where worms use a needle-like copulatory organ to inject sperm into the tissue of the partner. These contrasting mating behaviours are associated with striking differences in sperm and copulatory organ morphology. Here we expand the genomic resources within the genus to representatives of both behaviour types and investigate whether genes vary in their rate of evolution depending on their putative function. Results We present de novo assembled transcriptomes of three Macrostomum species, namely M. hystrix, a close relative of M. lignano that mates via hypodermic insemination, M. spirale, a more distantly related species that mates via reciprocal copulation, and finally M. pusillum, which represents a clade that is only distantly related to the other three species and also mates via hypodermic insemination. We infer 23,764 sets of homologous genes and annotate them using experimental evidence from M. lignano. Across the genus, we identify 521 gene families with conserved patterns of differential expression between juvenile vs. adult worms and 185 gene families with a putative expression in the testes that are restricted to the two reciprocally mating species. Further, we show that homologs of putative reproduction-related genes have a higher protein divergence across the four species than genes lacking such annotations and that they are more difficult to identify across the four species, indicating that these genes evolve more rapidly, while genes involved in neoblast function are more conserved. Conclusions This study improves the genus Macrostomum as a model system, by providing resources for the targeted investigation of gene function in a broad range of species. And we, for the first time, show that reproduction-related genes evolve at an accelerated rate in flatworms.

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.

Background: In December 2020, the United Kingdom (UK) reported a SARS-CoV-2 Variant of Concern (VoC) which is now named B.1.1.7. Based on initial data from the UK and later data from other countries, this variant was estimated to have a transmission fitness advantage of around 40–80 % (Volz et al., 2021; Leung et al., 2021; Davies et al., 2021). Aim: This study aims to estimate the transmission fitness advantage and the effective reproductive number of B.1.1.7 through time based on data from Switzerland. Methods: We generated whole genome sequences from 11.8 % of all confirmed SARS-CoV-2 cases in Switzerland between 14 December 2020 and 11 March 2021. Based on these data, we determine the daily frequency of the B.1.1.7 variant and quantify the variant’s transmission fitness advantage on a national and a regional scale. Results: We estimate B.1.1.7 had a transmission fitness advantage of 43–52 % compared to the other variants circulating in Switzerland during the study period. Further, we estimate B.1.1.7 had a reproductive number above 1 from 01 January 2021 until the end of the study period, compared to below 1 for the other variants. Specifically, we estimate the reproductive number for B.1.1.7 was 1.24 [1.07–1.41] from 01 January until 17 January 2021 and 1.18 [1.06–1.30] from 18 January until 01 March 2021 based on the whole genome sequencing data. From 10 March to 16 March 2021, once B.1.1.7 was dominant, we estimate the reproductive number was 1.14 [1.00–1.26] based on all confirmed cases. For reference, Switzerland applied more non-pharmaceutical interventions to combat SARS-CoV-2 on 18 January 2021 and lifted some measures again on 01 March 2021. Conclusion: The observed increase in B.1.1.7 frequency in Switzerland during the study period is as expected based on observations in the UK. In absolute numbers, B.1.1.7 increased exponentially with an estimated doubling time of around 2–3.5 weeks. To monitor the ongoing spread of B.1.1.7, our plots are available online.

Tumour progression is an evolutionary process in which different clones evolve over time, leading to intra-tumour heterogeneity. Interactions between clones can affect tumour evolution and hence disease progression and treatment outcome. Intra-tumoural pairs of mutations that are overrepresented in a co-occurring or clonally exclusive fashion over a cohort of patient samples may be suggestive of a synergistic effect between the different clones carrying these mutations. We therefore developed a novel statistical testing framework, called GeneAccord, to identify such gene pairs that are altered in distinct subclones of the same tumour. We analysed our framework for calibration and power. By comparing its performance to baseline methods, we demonstrate that to control type I errors, it is essential to account for the evolutionary dependencies among clones. In applying GeneAccord to the single-cell sequencing of a cohort of 123 acute myeloid leukaemia patients, we find 1 clonally co-occurring and 8 clonally exclusive gene pairs. The clonally exclusive pairs mostly involve genes of the key signalling pathways.

Summary: Background: South Africa has a high burden of rifampicin-resistant tuberculosis (including multidrug-resistant [MDR] tuberculosis), with increasing rifampicin-monoresistant (RMR) tuberculosis over time. Resistance acquisition during first-line tuberculosis treatment could be a key contributor to this burden, and HIV might increase the risk of acquiring rifampicin resistance. We assessed whether HIV during previous treatment was associated with RMR tuberculosis and resistance acquisition among a retrospective cohort of patients with MDR or rifampicin-resistant tuberculosis. Methods: In this retrospective cohort study, we included all patients routinely diagnosed with MDR or rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa, between Jan 1, 2008, and Dec 31, 2017. Patient-level data were obtained from a prospective database, complemented by data on previous tuberculosis treatment and HIV from a provincial health data exchange. Stored MDR or rifampicin-resistant tuberculosis isolates from patients underwent whole-genome sequencing (WGS). WGS data were used to infer resistance acquisition versus transmission, by identifying genomically unique isolates (single nucleotide polymorphism threshold of five). Logistic regression analyses were used to assess factors associated with RMR tuberculosis and genomic uniqueness. Findings: The cohort included 2041 patients diagnosed with MDR or rifampicin-resistant tuberculosis between Jan 1, 2008, and Dec 31, 2017; of those, 463 (22·7%) with RMR tuberculosis and 1354 (66·3%) with previous tuberculosis treatment. In previously treated patients, HIV positivity during previous tuberculosis treatment versus HIV negativity (adjusted odds ratio [OR] 2·07, 95% CI 1·35–3·18), and three or more previous tuberculosis treatment episodes versus one (1·96, 1·21–3·17) were associated with RMR tuberculosis. WGS data showing MDR or rifampicin-resistant tuberculosis were available for 1169 patients; 360 (30·8%) isolates were identified as unique. In previously treated patients, RMR tuberculosis versus MDR tuberculosis (adjusted OR 4·96, 3·40–7·23), HIV positivity during previous tuberculosis treatment (1·71, 1·03–2·84), and diagnosis in 2013–17 (1·42, 1·02–1·99) versus 2008–12, were associated with uniqueness. In previously treated patients with RMR tuberculosis, HIV positivity during previous treatment (adjusted OR 5·13, 1·61–16·32) was associated with uniqueness as was female sex (2·50 [1·18–5·26]). Interpretation: These data suggest that HIV contributes to rifampicin-resistance acquisition during first-line tuberculosis treatment and that this might be driving increasing RMR tuberculosis over time. Large-scale prospective cohort studies are required to further quantify this risk. Funding: Swiss National Science Foundation, South African National Research Foundation, and Wellcome Trust.

Each day, approximately 27,000 people become ill with tuberculosis (TB), and 4,000 die from this disease. Pulmonary TB is the main clinical form of TB, and affects the lungs with a considerably heterogeneous manifestation among patients. Immunomodulation by an interplay of host-, environment-, and pathogen-associated factors partially explains such heterogeneity. Microbial communities residing in the host's airways have immunomodulatory effects, but it is unclear if the inter-individual variability of these microbial communities is associated with the heterogeneity of pulmonary TB. Here, we investigated this possibility by characterizing the microbial composition in the sputum of 334 TB patients from Tanzania, and by assessing its association with three aspects of disease manifestations: sputum mycobacterial load, severe clinical findings, and chest x-ray (CXR) findings. Compositional data analysis of taxonomic profiles based on 16S-rRNA gene amplicon sequencing and on whole metagenome shotgun sequencing, and graph-based inference of microbial associations revealed that the airway microbiome of TB patients was shaped by inverse relationships between Streptococcus and two anaerobes: Selenomonas and Fusobacterium. Specifically, the strength of these microbial associations was negatively correlated with Faith's phylogenetic diversity (PD) and with the accumulation of transient genera. Furthermore, low body mass index (BMI) determined the association between abnormal CXRs and community diversity and composition. These associations were mediated by increased abundance of Selenomonas and Fusobacterium, relative to the abundance of Streptococcus, in underweight patients with lung parenchymal infiltrates and in comparison to those with normal chest x-rays. And last, the detection of herpesviruses and anelloviruses in sputum microbial assemblage was linked to co-infection with HIV. Given the anaerobic metabolism of Selenomonas and Fusobacterium, and the hypoxic environment of lung infiltrates, our results suggest that in underweight TB patients, lung tissue remodeling toward anaerobic conditions favors the growth of Selenomonas and Fusobacterium at the expense of Streptococcus. These new insights into the interplay among particular members of the airway microbiome, BMI, and lung parenchymal lesions in TB patients, add a new dimension to the long-known association between low BMI and pulmonary TB. Our results also drive attention to the airways virome in the context of HIV-TB coinfection.

Summary: High cell viability and recovered cell concentration are typical quality control requirements for single-cell processing and quality data. This protocol describes procedures for sampling, live-cell biobanking, preprocessing for single-cell RNA sequencing, and analysis of fine-needle aspiration (FNA) samples of the skin. The minimally invasive nature of FNA collection is more accepted by patients and allows for frequent longitudinal sampling, resulting in high-quality single-cell sequencing data that capture cellular heterogeneity in clinical samples.

Bi-allelic inactivation of the VHL gene on chromosome 3p is the characteristic feature in most clear cell renal cell carcinomas (ccRCC). Frequent gene alterations were also identified in SETD2, BAP1 and PBRM1, all of which are situated on chromosome 3p and encode histone/chromatin regulators. The relationship between gene mutation, loss of protein expression and the correlations with clinicopathological parameters is important for the understanding of renal cancer progression. We analyzed PBRM1 and BAP1 protein expression as well as the tri-methylation state of H3K36 as a surrogate marker for SETD2 activity in more than 700 RCC samples. In ccRCC loss of nuclear PBRM1 (68%), BAP1 (40%) and H3K36me3 (47%) expression was significantly correlated with each other, advanced tumor stage, poor tumor differentiation (P

Background: Lineage 1 (L1) and 3 (L3) are two lineages of the Mycobacterium tuberculosis complex (MTBC) causing tuberculosis (TB) in humans. L1 and L3 are prevalent around the rim of the Indian Ocean, the region that accounts for most of the world’s new TB cases. Despite their relevance for this region, L1 and L3 remain understudied. Methods: We analyzed 2,938 L1 and 2,030 L3 whole genome sequences originating from 69 countries. We reconstructed the evolutionary history of these two lineages and identified genes under positive selection. Results: We found a strongly asymmetric pattern of migration from South Asia toward neighboring regions, highlighting the historical role of South Asia in the dispersion of L1 and L3. Moreover, we found that several genes were under positive selection, including genes involved in virulence and resistance to antibiotics. For L1 we identified signatures of local adaptation at the esxH locus, a gene coding for a secreted effector that targets the human endosomal sorting complex, and is included in several vaccine candidates. Conclusions: Our study highlights the importance of genetic diversity in the MTBC, and sheds new light on two of the most important MTBC lineages affecting humans.

Background: Lineage 1 (L1) and 3 (L3) are two lineages of the Mycobacterium tuberculosis complex (MTBC) causing tuberculosis (TB) in humans. L1 and L3 are prevalent around the rim of the Indian Ocean, the region that accounts for most of the world’s new TB cases. Despite their relevance for this region, L1 and L3 remain understudied. Methods: We analyzed 2,938 L1 and 2,030 L3 whole genome sequences originating from 69 countries. We reconstructed the evolutionary history of these two lineages and identified genes under positive selection. Results: We found a strongly asymmetric pattern of migration from South Asia toward neighboring regions, highlighting the historical role of South Asia in the dispersion of L1 and L3. Moreover, we found that several genes were under positive selection, including genes involved in virulence and resistance to antibiotics . For L1 we identified signatures of local adaptation at the esxH locus, a gene coding for a secreted effector that targets the human endosomal sorting complex, and is included in several vaccine candidates. Conclusions: Our study highlights the importance of genetic diversity in the MTBC, and sheds new light on two of the most important MTBC lineages affecting humans.

Abstract Background Molecular precision oncology is an emerging practice to improve cancer therapy by decreasing the risk of choosing treatments that lack efficacy or cause adverse events. However, the challenges of integrating molecular profiling into routine clinical care are manifold. From a computational perspective these include the importance of a short analysis turnaround time, the interpretation of complex drug-gene and gene-gene interactions, and the necessity of standardized high-quality workflows. In addition, difficulties faced when integrating molecular diagnostics into clinical practice are ethical concerns, legal requirements, and limited availability of treatment options beyond standard of care as well as the overall lack of awareness of their existence. Methods To the best of our knowledge, we are the first group in Switzerland that established a workflow for personalized diagnostics based on comprehensive high-throughput sequencing of tumors at the clinic. Our workflow, named SwissMTB (Swiss Molecular Tumor Board), links genetic tumor alterations and gene expression to therapeutic options and clinical trial opportunities. The resulting treatment recommendations are summarized in a clinical report and discussed in a molecular tumor board at the clinic to support therapy decisions. Results Here we present results from an observational pilot study including 22 late-stage cancer patients. In this study we were able to identify actionable variants and corresponding therapies for 19 patients. Half of the patients were analyzed retrospectively. In two patients we identified resistance-associated variants explaining lack of therapy response. For five out of eleven patients analyzed before treatment the SwissMTB diagnostic influenced treatment decision. Conclusions SwissMTB enables the analysis and clinical interpretation of large numbers of potentially actionable molecular targets. Thus, our workflow paves the way towards a more frequent use of comprehensive molecular diagnostics in Swiss hospitals.

Abstract miRNAs are small RNAs that regulate gene expression post‐transcriptionally. By repressing the translation and promoting the degradation of target mRNAs, miRNAs may reduce the cell‐to‐cell variability in protein expression, induce correlations between target expression levels, and provide a layer through which targets can influence each other's expression as “competing RNAs” (ceRNAs). However, experimental evidence for these behaviors is limited. Combining mathematical modeling with RNA sequencing of individual human embryonic kidney cells in which the expression of two distinct miRNAs was induced over a wide range, we have inferred parameters describing the response of hundreds of miRNA targets to miRNA induction. Individual targets have widely different response dynamics, and only a small proportion of predicted targets exhibit high sensitivity to miRNA induction. Our data reveal for the first time the response parameters of the entire network of endogenous miRNA targets to miRNA induction, demonstrating that miRNAs correlate target expression and at the same time increase the variability in expression of individual targets across cells. The approach is generalizable to other miRNAs and post‐transcriptional regulators to improve the understanding of gene expression dynamics in individual cell types.

Abstract Background Large sequence datasets are difficult to visualize and handle. Additionally, they often do not represent a random subset of the natural diversity, but the result of uncoordinated and convenience sampling. Consequently, they can suffer from redundancy and sampling biases. Results Here we present Treemmer, a simple tool to evaluate the redundancy of phylogenetic trees and reduce their complexity by eliminating leaves that contribute the least to the tree diversity. Conclusions Treemmer can reduce the size of datasets with different phylogenetic structures and levels of redundancy while maintaining a sub-sample that is representative of the original diversity. Additionally, it is possible to fine-tune the behavior of Treemmer including any kind of meta-information, making Treemmer particularly useful for empirical studies.

The Lineage 2–Beijing (L2–Beijing) sub-lineage of Mycobacterium tuberculosis has received much attention due to its high virulence, fast disease progression, and association with antibiotic resistance. Despite several reports of the recent emergence of L2–Beijing in Africa, no study has investigated the evolutionary history of this sub-lineage on the continent. In this study, we used whole genome sequences of 781 L2 clinical strains from 14 geographical regions globally distributed to investigate the origins and onward spread of this lineage in Africa. Our results reveal multiple introductions of L2–Beijing into Africa linked to independent bacterial populations from East- and Southeast Asia. Bayesian analyses further indicate that these introductions occurred during the past 300 years, with most of these events pre-dating the antibiotic era. Hence, the success of L2–Beijing in Africa is most likely due to its hypervirulence and high transmissibility rather than drug resistance.

The Mycobacterium tuberculosis complex (MTBC) causes tuberculosis (TB) in humans and various other mammals. The human-adapted members of the MTBC comprise seven phylogenetic lineages that differ in their geographical distribution. There is growing evidence that this phylogeographic diversity modulates the outcome of TB infection and disease. For decades, TB research and development has focused on the two canonical MTBC laboratory strains H37Rv and Erdman, both of which belong to Lineage 4. Relying on only a few laboratory-adapted strains can be misleading as study results might not be directly transferrable to clinical settings where patients are infected with a diverse array of strains, including drug-resistant variants. Here, we argue for the need to expand TB research and development by incorporating the phylogenetic diversity of the MTBC. To facilitate such work, we have assembled a group of 20 genetically well-characterized clinical strains representing the seven known human-adapted MTBC lineages. With the "MTBC clinical strains reference set" we aim to provide a standardized resource for the TB community. We hope it will enable more direct comparisons between studies that explore the physiology of MTBC beyond the laboratory strains used thus far. We anticipate that detailed phenotypic analyses of this reference strain set will increase our understanding of TB biology and assist in the development of new control tools that are broadly effective.

Since the emergence of coronavirus disease (COVID-19) in late 2019, domestic cats have been demonstrated to be susceptible to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) under natural and experimental conditions. As pet cats often live in very close contact with their owners, it is essential to investigate SARS-CoV-2 infections in cats in a One-Health context. This study reports the first SARS-CoV-2 infection in a cat in a COVID-19-affected household in Switzerland. The cat (Cat 1) demonstrated signs of an upper respiratory tract infection, including sneezing, inappetence, and apathy, while the cohabiting cat (Cat 2) remained asymptomatic. Nasal, oral, fecal, fur, and environmental swab samples were collected twice from both cats and analyzed by RT-qPCR for the presence of SARS-CoV-2 viral RNA. Both nasal swabs from Cat 1 tested positive. In addition, the first oral swab from Cat 2 and fur and bedding swabs from both cats were RT-qPCR positive. The fecal swabs tested negative. The infection of Cat 1 was confirmed by positive SARS-CoV-2 S1 receptor binding domain (RBD) antibody testing and neutralizing activity in a surrogate assay. The viral genome sequence from Cat 1, obtained by next generation sequencing, showed the closest relation to a human sequence from the B.1.1.39 lineage, with one single nucleotide polymorphism (SNP) difference. This study demonstrates not only SARS-CoV-2 infection of a cat from a COVID-19-affected household but also contamination of the cats’ fur and bed with viral RNA. Our results are important to create awareness that SARS-CoV-2 infected people should observe hygienic measures to avoid infection and contamination of animal cohabitants.

Polycomb group (PcG) proteins are major determinants of gene silencing and epigenetic memory in higher eukaryotes. Here, we systematically mapped the human PcG complexome using a robust affinity purification mass spectrometry approach. Our high-density protein interaction network uncovered a diverse range of PcG complexes. Moreover, our analysis identified PcG interactors linking them to the PcG system, thus providing insight into the molecular function of PcG complexes and mechanisms of recruitment to target genes. We identified two human PRC2 complexes and two PR-DUB deubiquitination complexes, which contain the O-linked N-acetylglucosamine transferase OGT1 and several transcription factors. Finally, genome-wide profiling of PR-DUB components indicated that the human PR-DUB and PRC1 complexes bind distinct sets of target genes, suggesting differential impact on cellular processes in mammals.

Tuberculosis (TB) affects humans and other animals and is caused by bacteria from the Mycobacterium tuberculosis complex (MTBC). Previous studies have shown that there are at least nine members of the MTBC infecting animals other than humans; these have also been referred to as ecotypes. However, the ecology and the evolution of these animal-adapted MTBC ecotypes are poorly understood. Here we screened 12,886 publicly available MTBC genomes and newly sequenced 17 animal-adapted MTBC strains, gathering a total of 529 genomes of animal-adapted MTBC strains. Phylogenomic and comparative analyses confirm that the animal-adapted MTBC members are paraphyletic with some members more closely related to the human-adapted Mycobacterium africanum Lineage 6 than to other animal-adapted strains. Furthermore, we identified four main animal-adapted MTBC clades that might correspond to four main host shifts; two of these clades are hypothesized to reflect independent cattle domestication events. Contrary to what would be expected from an obligate pathogen, MTBC nucleotide diversity was not positively correlated with host phylogenetic distances, suggesting that host tropism in the animal-adapted MTBC seems to be driven by contact rates and demographic aspects of the host population rather by than host relatedness. By combining phylogenomics with ecological data, we propose an evolutionary scenario in which the ancestor of Lineage 6 and all animal-adapted MTBC ecotypes was a generalist pathogen that subsequently adapted to different host species. This study provides a new phylogenetic framework to better understand the evolution of the different ecotypes of the MTBC and guide future work aimed at elucidating the molecular mechanisms underlying host range.

Tumor initiation is often linked to a loss of cellular identity. Transcriptional programs determining cellular identity are preserved by epigenetically-acting chromatin factors. Although such regulators are among the most frequently mutated genes in cancer, it is not well understood how an abnormal epigenetic condition contributes to tumor onset. In this work, we investigated the gene signature of tumors caused by disruption of the Drosophila epigenetic regulator, polyhomeotic (ph). In larval tissue ph mutant cells show a shift towards an embryonic-like signature. Using loss- and gain-of-function experiments we uncovered the embryonic transcription factor knirps (kni) as a new oncogene. The oncogenic potential of kni lies in its ability to activate JAK/STAT signaling and block differentiation. Conversely, tumor growth in ph mutant cells can be substantially reduced by overexpressing a differentiation factor. This demonstrates that epigenetically derailed tumor conditions can be reversed when targeting key players in the transcriptional network.

At every cell cycle, faithful inheritance of metazoan genomes requires the concerted activation of thousands of DNA replication origins. However, the genetic and chromatin features defining metazoan replication start sites remain largely unknown. Here, we delineate the origin repertoire of the Drosophila genome at high resolution. We address the role of origin-proximal G-quadruplexes and suggest that they transiently stall replication forks in vivo. We dissect the chromatin configuration of replication origins and identify a rich spatial organization of chromatin features at initiation sites. DNA shape and chromatin configurations, not strict sequence motifs, mark and predict origins in higher eukaryotes. We further examine the link between transcription and origin firing and reveal that modulation of origin activity across cell types is intimately linked to cell-type-specific transcriptional programs. Our study unravels conserved origin features and provides unique insights into the relationship among DNA topology, chromatin, transcription, and replication initiation across metazoa.

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.

Genomic imprinting is an epigenetic phenomenon leading to parent-of-origin specific differential expression of maternally and paternally inherited alleles. In plants, genomic imprinting has mainly been observed in the endosperm, an ephemeral triploid tissue derived after fertilization of the diploid central cell with a haploid sperm cell. In an effort to identify novel imprinted genes in Arabidopsis thaliana, we generated deep sequencing RNA profiles of F1 hybrid seeds derived after reciprocal crosses of Arabidopsis Col-0 and Bur-0 accessions. Using polymorphic sites to quantify allele-specific expression levels, we could identify more than 60 genes with potential parent-of-origin specific expression. By analyzing the distribution of DNA methylation and epigenetic marks established by Polycomb group (PcG) proteins using publicly available datasets, we suggest that for maternally expressed genes (MEGs) repression of the paternally inherited alleles largely depends on DNA methylation or PcG-mediated repression, whereas repression of the maternal alleles of paternally expressed genes (PEGs) predominantly depends on PcG proteins. While maternal alleles of MEGs are also targeted by PcG proteins, such targeting does not cause complete repression. Candidate MEGs and PEGs are enriched for cis-proximal transposons, suggesting that transposons might be a driving force for the evolution of imprinted genes in Arabidopsis. In addition, we find that MEGs and PEGs are significantly faster evolving when compared to other genes in the genome. In contrast to the predominant location of mammalian imprinted genes in clusters, cluster formation was only detected for few MEGs and PEGs, suggesting that clustering is not a major requirement for imprinted gene regulation in Arabidopsis.

Background: Lineage 1 (L1) and 3 (L3) are two lineages of the Mycobacterium tuberculosis complex (MTBC) causing tuberculosis (TB) in humans. L1 and L3 are prevalent around the rim of the Indian Ocean, the region that accounts for most of the world’s new TB cases. Despite their relevance for this region, L1 and L3 remain understudied. Methods: We analyzed 2,938 L1 and 2,030 L3 whole genome sequences originating from 69 countries. We reconstructed the evolutionary history of these two lineages and identified genes under positive selection. Results: We found a strongly asymmetric pattern of migration from South Asia toward neighboring regions, highlighting the historical role of South Asia in the dispersion of L1 and L3. Moreover, we found that several genes were under positive selection, including genes involved in virulence and resistance to antibiotics . For L1 we identified signatures of local adaptation at the esxH locus, a gene coding for a secreted effector that targets the human endosomal sorting complex, and is included in several vaccine candidates. Conclusions: Our study highlights the importance of genetic diversity in the MTBC, and sheds new light on two of the most important MTBC lineages affecting humans.