Replication-associated inversions are the dominant form of bacterial chromosome structural variation.

The structural arrangements of bacterial chromosomes vary widely between closely related species and can result in significant phenotypic outcomes. The appearance of large-scale chromosomal inversions that are symmetric relative to markers for the origin of replication (OriC) has been previously observed; however, the overall prevalence of replication-associated structural rearrangements (RASRs) in bacteria and their causal mechanisms are currently unknown. Here, we systematically identify the locations of RASRs in species with multiple complete-sequenced genomes and investigate potential mediating biological mechanisms. We found that 247 of 313 species contained sequences with at least one large (>50 Kb) inversion in their sequence comparisons, and the aggregated inversion distances away from symmetry were normally distributed with a mean of zero. Many inversions that were offset from dnaA were found to be centered on a different marker for the OriC Instances of flanking repeats provide evidence that breaks formed during the replication process could be repaired to opposing positions. We also found a strong relationship between the later stages of replication and the range in distance variation from symmetry.

The Spruce Budworm Genome: Reconstructing the Evolutionary History of Antifreeze Proteins.

Insects have developed various adaptations to survive harsh winter conditions. Among freeze-intolerant species, some produce "antifreeze proteins" (AFPs) that bind to nascent ice crystals and inhibit further ice growth. Such is the case of the spruce budworm, Choristoneura fumiferana (Lepidoptera: Tortricidae), a destructive North American conifer pest that can withstand temperatures below -30°C. Despite the potential importance of AFPs in the adaptive diversification of Choristoneura, genomic tools to explore their origins have until now been limited. Here we present a chromosome-scale genome assembly for C. fumiferana, which we used to conduct comparative genomic analyses aimed at reconstructing the evolutionary history of tortricid AFPs. The budworm genome features 16 genes homologous to previously reported C. fumiferana AFPs (CfAFPs), 15 of which map to a single region on chromosome 18. Fourteen of these were also detected in five congeneric species, indicating Choristoneura AFP diversification occurred before the speciation event that led to C. fumiferana. Although budworm AFPs were previously considered unique to the genus Choristoneura, a search for homologs targeting recently sequenced tortricid genomes identified seven CfAFP-like genes in the distantly related Notocelia uddmanniana. High structural similarity between Notocelia and Choristoneura AFPs suggests a common origin, despite the absence of homologs in three related tortricids. Interestingly, one Notocelia AFP formed the C-terminus of a "zonadhesin-like" protein, possibly representing the ancestral condition from which tortricid AFPs evolved. Future work should clarify the evolutionary path of AFPs between Notocelia and Choristoneura and assess the role of the "zonadhesin-like" protein as precursor of tortricid AFPs.

Complete Genome Sequences of Klebsiella michiganensis and Citrobacter farmeri, KPC-2-Producers Serially Isolated from a Single Patient.

Carbapenemase-producing Enterobacterales, including KPC-2 producers, have become a major clinical problem. During an outbreak in Quebec City, Canada, KPC-2-producing Klebsiella michiganensis and Citrobacter farmeri were isolated from a patient six weeks apart. We determined their complete genome sequences. Both isolates carried nearly identical IncN2 plasmids with blaKPC-2 on a Tn4401b element. Both strains also carried IncP1 plasmids, but that of C. farmeri did not carry a Beta-lactamase gene, whereas that of K. michiganensis carried a second copy of blaKPC-2 on Tn4401b. These results suggest recent plasmid transfer between the two species and a recent transposition event.

Complete Genome Sequences of Two Marine Biofilm Isolates, Leisingera sp. nov. Strains 201A and 204H, Novel Representatives of the Roseobacter Group.

Here, we report the complete-genome assemblies of biofilm isolates 201A and 204H. They possess six and seven plasmids, respectively, with a size ranging from 44 kb to 159 kb. Genomic comparisons place the two strains into one new species belonging to the genus Leisingera as novel representatives of the Roseobacter group.

De novo assembly of the olive fruit fly (Bactrocera oleae) genome with linked-reads and long-read technologies minimizes gaps and provides exceptional Y chromosome assembly.

BACKGROUND: The olive fruit fly, Bactrocera oleae, is the most important pest in the olive fruit agribusiness industry. This is because female flies lay their eggs in the unripe fruits and upon hatching the larvae feed on the fruits thus destroying them. The lack of a high-quality genome and other genomic and transcriptomic data has hindered progress in understanding the fly's biology and proposing alternative control methods to pesticide use. RESULTS: Genomic DNA was sequenced from male and female Demokritos strain flies, maintained in the laboratory for over 45 years. We used short-, mate-pair-, and long-read sequencing technologies to generate a combined male-female genome assembly (GenBank accession GCA_001188975.2). Genomic DNA sequencing from male insects using 10x Genomics linked-reads technology followed by mate-pair and long-read scaffolding and gap-closing generated a highly contiguous 489 Mb genome with a scaffold N50 of 4.69 Mb and L50 of 30 scaffolds (GenBank accession GCA_001188975.4). RNA-seq data generated from 12 tissues and/or developmental stages allowed for genome annotation. Short reads from both males and females and the chromosome quotient method enabled identification of Y-chromosome scaffolds which were extensively validated by PCR. CONCLUSIONS: The high-quality genome generated represents a critical tool in olive fruit fly research. We provide an extensive RNA-seq data set, and genome annotation, critical towards gaining an insight into the biology of the olive fruit fly. In addition, elucidation of Y-chromosome sequences will advance our understanding of the Y-chromosome's organization, function and evolution and is poised to provide avenues for sterile insect technique approaches.

Dominant gut Prevotella copri in gastrectomised non-obese diabetic Goto-Kakizaki rats improves glucose homeostasis through enhanced FXR signalling.

AIMS/HYPOTHESIS: Drug and surgical-based therapies in type 2 diabetes are associated with altered gut microbiota architecture. Here we investigated the role of the gut microbiome in improved glucose homeostasis following bariatric surgery. METHODS: We carried out gut microbiome analyses in gastrectomised (by vertical sleeve gastrectomy [VSG]) rats of the Goto-Kakizaki (GK) non-obese model of spontaneously occurring type 2 diabetes, followed by physiological studies in the GK rat. RESULTS: VSG in the GK rat led to permanent improvement of glucose tolerance associated with minor changes in the gut microbiome, mostly characterised by significant enrichment of caecal Prevotella copri. Gut microbiota enrichment with P. copri in GK rats through permissive antibiotic treatment, inoculation of gut microbiota isolated from gastrectomised GK rats, and direct inoculation of P. copri, resulted in significant improvement of glucose tolerance, independent of changes in body weight. Plasma bile acids were increased in GK rats following inoculation with P. copri and P. copri-enriched microbiota from VSG-treated rats; the inoculated GK rats then showed increased liver glycogen and upregulated expression of Fxr (also known as Nr1h4), Srebf1c, Chrebp (also known as Mlxipl) and Il10 and downregulated expression of Cyp7a1. CONCLUSIONS: Our data underline the impact of intestinal P. copri on improved glucose homeostasis through enhanced bile acid metabolism and farnesoid X receptor (FXR) signalling, which may represent a promising opportunity for novel type 2 diabetes therapeutics.

Expansion of LINEs and species-specific DNA repeats drives genome expansion in Asian Gypsy Moths.

Two subspecies of Asian gypsy moth (AGM), Lymantria dispar asiatica and L. dispar japonica, pose a serious alien invasive threat to North American forests. Despite decades of research on the ecology and biology of this pest, limited AGM-specific genomic resources are currently available. Here, we report on the genome sequences and functional content of these AGM subspecies. The genomes of L.d. asiatica and L.d. japonica are the largest lepidopteran genomes sequenced to date, totaling 921 and 999 megabases, respectively. Large genome size in these subspecies is driven by the accumulation of specific classes of repeats. Genome-wide metabolic pathway reconstructions suggest strong genomic signatures of energy-related pathways in both subspecies, dominated by metabolic functions related to thermogenesis. The genome sequences reported here will provide tools for probing the molecular mechanisms underlying phenotypic traits that are thought to enhance AGM invasiveness.

Identification and characterization of OmpT-like proteases in uropathogenic Escherichia coli clinical isolates.

Bacterial colonization of the urogenital tract is limited by innate defenses, including the production of antimicrobial peptides (AMPs). Uropathogenic Escherichia coli (UPEC) resist AMP-killing to cause a range of urinary tract infections (UTIs) including asymptomatic bacteriuria, cystitis, pyelonephritis, and sepsis. UPEC strains have high genomic diversity and encode numerous virulence factors that differentiate them from non-UTI-causing strains, including ompT. As OmpT homologs cleave and inactivate AMPs, we hypothesized that UPEC strains from patients with symptomatic UTIs have high OmpT protease activity. Therefore, we measured OmpT activity in 58 clinical E. coli isolates. While heterogeneous OmpT activities were observed, OmpT activity was significantly greater in UPEC strains isolated from patients with symptomatic infections. Unexpectedly, UPEC strains exhibiting the greatest protease activities harbored an additional ompT-like gene called arlC (ompTp). The presence of two OmpT-like proteases in some UPEC isolates led us to compare the substrate specificities of OmpT-like proteases found in E. coli. While all three cleaved AMPs, cleavage efficiency varied on the basis of AMP size and secondary structure. Our findings suggest the presence of ArlC and OmpT in the same UPEC isolate may confer a fitness advantage by expanding the range of target substrates.

Altered microbiome composition in individuals with fibromyalgia.

Fibromyalgia (FM) is a prevalent syndrome, characterised by chronic widespread pain, fatigue, and impaired sleep, that is challenging to diagnose and difficult to treat. The microbiomes of 77 women with FM and that of 79 control participants were compared using 16S rRNA gene amplification and whole-genome sequencing. When comparing FM patients with unrelated controls using differential abundance analysis, significant differences were revealed in several bacterial taxa. Variance in the composition of the microbiomes was explained by FM-related variables more than by any other innate or environmental variable and correlated with clinical indices of FM. In line with observed alteration in butyrate-metabolising species, targeted serum metabolite analysis verified differences in the serum levels of butyrate and propionate in FM patients. Using machine-learning algorithms, the microbiome composition alone allowed for the classification of patients and controls (receiver operating characteristic area under the curve 87.8%). To the best of our knowledge, this is the first demonstration of gut microbiome alteration in nonvisceral pain. This observation paves the way for further studies, elucidating the pathophysiology of FM, developing diagnostic aids and possibly allowing for new treatment modalities to be explored.

A complete Leishmania donovani reference genome identifies novel genetic variations associated with virulence.

Leishmania donovani is responsible for visceral leishmaniasis, a neglected and lethal parasitic disease with limited treatment options and no vaccine. The study of L. donovani has been hindered by the lack of a high-quality reference genome and this can impact experimental outcomes including the identification of virulence genes, drug targets and vaccine development. We therefore generated a complete genome assembly by deep sequencing using a combination of second generation (Illumina) and third generation (PacBio) sequencing technologies. Compared to the current L. donovani assembly, the genome assembly reported within resulted in the closure over 2,000 gaps, the extension of several chromosomes up to telomeric repeats and the re-annotation of close to 15% of protein coding genes and the annotation of hundreds of non-coding RNA genes. It was possible to correctly assemble the highly repetitive A2 and Amastin virulence gene clusters. A comparative sequence analysis using the improved reference genome confirmed 70 published and identified 15 novel genomic differences between closely related visceral and atypical cutaneous disease-causing L. donovani strains providing a more complete map of genes associated with virulence and visceral organ tropism. Bioinformatic tools including protein variation effect analyzer and basic local alignment search tool were used to prioritize a list of potential virulence genes based on mutation severity, gene conservation and function. This complete genome assembly and novel information on virulence factors will support the identification of new drug targets and the development of a vaccine for L. donovani.

Publisher Correction: Ancient hybridization and strong adaptation to viruses across African vervet monkey populations.

In the version of this article published, in the Online Methods eight citations to supplementary material refer to the wrong supplementary items. See the correction notice for full details.

Canada-Wide Epidemic of emm74 Group A Streptococcus Invasive Disease.

BACKGROUND: The number of invasive group A Streptococcus (iGAS) infections due to hitherto extremely rare type emm74 strains has increased in several Canadian provinces since late 2015. We hypothesized that the cases recorded in the different provinces are linked and caused by strains of an emm74 clone that recently emerged and expanded explosively. METHODS: We analyzed both active and passive surveillance data for iGAS infections and used whole-genome sequencing to investigate the phylogenetic relationships of the emm74 strains responsible for these invasive infections country-wide. RESULTS: Genome analysis showed that highly clonal emm74 strains, genetically different from emm74 organisms previously circulating in Canada, were responsible for a country-wide epidemic of >160 invasive disease cases. The emerging clone belonged to multilocus sequence typing ST120. The analysis also revealed dissemination patterns of emm74 subclonal lineages across Canadian provinces. Clinical data analysis indicated that the emm74 epidemic disproportionally affected middle-aged or older male individuals. Homelessness, alcohol abuse, and intravenous drug usage were significantly associated with invasive emm74 infections. CONCLUSIONS: In a period of 20 months, an emm74 GAS clone emerged and rapidly spread across several Canadian provinces located more than 4500 km apart, causing invasive infections primarily among disadvantaged persons.

Ancient hybridization and strong adaptation to viruses across African vervet monkey populations.

Vervet monkeys are among the most widely distributed nonhuman primates, show considerable phenotypic diversity, and have long been an important biomedical model for a variety of human diseases and in vaccine research. Using whole-genome sequencing data from 163 vervets sampled from across Africa and the Caribbean, we find high diversity within and between taxa and clear evidence that taxonomic divergence was reticulate rather than following a simple branching pattern. A scan for diversifying selection across taxa identifies strong and highly polygenic selection signals affecting viral processes. Furthermore, selection scores are elevated in genes whose human orthologs interact with HIV and in genes that show a response to experimental simian immunodeficiency virus (SIV) infection in vervet monkeys but not in rhesus macaques, suggesting that part of the signal reflects taxon-specific adaptation to SIV.

Genetic variation and gene expression across multiple tissues and developmental stages in a nonhuman primate.

By analyzing multitissue gene expression and genome-wide genetic variation data in samples from a vervet monkey pedigree, we generated a transcriptome resource and produced the first catalog of expression quantitative trait loci (eQTLs) in a nonhuman primate model. This catalog contains more genome-wide significant eQTLs per sample than comparable human resources and identifies sex- and age-related expression patterns. Findings include a master regulatory locus that likely has a role in immune function and a locus regulating hippocampal long noncoding RNAs (lncRNAs), whose expression correlates with hippocampal volume. This resource will facilitate genetic investigation of quantitative traits, including brain and behavioral phenotypes relevant to neuropsychiatric disorders.

Comparative analysis of mitochondrial genomes of geographic variants of the gypsy moth, Lymantria dispar, reveals a previously undescribed genotypic entity.

The gypsy moth, Lymantria dispar L., is one of the most destructive forest pests in the world. While the subspecies established in North America is the European gypsy moth (L. dispar dispar), whose females are flightless, the two Asian subspecies, L. dispar asiatica and L. dispar japonica, have flight-capable females, enhancing their invasiveness and warranting precautionary measures to prevent their permanent establishment in North America. Various molecular tools have been developed to help distinguish European from Asian subspecies, several of which are based on the mitochondrial barcode region. In an effort to identify additional informative markers, we undertook the sequencing and analysis of the mitogenomes of 10 geographic variants of L. dispar, including two or more variants of each subspecies, plus the closely related L. umbrosa as outgroup. Several regions of the gypsy moth mitogenomes displayed nucleotide substitutions with potential usefulness for the identification of subspecies and/or geographic origins. Interestingly, the mitogenome of one geographic variant displayed significant divergence relative to the remaining variants, raising questions about its taxonomic status. Phylogenetic analyses placed this population from northern Iran as basal to the L. dispar clades. The present findings will help improve diagnostic tests aimed at limiting risks of AGM invasions.

A Syst-OMICS Approach to Ensuring Food Safety and Reducing the Economic Burden of Salmonellosis.

The Salmonella Syst-OMICS consortium is sequencing 4,500 Salmonella genomes and building an analysis pipeline for the study of Salmonella genome evolution, antibiotic resistance and virulence genes. Metadata, including phenotypic as well as genomic data, for isolates of the collection are provided through the Salmonella Foodborne Syst-OMICS database (SalFoS), at https://salfos.ibis.ulaval.ca/. Here, we present our strategy and the analysis of the first 3,377 genomes. Our data will be used to draw potential links between strains found in fresh produce, humans, animals and the environment. The ultimate goals are to understand how Salmonella evolves over time, improve the accuracy of diagnostic methods, develop control methods in the field, and identify prognostic markers for evidence-based decisions in epidemiology and surveillance.

Complete Genome of a Panresistant Pseudomonas aeruginosa Strain, Isolated from a Patient with Respiratory Failure in a Canadian Community Hospital.

We report here the complete genome sequence of a panresistant Pseudomonas aeruginosa strain, isolated from a patient with respiratory failure in Canada. No carbapenemase genes were identified. Carbapenem resistance is attributable to a frameshift in the oprD gene; the basis for colistin resistance remains undetermined.

Genome and Plasmid Analysis of blaIMP-4-Carrying Citrobacter freundii B38.

Sequencing of the blaIMP-4-carrying C. freundii B38 using the PacBio SMRT technique revealed that the genome contained a chromosome of 5,134,500 bp and three plasmids, pOZ172 (127,005 bp), pOZ181 (277,592 bp), and pOZ182 (18,467 bp). Plasmid pOZ172 was identified as IncFIIY, like pP10164-NDM and pNDM-EcGN174. It carries a class 1 integron with four cassettes (blaIMP-4-qacG2-aacA4-aphA15) and a complete hybrid tni module (tniR-tniQ-tniB-tniA). The recombination of tniR from Tn402 (identical) with tniQBA from Tn5053 (99%) occurred within the res site of Tn402/5053 The Tn402/5053-like integron, named Tn6017, was inserted into Tn1722 at the res II site. The replication, partitioning, and transfer systems of pOZ181 were similar to those of IncHI2 plasmids (e.g., R478) and contained a sul1-type class 1 integron with the cassette array orf-dfrA1-orf-gcu37-aadA5 linked to an upstream Tn1696 tnpA-tnpR and to a downstream 3' conserved sequence (3'-CS) and ISCR1 A Tn2 transposon encoding a blaTEM-1 ß-lactamase was identified on pOZ182. Other interesting resistance determinants encoded on the B38 chromosome included multidrug resistance (MDR) efflux pumps, an AmpC ß-lactamase, and resistances to Cu, Ag, As, and Zn. This is the first report of a complete tni module linked to a blaIMP-4-carrying class 1 integron, which, together with other recently reported non-sul1 integrons, represents the emergence of a distinct evolutionary lineage of class 1 integrons lacking a 3'-CS (qacEΔ1-sul1). The unique cassette array, complete tni module of Tn6017, and incompatibility group of pOZ172 suggest a blaIMP-4 evolutionary pathway in C. freundii B38 different from that for other blaIMP-4 genes found in Gram-negative bacteria in the Western Pacific region.

Genomic Recombination Leading to Decreased Virulence of Group B Streptococcus in a Mouse Model of Adult Invasive Disease.

Adult invasive disease caused by Group B Streptococcus (GBS) is increasing worldwide. Whole-genome sequencing (WGS) now permits rapid identification of recombination events, a phenomenon that occurs frequently in GBS. Using WGS, we described that strain NGBS375, a capsular serotype V GBS isolate of sequence type (ST)297, has an ST1 genomic background but has acquired approximately 300 kbp of genetic material likely from an ST17 strain. Here, we examined the virulence of this strain in an in vivo model of GBS adult invasive infection. The mosaic ST297 strain showed intermediate virulence, causing significantly less systemic infection and reduced mortality than a more virulent, serotype V ST1 isolate. Bacteremia induced by the ST297 strain was similar to that induced by a serotype III ST17 strain, which was the least virulent under the conditions tested. Yet, under normalized bacteremia levels, the in vivo intrinsic capacity to induce the production of pro-inflammatory cytokines was similar between the ST297 strain and the virulent ST1 strain. Thus, the diminished virulence of the mosaic strain may be due to reduced capacity to disseminate or multiply in blood during a systemic infection which could be mediated by regulatory factors contained in the recombined region.

Envelope-specific B-cell populations in African green monkeys chronically infected with simian immunodeficiency virus.

African green monkeys (AGMs) are natural primate hosts of simian immunodeficiency virus (SIV). Interestingly, features of the envelope-specific antibody responses in SIV-infected AGMs are distinct from that of HIV-infected humans and SIV-infected rhesus monkeys, including gp120-focused responses and rapid development of autologous neutralization. Yet, the lack of genetic tools to evaluate B-cell lineages hinders potential use of this unique non-human primate model for HIV vaccine development. Here we define features of the AGM Ig loci and compare the proportion of Env-specific memory B-cell populations to that of HIV-infected humans and SIV-infected rhesus monkeys. AGMs appear to have a higher proportion of Env-specific memory B cells that are mainly gp120 directed. Furthermore, AGM gp120-specific monoclonal antibodies display robust antibody-dependent cellular cytotoxicity and CD4-dependent virion capture activity. Our results support the use of AGMs to model induction of functional gp120-specific antibodies by HIV vaccine strategies.