Hereditary cancer syndrome-associated pathogenic variants are common in patients with hematologic malignancies subsequent to primary solid cancer.

Background: As the number of long-term survivors of solid cancers keeps increasing, risk assessment of secondary hematologic malignancies is important for the prognosis of the patient. Germline genetic predisposition to secondary hematologic malignancy has been studied widely in myeloid neoplasms and rarely in lymphoid neoplasms. This study aimed to profile the mutational spectrums of patients with subsequent lymphoid tissue neoplasm to shed some light on the understudied area. Methods: In total, 39 patients who had primary solid cancer and subsequent hematologic malignancies were enrolled. We performed two next-generation sequencing (NGS) panel tests encompassing hereditary cancer predisposition genes and genes related to clonal hematopoiesis of indeterminate potential (CHIP). All statistical analyses were performed using R 3.5.1. Results: We found 8 of 39 patients with germline mutations in cancer predisposition genes; 4 of 18 patients had therapy-related myeloid neoplasms (22.2%); and 4 of 15 patients had secondary lymphoid malignancies (26.7%). Notably, of 14 patients who initially suffered from thyroid cancer, 5 patients (35.7%) had germline mutations. Malignancy of lymphoid tissue showed no association with radioactive iodine therapy but was observed to a greater extent in germline mutation-positive thyroid cancer patients regardless of their history of treatment. We observed that 24 of 39 patients (61.5%) were CHIP carriers. Patients who had secondary lymphoid malignancy were less likely to have CHIP than those who had myeloid malignancy. Conclusions: In patients with primary solid cancer who are planning to undergo cytotoxic chemotherapy, radiotherapy, or radioactive iodine therapy, an initial assessment with germline mutation testing using an expanded NGS panel, including low, moderate, and high-risk cancer-associated genes, and somatic CHIP mutation testing can screen the patients who are at risk of developing therapy-related myeloid and lymphoid malignancies. Through careful screening and monitoring throughout the treatment process, patients can benefit from the early detection of secondary malignancies and receive proper treatment.

Quantitative and Qualitative QC of Next-Generation Sequencing for Detecting Somatic Variants: An Example of Detecting Clonal Hematopoiesis of Indeterminate Potential.

BACKGROUND: Because next-generation sequencing (NGS) for detecting somatic mutations has been adopted in clinical fields, both qualitative and quantitative QC of the somatic variants through whole coding regions detected by NGS is crucial. However, specific applications or guidelines, especially for quantitative QC, are currently insufficient. Our goal was to devise a practical approach for both quantitative and qualitative QC using an example of detecting clonal hematopoiesis of indeterminate potential (CHIP). METHODS: We applied the QC scheme using commercial reference materials and in-house QC materials (IQCM) composed of haplotype map and cancer cell lines for monitoring CHIP. RESULTS: This approach efficiently validated a customized CHIP NGS assay. Accuracy, analytical sensitivity, analytical specificity, qualitative precision (concordance), and limit of detection achieved were 99.87%, 98.53%, 100.00%, 100.00%, and 1.00%, respectively. The quantitative precision analysis also had a higher CV percentage at a lower alternative read depth (R2 = 0.749∼0.858). Use of IQCM ensured more than 100-fold reduction in the cost per run compared with that achieved using commercial reference materials. CONCLUSION: Our approach determined the general analytical performance of NGS for detecting CHIP and recognized limitations such as lower precision at a lower level of variant burden. This approach could also be theoretically expanded to a general NGS assay for detecting somatic variants. Considering the reliable NGS results and cost-effectiveness, we propose the use of IQCM for QC of NGS assays at clinical laboratories.

DeviCNV: detection and visualization of exon-level copy number variants in targeted next-generation sequencing data.

BACKGROUND: Targeted next-generation sequencing (NGS) is increasingly being adopted in clinical laboratories for genomic diagnostic tests. RESULTS: We developed a new computational method, DeviCNV, intended for the detection of exon-level copy number variants (CNVs) in targeted NGS data. DeviCNV builds linear regression models with bootstrapping for every probe to capture the relationship between read depth of an individual probe and the median of read depth values of all probes in the sample. From the regression models, it estimates the read depth ratio of the observed and predicted read depth with confidence interval for each probe which is applied to a circular binary segmentation (CBS) algorithm to obtain CNV candidates. Then, it assigns confidence scores to those candidates based on the reliability and strength of the CNV signals inferred from the read depth ratios of the probes within them. Finally, it also provides gene-centric plots with confidence levels of CNV candidates for visual inspection. We applied DeviCNV to targeted NGS data generated for newborn screening and demonstrated its ability to detect novel pathogenic CNVs from clinical samples. CONCLUSIONS: We propose a new pragmatic method for detecting CNVs in targeted NGS data with an intuitive visualization and a systematic method to assign confidence scores for candidate CNVs. Since DeviCNV was developed for use in clinical diagnosis, sensitivity is increased by the detection of exon-level CNVs.

A New Integrated Newborn Screening Workflow Can Provide a Shortcut to Differential Diagnosis and Confirmation of Inherited Metabolic Diseases.

PURPOSE: We developed a new workflow design which included results from both biochemical and targeted gene sequencing analysis interpreted comprehensively. We then conducted a pilot study to evaluate the benefit of this new approach in newborn screening (NBS) and demonstrated the efficiency of this workflow in detecting causative genetic variants. MATERIALS AND METHODS: Ten patients in Group 1 were diagnosed clinically using biochemical assays only, and 10 newborns in Group 2 were diagnosed with suspected inherited metabolic disease (IMD) in NBS. We applied NewbornDiscovery (SD Genomics), an integrated workflow design that encompasses analyte-phenotype-gene, single nucleotide variant/small insertion and deletion/copy number variation analyses along with clinical interpretation of genetic variants related to each participant's condition. RESULTS: A molecular genetic diagnosis was established in 95% (19/20) of individuals. In Group 1, 13 and 7 of 20 alleles were classified as pathogenic and likely pathogenic, respectively. In Group 2, 11 and 6 of 17 alleles with identified causative variants were pathogenic and likely pathogenic, respectively. There were no variants of uncertain significance. For each individual, the NewbornDiscovery and biochemical analysis results reached 100% concordance, since the single newborn testing negative for causative genetic variant in Group 2 showed a benign clinical course. CONCLUSION: This integrated diagnostic workflow resulted in a high yield. This approach not only enabled early confirmation of specific IMD, but also detected conditions not included in the current NBS.

Mutations in DDX58, which encodes RIG-I, cause atypical Singleton-Merten syndrome.

Singleton-Merten syndrome (SMS) is an autosomal-dominant multi-system disorder characterized by dental dysplasia, aortic calcification, skeletal abnormalities, glaucoma, psoriasis, and other conditions. Despite an apparent autosomal-dominant pattern of inheritance, the genetic background of SMS and information about its phenotypic heterogeneity remain unknown. Recently, we found a family affected by glaucoma, aortic calcification, and skeletal abnormalities. Unlike subjects with classic SMS, affected individuals showed normal dentition, suggesting atypical SMS. To identify genetic causes of the disease, we performed exome sequencing in this family and identified a variant (c.1118A>C [p.Glu373Ala]) of DDX58, whose protein product is also known as RIG-I. Further analysis of DDX58 in 100 individuals with congenital glaucoma identified another variant (c.803G>T [p.Cys268Phe]) in a family who harbored neither dental anomalies nor aortic calcification but who suffered from glaucoma and skeletal abnormalities. Cys268 and Glu373 residues of DDX58 belong to ATP-binding motifs I and II, respectively, and these residues are predicted to be located closer to the ADP and RNA molecules than other nonpathogenic missense variants by protein structure analysis. Functional assays revealed that DDX58 alterations confer constitutive activation and thus lead to increased interferon (IFN) activity and IFN-stimulated gene expression. In addition, when we transduced primary human trabecular meshwork cells with c.803G>T (p.Cys268Phe) and c.1118A>C (p.Glu373Ala) mutants, cytopathic effects and a significant decrease in cell number were observed. Taken together, our results demonstrate that DDX58 mutations cause atypical SMS manifesting with variable expression of glaucoma, aortic calcification, and skeletal abnormalities without dental anomalies.

Bacterial genome mapper: A comparative bacterial genome mapping tool.

UNLABELLED: Recently, next generation sequencing (NGS) technologies have led to a revolutionary increase in sequencing speed and costefficacy. Consequently, a vast number of contigs from many recently sequenced bacterial genomes remain to be accurately mapped and annotated, requiring the development of more convenient bioinformatics programs. In this paper, we present a newly developed web-based bioinformatics program, Bacterial Genome Mapper, which is suitable for mapping and annotating contigs that have been assembled from bacterial genome sequence raw data. By constructing a multiple alignment map between target contig sequences and two reference bacterial genome sequences, this program also provides very useful comparative genomics analysis of draft bacterial genomes. AVAILABILITY: The database is available for free at http://mbgm.kribb.re.kr.

CACG: a database for comparative analysis of conjoined genes.

A conjoined gene is defined as one formed at the time of transcription by combining at least part of one exon from each of two or more distinct genes that lie on the same chromosome, in the same or opposite orientation, which translate independently into different proteins. We comparatively studied the extent of conjoined genes in thirteen genomes by analyzing the public databases of expressed sequence tags and mRNA sequences using a set of computational tools designed to identify conjoined genes on the same DNA strand or opposite DNA strands of the same genomic locus. The CACG database, available at http://cgc.kribb.re.kr/map/, includes a number of conjoined genes (7131-human, 2-chimpanzee, 5-orangutan, 57-chicken, 4-rhesus monkey, 651-cow, 27-dog, 2512-mouse, 263-rat, 1482-zebrafish, 5-horse, 29-sheep, and 8-medaka) and is very effective and easy to use to analyze the evolutionary process of conjoined genes when comparing different species.

Genome sequence of Peptoniphilus rhinitidis 1-13T, an anaerobic coccus strain isolated from clinical specimens.

A new Peptoniphilus species has been isolated from samples from a patient who was scheduled for endoscopic sinus surgery for chronic rhinosinusitis. The isolate, Peptoniphilus rhinitidis 1-13(T) (KCTC 5985(T)), can use peptone as a sole carbon source and produce butyrate as a metabolic end product. This is the first report of the draft genome sequence of a novel species in the genus Peptoniphilus within the group of Gram-positive anaerobic cocci.

Genome analysis of the domestic dog (Korean Jindo) by massively parallel sequencing.

Although pioneering sequencing projects have shed light on the boxer and poodle genomes, a number of challenges need to be met before the sequencing and annotation of the dog genome can be considered complete. Here, we present the DNA sequence of the Jindo dog genome, sequenced to 45-fold average coverage using Illumina massively parallel sequencing technology. A comparison of the sequence to the reference boxer genome led to the identification of 4 675 437 single nucleotide polymorphisms (SNPs, including 3 346 058 novel SNPs), 71 642 indels and 8131 structural variations. Of these, 339 non-synonymous SNPs and 3 indels are located within coding sequences (CDS). In particular, 3 non-synonymous SNPs and a 26-bp deletion occur in the TCOF1 locus, implying that the difference observed in cranial facial morphology between Jindo and boxer dogs might be influenced by those variations. Through the annotation of the Jindo olfactory receptor gene family, we found 2 unique olfactory receptor genes and 236 olfactory receptor genes harbouring non-synonymous homozygous SNPs that are likely to affect smelling capability. In addition, we determined the DNA sequence of the Jindo dog mitochondrial genome and identified Jindo dog-specific mtDNA genotypes. This Jindo genome data upgrade our understanding of dog genomic architecture and will be a very valuable resource for investigating not only dog genetics and genomics but also human and dog disease genetics and comparative genomics.

Genome sequence of Myroides injenensis M09-0166(T), isolated from clinical specimens.

A new Myroides species has been isolated from the urine of a patient with fever in spite of multiple antibiotic treatments who had undergone a radical hysterectomy for cervical cancer and percutaneous nephrostomies for hydronephrosis in the past. The isolate, Myroides injenensis M09-0166(T) (KCTC 23367(T)), showed a high level of resistance to multiple antibiotic agents. Here we provide the first report of the draft genome sequence of a novel species in the genus Myroides within the nonfermenting Gram-negative group.

Genome sequence of the anaerobic bacterium Clostridium arbusti SL206(T).

A new Clostridium species has been isolated from pear orchard soil in Daejeon, Republic of Korea. The isolate, Clostridium arbusti SL206(T) (KCTC 5449(T)), showed a nitrogenase activity as well as an organic acid production. Here we first report the draft genome sequence of a novel species in the genus Clostridium within the largest Gram-positive group.

SpiroESTdb: a transcriptome database and online tool for sparganum expressed sequences tags.

BACKGROUND: Sparganum (plerocercoid of Spirometra erinacei) is a parasite that possesses the remarkable ability to survive by successfully modifying its physiology and morphology to suit various hosts and can be found in various tissues, even the nervous system. However, surprisingly little is known about the molecular function of genes that are expressed during the course of the parasite life cycle. To begin to decipher the molecular processes underlying gene function, we constructed a database of expressed sequence tags (ESTs) generated from sparganum. FINDINGS: SpiroESTdb is a web-based information resource that is built upon the annotation and curation of 5,655 ESTs data. SpiroESTdb provides an integrated platform for expressed sequence data, expression dynamics, functional genes, genetic markers including single nucleotide polymorphisms and tandem repeats, gene ontology and KEGG pathway information. Moreover, SpiroESTdb supports easy access to gene pages, such as (i) curation and query forms, (ii) in silico expression profiling and (iii) BLAST search tools. Comprehensive descriptions of the sparganum content of all sequenced data are available, including summary reports. The contents of SpiroESTdb can be viewed and downloaded from the web (http://pathod.cdc.go.kr/spiroestdb). CONCLUSIONS: This integrative web-based database of sequence data, functional annotations and expression profiling data will serve as a useful tool to help understand and expand the characterization of parasitic infections. It can also be used to identify potential industrial drug targets and vaccine candidate genes.

Genome sequence of Lactobacillus fructivorans KCTC 3543.

Lactobacillus fructivorans is important in the generation of particular flavors and in other ripening processes associated with fermented food. Here, we present the draft genome sequence of the type strain Lactobacillus fructivorans KCTC 3543 (1,373,326 bp, with a G+C content of 38.9%), which consists of 5 scaffolds. The genome sequence was obtained by using a whole-genome shotgun strategy with Roche 454 GS (FLX Titanium) pyrosequencing, and all of the reads were assembled using Newbler Assembler 2.3.

Novel mechanism of conjoined gene formation in the human genome.

Recently, conjoined genes (CGs) have emerged as important genetic factors necessary for understanding the human genome. However, their formation mechanism and precise structures have remained mysterious. Based on a detailed structural analysis of 57 human CG transcript variants (CGTVs, discovered in this study) and all (833) known CGs in the human genome, we discovered that the poly(A) signal site from the upstream parent gene region is completely removed via the skipping or truncation of the final exon; consequently, CG transcription is terminated at the poly(A) signal site of the downstream parent gene. This result led us to propose a novel mechanism of CG formation: the complete removal of the poly(A) signal site from the upstream parent gene is a prerequisite for the CG transcriptional machinery to continue transcribing uninterrupted into the intergenic region and downstream parent gene. The removal of the poly(A) signal sequence from the upstream gene region appears to be caused by a deletion or truncation mutation in the human genome rather than post-transcriptional trans-splicing events. With respect to the characteristics of CG sequence structures, we found that intergenic regions are hot spots for novel exon creation during CGTV formation and that exons farther from the intergenic regions are more highly conserved in the CGTVs. Interestingly, many novel exons newly created within the intergenic and intragenic regions originated from transposable element sequences. Additionally, the CGTVs showed tumor tissue-biased expression. In conclusion, our study provides novel insights into the CG formation mechanism and expands the present concepts of the genetic structural landscape, gene regulation, and gene formation mechanisms in the human genome.

Genome sequence of Leuconostoc carnosum KCTC 3525.

We announce the draft genome sequence of the type strain Leuconostoc carnosum KCTC 3525 (3,234,408 bp with a G+C content of 40.9%), one of the most prevalent lactic acid bacteria present during the manufacturing process of vacuum-packaged meats, which consists of 2,407 large contigs (>500 bp in size). The genome sequence was obtained by a whole-genome shotgun strategy using Roche 454 GS (FLX Titanium) pyrosequencing, and all of the reads were assembled using Newbler Assembler 2.3.

Genome sequence of Lactobacillus suebicus KCTC 3549.

Lactobacillus suebicus is important in the generation of particular flavors and in other ripening processes associated with apple mash. Here, we present the draft genome sequence of the type strain Lactobacillus suebicus KCTC 3549 (2,656,936 bp, with a G+C content of 39.0%), which consists of 143 large contigs (>100 bp).

Draft genome sequence of Lactobacillus malefermentans KCTC 3548.

We announce the draft genome sequence of the type strain Lactobacillus malefermentans KCTC 3548 (2,003,922 bp, with a G+C content of 41.1%), which is one of the most prevalent lactic acid bacteria present during the manufacturing process of beer; the genome consists of 172 large contigs (>100 bp in size). All of the contigs were assembled by using Newbler Assembler 2.3 (454 Life Science).

Genome sequence of Lactobacillus versmoldensis KCTC 3814.

Lactobacillus versmoldensis KCTC 3814 was isolated from raw fermented poultry salami. The species was present in high numbers and frequently dominated the lactic acid bacteria (LAB) populations of the products. Here, we announce the draft genome sequence of Lactobacillus versmoldensis KCTC 3814, isolated from poultry salami, and describe major findings from its annotation.

Draft genome sequence of Lactobacillus zeae KCTC 3804.

We announce the draft genome sequence of the type strain Lactobacillus zeae KCTC 3804 (3,110,326 bp, with a G+C content of 47.8%), which is one of the most prevalent lactic acid bacteria present during the processing of raw cow's milk. The genome consists of 113 large contigs (>100 bp). All of the contigs were assembled by Newbler Assembler 2.3 (454 Life Science).