A post-assembly genome-improvement toolkit (PAGIT) to obtain annotated genomes from contigs

Detecting Assembly Errors With Klumpy: Building Confidence in Your Daily Genomic Analysis

In the realm of genome assembly, even minor errors can send researchers down to rabbit holes of unintended misinterpretation. Enter Klumpy-a tool designed to help detecting these elusive mistakes before they cause significant problems. By providing detailed, region-specific assessments and an intuitive visualisation platform, Klumpy (Madrigal, et al. 2024) empowers researchers to pinpoint and resolve potential issues with precision, paving the way for more reliable downstream analyses and discoveries.

Exploring Pan-Genomes: An Overview of Resources and Tools for Unraveling Structure, Function, and Evolution of Crop Genes and Genomes

The availability of multiple sequenced genomes from a single species made it possible to explore intra- and inter-specific genomic comparisons at higher resolution and build clade-specific pan-genomes of several crops. The pan-genomes of crops constructed from various cultivars, accessions, landraces, and wild ancestral species represent a compendium of genes and structural variations and allow researchers to search for the novel genes and alleles that were inadvertently lost in domesticated crops during the historical process of crop domestication or in the process of extensive plant breeding. Fortunately, many valuable genes and alleles associated with desirable traits like disease resistance, abiotic stress tolerance, plant architecture, and nutrition qualities exist in landraces, ancestral species, and crop wild relatives. The novel genes from the wild ancestors and landraces can be introduced back to high-yielding varieties of modern crops by implementing classical plant breeding, genomic selection, and transgenic/gene editing approaches. Thus, pan-genomic represents a great leap in plant research and offers new avenues for targeted breeding to mitigate the impact of global climate change. Here, we summarize the tools used for pan-genome assembly and annotations, web-portals hosting plant pan-genomes, etc. Furthermore, we highlight a few discoveries made in crops using the pan-genomic approach and future potential of this emerging field of study.

From contigs towards chromosomes: automatic improvement of long read assemblies (ILRA)

Recent advances in long read technologies not only enable large consortia to aim to sequence all eukaryotes on Earth, but they also allow individual laboratories to sequence their species of interest with relatively low investment. Long read technologies embody the promise of overcoming scaffolding problems associated with repeats and low complexity sequences, but the number of contigs often far exceeds the number of chromosomes and they may contain many insertion and deletion errors around homopolymer tracts. To overcome these issues, we have implemented the ILRA pipeline to correct long read-based assemblies. Contigs are first reordered, renamed, merged, circularized, or filtered if erroneous or contaminated. Illumina short reads are used subsequently to correct homopolymer errors. We successfully tested our approach by improving the genome sequences of Homo sapiens, Trypanosoma brucei, and Leptosphaeria spp., and by generating four novel Plasmodium falciparum assemblies from field samples. We found that correcting homopolymer tracts reduced the number of genes incorrectly annotated as pseudogenes, but an iterative approach seems to be required to correct more sequencing errors. In summary, we describe and benchmark the performance of our new tool, which improved the quality of novel long read assemblies up to 1 Gbp. The pipeline is available at GitHub: https://github.com/ThomasDOtto/ILRA.

Repeat infections with chlamydia in women may be more transcriptionally active with lower responses from some immune genes

Chlamydia trachomatis, the most common bacterial sexually transmitted infection worldwide, is responsible for considerable health burden due to its significant sequelae. There are growing concerns about chlamydial treatment and management due to widely documented increasing burden of repeat infections. In the current study, a cohort study design of 305 women with urogenital chlamydial infections demonstrated that 11.8% of women experienced repeat infections after treatment with azithromycin. The chlamydial DNA load measured by quantitative PCR was higher in women who experienced a repeat infection (p = 0.0097) and repeat infection was associated with sexual contact. There was no genomic or phenotypic evidence of azithromycin resistance within the chlamydial isolates. During repeat infection, or repeat positive tests during follow up, vaginal chlamydial gene expression (ompA, euo, omcB, htrA, trpAB) was markedly higher compared to baseline, and two of the selected immune genes analyzed had significantly lower expression at the time of repeat infection. Overall, there are two implications of these results. The results could be generalized to all recent infections, or repeat positive events, and indicate that chlamydial infections are have higher transcriptional activity of select genes early in the infection in women. Alternatively, after azithromycin treatment, repeat infections of Chlamydia may be more transcriptionally active at certain genes, and there may be post-treatment immunological alterations that interplay into repeat exposures establishing an active infection. The potential that recent infections may involve a higher level of activity from the organism may have implications for management by more regular testing of the most at risk women to reduce the risk of sequelae.

Low carbon-to-nitrogen ratio digestate from high-rate anaerobic baffled reactor facilitates heterotrophic/autotrophic nitrifiers involved in nitrogen removal

In this study, baffled anaerobic-aerobic reactors (AOBRs) with modified basalt fiber (MBF) carriers and felt were used to treat domestic wastewater (DWW). The influent was first treated in anaerobic compartments, with the NH₄⁺-N containing digestate refluxed into aerobic compartment for nitrification. The nitrified liquid was channeled to the anaerobic compartments for further denitrification. Under optimal conditions, AOBR with MBF carriers could remove 91% chemical oxygen demand (COD) and 81% total nitrogen (TN), with biomass production increased by 7.6%, 4.5% and 8.7% in three successive anaerobic compartments compared to the control. Biological viability analysis showed that live cells outnumbered dead cells in bio-nests. Metagenomics analysis showed that multiple metabolic pathways accounted for nitrogen conversion in anaerobic and aerobic compartments. More importantly, low COD/TN ratio digestate facilitated heterotrophic nitrification-aerobic denitrification (HN-AD) species growth in aerobic compartment. This study provides a promising strategy to source treatment of DWW from urban communities.

Long-read assembly and comparative evidence-based reanalysis of Cryptosporidium genome sequences reveal expanded transporter repertoire and duplication of entire chromosome ends including subtelomeric regions

Cryptosporidiosis is a leading cause of waterborne diarrheal disease globally and an important contributor to mortality in infants and the immunosuppressed. Despite its importance, the Cryptosporidium community has only had access to a good, but incomplete, Cryptosporidium parvum IOWA reference genome sequence. Incomplete reference sequences hamper annotation, experimental design, and interpretation. We have generated a new C. parvum IOWA genome assembly supported by Pacific Biosciences (PacBio) and Oxford Nanopore long-read technologies and a new comparative and consistent genome annotation for three closely related species: C. parvum, Cryptosporidium hominis, and Cryptosporidium tyzzeri We made 1926 C. parvum annotation updates based on experimental evidence. They include new transporters, ncRNAs, introns, and altered gene structures. The new assembly and annotation revealed a complete Dnmt2 methylase ortholog. Comparative annotation between C. parvum, C. hominis, and C. tyzzeri revealed that most "missing" orthologs are found, suggesting that the biological differences between the species must result from gene copy number variation, differences in gene regulation, and single-nucleotide variants (SNVs). Using the new assembly and annotation as reference, 190 genes are identified as evolving under positive selection, including many not detected previously. The new C. parvum IOWA reference genome assembly is larger, gap free, and lacks ambiguous bases. This chromosomal assembly recovers all 16 chromosome ends, 13 of which are contiguously assembled. The three remaining chromosome ends are provisionally placed. These ends represent duplication of entire chromosome ends including subtelomeric regions revealing a new level of genome plasticity that will both inform and impact future research.

Microbial methylglyoxal metabolism contributes towards growth promotion and stress tolerance in plants

Plant growth promotion by microbes is a cumulative phenomenon involving multiple traits, many of which are not explored yet. Hence, to unravel microbial mechanisms underlying growth promotion, we have analysed the genomes of two potential growth-promoting microbes, viz., Pseudomonas sp. CK-NBRI-02 (P2) and Bacillus marisflavi CK-NBRI-03 (P3) for the presence of plant-beneficial traits. Besides known traits, we found that microbes differ in their ability to metabolize methylglyoxal (MG), a ubiquitous cytotoxin regarded as general consequence of stress in plants. P2 exhibited greater tolerance to MG and possessed better ability to sustain plant growth under dicarbonyl stress. However, under salinity, only P3 showed a dose-dependent induction in MG detoxification activity in accordance with concomitant increase in MG levels, contributing to enhanced salt tolerance. Furthermore, salt-stressed transcriptomes of both the strains showed differences with respect to MG, ion and osmolyte homeostasis, with P3 being more responsive to stress. Importantly, application of either strain altered MG levels and subsequently MG detoxification machinery in Arabidopsis, probably to strengthen plant defence response and growth. We therefore, suggest a crucial role of microbial MG resistance in plant growth promotion and that it should be considered as a beneficial trait while screening microbes for stress mitigation in plants.

Choice of assemblers has a critical impact on de novo assembly of SARS-CoV-2 genome and characterizing variants

Background

Coronavirus Disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has become a global pandemic following its initial emergence in China. SARS-CoV-2 has a positive-sense single-stranded RNA virus genome of around 30Kb. Using next-generation sequencing technologies, a large number of SARS-CoV-2 genomes are being sequenced at an unprecedented rate and being deposited in public repositories. For the de novo assembly of the SARS-CoV-2 genomes, a myriad of assemblers is being used, although their impact on the assembly quality has not been characterized for this virus. In this study, we aim to understand the variabilities on assembly qualities due to the choice of the assemblers.

Results

We performed 6648 de novo assemblies of 416 SARS-CoV-2 samples using eight different assemblers with different k-mer lengths. We used Illumina paired-end sequencing reads and compared the assembly quality of those assemblers. We showed that the choice of assembler plays a significant role in reconstructing the SARS-CoV-2 genome. Two metagenomic assemblers, e.g. MEGAHIT and metaSPAdes, performed better compared with others in most of the assembly quality metrics including, recovery of a larger fraction of the genome, constructing larger contigs and higher N50, NA50 values, etc. We showed that at least 09% (259/2873) of the variants present in the assemblies between MEGAHIT and metaSPAdes are unique to one of the assembly methods.

Conclusion

Our analyses indicate the critical role of assembly methods for assembling SARS-CoV-2 genome using short reads and their impact on variant characterization. This study could help guide future studies to determine the best-suited assembler for the de novo assembly of virus genomes.

ESREEM: Efficient Short Reads Error Estimation Computational Model for Next-generation Genome Sequencing

Aims:

To assess the error profile in NGS data, generated from high throughput sequencing machines.

Background:

Short-read sequencing data from Next Generation Sequencing (NGS) are currently being generated by a number of research projects. Depicting the errors produced by NGS platforms and expressing accurate genetic variation from reads are two inter-dependent phases. It has high significance in various analyses, such as genome sequence assembly, SNPs calling, evolutionary studies, and haplotype inference. The systematic and random errors show incidence profile for each of the sequencing platforms i.e. Illumina sequencing, Pacific Biosciences, 454 pyrosequencing, Complete Genomics DNA nanoball sequencing, Ion Torrent sequencing, and Oxford Nanopore sequencing. Advances in NGS deliver galactic data with the addition of errors. Some ratio of these errors may emulate genuine true biological signals i.e., mutation, and may subsequently negate the results. Various independent applications have been proposed to correct the sequencing errors. Systematic analysis of these algorithms shows that state-of-the-art models are missing.

Objective:

In this paper, an effcient error estimation computational model called ESREEM is proposed to assess the error rates in NGS data.

Methods:

The proposed model prospects the analysis that there exists a true linear regression association between the number of reads containing errors and the number of reads sequenced. The model is based on a probabilistic error model integrated with the Hidden Markov Model (HMM).

Result:

The proposed model is evaluated on several benchmark datasets and the results obtained are compared with state-of-the-art algorithms.

Conclusions:

Experimental results analyses show that the proposed model efficiently estimates errors and runs in less time as compared to others.

A novel terpene synthase controls differences in anti-aphrodisiac pheromone production between closely related Heliconius butterflies

Plants and insects often use the same compounds for chemical communication, but not much is known about the genetics of convergent evolution of chemical signals. The terpene (E)-β-ocimene is a common component of floral scent and is also used by the butterfly Heliconius melpomene as an anti-aphrodisiac pheromone. While the biosynthesis of terpenes has been described in plants and microorganisms, few terpene synthases (TPSs) have been identified in insects. Here, we study the recent divergence of 2 species, H. melpomene and Heliconius cydno, which differ in the presence of (E)-β-ocimene; combining linkage mapping, gene expression, and functional analyses, we identify 2 novel TPSs. Furthermore, we demonstrate that one, HmelOS, is able to synthesise (E)-β-ocimene in vitro. We find no evidence for TPS activity in HcydOS (HmelOS ortholog of H. cydno), suggesting that the loss of (E)-β-ocimene in this species is the result of coding, not regulatory, differences. The TPS enzymes we discovered are unrelated to previously described plant and insect TPSs, demonstrating that chemical convergence has independent evolutionary origins.

Plants and insects often use the same compounds for chemical communication, but little is known about the convergent evolution of such chemical signals. This study identifies a novel terpene synthase involved in production of an anti-aphrodisiac pheromone by the butterfly Heliconius melpomene. This enzyme is unrelated to other insect terpene synthases, providing evidence that the ability to synthesise terpenes has arisen multiple times independently within the insects.

Development of a PCR-based, genetic marker resource for the tomato-like nightshade relative, Solanum lycopersicoides using whole genome sequence analysis

Solanum lycopersicoides is a wild nightshade relative of tomato with known resistance to a wide range of pests and pathogens, as well as tolerance to cold, drought and salt stress. To effectively utilize S. lycopersicoides as a genetic resource in breeding for tomato improvement, the underlying basis of observable traits in the species needs to be understood. Molecular markers are important tools that can unlock the genetic underpinnings of phenotypic variation in wild crop relatives. Unfortunately, DNA markers that are specific to S. lycopersicoides are limited in number, distribution and polymorphism rate. In this study, we developed a suite of S. lycopersicoides-specific SSR and indel markers by sequencing, building and analyzing a draft assembly of the wild nightshade genome. Mapping of a total of 1.45 Gb of S. lycopersicoides contigs against the tomato reference genome assembled a moderate number of contiguous reads into longer scaffolds. Interrogation of the obtained draft yielded SSR information for more than 55,000 loci in S. lycopersicoides for which more than 35,000 primers pairs were designed. Additionally, indel markers were developed based on sequence alignments between S. lycopersicoides and tomato. Synthesis and experimental validation of 345 primer sets resulted in the amplification of single and multilocus targets in S. lycopersicoides and polymorphic loci between S. lycopersicoides and tomato. Cross-species amplification of the 345 markers in tomato, eggplant, silverleaf nightshade and pepper resulted in varying degrees of transferability that ranged from 55 to 83%. The markers reported in this study significantly expands the genetic marker resource for S. lycopersicoides, as well as for related Solanum spp. for applications in genetics and breeding studies.

Epstein-Barr Virus Genomes Reveal Population Structure and Type 1 Association with Endemic Burkitt Lymphoma

Endemic Burkitt lymphoma (eBL), the most prevalent pediatric cancer in sub-Saharan Africa, is distinguished by its inclusion of Epstein-Barr virus (EBV). In order to better understand the impact of EBV variation in eBL tumorigenesis, we improved viral DNA enrichment methods and generated a total of 98 new EBV genomes from both eBL cases (n = 58) and healthy controls (n = 40) residing in the same geographic region in Kenya. Using our unbiased methods, we found that EBV type 1 was significantly more prevalent in eBL patients (74.5%) than in healthy children (47.5%) (odds ratio = 3.24, 95% confidence interval = 1.36 to 7.71, P = 0.007), as opposed to similar proportions in both groups. Controlling for EBV type, we also performed a genome-wide association study identifying six nonsynonymous variants in the genes EBNA1, EBNA2, BcLF1, and BARF1 that were enriched in eBL patients. In addition, viruses isolated from plasma of eBL patients were identical to their tumor counterparts consistent with circulating viral DNA originating from the tumor. We also detected three intertypic recombinants carrying type 1 EBNA2 and type 2 EBNA3 regions, as well as one novel genome with a 20-kb deletion, resulting in the loss of multiple lytic and virion genes. Comparing EBV types, viral genes displayed differential variation rates as type 1 appeared to be more divergent, while type 2 demonstrated novel substructures. Overall, our findings highlight the complexities of the EBV population structure and provide new insight into viral variation, potentially deepening our understanding of eBL oncogenesis.IMPORTANCE Improved viral enrichment methods conclusively demonstrate EBV type 1 to be more prevalent in eBL patients than in geographically matched healthy controls, which previously underrepresented the prevalence of EBV type 2. Genome-wide association analysis between cases and controls identifies six eBL-associated nonsynonymous variants in EBNA1, EBNA2, BcLF1, and BARF1 genes. Analysis of population structure reveals that EBV type 2 exists as two genomic subgroups and was more commonly found in female than in male eBL patients.

Draft Genome Sequence of Bacillus marisflavi CK-NBRI-03, Isolated from Agricultural Soil

Here, we report the 4.34-Mb draft genome assembly of Bacillus marisflavi CK-NBRI-03 (or P3), a Gram-positive bacterium, with an average G+C content of 48.66%. P3 was isolated from agricultural soil from the Badaun (midwestern plain zone) region of Uttar Pradesh, India.

Here, we report the 4.34-Mb draft genome assembly of Bacillus marisflavi CK-NBRI-03 (or P3), a Gram-positive bacterium, with an average G+C content of 48.66%. P3 was isolated from agricultural soil from the Badaun (midwestern plain zone) region of Uttar Pradesh, India.

Draft Genome Sequences of Two Textile Azo Dye-Degrading Shewanella sp. Strains Isolated from a Textile Effluent in Peru

Here, we report the annotated genome sequences of two Shewanella sp. strains isolated from textile industry wastewater effluent in Peru. Potential genes for encoding enzymes that enable the strain to decolorize and degrade textile azo dyes were detected in both genomes.

Here, we report the annotated genome sequences of two Shewanella sp. strains isolated from textile industry wastewater effluent in Peru. Potential genes for encoding enzymes that enable the strain to decolorize and degrade textile azo dyes were detected in both genomes.

Complement Genome Annotation Lift Over Using a Weighted Sequence Alignment Strategy

With the broad application of high-throughput sequencing, more whole-genome resequencing data and de novo assemblies of natural populations are becoming available. For a particular species, in general, only the reference genome is well established and annotated. Computational tools based on sequence alignment have been developed to investigate the gene models of individuals belonging to the same or closely related species. During this process, inconsistent alignment often obscures genome annotation lift over and leads to improper functional impact prediction for a genomic variant, especially in plant species. Here, we proposed the zebraic striped dynamic programming algorithm, which provides different weights to genetic features to refine genome annotation lift over. Testing of our zebraic striped dynamic programming algorithm on both plant and animal genomic data showed complementation to standard sequence approach for highly diverse individuals. Using the lift over genome annotation as anchors, a base-pair resolution genome-wide sequence alignment and variant calling pipeline for de novo assembly has been implemented in the GEAN software. GEAN could be used to compare haplotype diversity, refine the genetic variant functional annotation, annotate de novo assembly genome sequence, detect homologous syntenic blocks, improve the quantification of gene expression levels using RNA-seq data, and unify genomic variants for population genetic analysis. We expect that GEAN will be a standard tool for the coming of age of de novo assembly population genetics.

Draft Genome Sequence of a Potential Plant Growth-Promoting Rhizobacterium, Pseudomonas sp. Strain CK-NBRI-02

Pseudomonas sp. strain CK-NBRI-02 is a potential plant growth-promoting Gram-negative rhizobacterium isolated from the rhizosphere of maize plants growing in fields in Srinagar, Jammu, and Kashmir, India. Here, we report a 5.25-Mb draft assembly of the genome sequence of Pseudomonas sp. strain CK-NBRI-02 with an average G+C content of 62.47%.

Pseudomonas sp. strain CK-NBRI-02 is a potential plant growth-promoting Gram-negative rhizobacterium isolated from the rhizosphere of maize plants growing in fields in Srinagar, Jammu, and Kashmir, India. Here, we report a 5.25-Mb draft assembly of the genome sequence of Pseudomonas sp. strain CK-NBRI-02 with an average G+C content of 62.47%.

BAUM: improving genome assembly by adaptive unique mapping and local overlap-layout-consensus approach

Motivation

It is highly desirable to assemble genomes of high continuity and consistency at low cost. The current bottleneck of draft genome continuity using the second generation sequencing (SGS) reads is primarily caused by uncertainty among repetitive sequences. Even though the single-molecule real-time sequencing technology is very promising to overcome the uncertainty issue, its relatively high cost and error rate add burden on budget or computation. Many long-read assemblers take the overlap-layout-consensus (OLC) paradigm, which is less sensitive to sequencing errors, heterozygosity and variability of coverage. However, current assemblers of SGS data do not sufficiently take advantage of the OLC approach.

Results

Aiming at minimizing uncertainty, the proposed method BAUM, breaks the whole genome into regions by adaptive unique mapping; then the local OLC is used to assemble each region in parallel. BAUM can (i) perform reference-assisted assembly based on the genome of a close species (ii) or improve the results of existing assemblies that are obtained based on short or long sequencing reads. The tests on two eukaryote genomes, a wild rice Oryza longistaminata and a parrot Melopsittacus undulatus, show that BAUM achieved substantial improvement on genome size and continuity. Besides, BAUM reconstructed a considerable amount of repetitive regions that failed to be assembled by existing short read assemblers. We also propose statistical approaches to control the uncertainty in different steps of BAUM.

Availability and implementation

http://www.zhanyuwang.xin/wordpress/index.php/2017/07/21/baum.

Supplementary information

Supplementary data are available at Bioinformatics online.

Whole-Genome Sequence of Brucella melitensis Strain B7, Isolated from a Blood Sample of a Brucellosis Patient from Hulunbuir, Inner Mongolia, China

We report here the draft genome sequence of Brucella melitensis strain B7, isolated from a blood sample of a brucellosis patient from Hulunbuir, Inner Mongolia, China. This report describes the genome sequence of a human B. melitensis isolate, which is endemic in this area of China.

We report here the draft genome sequence of Brucella melitensis strain B7, isolated from a blood sample of a brucellosis patient from Hulunbuir, Inner Mongolia, China. This report describes the genome sequence of a human B. melitensis isolate, which is endemic in this area of China.

Progression of the canonical reference malaria parasite genome from 2002-2019

Here we describe the ways in which the sequence and annotation of the Plasmodium falciparum reference genome has changed since its publication in 2002. As the malaria species responsible for the most deaths worldwide, the richness of annotation and accuracy of the sequence are important resources for the P. falciparum research community as well as the basis for interpreting the genomes of subsequently sequenced species. At the time of publication in 2002 over 60% of predicted genes had unknown functions. As of March 2019, this number has been significantly decreased to 33%. The reduction is due to the inclusion of genes that were subsequently characterised experimentally and genes with significant similarity to others with known functions. In addition, the structural annotation of genes has been significantly refined; 27% of gene structures have been changed since 2002, comprising changes in exon-intron boundaries, addition or deletion of exons and the addition or deletion of genes. The sequence has also undergone significant improvements. In addition to the correction of a large number of single-base and insertion or deletion errors, a major miss-assembly between the subtelomeres of chromosome 7 and 8 has been corrected. As the number of sequenced isolates continues to grow rapidly, a single reference genome will not be an adequate basis for interpreting intra-species sequence diversity. We therefore describe in this publication a population reference genome of P. falciparum, called Pfref1. This reference will enable the community to map to regions that are not present in the current assembly. P. falciparum 3D7 will continue to be maintained, with ongoing curation ensuring continual improvements in annotation quality.

Progression of the canonical reference malaria parasite genome from 2002-2019

Here we describe the ways in which the sequence and annotation of the Plasmodium falciparum reference genome has changed since its publication in 2002. As the malaria species responsible for the most deaths worldwide, the richness of annotation and accuracy of the sequence are important resources for the P. falciparum research community as well as the basis for interpreting the genomes of subsequently sequenced species. At the time of publication in 2002 over 60% of predicted genes had unknown functions. As of March 2019, this number has been significantly decreased to 33%. The reduction is due to the inclusion of genes that were subsequently characterised experimentally and genes with significant similarity to others with known functions. In addition, the structural annotation of genes has been significantly refined; 27% of gene structures have been changed since 2002, comprising changes in exon-intron boundaries, addition or deletion of exons and the addition or deletion of genes. The sequence has also undergone significant improvements. In addition to the correction of a large number of single-base and insertion or deletion errors, a major miss-assembly between the subtelomeres of chromosome 7 and 8 has been corrected. As the number of sequenced isolates continues to grow rapidly, a single reference genome will not be an adequate basis for interpreting intra-species sequence diversity. We therefore describe in this publication a population reference genome of P. falciparum, called Pfref1. This reference will enable the community to map to regions that are not present in the current assembly. P. falciparum 3D7 will continue to be maintained, with ongoing curation ensuring continual improvements in annotation quality.

New insights from Opisthorchis felineus genome: update on genomics of the epidemiologically important liver flukes

Background

The three epidemiologically important Opisthorchiidae liver flukes Opisthorchis felineus, O. viverrini, and Clonorchis sinensis, are believed to harbour similar potencies to provoke hepatobiliary diseases in their definitive hosts, although their populations have substantially different ecogeographical aspects including habitat, preferred hosts, population structure. Lack of O. felineus genomic data is an obstacle to the development of comparative molecular biological approaches necessary to obtain new knowledge about the biology of Opisthorchiidae trematodes, to identify essential pathways linked to parasite-host interaction, to predict genes that contribute to liver fluke pathogenesis and for the effective prevention and control of the disease.

Results

Here we present the first draft genome assembly of O. felineus and its gene repertoire accompanied by a comparative analysis with that of O. viverrini and Clonorchis sinensis. We observed both noticeably high heterozygosity of the sequenced individual and substantial genetic diversity in a pooled sample. This indicates that potency of O. felineus population for rapid adaptive response to control and preventive measures of opisthorchiasis is higher than in O. viverrini and C. sinensis. We also have found that all three species are characterized by more intensive involvement of trans-splicing in RNA processing compared to other trematodes.

Conclusion

All revealed peculiarities of structural organization of genomes are of extreme importance for a proper description of genes and their products in these parasitic species. This should be taken into account both in academic and applied research of epidemiologically important liver flukes. Further comparative genomics studies of liver flukes and non-carcinogenic flatworms allow for generation of well-grounded hypotheses on the mechanisms underlying development of cholangiocarcinoma associated with opisthorchiasis and clonorchiasis as well as species-specific mechanisms of these diseases.

Electronic supplementary material

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

Genomic analysis of Leptospira interrogans serovar Paidjan and Dadas isolates from carrier dogs and comparative genomic analysis to detect genes under positive selection

BACKGROUND

Leptospirosis is an emerging infectious disease worldwide that can cause high morbidity and mortality rates in humans and animals. The causative spirochetes have reservoirs in mammalian hosts, but there has been limited analysis of the genomes of isolates recovered from animals. The aims of this study were to characterize genomic features of two Leptospira interrogans strains recently isolated from asymptomatic dogs in Thailand (strains CUDO5 and CDUO8), and to perform comparative genome analyses with other strains. Molecular adaptive evolution in L. interrogans as signaled by positive selection also was analyzed.

RESULTS

Whole genome sequence analysis revealed that strains CUDO5 and CUDO8 had genome sizes of approximately 4.9 Mbp with 35.1% GC contents. Using monoclonal antibodies, strains CUDO5 and CUDO8 were identified as serovars Paidjan and Dadas, respectively. These strains harbored genes known to be associated with acute and chronic infections. Using Single Nucleotide Polymorphisms phylogeny (SNPs) with 97 L. interrogans strains, CUDO5 and CUDO8 had closest genetic relatedness with each other. Nevertheless, the serovar determinant region (rfb locus) showed variations in the genes encoding sugar biosynthesis. Amongst 13 representative L. interrogans strains examined for molecular adaptive evolution through positive selection under the site-model of Phylogenetic Analysis of Maximum Likelihood, genes responsible for iron acquisition (tlyA and hbpA), motility (fliN2, flgK, and flhB) and thermal adaptation (lpxD1) were under increased selective pressure.

CONCLUSIONS

L. interrogans serovar Paidjan strain CUDO5 and serovar Dadas strain CUDO8 had close genetic relatedness as analyzed by SNPs phylogeny. They contained genes with established roles in acute and chronic leptospirosis. The rfb locus in both serovars showed gene variation associated with sugar biosynthesis. Positive selection analysis indicated that genes encoding factors involved in motility, temperature adaptation, and iron acquisition were under strong positive selection in L. interrogans. These may be associated with adaptation in the early stages of infection.

Novel genetic code and record-setting AT-richness in the highly reduced plastid genome of the holoparasitic plant Balanophora

Plastid genomes (plastomes) vary enormously in size and gene content among the many lineages of nonphotosynthetic plants, but key lineages remain unexplored. We therefore investigated plastome sequence and expression in the holoparasitic and morphologically bizarre Balanophoraceae. The two Balanophora plastomes examined are remarkable, exhibiting features rarely if ever seen before in plastomes or in any other genomes. At 15.5 kb in size and with only 19 genes, they are among the most reduced plastomes known. They have no tRNA genes for protein synthesis, a trait found in only three other plastid lineages, and thus Balanophora plastids must import all tRNAs needed for translation. Balanophora plastomes are exceptionally compact, with numerous overlapping genes, highly reduced spacers, loss of all cis-spliced introns, and shrunken protein genes. With A+T contents of 87.8% and 88.4%, the Balanophora genomes are the most AT-rich genomes known save for a single mitochondrial genome that is merely bloated with AT-rich spacer DNA. Most plastid protein genes in Balanophora consist of ≥90% AT, with several between 95% and 98% AT, resulting in the most biased codon usage in any genome described to date. A potential consequence of its radical compositional evolution is the novel genetic code used by Balanophora plastids, in which TAG has been reassigned from stop to tryptophan. Despite its many exceptional properties, the Balanophora plastome must be functional because all examined genes are transcribed, its only intron is correctly trans-spliced, and its protein genes, although highly divergent, are evolving under various degrees of selective constraint.

Staphylococcus aureus from patients with chronic rhinosinusitis show minimal genetic association between polyp and non-polyp phenotypes

Background

Staphylococcus aureus has a high prevalence in chronic rhinosinusitis (CRS) patients and is suggested to play a more etiopathogenic role in CRS patients with nasal polyps (CRSwNP), a severe form of the CRS spectrum with poorer surgical outcomes. We performed a microbial genome-wide association study (mGWAS) to investigate whether S. aureus isolates from CRS patients have particular genetic markers associated with CRS with nasal polyps (CRSwNP) or CRS without nasal polyps (CRSsNP).

Methods

Whole genome sequencing was performed on S. aureus isolates collected from 28 CRSsNP and 30 CRSwNP patients. A mGWAS approach was employed using large-scale comparative genomics to identify genetic variation within our dataset.

Results

Considerable genetic variation was observed, with > 90,000 single nucleotide polymorphisms (SNPs) sites identified. There was little correlation with CRS subtype based on SNPs and Insertion/Delection (Indels). One indel was found to significantly correlate with CRSwNP and occurred in the promoter region of a bacitracin transport system ATP-binding protein. Additionally, two variants of the highly variable superantigen-like (SSL) proteins were found to significantly correlate with each CRS phenotype. No significant association with other virulence or antibiotic resistance genes were observed, consistent with previous studies.

Conclusion

To our knowledge this study is the first to use mGWAS to investigate the contribution of microbial genetic variation to CRS presentations. Utilising the most comprehensive genome-wide analysis methods available, our results suggest that CRS phenotype may be influenced by genetic factors other than specific virulence mechanisms within the S. aureus genome.

Electronic supplementary material

The online version of this article (10.1186/s12901-018-0064-1) contains supplementary material, which is available to authorized users.

An Introduced Crop Plant Is Driving Diversification of the Virulent Bacterial Pathogen Erwinia tracheiphila

Erwinia tracheiphila is the causal agent of bacterial wilt of cucurbits, an economically important phytopathogen affecting an economically important phytopathogen affecting few cultivated Cucurbitaceae few cultivated Cucurbitaceae host plant species in temperate eastern North America. However, essentially nothing is known about E. tracheiphila population structure or genetic diversity. To address this shortcoming, a representative collection of 88 E. tracheiphila isolates was gathered from throughout its geographic range, and their genomes were sequenced. Phylogenomic analysis revealed three genetic clusters with distinct hrpT3SS virulence gene repertoires, host plant association patterns, and geographic distributions. Low genetic heterogeneity within each cluster suggests a recent population bottleneck followed by population expansion. We showed that in the field and greenhouse, cucumber (Cucumis sativus), which was introduced to North America by early Spanish conquistadors, is the most susceptible host plant species and the only species susceptible to isolates from all three lineages. The establishment of large agricultural populations of highly susceptible C. sativus in temperate eastern North America may have facilitated the original emergence of E. tracheiphila into cucurbit agroecosystems, and this introduced plant species may now be acting as a highly susceptible reservoir host. Our findings have broad implications for agricultural sustainability by drawing attention to how worldwide crop plant movement, agricultural intensification, and locally unique environments may affect the emergence, evolution, and epidemic persistence of virulent microbial pathogens.IMPORTANCEErwinia tracheiphila is a virulent phytopathogen that infects two genera of cucurbit crop plants, Cucurbita spp. (pumpkin and squash) and Cucumis spp. (muskmelon and cucumber). One of the unusual ecological traits of this pathogen is that it is limited to temperate eastern North America. Here, we complete the first large-scale sequencing of an E. tracheiphila isolate collection. From phylogenomic, comparative genomic, and empirical analyses, we find that introduced Cucumis spp. crop plants are driving the diversification of E. tracheiphila into multiple lineages. Together, the results from this study show that locally unique biotic (plant population) and abiotic (climate) conditions can drive the evolutionary trajectories of locally endemic pathogens in unexpected ways.

Understanding Mucor circinelloides pathogenesis by comparative genomics and phenotypical studies

The increasing number of infections by species of Mucorales and their high mortality constitute an important concern for public health. This study aims to decipher the genetic basis of Mucor circinelloides pathogenicity, which displays virulence in a strain dependent manner. Assuming that genetic differences between strains may be linked to different pathotypes, we have conducted a study to explore genes responsible for virulence in M. circinelloides by whole genome sequencing of the avirulent strain NRRL3631 and comparison with the virulent strain CBS277.49. This genome analysis revealed 773 truncated, discontiguous and absent genes in the NRRL3631 strain. We also examined phenotypic traits resulting in reduced heat stress tolerance, chitosan content and lower susceptibility to toxic compounds (calcofluor white and sodium dodecyl sulphate) in the virulent strain, suggesting the influence of cell wall on pathogenesis. Based on these results, we focused on studying extracellular protein-coding genes by gene deletion and further pathotype characterization of mutants in murine models of pulmonary and systemic infection. Deletion of gene ID112092, which codes for a hypothetical extracellular protein of unknown function, resulted in significant reduction of virulence.

Although pathogenesis is a multifactorial process, these findings highlight the crucial role of surface and secreted proteins in M. circinelloides virulence and should promote further studies of other differential genes.

Genomic characterization of NDM-1 and 5, and OXA-181 carbapenemases in uropathogenic Escherichia coli isolates from Riyadh, Saudi Arabia

Urinary tract infections (UTIs) associated with Escherichia coli are a growing threat with an increase in the prevalence of multidrug resistant (MDR) strains, particularly ß-lactamase producers, occurring globally. We investigated the presence of carbapenem-resistant uropathogenic E. coli clones in community-acquired UTIs in Riyadh, Kingdom of Saudi Arabia (KSA) to identify the virulence and resistance structures of the resistant clones and relate the isolates to those circulating globally. A combination of comparative genomics and phenotypic approaches were used to characterize ten MDR-uropathogenic Escherichia coli isolates recovered from UTI patients in Riyadh between November 2014 and January 2015. We report the presence of NDM-1 and 5, and OXA-181 in carbapenem-resistant UPEC strains from Riyadh, KSA. Single nucleotide polymorphism analyses demonstrated that these ten isolates fell into four phylogenetically distinct clades within the UPEC phylogeny. Comparative genomic analyses indicate that these diverse clones could be distinguished according to their multilocus sequencing type (MLST), serology, and virulence and antimicrobial gene architectures. These clones include the blaNDM-1 carrying isolates of the globally predominant MDR ST131 and ST69 types, previously identified as one of the most common UPEC strains in KSA. This is in addition to clones of ST23Cplx (ST410) and ST448Cplx (ST448) that have likely evolved from common intestinal strains, carrying copies of ß-lactamase genes including blaNDM-5, blaCTX-M-15, blaTEM-1, blaCMY-42, blaOXA-1 and blaOXA-181. These data have identified an emerging public health concern and highlight the need to use comprehensive approaches to detect the structure of MDR E. coli populations associated with community-acquired UTIs in KSA.

The origin and adaptive evolution of domesticated populations of yeast from Far East Asia

The yeast Saccharomyces cerevisiae has been an essential component of human civilization because of its long global history of use in food and beverage fermentation. However, the diversity and evolutionary history of the domesticated populations of the yeast remain elusive. We show here that China/Far East Asia is likely the center of origin of the domesticated populations of the species. The domesticated populations form two major groups associated with solid- and liquid-state fermentation and appear to have originated from heterozygous ancestors, which were likely formed by outcrossing between diverse wild isolates primitively for adaptation to maltose-rich niches. We found consistent gene expansion and contraction in the whole domesticated population, as well as lineage-specific genome variations leading to adaptation to different environments. We show a nearly panoramic view of the diversity and life history of S. cerevisiae and provide new insights into the origin and evolution of the species.

Network-guided genomic and metagenomic analysis of the faecal microbiota of the critically endangered kakapo

The kakapo is a critically endangered, herbivorous parrot endemic to New Zealand. The kakapo hindgut hosts a dense microbial community of low taxonomic diversity, typically dominated by Escherichia fergusonii, and has proven to be a remarkably stable ecosystem, displaying little variation in core membership over years of study. To elucidate mechanisms underlying this robustness, we performed 16S rRNA gene-based co-occurrence network analysis to identify potential interactions between E. fergusonii and the wider bacterial community. Genomic and metagenomic sequencing were employed to facilitate interpretation of potential interactions observed in the network. E. fergusonii maintained very few correlations with other members of the microbiota, and isolates possessed genes for the generation of energy from a wide range of carbohydrate sources, including plant fibres such as cellulose. We surmise that this dominant microorganism is abundant not due to ecological interaction with other members of the microbiota, but its ability to metabolise a wide range of nutrients in the gut. This research represents the first concerted effort to understand the functional roles of the kakapo microbiota, and leverages metagenomic data to contextualise co-occurrence patterns. By combining these two techniques we provide a means for studying the diversity-stability hypothesis in the context of bacterial ecosystems.

Draft Genome Sequence of Staphylococcus microti DSM 22147, Isolated from the Common Vole

Staphylococcus microti DSM 22147 was isolated from viscera of common voles (Microtus arvalis Pallas) with generalized Brucella microti infection in the Czech Republic. To the best of our knowledge, the genome sequence of the species S. microti has not been previously studied. The complete genome sequence of strain DSM 22147 includes a genome of 2,381,859 bp (38.0% GC content) without any plasmids.

Complete Genome Sequence of Spiroplasma monobiae MQ-1T (ATCC 33825), a Bacterium Isolated from the Vespid Wasp (Monobia quadridens)

Spiroplasma monobiae MQ-1T (ATCC 33825) was isolated from the hemolymph of an adult vespid wasp (Monobia quadridens) collected in Maryland. Here, we report the complete genome sequence of this bacterium to facilitate the investigation of its biology and the comparative genomics among Spiroplasma species.

Genomic characterization of Nitrospirillum amazonense strain CBAmC, a nitrogen-fixing bacterium isolated from surface-sterilized sugarcane stems

Nitrospirillum amazonense is a nitrogen-fixing bacterium that shows potential to promote plant growth when inoculated into sugarcane and rice plants. This microorganism has been the subject of biochemical and genetic characterization to elucidate important functions related to host plant interaction and growth promotion, including the determination of draft genome sequences of two strains, Y2 and CBAmC, the second of which is the aim of the present study. CBAmC has been isolated from sugarcane (Saccharum spp.), and is currently used in a sugarcane consortium inoculant with four other nitrogen-fixing bacterial strains. The present paper describes a significant improvement in the genome sequence and assembly for the N. amazonense strain CBAmC, and determination for the first time of a complete genome sequence for this bacterial species, using PacBio technology. The analysis of the genomic data obtained allowed the discovery of genes coding for metabolic pathways and cellular structures that may be determinant for the success of the bacterial establishment and colonization into the host sugarcane plant, besides conferring important characteristics to the inoculant. These include genes for the use of sucrose and N-glycans, biosynthesis of autoinducer molecules, siderophore production and acquisition, auxin and polyamine biosynthesis, flagellum, σ-fimbriae, a variety of secretion systems, and a complete denitrification system. Concerning genes for nitrogenase and auxiliary proteins, it was possible to corroborate literature data that in N. amazonense these probably had originated from horizontal gene transfer, from bacteria of the Rhizobiales order. The complete genomic sequence of the CBAmC strain of N. amazonense revealed that the bacterium harbors four replicons, including three chromosomes and one chromid, a profile that coincides with that of other two strains, according to literature data, suggesting that as a replicon pattern for the species. Finally, results of phylogenomic analyses in this work support the recent reclassification of the species, separating it from the Azospirillum genus. More importantly, results of the present work shall guide subsequent studies on strain CBAmC as well as the development of a sugarcane inoculant.

Complete Genome Sequence of Agrobacterium tumefaciens 1D1609

Agrobacterium tumefaciens 1D1609 is a highly virulent strain isolated from a crown gall tumor of alfalfa (Medicago sativa L.). Compared to other well-characterized A. tumefaciens strains, such as C58 and Ach5, 1D1609 has a distinctive host range. Here, we report its complete genome sequence to facilitate future studies.

Complete Genome Sequence of Spiroplasma floricola 23-6 T (ATCC 29989), a Bacterium Isolated from a Tulip Tree (Liriodendron tulipifera L.)

Spiroplasma floricola 23-6^T (ATCC 29989) was isolated from the flower surface of a tulip tree (Liriodendron tulipifera L.). Here, we report the complete genome sequence of this bacterium to facilitate the investigation of its biology and the comparative genomics among Spiroplasma species.

Complete avian malaria parasite genomes reveal features associated with lineage-specific evolution in birds and mammals

Avian malaria parasites are prevalent around the world and infect a wide diversity of bird species. Here, we report the sequencing and analysis of high-quality draft genome sequences for two avian malaria species, Plasmodium relictum and Plasmodium gallinaceum We identify 50 genes that are specific to avian malaria, located in an otherwise conserved core of the genome that shares gene synteny with all other sequenced malaria genomes. Phylogenetic analysis suggests that the avian malaria species form an outgroup to the mammalian Plasmodium species, and using amino acid divergence between species, we estimate the avian- and mammalian-infective lineages diverged in the order of 10 million years ago. Consistent with their phylogenetic position, we identify orthologs of genes that had previously appeared to be restricted to the clades of parasites containing Plasmodium falciparum and Plasmodium vivax, the species with the greatest impact on human health. From these orthologs, we explore differential diversifying selection across the genus and show that the avian lineage is remarkable in the extent to which invasion-related genes are evolving. The subtelomeres of the P. relictum and P. gallinaceum genomes contain several novel gene families, including an expanded surf multigene family. We also identify an expansion of reticulocyte binding protein homologs in P. relictum, and within these proteins, we detect distinct regions that are specific to nonhuman primate, humans, rodent, and avian hosts. For the first time in the Plasmodium lineage, we find evidence of transposable elements, including several hundred fragments of LTR-retrotransposons in both species and an apparently complete LTR-retrotransposon in the genome of P. gallinaceum.

Comparative Genomics and Identification of an Enterotoxin-Bearing Pathogenicity Island, SEPI-1/SECI-1, in Staphylococcus epidermidis Pathogenic Strains

Staphylococcus epidermidis is a leading cause of nosocomial infections, majorly resistant to beta-lactam antibiotics, and may transfer several mobile genetic elements among the members of its own species, as well as to Staphylococcus aureus; however, a genetic exchange from S. aureus to S. epidermidis remains controversial. We recently identified two pathogenic clinical strains of S. epidermidis that produce a staphylococcal enterotoxin C3-like (SEC) similar to that by S. aureus pathogenicity islands. This study aimed to determine the genetic environment of the SEC-coding sequence and to identify the mobile genetic elements. Whole-genome sequencing and annotation of the S. epidermidis strains were performed using Illumina technology and a bioinformatics pipeline for assembly, which provided evidence that the SEC-coding sequences were located in a composite pathogenicity island that was previously described in the S. epidermidis strain FRI909, called SePI-1/SeCI-1, with 83.8–89.7% nucleotide similarity. Various other plasmids were identified, particularly p_3_95 and p_4_95, which carry antibiotic resistance genes (hsrA and dfrG, respectively), and share homologies with SAP085A and pUSA04-2-SUR11, two plasmids described in S. aureus. Eventually, one complete prophage was identified, ΦSE90, sharing 30 out of 52 coding sequences with the Acinetobacter phage vB_AbaM_IME200. Thus, the SePI-1/SeCI-1 pathogenicity island was identified in two pathogenic strains of S. epidermidis that produced a SEC enterotoxin causing septic shock. These findings suggest the existence of in vivo genetic exchange from S. aureus to S. epidermidis.

Whole genome sequencing of the monomorphic pathogen Mycobacterium bovis reveals local differentiation of cattle clinical isolates

Background

Bovine tuberculosis (bTB) poses serious risks to animal welfare and economy, as well as to public health as a zoonosis. Its etiological agent, Mycobacterium bovis, belongs to the Mycobacterium tuberculosis complex (MTBC), a group of genetically monomorphic organisms featured by a remarkably high overall nucleotide identity (99.9%). Indeed, this characteristic is of major concern for correct typing and determination of strain-specific traits based on sequence diversity. Due to its historical economic dependence on cattle production, Uruguay is deeply affected by the prevailing incidence of Mycobacterium bovis. With the world’s highest number of cattle per human, and its intensive cattle production, Uruguay represents a particularly suited setting to evaluate genomic variability among isolates, and the diversity traits associated to this pathogen.

Results

We compared 186 genomes from MTBC strains isolated worldwide, and found a highly structured population in M. bovis. The analysis of 23 new M. bovis genomes, belonging to strains isolated in Uruguay evidenced three groups present in the country. Despite presenting an expected highly conserved genomic structure and sequence, these strains segregate into a clustered manner within the worldwide phylogeny. Analysis of the non-pe/ppe differential areas against a reference genome defined four main sources of variability, namely: regions of difference (RD), variable genes, duplications and novel genes. RDs and variant analysis segregated the strains into clusters that are concordant with their spoligotype identities. Due to its high homoplasy rate, spoligotyping failed to reflect the true genomic diversity among worldwide representative strains, however, it remains a good indicator for closely related populations.

Conclusions

This study introduces a comprehensive population structure analysis of worldwide M. bovis isolates. The incorporation and analysis of 23 novel Uruguayan M. bovis genomes, sheds light onto the genomic diversity of this pathogen, evidencing the existence of greater genetic variability among strains than previously contemplated.

Electronic supplementary material

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

First insights into a type II toxin-antitoxin system from the clinical isolate Mycobacterium sp. MHSD3, similar to epsilon/zeta systems

A putative type II toxin-antitoxin (TA) system was found in the clinical isolate Mycobacterium sp. MHSD3, a strain closely related to Mycobacterium chelonae. Further analyses of the protein sequences of the two genes revealed the presence of domains related to a TA system. BLAST analyses indicated the presence of closely related proteins in the genomes of other recently published M. chelonae strains. The functionality of both elements of the TA system was demonstrated when expressed in Escherichia coli cells, and the predicted structure of the toxin is very similar to those of well-known zeta-toxins, leading to the definition of a type II TA system similar to epsilon/zeta TA systems in strains that are closely related to M. chelonae.

Whole-Genome Sequencing of Seven Strains of Staphylococcus lugdunensis Allows Identification of Mobile Genetic Elements

Coagulase negative staphylococci are normal inhabitant of the human skin flora that account for an increasing number of infections, particularly hospital-acquired infections. Staphylococcus lugdunensis has emerged as a most virulent species causing various infections with clinical characteristics close to what clinicians usually observe with Staphylococcus aureus and both bacteria share more than 70% of their genome. Virulence of S. aureus relies on a large repertoire of virulence factors, many of which are encoded on mobile genetic elements. S. lugdunensis also bears various putative virulence genes but only one complete genome with extensive analysis has been published with one prophage sequence (φSL2) and a unique plasmid was previously described. In this study, we performed de novo sequencing, whole genome assembly and annotation of seven strains of S. lugdunensis from VISLISI clinical trial. We searched for the presence of virulence genes and mobile genetics elements using bioinformatics tools. We identified four new prophages, named φSL2 to φSL4, belonging to the Siphoviridae class and five plasmids, named pVISLISI_1 to pVISLISI_5. Three plasmids are homologous to known plasmids that include, amongst others, one S. aureus plasmid. The two other plasmids were not described previously. This study provides a new context for the study of S. lugdunensis virulence suggesting the occurrence of several genetic recombination’ with other staphylococci.

Complete Genome Sequence of Spiroplasma corruscae EC-1 T (DSM 19793), a Bacterium Isolated from a Lampyrid Beetle (Ellychnia corrusca)

Spiroplasma corruscae EC-1^T (DSM 19793) was isolated from the gut of a lampryid beetle (Ellychnia corrusca) collected in Beltsville, MD, USA, in 1983. Here, we report the complete genome sequence of this bacterium to facilitate the investigation of its biology and the comparative genomics among Spiroplasma species.

Phylogeographic, genomic, and meropenem susceptibility analysis of Burkholderia ubonensis

The bacterium Burkholderia ubonensis is commonly co-isolated from environmental specimens harbouring the melioidosis pathogen, Burkholderia pseudomallei. B. ubonensis has been reported in northern Australia and Thailand but not North America, suggesting similar geographic distribution to B. pseudomallei. Unlike most other Burkholderia cepacia complex (Bcc) species, B. ubonensis is considered non-pathogenic, although its virulence potential has not been tested. Antibiotic resistance in B. ubonensis, particularly towards drugs used to treat the most severe B. pseudomallei infections, has also been poorly characterised. This study examined the population biology of B. ubonensis, and includes the first reported isolates from the Caribbean. Phylogenomic analysis of 264 B. ubonensis genomes identified distinct clades that corresponded with geographic origin, similar to B. pseudomallei. A small proportion (4%) of strains lacked the 920kb chromosome III replicon, with discordance of presence/absence amongst genetically highly related strains, demonstrating that the third chromosome of B. ubonensis, like other Bcc species, probably encodes for a nonessential pC3 megaplasmid. Multilocus sequence typing using the B. pseudomallei scheme revealed that one-third of strains lack the "housekeeping" narK locus. In comparison, all strains could be genotyped using the Bcc scheme. Several strains possessed high-level meropenem resistance (≥32 μg/mL), a concern due to potential transmission of this phenotype to B. pseudomallei. In silico analysis uncovered a high degree of heterogeneity among the lipopolysaccharide O-antigen cluster loci, with at least 35 different variants identified. Finally, we show that Asian B. ubonensis isolate RF23-BP41 is avirulent in the BALB/c mouse model via a subcutaneous route of infection. Our results provide several new insights into the biology of this understudied species.

A Case Study into Microbial Genome Assembly Gap Sequences and Finishing Strategies

This study characterized regions of DNA which remained unassembled by either PacBio and Illumina sequencing technologies for seven bacterial genomes. Two genomes were manually finished using bioinformatics and PCR/Sanger sequencing approaches and regions not assembled by automated software were analyzed. Gaps present within Illumina assemblies mostly correspond to repetitive DNA regions such as multiple rRNA operon sequences. PacBio gap sequences were evaluated for several properties such as GC content, read coverage, gap length, ability to form strong secondary structures, and corresponding annotations. Our hypothesis that strong secondary DNA structures blocked DNA polymerases and contributed to gap sequences was not accepted. PacBio assemblies had few limitations overall and gaps were explained as cumulative effect of lower than average sequence coverage and repetitive sequences at contig termini. An important aspect of the present study is the compilation of biological features that interfered with assembly and included active transposons, multiple plasmid sequences, phage DNA integration, and large sequence duplication. Our targeted genome finishing approach and systematic evaluation of the unassembled DNA will be useful for others looking to close, finish, and polish microbial genome sequences.

Draft Genome Sequence of Klebsiella michiganensis 3T412C, Harboring an Arsenic Resistance Genomic Island, Isolated from Mine Tailings in Peru

An arsenic resistance genomic island in the bacterium Klebsiella michiganensis 3T412C was isolated from mine tailings from Peru. This genomic island confers adaptation to extreme environments with high concentrations of arsenic. Isolate 3T412C contained a complete set of genes involved in resistance to arsenic. This operon is surrounded by putative genes for resistance to other heavy metals.

Approaches for in silico finishing of microbial genome sequences

The introduction of next-generation sequencing (NGS) had a significant effect on the availability of genomic information, leading to an increase in the number of sequenced genomes from a large spectrum of organisms. Unfortunately, due to the limitations implied by the short-read sequencing platforms, most of these newly sequenced genomes remained as "drafts", incomplete representations of the whole genetic content. The previous genome sequencing studies indicated that finishing a genome sequenced by NGS, even bacteria, may require additional sequencing to fill the gaps, making the entire process very expensive. As such, several in silico approaches have been developed to optimize the genome assemblies and facilitate the finishing process. The present review aims to explore some free (open source, in many cases) tools that are available to facilitate genome finishing.

An improved Plasmodium cynomolgi genome assembly reveals an unexpected methyltransferase gene expansion

Background: Plasmodium cynomolgi, a non-human primate malaria parasite species, has been an important model parasite since its discovery in 1907. Similarities in the biology of P. cynomolgi to the closely related, but less tractable, human malaria parasite P. vivax make it the model parasite of choice for liver biology and vaccine studies pertinent to P. vivax malaria. Molecular and genome-scale studies of P. cynomolgi have relied on the current reference genome sequence, which remains highly fragmented with 1,649 unassigned scaffolds and little representation of the subtelomeres. 

Methods: Using long-read sequence data (Pacific Biosciences SMRT technology), we assembled and annotated a new reference genome sequence, PcyM, sourced from an Indian rhesus monkey. We compare the newly assembled genome sequence with those of several other Plasmodium species, including a re-annotated P. coatneyi assembly.

Results: The new PcyM genome assembly is of significantly higher quality than the existing reference, comprising only 56 pieces, no gaps and an improved average gene length. Detailed manual curation has ensured a comprehensive annotation of the genome with 6,632 genes, nearly 1,000 more than previously attributed to P. cynomolgi. The new assembly also has an improved representation of the subtelomeric regions, which account for nearly 40% of the sequence. Within the subtelomeres, we identified more than 1300 Plasmodium interspersed repeat ( pir) genes, as well as a striking expansion of 36 methyltransferase pseudogenes that originated from a single copy on chromosome 9.

Conclusions: The manually curated PcyM reference genome sequence is an important new resource for the malaria research community. The high quality and contiguity of the data have enabled the discovery of a novel expansion of methyltransferase in the subtelomeres, and illustrates the new comparative genomics capabilities that are being unlocked by complete reference genomes.

The invasive Neisseria meningitidis MenC CC103 from Brazil is characterized by an accessory gene repertoire

Neisseria meningitidis infections are a major issue for global health. The invasive MenC ST-103 clonal complex (CC103) has been the most prevalent in meningococcal outbreaks in Brazil, occurring also in several countries worldwide. Here we have analysed the population structure and accessory genome of MenC CC103 strains from a global perspective. An in-depth phylogenomic analysis revealed a lineage of N. meningitidis causing meningitis in Brazil and the United Kingdom. This lineage was also characterized as harbouring a particular accessory genome composed of CRISPR/Cas and restriction modification systems. This lineage was also characterized by a genomic island resembling an integrative and conjugative element. This island carried genes potentially associated with virulence and fitness. We propose this accessory gene repertoire could be contributing to the spatial-temporal persistence of the invasive MenC CC103 lineage.

Plasmodium falciparum Maf1 Confers Survival upon Amino Acid Starvation

The target of rapamycin complex 1 (TORC1) pathway is a highly conserved signaling pathway across eukaryotes that integrates nutrient and stress signals to regulate the cellular growth rate and the transition into and maintenance of dormancy. The majority of the pathway's components, including the central TOR kinase, have been lost in the apicomplexan lineage, and it is unknown how these organisms detect and respond to nutrient starvation in its absence. Plasmodium falciparum encodes a putative ortholog of the RNA polymerase (Pol) III repressor Maf1, which has been demonstrated to modulate Pol III transcription in a TOR-dependent manner in a number of organisms. Here, we investigate the role of P. falciparum Maf1 (PfMaf1) in regulating RNA Pol III expression under conditions of nutrient starvation and other stresses. Using a transposon insertion mutant with an altered Maf1 expression profile, we demonstrated that proper Maf1 expression is necessary for survival of the dormancy-like state induced by prolonged amino acid starvation and is needed for full recovery from other stresses that slow or stall the parasite cell cycle. This Maf1 mutant is defective in the downregulation of pre-tRNA synthesis under nutrient-limiting conditions, indicating that the function of Maf1 as a stress-responsive regulator of structural RNA transcription is conserved in P. falciparum Recent work has demonstrated that parasites carrying artemisinin-resistant K13 alleles display an enhanced ability to recover from drug-induced growth retardation. We show that one such artemisinin-resistant line displays greater regulation of pre-tRNA expression and higher survival upon prolonged amino acid starvation, suggesting that overlapping, PfMaf1-associated pathways may regulate growth recovery from both artemisinin treatment and amino acid starvation.IMPORTANCE Eukaryote organisms sense changes in their environment and integrate this information through signaling pathways to activate response programs to ensure survival. The TOR pathway is a well-studied signaling pathway found throughout eukaryotes that is known to integrate a variety of signals to regulate organismal growth in response to starvation and other stresses. The human malaria parasite Plasmodium falciparum appears to have lost the TOR pathway over the course of evolution, and it is unclear how the parasite modulates its growth in response to starvation and drug treatment. Here, we show that Maf1, a protein regulated by TOR in other eukaryotes, plays an important role in maintaining the parasite's viability in the face of starvation and other forms of stress. This suggests that PfMaf1 is a component of a yet-to-be-described nutrient and stress response pathway.