De novo bacterial genome sequencing: millions of very short reads assembled on a desktop computer

Current status in whole genome sequencing and analysis of Ipomoea spp

The recent advances of next-generation sequencing have made it possible to construct reference genome sequences in divergent species. However, de novo assembly at the chromosome level remains challenging in polyploid species, due to the existence of more than two pairs of homoeologous chromosomes in one nucleus. Cultivated sweet potato (Ipomoea batatas (L.) Lam) is a hexaploid species with 90 chromosomes (2n = 6X = 90). Although the origin of sweet potato is also still under discussion, diploid relative species, I. trifida and I. triloba have been considered as one of the most possible progenitors. In this manuscript, we review the recent results and activities of whole-genome sequencing in the genus Ipomoea series Batatas, I. trifida, I. triloba and sweet potato (I. batatas). Most of the results of genome assembly suggest that the genomes of sweet potato consist of two pairs and four pairs of subgenomes, i.e., B1B1B2B2B2B2. The results also revealed the relation between sweet potato and other Ipomoea species. Together with the development of bioinformatics approaches, the large-scale publicly available genome and transcript sequence resources and international genome sequencing streams are expected to promote the genome sequence dissection in sweet potato.

Development and characterization of genomic microsatellite markers in the tree species, Rhodoleia championii, R. parvipetala, and R. forrestii (Hamamelidaceae)

Rhodoleia Champion ex Hooker is one of the most primitive relict genera of Hamamelidaceae, a key family exploited to understand the origin and early evolution of flowering plants. Genomic simple sequence repeats (SSRs) were developed for R. championii to perform genetic diversity, phylogeographical structure or even systematic evolution studies of the genus. Among the 278,743 contigs (105,758,242 bps) de novo assembled from the low-coverage whole genome sequencing of R. championii, a total of 9106 SSRs were detected in 8370 contigs, and SSR primer pairs were successfully designed for 6677 SSRs. Among the 110 selected primer pairs, 41 were amplified successfully in the preliminary test of SSR screening. Further amplification of these 41 primer pairs across the 122 individuals collected from six populations of the three Rhodoleia species showed that 32 and 40 SSR markers can be amplified in Vietnam and Jinping populations of R. parvipetala, 41, 33, and 41 SSR markers in Boluo, Hongkong and Xinyi populations of R. championii, 25 SSR markers in Fugong population of R. forrestii, and 20 SSR markers demonstrated to be polymorphic across the three species. Genetic analysis for these 20 polymorphic SSRs showed that Allele number (A) ranged from four to 13 and polymorphic information content (PIC) ranged from 0.479 to 0.876 across the three species. At the population level, observed heterozygosity (H_O) ranged from 0.000 to 1.000, and expected heterozygosity (H_E) ranged from 0.091 to 0.851. In the present study, we provided the first whole-genome sequencing database for the species R. championii, identified ample SSR loci with designed primers, and revealed that 20 of the 110 selected SSRs were polymorphic across three Rhodoleia species. These provide valuable resources for future studies on genetic study, species delimitation, phylogeography, and conservation of this genus.

Comprehensive genomic analysis of an indigenous Pseudomonas pseudoalcaligenes degrading phenolic compounds

Environmental contamination with aromatic compounds is a universal challenge. Aromatic-degrading microorganisms isolated from the same or similar polluted environments seem to be more suitable for bioremediation. Moreover, microorganisms adapted to contaminated environments are able to use toxic compounds as the sole sources of carbon and energy. An indigenous strain of Pseudomonas, isolated from the Mahshahr Petrochemical plant in the Khuzestan province, southwest of Iran, was studied genetically. It was characterized as a novel Gram-negative, aerobic, halotolerant, rod-shaped bacterium designated Pseudomonas YKJ, which was resistant to chloramphenicol and ampicillin. Genome of the strain was completely sequenced using Illumina technology to identify its genetic characteristics. MLST analysis revealed that the YKJ strain belongs to the genus Pseudomonas indicating the highest sequence similarity with Pseudomonas pseudoalcaligenes strain CECT 5344 (99% identity). Core- and pan-genome analysis indicated that P. pseudoalcaligenes contains 1,671 core and 3,935 unique genes for coding DNA sequences. The metabolic and degradation pathways for aromatic pollutants were investigated using the NCBI and KEGG databases. Genomic and experimental analyses showed that the YKJ strain is able to degrade certain aromatic compounds including bisphenol A, phenol, benzoate, styrene, xylene, benzene and chlorobenzene. Moreover, antibiotic resistance and chemotaxis properties of the YKJ strain were found to be controlled by two-component regulatory systems.

Novo&Stitch: accurate reconciliation of genome assemblies via optical maps

Motivation

De novo genome assembly is a challenging computational problem due to the high repetitive content of eukaryotic genomes and the imperfections of sequencing technologies (i.e. sequencing errors, uneven sequencing coverage and chimeric reads). Several assembly tools are currently available, each of which has strengths and weaknesses in dealing with the trade-off between maximizing contiguity and minimizing assembly errors (e.g. mis-joins). To obtain the best possible assembly, it is common practice to generate multiple assemblies from several assemblers and/or parameter settings and try to identify the highest quality assembly. Unfortunately, often there is no assembly that both maximizes contiguity and minimizes assembly errors, so one has to compromise one for the other.

Results

The concept of assembly reconciliation has been proposed as a way to obtain a higher quality assembly by merging or reconciling all the available assemblies. While several reconciliation methods have been introduced in the literature, we have shown in one of our recent papers that none of them can consistently produce assemblies that are better than the assemblies provided in input. Here we introduce Novo&Stitch, a novel method that takes advantage of optical maps to accurately carry out assembly reconciliation (assuming that the assembled contigs are sufficiently long to be reliably aligned to the optical maps, e.g. 50 Kbp or longer). Experimental results demonstrate that Novo&Stitch can double the contiguity (N50) of the input assemblies without introducing mis-joins or reducing genome completeness.

Availability and implementation

Novo&Stitch can be obtained from https://github.com/ucrbioinfo/Novo_Stitch.

Bioinformatics of nanopore sequencing

Nanopore sequencing is one of the most exciting new technologies that undergo dynamic development. With its development, a growing number of analytical tools are becoming available for researchers. To help them better navigate this ever changing field, we discuss a range of software available to analyze sequences obtained using nanopore technology.

GAAP: A Genome Assembly + Annotation Pipeline

Genomic analysis begins with de novo assembly of short-read fragments in order to reconstruct full-length base sequences without exploiting a reference genome sequence. Then, in the annotation step, gene locations are identified within the base sequences, and the structures and functions of these genes are determined. Recently, a wide range of powerful tools have been developed and published for whole-genome analysis, enabling even individual researchers in small laboratories to perform whole-genome analyses on their objects of interest. However, these analytical tools are generally complex and use diverse algorithms, parameter setting methods, and input formats; thus, it remains difficult for individual researchers to select, utilize, and combine these tools to obtain their final results. To resolve these issues, we have developed a genome analysis pipeline (GAAP) for semiautomated, iterative, and high-throughput analysis of whole-genome data. This pipeline is designed to perform read correction, de novo genome (transcriptome) assembly, gene prediction, and functional annotation using a range of proven tools and databases. We aim to assist non-IT researchers by describing each stage of analysis in detail and discussing current approaches. We also provide practical advice on how to access and use the bioinformatics tools and databases and how to implement the provided suggestions. Whole-genome analysis of Toxocara canis is used as case study to show intermediate results at each stage, demonstrating the practicality of the proposed method.

Complete Sequence and Annotation of the Mycoplasma phocicerebrale Strain 1049T Genome

The Mycoplasma phocicerebrale genome was analyzed to better understand this opportunistic pathogen. Amplification with ϕ29 polymerase was used to generate enough genomic DNA for large-insert library construction. Like other mycoplasmas from seals, M. phocicerebrale encodes an immunosuppressor that may predispose susceptibility to infection or influence intercurrent diseases of affected hosts.

The Mycoplasma phocicerebrale genome was analyzed to better understand this opportunistic pathogen. Amplification with ϕ29 polymerase was used to generate enough genomic DNA for large-insert library construction. Like other mycoplasmas from seals, M. phocicerebrale encodes an immunosuppressor that may predispose susceptibility to infection or influence intercurrent diseases of affected hosts.

Identifying similar transcripts in a related organism from de Bruijn graphs of RNA-Seq data, with applications to the study of salt and waterlogging tolerance in Melilotus

Background

A popular strategy to study alternative splicing in non-model organisms starts from sequencing the entire transcriptome, then assembling the reads by using de novo transcriptome assembly algorithms to obtain predicted transcripts. A similarity search algorithm is then applied to a related organism to infer possible function of these predicted transcripts. While some of these predictions may be inaccurate and transcripts with low coverage are often missed, we observe that it is possible to obtain a more complete set of transcripts to facilitate possible functional assignments by starting the search from the intermediate de Bruijn graph that contains all branching possibilities.

Results

We develop an algorithm to extract similar transcripts in a related organism by starting the search from the de Bruijn graph that represents the transcriptome instead of from predicted transcripts. We show that our algorithm is able to recover more similar transcripts than existing algorithms, with large improvements in obtaining longer transcripts and a finer resolution of isoforms. We apply our algorithm to study salt and waterlogging tolerance in two Melilotus species by constructing new RNA-Seq libraries.

Conclusions

We have developed an algorithm to identify paths in the de Bruijn graph that correspond to similar transcripts in a related organism directly. Our strategy bypasses the transcript prediction step in RNA-Seq data and makes use of support from evolutionary information.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5702-5) contains supplementary material, which is available to authorized users.

Identification and development of microsatellite (SSRs) makers of Exbucklandia (HAMAMELIDACEAE) by high-throughput sequencing

Hamamelidaceae (Saxifragales, previously Rosales) comprises approximately six subfamily, 30 genera and 140 species, most of which are Tertiary relicts. Exbucklandia is the only genus of the subfamily Exbucklandioideae, Hamelidaceae, containing only 2-4 species. Of them, the species E. longipetala H. T. Chang is endemic to China and listed as endangered in The Biodiversity Red List of China: Higher Plant, yet some taxonomists put forward that E. longipetala should be merged into E. tonkinensis (Lecomte) H. T. Chang. Currently, there was nearly no phylogeographic studies on this genus possibly due to the deficiency of efficient molecular markers. In this study, we sequenced the genome of E. tonkinensis based on high throughput sequencing technology, and obtained approximately 6 G raw data, which was further de novo assembled into 303,481 contigs. Based on them, 15,326 SSRs were identified from 13,596 contigs, and primers were successfully designed for 10,660 SSRs. A total of 139 paired primers were synthesized, 106 of them were successfully amplified in six Exbucklandia individuals with expected PCR product size, and 24 demonstrated to be polymorphic among three Exbucklandia populations. Accordingly, the expected and observed heterozygosity were between 0.097-0.717 and 0.098-0.583. Based on these efforts, future researches on genetic diversity and population structure of Exbucklandia can be performed to understand its phylogenetic origins and phylogeographic pattern.

Sequencing the Obligate Intracellular Rhabdochlamydia helvetica within Its Tick Host Ixodes ricinus to Investigate Their Symbiotic Relationship

The Rhabdochlamydiaceae family is one of the most widely distributed within the phylum Chlamydiae, but most of its members remain uncultivable. Rhabdochlamydia 16S rRNA was recently reported in more than 2% of 8,534 pools of ticks from Switzerland. Shotgun metagenomics was performed on a pool of five female Ixodes ricinus ticks presenting a high concentration of chlamydial DNA, allowing the assembly of a high-quality draft genome. About 60% of sequence reads originated from a single bacterial population that was named "Candidatus Rhabdochlamydia helvetica" whereas only few thousand reads mapped to the genome of "Candidatus Midichloria mitochondrii," a symbiont normally observed in all I. ricinus females. The 1.8 Mbp genome of R. helvetica is smaller than other Chlamydia-related bacteria. Comparative analyses with other chlamydial genomes identified transposases of the PD-(D/E)XK nuclease family that are unique to this new genome. These transposases show evidence of interphylum horizontal gene transfers between multiple arthropod endosymbionts, including Cardinium spp. (Bacteroidetes) and diverse proteobacteria such as Wolbachia, Rickettsia spp. (Rickettsiales), and Caedimonas varicaedens (Holosporales). Bacterial symbionts were previously suggested to provide B-vitamins to hematophagous hosts. However, incomplete metabolic capacities including for B-vitamin biosynthesis, high bacterial density and limited prevalence suggest that R. helvetica is parasitic rather than symbiotic to its host. The identification of novel Rhabdochlamydia strains in different hosts and their sequencing will help understanding if members of this genus have become highly specialized parasites with reduced genomes, like the Chlamydiaceae, or if they could be pathogenic to humans using ticks as a transmission vector.

Complete Sequence and Annotation of the Mycoplasma phocirhinis Strain 852T Genome

The genome of Mycoplasma phocirhinis strain 852^T was examined for determinants of tropism or virulence. It encodes multiple orthologs of an immunosuppressor that may predispose susceptibility to infection or influence outcomes of intercurrent diseases in marine mammals.

The genome of Mycoplasma phocirhinis strain 852^T was examined for determinants of tropism or virulence. It encodes multiple orthologs of an immunosuppressor that may predispose susceptibility to infection or influence outcomes of intercurrent diseases in marine mammals.

Genome and pan-genome analysis to classify emerging bacteria

Background

In the recent years, genomic and pan-genomic studies have become increasingly important. Culturomics allows to study human microbiota through the use of different culture conditions, coupled with a method of rapid identification by MALDI-TOF, or 16S rRNA. Bacterial taxonomy is undergoing many changes as a consequence. With the help of pan-genomic analyses, species can be redefined, and new species definitions generated.

Results

Genomics, coupled with culturomics, has led to the discovery of many novel bacterial species or genera, including Akkermansia muciniphila and Microvirga massiliensis. Using the genome to define species has been applied within the genus Klebsiella. A discontinuity or an abrupt break in the core/pan-genome ratio can uncover novel species.

Conclusions

Applying genomic and pan-genomic analyses to the reclassification of other bacterial species or genera will be important in the future of medical microbiology. The pan-genome is one of many new innovative tools in bacterial taxonomy.

Reviewers

This article was reviewed by William Martin, Eric Bapteste and James Mcinerney.

Open peer review

Reviewed by William Martin, Eric Bapteste and James Mcinerney.

Evolution of sex determination and heterogamety changes in section Otites of the genus Silene

Switches in heterogamety are known to occur in both animals and plants. Although plant sex determination systems probably often evolved more recently than those in several well-studied animals, including mammals, and have had less time for switches to occur, we previously detected a switch in heterogamety in the plant genus Silene: section Otites has both female and male heterogamety, whereas S. latifolia and its close relatives, in a different section of the genus, Melandrium (subgenus Behenantha), all have male heterogamety. Here we analyse the evolution of sex chromosomes in section Otites, which is estimated to have evolved only about 0.55 MYA. Our study confirms female heterogamety in S. otites and newly reveals female heterogamety in S. borysthenica. Sequence analyses and genetic mapping show that the sex-linked regions of these two species are the same, but the region in S. colpophylla, a close relative with male heterogamety, is different. The sex chromosome pairs of S. colpophylla and S. otites each correspond to an autosome of the other species, and both differ from the XY pair in S. latifolia. Silene section Otites species are suitable for detailed studies of the events involved in such changes, and our phylogenetic analysis suggests a possible change from female to male heterogamety within this section. Our analyses suggest a possibility that has so far not been considered, change in heterogamety through hybridization, in which a male-determining chromosome from one species is introgressed into another one, and over-rides its previous sex-determining system.

Shared genome analyses of notable listeriosis outbreaks, highlighting the critical importance of epidemiological evidence, input datasets and interpretation criteria

The persuasiveness of genomic evidence has pressured scientific agencies to supplement or replace well-established methodologies to inform public health and food safety decision-making. This study of 52 epidemiologically defined Listeria monocytogenes isolates, collected between 1981 and 2011, including nine outbreaks, was undertaken (1) to characterize their phylogenetic relationship at finished genome-level resolution, (2) to elucidate the underlying genetic diversity within an endemic subtype, CC8, and (3) to re-evaluate the genetic relationship and epidemiology of a CC8-delimited outbreak in Canada in 2008. Genomes representing Canadian Listeria outbreaks between 1981 and 2010 were closed and manually annotated. Single nucleotide variants (SNVs) and horizontally acquired traits were used to generate phylogenomic models. Phylogenomic relationships were congruent with classical subtyping and epidemiology, except for CC8 outbreaks, wherein the distribution of SNV and prophages revealed multiple co-evolving lineages. Chronophyletic reconstruction of CC8 evolution indicates that prophage-related genetic changes among CC8 strains manifest as PFGE subtype reversions, obscuring the relationship between CC8 isolates, and complicating the public health interpretation of subtyping data, even at maximum genome resolution. The size of the shared genome interrogated did not change the genetic relationship measured between highly related isolates near the tips of the phylogenetic tree, illustrating the robustness of these approaches for routine public health applications where the focus is recent ancestry. The possibility exists for temporally and epidemiologically distinct events to appear related even at maximum genome resolution, highlighting the continued importance of epidemiological evidence.

Prophage induction, but not production of phage particles, is required for lethal disease in a microbiome-replete murine model of enterohemorrhagic E. coli infection

Enterohemorrhagic Escherichia coli (EHEC) colonize intestinal epithelium by generating characteristic attaching and effacing (AE) lesions. They are lysogenized by prophage that encode Shiga toxin 2 (Stx2), which is responsible for severe clinical manifestations. As a lysogen, prophage genes leading to lytic growth and stx2 expression are repressed, whereas induction of the bacterial SOS response in response to DNA damage leads to lytic phage growth and Stx2 production both in vitro and in germ-free or streptomycin-treated mice. Some commensal bacteria diminish prophage induction and concomitant Stx2 production in vitro, whereas it has been proposed that phage-susceptible commensals may amplify Stx2 production by facilitating successive cycles of infection in vivo. We tested the role of phage induction in both Stx production and lethal disease in microbiome-replete mice, using our mouse model encompassing the murine pathogen Citrobacter rodentium lysogenized with the Stx2-encoding phage Φstx2dact. This strain generates EHEC-like AE lesions on the murine intestine and causes lethal Stx-mediated disease. We found that lethal mouse infection did not require that Φstx2dact infect or lysogenize commensal bacteria. In addition, we detected circularized phage genomes, potentially in the early stage of replication, in feces of infected mice, confirming that prophage induction occurs during infection of microbiota-replete mice. Further, C. rodentium (Φstx2dact) mutants that do not respond to DNA damage or express stx produced neither high levels of Stx2 in vitro or lethal infection in vivo, confirming that SOS induction and concomitant expression of phage-encoded stx genes are required for disease. In contrast, C. rodentium (Φstx2dact) mutants incapable of prophage genome excision or of packaging phage genomes retained the ability to produce Stx in vitro, as well as to cause lethal disease in mice. Thus, in a microbiome-replete EHEC infection model, lytic induction of Stx-encoding prophage is essential for lethal disease, but actual phage production is not.

Identification of the common biosynthetic gene cluster for both antimicrobial streptoaminals and antifungal 5-alkyl-1,2,3,4-tetrahydroquinolines

The new subfamily of type II PKS gene cluster is responsible for biosynthesis of structurally distinct streptoaminals (STAMs) and 5-alkyl-1,2,3,4-tetrahydroquinolines (5aTHQs).

Complete chloroplast genome sequence of Betula platyphylla: gene organization, RNA editing, and comparative and phylogenetic analyses

BACKGROUND

Betula platyphylla is a common tree species in northern China that has high economic and medicinal value. Our laboratory has been devoted to genome research on B. platyphylla for approximately 10 years. As primary organelle genomes, the complete genome sequences of chloroplasts are important to study the divergence of species, RNA editing and phylogeny. In this study, we sequenced and analyzed the complete chloroplast (cp) genome sequence of B. platyphylla.

RESULTS

The complete cp genome of B. platyphylla was 160,518 bp in length, which included a pair of inverted repeats (IRs) of 26,056 bp that separated a large single copy (LSC) region of 89,397 bp and a small single copy (SSC) region of 19,009 bp. The annotation contained a total of 129 genes, including 84 protein-coding genes, 37 tRNA genes and 8 rRNA genes. There were 3 genes using alternative initiation codons. Comparative genomics showed that the sequence of the Fagales species cp genome was relatively conserved, but there were still some high variation regions that could be used as molecular markers. The IR expansion event of B. platyphylla resulted in larger cp genomes and rps19 pseudogene formation. The simple sequence repeat (SSR) analysis showed that there were 105 SSRs in the cp genome of B. platyphylla. RNA editing sites recognition indicated that at least 80 RNA editing events occurred in the cp genome. Most of the substitutions were C to U, while a small proportion of them were not. In particular, three editing loci on the rRNA were converted to more than two other bases that had never been reported. For synonymous conversion, most of them increased the relative synonymous codon usage (RSCU) value of the codons. The phylogenetic analysis suggested that B. platyphylla had a closer evolutionary relationship with B. pendula than B. nana.

CONCLUSIONS

In this study, we not only obtained and annotated the complete cp genome sequence of B. platyphylla, but we also identified new RNA editing sites and predicted the phylogenetic relationships among Fagales species. These findings will facilitate genomic, genetic engineering and phylogenetic studies of this important species.

Virulotyping of Salmonella enterica serovar Typhi isolates from Pakistan: Absence of complete SPI-10 in Vi negative isolates

The pathogenesis of Salmonella enterica serovar Typhi (S. Typhi), the cause of typhoid fever in humans, is mainly attributed to the acquisition of horizontally acquired DNA elements. Salmonella pathogenicity islands (SPIs) are indubitably the most important form of horizontally acquired DNA with respect to pathogenesis of this bacterium. The insertion or deletion of any of these transferrable SPIs may have impact on the virulence potential of S. Typhi. In this study, the virulence potential and genetic relatedness of 35 S. Typhi isolates, collected from 2004 to 2013 was determined by identification of SPI and non-SPI virulence factors through a combination of techniques including virulotyping, Whole Genome Sequencing (WGS), and Variable Number of Tandem Repeats (VNTR) profiling. In order to determine the virulence potential of local S. Typhi isolates, 56 virulence related genes were studied by PCR. These genes are located in the core as well as accessory genome (SPIs and plasmid). Major variations among studied virulence determinants were found in case of SPI-7 and SPI-10 associated genes. On the basis of presence of virulence related genes, the studied S. Typhi isolates from Pakistan were clustered into two virulotypes Vi-positive and Vi-negative. Interestingly, SPI-7 and SPI-10 were collectively absent or present in Vi-negative and Vi-positive strains, respectively. Two Vi-negative and 11 Vi-positive S. Typhi strains were also analyzed by whole genome sequencing (WGS) and their results supported the PCR results. Genetic diversity was tested by VNTR-based molecular typing. All 35 isolates were clustered into five groups. Overall, all Vi-negative isolates were placed in a single group (T5) whereas Vi-positive isolates were grouped into four types. Vi-negative and Vi-positive isolates were mutually exclusive. This is the first report on the comparative distribution of SPI and non-SPI related virulence genes in Vi-negative and Vi-positive S. Typhi isolates with an important finding that SPI-10 is absent in all Vi-negative isolates.

The distribution of virulence factors in S. Typhi can vary in isolates from different geographical regions and can have significant effect on the disease control. In this study, we have checked the distribution of 56 reported virulence associated factors in 35 local isolates of S. Typhi to identify any variations that can help in designing effective control strategies for typhoid. We have identified four naturally occurring variants which are simultaneously lacking SPI-7 and SPI-10, two adjacently located pathogenicity islands on S. Typhi chromosome. These isolates are not producing Vi capsular antigen hence the Vi based vaccines will not be effective against them. These findings highlight the need to develop typhoid vaccines specifically effective in Pakistan.

Isolation and Genome Sequence Characterization of Bacteriophage vB_SalM_PM10, a Cba120virus, Concurrently Infecting Salmonella enterica Serovars Typhimurium, Typhi, and Enteritidis

The prevalence of multidrug-resistant Salmonella is ever increasing and calls for alternatives to antibiotics. The use of phages has been anticipated to reduce the multidrug-resistant human pathogens in food environment. Salmonella phage vB_SalM_PM10 (PM10) was isolated from sewage-polluted river in India. It shows an icosahedral head (94 ± 4 nm) along with a long contractile tail (106 ± 7 × 18 ± 2 nm), a morphotype of family Ackermannviridae. Additionally, the phage displayed the features resembling to existing Cba120viruses. Phage PM10 could infect S. enterica serovars Typhimurium, Typhi, and Enteritidis. The genome sequencing analysis of phage PM10 revealed circular 158.08 kb double-stranded DNA, with the GC content of 44.6%. Two hundred and nine ORFs, 171 putative promoters, 122 rho-independent terminators, and 5 transfer RNA encoding genes were found in the genome. The genome-wide comparisons and phylogenetic analyses showed that phage PM10 is closely related to Salmonella phage PhiSH19. Comparison of the tail-spike protein sequences encoded in PM10 and PhiSH19 genome showed the variation, which might have facilitated PM10's simultaneous infectivity to aforementioned S. enterica serovars. This is a varied host range than that of PhiSH19 or any other Cba120viruses.

Two Groups of Cocirculating, Epidemic Clostridiodes difficile Strains Microdiversify through Different Mechanisms

Clostridiodes difficile strains from the NAP_CR1/ST54 and NAP1/ST01 types have caused outbreaks despite of their notable differences in genome diversity. By comparing whole genome sequences of 32 NAP_CR1/ST54 isolates and 17 NAP1/ST01 recovered from patients infected with C. difficile we assessed whether mutation, homologous recombination (r) or nonhomologous recombination (NHR) through lateral gene transfer (LGT) have differentially shaped the microdiversification of these strains. The average number of single nucleotide polymorphisms (SNPs) in coding sequences (NAP_CR1/ST54 = 24; NAP1/ST01 = 19) and SNP densities (NAP_CR1/ST54 = 0.54/kb; NAP1/ST01 = 0.46/kb) in the NAP_CR1/ST54 and NAP1/ST01 isolates was comparable. However, the NAP1/ST01 isolates showed 3× higher average dN/dS rates (8.35) that the NAP_CR1/ST54 isolates (2.62). Regarding r, whereas 31 of the NAP_CR1/ST54 isolates showed 1 recombination block (3,301–8,226 bp), the NAP1/ST01 isolates showed no bases in recombination. As to NHR, the pangenome of the NAP_CR1/ST54 isolates was larger (4,802 gene clusters, 26% noncore genes) and more heterogeneous (644 ± 33 gene content changes) than that of the NAP1/ST01 isolates (3,829 gene clusters, ca. 6% noncore genes, 129 ± 37 gene content changes). Nearly 55% of the gene content changes seen among the NAP_CR1/ST54 isolates (355 ± 31) were traced back to MGEs with putative genes for antimicrobial resistance and virulence factors that were only detected in single isolates or isolate clusters. Congruently, the LGT/SNP rate calculated for the NAP_CR1/ST54 isolates (26.8 ± 2.8) was 4× higher than the one obtained for the NAP1/ST1 isolates (6.8 ± 2.0). We conclude that NHR-LGT has had a greater role in the microdiversification of the NAP_CR1/ST54 strains, opposite to the NAP1/ST01 strains, where mutation is known to play a more prominent role.

Identification and development of microsatellite markers in Hamamelis mollis (Hamamelidaceae)

PREMISE OF THE STUDY

Hamamelis mollis (Hamamelidaceae) is a Tertiary relict species endemic to southern China. Polymorphic microsatellite markers were developed to reveal the genetic diversity of this species.

METHODS AND RESULTS

The genome of H. mollis was sequenced and de novo assembled into 642,351 contigs. A total of 72,097 paired primers were successfully designed from 80,282 simple sequence repeat (SSR) markers identified in 63,419 contigs. PCR amplification showed that 96 of the 136 synthesized primers could be successfully amplified, and 22 demonstrated polymorphism. The mean number of alleles, levels of observed heterozygosity, and levels of expected heterozygosity were 4.602 ± 0.140, 0.632 ± 0.020, and 0.696 ± 0.010, respectively. The majority of the 96 primer pairs could be amplified in at least one other Hamamelidaceae species, including Distylium myricoides (60), Loropetalum chinense (39), Exbucklandia populnea (24), and E. tonkinensis (24).

CONCLUSIONS

These microsatellite loci provide abundant genomic SSR markers to evaluate genetic diversity of this woody ornamental plant.

The linear mitochondrial genome of the quarantine chytrid Synchytrium endobioticum; insights into the evolution and recent history of an obligate biotrophic plant pathogen

BACKGROUND

Chytridiomycota species (chytrids) belong to a basal lineage in the fungal kingdom. Inhabiting terrestrial and aquatic environments, most are free-living saprophytes but several species cause important diseases: e.g. Batrachochytrium dendrobatidis, responsible for worldwide amphibian decline; and Synchytrium endobioticum, causing potato wart disease. S. endobioticum has an obligate biotrophic lifestyle and isolates can be further characterized as pathotypes based on their virulence on a differential set of potato cultivars. Quarantine measures have been implemented globally to control the disease and prevent its spread. We used a comparative approach using chytrid mitogenomes to determine taxonomical relationships and to gain insights into the evolution and recent history of introductions of this plant pathogen.

RESULTS

We assembled and annotated the complete mitochondrial genome of 30 S. endobioticum isolates and generated mitochondrial genomes for five additional chytrid species. The mitochondrial genome of S. endobioticum is linear with terminal inverted repeats which was validated by tailing and PCR amplifying the telomeric ends. Surprisingly, no conservation in organisation and orientation of mitochondrial genes was observed among the Chytridiomycota except for S. endobioticum and its sister species Synchytrium microbalum. However, the mitochondrial genome of S. microbalum is circular and comprises only a third of the 72.9 Kbp found for S. endobioticum suggesting recent linearization and expansion. Four mitochondrial lineages were identified in the S. endobioticum mitochondrial genomes. Several pathotypes occur in different lineages, suggesting that these have emerged independently. In addition, variations for polymorphic sites in the mitochondrial genome of individual isolates were observed demonstrating that S. endobioticum isolates represent a community of different genotypes. Such communities were shown to be complex and stable over time, but we also demonstrate that the use of semi-resistant potato cultivars triggers a rapid shift in the mitochondrial haplotype associated with increased virulence.

CONCLUSIONS

Mitochondrial genomic variation shows that S. endobioticum has been introduced into Europe multiple times, that several pathotypes emerged multiple times, and that isolates represent communities of different genotypes. Our study represents the most comprehensive dataset of chytrid mitogenomes, which provides new insights into the extraordinary dynamics and evolution of mitochondrial genomes involving linearization, expansion and reshuffling.

Practical evaluation of 11 de novo assemblers in metagenome assembly

Next Generation Sequencing (NGS) technologies are revolutionizing the field of biology and metagenomic-based research. Since the volume of metagenomic data is typically very large, De novo metagenomic assembly can be effectively used to reduce the total amount of data and enhance quality of downstream analysis, such as annotation and binning. Although, there are many freely available assemblers, but selecting one suitable for a specific goal can be highly challenging. In this study, the performance of 11 well-known assemblers was evaluated in the assembly of three different metagenomes. The results obtained show that metaSPAdes is the best assembler and Megahit is a good choice for conservative assembly strategy. In addition, this research provides useful information regarding the pros and cons of each assembler and the effect of read length on assembly, thereby helping scholars to select the optimal assembler based on their objectives.

Genome Sequence of the Pseudomonas protegens Phage ΦGP100

We report here the complete annotated genome sequence of ΦGP100, a lytic bacteriophage of the Podoviridae family. ΦGP100 was isolated from rhizosphere soil in Switzerland and infects specifically strains of Pseudomonas protegens that are known for their plant-beneficial activities. Phage ΦGP100 has a 50,547-bp genome with 76 predicted open reading frames.

Genome comparison of different Zymomonas mobilis strains provides insights on conservation of the evolution

Zymomonas mobilis has the special Entner-Doudoroff (ED) pathway and it has excellent industrial characteristics, including low cell mass formation, high-specific productivity,ethanol yield, notable ethanol tolerance and wide pH range, a relatively small genome size. In this study, the genome sequences of NRRL B-14023 and NRRL B-12526 were sequenced and compared with other strains to explore their evolutionary relationships and the genetic basis of Z. mobilis. The comparative genomic analyses revealed that the 8 strains share a conserved core chromosomal backbone. ZM4, NRRL B-12526, NRRL B-14023, NCIMB 11163 and NRRL B-1960 share 98% sequence identity across the whole genome sequences. Highly similar plasmids and CRISPR repeats were detected in these strains. A whole-genome phylogenetic tree of the 8 strains indicated that NRRL B-12526, NRRL B-14023 and ATCC 10988 had a close evolutionary relationship with the strain ZM4. Furthermore, strains ATCC29191 and ATCC29192 had distinctive CRISPR with a far distant relationship. The size of the pan-genome was 1945 genes, including 1428 core genes and 517 accessory genes. The genomes of Z. mobilis were highly conserved; particularly strains ZM4, NRRL B-12526, NRRL B-14023, NCIMB 11163 and NRRL B-1960 had a close genomic relationship. This comparative study of Z. mobilis presents a foundation for future functional analyses and applications.

Genomic Structural Variations Affecting Virulence During Clonal Expansion of Pseudomonas syringae pv. actinidiae Biovar 3 in Europe

Pseudomonas syringae pv. actinidiae (Psa) biovar 3 caused pandemic bacterial canker of Actinidia chinensis and Actinidia deliciosa since 2008. In Europe, the disease spread rapidly in the kiwifruit cultivation areas from a single introduction. In this study, we investigated the genomic diversity of Psa biovar 3 strains during the primary clonal expansion in Europe using single molecule real-time (SMRT), Illumina and Sanger sequencing technologies. We recorded evidences of frequent mobilization and loss of transposon Tn6212, large chromosome inversions, and ectopic integration of IS sequences (remarkably ISPsy31, ISPsy36, and ISPsy37). While no phenotype change associated with Tn6212 mobilization could be detected, strains CRAFRU 12.29 and CRAFRU 12.50 did not elicit the hypersensitivity response (HR) on tobacco and eggplant leaves and were limited in their growth in kiwifruit leaves due to insertion of ISPsy31 and ISPsy36 in the hrpS and hrpR genes, respectively, interrupting the hrp cluster. Both strains had been isolated from symptomatic plants, suggesting coexistence of variant strains with reduced virulence together with virulent strains in mixed populations. The structural differences caused by rearrangements of self-genetic elements within European and New Zealand strains were comparable in number and type to those occurring among the European strains, in contrast with the significant difference in terms of nucleotide polymorphisms. We hypothesize a relaxation, during clonal expansion, of the selection limiting the accumulation of deleterious mutations associated with genome structural variation due to transposition of mobile elements. This consideration may be relevant when evaluating strategies to be adopted for epidemics management.

A Comprehensive Study of De Novo Genome Assemblers: Current Challenges and Future Prospective

BACKGROUND

Current advancements in next-generation sequencing technology have made possible to sequence whole genome but assembling a large number of short sequence reads is still a big challenge. In this article, we present the comparative study of seven assemblers, namely, ABySS, Velvet, Edena, SGA, Ray, SSAKE, and Perga, using prokaryotic and eukaryotic paired-end as well as single-end data sets from Illumina platform.

RESULTS

Results showed that in case of single-end data sets, Velvet and ABySS outperformed in all the seven assemblers with comparatively low assembling time and high genome fraction. Velvet consumed the least amount of memory than any other assembler. In case of paired-end data sets, Velvet consumed least amount of time and produced high genome fraction after ABySS and Ray. In terms of low memory usage, SGA and Edena outperformed in all the assemblers. Ray also showed good genome fraction; however, extremely high assembling time consumed by the Ray might make it prohibitively slow on larger data sets of single and paired-end data.

CONCLUSIONS

Our comparison study will provide assistance to the scientists for selecting the suitable assembler according to their data sets and will also assist the developers to upgrade or develop a new assembler for de novo assembling.

Inferences on specificity recognition at the Malus×domestica gametophytic self-incompatibility system

In Malus × domestica (Rosaceae) the product of each SFBB gene (the pollen component of the gametophytic self-incompatibility (GSI) system) of a S-haplotype (the combination of pistil and pollen genes that are linked) interacts with a sub-set of non-self S-RNases (the pistil component), but not with the self S-RNase. To understand how the Malus GSI system works, we identified 24 SFBB genes expressed in anthers, and determined their gene sequence in nine M. domestica cultivars. Expression of these SFBBs was not detected in the petal, sepal, filament, receptacle, style, stigma, ovary or young leaf. For all SFBBs (except SFBB15), identical sequences were obtained only in cultivars having the same S-RNase. Linkage with a particular S-RNase was further established using the progeny of three crosses. Such data is needed to understand how other genes not involved in GSI are affected by the S-locus region. To classify SFBBs specificity, the amino acids under positive selection obtained when performing intra-haplotypic analyses were used. Using this information and the previously identified S-RNase positively selected amino acid sites, inferences are made on the S-RNase amino acid properties (hydrophobicity, aromatic, aliphatic, polarity, and size), at these positions, that are critical features for GSI specificity determination.

A Novel Mechanism of Inactivating Antibacterial Nitro Compounds in the Human Pathogen Staphylococcus aureus by Overexpression of a NADH-Dependent Flavin Nitroreductase

Recently, the nitro-substituted bisquaternary bisnaphthalimides were reported to have substantial anti-infective activity against Gram-positive bacteria, including methicillin-resistant Staphylococcus aureus (MRSA). Here, we selected resistant S. aureus clones by cultivation in increasing concentrations of the most active compound, MT02. Interestingly, MT02-resistant variants induced a diffusible red color of the broth. Chromatographic and spectroscopic investigations revealed a stepwise reduction of the bisquaternary bisnaphthalimides' nitro groups to amino groups. The corresponding derivatives were completely inactive against staphylococci. RNA sequencing experiments revealed a strong overexpression of a novel oxidoreductase in MT02-resistant strains. Deletion mutants of this enzyme did not produce the red color and were not able to develop resistance against bisquaternary bisnaphthalimides. Biochemical reactions confirmed an NADH-dependent deactivation of the nitro-substituted compounds. Thus, this is the first report of a nitroreductase-based antibiotic resistance mechanism in the human pathogen S. aureus.

The Novel Phages phiCD5763 and phiCD2955 Represent Two Groups of Big Plasmidial Siphoviridae Phages of Clostridium difficile

Until recently, Clostridium difficile phages were limited to Myoviruses and Siphoviruses of medium genome length (32–57 kb). Here we report the finding of phiCD5763, a Siphovirus with a large extrachromosomal circular genome (132.5 kb, 172 ORFs) and a large capsid (205.6 ± 25.6 nm in diameter) infecting MLST Clade 1 strains of C. difficile. Two subgroups of big phage genomes similar to phiCD5763 were identified in 32 NAP_CR1/RT012/ST-54 C. difficile isolates from Costa Rica and in whole genome sequences (WGS) of 41 C. difficile isolates of Clades 1, 2, 3, and 4 from Canada, USA, UK, Belgium, Iraq, and China. Through comparative genomics we discovered another putative big phage genome in a non-NAP_CR1 isolate from Costa Rica, phiCD2955, which represents other big phage genomes found in 130 WGS of MLST Clade 1 and 2 isolates from Canada, USA, Hungary, France, Austria, and UK. phiCD2955 (131.6 kb, 172 ORFs) is related to a previously reported C. difficile phage genome, phiCD211/phiCDIF1296T. Detailed genome analyses of phiCD5763, phiCD2955, phiCD211/phiCDIF1296T, and seven other putative C. difficile big phage genome sequences of 131–136 kb reconstructed from publicly available WGS revealed a modular gene organization and high levels of sequence heterogeneity at several hotspots, suggesting that these genomes correspond to biological entities undergoing recombination. Compared to other C. difficile phages, these big phages have unique predicted terminase, capsid, portal, neck and tail proteins, receptor binding proteins (RBPs), recombinases, resolvases, primases, helicases, ligases, and hypothetical proteins. Moreover, their predicted gene load suggests a complex regulation of both phage and host functions. Overall, our results indicate that the prevalence of C. difficile big bacteriophages is more widespread than realized and open new avenues of research aiming to decipher how these viral elements influence the biology of this emerging pathogen.

Survey of gene splicing algorithms based on reads

Gene splicing is the process of assembling a large number of unordered short sequence fragments to the original genome sequence as accurately as possible. Several popular splicing algorithms based on reads are reviewed in this article, including reference genome algorithms and de novo splicing algorithms (Greedy-extension, Overlap-Layout-Consensus graph, De Bruijn graph). We also discuss a new splicing method based on the MapReduce strategy and Hadoop. By comparing these algorithms, some conclusions are drawn and some suggestions on gene splicing research are made.

Evidence of nuclei-encoded spliceosome mediating splicing of mitochondrial RNA

Mitochondria are thought to have originated as free-living prokaryotes. Mitochondria organelles have small circular genomes with substantial structural and genetic similarity to bacteria. Contrary to the prevailing concept of intronless mitochondria, here we present evidence that mitochondrial RNA transcripts (mtRNA) are not limited to policystronic molecules, but also processed as nuclei-like transcripts that are differentially spliced and expressed in a cell-type specific manner. The presence of canonical splice sites in the mtRNA introns and of core components of the nuclei-encoded spliceosome machinery within the mitochondrial organelle suggest that nuclei-encoded spliceosome can mediate splicing of mtRNA.

Development and characterization of genomic SSR markers for Anneslea fragrans (Pentaphylacaceae)

PREMISE OF THE STUDY

The genus Anneslea (Pentaphylacaceae) contains four species and six varieties, most of which are locally endemic. Here, simple sequence repeat (SSR) markers were developed for the conservation of these species.

METHODS AND RESULTS

The genome of A. fragrans was sequenced and de novo assembled into 445,162 contigs, of which 30,409 SSR loci were detected. Primers for 100 SSR loci were validated with PCR amplification in three populations of A. fragrans. Seventy-nine loci successfully amplified, and 30 were polymorphic. The mean number of alleles, observed heterozygosity, and expected heterozygosity were 7.01 ± 1.60, 0.817 ± 0.241, and 0.796 ± 0.145, respectively. Most primers could be amplified in Ternstroemia gymnanthera, T. kwangtungensis, and Cleyerapachyphylla.

CONCLUSIONS

Our study demonstrated that shotgun genome sequencing is an efficient way to develop genomic SSR markers for nonmodel species. These genomic SSR loci will be valuable in population genetic studies in Anneslea and its relatives.

Whole-Genome Sequencing of Bacterial Pathogens: the Future of Nosocomial Outbreak Analysis

Outbreaks of multidrug-resistant bacteria present a frequent threat to vulnerable patient populations in hospitals around the world. Intensive care unit (ICU) patients are particularly susceptible to nosocomial infections due to indwelling devices such as intravascular catheters, drains, and intratracheal tubes for mechanical ventilation. The increased vulnerability of infected ICU patients demonstrates the importance of effective outbreak management protocols to be in place. Understanding the transmission of pathogens via genotyping methods is an important tool for outbreak management. Recently, whole-genome sequencing (WGS) of pathogens has become more accessible and affordable as a tool for genotyping. Analysis of the entire pathogen genome via WGS could provide unprecedented resolution in discriminating even highly related lineages of bacteria and revolutionize outbreak analysis in hospitals. Nevertheless, clinicians have long been hesitant to implement WGS in outbreak analyses due to the expensive and cumbersome nature of early sequencing platforms. Recent improvements in sequencing technologies and analysis tools have rapidly increased the output and analysis speed as well as reduced the overall costs of WGS. In this review, we assess the feasibility of WGS technologies and bioinformatics analysis tools for nosocomial outbreak analyses and provide a comparison to conventional outbreak analysis workflows. Moreover, we review advantages and limitations of sequencing technologies and analysis tools and present a real-world example of the implementation of WGS for antimicrobial resistance analysis. We aimed to provide health care professionals with a guide to WGS outbreak analysis that highlights its benefits for hospitals and assists in the transition from conventional to WGS-based outbreak analysis.

Draft Genome Sequence of Klebsiella pneumoniae OK8, a Multidrug-Resistant Mouse and Human Pathogen

We report here the draft genome sequence of Klebsiella pneumoniae OK8, a multidrug-resistant strain which was isolated in 1976 from a human and is known to be a mouse pathogen.

Complete genome sequences of two novel autographiviruses infecting a bacterium from the Pseudomonas fluorescens group

In this paper, we describe two independent isolates of a new member of the subfamily Autographivirinae, Pseudomonas phage KNP. The type strain (KNP) has a linear, 40,491-bp-long genome with GC content of 57.3%, and 50 coding DNA sequences (CDSs). The genome of the second strain (WRT) contains one CDS less, encodes a significantly different tail fiber protein and is shorter (40,214 bp; GC content, 57.4%). Phylogenetic analysis indicates that both KNP and WRT belong to the genus T7virus. Together with genetically similar Pseudomonas phages (gh-1, phiPSA2, phiPsa17, PPPL-1, shl2, phi15, PPpW-4, UNO-SLW4, phiIBB-PF7A, Pf-10, and Phi-S1), they form a divergent yet coherent group that stands apart from the T7-like viruses (sensu lato). Analysis of the diversity of this group and its relatedness to other members of the subfamily Autographivirinae led us to the conclusion that this group might be considered as a candidate for a new genus.

Transcriptome analysis of primary aldosteronism in adrenal glands and controls

Primary aldosteronism (PA) is the most common form of endocrine hypertension. This study was to investigate the gene expression profile in PA adrenal glands and normal controls using RNA-Sequencing. By performing transcriptome analyses for 3 PA adrenal glands and 3 controls on Illumina platform, we identified 1,093 transcripts as significantly differently expressed genes (DEGs), which provided clues for further study of these transcript changes during PA pathogenesis. Further, Gene Set Enrichment Analysis (GSEA) identified 35 significant Kyoto Encyclopedia of Genes and Genomes (KEGG) biological pathways, including 'ribosome', 'oxidative phosphorylation', 'histidine metabolism', 'xenobiotics metabolism by Cytochrome P450', 'drug metabolism by Cytochrome P450', 'tyrosine metabolism' and 'glutathione metabolism'. In summary, we identified novel genes that are associated with PA phenotype, as well as differently regulated biological pathways relating to protein synthesis, energy acquisition and metabolism. Our study provides new candidates for further elucidation of the molecular mechanisms underlying PA pathogenesis.

Draft Genome Sequence of a Clinical Isolate of Streptococcus mutans Strain HM

We report the draft genome sequence of Streptococcus mutans strain HM isolated from a 4-year-old girl with infective endocarditis. The genomics information will provide information on the genetic diversity and virulence potential of S. mutans strain HM.

Draft Genome Sequence of the Nitrogen-Fixing Rhizobium sullae Type Strain IS123T Focusing on the Key Genes for Symbiosis with its Host Hedysarum coronarium L

The prominent feature of rhizobia is their molecular dialogue with plant hosts. Such interaction is enabled by the presence of a series of symbiotic genes encoding for the synthesis and export of signals triggering organogenetic and physiological responses in the plant. The genome of the Rhizobium sullae type strain IS123^T nodulating the legume Hedysarum coronarium, was sequenced and resulted in 317 scaffolds for a total assembled size of 7,889,576 bp. Its features were compared with those of genomes from rhizobia representing an increasing gradient of taxonomical distance, from a conspecific isolate (Rhizobium sullae WSM1592), to two congeneric cases (Rhizobium leguminosarum bv. viciae and Rhizobium etli) and up to different genera within the legume-nodulating taxa. The host plant is of agricultural importance, but, unlike the majority of other domesticated plant species, it is able to survive quite well in the wild. Data showed that that the type strain of R. sullae, isolated from a wild host specimen, is endowed with a richer array of symbiotic genes in comparison to other strains, species or genera of rhizobia that were rescued from domesticated plant ecotypes. The analysis revealed that the bacterium by itself is incapable of surviving in the extreme conditions that its host plant can tolerate. When exposed to drought or alkaline condition, the bacterium depends on its host to survive. Data are consistent with the view of the plant phenotype as the primary factor enabling symbiotic nitrogen fixing bacteria to survive in otherwise limiting environments.

Draft Genome Sequence of a Virulent Strain of Pasteurella Multocida Isolated From Alpaca

Pasteurella multocida is one of the most frequently isolated bacteria in acute pneumonia cases, being responsible for high mortality rates in Peruvian young alpacas, with consequent social and economic costs. Here we report the genome sequence of P. multocida strain UNMSM, isolated from the lung of an alpaca diagnosed with pneumonia, in Peru. The genome consists of 2,439,814 base pairs assembled into 82 contigs and 2,252 protein encoding genes, revealing the presence of known virulence-associated genes (ompH, ompA, tonB, tbpA, nanA, nanB, nanH, sodA, sodC, plpB and toxA). Further analysis could provide insights about bacterial pathogenesis and control strategies of this disease in Peruvian alpacas.

Insights into enterotoxigenic Escherichia coli diversity in Bangladesh utilizing genomic epidemiology

Enterotoxigenic Escherichia coli (ETEC) cause more than 500,000 deaths each year in the developing world and are characterized on a molecular level by the presence of genes that encode the heat-stable (ST) and/or heat-labile (LT) enterotoxins, as well as surface structures, known as colonization factors (CFs). Genome sequencing and comparative genomic analyses of 94 previously uncharacterized ETEC isolates demonstrated remarkable genomic diversity, with 28 distinct sequence types identified in three phylogenomic groups. Interestingly, there is a correlation between the genomic sequence type and virulence factor profiles based on prevalence of the isolate, suggesting that there is an optimal combination of genetic factors required for survival, virulence and transmission in the most successful clones. A large-scale BLAST score ratio (LS-BSR) analysis was further applied to identify ETEC-specific genomic regions when compared to non-ETEC genomes, as well as genes that are more associated with clinical presentations or other genotypic markers. Of the strains examined, 21 of 94 ETEC isolates lacked any previously identified CF. Homology searches with the structural subunits of known CFs identified 6 new putative CF variants. These studies provide a roadmap to exploit genomic analyses by directing investigations of pathogenesis, virulence regulation and vaccine development.

Monitoring of Newcastle disease virus in environmental samples

Newcastle disease virus (NDV) is an important pathogen in poultry. Waterfowl and a number of other avian species serve as the host for NDV. Severity of the disease is variable and infected animals mainly develop respiratory and neurological symptoms. Outbreaks of NDV in poultry are recorded regularly in the Republic of Kazakhstan despite the widespread use of vaccines. Here we present evidence that nucleic acid found in open water bodies in Kazakhstan can be detected by means of next-generation sequencing and belongs to at least three distinct genotypes of NDV.

A scalable and memory-efficient algorithm for de novo transcriptome assembly of non-model organisms

Background

With increased availability of de novo assembly algorithms, it is feasible to study entire transcriptomes of non-model organisms. While algorithms are available that are specifically designed for performing transcriptome assembly from high-throughput sequencing data, they are very memory-intensive, limiting their applications to small data sets with few libraries.

Results

We develop a transcriptome assembly algorithm that recovers alternatively spliced isoforms and expression levels while utilizing as many RNA-Seq libraries as possible that contain hundreds of gigabases of data. New techniques are developed so that computations can be performed on a computing cluster with moderate amount of physical memory.

Conclusions

Our strategy minimizes memory consumption while simultaneously obtaining comparable or improved accuracy over existing algorithms. It provides support for incremental updates of assemblies when new libraries become available.

Draft Genome Sequence of the Nylon Oligomer-Degrading Bacterium Arthrobacter sp. Strain KI72

We report here the 4.6-Mb genome sequence of a nylon oligomer-degrading bacterium, Arthrobacter sp. strain KI72. The draft genome sequence of strain KI72 consists of 4,568,574 bp, with a G+C content of 63.47%, 4,372 coding sequences (CDSs), 54 tRNAs, and six rRNAs.

Genome Sequence of Arcobacter sp. Strain LA11, Isolated from the Abalone Haliotis discus

Arcobacter sp. strain LA11 was isolated from the gut of the abalone Haliotis discus Here, we present the annotation and analysis of the draft genome of this strain, which is involved in nitrogen metabolism.

Analysis of de novo sequencing and transcriptome assembly and lignocellulolytic enzymes gene expression of Coriolopsis gallica HTC

White-rot basidiomycete Coriolopsis gallica HTC is one of the main biodegraders of poplar. In our previous study, we have shown the strong capacity of C. gallica HTC to degrade lignocellulose. In this study, equal amounts of total RNA fromC. Gallica HTC cultures grown in different conditions were pooled together. Illumina paired-end RNA sequencing was performed, and 13.2 million 90-bp paired-end reads were generated. We chose the Merged Assembly of Oases data-set for the following blast searches and gene ontology analyses. The reads were assembled de novo into 28,034 transcripts (≥ 100 bp) using combined assembly strategy MAO. The transcripts were annotated using Blast2GO. In all, 18,810 transcripts (≥100 bp) achieved BLASTX hits, of which, 7048 transcripts had GO term and 2074 had ECs. The expression level of 11 lignocellulolytic enzyme genes from the assembled C. gallica HTC transcriptome were detected by real-time quantitative polymerase chain reaction. The results showed that expression levels of these genes were affected by carbon source and nitrogen source at the level of transcription. The current abundant transcriptome data allowed the identification of many new transcripts in C. gallica HTC. Data provided here represent the most comprehensive and integrated genomic resources for cloning and identifying genes of interest from C. gallica HTC. Characterization of C. gallica HTC transcriptome provides an effective tool to understand mechanisms underlying cellular and molecular functions of C. gallica HTC.

Genome Sequence of Desulfurella amilsii Strain TR1 and Comparative Genomics of Desulfurellaceae Family

The acidotolerant sulfur reducer Desulfurella amilsii was isolated from sediments of Tinto River, an extremely acidic environment. Its ability to grow in a broad range of pH and to tolerate certain heavy metals offers potential for metal recovery processes. Here we report its high-quality draft genome sequence and compare it to the available genome sequences of other members of Desulfurellaceae family: D. acetivorans. D. multipotens, Hippea maritima. H. alviniae, H. medeae, and H. jasoniae. For most species, pairwise comparisons for average nucleotide identity (ANI) and in silico DNA-DNA hybridization (DDH) revealed ANI values from 67.5 to 80% and DDH values from 12.9 to 24.2%. D. acetivorans and D. multipotens, however, surpassed the estimated thresholds of species definition for both DDH (98.6%) and ANI (88.1%). Therefore, they should be merged to a single species. Comparative analysis of Desulfurellaceae genomes revealed different gene content for sulfur respiration between Desulfurella and Hippea species. Sulfur reductase is only encoded in D. amilsii, in which it is suggested to play a role in sulfur respiration, especially at low pH. Polysulfide reductase is only encoded in Hippea species; it is likely that this genus uses polysulfide as electron acceptor. Genes encoding thiosulfate reductase are present in all the genomes, but dissimilatory sulfite reductase is only present in Desulfurella species. Thus, thiosulfate respiration via sulfite is only likely in this genus. Although sulfur disproportionation occurs in Desulfurella species, the molecular mechanism behind this process is not yet understood, hampering a genome prediction. The metabolism of acetate in Desulfurella species can occur via the acetyl-CoA synthetase or via acetate kinase in combination with phosphate acetyltransferase, while in Hippea species, it might occur via the acetate kinase. Large differences in gene sets involved in resistance to acidic conditions were not detected among the genomes. Therefore, the regulation of those genes, or a mechanism not yet known, might be responsible for the unique ability of D. amilsii. This is the first report on comparative genomics of sulfur-reducing bacteria, which is valuable to give insight into this poorly understood metabolism, but of great potential for biotechnological purposes and of environmental significance.