Michael Zianni

Comparison of Commercially Available Target Enrichment Methods for Next-Generation Sequencing

Isolating high-priority segments of genomes greatly enhances the efficiency of next-generation sequencing (NGS) by allowing researchers to focus on their regions of interest. For the 2010-11 DNA Sequencing Research Group (DSRG) study, we compared outcomes from two leading companies, Agilent Technologies (Santa Clara, CA, USA) and Roche NimbleGen (Madison, WI, USA), which offer custom-targeted genomic enrichment methods. Both companies were provided with the same genomic sample and challenged to capture identical genomic locations for DNA NGS. The target region totaled 3.5 Mb and included 31 individual genes and a 2-Mb contiguous interval. Each company was asked to design its best assay, perform the capture in replicates, and return the captured material to the DSRG-participating laboratories. Sequencing was performed in two different laboratories on Genome Analyzer IIx systems (Illumina, San Diego, CA, USA). Sequencing data were analyzed for sensitivity, specificity, and coverage of the desired regions. The success of the enrichment was highly dependent on the design of the capture probes. Overall, coverage variability was higher for the Agilent samples. As variant discovery is the ultimate goal for a typical targeted sequencing project, we compared samples for their ability to sequence single-nucleotide polymorphisms (SNPs) as a test of the ability to capture both chromosomes from the sample. In the targeted regions, we detected 2546 SNPs with the NimbleGen samples and 2071 with Agilents. When limited to the regions that both companies included as baits, the number of SNPs was ∼1000 for each, with Agilent and NimbleGen finding a small number of unique SNPs not found by the other.

Borrelia burgdorferi oxidative stress regulator BosR directly represses lipoproteins primarily expressed in the tick during mammalian infection

Differential gene expression is a key strategy adopted by the Lyme disease spirochaete, Borrelia burgdorferi, for adaptation and survival in the mammalian host and the tick vector. Many B. burgdorferi surface lipoproteins fall into two distinct groups according to their expression patterns: one group primarily expressed in the tick and the other group primarily expressed in the mammal. Here, we show that the Fur homologue in this bacterium, also known as Borrelia oxidative stress regulator (BosR), is required for repression of outer surface protein A (OspA) and OspD in the mammal. Furthermore, BosR binds directly to sequences upstream of the ospAB operon and the ospD gene through recognition of palindromic motifs similar to those recognized by other Fur homologues but with a 1 bp variation in the spacer length. Putative BosR binding sites have been identified upstream of 156 B. burgdorferi genes. Some of these genes share the same expression pattern as ospA and ospD. Most notably, 12 (67%) of the 18 genes previously identified in a genome‐wide microarray study to be most significantly repressed in the mammal are among the putative BosR regulon. These data indicate that BosR may directly repress transcription of many genes that are downregulated in the mammal.

Molecular Characterization of Pantoea stewartii subsp. stewartii HrpY, a Conserved Response Regulator of the Hrp Type III Secretion System, and its Interaction with the hrpS Promoter

Pantoea stewartii subsp. stewartii is a bacterial pathogen of corn. Its pathogenicity depends on the translocation of effector proteins into host cells by the Hrp type III secretion system. We previously showed by genetic analysis that the HrpX sensor kinase and the HrpY response regulator are at the head of a complex cascade of regulators controlling hrp/hrc secretion and wts effector genes. This cascade also includes the HrpS response regulator and the HrpL alternative sigma factor. These regulators are shared among many important plant pathogens in the genera Pantoea, Erwinia, and Pseudomonas. In this study, we dissect the regulatory elements in the hrpS promoter region, using genetic and biochemical approaches, and show how it integrates various environmental signals, only some of which are dependent on phosphorylation of HrpY. Primer extension located the transcriptional start site of hrpS at a sigma70 promoter 601 bp upstream of the open reading frame. Electrophoretic mobility shift assays and DNase I footprinting analysis demonstrated that HrpY binds to conserved regulatory elements immediately adjacent to this promoter, and its binding affinity was increased by phosphorylation at D57. A consensus sequence for the two direct repeats bound by HrpY is proposed. Deletion analysis of the promoter region revealed that both the HrpY binding site and additional sequences farther upstream, including a putative integration host factor binding site, are required for hrpS expression. This finding suggests that other unknown regulatory proteins may act cooperatively with HrpY.

Sequencing of the Dutch elm disease fungus genome using the Roche/454 GS-FLX Titanium System in a comparison of multiple genomics core facilities.

As part of the DNA Sequencing Research Group of the Association of Biomolecular Resource Facilities, we have tested the reproducibility of the Roche/454 GS-FLX Titanium System at five core facilities. Experience with the Roche/454 system ranged from <10 to >340 sequencing runs performed. All participating sites were supplied with an aliquot of a common DNA preparation and were requested to conduct sequencing at a common loading condition. The evaluation of sequencing yield and accuracy metrics was assessed at a single site. The study was conducted using a laboratory strain of the Dutch elm disease fungus Ophiostoma novo-ulmi strain H327, an ascomycete, vegetatively haploid fungus with an estimated genome size of 30-50 Mb. We show that the Titanium System is reproducible, with some variation detected in loading conditions, sequencing yield, and homopolymer length accuracy. We demonstrate that reads shorter than the theoretical minimum length are of lower overall quality and not simply truncated reads. The O. novo-ulmi H327 genome assembly is 31.8 Mb and is comprised of eight chromosome-length linear scaffolds, a circular mitochondrial conti of 66.4 kb, and a putative 4.2-kb linear plasmid. We estimate that the nuclear genome encodes 8613 protein coding genes, and the mitochondrion encodes 15 genes and 26 tRNAs.

A variable DNA recognition site organization establishes the LiaR-mediated cell envelope stress response of enterococci to daptomycin

LiaR is a ‘master regulator’ of the cell envelope stress response in enterococci and many other Gram-positive organisms. Mutations to liaR can lead to antibiotic resistance to a variety of antibiotics including the cyclic lipopeptide daptomycin. LiaR is phosphorylated in response to membrane stress to regulate downstream target operons. Using DNA footprinting of the regions upstream of the liaXYZ and liaFSR operons we show that LiaR binds an extended stretch of DNA that extends beyond the proposed canonical consensus sequence suggesting a more complex level of regulatory control of target operons. We go on to determine the biochemical and structural basis for increased resistance to daptomycin by the adaptive mutation to LiaR (D191N) first identified from the pathogen Enterococcus faecalis S613. LiaRD191N increases oligomerization of LiaR to form a constitutively activated tetramer that has high affinity for DNA even in the absence of phosphorylation leading to increased resistance. Crystal structures of the LiaR DNA binding domain complexed to the putative consensus sequence as well as an adjoining secondary sequence show that upon binding, LiaR induces DNA bending that is consistent with increased recruitment of RNA polymerase to the transcription start site and upregulation of target operons.

Application of physicochemically modified silicon substrates as reverse phase protein microarrays

Physicochemically modified silicon substrates can provide a high quality alternative to nitrocellulose-coated glass slides for use in reverse-phase protein microarrays. Enhancement of protein microarray sensitivities is an important goal, especially because molecular targets within patient tissues exist in low abundance. The ideal array substrate has a high protein binding affinity and low intrinsic background signal. Silicon, which has low intrinsic autofluorescence, is being explored as a potential microarray surface. In a previous paper ( Nijdam , A. J. ; Cheng , M. M.-C. ; Fedele , R. ; Geho , D. H. ; Herrmann , P. ; Killian , K. ; Espina , V. ; Petricoin , E. F. ; Liotta , L. A. ; Ferrari , M. Physicochemically Modified Silicon as Substrate for Protein Microarrays . Biomaterials 2007 , 28 , 550 - 558 ), it is shown that physicochemical modification of silicon substrates increases the binding of protein to silicon to a level comparable with that of nitrocellulose. Here, we apply such substrates in a reverse-phase protein microarray setting in two model systems.

Further Unraveling the Regulatory Twist by Elucidating Metabolic Coinducer-Mediated CbbR-cbbI Promoter Interactions in Rhodopseudomonas palustris CGA010

The cbb(I) region of Rhodopseudomonas palustris (Rp. palustris) contains the cbbLS genes encoding form I ribulose-1,5-bisphosphate (RuBP) carboxylase oxygenase (RubisCO) along with a divergently transcribed regulator gene, cbbR. Juxtaposed between cbbR and cbbLS are the cbbRRS genes, encoding an unusual three-protein two-component (CbbRRS) system that modulates the ability of CbbR to influence cbbLS expression. The nature of the metabolic signals that Rp. palustris CbbR perceives to regulate cbbLS transcription is not known. Thus, in this study, the CbbR binding region was first mapped within the cbbLS promoter by the use of gel mobility shift assays and DNase I footprinting. In addition, potential metabolic coinducers (metabolites) were tested for their ability to alter the cbbLS promoter binding properties of CbbR. Gel mobility shift assays and surface plasmon resonance analyses together indicated that biosynthetic intermediates such as RuBP, ATP, fructose 1,6-bisphosphate, and NADPH enhanced DNA binding by CbbR. These coinducers did not yield identical CbbR-dependent DNase I footprints, indicating that the coinducers caused significant changes in DNA structure. These in vitro studies suggest that cellular signals such as fluctuating metabolite concentrations are perceived by and transduced to the cbbLS promoter via the master regulator CbbR.

Regulatory Twist and Synergistic Role of Metabolic Coinducer- and Response Regulator-Mediated CbbR-cbbI Interactions in Rhodopseudomonas palustris CGA010

Rhodopseudomonas palustris assimilates CO2 by the Calvin-Benson-Bassham (CBB) reductive pentose phosphate pathway. Most genes required for a functional CBB pathway are clustered into the cbbI and cbbII operons, with the cbbI operon subject to control by a LysR transcriptional activator, CbbR, encoded by cbbR, which is divergently transcribed from the cbbLS genes (encoding form I RubisCO) of the cbbI operon. Juxtaposed between the genes encoding CbbR and CbbLS are genes that encode a three-protein two-component system (CbbRRS system) that functions to modify the ability of CbbR to regulate cbbLS expression. Previous studies indicated that the response regulators, as well as various coinducers (effectors), specifically influence CbbR-promoter interactions. In the current study, it was shown via several experimental approaches that the response regulators and coinducers act synergistically on CbbR to influence cbbLS transcription. Synergistic effects on the formation of specific CbbR-DNA complexes were quantified using surface plasmon resonance (SPR) procedures. Gel mobility shift and DNA footprint analyses further indicated structural changes in the DNA arising from the presence of response regulators and coinducer molecules binding to CbbR. Based on previous studies, and especially emphasized by the current investigation, it is clear that protein complexes influence promoter activity and the cbbLS transcription machinery.

Rapid Discrimination between Anopheles gambiae s.s. and Anopheles arabiensis by High-Resolution Melt (HRM) Analysis

There is a need for more cost-effective options to more accurately discriminate among members of the Anopheles gambiae complex, particularly An. gambiae and Anopheles arabiensis. These species are morphologically indistinguishable in the adult stage, have overlapping distributions, but are behaviorally and ecologically different, yet both are efficient vectors of malaria in equatorial Africa. The method described here, High-Resolution Melt (HRM) analysis, takes advantage of minute differences in DNA melting characteristics, depending on the number of incongruent single nucleotide polymorphisms in an intragenic spacer region of the X-chromosome-based ribosomal DNA. The two species in question differ by an average of 13 single-nucleotide polymorphisms giving widely divergent melting curves. A real-time PCR system, Bio-Rad CFX96, was used in combination with a dsDNA-specific dye, EvaGreen, to detect and measure the melting properties of the amplicon generated from leg-extracted DNA of selected mosquitoes. Results with seven individuals from pure colonies of known species, as well as 10 field-captured individuals unambiguously identified by DNA sequencing, demonstrated that the method provided a high level of accuracy. The method was used to identify 86 field mosquitoes through the assignment of each to the two common clusters with a high degree of certainty. Each cluster was defined by individuals from pure colonies. HRM analysis is simpler to use than most other methods and provides comparable or more accurate discrimination between the two sibling species but requires a specialized melt-analysis instrument and software.

Altered residues in key proteins influence the expression and activity of the nitrogenase complex in an adaptive CO2 fixation-deficient mutant strain of Rhodobacter sphaeroides.

Previously, the RubisCO-compromised spontaneous adaptive Rhodobacter sphaeroides mutant, strain 16PHC, was shown to derepress the expression of genes that encode the nitrogenase complex under normal repressive conditions. As a result of this adaptation, the active nitrogenase complex restored redox balance, thus allowing strain 16PHC to grow under photoheterotrophic conditions in the absence of an exogenous electron acceptor. A combination of whole genome pyrosequencing and whole genome microarray analyses was employed to identify possible loci responsible for the observed phenotype. Mutations were found in two genes, glnA and nifA, whose products are involved in the regulatory cascade that controls nitrogenase complex gene expression. In addition, a nucleotide reversion within the nifK gene, which encodes a subunit of the nitrogenase complex, was also identified. Subsequent genetic, physiological and biochemical studies revealed alterations that led to derepression of the synthesis of an active nitrogenase complex in strain 16PHC.