AlignGraph: algorithm for secondary de novo genome assembly guided by closely related references

Genomic study of Acinetobacter baumannii strains co-harboring bla OXA-58 and bla NDM-1 reveals a large multidrug-resistant plasmid encoding these carbapenemases in Brazil

Introduction

Acinetobacter baumannii contributes significantly to the global issue of multidrug-resistant (MDR) nosocomial infections. Often, these strains demonstrate resistance to carbapenems (MDR-CRAB), the first-line treatment for infections instigated by MDR A. baumannii. Our study focused on the antimicrobial susceptibility and genomic sequences related to plasmids from 12 clinical isolates of A. baumannii that carry both the blaOXA-58 and bla NDM-1 carbapenemase genes.

Methods

Whole-genome sequencing with long-read technology was employed for the characterization of an A. baumannii plasmid that harbors the bla OXA-58 and blaNDM-1 genes. The location of the bla OXA-58 and bla NDM-1 genes was confirmed through Southern blot hybridization assays. Antimicrobial susceptibility tests were conducted, and molecular characterization was performed using PCR and PFGE.

Results

Multilocus Sequence Typing analysis revealed considerable genetic diversity among bla OXA-58 and bla NDM-1 positive strains in Brazil. It was confirmed that these genes were located on a plasmid larger than 300 kb in isolates from the same hospital, which also carry other antimicrobial resistance genes. Different genetic contexts were observed for the co-occurrence of these carbapenemase-encoding genes in Brazilian strains.

Discussion

The propagation of bla OXA-58 and bla NDM-1 genes on the same plasmid, which also carries other resistance determinants, could potentially lead to the emergence of bacterial strains resistant to multiple classes of antimicrobials. Therefore, the characterization of these strains is of paramount importance for monitoring resistance evolution, curbing their rapid global dissemination, averting outbreaks, and optimizing therapy.

GeneMiner: A tool for extracting phylogenetic markers from next-generation sequencing data

The advancement of next-generation sequencing (NGS) technologies has been revolutionary for the field of evolutionary biology. This technology has led to an abundance of available genomes and transcriptomes for researchers to mine. Specifically, researchers can mine for various types of molecular markers that are vital for phylogenetic, evolutionary and ecological studies. Numerous tools have been developed to extract these molecular markers from NGS data. However, due to an insufficient number of well-annotated reference genomes for non-model organisms, it remains challenging to obtain these markers accurately and efficiently. Here, we present GeneMiner, an improved and expanded version of our previous tool, Easy353. GeneMiner combines the reference-guided de Bruijn graph assembly with seed self-discovery and greedy extension. Additionally, it includes a verification step using a parameter-bootstrap method to reduce the pitfalls associated with using a relatively distant reference. Our results, using both experimental and simulation data, showed GeneMiner can accurately acquire phylogenetic molecular markers for plants using transcriptomic, genomic and other NGS data. GeneMiner is designed to be user-friendly, fast and memory-efficient. Further, it is compatible with Linux, Windows and macOS. All source codes are publicly available on GitHub (https://github.com/sculab/GeneMiner) and Gitee (https://gitee.com/sculab/GeneMiner) for easy accessibility and transparency.

The Utilization of Reference-Guided Assembly and In Silico Libraries Improves the Draft Genome of Clarias batrachus and Culter alburnus

Long-read sequencing technologies can generate highly contiguous genome assemblies compared to short-read methods. However, their higher cost often poses a significant barrier. To address this, we explore the utilization of mapping-based genome assembly and reference-guided assembly as cost-effective alternative approaches. We assess the efficacy of these approaches in improving the contiguity of Clarias batrachus and Culter alburnus draft genomes. Our findings demonstrate that employing an iterative mapping strategy leads to a reduction in assembly errors. Specifically, after three iterations, the Mismatches per 100 kbp value for the C. batrachus genome decreased from 2447.20 to 2432.67, reaching a minimum of 2422.67 after two iterations. Additionally, the N50 value for the C. batrachus genome increased from 362,143 to 1,315,126 bp, with a maximum of 1,315,403 bp after two iterations. Furthermore, we achieved Mismatches per 100 kbp values of 3.70 for the reference-guided assembly of C. batrachus and 0.34 for C. alburnus. Correspondingly, the N50 value for the C. batrachus and C. alburnus genomes increased from 362,143 bp and 3,686,385 bp to 2,026,888 bp and 43,735,735 bp, respectively. Finally, we successfully utilized the improved C. batrachus and C. alburnus genomes to compare genome studies using the combined approach of Ragout and Ragtag. Through a comprehensive comparative analysis of mapping-based and reference-guided genome assembly methods, we shed light on the specific contributions of reference-guided assembly in reducing assembly errors and improving assembly continuity and integrity. These advancements establish reference-guided assembly and the utilization of in silico libraries as a promising and suitable approach for comparative genomics studies.

Genomic analysis of phylogroup D Escherichia coli strains using novel de-novo reference-based guided assembly

Escherichia coli are highly diverse bacteria with different pathogenic types, serotypes and phylogenetic types/phylotypes. In recent years, infections with E. coli have increased worldwide and so has the emergence of antibiotic resistant strains. In the present study we have assembled, annotated and analysed genome sequences of three strains of the phylogroup D of E. coli. These strains were isolated from the river Yamuna, a prominent anthropogenic urban river of northern India. These strains showed varied antibiotic susceptibilities, one was susceptible to all the antibiotics tested except ampicillin while of the other two, one was multi-β-lactam resistant and the other was multi-drug resistant (resistant to multiple β-lactams, fluoroquinolones and kanamycin). The short-sequence reads were assembled into contigs using the de-novo approach and further, scaffolding of contigs was performed by using the best reference genome for a particular isolate which resulted in a significant increase in the N50 value of each assembly. The bioinformatics assembly approach used in this study could be easily applied to study other bacterial genomes.

Optimization of the "in-silico" mate-pair method improves contiguity and accuracy of genome assembly

A combination of short-insert paired-ended and mate-pair libraries of large insert sizes is used as a standard method to generate genome assemblies with high contiguity. The third-generation sequencing techniques also are used to improve the quality of assembled genomes. However, both mate-pair libraries and the third-generation libraries require high-molecular-weight DNA, making the use of these libraries inappropriate for samples with only degraded DNA. An in silico method that generates mate-pair libraries using a reference genome was devised for the task of assembling target genomes. Although the contiguity and completeness of assembled genomes were significantly improved by this method, a high level of errors manifested in the assembly, further to which the methods for using reference genomes, was not optimized. Here, we tested different strategies for using reference genomes to generate in silico mate-pairs. The results showed that using a closely related reference genome from the same genus was more effective than using divergent references. Conservation of in silico mate-pairs by comparing two references and using those to guide genome assembly reduced the number of misassemblies (18.6%-46.1%) and increased the contiguity of assembled genomes (9.7%-70.7%), while maintaining gene completeness at a level that was either similar or marginally lower than that obtained via the current method. Finally, we developed a pipeline of the optimized in silico method and compared it with another reference-guided assembler, RagTag. We found that RagTag produced longer scaffolds (17.8 Mbp vs 3.0 Mbp), but resulted in a much higher misassembly rate (85.68%) than our optimized in silico mate-pair method. This optimized in silico pipeline developed in this study should facilitate further studies on genomics, population genetics, and conservation of endangered species.

Easy353: A Tool to Get Angiosperms353 Genes for Phylogenomic Research

The Angiosperms353 gene set (AGS) consists of a set of 353 universal low-copy nuclear genes that were selected by examining more than 600 angiosperm species. These genes can be used for phylogenetic studies and population genetics at multiple taxonomic scales. However, current pipelines are not able to recover Angiosperms353 genes efficiently and accurately from high-throughput sequences. Here, we developed Easy353, a reference-guided assembly tool to recover the AGS from high-throughput sequencing (HTS) data (including genome skimming, RNA-seq, and target enrichment). Easy353 is an open-source user-friendly assembler for diverse types of high-throughput data. It has a graphical user interface and a command-line interface that is compatible with all widely-used computer systems. Evaluations, based on both simulated and empirical data, suggest that Easy353 yields low rates of assembly errors.

Validation of reference-assisted assembly using existing and novel Heliothine genomes

Chloridea subflexa and Chloridea virescens are a pair of closely related noctuid species exhibiting pheromone-based sexual isolation and divergent host plant preferences. We produced a novel Illumina short read C. subflexa genome assembly and an improved C. virescens genome assembly, which offer opportunities to study the genomic basis for evolutionarily important traits in this lepidopteran family with few genomic resources. We then examined the feasibility of reference-assisted assembly, an approach that leverages existing high quality genomic resources for genome improvement in closely related taxa and applied it to our Heliothine genomes. Our work demonstrates that reference-assisted assembly has the potential to enhance contiguity and completeness of existing insect genomic resources with minimal additional laboratory costs. We conclude by discussing both the potential and pitfalls of reference-assisted assembly according to the intended downstream assembly application.

The Visayan warty pig (Sus cebifrons) genome provides insight into chromosome evolution and sensory adaptation in pigs

It is largely unknown how mammalian genomes evolve under rapid speciation and environmental adaptation. An excellent model for understanding fast evolution is provided by the genus Sus, which diverged relatively recently and lacks post-zygotic isolation. Here, we present a high-quality reference genome of the Visayan warty pig, which is specialized to a tropical island environment. Comparing the genome sequences and chromatin contact maps of the Visayan warty pig (Sus cebifrons) and domestic pig (Sus scrofa), we characterized the dynamics of chromosomal structure evolution during Sus speciation, revealing the similar chromosome conformation as the potential biological mechanism of frequent post-divergence hybridization among Suidae. We further investigated the different signatures of adaptive selection and domestication in Visayan warty pig and domestic pig with specific emphasize on the evolution of olfactory and gustatory genes, elucidating higher olfactory diversity in Visayan warty pig and positive and relaxed evolution of bitter and fat taste receptors, respectively, in domestic pig. Our comprehensive evolutionary and comparative genome analyses provide insight into the dynamics of genomes and how these change over relative short evolutionary times, as well as how these genomic differences encode for differences in the phenotypes.

HGGA: hierarchical guided genome assembler

BACKGROUND

De novo genome assembly typically produces a set of contigs instead of the complete genome. Thus additional data such as genetic linkage maps, optical maps, or Hi-C data is needed to resolve the complete structure of the genome. Most of the previous work uses the additional data to order and orient contigs.

RESULTS

Here we introduce a framework to guide genome assembly with additional data. Our approach is based on clustering the reads, such that each read in each cluster originates from nearby positions in the genome according to the additional data. These sets are then assembled independently and the resulting contigs are further assembled in a hierarchical manner. We implemented our approach for genetic linkage maps in a tool called HGGA.

CONCLUSIONS

Our experiments on simulated and real Pacific Biosciences long reads and genetic linkage maps show that HGGA produces a more contiguous assembly with less contigs and from 1.2 to 9.8 times higher NGA50 or N50 than a plain assembly of the reads and 1.03 to 6.5 times higher NGA50 or N50 than a previous approach integrating genetic linkage maps with contig assembly. Furthermore, also the correctness of the assembly remains similar or improves as compared to an assembly using only the read data.

Genomic Analyses of Penicillium Species Have Revealed Patulin and Citrinin Gene Clusters and Novel Loci Involved in Oxylipin Production

Blue mold of apple is caused by several different Penicillium species, among which P. expansum and P. solitum are the most frequently isolated. P. expansum is the most aggressive species, and P. solitum is very weak when infecting apple fruit during storage. In this study, we report complete genomic analyses of three different Penicillium species: P. expansum R21 and P. crustosum NJ1, isolated from stored apple fruit; and P. maximae 113, isolated in 2013 from a flooded home in New Jersey, USA, in the aftermath of Hurricane Sandy. Patulin and citrinin gene cluster analyses explained the lack of patulin production in NJ1 compared to R21 and lack of citrinin production in all three strains. A Drosophila bioassay demonstrated that volatiles emitted by P. solitum SA and P. polonicum RS1 were more toxic than those from P. expansum and P. crustosum strains (R27, R11, R21, G10, and R19). The toxicity was hypothesized to be related to production of eight-carbon oxylipins. Putative lipoxygenase genes were identified in P. expansum and P. maximae strains, but not in P. crustosum. Our data will provide a better understanding of Penicillium spp. complex secondary metabolic capabilities, especially concerning the genetic bases of mycotoxins and toxic VOCs.

Whole-genome assembly and annotation of northern wild rice, Zizania palustris L., supports a whole-genome duplication in the Zizania genus

Zizania palustris L. (northern wild rice, NWR) is an aquatic grass native to North America that is notable for its nutritious grain. This is an important species with ecological, cultural and agricultural significance, specifically in the Great Lakes region of the USA. Using flow cytometry, we first estimated the NWR genome size to be 1.8 Gb. Using long- and short-range sequencing, Hi-C scaffolding and RNA-seq data from eight tissues, we generated an annotated whole-genome de novo assembly of NWR. The assembly was 1.29 Gb in length, highly repetitive (approx. 76.0%) and contained 46 421 putative protein-coding genes. The expansion of retrotransposons within the genome and a whole-genome duplication (WGD) after the Zizania-Oryza speciation event have both led to an increase in the genome size of NWR in comparison with Oryza sativa L. and Zizania latifolia. Both events depict a genome rapidly undergoing change over a short evolutionary time. Comparative analyses revealed the conservation of large syntenic blocks between NWR and O. sativa, which were used to identify putative seed-shattering genes. Estimates of divergence times revealed that the Zizania genus diverged from Oryza approximately 26-30 million years ago (26-30 MYA), whereas NWR and Z. latifolia diverged from one another approximately 6-8 MYA. Comparative genomics confirmed evidence of a WGD in the Zizania genus and provided support that the event occurred prior to the NWR-Z. latifolia speciation event. This genome assembly and annotation provides a valuable resource for comparative genomics in the Oryzeae tribe and provides an important resource for future conservation and breeding efforts of NWR.

A Novel Polyvalent Bacteriophage vB_EcoM_swi3 Infects Pathogenic Escherichia coli and Salmonella enteritidis

A novel virulent bacteriophage vB_EcoM_swi3 (swi3), isolated from swine feces, lyzed 9% (6/65) of Escherichia coli and isolates 54% (39/72) of Salmonella enteritidis isolates, which were all clinically pathogenic multidrug-resistant strains. Morphological observation showed that phage swi3 belonged to the Myoviridae family with an icosahedral head (80 nm in diameter) and a contractile sheathed tail (120 nm in length). At the optimal multiplicity of infection of 1, the one-step growth analysis of swi3 showed a 25-min latent period with a burst size of 25-plaque-forming units (PFU)/infected cell. Phage swi3 remained stable both at pH 6.0–8.0 and at less than 50°C for at least 1 h. Genomic sequencing and bioinformatics analysis based on genomic sequences and the terminase large subunit showed that phage swi3 was a novel member that was most closely related to Salmonella phages and belonged to the Rosemountvirus genus. Phage swi3 harbored a 52-kb double-stranded DNA genome with 46.02% GC content. Seventy-two potential open reading frames were identified and annotated, only 15 of which had been assigned to functional genes. No gene associated with pathogenicity and virulence was identified. The effects of phage swi3 in treating pathologic E. coli infections in vivo were evaluated using a mouse model. The administration of a single intraperitoneal injection of swi3 (10⁶ PFU) at 2 h after challenge with the E. coli strain (serotype K88) (10⁸ colony-forming units) sufficiently protected all mice without toxic side effects. This finding highlighted that phage swi3 might be used as an effective antibacterial agent to prevent E. coli and S. enteritidis infection.

Characterization of a Novel Bacteriophage swi2 Harboring Two Lysins Can Naturally Lyse Escherichia coli

The novel virulent Siphoviridae bacteriophage swi2 was isolated from a pig farm, and its biological characteristics, genome architecture, and infection-related properties were characterized. Phage swi2 has a high titer of 1.01 × 10¹² PFU/mL with good tolerance to UV rays and remains stable in the pH range of 6–10 and at temperatures less than 50°C. One-step growth analysis revealed that phage swi2 had a 25 min latent period with a large burst size (1,000 PFU/cell). The biological characteristics indicated that swi2 had good host infectivity and effective lytic activities. The genome of phage swi2 is composed of 47,611 bp with a G + C content of 46.50%. Eighty-nine orfs were predicted, and only 18 of them have known functions. No virulence genes or drug resistance genes were found in the genome. Genome sequence comparison of phage swi2 showed that there were a total of 10 homologous phages in the database with low similarity (less than 92.51% nucleotide identity and 66% query coverage). The predicted host lysis-related genes of phage swi2 consist of one holin, two endolysins, and Rz/Rz1 equivalents. Antibacterial activity assays showed that both endolysins could naturally reduce the host Escherichia coli 51 titers by -1 log unit both in vitro and in vivo, EDTA showed no obvious synergistic action, and holin had no lytic effects on the host cell. These results provide necessary information for the development of antibiotic alternatives for the treatment of multidrug-resistant Escherichia coli infection.

Draft Genome Sequences of the Clinical Isolates Kp 23 and KPM 20

Klebsiella pneumoniae strains are capable of becoming resistant through multiple mechanisms. Here, we announce draft sequences for Kp 23, a clinical isolate with no plasmid-encoded β-lactamases, and KPM 20, a clinical isolate with no plasmid-encoded β-lactamases and no detectable OmpK35, OmpK36, or PhoE in the outer membrane.

Klebsiella pneumoniae strains are capable of becoming resistant through multiple mechanisms. Here, we announce draft sequences for Kp 23, a clinical isolate with no plasmid-encoded β-lactamases, and KPM 20, a clinical isolate with no plasmid-encoded β-lactamases and no detectable OmpK35, OmpK36, or PhoE in the outer membrane.

AlignGraph2: similar genome-assisted reassembly pipeline for PacBio long reads

Contigs assembled from the third-generation sequencing long reads are usually more complete than the second-generation short reads. However, the current algorithms still have difficulty in assembling the long reads into the ideal complete and accurate genome, or the theoretical best result [1]. To improve the long read contigs and with more and more fully sequenced genomes available, it could still be possible to use the similar genome-assisted reassembly method [2], which was initially proposed for the short reads making use of a closely related genome (similar genome) to the sequencing genome (target genome). The method aligns the contigs and reads to the similar genome, and then extends and refines the aligned contigs with the aligned reads. Here, we introduce AlignGraph2, a similar genome-assisted reassembly pipeline for the PacBio long reads. The AlignGraph2 pipeline is the second version of AlignGraph algorithm proposed by us but completely redesigned, can be inputted with either error-prone or HiFi long reads, and contains four novel algorithms: similarity-aware alignment algorithm and alignment filtration algorithm for alignment of the long reads and preassembled contigs to the similar genome, and reassembly algorithm and weight-adjusted consensus algorithm for extension and refinement of the preassembled contigs. In our performance tests on both error-prone and HiFi long reads, AlignGraph2 can align 5.7-27.2% more long reads and 7.3-56.0% more bases than some current alignment algorithm and is more efficient or comparable to the others. For contigs assembled with various de novo algorithms and aligned to similar genomes (aligned contigs), AlignGraph2 can extend 8.7-94.7% of them (extendable contigs), and obtain contigs of 7.0-249.6% larger N50 value and 5.2-87.7% smaller number of indels per 100 kbp (extended contigs). With genomes of decreased similarities, AlignGraph2 also has relatively stable performance. The AlignGraph2 software can be downloaded for free from this site: https://github.com/huangs001/AlignGraph2.

Letter to the Editor: Escherichia fergusonii Harboring IncHI2 Plasmid Containing mcr-1 Gene-A Novel Reservoir for Colistin Resistance in Brazil

Emergence of colistin-resistant bacteria harboring mobile colistin resistance genes (mcr genes) pose a threat for food-producing animals and humans. In this article, we aim to highlight the emergence of Escherichia fergusonii as an important new reservoir to mcr-1-harboring plasmid in poultry production. Three strains closely related were isolated from cloacal swabs. Their genome contains four plasmids, including a 182,869 bp IncHI2 plasmid harboring the colistin resistance gene mcr-1. These results will contribute to our understanding of plasmid-mediated mcr-1 gene presence and transmission in E. fergusonii.

Harnessing Pseudomonas protegens to Control Bacterial Panicle Blight of Rice

Bacterial panicle blight of rice is a seedborne disease caused by the bacterium Burkholderia glumae. This disease has affected rice production worldwide and its effects are likely to become more devastating with the continuous increase in global temperatures, especially during the growing season. The bacterium can cause disease symptoms in different tissues and at different developmental stages. In reproductive stages, the bacterium interferes with grain development in the panicles and, as a result, directly affects rice yield. Currently, there are no methods to control the disease because chemical control is not effective and completely resistant cultivars are not available. Thus, a promising approach is the use of antagonistic microorganisms. In this work, we identified one strain of Pseudomonas protegens and one strain of B. cepacia with antimicrobial activity against B. glumae in vitro and in planta. We further characterized the antimicrobial activity of P. protegens and found that this activity is associated with bacterial secretions. Cell-free secretions from P. protegens inhibited the growth of B. glumae in vitro and also prevented B. glumae from causing disease in rice. Although the specific molecules associated with these activities have not been identified, these findings suggest that the secreted fractions from P. protegens could be harnessed as biopesticides to control bacterial panicle blight of rice.

The functional determinants in the organization of bacterial genomes

Bacterial genomes are now recognized as interacting intimately with cellular processes. Uncovering organizational mechanisms of bacterial genomes has been a primary focus of researchers to reveal the potential cellular activities. The advances in both experimental techniques and computational models provide a tremendous opportunity for understanding these mechanisms, and various studies have been proposed to explore the organization rules of bacterial genomes associated with functions recently. This review focuses mainly on the principles that shape the organization of bacterial genomes, both locally and globally. We first illustrate local structures as operons/transcription units for facilitating co-transcription and horizontal transfer of genes. We then clarify the constraints that globally shape bacterial genomes, such as metabolism, transcription and replication. Finally, we highlight challenges and opportunities to advance bacterial genomic studies and provide application perspectives of genome organization, including pathway hole assignment and genome assembly and understanding disease mechanisms.

Plant Growth-Promoting Methylobacteria Selectively Increase the Biomass of Biotechnologically Relevant Microalgae

Microalgae, a diverse group of single-celled organisms exhibiting versatile traits, find broad applications in industry. However, high production costs require further efforts to optimize their production and to enhance biomass yields. In the present study, co-occurrence of algae and methylobacteria was observed when naturally occurring microalgae biofilms were subjected to 16S rRNA gene fragment amplicon sequencing. This bacterial group is so far less explored than other microalgae-associated bacteria in terms of mutualistic relationships that might be exploitable for biotechnological applications. In order to assess the potential of four plant growth-promoting strains from the genus Methylobacterium for increased algae biomass production, co-cultivation experiments were conducted with three industrially relevant microalgae (Chlorella vulgaris, Scenedesmus vacuolatus, and Haematococcus lacustris). For S. vacuolatus and H. lacustris, a significant increase in algal biomass formation of 1.3-fold to up to 14-fold was observed after 7 days of co-incubation. Visualization of mixed cultures using confocal laser scanning microscopy revealed a high abundance of methylobacteria in the phycosphere of H. lacustris and S. vacuolatus, visually attached to the algae’s surface forming a biofilm-like assemblage. Genome analyses revealed that features attributable to enhanced algal growth include genes involved in the synthesis of vitamins, siderophores and plant hormones. Our results provide evidence for the constructability of novel symbiotic algae-bacteria relationships with inter-kingdom supportive capacities, underlining the potential of microbial consortia as promising tool for sustainable biotechnology and agriculture.

Vanillin Production in Pseudomonas: Whole-Genome Sequencing of Pseudomonas sp. Strain 9.1 and Reannotation of Pseudomonas putida CalA as a Vanillin Reductase

Microbial degradation of lignin and its related aromatic compounds has great potential for the sustainable production of chemicals and bioremediation of contaminated soils. We previously isolated Pseudomonas sp. strain 9.1 from historical waste deposits (forming so-called fiber banks) released from pulp and paper mills along the Baltic Sea coast. The strain accumulated vanillyl alcohol during growth on vanillin, and while reported in other microbes, this phenotype is less common in wild-type pseudomonads. As the reduction of vanillin to vanillyl alcohol is an undesired trait in Pseudomonas strains engineered to accumulate vanillin, connecting the strain 9.1 phenotype with a genotype would increase the fundamental understanding and genetic engineering potential of microbial vanillin metabolism. The genome of Pseudomonas sp. 9.1 was sequenced and assembled. Annotation identified oxidoreductases with homology to Saccharomyces cerevisiae alcohol dehydrogenase ScADH6p, known to reduce vanillin to vanillyl alcohol, in both the 9.1 genome and the model strain Pseudomonas putida KT2440. Recombinant expression of the Pseudomonas sp. 9.1 FEZ21_09870 and P. putida KT2440 PP_2426 (calA) genes in Escherichia coli revealed that these open reading frames encode aldehyde reductases that convert vanillin to vanillyl alcohol, and that P. putida KT2440 PP_3839 encodes a coniferyl alcohol dehydrogenase that oxidizes coniferyl alcohol to coniferyl aldehyde (i.e., the function previously assigned to calA). The deletion of PP_2426 in P. putida GN442 engineered to accumulate vanillin resulted in a decrease in by-product (vanillyl alcohol) yield from 17% to ∼1%. Based on these results, we propose the reannotation of PP_2426 and FEZ21_09870 as areA and PP_3839 as calA-II IMPORTANCE Valorization of lignocellulose (nonedible plant matter) is of key interest for the sustainable production of chemicals from renewable resources. Lignin, one of the main constituents of lignocellulose, is a heterogeneous aromatic biopolymer that can be chemically depolymerized into a heterogeneous mixture of aromatic building blocks; those can be further converted by certain microbes into value-added aromatic chemicals, e.g., the flavoring agent vanillin. We previously isolated a Pseudomonas sp. strain with the (for the genus) unusual trait of vanillyl alcohol production during growth on vanillin. Whole-genome sequencing of the isolate led to the identification of a vanillin reductase candidate gene whose deletion in a recombinant vanillin-accumulating P. putida strain almost completely alleviated the undesired vanillyl alcohol by-product yield. These results represent an important step toward biotechnological production of vanillin from lignin using bacterial cell factories.

Mapping Ethanol Tolerance in Budding Yeast Reveals High Genetic Variation in a Wild Isolate

Ethanol tolerance, a polygenic trait of the yeast Saccharomyces cerevisiae, is the primary factor determining industrial bioethanol productivity. Until now, genomic elements affecting ethanol tolerance have been mapped only at low resolution, hindering their identification. Here, we explore the genetic architecture of ethanol tolerance, in the F6 generation of an Advanced Intercrossed Line (AIL) mapping population between two phylogenetically distinct, but phenotypically similar, S. cerevisiae strains (a common laboratory strain and a wild strain isolated from nature). Under ethanol stress, 51 quantitative trait loci (QTLs) affecting growth and 96 QTLs affecting survival, most of them novel, were identified, with high resolution, in some cases to single genes, using a High-Resolution Mapping Package of methodologies that provided high power and high resolution. We confirmed our results experimentally by showing the effects of the novel mapped genes: MOG1, MGS1, and YJR154W. The mapped QTLs explained 34% of phenotypic variation for growth and 72% for survival. High statistical power provided by our analysis allowed detection of many loci with small, but mappable effects, uncovering a novel “quasi-infinitesimal” genetic architecture. These results are striking demonstration of tremendous amounts of hidden genetic variation exposed in crosses between phylogenetically separated strains with similar phenotypes; as opposed to the more common design where strains with distinct phenotypes are crossed. Our findings suggest that ethanol tolerance is under natural evolutionary fitness-selection for an optimum phenotype that would tend to eliminate alleles of large effect. The study provides a platform for development of superior ethanol-tolerant strains using genome editing or selection.

Genetic and Genome Analyses Reveal Genetically Distinct Populations of the Bee Pathogen Nosema ceranae from Thailand

The recent global decline in Western honeybee (Apis mellifera) populations is of great concern for pollination and honey production worldwide. Declining honeybee populations are frequently infected by the microsporidian pathogen Nosema ceranae. This species was originally described in the Asiatic honeybee (Apis cerana), and its identification in global A. mellifera hives could result from a recent host transfer. Recent genome studies have found that global populations of this parasite are polyploid and that humans may have fueled their global expansion. To better understand N. ceranae biology, we investigated its genetic diversity within part of their native range (Thailand) and among different hosts (A. mellifera, A. cerana) using both PCR and genome-based methods. We find that Thai N. ceranae populations share many SNPs with other global populations and appear to be clonal. However, in stark contrast with previous studies, we found that these populations also carry many SNPs not found elsewhere, indicating that these populations have evolved in their current geographic location for some time. Our genome analyses also indicate the potential presence of diploidy within Thai populations of N. ceranae.

lptG contributes to changes in membrane permeability and the emergence of multidrug hypersusceptibility in a cystic fibrosis isolate of Pseudomonas aeruginosa

PURPOSE

In the lungs of cystic fibrosis patients, Pseudomonas aeruginosa is exposed to a myriad of antibiotics leading to alterations in antibiotic susceptibility. This study identifies mutations resulting in hypersusceptibility in isogenic mutants of a P. aeruginosa clinical isolate, PA34.

METHODS

PA34 was exposed to subinhibitory concentrations of doripenem or meropenem during growth to mid-log phase. Antibiotic susceptibility of surviving colonies was determined by agar dilution. Two carbapenem-resistant colonies hypersusceptible to non-carbapenem antibiotics were selected for further analysis. Antibiotic resistance gene expression was evaluated by RT-rtPCR and OprD production by SDS-PAGE. PA34 and isogenic mutants were evaluated with whole genome sequencing. Sequence variants were confirmed by Sanger sequencing, and cognate genes in eight carbapenem-resistant clinical isolates hypersusceptible to non-carbapenem antibiotics were sequenced. Lipopolysaccharide preparations of PA34 and hypersusceptible mutants were evaluated with ProQ-Emerald stain.

RESULTS

Isogenic mutants showed 4- to 8-fold MIC increase for imipenem, meropenem, and doripenem. However, they were hypersusceptible (≥4-fold MIC decrease) to aminoglycosides, fluoroquinolones, and non-carbapenem β-lactams. Expression of ampC or mex-opr efflux pumps was unchanged, but OprD production was decreased. Mutations causing Q86H AlgU and G77C LptG amino acid substitutions and nonsense mutations within OprD were observed in both mutants. Lipopolysaccharide modifications were observed between isogenic mutants and PA34. Non-synonymous mutations in LptF or LptG were observed in 6/8 hypersusceptible clinical isolates resistant to carbapenem antibiotics.

CONCLUSION

Evaluation of hypersusceptible mutants identified the association between lptG and a hypersusceptible phenotype. Modifications in lipopolysaccharide profiles suggests LptG modification interferes with lipopolysaccharide transport and contributes to hypersusceptibility.

Draft Genome Assemblies of Clinical Isolates of Klebsiella pneumoniae V9011662 and Enterobacter hormaechei Entb306

Enterobacter hormaechei and Klebsiella pneumoniae are pathogenic Enterobacteriaceae that have been associated with the spread of antibiotic resistance. Here, we report draft genome assemblies of an Enterobacter hormaechei clinical isolate and a multidrug-resistant clinical isolate of Klebsiella pneumoniae.

Enterobacter hormaechei and Klebsiella pneumoniae are pathogenic Enterobacteriaceae that have been associated with the spread of antibiotic resistance. Here, we report draft genome assemblies of an Enterobacter hormaechei clinical isolate and a multidrug-resistant clinical isolate of Klebsiella pneumoniae.

Kermit: linkage map guided long read assembly

BACKGROUND

With long reads getting even longer and cheaper, large scale sequencing projects can be accomplished without short reads at an affordable cost. Due to the high error rates and less mature tools, de novo assembly of long reads is still challenging and often results in a large collection of contigs. Dense linkage maps are collections of markers whose location on the genome is approximately known. Therefore they provide long range information that has the potential to greatly aid in de novo assembly. Previously linkage maps have been used to detect misassemblies and to manually order contigs. However, no fully automated tools exist to incorporate linkage maps in assembly but instead large amounts of manual labour is needed to order the contigs into chromosomes.

RESULTS

We formulate the genome assembly problem in the presence of linkage maps and present the first method for guided genome assembly using linkage maps. Our method is based on an additional cleaning step added to the assembly. We show that it can simplify the underlying assembly graph, resulting in more contiguous assemblies and reducing the amount of misassemblies when compared to de novo assembly.

CONCLUSIONS

We present the first method to integrate linkage maps directly into genome assembly. With a modest increase in runtime, our method improves contiguity and correctness of genome assembly.

Nearly Complete Genome Sequence of a Human Norovirus GII.P17-GII.17 Strain Isolated from Brazil in 2015

Human noroviruses are the most common cause of nonbacterial acute gastroenteritis worldwide. We report here the nearly complete genome sequence (7,551 nucleotides) of a human norovirus GII.P17-GII.17 strain detected in July 2015 in the stool sample from an adult with acute gastroenteritis in Brazil.

Human noroviruses are the most common cause of nonbacterial acute gastroenteritis worldwide. We report here the nearly complete genome sequence (7,551 nucleotides) of a human norovirus GII.P17-GII.17 strain detected in July 2015 in the stool sample from an adult with acute gastroenteritis in Brazil.

Improving draft genome contiguity with reference-derived in silico mate-pair libraries

Background

Contiguous genome assemblies are a highly valued biological resource because of the higher number of completely annotated genes and genomic elements that are usable compared to fragmented draft genomes. Nonetheless, contiguity is difficult to obtain if only low coverage data and/or only distantly related reference genome assemblies are available.

Findings

In order to improve genome contiguity, we have developed Cross-Species Scaffolding—a new pipeline that imports long-range distance information directly into the de novo assembly process by constructing mate-pair libraries in silico.

Conclusions

We show how genome assembly metrics and gene prediction dramatically improve with our pipeline by assembling two primate genomes solely based on ∼30x coverage of shotgun sequencing data.

Trichoplax genomes reveal profound admixture and suggest stable wild populations without bisexual reproduction

The phylum Placozoa officially consists of only a single described species, Trichoplax adhaerens, although several lineages can be separated by molecular markers, geographical distributions and environmental demands. The placozoan 16S haplotype H2 (Trichoplax sp. H2) is the most robust and cosmopolitan lineage of placozoans found to date. In this study, its genome was found to be distinct but highly related to the Trichoplax adhaerens reference genome, for remarkably unique reasons. The pattern of variation and allele distribution between the two lineages suggests that both originate from a single interbreeding event in the wild, dating back at least several decades ago, and both seem not to have engaged in sexual reproduction since. We conclude that populations of certain placozoan haplotypes remain stable for long periods without bisexual reproduction. Furthermore, allelic variation within and between the two Trichoplax lineages indicates that successful bisexual reproduction between related placozoan lineages might serve to either counter accumulated negative somatic mutations or to cope with changing environmental conditions. On the other hand, enrichment of neutral or beneficial somatic mutations by vegetative reproduction, combined with rare sexual reproduction, could instantaneously boost genetic variation, generating novel ecotypes and eventually species.

Omics approaches to study gene regulatory networks for development in echinoderms

Gene regulatory networks (GRNs) describe the interactions for a developmental process at a given time and space. Historically, perturbation experiments represent one of the key methods for analyzing and reconstructing a GRN, and the GRN governing early development in the sea urchin embryo stands as one of the more deeply dissected so far. As technology progresses, so do the methods used to address different biological questions. Next-generation sequencing (NGS) has become a standard experimental technique for genome and transcriptome sequencing and studies of protein-DNA interactions and DNA accessibility. While several efforts have been made toward the integration of different omics approaches for the study of the regulatory genome in many animals, in a few cases, these are applied with the purpose of reconstructing and experimentally testing developmental GRNs. Here, we review emerging approaches integrating multiple NGS technologies for the prediction and validation of gene interactions within echinoderm GRNs. These approaches can be applied to both 'model' and 'non-model' organisms. Although a number of issues still need to be addressed, advances in NGS applications, such as assay for transposase-accessible chromatin sequencing, combined with the availability of embryos belonging to different species, all separated by various evolutionary distances and accessible to experimental regulatory biology, place echinoderms in an unprecedented position for the reconstruction and evolutionary comparison of developmental GRNs. We conclude that sequencing technologies and integrated omics approaches allow the examination of GRNs on a genome-wide scale only if biological perturbation and cis-regulatory analyses are experimentally accessible, as in the case of echinoderm embryos.

Phylogenetic signal from rearrangements in 18 Anopheles species by joint scaffolding extant and ancestral genomes

Background

Genomes rearrangements carry valuable information for phylogenetic inference or the elucidation of molecular mechanisms of adaptation. However, the detection of genome rearrangements is often hampered by current deficiencies in data and methods: Genomes obtained from short sequence reads have generally very fragmented assemblies, and comparing multiple gene orders generally leads to computationally intractable algorithmic questions.

Results

We present a computational method, ADseq, which, by combining ancestral gene order reconstruction, comparative scaffolding and de novo scaffolding methods, overcomes these two caveats. ADseq provides simultaneously improved assemblies and ancestral genomes, with statistical supports on all local features. Compared to previous comparative methods, it runs in polynomial time, it samples solutions in a probabilistic space, and it can handle a significantly larger gene complement from the considered extant genomes, with complex histories including gene duplications and losses. We use ADseq to provide improved assemblies and a genome history made of duplications, losses, gene translocations, rearrangements, of 18 complete Anopheles genomes, including several important malaria vectors. We also provide additional support for a differentiated mode of evolution of the sex chromosome and of the autosomes in these mosquito genomes.

Conclusions

We demonstrate the method’s ability to improve extant assemblies accurately through a procedure simulating realistic assembly fragmentation. We study a debated issue regarding the phylogeny of the Gambiae complex group of Anopheles genomes in the light of the evolution of chromosomal rearrangements, suggesting that the phylogenetic signal they carry can differ from the phylogenetic signal carried by gene sequences, more prone to introgression.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4466-7) contains supplementary material, which is available to authorized users.

A hybrid reference-guided de novo assembly approach for generating Cyclospora mitochondrion genomes

Cyclospora cayetanensis is a coccidian parasite associated with large and complex foodborne outbreaks worldwide. Linking samples from cyclosporiasis patients during foodborne outbreaks with suspected contaminated food sources, using conventional epidemiological methods, has been a persistent challenge. To address this issue, development of new methods based on potential genomically-derived markers for strain-level identification has been a priority for the food safety research community. The absence of reference genomes to identify nucleotide and structural variants with a high degree of confidence has limited the application of using sequencing data for source tracking during outbreak investigations. In this work, we determined the quality of a high resolution, curated, public mitochondrial genome assembly to be used as a reference genome by applying bioinformatic analyses. Using this reference genome, three new mitochondrial genome assemblies were built starting with metagenomic reads generated by sequencing DNA extracted from oocysts present in stool samples from cyclosporiasis patients. Nucleotide variants were identified in the new and other publicly available genomes in comparison with the mitochondrial reference genome. A consolidated workflow, presented here, to generate new mitochondrion genomes using our reference-guided de novo assembly approach could be useful in facilitating the generation of other mitochondrion sequences, and in their application for subtyping C. cayetanensis strains during foodborne outbreak investigations.

Genomic repeats, misassembly and reannotation: a case study with long-read resequencing of Porphyromonas gingivalis reference strains

BACKGROUND

Without knowledge of their genomic sequences, it is impossible to make functional models of the bacteria that make up human and animal microbiota. Unfortunately, the vast majority of publicly available genomes are only working drafts, an incompleteness that causes numerous problems and constitutes a major obstacle to genotypic and phenotypic interpretation. In this work, we began with an example from the class Bacteroidia in the phylum Bacteroidetes, which is preponderant among human orodigestive microbiota. We successfully identify the genetic loci responsible for assembly breaks and misassemblies and demonstrate the importance and usefulness of long-read sequencing and curated reannotation.

RESULTS

We showed that the fragmentation in Bacteroidia draft genomes assembled from massively parallel sequencing linearly correlates with genomic repeats of the same or greater size than the reads. We also demonstrated that some of these repeats, especially the long ones, correspond to misassembled loci in three reference Porphyromonas gingivalis genomes marked as circularized (thus complete or finished). We prove that even at modest coverage (30X), long-read resequencing together with PCR contiguity verification (rrn operons and an integrative and conjugative element or ICE) can be used to identify and correct the wrongly combined or assembled regions. Finally, although time-consuming and labor-intensive, consistent manual biocuration of three P. gingivalis strains allowed us to compare and correct the existing genomic annotations, resulting in a more accurate interpretation of the genomic differences among these strains.

CONCLUSIONS

In this study, we demonstrate the usefulness and importance of long-read sequencing in verifying published genomes (even when complete) and generating assemblies for new bacterial strains/species with high genomic plasticity. We also show that when combined with biological validation processes and diligent biocurated annotation, this strategy helps reduce the propagation of errors in shared databases, thus limiting false conclusions based on incomplete or misleading information.

Draft Genome Sequence of USA100 Methicillin-Resistant Staphylococcus aureus Strain 209

USA100 strains are significant contributors to the overall burden of health care-associated methicillin-resistant Staphylococcus aureus (MRSA) infections. Strain 209 is a representative MRSA isolate that serves as a model organism for agr type II studies and USA100 virulence assessments. We present a draft genome sequence of this strain.

Draft Genome Sequence of the Mucoid Pseudomonas aeruginosa Clinical Isolate PA34

Pseudomonas aeruginosa is a serious threat to patients suffering from cystic fibrosis. These organisms are exposed to a unique set of selective pressures within the lung. Here, we report the draft genome sequence of a mucoid P. aeruginosa clinical isolate obtained from a cystic fibrosis patient colonized with P. aeruginosa.

Reference-guided de novo assembly approach improves genome reconstruction for related species

Background

The development of next-generation sequencing has made it possible to sequence whole genomes at a relatively low cost. However, de novo genome assemblies remain challenging due to short read length, missing data, repetitive regions, polymorphisms and sequencing errors. As more and more genomes are sequenced, reference-guided assembly approaches can be used to assist the assembly process. However, previous methods mostly focused on the assembly of other genotypes within the same species. We adapted and extended a reference-guided de novo assembly approach, which enables the usage of a related reference sequence to guide the genome assembly. In order to compare and evaluate de novo and our reference-guided de novo assembly approaches, we used a simulated data set of a repetitive and heterozygotic plant genome.

Results

The extended reference-guided de novo assembly approach almost always outperforms the corresponding de novo assembly program even when a reference of a different species is used. Similar improvements can be observed in high and low coverage situations. In addition, we show that a single evaluation metric, like the widely used N50 length, is not enough to properly rate assemblies as it not always points to the best assembly evaluated with other criteria. Therefore, we used the summed z-scores of 36 different statistics to evaluate the assemblies.

Conclusions

The combination of reference mapping and de novo assembly provides a powerful tool to improve genome reconstruction by integrating information of a related genome. Our extension of the reference-guided de novo assembly approach enables the application of this strategy not only within but also between related species. Finally, the evaluation of genome assemblies is often not straight forward, as the truth is not known. Thus one should always use a combination of evaluation metrics, which not only try to assess the continuity but also the accuracy of an assembly.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1911-6) contains supplementary material, which is available to authorized users.

Pathogenicity Locus, Core Genome, and Accessory Gene Contributions to Clostridium difficile Virulence

Clostridium difficile is a spore-forming anaerobic bacterium that causes colitis in patients with disrupted colonic microbiota. While some individuals are asymptomatic C. difficile carriers, symptomatic disease ranges from mild diarrhea to potentially lethal toxic megacolon. The wide disease spectrum has been attributed to the infected host’s age, underlying diseases, immune status, and microbiome composition. However, strain-specific differences in C. difficile virulence have also been implicated in determining colitis severity. Because patients infected with C. difficile are unique in terms of medical history, microbiome composition, and immune competence, determining the relative contribution of C. difficile virulence to disease severity has been challenging, and conclusions regarding the virulence of specific strains have been inconsistent. To address this, we used a mouse model to test 33 clinical C. difficile strains isolated from patients with disease severities ranging from asymptomatic carriage to severe colitis, and we determined their relative in vivo virulence in genetically identical, antibiotic-pretreated mice. We found that murine infections with C. difficile clade 2 strains (including multilocus sequence type 1/ribotype 027) were associated with higher lethality and that C. difficile strains associated with greater human disease severity caused more severe disease in mice. While toxin production was not strongly correlated with in vivo colonic pathology, the ability of C. difficile strains to grow in the presence of secondary bile acids was associated with greater disease severity. Whole-genome sequencing and identification of core and accessory genes identified a subset of accessory genes that distinguish high-virulence from lower-virulence C. difficile strains.

Clostridium difficile is an important cause of hospital-associated intestinal infections, and recent years have seen an increase in the number and severity of cases in the United States. A patient’s antibiotic history, immune status, and medical comorbidities determine, in part, the severity of C. difficile infection. The relative virulence of different clinical C. difficile strains, although postulated to determine disease severity in patients, has been more difficult to consistently associate with mild versus severe colitis. We tested 33 distinct clinical C. difficile isolates for their ability to cause disease in genetically identical mice and found that C. difficile strains belonging to clade 2 were associated with higher mortality. Differences in survival were not attributed to differences in toxin production but likely resulted from the distinct gene content in the various clinical isolates.

GAAP: Genome-organization-framework-Assisted Assembly Pipeline for prokaryotic genomes

BACKGROUND

Next-generation sequencing (NGS) technologies have greatly promoted the genomic study of prokaryotes. However, highly fragmented assemblies due to short reads from NGS are still a limiting factor in gaining insights into the genome biology. Reference-assisted tools are promising in genome assembly, but tend to result in false assembly when the assigned reference has extensive rearrangements.

RESULTS

Herein, we present GAAP, a genome assembly pipeline for scaffolding based on core-gene-defined Genome Organizational Framework (cGOF) described in our previous study. Instead of assigning references, we use the multiple-reference-derived cGOFs as indexes to assist in order and orientation of the scaffolds and build a skeleton structure, and then use read pairs to extend scaffolds, called local scaffolding, and distinguish between true and chimeric adjacencies in the scaffolds. In our performance tests using both empirical and simulated data of 15 genomes in six species with diverse genome size, complexity, and all three categories of cGOFs, GAAP outcompetes or achieves comparable results when compared to three other reference-assisted programs, AlignGraph, Ragout and MeDuSa.

CONCLUSIONS

GAAP uses both cGOF and pair-end reads to create assemblies in genomic scale, and performs better than the currently available reference-assisted assembly tools as it recovers more assemblies and makes fewer false locations, especially for species with extensive rearranged genomes. Our method is a promising solution for reconstruction of genome sequence from short reads of NGS.

A Bacillus anthracis Genome Sequence from the Sverdlovsk 1979 Autopsy Specimens

Anthrax is a zoonotic disease that occurs naturally in wild and domestic animals but has been used by both state-sponsored programs and terrorists as a biological weapon. A Soviet industrial production facility in Sverdlovsk, USSR, proved deficient in 1979 when a plume of spores was accidentally released and resulted in one of the largest known human anthrax outbreaks. In order to understand this outbreak and others, we generated a Bacillus anthracis population genetic database based upon whole-genome analysis to identify all single-nucleotide polymorphisms (SNPs) across a reference genome. Phylogenetic analysis has defined three major clades (A, B, and C), B and C being relatively rare compared to A. The A clade has numerous subclades, including a major polytomy named the trans-Eurasian (TEA) group. The TEA radiation is a dominant evolutionary feature of B. anthracis, with many contemporary populations having resulted from a large spatial dispersal of spores from a single source. Two autopsy specimens from the Sverdlovsk outbreak were deep sequenced to produce draft B. anthracis genomes. This allowed the phylogenetic placement of the Sverdlovsk strain into a clade with two Asian live vaccine strains, including the Russian Tsiankovskii strain. The genome was examined for evidence of drug resistance manipulation or other genetic engineering, but none was found. The Soviet Sverdlovsk strain genome is consistent with a wild-type strain from Russia that had no evidence of genetic manipulation during its industrial production. This work provides insights into the world’s largest biological weapons program and provides an extensive B. anthracis phylogenetic reference.

The 1979 Russian anthrax outbreak resulted from an industrial accident at the Soviet anthrax spore production facility in the city of Sverdlovsk. Deep genomic sequencing of two autopsy specimens generated a draft genome and phylogenetic placement of the Soviet Sverdlovsk anthrax strain. While it is known that Soviet scientists had genetically manipulated Bacillus anthracis with the potential to evade vaccine prophylaxis and antibiotic therapeutics, there was no genomic evidence of this from the Sverdlovsk production strain genome. The whole-genome SNP genotype of the Sverdlovsk strain was used to precisely identify it and its close relatives in the context of an extensive global B. anthracis strain collection. This genomic identity can now be used for forensic tracking of this weapons material on a global scale and for future anthrax investigations.

Chromosomer: a reference-based genome arrangement tool for producing draft chromosome sequences

BACKGROUND

As the number of sequenced genomes rapidly increases, chromosome assembly is becoming an even more crucial step of any genome study. Since de novo chromosome assemblies are confounded by repeat-mediated artifacts, reference-assisted assemblies that use comparative inference have become widely used, prompting the development of several reference-assisted assembly programs for prokaryotic and eukaryotic genomes.

FINDINGS

We developed Chromosomer - a reference-based genome arrangement tool, which rapidly builds chromosomes from genome contigs or scaffolds using their alignments to a reference genome of a closely related species. Chromosomer does not require mate-pair libraries and it offers a number of auxiliary tools that implement common operations accompanying the genome assembly process.

CONCLUSIONS

Despite implementing a straightforward alignment-based approach, Chromosomer is a useful tool for genomic analysis of species without chromosome maps. Putative chromosome assemblies by Chromosomer can be used in comparative genomic analysis, genomic variation assessment, potential linkage group inference and other kinds of analysis involving contig or scaffold mapping to a high-quality assembly.

Chapter 6. Comparative genomic approaches to investigate molecular traits specific to social insects

Ageing is a feature of nearly all known organisms and, by its connection to survival, appears to trade off with fecundity. However, in some organisms such as in queens of social insects, this negative relation appears reversed and individuals live long and reproduce much. Since new experimental techniques, transcriptomes and genomes of many social insects have recently become available, a comparison of these data in a phylogenetic framework becomes feasible. This allows the study of general trends, species specific oddities and evolutionary dynamics of the molecular properties and changes which underlie ageing, fecundity and the reversal of this negative association. In the framework of social insect evolution, we review the most important recent insights, computational methods, their applications and data resources which are available.

In silico serotyping of E. coli from short read data identifies limited novel O-loci but extensive diversity of O:H serotype combinations within and between pathogenic lineages

The lipopolysaccharide (O) and flagellar (H) surface antigens of Escherichia coli are targets for serotyping that have traditionally been used to identify pathogenic lineages. These surface antigens are important for the survival of E. coli within mammalian hosts. However, traditional serotyping has several limitations, and public health reference laboratories are increasingly moving towards whole genome sequencing (WGS) to characterize bacterial isolates. Here we present a method to rapidly and accurately serotype E. coli isolates from raw, short read WGS data. Our approach bypasses the need for de novo genome assembly by directly screening WGS reads against a curated database of alleles linked to known and novel E. coli O-groups and H-types (the EcOH database) using the software package srst2. We validated the approach by comparing in silico results for 197 enteropathogenic E. coli isolates with those obtained by serological phenotyping in an independent laboratory. We then demonstrated the utility of our method to characterize isolates in public health and clinical settings, and to explore the genetic diversity of >1500 E. coli genomes from multiple sources. Importantly, we showed that transfer of O- and H-antigen loci between E. coli chromosomal backbones is common, with little evidence of constraints by host or pathotype, suggesting that E. coli ‘strain space’ may be virtually unlimited, even within specific pathotypes. Our findings show that serotyping is most useful when used in combination with strain genotyping to characterize microevolution events within an inferred population structure.

Transcriptome analysis reveals differentially expressed genes associated with germ cell and gonad development in the Southern bluefin tuna (Thunnus maccoyii)

BACKGROUND

Controlling and managing the breeding of bluefin tuna (Thunnus spp.) in captivity is an imperative step towards obtaining a sustainable supply of these fish in aquaculture production systems. Germ cell transplantation (GCT) is an innovative technology for the production of inter-species surrogates, by transplanting undifferentiated germ cells derived from a donor species into larvae of a host species. The transplanted surrogates will then grow and mature to produce donor-derived seed, thus providing a simpler alternative to maintaining large-bodied broodstock such as the bluefin tuna. Implementation of GCT for new species requires the development of molecular tools to follow the fate of the transplanted germ cells. These tools are based on key reproductive and germ cell-specific genes. RNA-Sequencing (RNA-Seq) provides a rapid, cost-effective method for high throughput gene identification in non-model species. This study utilized RNA-Seq to identify key genes expressed in the gonads of Southern bluefin tuna (Thunnus maccoyii, SBT) and their specific expression patterns in male and female gonad cells.

RESULTS

Key genes involved in the reproductive molecular pathway and specifically, germ cell development in gonads, were identified using analysis of RNA-Seq transcriptomes of male and female SBT gonad cells. Expression profiles of transcripts from ovary and testis cells were compared, as well as testis germ cell-enriched fraction prepared with Percoll gradient, as used in GCT studies. Ovary cells demonstrated over-expression of genes related to stem cell maintenance, while in testis cells, transcripts encoding for reproduction-associated receptors, sex steroids and hormone synthesis and signaling genes were over-expressed. Within the testis cells, the Percoll-enriched fraction showed over-expression of genes that are related to post-meiosis germ cell populations.

CONCLUSIONS

Gonad development and germ cell related genes were identified from SBT gonads and their expression patterns in ovary and testis cells were determined. These expression patterns correlate with the reproductive developmental stage of the sampled fish. The majority of the genes described in this study were sequenced for the first time in T. maccoyii. The wealth of SBT gonadal and germ cell-related gene sequences made publicly available by this study provides an extensive resource for further GCT and reproductive molecular biology studies of this commercially valuable fish.

RECORD: Reference-Assisted Genome Assembly for Closely Related Genomes

Background. Next-generation sequencing technologies are now producing multiple times the genome size in total reads from a single experiment. This is enough information to reconstruct at least some of the differences between the individual genome studied in the experiment and the reference genome of the species. However, in most typical protocols, this information is disregarded and the reference genome is used. Results. We provide a new approach that allows researchers to reconstruct genomes very closely related to the reference genome (e.g., mutants of the same species) directly from the reads used in the experiment. Our approach applies de novo assembly software to experimental reads and so-called pseudoreads and uses the resulting contigs to generate a modified reference sequence. In this way, it can very quickly, and at no additional sequencing cost, generate new, modified reference sequence that is closer to the actual sequenced genome and has a full coverage. In this paper, we describe our approach and test its implementation called RECORD. We evaluate RECORD on both simulated and real data. We made our software publicly available on sourceforge. Conclusion. Our tests show that on closely related sequences RECORD outperforms more general assisted-assembly software.

Draft Genome Sequence of a Clinical Isolate of Serratia marcescens, Strain AH0650_Sm1

Serratia marcescens strain AH0650_Sm1 is a clinical multidrug-resistant isolate from Australia. Here, we report its annotated draft genome comprising 20 contigs. We identified chromosomal antimicrobial resistance genes including a tet(41) variant, an aac(6′)-Ic variant, ampC, a metallo-beta-lactamase, and several putative multidrug efflux pumps, as well as a novel prophage.