Search

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., Competing Interests: Declaration of interests D.B.J., B.A.A., M.J.T.S., and W.J.M. were employees and shareholders of 10× Genomics, Inc., at the time this work was performed. D.B.J., B.A.A., M.J.T.S., and W.J.M. are inventors on patents in relation to algorithms, therapeutic candidates, and other components of this work. W.J.M. and D.B.J. are founders, employees, and shareholders of Infinimmune. B.B. is a shareholder of Infinimmune and a member of its scientific advisory board., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

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., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Hao et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The vertebrate adaptive immune system modifies the genome of individual B cells to encode antibodies that bind particular antigens 1 . 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 antigens 2-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 clonotypes 23,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., (© 2022. The Author(s).)

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 Ca 2+ channel is suggested) that, in turn, boosts rebound depolarization at the soma. The model suggests that the fAHP may both reduce Ca 2+ 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.

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.

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., (© 2017 Weisenfeld et al.; Published by Cold Spring Harbor Laboratory Press.)

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., (Copyright © 2016 the American Physiological Society.)

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.

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., (Copyright © 2015 the American Physiological Society.)

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.

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., (Copyright © 2014 The Authors. Published by Elsevier B.V. All rights reserved.)

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.

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.

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.

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.

The discovery of a living coelacanth specimen in 1938 was remarkable, as this lineage of lobe-finned fish was thought to have become extinct 70 million years ago. The modern coelacanth looks remarkably similar to many of its ancient relatives, and its evolutionary proximity to our own fish ancestors provides a glimpse of the fish that first walked on land. Here we report the genome sequence of the African coelacanth, Latimeria chalumnae. Through a phylogenomic analysis, we conclude that the lungfish, and not the coelacanth, is the closest living relative of tetrapods. Coelacanth protein-coding genes are significantly more slowly evolving than those of tetrapods, unlike other genomic features. Analyses of changes in genes and regulatory elements during the vertebrate adaptation to land highlight genes involved in immunity, nitrogen excretion and the development of fins, tail, ear, eye, brain and olfaction. Functional assays of enhancers involved in the fin-to-limb transition and in the emergence of extra-embryonic tissues show the importance of the coelacanth genome as a blueprint for understanding tetrapod evolution.

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.

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.

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.

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.

Advances in next-generation sequencing (NGS) promise to facilitate diagnosis of inherited disorders. Although in research settings NGS has pinpointed causal alleles using segregation in large families, the key challenge for clinical diagnosis is application to single individuals. To explore its diagnostic use, we performed targeted NGS in 42 unrelated infants with clinical and biochemical evidence of mitochondrial oxidative phosphorylation disease. These devastating mitochondrial disorders are characterized by phenotypic and genetic heterogeneity, with more than 100 causal genes identified to date. We performed "MitoExome" sequencing of the mitochondrial DNA (mtDNA) and exons of ~1000 nuclear genes encoding mitochondrial proteins and prioritized rare mutations predicted to disrupt function. Because patients and healthy control individuals harbored a comparable number of such heterozygous alleles, we could not prioritize dominant-acting genes. However, patients showed a fivefold enrichment of genes with two such mutations that could underlie recessive disease. In total, 23 of 42 (55%) patients harbored such recessive genes or pathogenic mtDNA variants. Firm diagnoses were enabled in 10 patients (24%) who had mutations in genes previously linked to disease. Thirteen patients (31%) had mutations in nuclear genes not previously linked to disease. The pathogenicity of two such genes, NDUFB3 and AGK, was supported by complementation studies and evidence from multiple patients, respectively. The results underscore the potential and challenges of deploying NGS in clinical settings.

Low-cost short read sequencing technology has revolutionized genomics, though it is only just becoming practical for the high-quality de novo assembly of a novel large genome. We describe the Assemblathon 1 competition, which aimed to comprehensively assess the state of the art in de novo assembly methods when applied to current sequencing technologies. In a collaborative effort, teams were asked to assemble a simulated Illumina HiSeq data set of an unknown, simulated diploid genome. A total of 41 assemblies from 17 different groups were received. Novel haplotype aware assessments of coverage, contiguity, structure, base calling, and copy number were made. We establish that within this benchmark: (1) It is possible to assemble the genome to a high level of coverage and accuracy, and that (2) large differences exist between the assemblies, suggesting room for further improvements in current methods. The simulated benchmark, including the correct answer, the assemblies, and the code that was used to evaluate the assemblies is now public and freely available from http://www.assemblathon.org/.

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.

The BK channel is a Ca(2+) and voltage-gated conductance responsible for shaping action potential waveforms in many types of neurons. Type II BK channels are differentiated from type I channels by their pharmacology and slow gating kinetics. The β4 accessory subunit confers type II properties on BK α subunits. Empirically derived properties of BK channels, with and without the β4 accessory subunit, were obtained using a heterologous expression system under physiological ionic conditions. These data were then used to study how BK channels alone (type I) and with the accessory β4 subunit (type II) modulate action potential properties in biophysical neuron models. Overall, the models support the hypothesis that it is the slower kinetics provided by the β4 subunit that endows the BK channel with type II properties, which leads to broadening of action potentials and, secondarily, to greater recruitment of SK channels reducing neuronal excitability. Two regions of parameter space distinguished type II and type I effects; one where the range of BK-activating Ca(2+) was high (>20 μM) and the other where BK-activating Ca(2+) was low (∼0.4-1.2 μM). The latter required an elevated BK channel density, possibly beyond a likely physiological range. BK-mediated sharpening of the spike waveform associated with the lack of the β4 subunit was sensitive to the properties of voltage-gated Ca(2+) channels due to electrogenic effects on spike duration. We also found that depending on Ca(2+) dynamics, type II BK channels may have the ability to contribute to the medium AHP, a property not generally ascribed to BK channels, influencing the frequency-current relationship. Finally, we show how the broadening of action potentials conferred by type II BK channels can also indirectly increase the recruitment of SK-type channels decreasing the excitability of the neuron., (Copyright © 2011 IBRO. Published by Elsevier Ltd. All rights reserved.)

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.

Despite the ever-increasing output of Illumina sequencing data, loci with extreme base compositions are often under-represented or absent. To evaluate sources of base-composition bias, we traced genomic sequences ranging from 6% to 90% GC through the process by quantitative PCR. We identified PCR during library preparation as a principal source of bias and optimized the conditions. Our improved protocol significantly reduces amplification bias and minimizes the previously severe effects of PCR instrument and temperature ramp rate.

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.

We present an automated, high throughput library construction process for 454 technology. Sample handling errors and cross-contamination are minimized via end-to-end barcoding of plasticware, along with molecular DNA barcoding of constructs. Automation-friendly magnetic bead-based size selection and cleanup steps have been devised, eliminating major bottlenecks and significant sources of error. Using this methodology, one technician can create 96 sequence-ready 454 libraries in 2 days, a dramatic improvement over the standard method.

Pyramidal neurons in the CA3 region of the hippocampal formation integrate synaptic information arriving in the dendrites within discrete laminar regions. At potentials near or below the resting potential integration of synaptic signals is most affected by the passive properties of the cell and hyperpolarization-activated currents (I(h)). Here we focused specifically on a subset of neurons within the CA3b subregion of the rat hippocampus in order to better understand their membrane response within subthreshold voltage ranges. Using a combined experimental and computational approach we found that the passive properties of these neurons varied up to fivefold between cells. Likewise, there was a large variance in the expression of I(h) channels. However, the contribution of I(h) was minimal at resting potentials endowing the membrane with an apparent linear response to somatic current injection within +/-10 mV. Unlike in CA1 pyramidal neurons, however, I(h) activation was not potentiated in an activity-dependent manner. Computer modeling, based on a combination of voltage- and current-clamp data, suggested that an increasing density of these channels with distance from the soma, compared with a uniform distribution, would have no significant effect on the general properties of the cell because of their relatively lower expression. Nonetheless, temporal summation of excitatory inputs was affected by the presence of I(h) in the dendrites in a frequency- and distance-dependent fashion.

Targeting genomic loci by massively parallel sequencing requires new methods to enrich templates to be sequenced. We developed a capture method that uses biotinylated RNA 'baits' to fish targets out of a 'pond' of DNA fragments. The RNA is transcribed from PCR-amplified oligodeoxynucleotides originally synthesized on a microarray, generating sufficient bait for multiple captures at concentrations high enough to drive the hybridization. We tested this method with 170-mer baits that target >15,000 coding exons (2.5 Mb) and four regions (1.7 Mb total) using Illumina sequencing as read-out. About 90% of uniquely aligning bases fell on or near bait sequence; up to 50% lay on exons proper. The uniformity was such that approximately 60% of target bases in the exonic 'catch', and approximately 80% in the regional catch, had at least half the mean coverage. One lane of Illumina sequence was sufficient to call high-confidence genotypes for 89% of the targeted exon space.

Cancer results from somatic alterations in key genes, including point mutations, copy-number alterations and structural rearrangements. A powerful way to discover cancer-causing genes is to identify genomic regions that show recurrent copy-number alterations (gains and losses) in tumor genomes. Recent advances in sequencing technologies suggest that massively parallel sequencing may provide a feasible alternative to DNA microarrays for detecting copy-number alterations. Here we present: (i) a statistical analysis of the power to detect copy-number alterations of a given size; (ii) SegSeq, an algorithm to segment equal copy numbers from massively parallel sequence data; and (iii) analysis of experimental data from three matched pairs of tumor and normal cell lines. We show that a collection of approximately 14 million aligned sequence reads from human cell lines has comparable power to detect events as the current generation of DNA microarrays and has over twofold better precision for localizing breakpoints (typically, to within approximately 1 kilobase).

We demonstrate that genome sequences approaching finished quality can be generated from short paired reads. Using 36 base (fragment) and 26 base (jumping) reads from five microbial genomes of varied GC composition and sizes up to 40 Mb, ALLPATHS2 generated assemblies with long, accurate contigs and scaffolds. Velvet and EULER-SR were less accurate. For example, for Escherichia coli, the fraction of 10-kb stretches that were perfect was 99.8% (ALLPATHS2), 68.7% (Velvet), and 42.1% (EULER-SR).

We describe a new assembly algorithm, where a genome assembly with low sequence coverage, either throughout the genome or locally, due to cloning bias, is considerably improved through an assisting process via a related genome. We show that the information provided by aligning the whole-genome shotgun reads of the target against a reference genome can be used to substantially improve the quality of the resulting assembly.

Our variant ascertainment algorithm, VAAL, uses massively parallel DNA sequence data to identify differences between bacterial genomes with high sensitivity and specificity. VAAL detected approximately 98% of differences (including large insertion-deletions) between pairs of strains from three species while calling no false positives. VAAL also pinpointed a single mutation between Vibrio cholerae genomes, identifying an antibiotic's site of action by identifying sequence differences between drug-sensitive strains and drug-resistant derivatives.

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.

Promising new sequencing technologies, based on sequencing-by-synthesis (SBS), are starting to deliver large amounts of DNA sequence at very low cost. Polymorphism detection is a key application. We describe general methods for improved quality scores and accurate automated polymorphism detection, and apply them to data from the Roche (454) Genome Sequencer 20. We assess our methods using known-truth data sets, which is critical to the validity of the assessments. We developed informative, base-by-base error predictors for this sequencer and used a variant of the phred binning algorithm to combine them into a single empirically derived quality score. These quality scores are more useful than those produced by the system software: They both better predict actual error rates and identify many more high-quality bases. We developed a SNP detection method, with variants for low coverage, high coverage, and PCR amplicon applications, and evaluated it on known-truth data sets. We demonstrate good specificity in single reads, and excellent specificity (no false positives in 215 kb of genome) in high-coverage data.

New DNA sequencing technologies deliver data at dramatically lower costs but demand new analytical methods to take full advantage of the very short reads that they produce. We provide an initial, theoretical solution to the challenge of de novo assembly from whole-genome shotgun "microreads." For 11 genomes of sizes up to 39 Mb, we generated high-quality assemblies from 80x coverage by paired 30-base simulated reads modeled after real Illumina-Solexa reads. The bacterial genomes of Campylobacter jejuni and Escherichia coli assemble optimally, yielding single perfect contigs, and larger genomes yield assemblies that are highly connected and accurate. Assemblies are presented in a graph form that retains intrinsic ambiguities such as those arising from polymorphism, thereby providing information that has been absent from previous genome assemblies. For both C. jejuni and E. coli, this assembly graph is a single edge encompassing the entire genome. Larger genomes produce more complicated graphs, but the vast majority of the bases in their assemblies are present in long edges that are nearly always perfect. We describe a general method for genome assembly that can be applied to all types of DNA sequence data, not only short read data, but also conventional sequence reads.

It is thought that CA3 pyramidal neurons communicate mainly through bursts of spikes rather than so-called trains of regular firing action potentials. Reports of both burst firing and nonburst firing CA3 cells suggest that they may fire with more than one output pattern. With the use of whole-cell recording methods we studied the firing properties of rat hippocampal pyramidal neurons in vitro within the CA3b subregion and found three distinct types of firing patterns. Approximately 37% of cells were regular firing where spikes generated by minimal current injection (rheobase) were elicited with a short latency and with stronger current intensities trains of spikes exhibited spike frequency adaptation (SFA). Another 46% of neurons exhibited a delayed onset at rheobase with a weakly-adapting firing pattern upon stronger stimulation. The remaining 17% of cells showed a burst-firing pattern, though only elicited in response to strong current injection and spontaneous bursts were never observed. Control experiments indicated that the distinct firing patterns were not due to our particular slicing methods or recording techniques. Finally, computer modeling was used to identify how relative differences in K+ conductances, specifically K(C), K(M), and K(D), between cells contribute to the different characteristics of the three types of firing patterns observed experimentally., ((c) 2008 Wiley-Liss, Inc.)

Comparative analysis of multiple genomes in a phylogenetic framework dramatically improves the precision and sensitivity of evolutionary inference, producing more robust results than single-genome analyses can provide. The genomes of 12 Drosophila species, ten of which are presented here for the first time (sechellia, simulans, yakuba, erecta, ananassae, persimilis, willistoni, mojavensis, virilis and grimshawi), illustrate how rates and patterns of sequence divergence across taxa can illuminate evolutionary processes on a genomic scale. These genome sequences augment the formidable genetic tools that have made Drosophila melanogaster a pre-eminent model for animal genetics, and will further catalyse fundamental research on mechanisms of development, cell biology, genetics, disease, neurobiology, behaviour, physiology and evolution. Despite remarkable similarities among these Drosophila species, we identified many putatively non-neutral changes in protein-coding genes, non-coding RNA genes, and cis-regulatory regions. These may prove to underlie differences in the ecology and behaviour of these diverse species.

Purinergic P2Y receptors are one of four types of G(q/11)-coupled receptors in rat superior cervical ganglia (SCG) sympathetic neurons. In cultured SCG neurons, purinergic and bradykinin suppression of I(M) were similar in magnitude and somewhat less than that by muscarinic agonists. The effects of the P2Y receptor agonist UTP on neuronal excitability and discharge properties were studied. Under current clamp, UTP increased action potential (AP) firing in response to depolarizing current steps, depolarized the resting potential, decreased the threshold current required to fire an AP, and decreased spike-frequency adaptation. These effects were very similar to those resulting from bradykinin stimulation and not as profound as from muscarinic stimulation or full M-current blockade. We then examined the P2Y mechanism of action. Like bradykinin, but unlike muscarinic, purinergic stimulation induced rises in intracellular [Ca(2+)](i). Tests using expression of IP(3)"sponge" or IP(3) phosphatase constructs implicated IP(3) accumulation as necessary for purinergic suppression of I(M). Overexpression of wild-type or dominant-negative calmodulin (CaM) implicated Ca(2+)/CaM in the purinergic action. Both sets of results were similar to bradykinin, and opposite to muscarinic, suppression. We also examined modulation of Ca(2+) channels. As for bradykinin, purinergic stimulation did not suppress I(Ca), unless neuronal calcium sensor-1 (NCS-1) activity was blocked by a dominant-negative NCS-1 construct. Our results indicate that P2Y receptors modulate M-type channels in SCG cells via IP(3)-mediated [Ca(2+)](i) signals in concert with CaM and not by depletion of phosphatidylinositol-4, 5-biphosphate. We group purinergic P2Y and bradykinin B(2) receptors together as having a common mode of action.

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.

We report the generation and analysis of functional data from multiple, diverse experiments performed on a targeted 1% of the human genome as part of the pilot phase of the ENCODE Project. These data have been further integrated and augmented by a number of evolutionary and computational analyses. Together, our results advance the collective knowledge about human genome function in several major areas. First, our studies provide convincing evidence that the genome is pervasively transcribed, such that the majority of its bases can be found in primary transcripts, including non-protein-coding transcripts, and those that extensively overlap one another. Second, systematic examination of transcriptional regulation has yielded new understanding about transcription start sites, including their relationship to specific regulatory sequences and features of chromatin accessibility and histone modification. Third, a more sophisticated view of chromatin structure has emerged, including its inter-relationship with DNA replication and transcriptional regulation. Finally, integration of these new sources of information, in particular with respect to mammalian evolution based on inter- and intra-species sequence comparisons, has yielded new mechanistic and evolutionary insights concerning the functional landscape of the human genome. Together, these studies are defining a path for pursuit of a more comprehensive characterization of human genome function.

A key component of the ongoing ENCODE project involves rigorous comparative sequence analyses for the initially targeted 1% of the human genome. Here, we present orthologous sequence generation, alignment, and evolutionary constraint analyses of 23 mammalian species for all ENCODE targets. Alignments were generated using four different methods; comparisons of these methods reveal large-scale consistency but substantial differences in terms of small genomic rearrangements, sensitivity (sequence coverage), and specificity (alignment accuracy). We describe the quantitative and qualitative trade-offs concomitant with alignment method choice and the levels of technical error that need to be accounted for in applications that require multisequence alignments. Using the generated alignments, we identified constrained regions using three different methods. While the different constraint-detecting methods are in general agreement, there are important discrepancies relating to both the underlying alignments and the specific algorithms. However, by integrating the results across the alignments and constraint-detecting methods, we produced constraint annotations that were found to be robust based on multiple independent measures. Analyses of these annotations illustrate that most classes of experimentally annotated functional elements are enriched for constrained sequences; however, large portions of each class (with the exception of protein-coding sequences) do not overlap constrained regions. The latter elements might not be under primary sequence constraint, might not be constrained across all mammals, or might have expendable molecular functions. Conversely, 40% of the constrained sequences do not overlap any of the functional elements that have been experimentally identified. Together, these findings demonstrate and quantify how many genomic functional elements await basic molecular characterization.

Communication between the dentate gyrus (DG) and area CA3 of the hippocampus proper is transmitted via axons of granule cells--the mossy fiber (MF) pathway. In this review we discuss and compare the properties of transmitter release from the MFs onto pyramidal neurons and interneurons. An examination of the anatomical connectivity from DG to CA3 reveals a surprising interplay between excitation and inhibition for this circuit. In this respect it is particularly relevant that the major targets of the MFs are interneurons and that the consequence of MF input into CA3 may be inhibitory or excitatory, conditionally dependent on the frequency of input and modulatory regulation. This is further complicated by the properties of transmitter release from the MFs where a large number of co-localized transmitters, including GABAergic inhibitory transmitter release, and the effects of presynaptic modulation finely tune transmitter release. A picture emerges that extends beyond the hypothesis that the MFs are simply "detonators" of CA3 pyramidal neurons; the properties of synaptic information flow from the DG have more subtle and complex influences on the CA3 network.

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.

Ustilago maydis is a ubiquitous pathogen of maize and a well-established model organism for the study of plant-microbe interactions. This basidiomycete fungus does not use aggressive virulence strategies to kill its host. U. maydis belongs to the group of biotrophic parasites (the smuts) that depend on living tissue for proliferation and development. Here we report the genome sequence for a member of this economically important group of biotrophic fungi. The 20.5-million-base U. maydis genome assembly contains 6,902 predicted protein-encoding genes and lacks pathogenicity signatures found in the genomes of aggressive pathogenic fungi, for example a battery of cell-wall-degrading enzymes. However, we detected unexpected genomic features responsible for the pathogenicity of this organism. Specifically, we found 12 clusters of genes encoding small secreted proteins with unknown function. A significant fraction of these genes exists in small gene families. Expression analysis showed that most of the genes contained in these clusters are regulated together and induced in infected tissue. Deletion of individual clusters altered the virulence of U. maydis in five cases, ranging from a complete lack of symptoms to hypervirulence. Despite years of research into the mechanism of pathogenicity in U. maydis, no 'true' virulence factors had been previously identified. Thus, the discovery of the secreted protein gene clusters and the functional demonstration of their decisive role in the infection process illuminate previously unknown mechanisms of pathogenicity operating in biotrophic fungi. Genomic analysis is, similarly, likely to open up new avenues for the discovery of virulence determinants in other pathogens.