Functional antibodies exhibit light chain coherence.

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.

Dorsal root ganglia control nociceptive input to the central nervous system.

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.

Direct determination of diploid genome sequences.

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.

The genomic substrate for adaptive radiation in African cichlid fish.

Cichlid fishes are famous for large, diverse and replicated adaptive radiations in the Great Lakes of East Africa. To understand the molecular mechanisms underlying cichlid phenotypic diversity, we sequenced the genomes and transcriptomes of five lineages of African cichlids: the Nile tilapia (Oreochromis niloticus), an ancestral lineage with low diversity; and four members of the East African lineage: Neolamprologus brichardi/pulcher (older radiation, Lake Tanganyika), Metriaclima zebra (recent radiation, Lake Malawi), Pundamilia nyererei (very recent radiation, Lake Victoria), and Astatotilapia burtoni (riverine species around Lake Tanganyika). We found an excess of gene duplications in the East African lineage compared to tilapia and other teleosts, an abundance of non-coding element divergence, accelerated coding sequence evolution, expression divergence associated with transposable element insertions, and regulation by novel microRNAs. In addition, we analysed sequence data from sixty individuals representing six closely related species from Lake Victoria, and show genome-wide diversifying selection on coding and regulatory variants, some of which were recruited from ancient polymorphisms. We conclude that a number of molecular mechanisms shaped East African cichlid genomes, and that amassing of standing variation during periods of relaxed purifying selection may have been important in facilitating subsequent evolutionary diversification.

A computational model for how the fast afterhyperpolarization paradoxically increases gain in regularly firing neurons.

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.

The genomic basis of adaptive evolution in threespine sticklebacks.

Marine stickleback fish have colonized and adapted to thousands of streams and lakes formed since the last ice age, providing an exceptional opportunity to characterize genomic mechanisms underlying repeated ecological adaptation in nature. Here we develop a high-quality reference genome assembly for threespine sticklebacks. By sequencing the genomes of twenty additional individuals from a global set of marine and freshwater populations, we identify a genome-wide set of loci that are consistently associated with marine-freshwater divergence. Our results indicate that reuse of globally shared standing genetic variation, including chromosomal inversions, has an important role in repeated evolution of distinct marine and freshwater sticklebacks, and in the maintenance of divergent ecotypes during early stages of reproductive isolation. Both coding and regulatory changes occur in the set of loci underlying marine-freshwater evolution, but regulatory changes appear to predominate in this well known example of repeated adaptive evolution in nature.

A high-resolution map of human evolutionary constraint using 29 mammals.

The comparison of related genomes has emerged as a powerful lens for genome interpretation. Here we report the sequencing and comparative analysis of 29 eutherian genomes. We confirm that at least 5.5% of the human genome has undergone purifying selection, and locate constrained elements covering ∼4.2% of the genome. We use evolutionary signatures and comparisons with experimental data sets to suggest candidate functions for ∼60% of constrained bases. These elements reveal a small number of new coding exons, candidate stop codon readthrough events and over 10,000 regions of overlapping synonymous constraint within protein-coding exons. We find 220 candidate RNA structural families, and nearly a million elements overlapping potential promoter, enhancer and insulator regions. We report specific amino acid residues that have undergone positive selection, 280,000 non-coding elements exapted from mobile elements and more than 1,000 primate- and human-accelerated elements. Overlap with disease-associated variants indicates that our findings will be relevant for studies of human biology, health and disease.

Evaluation of DISCOVAR de novo using a mosquito sample for cost-effective short-read genome assembly.

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.

Knockout of the BK ß4-subunit promotes a functional coupling of BK channels and ryanodine receptors that mediate a fAHP-induced increase in excitability.

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.

A genome-wide map of diversity in Plasmodium falciparum.

Genetic variation allows the malaria parasite Plasmodium falciparum to overcome chemotherapeutic agents, vaccines and vector control strategies and remain a leading cause of global morbidity and mortality. Here we describe an initial survey of genetic variation across the P. falciparum genome. We performed extensive sequencing of 16 geographically diverse parasites and identified 46,937 SNPs, demonstrating rich diversity among P. falciparum parasites (pi = 1.16 x 10(-3)) and strong correlation with gene function. We identified multiple regions with signatures of selective sweeps in drug-resistant parasites, including a previously unidentified 160-kb region with extremely low polymorphism in pyrimethamine-resistant parasites. We further characterized 54 worldwide isolates by genotyping SNPs across 20 genomic regions. These data begin to define population structure among African, Asian and American groups and illustrate the degree of linkage disequilibrium, which extends over relatively short distances in African parasites but over longer distances in Asian parasites. We provide an initial map of genetic diversity in P. falciparum and demonstrate its potential utility in identifying genes subject to recent natural selection and in understanding the population genetics of this parasite.

Spike propagation through the dorsal root ganglia in an unmyelinated sensory neuron: a modeling study.

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 Ca(2+)-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 Ca(2+)-dependent K(+) current. These simulations suggest how specific ion channels within the DRG may contribute toward therapeutic treatments for chronic pain.

Finished bacterial genomes from shotgun sequence data.

Exceptionally accurate genome reference sequences have proven to be of great value to microbial researchers. Thus, to date, about 1800 bacterial genome assemblies have been "finished" at great expense with the aid of manual laboratory and computational processes that typically iterate over a period of months or even years. By applying a new laboratory design and new assembly algorithm to 16 samples, we demonstrate that assemblies exceeding finished quality can be obtained from whole-genome shotgun data and automated computation. Cost and time requirements are thus dramatically reduced.

Paired-end sequencing of Fosmid libraries by Illumina.

Eliminating the bacterial cloning step has been a major factor in the vastly improved efficiency of massively parallel sequencing approaches. However, this also has made it a technical challenge to produce the modern equivalent of the Fosmid- or BAC-end sequences that were crucial for assembling and analyzing complex genomes during the Sanger-based sequencing era. To close this technology gap, we developed Fosill, a method for converting Fosmids to Illumina-compatible jumping libraries. We constructed Fosmid libraries in vectors with Illumina primer sequences and specific nicking sites flanking the cloning site. Our family of pFosill vectors allows multiplex Fosmid cloning of end-tagged genomic fragments without physical size selection and is compatible with standard and multiplex paired-end Illumina sequencing. To excise the bulk of each cloned insert, we introduced two nicks in the vector, translated them into the inserts, and cleaved them. Recircularization of the vector via coligation of insert termini followed by inverse PCR generates a jumping library for paired-end sequencing with 101-base reads. The yield of unique Fosmid-sized jumps is sufficiently high, and the background of short, incorrectly spaced and chimeric artifacts sufficiently low, to enable applications such as mapping of structural variation and scaffolding of de novo assemblies. We demonstrate the power of Fosill to map genome rearrangements in a cancer cell line and identified three fusion genes that were corroborated by RNA-seq data. Our Fosill-powered assembly of the mouse genome has an N50 scaffold length of 17.0 Mb, rivaling the connectivity (16.9 Mb) of the Sanger-sequencing based draft assembly.

Genome-scale DNA methylation maps of pluripotent and differentiated cells.

DNA methylation is essential for normal development and has been implicated in many pathologies including cancer. Our knowledge about the genome-wide distribution of DNA methylation, how it changes during cellular differentiation and how it relates to histone methylation and other chromatin modifications in mammals remains limited. Here we report the generation and analysis of genome-scale DNA methylation profiles at nucleotide resolution in mammalian cells. Using high-throughput reduced representation bisulphite sequencing and single-molecule-based sequencing, we generated DNA methylation maps covering most CpG islands, and a representative sampling of conserved non-coding elements, transposons and other genomic features, for mouse embryonic stem cells, embryonic-stem-cell-derived and primary neural cells, and eight other primary tissues. Several key findings emerge from the data. First, DNA methylation patterns are better correlated with histone methylation patterns than with the underlying genome sequence context. Second, methylation of CpGs are dynamic epigenetic marks that undergo extensive changes during cellular differentiation, particularly in regulatory regions outside of core promoters. Third, analysis of embryonic-stem-cell-derived and primary cells reveals that 'weak' CpG islands associated with a specific set of developmentally regulated genes undergo aberrant hypermethylation during extended proliferation in vitro, in a pattern reminiscent of that reported in some primary tumours. More generally, the results establish reduced representation bisulphite sequencing as a powerful technology for epigenetic profiling of cell populations relevant to developmental biology, cancer and regenerative medicine.

High-quality draft assemblies of mammalian genomes from massively parallel sequence data.

Massively parallel DNA sequencing technologies are revolutionizing genomics by making it possible to generate billions of relatively short (~100-base) sequence reads at very low cost. Whereas such data can be readily used for a wide range of biomedical applications, it has proven difficult to use them to generate high-quality de novo genome assemblies of large, repeat-rich vertebrate genomes. To date, the genome assemblies generated from such data have fallen far short of those obtained with the older (but much more expensive) capillary-based sequencing approach. Here, we report the development of an algorithm for genome assembly, ALLPATHS-LG, and its application to massively parallel DNA sequence data from the human and mouse genomes, generated on the Illumina platform. The resulting draft genome assemblies have good accuracy, short-range contiguity, long-range connectivity, and coverage of the genome. In particular, the base accuracy is high (≥99.95%) and the scaffold sizes (N50 size = 11.5 Mb for human and 7.2 Mb for mouse) approach those obtained with capillary-based sequencing. The combination of improved sequencing technology and improved computational methods should now make it possible to increase dramatically the de novo sequencing of large genomes. The ALLPATHS-LG program is available at http://www.broadinstitute.org/science/programs/genome-biology/crd.

Deep repertoire mining uncovers ultra-broad coronavirus neutralizing antibodies targeting multiple spike epitopes.

The development of vaccines and therapeutics that are broadly effective against known and emergent coronaviruses is an urgent priority. We screened the circulating B cell repertoires of COVID-19 survivors and vaccinees to isolate over 9,000 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-specific monoclonal antibodies (mAbs), providing an expansive view of the SARS-CoV-2-specific Ab repertoire. Among the recovered antibodies was TXG-0078, an N-terminal domain (NTD)-specific neutralizing mAb that recognizes diverse 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 CC24.2, a pan-sarbecovirus neutralizing antibody that targets a unique receptor-binding domain (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 shows protection in vivo, suggesting their potential use in variant-resistant therapeutic Ab cocktails and as templates for pan-coronavirus vaccine design.

Integrative analysis of the melanoma transcriptome.

Global studies of transcript structure and abundance in cancer cells enable the systematic discovery of aberrations that contribute to carcinogenesis, including gene fusions, alternative splice isoforms, and somatic mutations. We developed a systematic approach to characterize the spectrum of cancer-associated mRNA alterations through integration of transcriptomic and structural genomic data, and we applied this approach to generate new insights into melanoma biology. Using paired-end massively parallel sequencing of cDNA (RNA-seq) together with analyses of high-resolution chromosomal copy number data, we identified 11 novel melanoma gene fusions produced by underlying genomic rearrangements, as well as 12 novel readthrough transcripts. We mapped these chimeric transcripts to base-pair resolution and traced them to their genomic origins using matched chromosomal copy number information. We also used these data to discover and validate base-pair mutations that accumulated in these melanomas, revealing a surprisingly high rate of somatic mutation and lending support to the notion that point mutations constitute the major driver of melanoma progression. Taken together, these results may indicate new avenues for target discovery in melanoma, while also providing a template for large-scale transcriptome studies across many tumor types.

Genome-wide maps of chromatin state in pluripotent and lineage-committed cells.

We report the application of single-molecule-based sequencing technology for high-throughput profiling of histone modifications in mammalian cells. By obtaining over four billion bases of sequence from chromatin immunoprecipitated DNA, we generated genome-wide chromatin-state maps of mouse embryonic stem cells, neural progenitor cells and embryonic fibroblasts. We find that lysine 4 and lysine 27 trimethylation effectively discriminates genes that are expressed, poised for expression, or stably repressed, and therefore reflect cell state and lineage potential. Lysine 36 trimethylation marks primary coding and non-coding transcripts, facilitating gene annotation. Trimethylation of lysine 9 and lysine 20 is detected at satellite, telomeric and active long-terminal repeats, and can spread into proximal unique sequences. Lysine 4 and lysine 9 trimethylation marks imprinting control regions. Finally, we show that chromatin state can be read in an allele-specific manner by using single nucleotide polymorphisms. This study provides a framework for the application of comprehensive chromatin profiling towards characterization of diverse mammalian cell populations.

Deep repertoire mining uncovers ultra-broad coronavirus neutralizing antibodies targeting multiple spike epitopes.

Development 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 against in vivo challenge with SARS-CoV-2, suggesting potential future use in variant-resistant therapeutic Ab cocktails and as templates for pan-coronavirus vaccine design.