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.

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.

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.

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.

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.

Downregulation of KCNMB4 expression and changes in BK channel subtype in hippocampal granule neurons following seizure activity.

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.

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.

Local GABAergic signaling within sensory ganglia controls peripheral nociceptive transmission.

The integration of somatosensory information is generally assumed to be a function of the central nervous system (CNS). Here we describe fully functional GABAergic communication within rodent peripheral sensory ganglia and show that it can modulate transmission of pain-related signals from the peripheral sensory nerves to the CNS. We found that sensory neurons express major proteins necessary for GABA synthesis and release and that sensory neurons released GABA in response to depolarization. In vivo focal infusion of GABA or GABA reuptake inhibitor to sensory ganglia dramatically reduced acute peripherally induced nociception and alleviated neuropathic and inflammatory pain. In addition, focal application of GABA receptor antagonists to sensory ganglia triggered or exacerbated peripherally induced nociception. We also demonstrated that chemogenetic or optogenetic depolarization of GABAergic dorsal root ganglion neurons in vivo reduced acute and chronic peripherally induced nociception. Mechanistically, GABA depolarized the majority of sensory neuron somata, yet produced a net inhibitory effect on the nociceptive transmission due to the filtering effect at nociceptive fiber T-junctions. Our findings indicate that peripheral somatosensory ganglia represent a hitherto underappreciated site of somatosensory signal integration and offer a potential target for therapeutic intervention.

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.

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.

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.

Comprehensive variation discovery in single human genomes.

Complete knowledge of the genetic variation in individual human genomes is a crucial foundation for understanding the etiology of disease. Genetic variation is typically characterized by sequencing individual genomes and comparing reads to a reference. Existing methods do an excellent job of detecting variants in approximately 90% of the human genome; however, calling variants in the remaining 10% of the genome (largely low-complexity sequence and segmental duplications) is challenging. To improve variant calling, we developed a new algorithm, DISCOVAR, and examined its performance on improved, low-cost sequence data. Using a newly created reference set of variants from the finished sequence of 103 randomly chosen fosmids, we find that some standard variant call sets miss up to 25% of variants. We show that the combination of new methods and improved data increases sensitivity by several fold, with the greatest impact in challenging regions of the human genome.

Current understanding of iberiotoxin-resistant BK channels in the nervous system.

While most large-conductance, calcium-, and voltage-activated potassium channels (BK or Maxi-K type) are blocked by the scorpion venom iberiotoxin, the so-called "type II" subtype has the property of toxin resistance. This property is uniquely mediated by channel assembly with one member of the BK accessory ß subunit family, the neuron-enriched ß4 subunit. This review will focus on current understanding of iberiotoxin-resistant, ß4-containing BK channel properties and their function in the CNS. Studies have shown that ß4 dramatically promotes BK channel opening by shifting voltage sensor activation to more negative voltage ranges, but also slows activation to timescales that theoretically preclude BK ability to shape action potentials (APs). In addition, ß4 membrane trafficking is regulated through an endoplasmic retention signal and palmitoylation. More recently, the challenge has been to understand the functional role of the iberiotoxin-resistant BK subtype utilizing computational modeling of neurons and neurophysiological approaches. Utilizing iberiotoxin-resistance as a footprint for these channels, they have been identified in dentate gyrus granule neurons and in purkinje neurons of the cerebellum. In these neurons, the role of these channels is largely consistent with slow-gated channels that reduce excitability either through an interspike conductance, such as in purkinje neurons, or by replacing fast-gating BK channels that otherwise facilitate high frequency AP firing, such as in dentate gyrus neurons. They are also observed in presynaptic mossy fiber terminals of the dentate gyrus and posterior pituitary terminals. More recent studies suggest that ß4 subunits may also be expressed in some neurons lacking iberiotoxin-resistant BK channels, such as in CA3 hippocampus neurons. Ongoing research using novel, specific blockers and agonists of BK/ß4, and ß4 knockout mice, will continue to move the field forward in understanding the function of these channels.

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.

Control of somatic membrane potential in nociceptive neurons and its implications for peripheral nociceptive transmission.

Peripheral sensory ganglia contain somata of afferent fibres conveying somatosensory inputs to the central nervous system. Growing evidence suggests that the somatic/perisomatic region of sensory neurons can influence peripheral sensory transmission. Control of resting membrane potential (Erest) is an important mechanism regulating excitability, but surprisingly little is known about how Erest is regulated in sensory neuron somata or how changes in somatic/perisomatic Erest affect peripheral sensory transmission. We first evaluated the influence of several major ion channels on Erest in cultured small-diameter, mostly capsaicin-sensitive (presumed nociceptive) dorsal root ganglion (DRG) neurons. The strongest and most prevalent effect on Erest was achieved by modulating M channels, K2P and 4-aminopiridine-sensitive KV channels, while hyperpolarization-activated cyclic nucleotide-gated, voltage-gated Na(+), and T-type Ca(2+) channels to a lesser extent also contributed to Erest. Second, we investigated how varying somatic/perisomatic membrane potential, by manipulating ion channels of sensory neurons within the DRG, affected peripheral nociceptive transmission in vivo. Acute focal application of M or KATP channel enhancers or a hyperpolarization-activated cyclic nucleotide-gated channel blocker to L5 DRG in vivo significantly alleviated pain induced by hind paw injection of bradykinin. Finally, we show with computational modelling how somatic/perisomatic hyperpolarization, in concert with the low-pass filtering properties of the t-junction within the DRG, can interfere with action potential propagation. Our study deciphers a complement of ion channels that sets the somatic Erest of nociceptive neurons and provides strong evidence for a robust filtering role of the somatic and perisomatic compartments of peripheral nociceptive neuron.

Assemblathon 2: evaluating de novo methods of genome assembly in three vertebrate species.

BACKGROUND: The process of generating raw genome sequence data continues to become cheaper, faster, and more accurate. However, assembly of such data into high-quality, finished genome sequences remains challenging. Many genome assembly tools are available, but they differ greatly in terms of their performance (speed, scalability, hardware requirements, acceptance of newer read technologies) and in their final output (composition of assembled sequence). More importantly, it remains largely unclear how to best assess the quality of assembled genome sequences. The Assemblathon competitions are intended to assess current state-of-the-art methods in genome assembly. RESULTS: In Assemblathon 2, we provided a variety of sequence data to be assembled for three vertebrate species (a bird, a fish, and snake). This resulted in a total of 43 submitted assemblies from 21 participating teams. We evaluated these assemblies using a combination of optical map data, Fosmid sequences, and several statistical methods. From over 100 different metrics, we chose ten key measures by which to assess the overall quality of the assemblies. CONCLUSIONS: Many current genome assemblers produced useful assemblies, containing a significant representation of their genes and overall genome structure. However, the high degree of variability between the entries suggests that there is still much room for improvement in the field of genome assembly and that approaches which work well in assembling the genome of one species may not necessarily work well for another.

Characterizing and measuring bias in sequence data.

BACKGROUND: DNA sequencing technologies deviate from the ideal uniform distribution of reads. These biases impair scientific and medical applications. Accordingly, we have developed computational methods for discovering, describing and measuring bias. RESULTS: We applied these methods to the Illumina, Ion Torrent, Pacific Biosciences and Complete Genomics sequencing platforms, using data from human and from a set of microbes with diverse base compositions. As in previous work, library construction conditions significantly influence sequencing bias. Pacific Biosciences coverage levels are the least biased, followed by Illumina, although all technologies exhibit error-rate biases in high- and low-GC regions and at long homopolymer runs. The GC-rich regions prone to low coverage include a number of human promoters, so we therefore catalog 1,000 that were exceptionally resistant to sequencing. Our results indicate that combining data from two technologies can reduce coverage bias if the biases in the component technologies are complementary and of similar magnitude. Analysis of Illumina data representing 120-fold coverage of a well-studied human sample reveals that 0.20% of the autosomal genome was covered at less than 10% of the genome-wide average. Excluding locations that were similar to known bias motifs or likely due to sample-reference variations left only 0.045% of the autosomal genome with unexplained poor coverage. CONCLUSIONS: The assays presented in this paper provide a comprehensive view of sequencing bias, which can be used to drive laboratory improvements and to monitor production processes. Development guided by these assays should result in improved genome assemblies and better coverage of biologically important loci.

Mutations causing medullary cystic kidney disease type 1 lie in a large VNTR in MUC1 missed by massively parallel sequencing.

Although genetic lesions responsible for some mendelian disorders can be rapidly discovered through massively parallel sequencing of whole genomes or exomes, not all diseases readily yield to such efforts. We describe the illustrative case of the simple mendelian disorder medullary cystic kidney disease type 1 (MCKD1), mapped more than a decade ago to a 2-Mb region on chromosome 1. Ultimately, only by cloning, capillary sequencing and de novo assembly did we find that each of six families with MCKD1 harbors an equivalent but apparently independently arising mutation in sequence markedly under-represented in massively parallel sequencing data: the insertion of a single cytosine in one copy (but a different copy in each family) of the repeat unit comprising the extremely long (∼1.5-5 kb), GC-rich (>80%) coding variable-number tandem repeat (VNTR) sequence in the MUC1 gene encoding mucin 1. These results provide a cautionary tale about the challenges in identifying the genes responsible for mendelian, let alone more complex, disorders through massively parallel sequencing.

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.