A 2.91-billion base pair (bp) consensus sequence of the euchromatic portion of the human genome was generated by the whole-genome shotgun sequencing method. The 14.8-billion bp DNA sequence was generated over 9 months from 27,271,853 high-quality sequence reads (5.11-fold coverage of the genome) from both ends of plasmid clones made from the DNA of five individuals. Two assembly strategies—a whole-genome assembly and a regional chromosome assembly—were used, each combining sequence data from Celera and the publicly funded genome effort. The public data were shredded into 550-bp segments to create a 2.9-fold coverage of those genome regions that had been sequenced, without including biases inherent in the cloning and assembly procedure used by the publicly funded group. This brought the effective coverage in the assemblies to eightfold, reducing the number and size of gaps in the final assembly over what would be obtained with 5.11-fold coverage. The two assembly strategies yielded very similar results that largely agree with independent mapping data. The assemblies effectively cover the euchromatic regions of the human chromosomes. More than 90% of the genome is in scaffold assemblies of 100,000 bp or more, and 25% of the genome is in scaffolds of 10 million bp or larger. Analysis of the genome sequence revealed 26,588 protein-encoding transcripts for which there was strong corroborating evidence and an additional ∼12,000 computationally derived genes with mouse matches or other weak supporting evidence. Although gene-dense clusters are obvious, almost half the genes are dispersed in low G+C sequence separated by large tracts of apparently noncoding sequence. Only 1.1% of the genome is spanned by exons, whereas 24% is in introns, with 75% of the genome being intergenic DNA. Duplications of segmental blocks, ranging in size up to chromosomal lengths, are abundant throughout the genome and reveal a complex evolutionary history. Comparative genomic analysis indicates vertebrate expansions of genes associated with neuronal function, with tissue-specific developmental regulation, and with the hemostasis and immune systems. DNA sequence comparisons between the consensus sequence and publicly funded genome data provided locations of 2.1 million single-nucleotide polymorphisms (SNPs). A random pair of human haploid genomes differed at a rate of 1 bp per 1250 on average, but there was marked heterogeneity in the level of polymorphism across the genome. Less than 1% of all SNPs resulted in variation in proteins, but the task of determining which SNPs have functional consequences remains an open challenge.

The high degree of similarity between the mouse and human genomes is demonstrated through analysis of the sequence of mouse chromosome 16 (Mmu 16), which was obtained as part of a whole-genome shotgun assembly of the mouse genome. The mouse genome is about 10% smaller than the human genome, owing to a lower repetitive DNA content. Comparison of the structure and protein-coding potential of Mmu 16 with that of the homologous segments of the human genome identifies regions of conserved synteny with human chromosomes (Hsa) 3, 8, 12, 16, 21, and 22. Gene content and order are highly conserved between Mmu 16 and the syntenic blocks of the human genome. Of the 731 predicted genes on Mmu 16, 509 align with orthologs on the corresponding portions of the human genome, 44 are likely paralogous to these genes, and 164 genes have homologs elsewhere in the human genome; there are 14 genes for which we could find no human counterpart.

We report a whole-genome shotgun assembly (called WGSA) of the human genome generated at Celera in 2001. The Celeragenerated shotgun data set consisted of 27 million sequencing reads organized in pairs by virtue of end-sequencing 2-kbp, 10-kbp, and 50-kbp inserts from shotgun clone libraries.

Escherichia coli, Escherichia albertii, and Escherichia fergusonii are closely related bacteria that can cause illness in humans, such as bacteremia, urinary tract infections and diarrhea. Current identification strategies for these three species vary in complexity and typically rely on the use of multiple phenotypic and genetic tests. To facilitate their rapid identification, we developed a multiplex PCR assay targeting conserved, species-specific genes. We used the Daydreamer™ (Pattern Genomics, USA) software platform to concurrently analyze whole genome sequence assemblies (WGS) from 150 Enterobacteriaceae genomes (107 E. coli, 5 Shigella spp., 21 E. albertii, 12 E. fergusonii and 5 other species) and design primers for the following species-specific regions: a 212 bp region of the cyclic di-GMP regulator gene (cdgR, AW869_22935 from genome K-12 MG1655, CP014225) for E. coli/Shigella; a 393 bp region of the DNA-binding transcriptional activator of cysteine biosynthesis gene (EAKF1_ch4033 from genome KF1, CP007025) for E. albertii; and a 575 bp region of the palmitoleoyl-acyl carrier protein (ACP)-dependent acyltransferase (EFER_0790 from genome ATCC 35469, CU928158) for E. fergusonii. We incorporated the species-specific primers into a conventional multiplex PCR assay and assessed its performance with a collection of 97 Enterobacteriaceae strains. The assay was 100% sensitive and specific for detecting the expected species and offers a quick and accurate strategy for identifying E. coli, E. albertii, and E. fergusonii in either a single reaction or by in silico PCR with sequence assemblies.

Identifying and tracking microbial strains as microbiomes evolve are major challenges in the field of microbiome research. We utilized a new sequencing kit that combines DNA extraction with PCR amplification of a large region of the rRNA operon and downstream bioinformatic data analysis. Longitudinal microbiome samples of coadmitted twins from two different neonatal intensive care units (NICUs) were analyzed using an ∼2,500-base amplicon that spans the 16S and 23S rRNA genes and mapped to a new, custom 16S-23S rRNA database. Amplicon sequence variants (ASVs) inferred using DADA2 provided sufficient resolution for the differentiation of rRNA variants from closely related but not previously sequenced Klebsiella, Escherichia coli, and Enterobacter strains, among the first bacteria colonizing the gut of these infants after admission to the NICU. Distinct ASV groups (fingerprints) were monitored between coadmitted twins over time, demonstrating the potential to track the source and spread of both commensals and pathogens. The high-resolution taxonomy obtained from long amplicon sequencing enables the tracking of strains temporally and spatially as microbiomes are established in infants in the hospital environment. IMPORTANCE Achieving strain-level resolution is a major obstacle for source tracking and temporal studies of microbiomes. In this study, we describe a novel deep-sequencing approach that provides species- and strain-level resolution of the neonatal microbiome. Using Klebsiella, E. coli, and Enterobacter as examples, we could monitor their temporal dynamics after antibiotic treatment and in pairs of twins. The strain-level resolution, combined with the greater sequencing depth and decreased cost per read of PacBio Sequel 2, enables this advantageous source- and strain-tracking analysis method to be implemented widely across more complex microbiomes.