Michael Egholm

Genome sequencing in microfabricated high-density picolitre reactors

The proliferation of large-scale DNA-sequencing projects in recent years has driven a search for alternative methods to reduce time and cost. Here we describe a scalable, highly parallel sequencing system with raw throughput significantly greater than that of state-of-the-art capillary electrophoresis instruments. The apparatus uses a novel fibre-optic slide of individual wells and is able to sequence 25 million bases, at 99% or better accuracy, in one four-hour run. To achieve an approximately 100-fold increase in throughput over current Sanger sequencing technology, we have developed an emulsion method for DNA amplification and an instrument for sequencing by synthesis using a pyrosequencing protocol optimized for solid support and picolitre-scale volumes. Here we show the utility, throughput, accuracy and robustness of this system by shotgun sequencing and de novo assembly of the Mycoplasma genitalium genome with 96% coverage at 99.96% accuracy in one run of the machine.

A core gut microbiome in obese and lean twins

The human distal gut harbours a vast ensemble of microbes (the microbiota) that provide important metabolic capabilities, including the ability to extract energy from otherwise indigestible dietary polysaccharides. Studies of a few unrelated, healthy adults have revealed substantial diversity in their gut communities, as measured by sequencing 16S rRNA genes, yet how this diversity relates to function and to the rest of the genes in the collective genomes of the microbiota (the gut microbiome) remains obscure. Studies of lean and obese mice suggest that the gut microbiota affects energy balance by influencing the efficiency of calorie harvest from the diet, and how this harvested energy is used and stored. Here we characterize the faecal microbial communities of adult female monozygotic and dizygotic twin pairs concordant for leanness or obesity, and their mothers, to address how host genotype, environmental exposure and host adiposity influence the gut microbiome. Analysis of 154 individuals yielded 9,920 near full-length and 1,937,461 partial bacterial 16S rRNA sequences, plus 2.14 gigabases from their microbiomes. The results reveal that the human gut microbiome is shared among family members, but that each person’s gut microbial community varies in the specific bacterial lineages present, with a comparable degree of co-variation between adult monozygotic and dizygotic twin pairs. However, there was a wide array of shared microbial genes among sampled individuals, comprising an extensive, identifiable ‘core microbiome’ at the gene, rather than at the organismal lineage, level. Obesity is associated with phylum-level changes in the microbiota, reduced bacterial diversity and altered representation of bacterial genes and metabolic pathways. These results demonstrate that a diversity of organismal assemblages can nonetheless yield a core microbiome at a functional level, and that deviations from this core are associated with different physiological states (obese compared with lean).

A draft sequence of the Neandertal genome.

Kissing Cousins Neandertals, our closest relatives, ranged across Europe and Southwest Asia before their extinction approximately 30,000 years ago. Green et al. (p. 710) report a draft sequence of the Neandertal genome, created from three individuals, and compare it with genomes of five modern humans. The results suggest that ancient genomes of human relatives can be recovered with acceptably low contamination from modern human DNA. Because ancient DNA can be contaminated with microbial DNA, Burbano et al. (p. 723) developed a target sequence capture approach to obtain 14 kilobases of Neandertal DNA from a fairly poorly preserved sample with a high microbial load. A number of genomic regions and genes were revealed as candidates for positive selection early in modern human history. The genomic data suggest that Neandertals mixed with modern human ancestors some 120,000 years ago, leaving traces of Neandertal DNA in contemporary humans. Gene flow has occurred from Neandertals to humans of Eurasian descent, but not to Africans. Neandertals, the closest evolutionary relatives of present-day humans, lived in large parts of Europe and western Asia before disappearing 30,000 years ago. We present a draft sequence of the Neandertal genome composed of more than 4 billion nucleotides from three individuals. Comparisons of the Neandertal genome to the genomes of five present-day humans from different parts of the world identify a number of genomic regions that may have been affected by positive selection in ancestral modern humans, including genes involved in metabolism and in cognitive and skeletal development. We show that Neandertals shared more genetic variants with present-day humans in Eurasia than with present-day humans in sub-Saharan Africa, suggesting that gene flow from Neandertals into the ancestors of non-Africans occurred before the divergence of Eurasian groups from each other.

The complete genome of an individual by massively parallel DNA sequencing.

The association of genetic variation with disease and drug response, and improvements in nucleic acid technologies, have given great optimism for the impact of ‘genomic medicine’. However, the formidable size of the diploid human genome, approximately 6 gigabases, has prevented the routine application of sequencing methods to deciphering complete individual human genomes. To realize the full potential of genomics for human health, this limitation must be overcome. Here we report the DNA sequence of a diploid genome of a single individual, James D. Watson, sequenced to 7.4-fold redundancy in two months using massively parallel sequencing in picolitre-size reaction vessels. This sequence was completed in two months at approximately one-hundredth of the cost of traditional capillary electrophoresis methods. Comparison of the sequence to the reference genome led to the identification of 3.3 million single nucleotide polymorphisms, of which 10,654 cause amino-acid substitution within the coding sequence. In addition, we accurately identified small-scale (2–40,000 base pair (bp)) insertion and deletion polymorphism as well as copy number variation resulting in the large-scale gain and loss of chromosomal segments ranging from 26,000 to 1.5 million base pairs. Overall, these results agree well with recent results of sequencing of a single individual by traditional methods. However, in addition to being faster and significantly less expensive, this sequencing technology avoids the arbitrary loss of genomic sequences inherent in random shotgun sequencing by bacterial cloning because it amplifies DNA in a cell-free system. As a result, we further demonstrate the acquisition of novel human sequence, including novel genes not previously identified by traditional genomic sequencing. This is the first genome sequenced by next-generation technologies. Therefore it is a pilot for the future challenges of ‘personalized genome sequencing’.

Comprehensive human genome amplification using multiple displacement amplification

Fundamental to most genetic analysis is availability of genomic DNA of adequate quality and quantity. Because DNA yield from human samples is frequently limiting, much effort has been invested in developing methods for whole genome amplification (WGA) by random or degenerate oligonucleotide-primed PCR. However, existing WGA methods like degenerate oligonucleotide-primed PCR suffer from incomplete coverage and inadequate average DNA size. We describe a method, termed multiple displacement amplification (MDA), which provides a highly uniform representation across the genome. Amplification bias among eight chromosomal loci was less than 3-fold in contrast to 4–6 orders of magnitude for PCR-based WGA methods. Average product length was >10 kb. MDA is an isothermal, strand-displacing amplification yielding about 20–30 μg product from as few as 1–10 copies of human genomic DNA. Amplification can be carried out directly from biological samples including crude whole blood and tissue culture cells. MDA-amplified human DNA is useful for several common methods of genetic analysis, including genotyping of single nucleotide polymorphisms, chromosome painting, Southern blotting and restriction fragment length polymorphism analysis, subcloning, and DNA sequencing. MDA-based WGA is a simple and reliable method that could have significant implications for genetic studies, forensics, diagnostics, and long-term sample storage.

Peptide nucleic acids

Peptide nucleic acids (PNAs) are attractive, as compared to other classes of oligonucleotides that have been developed to date, in that they are relatively easy to synthesize and modify, hybridize to DNA and RNA with high affi nity and sequence selectivity, and are resistant to enzymatic degradation by proteases and nucleases; however, the downside is that they are only moderately soluble in aqueous solution. Herein we describe the protocols for synthesizing the second-generation γPNAs, both the monomers and oligomers, containing MiniPEG side chain with considerable improvements in water solubility, biocompatibility, and hybridization properties.

Rare Variants of IFIH1, a Gene Implicated in Antiviral Responses, Protect Against Type 1 Diabetes

Genome-wide association studies (GWASs) are regularly used to map genomic regions contributing to common human diseases, but they often do not identify the precise causative genes and sequence variants. To identify causative type 1 diabetes (T1D) variants, we resequenced exons and splice sites of 10 candidate genes in pools of DNA from 480 patients and 480 controls and tested their disease association in over 30,000 participants. We discovered four rare variants that lowered T1D risk independently of each other (odds ratio = 0.51 to 0.74; P = 1.3 × 10–3 to 2.1 × 10–16) in IFIH1 (interferon induced with helicase C domain 1), a gene located in a region previously associated with T1D by GWASs. These variants are predicted to alter the expression and structure of IFIH1 [MDA5 (melanoma differentiation-associated protein 5)], a cytoplasmic helicase that mediates induction of interferon response to viral RNA. This finding firmly establishes the role of IFIH1 in T1D and demonstrates that resequencing studies can pinpoint disease-causing genes in genomic regions initially identified by GWASs.

Analysis of one million base pairs of Neanderthal DNA

Neanderthals are the extinct hominid group most closely related to contemporary humans, so their genome offers a unique opportunity to identify genetic changes specific to anatomically fully modern humans. We have identified a 38,000-year-old Neanderthal fossil that is exceptionally free of contamination from modern human DNA. Direct high-throughput sequencing of a DNA extract from this fossil has thus far yielded over one million base pairs of hominoid nuclear DNA sequences. Comparison with the human and chimpanzee genomes reveals that modern human and Neanderthal DNA sequences diverged on average about 500,000 years ago. Existing technology and fossil resources are now sufficient to initiate a Neanderthal genome-sequencing effort.

Stability of peptide nucleic acids in human serum and cellular extracts

The stability of a new type of DNA mimic, peptide nucleic acid (PNA) in human blood serum, Eschericia coli and Micrococcus luteus extracts and nuclear and cytoplasmic extracts from mouse Ehrlich ascites tumor cells was investigated using HPLC analysis. Under conditions that caused complete cleavage of a control peptide, adrenocorticotropic hormone fragment 4-10, no significant degradation of the PNAs, H-T10-LysNH2 and H-TGTACGTCACAACTA-NH2, could be detected. Similarly, PNA H-T5-LysNH2 was found to resist attack by fungal proteinase K or porcine intestinal mucosa peptidase at concentrations exceeding those necessary to completely degrade a control peptide, H-Phe-Trp-Tyr-Cys-Phe-Trp-Tyr-Lys-Phe-Trp-Tyr-Lys-OH, by at least 1000- and 30-fold, respectively. Thus PNA is expected to have sufficient biostability to be used as a drug.

A New Arenavirus in a Cluster of Fatal Transplant-Associated Diseases

BACKGROUND Three patients who received visceral-organ transplants from a single donor on the same day died of a febrile illness 4 to 6 weeks after transplantation. Culture, polymerase-chain-reaction (PCR) and serologic assays, and oligonucleotide microarray analysis for a wide range of infectious agents were not informative. METHODS We evaluated RNA obtained from the liver and kidney transplant recipients. Unbiased high-throughput sequencing was used to identify microbial sequences not found by means of other methods. The specificity of sequences for a new candidate pathogen was confirmed by means of culture and by means of PCR, immunohistochemical, and serologic analyses. RESULTS High-throughput sequencing yielded 103,632 sequences, of which 14 represented an Old World arenavirus. Additional sequence analysis showed that this new arenavirus was related to lymphocytic choriomeningitis viruses. Specific PCR assays based on a unique sequence confirmed the presence of the virus in the kidneys, liver, blood, and cerebrospinal fluid of the recipients. Immunohistochemical analysis revealed arenavirus antigen in the liver and kidney transplants in the recipients. IgM and IgG antiviral antibodies were detected in the serum of the donor. Seroconversion was evident in serum specimens obtained from one recipient at two time points. CONCLUSIONS Unbiased high-throughput sequencing is a powerful tool for the discovery of pathogens. The use of this method during an outbreak of disease facilitated the identification of a new arenavirus transmitted through solid-organ transplantation.