Yoo-Bok Cho

Comprehensive Analysis to Improve the Validation Rate for Single Nucleotide Variants Detected by Next-Generation Sequencing

Next-generation sequencing (NGS) has enabled the high-throughput discovery of germline and somatic mutations. However, NGS-based variant detection is still prone to errors, resulting in inaccurate variant calls. Here, we categorized the variants detected by NGS according to total read depth (TD) and SNP quality (SNPQ), and performed Sanger sequencing with 348 selected non-synonymous single nucleotide variants (SNVs) for validation. Using the SAMtools and GATK algorithms, the validation rate was positively correlated with SNPQ but showed no correlation with TD. In addition, common variants called by both programs had a higher validation rate than caller-specific variants. We further examined several parameters to improve the validation rate, and found that strand bias (SB) was a key parameter. SB in NGS data showed a strong difference between the variants passing validation and those that failed validation, showing a validation rate of more than 92% (filtering cutoff value: alternate allele forward [AF]≥20 and AF<80 in SAMtools, SB<–10 in GATK). Moreover, the validation rate increased significantly (up to 97–99%) when the variant was filtered together with the suggested values of mapping quality (MQ), SNPQ and SB. This detailed and systematic study provides comprehensive recommendations for improving validation rates, saving time and lowering cost in NGS analyses.

The architecture of ArgR-DNA complexes at the genome-scale in Escherichia coli

DNA-binding motifs that are recognized by transcription factors (TFs) have been well studied; however, challenges remain in determining the in vivo architecture of TF-DNA complexes on a genome-scale. Here, we determined the in vivo architecture of Escherichia coli arginine repressor (ArgR)-DNA complexes using high-throughput sequencing of exonuclease-treated chromatin-immunoprecipitated DNA (ChIP-exo). The ChIP-exo has a unique peak-pair pattern indicating 5′ and 3′ ends of ArgR-binding region. We identified 62 ArgR-binding loci, which were classified into three groups, comprising single, double and triple peak-pairs. Each peak-pair has a unique 93 base pair (bp)-long (±2 bp) ArgR-binding sequence containing two ARG boxes (39 bp) and residual sequences. Moreover, the three ArgR-binding modes defined by the position of the two ARG boxes indicate that DNA bends centered between the pair of ARG boxes facilitate the non-specific contacts between ArgR subunits and the residual sequences. Additionally, our approach may also reveal other fundamental structural features of TF-DNA interactions that have implications for studying genome-scale transcriptional regulatory networks.

Current challenges in bacterial transcriptomics.

Over the past decade or so, dramatic developments in our ability to experimentally determine the content and function of genomes have taken place. In particular, next-generation sequencing technologies are now inspiring a new understanding of bacterial transcriptomes on a global scale. In bacterial cells, whole-transcriptome studies have not received attention, owing to the general view that bacterial genomes are simple. However, several recent RNA sequencing results are revealing unexpected levels of complexity in bacterial transcriptomes, indicating that the transcribed regions of genomes are much larger and complex than previously anticipated. In particular, these data show a wide array of small RNAs, antisense RNAs, and alternative transcripts. Here, we review how current transcriptomics are now revolutionizing our understanding of the complexity and regulation of bacterial transcriptomes.

Determination of single nucleotide variants in Escherichia coli DH5α by using short-read sequencing

Escherichia coli DH5α is a common laboratory strain that provides an important platform for routine use in cloning and synthetic biology applications. Many synthetic circuits have been constructed and successfully expressed in E. coli DH5α; however, its genome sequence has not been determined yet. Here, we determined E. coli DH5α genome sequence and identified genetic mutations that affect its phenotypic functions by using short-read sequencing. The sequencing results clearly described the genotypes of E. coli DH5α, which aid in further studies using the strain. Additionally, we observed 105 single nucleotide variants (SNVs), 83% of which were detected in protein-coding regions compared to the parental strain E. coli DH1. Interestingly, 23% of the protein-coding regions have mutations in their amino acid residues, whose biological functions were categorized into two-component systems, peptidoglycan biosynthesis and lipopolysaccharide biosynthesis. These results underscore the advantages of E. coli DH5α, which tolerates the components of transformation buffer and expresses foreign plasmids efficiently. Moreover, these SNVs were also observed in the commercially available strain. These data provide the genetic information of E. coli DH5α for its future application in metabolic engineering and synthetic biology.

Exploring the Functional Residues in a Flavin-Binding Fluorescent Protein Using Deep Mutational Scanning

Flavin mononucleotide (FMN)-based fluorescent proteins are versatile reporters that can monitor various cellular processes in both aerobic and anaerobic conditions. However, the understanding of the role of individual amino acid residues on the protein function has been limited and has restricted the development of better functional variants. Here we examine the functional amino acid residues of Escherichia coli flavin mononucleotide binding fluorescent protein (EcFbFP) using the application of high-throughput sequencing of functional variants, termed deep mutational scanning. The variants were classified into 329 function-retained (FR) and 259 function-loss (FL) mutations, and further the mutational enrichment in each amino acid residues was weighed to find the functionally important residues of EcFbFP. We show that the crucial amino acid residues of EcFbFP lie among the FMN-binding pocket, turns and loops of the protein where conformation changes occur, and spatially clustered residues near the E56-K97 salt bridges. In addition, the mutational sensitivity of the critical residues was confirmed by site-directed mutagenesis. The deep mutational scanning of EcFbFP has demonstrated important implications for constructing better functioning protein variants.