Mugsy: fast multiple alignment of closely related whole genomes

Efficient merging of genome profile alignments

Motivation

Whole-genome alignment (WGA) methods show insufficient scalability toward the generation of large-scale WGAs. Profile alignment-based approaches revolutionized the fields of multiple sequence alignment construction methods by significantly reducing computational complexity and runtime. However, WGAs need to consider genomic rearrangements between genomes, which make the profile-based extension of several whole-genomes challenging. Currently, none of the available methods offer the possibility to align or extend WGA profiles.

Results

Here, we present genome profile alignment, an approach that aligns the profiles of WGAs and that is capable of producing large-scale WGAs many times faster than conventional methods. Our concept relies on already available whole-genome aligners, which are used to compute several smaller sets of aligned genomes that are combined to a full WGA with a divide and conquer approach. To align or extend WGA profiles, we make use of the SuperGenome data structure, which features a bidirectional mapping between individual sequence and alignment coordinates. This data structure is used to efficiently transfer different coordinate systems into a common one based on the principles of profiles alignments. The approach allows the computation of a WGA where alignments are subsequently merged along a guide tree. The current implementation uses progressiveMauve and offers the possibility for parallel computation of independent genome alignments. Our results based on various bacterial datasets up to several hundred genomes show that we can reduce the runtime from months to hours with a quality that is negligibly worse than the WGA computed with the conventional progressiveMauve tool.

Availability and implementation

GPA is freely available at https://lambda.informatik.uni-tuebingen.de/gitlab/ahennig/GPA. GPA is implemented in Java, uses progressiveMauve and offers a parallel computation of WGAs.

Supplementary information

Supplementary data are available at Bioinformatics online.

Nanopore sequencing data analysis: state of the art, applications and challenges

The nanopore sequencing process is based on the transit of a DNA molecule through a nanoscopic pore, and since the 90s is considered as one of the most promising approaches to detect polymeric molecules. In 2014, Oxford Nanopore Technologies (ONT) launched a beta-testing program that supplied the scientific community with the first prototype of a nanopore sequencer: the MinION. Thanks to this program, several research groups had the opportunity to evaluate the performance of this novel instrument and develop novel computational approaches for analyzing this new generation of data. Despite the short period of time from the release of the MinION, a large number of algorithms and tools have been developed for base calling, data handling, read mapping, de novo assembly and variant discovery. Here, we face the main computational challenges related to the analysis of nanopore data, and we carry out a comprehensive and up-to-date survey of the algorithmic solutions adopted by the bioinformatic community comparing performance and reporting limits and advantages of using this new generation of sequences for genomic analyses. Our analyses demonstrate that the use of nanopore data dramatically improves the de novo assembly of genomes and allows for the exploration of structural variants with an unprecedented accuracy and resolution. However, despite the impressive improvements reached by ONT in the past 2 years, the use of these data for small-variant calling is still challenging, and at present, it needs to be coupled with complementary short sequences for mitigating the intrinsic biases of nanopore sequencing technology.

OMMA enables population-scale analysis of complex genomic features and phylogenomic relationships from nanochannel-based optical maps

BACKGROUND

Optical mapping is an emerging technology that complements sequencing-based methods in genome analysis. It is widely used in improving genome assemblies and detecting structural variations by providing information over much longer (up to 1 Mb) reads. Current standards in optical mapping analysis involve assembling optical maps into contigs and aligning them to a reference, which is limited to pairwise comparison and becomes bias-prone when analyzing multiple samples.

FINDINGS

We present a new method, OMMA, that extends optical mapping to the study of complex genomic features by simultaneously interrogating optical maps across many samples in a reference-independent manner. OMMA captures and characterizes complex genomic features, e.g., multiple haplotypes, copy number variations, and subtelomeric structures when applied to 154 human samples across the 26 populations sequenced in the 1000 Genomes Project. For small genomes such as pathogenic bacteria, OMMA accurately reconstructs the phylogenomic relationships and identifies functional elements across 21 Acinetobacter baumannii strains.

CONCLUSIONS

With the increasing data throughput of optical mapping system, the use of this technology in comparative genome analysis across many samples will become feasible. OMMA is a timely solution that can address such computational need. The OMMA software is available at https://github.com/TF-Chan-Lab/OMTools.

Lineage specific histories of Mycobacterium tuberculosis dispersal in Africa and Eurasia

Mycobacterium tuberculosis (M.tb) is a globally distributed, obligate pathogen of humans that can be divided into seven clearly defined lineages. An emerging consensus places the origin and global dispersal of M.tb within the past 6,000 years: identifying how the ancestral clone of M.tb spread and differentiated within this timeframe is important for identifying the ecological drivers of the current pandemic. We used Bayesian phylogeographic inference to reconstruct the migratory history of M.tb in Africa and Eurasia and to investigate lineage specific patterns of spread from a geographically diverse sample of 552 M.tb genomes. Applying evolutionary rates inferred with ancient M.tb genome calibration, we estimated the timing of major events in the migratory history of the pathogen. Inferred timings contextualize M.tb dispersal within historical phenomena that altered patterns of connectivity throughout Africa and Eurasia: trans-Indian Ocean trade in spices and other goods, the Silk Road and its predecessors, the expansion of the Roman Empire, and the European Age of Exploration. We found that Eastern Africa and Southeast Asia have been critical in the dispersal of M.tb. Our results further reveal that M.tb populations have grown through range expansion, as well as in situ, and delineate the independent evolutionary trajectories of bacterial subpopulations underlying the current pandemic.

Genetic changes associated with the temporal shift in invasive non-typhoidal Salmonella serovars in Bamako Mali

BACKGROUND

Invasive non-typhoidal Salmonella (iNTS) serovars S. Typhimurium and S. Enteritidis are major etiologic agents of invasive bacterial disease among infants and young children in sub-Saharan Africa, including in Mali. Early studies of iNTS serovars in several countries indicated that S. Typhimurium was more prevalent than S. Enteritidis, including in Mali before 2008. We investigated genomic and associated phenotypic changes associated with an increase in the relative proportion of iNTS caused by S. Enteritidis versus S. Typhimurium in Bamako, Mali, during the period 2002-2012.

METHODOLOGY/PRINCIPAL FINDINGS

Comparative genomics studies identified homologs of tetracycline resistance and arsenic utilization genes that were associated with the temporal shift of serovars causing iNTS shift, along with several hypothetical proteins. These findings, validated through PCR screening and phenotypic assays, provide initial steps towards characterizing the genomic changes consequent to unknown evolutionary pressures associated with the shift in serovar prevalence.

CONCLUSIONS/SIGNIFICANCE

This work identified a shift to S. Enteritidis from the more classic S. Typhimurium, associated with iNTS in Bamako, Mali, during the period 2002-2012. This type of shift in underlying iNTS pathogens are of great importance to pediatric public health in endemic regions of sub-Saharan Africa. Additionally, this work demonstrates the utility of combining epidemiologic data, whole genome sequencing, and functional characterization in the laboratory to identify and characterize genomic changes in the isolates that may be involved with the observed shift in circulating iNTS agents.

Clustering of Vibrio parahaemolyticus Isolates Using MLST and Whole-Genome Phylogenetics and Protein Motif Fingerprinting

Vibrio parahaemolyticus is a ubiquitous and abundant member of native microbial assemblages in coastal waters and shellfish. Though V. parahaemolyticus is predominantly environmental, some strains have infected human hosts and caused outbreaks of seafood-related gastroenteritis. In order to understand differences among clinical and environmental V. parahaemolyticus strains, we used high quality DNA sequencing data to compare the genomes of V. parahaemolyticus isolates (n = 43) from a variety of geographic locations and clinical and environmental sample matrices. We used phylogenetic trees inferred from multilocus sequence typing (MLST) and whole-genome (WG) alignments, as well as a novel classification and genome clustering approach that relies on protein motif fingerprints (MFs), to assess relationships between V. parahaemolyticus strains and identify novel molecular targets associated with virulence. Differences in strain clustering at more than one position were observed between the MLST and WG phylogenetic trees. The WG phylogeny had higher support values and strain resolution since isolates of the same sequence type could be differentiated. The MF analysis revealed groups of protein motifs that were associated with the pathogenic MLST type ST36 and a large group of clinical strains isolated from human stool. A subset of the stool and ST36-associated protein motifs were selected for further analysis and the motif sequences were found in genes with a variety of functions, including transposases, secretion system components and effectors, and hypothetical proteins. DNA sequences associated with these protein motifs are candidate targets for future molecular assays in order to improve surveys of pathogenic V. parahaemolyticus in the environment and seafood.

Insight into the evolution and functional characteristics of the pan-genome assembly from sesame landraces and modern cultivars

Sesame (Sesamum indicum L.) is an important oil crop renowned for its high oil content and quality. Recently, genome assemblies for five sesame varieties including two landraces (S. indicum cv. Baizhima and Mishuozhima) and three modern cultivars (S. indicum var. Zhongzhi13, Yuzhi11 and Swetha), have become available providing a rich resource for comparative genomic analyses and gene discovery. Here, we employed a reference-assisted assembly approach to improve the draft assemblies of four of the sesame varieties. We then constructed a sesame pan-genome of 554.05 Mb. The pan-genome contained 26 472 orthologous gene clusters; 15 409 (58.21%) of them were core (present across all five sesame genomes), whereas the remaining 41.79% (11 063) clusters and the 15 890 variety-specific genes were dispensable. Comparisons between varieties suggest that modern cultivars from China and India display significant genomic variation. The gene families unique to the sesame modern cultivars contain genes mainly related to yield and quality, while those unique to the landraces contain genes involved in environmental adaptation. Comparative evolutionary analysis indicates that several genes involved in plant-pathogen interaction and lipid metabolism are under positive selection, which may be associated with sesame environmental adaption and selection for high seed oil content. This study of the sesame pan-genome provides insights into the evolution and genomic characteristics of this important oilseed and constitutes a resource for further sesame crop improvement.

Amplification Dynamics of Platy-1 Retrotransposons in the Cebidae Platyrrhine Lineage

Platy-1 elements are Platyrrhine-specific, short interspersed elements originally discovered in the Callithrix jacchus (common marmoset) genome. To date, only the marmoset genome has been analyzed for Platy-1 repeat content. Here, we report full-length Platy-1 insertions in other New World monkey (NWM) genomes (Saimiri boliviensis, squirrel monkey; Cebus imitator, capuchin monkey; and Aotus nancymaae, owl monkey) and analyze the amplification dynamics of lineage-specific Platy-1 insertions. A relatively small number of full-length and lineage-specific Platy-1 elements were found in the squirrel, capuchin, and owl monkey genomes compared with the marmoset genome. In addition, only a few older Platy-1 subfamilies were recovered in this study, with no Platy-1 subfamilies younger than Platy-1-6. By contrast, 62 Platy-1 subfamilies were discovered in the marmoset genome. All of the lineage-specific insertions found in the squirrel and capuchin monkeys were fixed present. However, ∼15% of the lineage-specific Platy-1 loci in Aotus were polymorphic for insertion presence/absence. In addition, two new Platy-1 subfamilies were identified in the owl monkey genome with low nucleotide divergences compared with their respective consensus sequences, suggesting minimal ongoing retrotransposition in the Aotus genus and no current activity in the Saimiri, Cebus, and Sapajus genera. These comparative analyses highlight the finding that the high number of Platy-1 elements discovered in the marmoset genome is an exception among NWM analyzed thus far, rather than the rule. Future studies are needed to expand upon our knowledge of Platy-1 amplification in other NWM genomes.

Genome sequence of Xanthomonas fuscans subsp. fuscans strain Xff49: a new isolate obtained from common beans in Southern Brazil

The genus Xanthomonas comprises Gram-negative bacteria, many of which are phytopathogens. Xanthomonas fuscans subsp. fuscans is one of the most devastating pathogens affecting the bean plant, resulting in the common bacterial blight of bean (CBB). The disease is mainly foliar and affects a wide variety of bean species, thus acting as the yield-limiting factor for the bean crop. Here, we report the whole-genome sequencing of a new strain of X. fuscans subsp. fuscans, named Xff49, isolated from the infected and symptomatic beans from Capão do Leão, Southern Brazil. The genetic analysis demonstrated the presence of single-nucleotide variants (SNVs) in this strain, potentially affecting the mobilome, cell mobility, and inorganic ion metabolism. In addition, the analysis resulted in the identification of a new plasmid similar to the pAX22 derived from Achromobacter denitrificans, which was named plX, along with plA and plC, previously reported in other strains of X. fuscans subsp. fuscans. Xff49 represents the first Brazilian genome of X. fuscans subsp. fuscans and might provide useful information applicable to the studies of phylogenetics, evolution, and pathogenomics, thereby allowing a better understanding of the genomic features present in the Brazilian strains.

Comparative genomic analysis of enterotoxigenic Escherichia coli O159 strains isolated from diarrheal patients in Korea

BACKGROUND

Enterotoxigenic Escherichia coli (ETEC) is a common cause of bacterial infection that leads to diarrhea. Although some studies have proposed a potential association between the toxic profile and genetic background, association between toxin of ETEC and phylo-group has not been reported yet. The objective of this study was to examine genomic and phylogenetic characteristics of ETEC strain NCCP15731 and NCCP15733 by whole genome sequencing and comparative genomic analysis of two phylo-groups of O159 reference strains.

RESULTS

Whole genome sequencing showed that genome size of NCCP15731 strain was 4,663,459 bp, containing 4435 CDS and 19 RNAs. The genome size of NCCP15733 was 4,645,336 bp, containing 4369 CDS and 23 RNAs. Both NCCP15731 and NCCP15733 were classified in the phylo-group A, which is one of major E. coli phylogenetic groups. Their serotype was O159:H34. They possessed the virulence factor such as adherence systems, auto transporter systems, and flagella segments of major driving force for ETEC pathogenicity. They also harbored STh enterotoxin. Hierarchical clustering result based on the presence or absence of a total of 108 major virulence factors of 14 O159 ETEC strains showed that seven strains in phylo-group A and seven strains in phylo-group B1 were clustered each other, respectively. Colonization factors (CFs) of NCCP15731 or NCCP15733 were not detected.

CONCLUSIONS

Serotype of NCCP15731 and NCCP15733, representing major types of ETEC in Korea, was O159:H34 and their MLST type was ST218. Comparison with other O159 strains revealed that NCCP15731 was specialized for transporter system and secretion system whereas NCCP15733 had unique genes related to capsular polysaccharide. Compared with E159, the most recent common ancestor, these two strains had different toxin type and virulence factors. These results will improve our understanding of ETEC O159 strains to prevent ETEC disease.

Whole-Genome Alignment and Comparative Annotation

Rapidly improving sequencing technology coupled with computational developments in sequence assembly are making reference-quality genome assembly economical. Hundreds of vertebrate genome assemblies are now publicly available, and projects are being proposed to sequence thousands of additional species in the next few years. Such dense sampling of the tree of life should give an unprecedented new understanding of evolution and allow a detailed determination of the events that led to the wealth of biodiversity around us. To gain this knowledge, these new genomes must be compared through genome alignment (at the sequence level) and comparative annotation (at the gene level). However, different alignment and annotation methods have different characteristics; before starting a comparative genomics analysis, it is important to understand the nature of, and biases and limitations inherent in, the chosen methods. This review is intended to act as a technical but high-level overview of the field that should provide this understanding. We briefly survey the state of the genome alignment and comparative annotation fields and potential future directions for these fields in a new, large-scale era of comparative genomics.

Draft Genome Sequences of Eight Vibrio sp. Clinical Isolates from across the United States That Form a Basal Sister Clade to Vibrio cholerae

We sequenced the genomes of eight isolates from various regions of the United States. These isolates form a monophyletic cluster clearly related to but distinct from Vibrio cholerae. Phylogenetic and genomic analyses suggest that they represent a basal lineage highly divergent from Vibrio cholerae or a novel species.

Highly flexible infection programs in a specialized wheat pathogen

Many filamentous plant pathogens exhibit high levels of genomic variability, yet the impact of this variation on host-pathogen interactions is largely unknown. We have addressed host specialization in the wheat pathogen Zymoseptoria tritici. Our study builds on comparative analyses of infection and gene expression phenotypes of three isolates and reveals the extent to which genomic variation translates into phenotypic variation. The isolates exhibit genetic and genomic variation but are similarly virulent. By combining confocal microscopy, disease monitoring, staining of ROS, and comparative transcriptome analyses, we conducted a detailed comparison of the infection processes of these isolates in a susceptible wheat cultivar. We characterized four core infection stages: establishment, biotrophic growth, lifestyle transition, and necrotrophic growth and asexual reproduction that are shared by the three isolates. However, we demonstrate differentiated temporal and spatial infection development and significant differences in the expression profiles of the three isolates during the infection stages. More than 20% of the genes were differentially expressed and these genes were located significantly closer to transposable elements, suggesting an impact of epigenetic regulation. Further, differentially expressed genes were enriched in effector candidates suggesting that isolate-specific strategies for manipulating host defenses are present in Z. tritici. We demonstrate that individuals of a host-specialized pathogen have highly differentiated infection programs characterized by flexible infection development and functional redundancy. This illustrates how high genetic diversity in pathogen populations results in highly differentiated infection phenotypes, which fact needs to be acknowledged to understand host-pathogen interactions and pathogen evolution.

Whole-Genome Alignment

Whole-genome alignment (WGA) is the prediction of evolutionary relationships at the nucleotide level between two or more genomes. It combines aspects of both colinear sequence alignment and gene orthology prediction and is typically more challenging to address than either of these tasks due to the size and complexity of whole genomes. Despite the difficulty of this problem, numerous methods have been developed for its solution because WGAs are valuable for genome-wide analyses such as phylogenetic inference, genome annotation, and function prediction. In this chapter, we discuss the meaning and significance of WGA and present an overview of the methods that address it. We also examine the problem of evaluating whole-genome aligners and offer a set of methodological challenges that need to be tackled in order to make most effective use of our rapidly growing databases of whole genomes.

Clinical, immunological and bacteriological characteristics of H7N9 patients nosocomially co-infected by Acinetobacter Baumannii: a case control study

BACKGROUND

Bacterial co-infection of patients suffering from influenza pneumonia is a key element that increases morbidity and mortality. The occurrence of Acinetobacter baumannii co-infection in patients with avian influenza A (H7N9) virus infection has been described as one of the most prevalent bacterial co-infections. However, the clinical and laboratory features of this entity of H7N9 and A. baumannii co-infection have not been systematically investigated.

METHODS

We collected clinical and laboratory data from laboratory-confirmed H7N9 cases co-infected by A. baumannii. H7N9 patients without bacterial co-infection and patients with A. baumannii-related pneumonia in the same hospital during the same period were recruited as controls. The antibiotic resistance features and the corresponding genome determinants of A. baumannii and the immune responses of the patients were tested through the respiratory and peripheral blood specimens.

RESULTS

Invasive mechanical ventilation was the most significant risk factor for the nosocomial A. baumannii co-infection in H7N9 patients. The co-infection resulted in severe clinical manifestation which was associated with the dysregulation of immune responses including deranged T-cell counts, antigen-specific T-cell responses and plasma cytokines. The emergence of genome variations of extensively drug-resistant A. baumannii associated with acquired polymyxin resistance contributed to the fatal outcome of a co-infected patient.

CONCLUSIONS

The co-infection of H7N9 patients by extensively drug-resistant A. baumannii with H7N9 infection is an important issue which deserves attention. The dysfunctions of immune responses were associated with the co-infection and were correlated with the disease severity. These data provide useful reference for the diagnosis and treatment of H7N9 infection.

Genome-wide identification and expression analysis of StTCP transcription factors of potato (Solanum tuberosum L.)

The plant-specific TCP transcription factors, which play critical roles in diverse aspects of biological processes, have been identified and analyzed in various plant species. However, no systematical study of TCP family genes in potato (Solanum tuberosum L.) has been undertaken. In this study, a total of 31 non-redundant TCP transcription factors of potato were identified and divided into two subfamilies including three distinct subclades. The various orthologous TCP genes in Arabidopsis, rice, potato and tomato were identified using synteny and phylogenetic analysis. Protein motif analysis demonstrated that StTCPs in the same subclade shared similar conserved motif structures. Gene structure analysis showed that almost all StTCPs displayed highly conserved exon-intron organization. The analysis of StTCP gene promoter regions revealed that multiple cis-acting elements were involved in plant growth, development, hormone responses as well as stress responses. The result of StTCP gene expression profiles showed they had tissue-specific expression patterns which implied their differentiated functions. According to the results of quantitative RT-PCR (qRT-PCR), 7 StTCP genes were dramatically up-regulated during the release of tuber dormancy and some specific StTCP genes were strongly responding to different abiotic stresses and multiple hormones, which suggested they had important roles in potato growth and development processes. The results of our findings could provide comprehensive insights in StTCP family genes of potato for further functional investigations.

Draft Genome Sequences of Nine Vibrio sp. Isolates from across the United States Closely Related to Vibrio cholerae

We are reporting whole-genome sequences of nine Vibrio sp. isolates closely related to the waterborne human pathogen Vibrio cholerae.

We are reporting whole-genome sequences of nine Vibrio sp. isolates closely related to the waterborne human pathogen Vibrio cholerae. These isolates were recovered from sources, including human samples, from different regions of the United States. Genome analysis suggests that this group of isolates represents a highly divergent basal V. cholerae lineage or a closely related novel species.

Comprehensive study of instable regions in Pseudomonas aeruginosa and Mycobacterium tuberculosis

Background

Pseudomonas aeruginosa is a common bacterium which is recognized for its association with hospital-acquired infections and its advanced antibiotic resistance mechanisms. Tuberculosis, one of the major causes of mortality, is initiated by the deposition of Mycobacterium tuberculosis. Accessory sequences shared by a subset of strains of a species play an important role in a species’ evolution, antibiotic resistance and infectious potential.

Results

Here, with a multiple sequence aligner, we segmented 25 P. aeruginosa genomes and 28 M. tuberculosis genomes into core blocks (include sequences shared by all the input genomes) and dispensable blocks (include sequences shared by a subset of the input genomes), respectively. For each input genome, we then constructed a scaffold consisting of its core and dispensable blocks sorted by blocks’ locations on the chromosomes. Consecutive dispensable blocks on these scaffold formed instable regions. After a comprehensive study of these instable regions, three characteristics of instable regions are summarized: instable regions were short, site specific and varied in different strains. Three DNA elements (directed repeats (DRs), transposons and integrons) were then studied to see whether these DNA elements are associated with the variation of instable regions. A pipeline was developed to search for DR pairs on the flank of every instable sequence. 27 DR pairs in P. aeruginosa strains and 6 pairs in M. tuberculosis strains were found to exist in the instable regions. On the average, 14% and 12% of instable regions in P. aeruginosa strains covered transposase genes and integrase genes, respectively. In M. tuberculosis strains, an average of 43% and 8% of instable regions contain transposase genes and integrase genes, respectively.

Conclusions

Instable regions were short, site specific and varied in different strains for both P. aeruginosa and M. tuberculosis. Our experimental results showed that DRs, transposons and integrons may be associated with variation of instable regions.

Electronic supplementary material

The online version of this article (10.1186/s12938-018-0563-8) contains supplementary material, which is available to authorized users.

Assembly of a Complete Mitogenome of Chrysanthemum nankingense Using Oxford Nanopore Long Reads and the Diversity and Evolution of Asteraceae Mitogenomes

Diversity in structure and organization is one of the main features of angiosperm mitochondrial genomes (mitogenomes). The ultra-long reads of Oxford Nanopore Technology (ONT) provide an opportunity to obtain a complete mitogenome and investigate the structural variation in unprecedented detail. In this study, we compared mitogenome assembly methods using Illumina and/or ONT sequencing data and obtained the complete mitogenome (208 kb) of Chrysanthemum nankingense based on the hybrid assembly method. The mitogenome encoded 19 transfer RNA genes, three ribosomal RNA genes, and 34 protein-coding genes with 21 group II introns disrupting eight intron-contained genes. A total of seven medium repeats were related to homologous recombination at different frequencies as supported by the long ONT reads. Subsequently, we investigated the variations in gene content and constitution of 28 near-complete mitogenomes from Asteraceae. A total of six protein-coding genes were missing in all Asteraceae mitogenomes, while four other genes were not detected in some lineages. The core fragments (~88 kb) of the Asteraceae mitogenomes had a higher GC content (~46.7%) than the variable and specific fragments. The phylogenetic topology based on the core fragments of the Asteraceae mitogenomes was highly consistent with the topologies obtained from the corresponding plastid datasets. Our results highlighted the advantages of the complete assembly of the C. nankingense mitogenome and the investigation of its structural variation based on ONT sequencing data. Moreover, the method based on local collinear blocks of the mitogenomes could achieve the alignment of highly rearrangeable and variable plant mitogenomes as well as construct a robust phylogenetic topology.

Genome sequence of the potato pathogenic fungus Alternaria solani HWC-168 reveals clues for its conidiation and virulence

BACKGROUND

Alternaria solani is a known air-born deuteromycete fungus with a polycyclic life cycle and is the causal agent of early blight that causes significant yield losses of potato worldwide. However, the molecular mechanisms underlying the conidiation and pathogenicity remain largely unknown.

RESULTS

We produced a high-quality genome assembly of A. solani HWC-168 that was isolated from a major potato-producing region of Northern China, which facilitated a comprehensive gene annotation, the accurate prediction of genes encoding secreted proteins and identification of conidiation-related genes. The assembled genome of A. solani HWC-168 has a genome size 32.8 Mb and encodes 10,358 predicted genes that are highly similar with related Alternaria species including Alternaria arborescens and Alternaria brassicicola. We identified conidiation-related genes in the genome of A. solani HWC-168 by searching for sporulation-related homologues identified from Aspergillus nidulans. A total of 975 secreted protein-encoding genes, which might act as virulence factors, were identified in the genome of A. solani HWC-168. The predicted secretome of A. solani HWC-168 possesses 261 carbohydrate-active enzymes (CAZy), 119 proteins containing RxLx[EDQ] motif and 27 secreted proteins unique to A. solani.

CONCLUSIONS

Our findings will facilitate the identification of conidiation- and virulence-related genes in the genome of A. solani. This will permit new insights into understanding the molecular mechanisms underlying the A. solani-potato pathosystem and will add value to the global fungal genome database.

Using Core Genome Alignments To Assign Bacterial Species

With the exponential increase in the number of bacterial taxa with genome sequence data, a new standardized method to assign species designations is needed that is consistent with classically obtained taxonomic analyses. This is particularly acute for unculturable, obligate intracellular bacteria with which classically defined methods, like DNA-DNA hybridization, cannot be used, such as those in the Rickettsiales. In this study, we generated nucleotide-based core genome alignments for a wide range of genera with classically defined species, as well as those within the Rickettsiales. We created a workflow that uses the length, sequence identity, and phylogenetic relationships inferred from core genome alignments to assign genus and species designations that recapitulate classically obtained results. Using this method, most classically defined bacterial genera have a core genome alignment that is ≥10% of the average input genome length. Both Anaplasma and Neorickettsia fail to meet this criterion, indicating that the taxonomy of these genera should be reexamined. Consistently, genomes from organisms with the same species epithet have ≥96.8% identity of their core genome alignments. Additionally, these core genome alignments can be used to generate phylogenomic trees to identify monophyletic clades that define species and neighbor-network trees to assess recombination across different taxa. By these criteria, Wolbachia organisms are delineated into species different from the currently used supergroup designations, while Rickettsia organisms are delineated into 9 distinct species, compared to the current 27 species. By using core genome alignments to assign taxonomic designations, we aim to provide a high-resolution, robust method to guide bacterial nomenclature that is aligned with classically obtained results. IMPORTANCE With the increasing availability of genome sequences, we sought to develop and apply a robust, portable, and high-resolution method for the assignment of genera and species designations that can recapitulate classically defined taxonomic designations. Using cutoffs derived from the lengths and sequence identities of core genome alignments along with phylogenetic analyses, we sought to evaluate or reevaluate genus- and species-level designations for diverse taxa, with an emphasis on the order Rickettsiales, where species designations have been applied inconsistently. Our results indicate that the Rickettsia genus has an overabundance of species designations, that the current Anaplasma and Neorickettsia genus designations are both too broad and need to be divided, and that there are clear demarcations of Wolbachia species that do not align precisely with the existing supergroup designations.

Evolution of Olfactory Functions on the Fire Ant Social Chromosome

Understanding the molecular evolutionary basis of social behavior is a major challenge in evolutionary biology. Social insects evolved a complex language of chemical signals to coordinate thousands of individuals. In the fire ant Solenopsis invicta, chemical signals are involved in the determination of a polymorphic social organization. Single-queen (monogyne) or multiqueen (polygyne) social structure is determined by the "social chromosome," a nonrecombining region containing ∼504 genes with two distinct haplotypes, SB and Sb. Monogyne queens are always SBB, while polygyne queens are always SBb. Workers discriminate monogyne from polygyne queens based on olfactory cues. Here, we took an evolutionary genomics approach to search for candidate genes in the social chromosome that could be responsible for this discrimination. We compared the SB and Sb haplotypes and analyzed the evolutionary rates since their divergence. Notably, we identified a cluster of 23 odorant receptors in the nonrecombining region of the social chromosome that stands out in terms of nonsynonymous changes in both haplotypes. The cluster includes twelve genes formed by recent Solenopsis-specific duplications. We found evidence for positive selection on several tree branches and significant differences between the SB and Sb haplotypes of these genes. The most dramatic difference is the complete deletion of two of these genes in Sb. These results suggest that the evolution of polygyne social organization involved adaptations in olfactory genes and opens the way for functional studies of the molecular mechanisms underlying social behavior.

Comparative genomics of the Erwinia and Enterobacter olive fly endosymbionts

The pestivorous tephritid olive fly has long been known as a frequent host of the obligately host-associated bacterial endosymbiont, Erwinia dacicola, as well as other facultative endosymbionts. The genomes of Erwinia dacicola and Enterobacter sp. OLF, isolated from a California olive fly, encode the ability to supplement amino acids and vitamins missing from the olive fruit on which the larvae feed. The Enterobacter sp. OLF genome encodes both uricase and ureases, and the Er. dacicola genome encodes an allantoate transport pathway, suggesting that bird feces or recycling the fly's waste products may be important sources of nitrogen. No homologs to known nitrogenases were identified in either bacterial genome, despite suggestions of their presence from experiments with antibiotic-treated flies. Comparisons between the olive fly endosymbionts and their free-living relatives revealed similar GC composition and genome size. The Er. dacicola genome has fewer genes for amino acid metabolism, cell motility, and carbohydrate transport and metabolism than free-living Erwinia spp. while having more genes for cell division, nucleotide metabolism and replication as well as mobile elements. A 6,696 bp potential lateral gene transfer composed primarily of amino acid synthesis and transport genes was identified that is also observed in Pseudomonas savastanoii pv savastanoii, the causative agent of olive knot disease.

Recombinational DSBs-intersected genes converge on specific disease- and adaptability-related pathways

Motivation

The budding yeast Saccharomyces cerevisiae is a model species powerful for studying the recombination of eukaryotes. Although many recombination studies have been performed for this species by experimental methods, the population genomic study based on bioinformatics analyses is urgently needed to greatly increase the range and accuracy of recombination detection. Here, we carry out the population genomic analysis of recombination in S.cerevisiae to reveal the potential rules between recombination and evolution in eukaryotes.

Results

By population genomic analysis, we discover significantly more and longer recombination events in clinical strains, which indicates that adverse environmental conditions create an obviously wider range of genetic combination in response to the selective pressure. Based on the analysis of recombinational double strand breaks (DSBs)-intersected genes (RDIGs), we find that RDIGs significantly converge on specific disease- and adaptability-related pathways, indicating that recombination plays a biologically key role in the repair of DSBs related to diseases and environmental adaptability, especially the human neurological disorders. By evolutionary analysis of RDIGs, we find that the RDIGs highly prevailing in populations of yeast tend to be more evolutionarily conserved, indicating the accurate repair of DSBs in these RDIGs is critical to ensure the eukaryotic survival or fitness.

Supplementary information

Supplementary data are available at Bioinformatics online.

Determination and genome-wide analysis of Epstein-Barr virus (EBV) sequences in EBV-associated gastric carcinoma from Guangdong, an endemic area of nasopharyngeal carcinoma

About 10 % of gastric carcinoma worldwide is associated with EBV, which is defined as EBV-associated gastric carcinoma (EBVaGC). To date, EBV sequence data from EBVaGC in Guangdong, China, an endemic area of nasopharyngeal carcinoma (NPC), are not available. In the present study, two EBV genomes from EBVaGC specimens from Guangdong (designated as GDGC1 and GDGC2) were determined by next-generation sequencing, de novo assembly and joining of contigs by Sanger sequencing. In addition, we sequenced EBV from two Korean EBVaGC cell lines, YCCEL1 and SNU-719. Genomic diversity, including single nucleotide polymorphisms (SNPs) and insertions and deletions (indels), phylogenetic analysis and rates of protein evolution, was performed using bioinformatics software. The four gastric carcinoma-derived EBV (GC-EBV) were all type I. Compared with the reference EBV genome, a total of 1815 SNPs (146 indels), 1519 SNPs (106 indels), 1812 SNPs (126 indels) and 1484 SNPs (106 indels) were found in GDGC1, GDGC2, YCCEL1 and SNU-719, respectively. These variations were distributed across the entire genome, especially in latent genes. In contrast, the sequences of promoters and non-coding RNAs were strictly conserved. Phylogenetic analyses suggested the presence of at least two parental lineages of EBV among the GC-EBV genomes. Rates of protein evolution analyses showed that lytic genes were under purifying selection; in contrast, latency genes were under positive selection. In conclusion, this study determined the EBV genomes in EBVaGC from Guangdong and performed a detailed genome-wide analysis of GC-EBV, which would be helpful for further understanding of the relationship between EBV genomic variation and EBVaGC carcinogenesis.

Genome Plasticity of agr-Defective Staphylococcus aureus during Clinical Infection

Therapy for bacteremia caused by Staphylococcus aureus is often ineffective, even when treatment conditions are optimal according to experimental protocols. Adapted subclones, such as those bearing mutations that attenuate agr-mediated virulence activation, are associated with persistent infection and patient mortality.

Therapy for bacteremia caused by Staphylococcus aureus is often ineffective, even when treatment conditions are optimal according to experimental protocols. Adapted subclones, such as those bearing mutations that attenuate agr-mediated virulence activation, are associated with persistent infection and patient mortality. To identify additional alterations in agr-defective mutants, we sequenced and assembled the complete genomes of clone pairs from colonizing and infected sites of several patients in whom S. aureus demonstrated a within-host loss of agr function. We report that events associated with agr inactivation result in agr-defective blood and nares strain pairs that are enriched in mutations compared to pairs from wild-type controls. The random distribution of mutations between colonizing and infecting strains from the same patient, and between strains from different patients, suggests that much of the genetic complexity of agr-defective strains results from prolonged infection or therapy-induced stress. However, in one of the agr-defective infecting strains, multiple genetic changes resulted in increased virulence in a murine model of bloodstream infection, bypassing the mutation of agr and raising the possibility that some changes were selected. Expression profiling correlated the elevated virulence of this agr-defective mutant to restored expression of the agr-regulated ESAT6-like type VII secretion system, a known virulence factor. Thus, additional mutations outside the agr locus can contribute to diversification and adaptation during infection by S. aureus agr mutants associated with poor patient outcomes.

Complete genome sequence of soil actinobacteria Streptomyces cavourensis TJ430

A new actinobacteria Streptomyces cavourensis TJ430 was isolated from the mountain soil collected from the southwest of China. In previous study, TJ430 showed striking bactericidal activities and strong ability of antibiotic production. Here, we report complete genome of this bacterium, consisting of 7.6 Mb linear chromosome and 0.2 Mb plasmids. It was predicted 6450 genes in chromosome and 225 genes in plasmids, as well as 12 gene islands in chromosome. Abundant genes have predicted functions in antibiotic metabolism and stress resistance. A whole-genome comparison of S. cavourensis TJ430, S. coelicolor A3(2), and S. lividans 66 indicates that TJ430 has a relatively high degree of strain specificity. The 16S rRNA phylogenetic tree shows the high identities (99.79%) of TJ430 with S. cavourensis DSM40300. TJ430 is a new and rare Streptomyces species, and analysis of its genome helps us to better understand primary metabolism mechanism of this isolate, as well as the evolutionary biology.

NovoGraph: Human genome graph construction from multiple long-read de novo assemblies

Genome graphs are emerging as an important novel approach to the analysis of high-throughput human sequencing data. By explicitly representing genetic variants and alternative haplotypes in a mappable data structure, they can enable the improved analysis of structurally variable and hyperpolymorphic regions of the genome. In most existing approaches, graphs are constructed from variant call sets derived from short-read sequencing. As long-read sequencing becomes more cost-effective and enables de novo assembly for increasing numbers of whole genomes, a method for the direct construction of a genome graph from sets of assembled human genomes would be desirable. Such assembly-based genome graphs would encompass the wide spectrum of genetic variation accessible to long-read-based de novo assembly, including large structural variants and divergent haplotypes. Here we present NovoGraph, a method for the construction of a human genome graph directly from a set of de novo assemblies. NovoGraph constructs a genome-wide multiple sequence alignment of all input contigs and creates a graph by merging the input sequences at positions that are both homologous and sequence-identical. NovoGraph outputs resulting graphs in VCF format that can be loaded into third-party genome graph toolkits. To demonstrate NovoGraph, we construct a genome graph with 23,478,835 variant sites and 30,582,795 variant alleles from de novo assemblies of seven ethnically diverse human genomes (AK1, CHM1, CHM13, HG003, HG004, HX1, NA19240). Initial evaluations show that mapping against the constructed graph reduces the average mismatch rate of reads from sample NA12878 by approximately 0.2%, albeit at a slightly increased rate of reads that remain unmapped.

Closed Genome and Comparative Phylogenetic Analysis of the Clinical Multidrug Resistant Shigella sonnei Strain 866

Shigella sonnei is responsible for the majority of shigellosis infections in the US with over 500,000 cases reported annually. Here, we present the complete genome of the clinical multidrug resistant (MDR) strain 866, which is highly susceptible to bacteriophage infections. The strain has a circular chromosome of 4.85 Mb and carries a 113 kb MDR plasmid. This IncB/O/K/Z-type plasmid, termed p866, confers resistance to five different classes of antibiotics including ß-lactamase, sulfonamide, tetracycline, aminoglycoside, and trimethoprim. Comparative analysis of the plasmid architecture and gene inventory revealed that p866 shares its plasmid backbone with previously described IncB/O/K/Z-type Shigella spp. and Escherichia coli plasmids, but is differentiated by the insertion of antibiotic resistance cassettes, which we found associated with mobile genetic elements such as Tn3, Tn7, and Tn10. A whole genome-derived phylogenetic reconstruction showed the evolutionary relationships of S. sonnei strain 866 and the four established Shigella species, highlighting the clonal nature of S. sonnei.

Isolation of a divergent strain of Rickettsia japonica from Dew's Australian bat Argasid ticks (Argas (Carios) dewae) in Victoria, Australia

A divergent strain of Rickettsia japonica was isolated from a Dew's Australian bat argasid tick, Argas (Carios) dewae, collected in southern Victoria, Australia and a full-genome analysis along with sequencing of 5 core gene fragments was undertaken. This isolate was designated Rickettsia japonica str. argasii (ATCC VR-1665, CSUR R179).

Approximate, simultaneous comparison of microbial genome architectures via syntenic anchoring of quiver representations

Motivation

A long-standing limitation in comparative genomic studies is the dependency on a reference genome, which hinders the spectrum of genetic diversity that can be identified across a population of organisms. This is especially true in the microbial world where genome architectures can significantly vary. There is therefore a need for computational methods that can simultaneously analyze the architectures of multiple genomes without introducing bias from a reference.

Results

In this article, we present Ptolemy: a novel method for studying the diversity of genome architectures-such as structural variation and pan-genomes-across a collection of microbial assemblies without the need of a reference. Ptolemy is a 'top-down' approach to compare whole genome assemblies. Genomes are represented as labeled multi-directed graphs-known as quivers-which are then merged into a single, canonical quiver by identifying 'gene anchors' via synteny analysis. The canonical quiver represents an approximate, structural alignment of all genomes in a given collection encoding structural variation across (sub-) populations within the collection. We highlight various applications of Ptolemy by analyzing structural variation and the pan-genomes of different datasets composing of Mycobacterium, Saccharomyces, Escherichia and Shigella species. Our results show that Ptolemy is flexible and can handle both conserved and highly dynamic genome architectures. Ptolemy is user-friendly-requires only FASTA-formatted assembly along with a corresponding GFF-formatted file-and resource-friendly-can align 24 genomes in ∼10 mins with four CPUs and <2 GB of RAM.

Availability and implementation

Github: https://github.com/AbeelLab/ptolemy.

Supplementary information

Supplementary data are available at Bioinformatics online.

Comparative Genomics and Description of Putative Virulence Factors of Melissococcus plutonius, the Causative Agent of European Foulbrood Disease in Honey Bees

In Europe, approximately 84% of cultivated crop species depend on insect pollinators, mainly bees. Apis mellifera (the Western honey bee) is the most important commercial pollinator worldwide. The Gram-positive bacterium Melissococcus plutonius is the causative agent of European foulbrood (EFB), a global honey bee brood disease. In order to detect putative virulence factors, we sequenced and analyzed the genomes of 14 M. plutonius strains, including two reference isolates. The isolates do not show a high diversity in genome size or number of predicted protein-encoding genes, ranging from 2.021 to 2.101 Mbp and 1589 to 1686, respectively. Comparative genomics detected genes that might play a role in EFB pathogenesis and ultimately in the death of the honey bee larvae. These include bacteriocins, bacteria cell surface- and host cell adhesion-associated proteins, an enterococcal polysaccharide antigen, an epsilon toxin, proteolytic enzymes, and capsule-associated proteins. In vivo expression of three putative virulence factors (endo-alpha-N-acetylgalactosaminidase, enhancin and epsilon toxin) was verified using naturally infected larvae. With our strain collection, we show for the first time that genomic differences exist between non-virulent and virulent typical strains, as well as a highly virulent atypical strain, that may contribute to the virulence of M. plutonius. Finally, we also detected a high number of conserved pseudogenes (75 to 156) per genome, which indicates genomic reduction during evolutionary host adaptation.

GRSR: a tool for deriving genome rearrangement scenarios from multiple unichromosomal genome sequences

Background

Genome rearrangements describe changes in the genetic linkage relationship of large chromosomal regions, involving reversals, transpositions, block interchanges, deletions, insertions, fissions, fusions and translocations etc. Many algorithms for calculating rearrangement scenarios between two genomes have been proposed. Very often, the calculated rearrangement scenario is not unique for the same pair of permutations. Hence, how to decide which calculated rearrangement scenario is more biologically meaningful becomes an essential task. Up to now, several mechanisms for genome rearrangements have been studied. One important theory is that genome rearrangement may be mediated by repeats, especially for reversal events. Many reversal regions are found to be flanked by a pair of inverted repeats. As a result, whether there are repeats at the breakpoints of the calculated rearrangement events can shed a light on deciding whether the calculated rearrangement events is biologically meaningful. To our knowledge, there is no tool which can automatically identify rearrangement events and check whether there exist repeats at the breakpoints of each calculated rearrangement event.

Results

In this paper, we describe a new tool named GRSR which allows us to compare multiple unichromosomal genomes to identify “independent” (obvious) rearrangement events such as reversals, (inverted) block interchanges and (inverted) transpositions and automatically searches for repeats at the breakpoints of each rearrangement event. We apply our tool on the complete genomes of 28 Mycobacterium tuberculosis strains and 24 Shewanella strains respectively. In both Mycobacterium tuberculosis and Shewanella strains, our tool finds many reversal regions flanked by a pair of inverted repeats. In particular, the GRSR tool also finds an inverted transposition and an inverted block interchange in Shewanella, where the repeats at the ends of rearrangement regions remain unchanged after the rearrangement event. To our knowledge, this is the first time such a phenomenon for inverted transposition and inverted block interchange is reported in Shewanella.

Conclusions

From the calculated results, there are many examples supporting the theory that the existence of repeats at the breakpoints of a rearrangement event can make the sequences at the breakpoints remain unchanged before and after the rearrangement events, suggesting that the conservation of ends could possibly be a popular phenomenon in many types of genome rearrangement events.

Draft genome sequence of the Wolbachia endosymbiont of Wuchereria bancrofti wWb

The draft genome assembly of the Wolbachia endosymbiont of Wuchereria bancrofti (wWb) consists of 1060 850 bp in 100 contigs and contains 961 ORFs, with a single copy of the 5S rRNA, 16S rRNA and 23S rRNA and each of the 34 tRNA genes. Phylogenetic core genome analyses show wWb to cluster with other strains in supergroup D of the Wolbachia phylogeny, while being most closely related to the Wolbachia endosymbiont of Brugia malayi strain TRS (wBm). The wWb and wBm genomes share 779 orthologous clusters with wWb having 101 unclustered genes and wBm having 23 unclustered genes. The higher number of unclustered genes in the wWb genome likely reflects the fragmentation of the draft genome.

Global Scale Dissemination of ST93: A Divergent Staphylococcus aureus Epidemic Lineage That Has Recently Emerged From Remote Northern Australia

Background: In Australia, community-associated methicillin-resistant Staphylococcus aureus (MRSA) lineage sequence type (ST) 93 has rapidly risen to dominance since being described in the early 1990s. We examined 459 ST93 genome sequences from Australia, New Zealand, Samoa, and Europe to investigate the evolutionary history of ST93, its emergence in Australia and subsequent spread overseas. Results: Comparisons with other S. aureus genomes indicate that ST93 is an early diverging and recombinant lineage, comprising of segments from the ST59/ST121 lineage and from a divergent but currently unsampled Staphylococcal population. However, within extant ST93 strains limited genetic diversity was apparent with the most recent common ancestor dated to 1977 (95% highest posterior density 1973-1981). An epidemic ST93 population arose from a methicillin-susceptible progenitor in remote Northern Australia, which has a proportionally large Indigenous population, with documented overcrowded housing and a high burden of skin infection. Methicillin-resistance was acquired three times in these regions, with a clade harboring a staphylococcal cassette chromosome mec (SCCmec) IVa expanding and spreading to Australia's east coast by 2000. We observed sporadic and non-sustained introductions of ST93-MRSA-IVa to the United Kingdom. In contrast, in New Zealand, ST93-MRSA-IVa was sustainably transmitted with clonal expansion within the Pacific Islander population, who experience similar disadvantages as Australian Indigenous populations. Conclusion: ST93 has a highly recombinant genome including portions derived from an early diverging S. aureus population. Our findings highlight the need to understand host population factors in the emergence and spread of antimicrobial resistant community pathogens.

Fern genomes elucidate land plant evolution and cyanobacterial symbioses

Ferns are the closest sister group to all seed plants, yet little is known about their genomes other than that they are generally colossal. Here, we report on the genomes of Azolla filiculoides and Salvinia cucullata (Salviniales) and present evidence for episodic whole-genome duplication in ferns-one at the base of 'core leptosporangiates' and one specific to Azolla. One fern-specific gene that we identified, recently shown to confer high insect resistance, seems to have been derived from bacteria through horizontal gene transfer. Azolla coexists in a unique symbiosis with N2-fixing cyanobacteria, and we demonstrate a clear pattern of cospeciation between the two partners. Furthermore, the Azolla genome lacks genes that are common to arbuscular mycorrhizal and root nodule symbioses, and we identify several putative transporter genes specific to Azolla-cyanobacterial symbiosis. These genomic resources will help in exploring the biotechnological potential of Azolla and address fundamental questions in the evolution of plant life.

Ultrasensitive Capture of Human Herpes Simplex Virus Genomes Directly from Clinical Samples Reveals Extraordinarily Limited Evolution in Cell Culture

Herpes simplex viruses (HSVs) are difficult to sequence due to their large DNA genome, high GC content, and the presence of repeats. To date, most HSV genomes have been recovered from culture isolates, raising concern that these genomes may not accurately represent circulating clinical strains. We report the development and validation of a DNA oligonucleotide hybridization panel to recover nearly complete HSV genomes at abundances up to 50,000-fold lower than previously reported. Using copy number information on herpesvirus and host DNA background via quantitative PCR, we developed a protocol for pooling for cost-effective recovery of more than 50 HSV-1 or HSV-2 genomes per MiSeq run. We demonstrate the ability to recover >99% of the HSV genome at >100× coverage in 72 h at viral loads that allow whole-genome recovery from latently infected ganglia. We also report a new computational pipeline for rapid HSV genome assembly and annotation. Using the above tools and a series of 17 HSV-1-positive clinical swabs sent to our laboratory for viral isolation, we show limited evolution of HSV-1 during viral isolation in human fibroblast cells compared to the original clinical samples. Our data indicate that previous studies using low-passage-number clinical isolates of herpes simplex viruses are reflective of the viral sequences present in the lesion and thus can be used in phylogenetic analyses. We also detect superinfection within a single sample with unrelated HSV-1 strains recovered from separate oral lesions in an immunosuppressed patient during a 2.5-week period, illustrating the power of direct-from-specimen sequencing of HSV.IMPORTANCE Herpes simplex viruses affect more than 4 billion people across the globe, constituting a large burden of disease. Understanding the global diversity of herpes simplex viruses is important for diagnostics and therapeutics as well as cure research and tracking transmission among humans. To date, most HSV genomics has been performed on culture isolates and DNA swabs with high quantities of virus. We describe the development of wet-lab and computational tools that enable the accurate sequencing of near-complete genomes of HSV-1 and HSV-2 directly from clinical specimens at abundances >50,000-fold lower than previously sequenced and at significantly reduced cost. We use these tools to profile circulating HSV-1 strains in the community and illustrate limited changes to the viral genome during the viral isolation process. These techniques enable cost-effective, rapid sequencing of HSV-1 and HSV-2 genomes that will help enable improved detection, surveillance, and control of this human pathogen.

Genomic analysis of a Raoultella ornithinolytica strain causing prosthetic joint infection in an immunocompetent patient

We sequenced the genome of Raoultella ornithinolytica strain Marseille-P1025 that caused a rare case of prosthetic joint infection in a 67-year-old immunocompetent male. The 6.7-Mb genome exhibited a genomic island (RoGI) that was unique among R. ornithinolytica strains. RoGI was likely acquired by lateral gene transfer from a member of the Pectobacterium genus and coded for a type IVa secretion system found in other pathogenic bacteria and that may have conferred strain Marseille-P1025 an increased virulence. Strain Marseille-P1025 was also able to infect, multiply within, and kill Acanthamoaeba castellanii amoebae.

Phylogenetic conservatism of thermal traits explains dispersal limitation and genomic differentiation of Streptomyces sister-taxa

The latitudinal diversity gradient is a pattern of biogeography observed broadly in plants and animals but largely undocumented in terrestrial microbial systems. Although patterns of microbial biogeography across broad taxonomic scales have been described in a range of contexts, the mechanisms that generate biogeographic patterns between closely related taxa remain incompletely characterized. Adaptive processes are a major driver of microbial biogeography, but there is less understanding of how microbial biogeography and diversification are shaped by dispersal limitation and drift. We recently described a latitudinal diversity gradient of species richness and intraspecific genetic diversity in Streptomyces by using a geographically explicit culture collection. Within this geographically explicit culture collection, we have identified Streptomyces sister-taxa whose geographic distribution is delimited by latitude. These sister-taxa differ in geographic distribution, genomic diversity, and ecological traits despite having nearly identical SSU rRNA gene sequences. Comparative genomic analysis reveals genomic differentiation of these sister-taxa consistent with restricted gene flow across latitude. Furthermore, we show phylogenetic conservatism of thermal traits between the sister-taxa suggesting that thermal trait adaptation limits dispersal and gene flow across climate regimes as defined by latitude. Such phylogenetic conservatism of thermal traits is commonly associated with latitudinal diversity gradients for plants and animals. These data provide further support for the hypothesis that the Streptomyces latitudinal diversity gradient was formed as a result of historical demographic processes defined by dispersal limitation and driven by paleoclimate dynamics.

ClustAGE: a tool for clustering and distribution analysis of bacterial accessory genomic elements

BACKGROUND

The non-conserved accessory genome of bacteria can be associated with important adaptive characteristics that can contribute to niche specificity or pathogenicity of strains. High degrees of structural and compositional diversity in genomic islands and other elements of the accessory genome can complicate characterization of accessory genome contents among populations of strains. Methods for easily and effectively defining the distributions of discrete elements of the accessory genome among bacterial strains in a population are needed to explore the relationships between the flexible genome and bacterial adaptive traits.

RESULTS

We have developed the open-source software package ClustAGE. This program, written in Perl, uses BLAST to cluster nucleotide accessory genomic elements from the genomes of multiple bacterial strains and to identify their distribution within the study population. The program output can be used in combination with strain phenotype data or other characteristics to detect associations. Optional graphical output is available for visualizing accessory genome gene content and distribution patterns. The capabilities of the software are demonstrated on a collection of 14 Pseudomonas aeruginosa genome sequences.

CONCLUSIONS

The ClustAGE software and utilities are effective for identifying characteristics and distributions of accessory genomic elements among groups of bacterial genomes. The ability to easily and effectively characterize the accessory genome of a sequence collection may provide a better understanding of the accessory genome's contribution to a species' adaptation and pathogenesis. The ClustAGE source code can be downloaded from https://clustage.sourceforge.io and a limited web-based implementation is available at http://vfsmspineagent.fsm.northwestern.edu/cgi-bin/clustage.cgi .

Natural selection and the predictability of evolution in Timema stick insects

Predicting evolution remains difficult. We studied the evolution of cryptic body coloration and pattern in a stick insect using 25 years of field data, experiments, and genomics. We found that evolution is more difficult to predict when it involves a balance between multiple selective factors and uncertainty in environmental conditions than when it involves feedback loops that cause consistent back-and-forth fluctuations. Specifically, changes in color-morph frequencies are modestly predictable through time (r² = 0.14) and driven by complex selective regimes and yearly fluctuations in climate. In contrast, temporal changes in pattern-morph frequencies are highly predictable due to negative frequency-dependent selection (r² = 0.86). For both traits, however, natural selection drives evolution around a dynamic equilibrium, providing some predictability to the process.

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.

Primordial origin and diversification of plasmids in Lyme disease agent bacteria

BACKGROUND

With approximately one-third of their genomes consisting of linear and circular plasmids, the Lyme disease agent cluster of species has the most complex genomes among known bacteria. We report here a comparative analysis of plasmids in eleven Borreliella (also known as Borrelia burgdorferi sensu lato) species.

RESULTS

We sequenced the complete genomes of two B. afzelii, two B. garinii, and individual B. spielmanii, B. bissettiae, B. valaisiana and B. finlandensis isolates. These individual isolates carry between seven and sixteen plasmids, and together harbor 99 plasmids. We report here a comparative analysis of these plasmids, along with 70 additional Borreliella plasmids available in the public sequence databases. We identify only one new putative plasmid compatibility type (the 30th) among these 169 plasmid sequences, suggesting that all or nearly all such types have now been discovered. We find that the linear plasmids in the non-B. burgdorferi species have undergone the same kinds of apparently random, chaotic rearrangements mediated by non-homologous recombination that we previously discovered in B. burgdorferi. These rearrangements occurred independently in the different species lineages, and they, along with an expanded chromosomal phylogeny reported here, allow the identification of several whole plasmid transfer events among these species. Phylogenetic analyses of the plasmid partition genes show that a majority of the plasmid compatibility types arose early, most likely before separation of the Lyme agent Borreliella and relapsing fever Borrelia clades, and this, with occasional cross species plasmid transfers, has resulted in few if any species-specific or geographic region-specific Borreliella plasmid types.

CONCLUSIONS

The primordial origin and persistent maintenance of the Borreliella plasmid types support their functional indispensability as well as evolutionary roles in facilitating genome diversity. The improved resolution of Borreliella plasmid phylogeny based on conserved partition-gene clusters will lead to better determination of gene orthology which is essential for prediction of biological function, and it will provide a basis for inferring detailed evolutionary mechanisms of Borreliella genomic variability including homologous gene and plasmid exchanges as well as non-homologous rearrangements.

Whole-genome sequencing of Leptospira interrogans from southern Brazil: genetic features of a highly virulent strain

BACKGROUND

Leptospirosis is the most widespread zoonotic disease. It is caused by infection with pathogenic Leptospira species, of which over 300 serovars have been described. The accurate identification of the causative Leptospira spp. is required to ascertain the pathogenic status of the local isolates.

OBJECTIVES

This study aimed to obtain the complete genome sequence of a virulent Leptospira interrogans strain isolated from southern Brazil and to describe its genetic features.

METHODS

The whole genome was sequenced by next-generation sequencing (Ion Torrent). The genome was assembled, scaffolded, annotated, and manually reviewed. Mutations were identified based on a variant calling analysis using the genome of L. interrogans strain Fiocruz L1-130 as a reference.

FINDINGS

The entire genome had an average GC content of 35%. The variant calling analysis identified 119 single nucleotide polymorphisms (SNPs), from which 30 led to a missense mutation. The structural analyses identified potential evidence of genomic inversions, translocations, and deletions in both the chromosomes.

MAIN CONCLUSIONS

The genome properties provide comprehensive information about the local isolates of Leptospira spp., and thereby, could facilitate the identification of new targets for the development of diagnostic kits and vaccines.

Diversification of Secondary Metabolite Biosynthetic Gene Clusters Coincides with Lineage Divergence in Streptomyces

We have identified Streptomyces sister-taxa which share a recent common ancestor and nearly identical small subunit (SSU) rRNA gene sequences, but inhabit distinct geographic ranges demarcated by latitude and have sufficient genomic divergence to represent distinct species. Here, we explore the evolutionary dynamics of secondary metabolite biosynthetic gene clusters (SMGCs) following lineage divergence of these sister-taxa. These sister-taxa strains contained 310 distinct SMGCs belonging to 22 different gene cluster classes. While there was broad conservation of these 22 gene cluster classes among the genomes analyzed, each individual genome harbored a different number of gene clusters within each class. A total of nine SMGCs were conserved across nearly all strains, but the majority (57%) of SMGCs were strain-specific. We show that while each individual genome has a unique combination of SMGCs, this diversity displays lineage-level modularity. Overall, the northern-derived (NDR) clade had more SMGCs than the southern-derived (SDR) clade (40.7 ± 3.9 and 33.8 ± 3.9, mean and S.D., respectively). This difference in SMGC content corresponded with differences in the number of predicted open reading frames (ORFs) per genome (7775 ± 196 and 7093 ± 205, mean and S.D., respectively) such that the ratio of SMGC:ORF did not differ between sister-taxa genomes. We show that changes in SMGC diversity between the sister-taxa were driven primarily by gene acquisition and deletion events, and these changes were associated with an overall change in genome size which accompanied lineage divergence.

ReprDB and panDB: minimalist databases with maximal microbial representation

BACKGROUND

Profiling of shotgun metagenomic samples is hindered by a lack of unified microbial reference genome databases that (i) assemble genomic information from all open access microbial genomes, (ii) have relatively small sizes, and (iii) are compatible to various metagenomic read mapping tools. Moreover, computational tools to rapidly compile and update such databases to accommodate the rapid increase in new reference genomes do not exist. As a result, database-guided analyses often fail to profile a substantial fraction of metagenomic shotgun sequencing reads from complex microbiomes.

RESULTS

We report pipelines that efficiently traverse all open access microbial genomes and assemble non-redundant genomic information. The pipelines result in two species-resolution microbial reference databases of relatively small sizes: reprDB, which assembles microbial representative or reference genomes, and panDB, for which we developed a novel iterative alignment algorithm to identify and assemble non-redundant genomic regions in multiple sequenced strains. With the databases, we managed to assign taxonomic labels and genome positions to the majority of metagenomic reads from human skin and gut microbiomes, demonstrating a significant improvement over a previous database-guided analysis on the same datasets.

CONCLUSIONS

reprDB and panDB leverage the rapid increases in the number of open access microbial genomes to more fully profile metagenomic samples. Additionally, the databases exclude redundant sequence information to avoid inflated storage or memory space and indexing or analyzing time. Finally, the novel iterative alignment algorithm significantly increases efficiency in pan-genome identification and can be useful in comparative genomic analyses.

seq-seq-pan: building a computational pan-genome data structure on whole genome alignment

Background

The increasing application of next generation sequencing technologies has led to the availability of thousands of reference genomes, often providing multiple genomes for the same or closely related species. The current approach to represent a species or a population with a single reference sequence and a set of variations cannot represent their full diversity and introduces bias towards the chosen reference. There is a need for the representation of multiple sequences in a composite way that is compatible with existing data sources for annotation and suitable for established sequence analysis methods. At the same time, this representation needs to be easily accessible and extendable to account for the constant change of available genomes.

Results

We introduce seq-seq-pan, a framework that provides methods for adding or removing new genomes from a set of aligned genomes and uses these to construct a whole genome alignment. Throughout the sequential workflow the alignment is optimized for generating a representative linear presentation of the aligned set of genomes, that enables its usage for annotation and in downstream analyses.

Conclusions

By providing dynamic updates and optimized processing, our approach enables the usage of whole genome alignment in the field of pan-genomics. In addition, the sequential workflow can be used as a fast alternative to existing whole genome aligners for aligning closely related genomes.

seq-seq-pan is freely available at https://gitlab.com/rki_bioinformatics

Electronic supplementary material

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

An Analysis of the Epidemic of Klebsiella pneumoniae Carbapenemase-Producing K. pneumoniae: Convergence of Two Evolutionary Mechanisms Creates the "Perfect Storm"

Background

Carbapenem resistance is a critical healthcare challenge worldwide. Particularly concerning is the widespread dissemination of Klebsiella pneumoniae carbapenemase (KPC). Klebsiella pneumoniae harboring blaKPC (KPC-Kpn) is endemic in many areas including the United States, where the epidemic was primarily mediated by the clonal dissemination of Kpn ST258. We postulated that the spread of blaKPC in other regions occurs by different and more complex mechanisms. To test this, we investigated the evolution and dynamics of spread of KPC-Kpn in Colombia, where KPC became rapidly endemic after emerging in 2005.

Methods

We sequenced the genomes of 133 clinical isolates recovered from 24 tertiary care hospitals located in 10 cities throughout Colombia, between 2002 (before the emergence of KPC-Kpn) and 2014. Phylogenetic reconstructions and evolutionary mapping were performed to determine temporal and genetic associations between the isolates.

Results

Our results indicate that the start of the epidemic was driven by horizontal dissemination of mobile genetic elements carrying blaKPC-2, followed by the introduction and subsequent spread of clonal group 258 (CG258) isolates containing blaKPC-3.

Conclusions

The combination of 2 evolutionary mechanisms of KPC-Kpn within a challenged health system of a developing country created the "perfect storm" for sustained endemicity of these multidrug-resistant organisms in Colombia.