Explore the vast range of titles available.

Background Next generation sequencing (NGS) technology has revolutionized genomic and genetic research. The pace of change in this area is rapid with three major new sequencing platforms having been released in 2011: Ion Torrent’s PGM, Pacific Biosciences’ RS and the Illumina MiSeq. Here we compare the results obtained with those platforms to the performance of the Illumina HiSeq, the current market leader. In order to compare these platforms, and get sufficient coverage depth to allow meaningful analysis, we have sequenced a set of 4 microbial genomes with mean GC content ranging from 19.3 to 67.7%. Together, these represent a comprehensive range of genome content. Here we report our analysis of that sequence data in terms of coverage distribution, bias, GC distribution, variant detection and accuracy. Results Sequence generated by Ion Torrent, MiSeq and Pacific Biosciences technologies displays near perfect coverage behaviour on GC-rich, neutral and moderately AT-rich genomes, but a profound bias was observed upon sequencing the extremely AT-rich genome of Plasmodium falciparum on the PGM, resulting in no coverage for approximately 30% of the genome. We analysed the ability to call variants from each platform and found that we could call slightly more variants from Ion Torrent data compared to MiSeq data, but at the expense of a higher false positive rate. Variant calling from Pacific Biosciences data was possible but higher coverage depth was required. Context specific errors were observed in both PGM and MiSeq data, but not in that from the Pacific Biosciences platform. Conclusions All three fast turnaround sequencers evaluated here were able to generate usable sequence. However there are key differences between the quality of that data and the applications it will support.

The assembly of DNA sequence data is undergoing a renaissance thanks to emerging technologies capable of producing reads tens of kilobases long. Assembling complete bacterial and small eukaryotic genomes is now possible, but the final step of circularizing sequences remains unsolved. Here we present Circlator, the first tool to automate assembly circularization and produce accurate linear representations of circular sequences. Using Pacific Biosciences and Oxford Nanopore data, Circlator correctly circularized 26 of 27 circularizable sequences, comprising 11 chromosomes and 12 plasmids from bacteria, the apicoplast and mitochondrion of Plasmodium falciparum and a human mitochondrion. Circlator is available at http://sanger-pathogens.github.io/circlator/ .

Immune-regulatory mechanisms of drug-free remission in rheumatoid arthritis (RA) are unknown. We hypothesized that synovial tissue macrophages (STM), which persist in remission, contribute to joint homeostasis. We used single-cell transcriptomics to profile 32,000 STMs and identified phenotypic changes in patients with early/active RA, treatment-refractory/active RA and RA in sustained remission. Each clinical state was characterized by different frequencies of nine discrete phenotypic clusters within four distinct STM subpopulations with diverse homeostatic, regulatory and inflammatory functions. This cellular atlas, combined with deep-phenotypic, spatial and functional analyses of synovial biopsy fluorescent activated cell sorted STMs, revealed two STM subpopulations (MerTK pos TREM2 high and MerTK pos LYVE1 pos ) with unique remission transcriptomic signatures enriched in negative regulators of inflammation. These STMs were potent producers of inflammation-resolving lipid mediators and induced the repair response of synovial fibroblasts in vitro. A low proportion of MerTK pos STMs in remission was associated with increased risk of disease flare after treatment cessation. Therapeutic modulation of MerTK pos STM subpopulations could therefore be a potential treatment strategy for RA. Multiple subpopulations of synovial tissue macrophages with varied transcriptional, phenotypic and functional features may contribute to disease flare and tissue repair in patients with active rheumatoid arthritis and patients in clinical remission.

Plasmodium falciparum erythrocyte membrane protein 1 ( Pf EMP1), encoded by the hypervariable var gene family, is central to malaria pathogenesis, influencing both disease severity and immune evasion. Classifying var genes into upstream groups (upsA, upsB, upsC, upsE) is important for understanding parasite biology and clinical outcomes, but remains challenging, especially with partial sequences, such as the DBLα tag or RNA-Seq assemblies. We developed upsML, a machine-learning-based classifier trained on 2,530 curated var genes, to accurately assign upstream groups based on sequence features from different partial gene regions. We compared seven methods, including support vector machines, random forests, XGBoost, and HMMER models. Several models in upsML achieve accuracies of 83% for DBLα-tag sequences and 92% for full-length Pf EMP1 sequences, thereby significantly outperforming existing tools. Additionally, we developed a model to distinguish internal from subtelomeric var genes, which we applied to a global collection of P. falciparum genomes, revealing a higher frequency of internal var genes in Asia. upsML is available at https://github.com/sii-scRNA-Seq/upsML , providing a robust and efficient resource for large-scale var gene analysis. It can classify var genes from 20 genomes in under one second.

Developmental steps in the trypanosome life-cycle involve transition between replicative and non-replicative forms specialised for survival in, and transmission between, mammalian and tsetse fly hosts. Here, using oligopeptide-induced differentiation in vitro, we model the progressive development of replicative ‘slender’ to transmissible ‘stumpy’ bloodstream form Trypanosoma brucei and capture the transcriptomes of 8,599 parasites using single cell transcriptomics (scRNA-seq). Using this framework, we detail the relative order of biological events during asynchronous development, profile dynamic gene expression patterns and identify putative regulators. We additionally map the cell cycle of proliferating parasites and position stumpy cell-cycle exit at early G1 before progression to a distinct G0 state. A null mutant for one transiently elevated developmental regulator, ZC3H20 is further analysed by scRNA-seq, identifying its point of failure in the developmental atlas. This approach provides a paradigm for the dissection of differentiation events in parasites, relevant to diverse transitions in pathogen biology. Trypanosoma brucei undergoes developmental steps during host infection. Here, using oligopeptide-induced differentiation in vitro, authors model replicative ‘slender’ to transmissible ‘stumpy’ bloodstream forms and identify developmental and cell cycle regulators by single cell transcriptomics.

Background Plasmodium parasites undergo several major developmental transitions during their complex lifecycle, which are enabled by precisely ordered gene expression programs. Transcriptomes from the 48-h blood stages of the major human malaria parasite Plasmodium falciparum have been described using cDNA microarrays and RNA-seq, but these assays have not always performed well within non-coding regions, where the AT-content is often 90–95%. Results We developed a directional, amplification-free RNA-seq protocol (DAFT-seq) to reduce bias against AT-rich cDNA, which we have applied to three strains of P. falciparum (3D7, HB3 and IT). While strain-specific differences were detected, overall there is strong conservation between the transcriptional profiles. For the 3D7 reference strain, transcription was detected from 89% of the genome, with over 78% of the genome transcribed into mRNAs. We also find that transcription from bidirectional promoters frequently results in non-coding, antisense transcripts. These datasets allowed us to refine the 5′ and 3′ untranslated regions (UTRs), which can be variable, long (> 1000 nt), and often overlap those of adjacent transcripts. Conclusions The approaches applied in this study allow a refined description of the transcriptional landscape of P. falciparum and demonstrate that very little of the densely packed P. falciparum genome is inactive or redundant. By capturing the 5′ and 3′ ends of mRNAs, we reveal both constant and dynamic use of transcriptional start sites across the intraerythrocytic developmental cycle that will be useful in guiding the definition of regulatory regions for use in future experimental gene expression studies.

The most polymorphic gene family in P. falciparum is the ∼60 var genes distributed across parasite chromosomes, both in the subtelomeres and in internal regions. They encode hypervariable surface proteins known as P. falciparum erythrocyte membrane protein 1 (PfEMP1) that are critical for pathogenesis and immune evasion in Plasmodium falciparum. How var gene sequence diversity is generated is not currently completely understood. To address this, we constructed large clone trees and performed whole genome sequence analysis to study the generation of novel var gene sequences in asexually replicating parasites. While single nucleotide polymorphisms (SNPs) were scattered across the genome, structural variants (deletions, duplications, translocations) were focused in and around var genes, with considerable variation in frequency between strains. Analysis of more than 100 recombination events involving var exon 1 revealed that the average nucleotide sequence identity of two recombining exons was only 63% (range: 52.7-72.4%) yet the crossovers were error-free and occurred in such a way that the resulting sequence was in frame and domain architecture was preserved. Var exon 1, which encodes the immunologically exposed part of the protein, recombined in up to 0.2% of infected erythrocytes in vitro per life cycle. The high rate of var exon 1 recombination indicates that millions of new antigenic structures could potentially be generated each day in a single infected individual. We propose a model whereby var gene sequence polymorphism is mainly generated during the asexual part of the life cycle.

Many evolutionarily distant pathogenic organisms have evolved similar survival strategies to evade the immune responses of their hosts. These include antigenic variation, through which an infecting organism prevents clearance by periodically altering the identity of proteins that are visible to the immune system of the host 1 . Antigenic variation requires large reservoirs of immunologically diverse antigen genes, which are often generated through homologous recombination, as well as mechanisms to ensure the expression of one or very few antigens at any given time. Both homologous recombination and gene expression are affected by three-dimensional genome architecture and local DNA accessibility 2 , 3 . Factors that link three-dimensional genome architecture, local chromatin conformation and antigenic variation have, to our knowledge, not yet been identified in any organism. One of the major obstacles to studying the role of genome architecture in antigenic variation has been the highly repetitive nature and heterozygosity of antigen-gene arrays, which has precluded complete genome assembly in many pathogens. Here we report the de novo haplotype-specific assembly and scaffolding of the long antigen-gene arrays of the model protozoan parasite Trypanosoma brucei , using long-read sequencing technology and conserved features of chromosome folding 4 . Genome-wide chromosome conformation capture (Hi-C) reveals a distinct partitioning of the genome, with antigen-encoding subtelomeric regions that are folded into distinct, highly compact compartments. In addition, we performed a range of analyses—Hi-C, fluorescence in situ hybridization, assays for transposase-accessible chromatin using sequencing and single-cell RNA sequencing—that showed that deletion of the histone variants H3.V and H4.V increases antigen-gene clustering, DNA accessibility across sites of antigen expression and switching of the expressed antigen isoform, via homologous recombination. Our analyses identify histone variants as a molecular link between global genome architecture, local chromatin conformation and antigenic variation. Long-read sequencing allows the assembly of antigen-gene arrays in Trypanosoma brucei and, coupled with deletion experiments, demonstrates that histone variants act as a molecular link between genome architecture, chromatin conformation and antigen variation.

Background Parasites use polymorphic gene families to evade the immune system or interact with the host. Assessing the diversity and expression of such gene families in pathogens can inform on the repertoire or host interaction phenotypes of clinical relevance. However, obtaining the sequences and quantifying their expression is a challenge. In Plasmodium falciparum , the highly polymorphic var genes encode the major virulence protein, PfEMP1, which bind a range of human receptors through varying combinations of DBL and CIDR domains. Here we present a tool, Varia, to predict near full-length gene sequences and domain compositions of query genes from database genes sharing short sequence tags. Varia generates output through two complementary pipelines. Varia_VIP returns all putative gene sequences and domain compositions of the query gene from any partial sequence provided, thereby enabling experimental validation of specific genes of interest and detailed assessment of their putative domain structure. Varia_GEM accommodates rapid profiling of  var  gene expression in complex patient samples from DBLα expression sequence tags (EST), by computing a sample overall transcript profile stratified by PfEMP1 domain types. Results Varia_VIP was tested querying sequence tags from all DBL domain types using different search criteria. On average 92% of query tags had one or more 99% identical database hits, resulting in the full-length query gene sequence being identified (> 99% identical DNA > 80% of query gene) among the five most prominent database hits, for ~ 33% of the query genes. Optimized Varia_GEM settings allowed correct prediction of > 90% of domains placed among the four most N-terminal domains, including the DBLα domain, and > 70% of C-terminal domains. With this accuracy, N-terminal domains could be predicted for > 80% of queries, whereas prediction rates of C-terminal domains dropped with the distance from the DBLα from 70 to 40%. Conclusion Prediction of var sequence and domain composition is possible from short sequence tags. Varia can be used to guide experimental validation of PfEMP1 sequences of interest and conduct high-throughput analysis of var type expression in patient samples.

Human African trypanosomiasis, or sleeping sickness, is caused by the protozoan parasite Trypanosoma brucei and induces profound reactivity of glial cells and neuroinflammation when the parasites colonise the central nervous system. However, the transcriptional and functional responses of the brain to chronic T. brucei infection remain poorly understood. By integrating single cell and spatial transcriptomics of the mouse brain, we identify that glial responses triggered by infection are readily detected in the proximity to the circumventricular organs, including the lateral and 3 rd ventricle. This coincides with the spatial localisation of both slender and stumpy forms of T. brucei . Furthermore, in silico predictions and functional validations led us to identify a previously unknown crosstalk between homeostatic microglia and Cd138 + plasma cells mediated by IL-10 and B cell activating factor (BAFF) signalling. This study provides important insights and resources to improve understanding of the molecular and cellular responses in the brain during infection with African trypanosomes. Detailed insight into how the brain responds to Trypanosoma brucei infection is lacking. Here, single cell and spatial transcriptomics are integrated to characterise this response, identifying a unique crosstalk between microglia and plasma cells.