Michael E. Zwick

Microarray-based genomic selection for high-throughput resequencing

We developed a general method, microarray-based genomic selection (MGS), capable of selecting and enriching targeted sequences from complex eukaryotic genomes without the repeat blocking steps necessary for bacterial artificial chromosome (BAC)-based genomic selection. We demonstrate that large human genomic regions, on the order of hundreds of kilobases, can be enriched and resequenced with resequencing arrays. MGS, when combined with a next-generation resequencing technology, can enable large-scale resequencing in single-investigator laboratories.

Structure and Complexity of a Bacterial Transcriptome

Although gene expression has been studied in bacteria for decades, many aspects of the bacterial transcriptome remain poorly understood. Transcript structure, operon linkages, and information on absolute abundance all provide valuable insights into gene function and regulation, but none has ever been determined on a genome-wide scale for any bacterium. Indeed, these aspects of the prokaryotic transcriptome have been explored on a large scale in only a few instances, and consequently little is known about the absolute composition of the mRNA population within a bacterial cell. Here we report the use of a high-throughput sequencing-based approach in assembling the first comprehensive, single-nucleotide resolution view of a bacterial transcriptome. We sampled the Bacillus anthracis transcriptome under a variety of growth conditions and showed that the data provide an accurate and high-resolution map of transcript start sites and operon structure throughout the genome. Further, the sequence data identified previously nonannotated regions with significant transcriptional activity and enhanced the accuracy of existing genome annotations. Finally, our data provide estimates of absolute transcript abundance and suggest that there is significant transcriptional heterogeneity within a clonal, synchronized bacterial population. Overall, our results offer an unprecedented view of gene expression and regulation in a bacterial cell.

Haplotype Inference in Random Population Samples

Contemporary genotyping and sequencing methods do not provide information on linkage phase in diploid organisms. The application of statistical methods to infer and reconstruct linkage phase in samples of diploid sequences is a potentially time- and labor-saving method. The Stephens-Smith-Donnelly (SSD) algorithm is one such method, which incorporates concepts from population genetics theory in a Markov chain-Monte Carlo technique. We applied a modified SSD method, as well as the expectation-maximization and partition-ligation algorithms, to sequence data from eight loci spanning >1 Mb on the human X chromosome. We demonstrate that the accuracy of the modified SSD method is better than that of the other algorithms and is superior in terms of the number of sites that may be processed. Also, we find phase reconstructions by the modified SSD method to be highly accurate over regions with high linkage disequilibrium (LD). If only polymorphisms with a minor allele frequency >0.2 are analyzed and scored according to the fraction of neighbor relations correctly called, reconstructions are 95.2% accurate over entire 100-kb stretches and are 98.6% accurate within blocks of high LD.

Genomic characterization of the Bacillus cereus sensu lato species: Backdrop to the evolution of Bacillus anthracis

The key genes required for Bacillus anthracis to cause anthrax have been acquired recently by horizontal gene transfer. To understand the genetic background for the evolution of B. anthracis virulence, we obtained high-redundancy genome sequences of 45 strains of the Bacillus cereus sensu lato (s.l.) species that were chosen for their genetic diversity within the species based on the existing multilocus sequence typing scheme. From the resulting data, we called more than 324,000 new genes representing more than 12,333 new gene families for this group. The core genome size for the B. cereus s.l. group was ∼1750 genes, with another 2150 genes found in almost every genome constituting the extended core. There was a paucity of genes specific and conserved in any clade. We found no evidence of recent large-scale gene loss in B. anthracis or for unusual accumulation of nonsynonymous DNA substitutions in the chromosome; however, several B. cereus genomes isolated from soil and not previously associated with human disease were degraded to various degrees. Although B. anthracis has undergone an ecological shift within the species, its chromosome does not appear to be exceptional on a macroscopic scale compared with close relatives.

Microarray‐based mutation detection in the dystrophin gene

Duchenne and Becker muscular dystrophies (DMD and BMD) are X‐linked recessive neuromuscular disorders caused by mutations in the dystrophin gene affecting approximately 1 in 3,500 males. The human dystrophin gene spans>2,200 kb, or roughly 0.1% of the genome, and is composed of 79 exons. The mutational spectrum of disease‐causing alleles, including exonic copy number variations (CNVs), is complex. Deletions account for approximately 65% of DMD mutations and 85% of BMD mutations. Duplications occur in approximately 6 to 10% of males with either DMD or BMD. The remaining 30 to 35% of mutations consist of small deletions, insertions, point mutations, or splicing mutations, most of which introduce a premature stop codon. Laboratory analysis of dystrophin can be used to confirm a clinical diagnosis of DMD, characterize the type of dystrophin mutation, and perform prenatal testing and carrier testing for females. Current dystrophin diagnostic assays involve a variety of methodologies, including multiplex PCR, Southern blot analysis, multiplex ligation‐dependent probe amplification (MLPA), detection of virtually all mutations‐SSCP (DOVAM‐S), and single condition amplification/internal primer sequencing (SCAIP); however, these methods are time‐consuming, laborious, and do not accurately detect duplication mutations in the dystrophin gene. Furthermore, carrier testing in females is often difficult when a related affected male is unavailable. Here we describe the development, design, validation, and implementation of a high‐resolution comparative genomic hybridization (CGH) microarray‐based approach capable of accurately detecting both deletions and duplications in the dystrophin gene. This assay can be readily adopted by clinical molecular testing laboratories and represents a rapid, cost‐effective approach for screening a large gene, such as dystrophin. Hum Mutat 0, 1–9, 2008.

Genomic characterization of the Yersinia genus

BackgroundNew DNA sequencing technologies have enabled detailed comparative genomic analyses of entire genera of bacterial pathogens. Prior to this study, three species of the enterobacterial genus Yersinia that cause invasive human diseases (Yersinia pestis, Yersinia pseudotuberculosis, and Yersinia enterocolitica) had been sequenced. However, there were no genomic data on the Yersinia species with more limited virulence potential, frequently found in soil and water environments.ResultsWe used high-throughput sequencing-by-synthesis instruments to obtain 25- to 42-fold average redundancy, whole-genome shotgun data from the type strains of eight species: Y. aldovae, Y. bercovieri, Y. frederiksenii, Y. kristensenii, Y. intermedia, Y. mollaretii, Y. rohdei, and Y. ruckeri. The deepest branching species in the genus, Y. ruckeri, causative agent of red mouth disease in fish, has the smallest genome (3.7 Mb), although it shares the same core set of approximately 2,500 genes as the other members of the species, whose genomes range in size from 4.3 to 4.8 Mb. Yersinia genomes had a similar global partition of protein functions, as measured by the distribution of Cluster of Orthologous Groups families. Genome to genome variation in islands with genes encoding functions such as ureases, hydrogeneases and B-12 cofactor metabolite reactions may reflect adaptations to colonizing specific host habitats.ConclusionsRapid high-quality draft sequencing was used successfully to compare pathogenic and non-pathogenic members of the Yersinia genus. This work underscores the importance of the acquisition of horizontally transferred genes in the evolution of Y. pestis and points to virulence determinants that have been gained and lost on multiple occasions in the history of the genus.

Identification of Novel FMR1 Variants by Massively Parallel Sequencing in Developmentally Delayed Males

Fragile X syndrome (FXS), the most common inherited form of developmental delay, is typically caused by CGG‐repeat expansion in FMR1. However, little attention has been paid to sequence variants in FMR1. Through the use of pooled‐template massively parallel sequencing, we identified 130 novel FMR1 sequence variants in a population of 963 developmentally delayed males without CGG‐repeat expansion mutations. Among these, we identified a novel missense change, p.R138Q, which alters a conserved residue in the nuclear localization signal of FMRP. We have also identified three promoter mutations in this population, all of which significantly reduce in vitro levels of FMR1 transcription. Additionally, we identified 10 noncoding variants of possible functional significance in the introns and 3′‐untranslated region of FMR1, including two predicted splice site mutations. These findings greatly expand the catalog of known FMR1 sequence variants and suggest that FMR1 sequence variants may represent an important cause of developmental delay.

Microarray-based resequencing of multiple Bacillus anthracis isolates

We used custom-designed resequencing arrays to generate 3.1 Mb of genomic sequence from a panel of 56 Bacillus anthracis strains. Sequence quality was shown to be very high by replication (discrepancy rate of 7.4 × 10-7) and by comparison to independently generated shotgun sequence (discrepancy rate < 2.5 × 10-6). Population genomics studies of microbial pathogens using rapid resequencing technologies such as resequencing arrays are critical for recognizing newly emerging or genetically engineered strains.

Mutations in the human SC4MOL gene encoding a methyl sterol oxidase cause psoriasiform dermatitis, microcephaly, and developmental delay

Defects in cholesterol synthesis result in a wide variety of symptoms, from neonatal lethality to the relatively mild dysmorphic features and developmental delay found in individuals with Smith-Lemli-Opitz syndrome. We report here the identification of mutations in sterol-C4-methyl oxidase–like gene (SC4MOL) as the cause of an autosomal recessive syndrome in a human patient with psoriasiform dermatitis, arthralgias, congenital cataracts, microcephaly, and developmental delay. This gene encodes a sterol-C4-methyl oxidase (SMO), which catalyzes demethylation of C4-methylsterols in the cholesterol synthesis pathway. C4-Methylsterols are meiosis-activating sterols (MASs). They exist at high concentrations in the testis and ovary and play roles in meiosis activation. In this study, we found that an accumulation of MASs in the patient led to cell overproliferation in both skin and blood. SMO deficiency also substantially altered immunocyte phenotype and in vitro function. MASs serve as ligands for liver X receptors α and β(LXRα and LXRβ), which are important in regulating not only lipid transport in the epidermis, but also innate and adaptive immunity. Deficiency of SMO represents a biochemical defect in the cholesterol synthesis pathway, the clinical spectrum of which remains to be defined.

SeqAnt: A web service to rapidly identify and annotate DNA sequence variations

BackgroundThe enormous throughput and low cost of second-generation sequencing platforms now allow research and clinical geneticists to routinely perform single experiments that identify tens of thousands to millions of variant sites. Existing methods to annotate variant sites using information from publicly available databases via web browsers are too slow to be useful for the large sequencing datasets being routinely generated by geneticists. Because sequence annotation of variant sites is required before functional characterization can proceed, the lack of a high-throughput pipeline to efficiently annotate variant sites can act as a significant bottleneck in genetics research.ResultsSeqAnt (Seq uence An notator) is an open source web service and software package that rapidly annotates DNA sequence variants and identifies recessive or compound heterozygous loci in human, mouse, fly, and worm genome sequencing experiments. Variants are characterized with respect to their functional type, frequency, and evolutionary conservation. Annotated variants can be viewed on a web browser, downloaded in a tab-delimited text file, or directly uploaded in a BED format to the UCSC genome browser. To demonstrate the speed of SeqAnt, we annotated a series of publicly available datasets that ranged in size from 37 to 3,439,107 variant sites. The total time to completely annotate these data completely ranged from 0.17 seconds to 28 minutes 49.8 seconds.ConclusionSeqAnt is an open source web service and software package that overcomes a critical bottleneck facing research and clinical geneticists using second-generation sequencing platforms. SeqAnt will prove especially useful for those investigators who lack dedicated bioinformatics personnel or infrastructure in their laboratories.