Search

The vertebrate adaptive immune system modifies the genome of individual B cells to encode antibodies that bind particular antigens1. In most mammals, antibodies are composed of heavy and light chains that are generated sequentially by recombination of V, D (for heavy chains), J and C gene segments. Each chain contains three complementarity-determining regions (CDR1–CDR3), which contribute to antigen specificity. Certain heavy and light chains are preferred for particular antigens2–22. Here we consider pairs of B cells that share the same heavy chain V gene and CDRH3 amino acid sequence and were isolated from different donors, also known as public clonotypes23,24. We show that for naive antibodies (those not yet adapted to antigens), the probability that they use the same light chain V gene is around 10%, whereas for memory (functional) antibodies, it is around 80%, even if only one cell per clonotype is used. This property of functional antibodies is a phenomenon that we call light chain coherence. We also observe this phenomenon when similar heavy chains recur within a donor. Thus, although naive antibodies seem to recur by chance, the recurrence of functional antibodies reveals surprising constraint and determinism in the processes of V(D)J recombination and immune selection. For most functional antibodies, the heavy chain determines the light chain.

SummaryDevelopment of vaccines and therapeutics that are broadly effective against known and emergent coronaviruses is an urgent priority. Current strategies for developing pan-coronavirus countermeasures have largely focused on the receptor binding domain (RBD) and S2 regions of the coronavirus Spike protein; it has been unclear whether the N-terminal domain (NTD) is a viable target for universal vaccines and broadly neutralizing antibodies (Abs). Additionally, many RBD-targeting Abs have proven susceptible to viral escape. We screened the circulating B cell repertoires of COVID-19 survivors and vaccinees using multiplexed panels of uniquely barcoded antigens in a high-throughput single cell workflow to isolate over 9,000 SARS-CoV-2-specific monoclonal Abs (mAbs), providing an expansive view of the SARS-CoV-2-specific Ab repertoire. We observed many instances of clonal coalescence between individuals, suggesting that Ab responses frequently converge independently on similar genetic solutions. Among the recovered antibodies was TXG-0078, a public neutralizing mAb that binds the NTD supersite region of the coronavirus Spike protein and recognizes a diverse collection of alpha- and beta-coronaviruses. TXG-0078 achieves its exceptional binding breadth while utilizing the same VH1-24 variable gene signature and heavy chain-dominant binding pattern seen in other NTD supersite-specific neutralizing Abs with much narrower specificity. We also report the discovery of CC24.2, a pan-sarbecovirus neutralizing mAb that targets a novel RBD epitope and shows similar neutralization potency against all tested SARS-CoV-2 variants, including BQ.1.1 and XBB.1.5. A cocktail of TXG-0078 and CC24.2 provides protection againstin vivochallenge with SARS-CoV-2, suggesting potential future use in variant-resistant therapeutic Ab cocktails and as templates for pan-coronavirus vaccine design.

Accumulating observations suggest that peripheral somatosensory ganglia may regulate nociceptive transmission, yet direct evidence is sparse. Here, in experiments on rats and mice, we show that the peripheral afferent nociceptive information undergoes dynamic filtering within the dorsal root ganglion (DRG) and suggest that this filtering occurs at the axonal bifurcations (t-junctions). Using synchronous in vivo electrophysiological recordings from the peripheral and central processes of sensory neurons (in the spinal nerve and dorsal root), ganglionic transplantation of GABAergic progenitor cells, and optogenetics, we demonstrate existence of tonic and dynamic filtering of action potentials traveling through the DRG. Filtering induced by focal application of GABA or optogenetic GABA release from the DRG-transplanted GABAergic progenitor cells was specific to nociceptive fibers. Light-sheet imaging and computer modeling demonstrated that, compared to other somatosensory fiber types, nociceptors have shorter stem axons, making somatic control over t-junctional filtering more efficient. Optogenetically induced GABA release within DRG from the transplanted GABAergic cells enhanced filtering and alleviated hypersensitivity to noxious stimulation produced by chronic inflammation and neuropathic injury in vivo. These findings support “gating” of pain information by DRGs and suggest new therapeutic approaches for pain relief.

Half a billion years of evolutionary battle forged the vertebrate adaptive immune system, an astonishingly versatile factory for molecules that can adapt to arbitrary attacks. The history of an individual encounter is chronicled within a clonotype: the descendants of a single fully rearranged adaptive immune cell. For B cells, reading this immune history for an individual remains a fundamental challenge of modern immunology. Identification of such clonotypes is a magnificently challenging problem for three reasons:The cell history is inferred rather than directly observed: the only available data are the sequences of V(D)J molecules occurring in a sample of cells.Each immune receptor is a pair of V(D)J molecules. Identifying these pairs at scale is a technological challenge and cannot be done with perfect accuracy—real samples are mixtures of cells and fragments thereof.These molecules can be intensely mutatedduring the optimization of the response to particular antigens, blurring distinctions between kindred molecules.It is thus impossible to determine clonotypes exactly. All solutions to this problem make a trade-off between sensitivity and specificity; useful solutions must address actual artifacts found in real data.We presentenclone1, a system for computing approximate clonotypes from single cell data, and demonstrate its use and value with the 10x Genomics Immune Profiling Solution. To test it, we generate data for 1.6 million individual B cells, from four humans, including deliberately enriched memory cells, to tax the algorithm and provide a resource for the community. We analytically determine the specificity ofenclone’s clonotyping algorithm, showing that on this dataset the probability of co-clonotyping two unrelated B cells is around 10−9. We prove that using only heavy chains increases the error rate by two orders of magnitude.enclonecomprises a comprehensive toolkit for the analysis and display of immune receptor data. It is ultra-fast, easy to install, has public source code, comes with public data, and is documented atbit.ly/enclone. It has three “flavors” of use: (1) as a command-line tool run from a terminal window, that yields visual output; (2) as a command-line tool that yields parseable output that can be fed to other programs; and (3) as a graphical version (GUI).

Accumulating observations suggest that peripheral somatosensory ganglia may regulate pain transmission, yet direct evidence is sparse. Here we show that the peripheral afferent nociceptive information undergoes dynamic filtering within dorsal root ganglia (DRG) and suggest that this filtering occurs at the axonal bifurcations (t-junctions). Using simultaneous in vivo electrophysiological recordings from the peripheral (spinal nerve) and central (dorsal root) aspects of rodent spinal nerves, ganglionic transplantation of GABAergic progenitor cells, and optogenetics we demonstrate tonic and dynamic filtering of action potentials traveling through the DRG. Filtering induced by focal application of GABA or optogenetic GABA release from the DRG-transplanted GABAergic progenitor cells was specific to nociceptive fibers. Light-sheet imaging and computer modeling demonstrated that, compared to other somatosensory fiber types, nociceptors have shorter stem axons, making somatic control over t-junctional filtering more efficient. Optogenetically-induced GABA release within DRG enhanced filtering and reduced both acute and chronic inflammatory and neuropathic pain in vivo. These findings support the potential gating of pain information within the somatosensory system, and suggests new therapeutic approaches for pain relief. Significance Statement We show that the ascending somatosensory information undergoes dynamic filtering within the dorsal root ganglia (DRGs). Filtering (i.e. failure of action potentials to propagate through the axonal bifurcation or ‘t-junction’) is predominant in nociceptive nerve fibers and can be scaled up or down by the activity of GABAergic inhibitory system within the ganglia. Thus, spinal ganglion represents an intrinsic filtering device within the peripheral branch of somatosensory system. These findings may identify new strategies for pain relief.

We here pioneer a low-cost assembly strategy for 20 Heliconiini genomes to characterize the evolutionary history of the rapidly radiating genusHeliconius. A bifurcating tree provides a poor fit to the data, and we therefore explore a reticulate phylogeny forHeliconius. We probe the genomic architecture of gene flow, and develop a new method to distinguish incomplete lineage sorting from introgression. We find that most loci with non-canonical histories arose through introgression, and are strongly underrepresented in regions of low recombination and high gene density. This is expected if introgressed alleles are more likely to be purged in such regions due to tighter linkage with incompatibility loci. Finally, we identify a hitherto unrecognized inversion, and show it is a convergent structural rearrangement that captures a known color pattern switch locus within the genus. Our multi-genome assembly approach enables an improved understanding of adaptive radiation.

Continued advances in single cell gene detection sensitivity and the ability to analyze data quickly with biological context are critical for discoveries in therapeutic research. The second version of the Chromium Single Cell Immune Profiling Solution by 10x Genomics enables highly sensitive detection of gene expression, full-length paired T-cell alpha- and beta-chain and immunoglobulin sequences, T-cell antigen specificity, and cell surface protein expression from the same single cells, allowing a comprehensive view of the immune response at the cellular level. The new workflow provides a 45% increase in the number of genes detected per cell and up to a 25% increase in the cells detected with paired full-length V(D)J receptor sequences in melanoma tumor derived cells. Using a new version of Cell Ranger (v5.0), clonotypes were grouped from >30,000 cells, increasing the power to detect small clonotype expansions with fewer total cells. In addition, the new software enabled comparison of pre- and post-influenza vaccination B-cell receptor sequences from a single donor, identifying post-vaccination specific clonotypes. These technological and informatics advancements enhance a researcher’s ability to perform a broad characterization of immune cell populations at unprecedented throughput and resolution.

BK channels are large-conductance calcium- and voltage-activated potassium channels with diverse properties. Knockout of the accessory BK β4-subunit in hippocampus dentate gyrus granule neurons causes BK channels to change properties from slow-gated type II channels to fast-gated type I channels that sharpen the action potential, increase the fast afterhyperpolarization (fAHP) amplitude, and increase spike frequency. Here we studied the calcium channels that contribute to fast-gated BK channel activation and increased excitability of β4 knockout neurons. By using pharmacological blockers during current-clamp recording, we find that BK channel activation during the fAHP is dependent on ryanodine receptor activation. In contrast, L-type calcium channel blocker (nifedipine) affects the BK channel-dependent repolarization phase of the action potential but has no effect on the fAHP. Reducing BK channel activation during the repolarization phase with nifedipine, or during the fAHP with ryanodine, indicated that it is the BK-mediated increase of the fAHP that confers proexcitatory effects. The proexcitatory role of the fAHP was corroborated using dynamic current clamp. Increase or decrease of the fAHP amplitude during spiking revealed an inverse relationship between fAHP amplitude and interspike interval. Finally, we show that the seizure-prone ryanodine receptor gain-of-function (R2474S) knockin mice have an unaltered repolarization phase but larger fAHP and increased AP frequency compared with their control littermates. In summary, these results indicate that an important role of the β4-subunit is to reduce ryanodine receptor-BK channel functional coupling during the fAHP component of the action potential, thereby decreasing excitability of dentate gyrus neurons.

Profiling the complex interactions of immune infiltrate with tumor cells in a tumor microenvironment is critical for advancing our understanding of tumor biology for developing personalized cancer therapies. Using a droplet-based single cell RNA sequencing (scRNA-seq) platform, we profiled the transcriptome and immune repertoire of gastric, kidney, and lung cancer cells that are primarily immune cells from three different donors. scRNA-seq analysis also identified similar fractions of CD45+ immune cells, CD4+ and CD8+ T cells, CD19+ B cells, and myeloid cells compared to flow. Targeted scRNA-seq was also used to identify paired, full length B-cell (BCR) and T-cell (TCR) receptors. In the gastric adenocarcinoma tumor cells, gene expression analysis identified a large population of B cells, but with no clonal expansion. T cells expressing CD8A constituted about 10% of cells with the top clonotype being 1% of all clones. The clear cell RCC sample had a modest fraction of infiltrating T cells with the top clonotype representing 7.4% of all clones, and no B cell infiltrate. The NSCLC cells were mainly T and B cells but limited clonal expansion was observed; the top TCR clonotype was present on 2.2% of all T cells, no expansion was seen in any of the B cell clonotypes. These findings highlight the value of profiling of tumor immune cells holistically and not relying on the presence of B or T cells alone to understand the immune dynamics of the tumor microenvironment. The presence of tumor-infiltrating lymphocytes is associated with favorable clinical outcomes in some cancers, but understanding their cellular subtype and clonality by high resolution profiling is key in the development of immune-based cancer treatments

The gain of a neuron, the number and frequency of action potentials triggered in response to a given amount of depolarizing injection, is an important behavior underlying a neuron’s function. Variations in action potential waveform can influence neuronal discharges by the differential activation of voltage- and ion-gated channels long after the end of a spike. One component of the action potential waveform, the afterhyperpolarization (AHP), is generally considered an inhibitory mechanism for limiting firing rates. In dentate gyrus granule cells (DGCs) expressing fast-gated BK channels, large fast AHPs (fAHP) are paradoxically associated with increased gain. In this article, we describe a mechanism for this behavior using a computational model. Hyperpolarization provided by the fAHP enhances activation of a dendritic inward current (a T-type Ca2+channel is suggested) that, in turn, boosts rebound depolarization at the soma. The model suggests that the fAHP may both reduce Ca2+channel inactivation and, counterintuitively, enhance its activation. The magnitude of the rebound depolarization, in turn, determines the activation of a subsequent, slower inward current (a persistent Na+current is suggested) limiting the interspike interval. Simulations also show that the effect of AHP on gain is also effective for physiologically relevant stimulation; varying AHP amplitude affects interspike interval across a range of “noisy” stimulus frequency and amplitudes. The mechanism proposed suggests that small fAHPs in DGCs may contribute to their limited excitability.NEW & NOTEWORTHY The afterhyperpolarization (AHP) is canonically viewed as a major factor underlying the refractory period, serving to limit neuronal firing rate. We recently reported that enhancing the amplitude of the fast AHP (fAHP) in a relatively slowly firing neuron (vs. fast spiking neurons) expressing fast-gated BK channels augments neuronal excitability. In this computational study, we present a novel, quantitative hypothesis for how varying the amplitude of the fAHP can, paradoxically, influence a subsequent spike tens of milliseconds later.

Linked-Read sequencing technology has recently been employed successfully for de novo assembly of human genomes, however, the utility of this technology for complex plant genomes is unproven. We evaluated the technology for this purpose by sequencing the 3.5-gigabase (Gb) diploid pepper (Capsicum annuum) genome with a single Linked-Read library. Plant genomes, including pepper, are characterized by long, highly similar repetitive sequences. Accordingly, significant effort is used to ensure that the sequenced plant is highly homozygous and the resulting assembly is a haploid consensus. With a phased assembly approach, we targeted a heterozygous F1 derived from a wide cross to assess the ability to derive both haplotypes and characterize a pungency gene with a large insertion/deletion. The Supernova software generated a highly ordered, more contiguous sequence assembly than all currently available C. annuum reference genomes. Over 83% of the final assembly was anchored and oriented using four publicly available de novo linkage maps. A comparison of the annotation of conserved eukaryotic genes indicated the completeness of assembly. The validity of the phased assembly is further demonstrated with the complete recovery of both 2.5-Kb insertion/deletion haplotypes of the PUN1 locus in the F1 sample that represents pungent and nonpungent peppers, as well as nearly full recovery of the BUSCO2 gene set within each of the two haplotypes. The most contiguous pepper genome assembly to date has been generated which demonstrates that Linked-Read library technology provides a tool to de novo assemble complex highly repetitive heterozygous plant genomes. This technology can provide an opportunity to cost-effectively develop high-quality genome assemblies for other complex plants and compare structural and gene differences through accurate haplotype reconstruction.

A major challenge is to understand maladaptive changes in ion channels that sets neurons on a course towards epilepsy development. Voltage- and calcium-activated K+ (BK) channels contribute to early spike timing in neurons, and studies indicate that the BK channel plays a pathological role in increasing excitability early after a seizure. Here, we have investigated changes in BK channels and their accessory β4 subunit (KCNMB4) in dentate gyrus (DG) granule neurons of the hippocampus, key neurons that regulate excitability of the hippocampus circuit. Two days after pilocarpine-induced seizures, we found that the predominant effect is a downregulation of the β4 accessory subunit mRNA. Consistent with reduced expression, single channel recording and pharmacology indicate a switch in the subtype of channels expressed; from iberiotoxin-resistant, type II BK channels (BK α/β4) that have higher channel open probability and slow gating, to iberiotoxin-sensitive type I channels (BK α alone) with low open probability and faster gating. The switch to a majority of type I channel expression following seizure activity is correlated with a loss of BK channel function on spike threshold while maintaining the channel's contribution to increased early spike frequency. Using heterozygous β4 knockout mice, we find reduced expression is sufficient to increase seizure sensitivity. We conclude that seizure-induced downregulation of KCNMB4 is an activity dependent mechanism that increases the excitability of DG neurons. These novel findings indicate that BK channel subtypes are not only defined by cell-specific expression, but can also be plastic depending on the recent history of neuronal excitability.

BackgroundLinked-Read sequencing technology has recently been employed successfully for de novo assembly of multiple human genomes, however the utility of this technology for complex plant genomes is unproven. We evaluated the technology for this purpose by sequencing the 3.5 gigabase (Gb) diploid pepper (Capsicum annuum) genome with a single Linked-Read library. Plant genomes, including pepper, are characterized by long, highly similar repetitive sequences. Accordingly, significant effort is used to ensure the sequenced plant is highly homozygous and the resulting assembly is a haploid consensus. With a phased assembly approach, we targeted a heterozygous F1 derived from a wide cross to assess the ability to derive both haplotypes for a pungency gene characterized by a large insertion/deletion.ResultsThe Supernova software generated a highly ordered, more contiguous sequence assembly than all currently available C. annuum reference genomes. Eighty-four percent of the final assembly was anchored and oriented using four de novo linkage maps. A comparison of the annotation of conserved eukaryotic genes indicated the completeness of assembly. The validity of the phased assembly is further demonstrated with the complete recovery of both 2.5 kb insertion/deletion haplotypes of the PUN1 locus in the F1 sample that represents pungent and non-pungent peppers.ConclusionsThe most contiguous pepper genome assembly to date has been generated through this work which demonstrates that Linked-Read library technology provides a rapid tool to assemble de novo complex highly repetitive heterozygous plant genomes. This technology can provide an opportunity to cost-effectively develop high-quality reference genome assemblies for other complex plants and compare structural and gene differences through accurate haplotype reconstruction.

Determining the genome sequence of an organism is challenging, yet fundamental to understanding its biology. Over the past decade, thousands of human genomes have been sequenced, contributing deeply to biomedical research. In the vast majority of cases, these have been analyzed by aligning sequence reads to a single reference genome, biasing the resulting analyses and, in general, failing to capture sequences novel to a given genome.Some de novo assemblies have been constructed, free of reference bias, but nearly all were constructed by merging homologous loci into single ‘consensus’ sequences, generally absent from nature. These assemblies do not correctly represent the diploid biology of an individual. In exactly two cases, true diploid de novo assemblies have been made, at great expense. One was generated using Sanger sequencing and one using thousands of clone pools.Here we demonstrate a straightforward and low-cost method for creating true diploid de novo assemblies. We make a single library from ~1 ng of high molecular weight DNA, using the 10x Genomics microfluidic platform to partition the genome. We applied this technique to seven human samples, generating low-cost HiSeq X data, then assembled these using a new ‘pushbutton’ algorithm, Supernova. Each computation took two days on a single server. Each yielded contigs longer than 100 kb, phase blocks longer than 2.5 Mb, and scaffolds longer than 15 Mb. Our method provides a scalable capability for determining the actual diploid genome sequence in a sample, opening the door to new approaches in genomic biology and medicine.

Background De novo reference assemblies that are affordable, practical to produce, and of sufficient quality for most downstream applications, remain an unattained goal for many taxa. Insects, which may yield too little DNA from individual specimens for long-read sequencing library construction and often have highly heterozygous genomes, can be particularly hard to assemble using inexpensive short-read sequencing data. The large number of insect species with medical or economic importance makes this a critical problem to address. Results Using the assembler DISCOVAR de novo, we assembled the genome of the African malaria mosquito Anopheles arabiensis using 250 bp reads from a single library. The resulting assembly had a contig N50 of 22,433 bp, and recovered the gene set nearly as well as the ALLPATHS-LG AaraD1 An. arabiensis assembly produced with reads from three sequencing libraries and much greater resources. DISCOVAR de novo appeared to perform better than ALLPATHS-LG in regions of low complexity. Conclusions DISCOVAR de novo performed well assembling the genome of an insect of medical importance, using simpler sequencing input than previous anopheline assemblies. We have shown that this program is a viable tool for cost-effective assembly of a modestly-sized insect genome. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-2531-7) contains supplementary material, which is available to authorized users.

Unmyelinated C-fibers are a major type of sensory neurons conveying pain information. Action potential conduction is regulated by the bifurcation (T-junction) of sensory neuron axons within the dorsal root ganglia (DRG). Understanding how C-fiber signaling is influenced by the morphology of the T-junction and the local expression of ion channels is important for understanding pain signaling. In this study we used biophysical computer modeling to investigate the influence of axon morphology within the DRG and various membrane conductances on the reliability of spike propagation. As expected, calculated input impedance and the amplitude of propagating action potentials were both lowest at the T-junction. Propagation reliability for single spikes was highly sensitive to the diameter of the stem axon and the density of voltage-gated Na+ channels. A model containing only fast voltage-gated Na+ and delayed-rectifier K+ channels conducted trains of spikes up to frequencies of 110 Hz. The addition of slowly activating KCNQ channels (i.e., KV7 or M-channels) to the model reduced the following frequency to 30 Hz. Hyperpolarization produced by addition of a much slower conductance, such as a Ca2+-dependent K+ current, was needed to reduce the following frequency to 6 Hz. Attenuation of driving force due to ion accumulation or hyperpolarization produced by a Na+-K+ pump had no effect on following frequency but could influence the reliability of spike propagation mutually with the voltage shift generated by a Ca2+-dependent K+ current. These simulations suggest how specific ion channels within the DRG may contribute toward therapeutic treatments for chronic pain.

Determining the genome sequence of an organism is challenging, yet fundamental to understanding its biology. Over the past decade, thousands of human genomes have been sequenced, contributing deeply to biomedical research. In the vast majority of cases, these have been analyzed by aligning sequence reads to a single reference genome, biasing the resulting analyses, and in general, failing to capture sequences novel to a given genome. Some de novo assemblies have been constructed free of reference bias, but nearly all were constructed by merging homologous loci into single “consensus” sequences, generally absent from nature. These assemblies do not correctly represent the diploid biology of an individual. In exactly two cases, true diploid de novo assemblies have been made, at great expense. One was generated using Sanger sequencing, and one using thousands of clone pools. Here, we demonstrate a straightforward and low-cost method for creating true diploid de novo assemblies. We make a single library from ∼1 ng of high molecular weight DNA, using the 10x Genomics microfluidic platform to partition the genome. We applied this technique to seven human samples, generating low-cost HiSeq X data, then assembled these using a new “pushbutton” algorithm, Supernova. Each computation took 2 d on a single server. Each yielded contigs longer than 100 kb, phase blocks longer than 2.5 Mb, and scaffolds longer than 15 Mb. Our method provides a scalable capability for determining the actual diploid genome sequence in a sample, opening the door to new approaches in genomic biology and medicine.