SeqKit2: A Swiss army knife for sequence and alignment processing

Deep mutational scanning reveals transmembrane features governing surface expression of the B cell antigen receptor

B cells surveil the body for foreign matter using their surface-expressed B cell antigen receptor (BCR), a tetrameric complex comprising a membrane-tethered antibody (mIg) that binds antigens and a signaling dimer (CD79AB) that conveys this interaction to the B cell. Recent cryogenic electron microscopy (cryo-EM) structures of IgM and IgG isotype BCRs provide the first complete views of their architecture, revealing that the largest interaction surfaces between the mIg and CD79AB are in their transmembrane domains (TMDs). These structures support decades of biochemical work interrogating the requirements for assembly of a functional BCR and provide the basis for explaining the effects of mutations. Here we report a focused saturating mutagenesis to comprehensively characterize the nature of the interactions in the mIg TMD that are required for BCR surface expression. We examined the effects of 600 single-amino-acid changes simultaneously in a pooled competition assay and quantified their effects by next-generation sequencing. Our deep mutational scanning results reflect a feature-rich TMD sequence, with some positions completely intolerant to mutation and others requiring specific biochemical properties such as charge, polarity or hydrophobicity, emphasizing the high value of saturating mutagenesis over, for example, alanine scanning. The data agree closely with published mutagenesis and the cryo-EM structures, while also highlighting several positions and surfaces that have not previously been characterized or have effects that are difficult to rationalize purely based on structure. This unbiased and complete mutagenesis dataset serves as a reference and framework for informed hypothesis testing, design of therapeutics to regulate BCR surface expression and to annotate patient mutations.

Evidence of RNA polymerase III recruitment and transcription at protein-coding gene promoters

RNA polymerase (Pol) I, II, and III are most commonly described as having distinct roles in synthesizing ribosomal RNA (rRNA), messenger RNA (mRNA), and specific small noncoding (nc)RNAs, respectively. This delineation of transcriptional responsibilities is not definitive, however, as evidenced by instances of Pol II recruitment to genes conventionally transcribed by Pol III, including the co-transcription of RPPH1 - the catalytic RNA component of RNase P. A comprehensive understanding of the interplay between RNA polymerase complexes remains lacking, however, due to limited comparative analyses for all three enzymes. To address this gap, we applied a uniform framework for quantifying global Pol I, II, and III occupancies that integrates currently available human RNA polymerase ChIP-seq datasets. Occupancy maps are combined with a comprehensive multi-class promoter set that includes protein-coding genes, noncoding genes, and repetitive elements. While our genomic survey appropriately identifies recruitment of Pol I, II, and III to canonical target genes, we unexpectedly discover widespread recruitment of the Pol III machinery to promoters of specific protein-coding genes, supported by colocalization patterns observed for several Pol III-specific subunits. We show that Pol III-occupied Pol II promoters are enriched for small, nascent RNA reads terminating in a run of 4 Ts, a unique hallmark of Pol III transcription termination and evidence of active Pol III activity at these sites. Pol III disruption differentially modulates the expression of Pol III-occupied coding genes, which are functionally enriched for ribosomal proteins and genes broadly linked to unfavorable outcomes in cancer. Our map also identifies additional, currently unannotated genomic elements occupied by Pol III with clear signatures of nascent RNA species that are sensitive to disruption of La (SSB) - a Pol III-related RNA chaperone protein. These findings reshape our current understanding of the interplay between Pols II and III and identify potentially novel small ncRNAs with broad implications for gene regulatory paradigms and RNA biology.

How low can you go? Short-read polishing of Oxford Nanopore bacterial genome assemblies

It is now possible to assemble near-perfect bacterial genomes using Oxford Nanopore Technologies (ONT) long reads, but short-read polishing is usually required for perfection. However, the effect of short-read depth on polishing performance is not well understood. Here, we introduce Pypolca (with default and careful parameters) and Polypolish v0.6.0 (with a new careful parameter). We then show that: (1) all polishers other than Pypolca-careful, Polypolish-default and Polypolish-careful commonly introduce false-positive errors at low read depth; (2) most of the benefit of short-read polishing occurs by 25× depth; (3) Polypolish-careful almost never introduces false-positive errors at any depth; and (4) Pypolca-careful is the single most effective polisher. Overall, we recommend the following polishing strategies: Polypolish-careful alone when depth is very low (<5×), Polypolish-careful and Pypolca-careful when depth is low (5-25×), and Polypolish-default and Pypolca-careful when depth is sufficient (>25×).