Sequencing artifacts derived from a library preparation method using enzymatic fragmentation.

DNA fragmentation is a fundamental step during library preparation in hybridization capture-based, short-read sequencing. Ultra-sonication has been used thus far to prepare DNA of an appropriate size, but this method is associated with a considerable loss of DNA sample. More recently, studies have employed library preparation methods that rely on enzymatic fragmentation with DNA endonucleases to minimize DNA loss, particularly in nano-quantity samples. Yet, despite their wide use, the effect of enzymatic fragmentation on the resultant sequences has not been carefully assessed. Here, we used pairwise comparisons of somatic variants of the same tumor DNA samples prepared using ultrasonic and enzymatic fragmentation methods. Our analysis revealed a substantially larger number of recurrent artifactual SNVs/indels in endonuclease-treated libraries as compared with those created through ultrasonication. These artifacts were marked by palindromic structure in the genomic context, positional bias in sequenced reads, and multi-nucleotide substitutions. Taking advantage of these distinctive features, we developed a filtering algorithm to distinguish genuine somatic mutations from artifactual noise with high specificity and sensitivity. Noise cancelling recovered the composition of the mutational signatures in the tumor samples. Thus, we provide an informatics algorithm as a solution to the sequencing errors produced as a consequence of endonuclease-mediated fragmentation, highlighted for the first time in this study.

Genome sequence of a mesophilic hydrogenotrophic methanogen Methanocella paludicola, the first cultivated representative of the order Methanocellales.

We report complete genome sequence of a mesophilic hydrogenotrophic methanogen Methanocella paludicola, the first cultured representative of the order Methanocellales once recognized as an uncultured key archaeal group for methane emission in rice fields. The genome sequence of M. paludicola consists of a single circular chromosome of 2,957,635 bp containing 3004 protein-coding sequences (CDS). Genes for most of the functions known in the methanogenic archaea were identified, e.g. a full complement of hydrogenases and methanogenesis enzymes. The mixotrophic growth of M. paludicola was clarified by the genomic characterization and re-examined by the subsequent growth experiments. Comparative genome analysis with the previously reported genome sequence of RC-I(MRE50), which was metagenomically reconstructed, demonstrated that about 70% of M. paludicola CDSs were genetically related with RC-I(MRE50) CDSs. These CDSs included the genes involved in hydrogenotrophic methane production, incomplete TCA cycle, assimilatory sulfate reduction and so on. However, the genetic components for the carbon and nitrogen fixation and antioxidant system were different between the two Methanocellales genomes. The difference is likely associated with the physiological variability between M. paludicola and RC-I(MRE50), further suggesting the genomic and physiological diversity of the Methanocellales methanogens. Comparative genome analysis among the previously determined methanogen genomes points to the genome-wide relatedness of the Methanocellales methanogens to the orders Methanosarcinales and Methanomicrobiales methanogens in terms of the genetic repertoire. Meanwhile, the unique evolutionary history of the Methanocellales methanogens is also traced in an aspect by the comparative genome analysis among the methanogens.

Whole-genome analysis of Salmonella enterica serovar Typhimurium T000240 reveals the acquisition of a genomic island involved in multidrug resistance via IS1 derivatives on the chromosome.

Salmonella enterica serovar Typhimurium is frequently associated with life-threatening systemic infections, and the recent global emergence of multidrug resistance in S. enterica isolates from agricultural and clinical settings has raised concerns. In this study, we determined the whole-genome sequence of fluoroquinolone-resistant S. enterica serovar Typhimurium T000240 strain (DT12) isolated from human gastroenteritis in 2000. Comparative genome analysis revealed that T000240 displays high sequence similarity to strain LT2, which was originally isolated in 1940, indicating that progeny of LT2 might be reemerging. T000240 possesses a unique 82-kb genomic island, designated as GI-DT12, which is composed of multidrug resistance determinants, including a Tn2670-like composite transposon (class 1 integron [intI1, bla(oxa-30), aadA1, qacEΔ1, and sul1], mercury resistance proteins, and chloramphenicol acetyltransferase), a Tn10-like tetracycline resistance protein (tetA), the aerobactin iron-acquisition siderophore system (lutA and lucABC), and an iron transporter (sitABCD). Since GI-DT12 is flanked by IS1 derivatives, IS1-mediated recombination likely played a role in the acquisition of this genomic island through horizontal gene transfer. The aminoglycoside-(3)-N-acetyltransferase (aac(3)) gene and a class 1 integron harboring the dfrA1 gene cassette responsible for gentamicin and trimethoprim resistance, respectively, were identified on plasmid pSTMDT12_L and appeared to have been acquired through homologous recombination with IS26. This study represents the first characterization of the unique genomic island GI-DT12 that appears to be associated with possible IS1-mediated recombination in S. enterica serovar Typhimurium. It is expected that future whole-genome studies will aid in the characterization of the horizontal gene transfer events for the emerging S. enterica serovar Typhimurium strains.

Sequence analysis of three plasmids harboured in Rhodococcus erythropolis strain PR4.

Rhodococcus erythropolis strain PR4 has been isolated as an alkane-degrading bacterium. The strain harbours one linear plasmid, pREL1 (271 577 bp) and two circular plasmids, pREC1 (104 014 bp) and pREC2 (3637 bp), all with some sequence similarities to other Rhodococcus plasmids. For pREL1, pREC1 and pREC2, 298, 102 and 3 open reading frames, respectively, were predicted. Linear plasmid pREL1 has several regions homologous to plasmid pBD2 found in R. erythropolis BD2. Sequence analysis of pREL1 and pBD2 identified common metal-resistance genes on both, but pREL1 also encodes alkane-degradation genes not found on pBD2, with enzyme constituents some of which are quite different from those of other organisms. The alkane hydroxylase consisted of a cytochrome P450 monooxygenase, a 2Fe-2S ferredoxin, and a ferredoxin reductase. The ferredoxin reductase amino acid sequence resembles the AlkT (rubredoxin reductase) sequence. A zinc-containing alcohol dehydrogenase further oxydizes alkanols, alkane oxidation products catalysed by alkane hydroxylase. Of the circular plasmids, the pREC1 sequence is partially similar to the sequence of pREAT701, the virulence plasmid found in Rhodococcus equi. pREC1 has no pREAT701 virulence genes and encodes genes for beta-oxidation of fatty acids. Thus, joint actions of enzymes encoded by pREL1 and pREC1 may enable efficient mineralization of alkanes.