Reconstructing B-cell receptor sequences from short-read single-cell RNA sequencing with BRAPeS.

RNA sequencing of single B cells provides simultaneous measurements of the cell state and its antigen specificity as determined by the B-cell receptor (BCR). However, to uncover the latter, further reconstruction of the BCR sequence is needed. We present BRAPeS ("BCR Reconstruction Algorithm for Paired-end Single cells" ), an algorithm for reconstructing BCRs from short-read paired-end single-cell RNA sequencing. BRAPeS is accurate and achieves a high success rate even at very short (25 bp) read length, which can decrease the cost and increase the number of cells that can be analyzed compared with long reads. BRAPeS is publicly available at the following link: https://github.com/YosefLab/BRAPeS.

Comparative Analysis of Immune Cells Reveals a Conserved Regulatory Lexicon.

Most well-characterized enhancers are deeply conserved. In contrast, genome-wide comparative studies of steady-state systems showed that only a small fraction of active enhancers are conserved. To better understand conservation of enhancer activity, we used a comparative genomics approach that integrates temporal expression and epigenetic profiles in an innate immune system. We found that gene expression programs diverge among mildly induced genes, while being highly conserved for strongly induced genes. The fraction of conserved enhancers varies greatly across gene expression programs, with induced genes and early-response genes, in particular, being regulated by a higher fraction of conserved enhancers. Clustering of conserved accessible DNA sequences within enhancers resulted in over 60 sequence motifs including motifs for known factors, as well as many with unknown function. We further show that the number of instances of these motifs is a strong predictor of the responsiveness of a gene to pathogen detection.

Targeted reconstruction of T cell receptor sequence from single cell RNA-seq links CDR3 length to T cell differentiation state.

The T cell compartment must contain diversity in both T cell receptor (TCR) repertoire and cell state to provide effective immunity against pathogens. However, it remains unclear how differences in the TCR contribute to heterogeneity in T cell state. Single cell RNA-sequencing (scRNA-seq) can allow simultaneous measurement of TCR sequence and global transcriptional profile from single cells. However, current methods for TCR inference from scRNA-seq are limited in their sensitivity and require long sequencing reads, thus increasing the cost and decreasing the number of cells that can be feasibly analyzed. Here we present TRAPeS, a publicly available tool that can efficiently extract TCR sequence information from short-read scRNA-seq libraries. We apply it to investigate heterogeneity in the CD8+ T cell response in humans and mice, and show that it is accurate and more sensitive than existing approaches. Coupling TRAPeS with transcriptome analysis of CD8+ T cells specific for a single epitope from Yellow Fever Virus (YFV), we show that the recently described 'naive-like' memory population have significantly longer CDR3 regions and greater divergence from germline sequence than do effector-memory phenotype cells. This suggests that TCR usage is associated with the differentiation state of the CD8+ T cell response to YFV.

Defining the 5΄ and 3΄ landscape of the Drosophila transcriptome with Exo-seq and RNaseH-seq.

Cells regulate biological responses in part through changes in transcription start sites (TSS) or cleavage and polyadenylation sites (PAS). To fully understand gene regulatory networks, it is therefore critical to accurately annotate cell type-specific TSS and PAS. Here we present a simple and straightforward approach for genome-wide annotation of 5΄- and 3΄-RNA ends. Our approach reliably discerns bona fide PAS from false PAS that arise due to internal poly(A) tracts, a common problem with current PAS annotation methods. We applied our methodology to study the impact of temperature on the Drosophila melanogaster head transcriptome. We found hundreds of previously unidentified TSS and PAS which revealed two interesting phenomena: first, genes with multiple PASs tend to harbor a motif near the most proximal PAS, which likely represents a new cleavage and polyadenylation signal. Second, motif analysis of promoters of genes affected by temperature suggested that boundary element association factor of 32 kDa (BEAF-32) and DREF mediates a transcriptional program at warm temperatures, a result we validated in a fly line where beaf-32 is downregulated. These results demonstrate the utility of a high-throughput platform for complete experimental and computational analysis of mRNA-ends to improve gene annotation.

Clk post-transcriptional control denoises circadian transcription both temporally and spatially.

The transcription factor CLOCK (CLK) is essential for the development and maintenance of circadian rhythms in Drosophila. However, little is known about how CLK levels are controlled. Here we show that Clk mRNA is strongly regulated post-transcriptionally through its 3' UTR. Flies expressing Clk transgenes without normal 3' UTR exhibit variable CLK-driven transcription and circadian behaviour as well as ectopic expression of CLK-target genes in the brain. In these flies, the number of the key circadian neurons differs stochastically between individuals and within the two hemispheres of the same brain. Moreover, flies carrying Clk transgenes with deletions in the binding sites for the miRNA bantam have stochastic number of pacemaker neurons, suggesting that this miRNA mediates the deterministic expression of CLK. Overall our results demonstrate a key role of Clk post-transcriptional control in stabilizing circadian transcription, which is essential for proper development and maintenance of circadian rhythms in Drosophila.

Novel Observations From Next-Generation RNA Sequencing of Highly Purified Human Adult and Fetal Islet Cell Subsets.

Understanding distinct gene expression patterns of normal adult and developing fetal human pancreatic α- and ß-cells is crucial for developing stem cell therapies, islet regeneration strategies, and therapies designed to increase ß-cell function in patients with diabetes (type 1 or 2). Toward that end, we have developed methods to highly purify α-, ß-, and δ-cells from human fetal and adult pancreata by intracellular staining for the cell-specific hormone content, sorting the subpopulations by flow cytometry, and, using next-generation RNA sequencing, we report the detailed transcriptomes of fetal and adult α- and ß-cells. We observed that human islet composition was not influenced by age, sex, or BMI, and transcripts for inflammatory gene products were noted in fetal ß-cells. In addition, within highly purified adult glucagon-expressing α-cells, we observed surprisingly high insulin mRNA expression, but not insulin protein expression. This transcriptome analysis from highly purified islet α- and ß-cell subsets from fetal and adult pancreata offers clear implications for strategies that seek to increase insulin expression in type 1 and type 2 diabetes.

Differentiation of human parthenogenetic pluripotent stem cells reveals multiple tissue- and isoform-specific imprinted transcripts.

Parental imprinting results in monoallelic parent-of-origin-dependent gene expression. However, many imprinted genes identified by differential methylation do not exhibit complete monoallelic expression. Previous studies demonstrated complex tissue-dependent expression patterns for some imprinted genes. Still, the complete magnitude of this phenomenon remains largely unknown. By differentiating human parthenogenetic induced pluripotent stem cells into different cell types and combining DNA methylation with a 5' RNA sequencing methodology, we were able to identify tissue- and isoform-dependent imprinted genes in a genome-wide manner. We demonstrate that nearly half of all imprinted genes express both biallelic and monoallelic isoforms that are controlled by tissue-specific alternative promoters. This study provides a global analysis of tissue-specific imprinting in humans and suggests that alternative promoters are central in the regulation of imprinted genes.

BREX is a novel phage resistance system widespread in microbial genomes.

The perpetual arms race between bacteria and phage has resulted in the evolution of efficient resistance systems that protect bacteria from phage infection. Such systems, which include the CRISPR-Cas and restriction-modification systems, have proven to be invaluable in the biotechnology and dairy industries. Here, we report on a six-gene cassette in Bacillus cereus which, when integrated into the Bacillus subtilis genome, confers resistance to a broad range of phages, including both virulent and temperate ones. This cassette includes a putative Lon-like protease, an alkaline phosphatase domain protein, a putative RNA-binding protein, a DNA methylase, an ATPase-domain protein, and a protein of unknown function. We denote this novel defense system BREX (Bacteriophage Exclusion) and show that it allows phage adsorption but blocks phage DNA replication. Furthermore, our results suggest that methylation on non-palindromic TAGGAG motifs in the bacterial genome guides self/non-self discrimination and is essential for the defensive function of the BREX system. However, unlike restriction-modification systems, phage DNA does not appear to be cleaved or degraded by BREX, suggesting a novel mechanism of defense. Pan genomic analysis revealed that BREX and BREX-like systems, including the distantly related Pgl system described in Streptomyces coelicolor, are widely distributed in ~10% of all sequenced microbial genomes and can be divided into six coherent subtypes in which the gene composition and order is conserved. Finally, we detected a phage family that evades the BREX defense, implying that anti-BREX mechanisms may have evolved in some phages as part of their arms race with bacteria.