Variation In Gene Content Research Articles

Time-calibrated genomic evolution of a monomorphic bacterium during its establishment as an endemic crop pathogen.

Horizontal gene transfer is of major evolutionary importance as it allows for the redistribution of phenotypically important genes among lineages. Such genes with essential functions include those involved in resistance to antimicrobial compounds and virulence factors in pathogenic bacteria. Understanding gene turnover at microevolutionary scales is critical to assess the pace of this evolutionary process. Here, we characterized and quantified gene turnover for the epidemic lineage of a bacterial plant pathogen of major agricultural importance worldwide. Relying on a dense geographic sampling spanning 39years of evolution, we estimated both the dynamics of single nucleotide polymorphism accumulation and gene content turnover. We identified extensive gene content variation among lineages even at the smallest phylogenetic and geographic scales. Gene turnover rate exceeded nucleotide substitution rate by three orders of magnitude. Accessory genes were found preferentially located on plasmids, but we identified a highly plastic chromosomal region hosting ecologically important genes such as transcription activator-like effectors. Whereas most changes in the gene content are probably transient, the rapid spread of a mobile element conferring resistance to copper compounds widely used for the management of plant bacterial pathogens illustrates how some accessory genes can become ubiquitous within a population over short timeframes.

Carriage and Gene Content Variability of the pESI-Like Plasmid Associated with Salmonella Infantis Recently Established in United States Poultry Production.

Salmonella Infantis carrying extended spectrum β-lactamase blaCTX-M-65 on a pESI-like megaplasmid has recently emerged in United States poultry. In order to determine the carriage rate and gene content variability of this plasmid in U.S. Salmonella Infantis, whole genome sequences of Salmonella isolates from humans and animals in the U.S. and internationally containing the pESI-like plasmid were analyzed. The U.S. Department of Agriculture Food Safety and Inspection Service (FSIS) identified 654 product sampling isolates containing pESI-like plasmids through hazard analysis and critical control point (HACCP) verification testing in 2017 and 2018. The Centers for Disease Control and Prevention identified 55 isolates with pESI-like plasmids in 2016–2018 through the National Antimicrobial Resistance Monitoring System. Approximately 49% of pESI-like plasmids from FSIS verification isolates and 71% from CDC NARMS contained blaCTX-M-65. Pan-plasmid genome analysis was also performed. All plasmids contained traN and more than 95% contained 172 other conserved genes; 61% contained blaCTX-M-65. In a hierarchical clustering analysis, some plasmids from U.S. animal sources clustered together and some plasmids from South America clustered together, possibly indicating multiple plasmid lineages. However, most plasmids contained similar genes regardless of origin. Carriage of the pESI-like plasmid in U.S. appears to be limited to Salmonella Infantis and carriage rates increased from 2017 to 2018.

Evolution along the parasitism-mutualism continuum determines the genetic repertoire of prophages.

Integrated into their bacterial hosts’ genomes, prophage sequences exhibit a wide diversity of length and gene content, from highly degraded cryptic sequences to intact, functional prophages that retain a full complement of lytic-function genes. We apply three approaches—bioinformatics, analytical modelling and computational simulation—to understand the diverse gene content of prophages. In the bioinformatics work, we examine the distributions of over 50,000 annotated prophage genes identified in 1384 prophage sequences, comparing the gene repertoires of intact and incomplete prophages. These data indicate that genes involved in the replication, packaging, and release of phage particles have been preferentially lost in incomplete prophages, while tail fiber, transposase and integrase genes are significantly enriched. Consistent with these results, our mathematical and computational approaches predict that genes involved in phage lytic function are preferentially lost, resulting in shorter prophages that often retain genes that benefit the host. Informed by these models, we offer novel hypotheses for the enrichment of integrase and transposase genes in cryptic prophages. Overall, we demonstrate that functional and cryptic prophages represent a diversity of genetic sequences that evolve along a parasitism-mutualism continuum.

Genomic epidemiology of nontoxigenic Corynebacterium diphtheriae from King County, Washington State, USA between July 2018 and May 2019.

Between July 2018 and May 2019, Corynebacterium diphtheriae was isolated from eight patients with non-respiratory infections, seven of whom experienced homelessness and had stayed at shelters in King County, WA, USA. All isolates were microbiologically identified as nontoxigenic C. diphtheriae biovar mitis. Whole-genome sequencing confirmed that all case isolates were genetically related, associated with sequence type 445 and differing by fewer than 24 single-nucleotide polymorphisms (SNPs). Compared to publicly available C. diphtheriae genomic data, these WA isolates formed a discrete cluster with SNP variation consistent with previously reported outbreaks. Virulence-related gene content variation within the highly related WA cluster isolates was also observed. These results indicated that genome characterization can readily support epidemiology of nontoxigenic C. diphtheriae .

Diversification of methanogens into hyperalkaline serpentinizing environments through adaptations to minimize oxidant limitation

Metagenome assembled genomes (MAGs) and single amplified genomes (SAGs) affiliated with two distinct Methanobacterium lineages were recovered from subsurface fracture waters of the Samail Ophiolite, Sultanate of Oman. Lineage Type I was abundant in waters with circumneutral pH, whereas lineage Type II was abundant in hydrogen rich, hyperalkaline waters. Type I encoded proteins to couple hydrogen oxidation to CO2 reduction, typical of hydrogenotrophic methanogens. Surprisingly, Type II, which branched from the Type I lineage, lacked homologs of two key oxidative [NiFe]-hydrogenases. These functions were presumably replaced by formate dehydrogenases that oxidize formate to yield reductant and cytoplasmic CO2 via a pathway that was unique among characterized Methanobacteria, allowing cells to overcome CO2/oxidant limitation in high pH waters. This prediction was supported by microcosm-based radiotracer experiments that showed significant biological methane generation from formate, but not bicarbonate, in waters where the Type II lineage was detected in highest relative abundance. Phylogenetic analyses and variability in gene content suggested that recent and ongoing diversification of the Type II lineage was enabled by gene transfer, loss, and transposition. These data indicate that selection imposed by CO2/oxidant availability drove recent methanogen diversification into hyperalkaline waters that are heavily impacted by serpentinization.

Complete genome sequence and analysis of Alcaligenes faecalis strain Mc250, a new potential plant bioinoculant

Here we present and analyze the complete genome of Alcaligenes faecalis strain Mc250 (Mc250), a bacterium isolated from the roots of Mimosa calodendron, an endemic plant growing in ferruginous rupestrian grasslands in Minas Gerais State, Brazil. The genome has 4,159,911 bp and 3,719 predicted protein-coding genes, in a single chromosome. Comparison of the Mc250 genome with 36 other Alcaligenes faecalis genomes revealed that there is considerable gene content variation among these strains, with the core genome representing only 39% of the protein-coding gene repertoire of Mc250. Mc250 encodes a complete denitrification pathway, a network of pathways associated with phenolic compounds degradation, and genes associated with HCN and siderophores synthesis; we also found a repertoire of genes associated with metal internalization and metabolism, sulfate/sulfonate and cysteine metabolism, oxidative stress and DNA repair. These findings reveal the genomic basis for the adaptation of this bacterium to the harsh environmental conditions from where it was isolated. Gene clusters associated with ectoine, terpene, resorcinol, and emulsan biosynthesis that can confer some competitive advantage were also found. Experimental results showed that Mc250 was able to reduce (~60%) the virulence phenotype of the plant pathogen Xanthomonas citri subsp. citri when co-inoculated in Citrus sinensis, and was able to eradicate 98% of juveniles and stabilize the hatching rate of eggs to 4% in two species of agricultural nematodes. These results reveal biotechnological potential for the Mc250 strain and warrant its further investigation as a biocontrol and plant growth-promoting bacterium.

Genomic Analysis of SXT/R391 Integrative Conjugative Elements From Proteus mirabilis Isolated in Brazil.

Integrative conjugative elements (ICEs) are widespread in many bacterial species, often carrying antibiotic resistance determinants. In the present work, we screened a collection of Proteus mirabilis clinical isolates for the presence of type 1 SXT/R391 ICEs. Among the 76 isolates analyzed, 5 of them carry such elements. The complete sequences of these elements were obtained. One of the isolates carried the CMY-2 beta-lactamase gene in a transposon and is nearly identical to the element ICEPmiJpn1 previously described in Japan, and later shown to be present in other parts of the world, indicating global spread of this element. Nevertheless, the Brazilian isolate carrying ICEPmiJpn1 is not clonally related to the other lineages carrying the same element around the world. The other ICEs identified in this work do not carry known antibiotic resistance markers and are diverse in variable gene content and size, suggesting that these elements may be responsible for the acquisition of other advantageous traits by bacteria. Some sequences carried by these elements in Brazilian strains were not previously found in other SXT/R391 variants.

Evolutionary Mechanisms of Long-Term Genome Diversification Associated With Niche Partitioning in Marine Picocyanobacteria.

Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus are the most abundant photosynthetic organisms on Earth, an ecological success thought to be linked to the differential partitioning of distinct ecotypes into specific ecological niches. However, the underlying processes that governed the diversification of these microorganisms and the appearance of niche-related phenotypic traits are just starting to be elucidated. Here, by comparing 81 genomes, including 34 new Synechococcus, we explored the evolutionary processes that shaped the genomic diversity of picocyanobacteria. Time-calibration of a core-protein tree showed that gene gain/loss occurred at an unexpectedly low rate between the different lineages, with for instance 5.6 genes gained per million years (My) for the major Synechococcus lineage (sub-cluster 5.1), among which only 0.71/My have been fixed in the long term. Gene content comparisons revealed a number of candidates involved in nutrient adaptation, a large proportion of which are located in genomic islands shared between either closely or more distantly related strains, as identified using an original network construction approach. Interestingly, strains representative of the different ecotypes co-occurring in phosphorus-depleted waters (Synechococcus clades III, WPC1, and sub-cluster 5.3) were shown to display different adaptation strategies to this limitation. In contrast, we found few genes potentially involved in adaptation to temperature when comparing cold and warm thermotypes. Indeed, comparison of core protein sequences highlighted variants specific to cold thermotypes, notably involved in carotenoid biosynthesis and the oxidative stress response, revealing that long-term adaptation to thermal niches relies on amino acid substitutions rather than on gene content variation. Altogether, this study not only deciphers the respective roles of gene gains/losses and sequence variation but also uncovers numerous gene candidates likely involved in niche partitioning of two key members of the marine phytoplankton.

Comparative analyses of saprotrophy in Salisapilia sapeloensis and diverse plant pathogenic oomycetes reveal lifestyle-specific gene expression.

ABSTRACTOomycetes include many devastating plant pathogens. Across oomycete diversity, plant-infecting lineages are interspersed by non-pathogenic ones. Unfortunately, our understanding of the evolution of lifestyle switches is hampered by a scarcity of data on the molecular biology of saprotrophic oomycetes, ecologically important primary colonizers of dead tissue that can serve as informative reference points for understanding the evolution of pathogens. Here, we established Salisapilia sapeloensis as a tractable system for the study of saprotrophic oomycetes. We generated multiple transcriptomes from S. sapeloensis and compared them with (i) 22 oomycete genomes and (ii) the transcriptomes of eight pathogenic oomycetes grown under 13 conditions. We obtained a global perspective on gene expression signatures of oomycete lifestyles. Our data reveal that oomycete saprotrophs and pathogens use similar molecular mechanisms for colonization but exhibit distinct expression patterns. We identify a S. sapeloensis-specific array and expression of carbohydrate-active enzymes and putative regulatory differences, highlighted by distinct expression levels of transcription factors. Salisapilia sapeloensis expresses only a small repertoire of candidates for virulence-associated genes. Our analyses suggest lifestyle-specific gene regulatory signatures and that, in addition to variation in gene content, shifts in gene regulatory networks underpin the evolution of oomycete lifestyles.

Multiple Intraspecific Variations of Mitochondrial Genomes in the Green-Tide Forming Alga, Ulva compressa Linnaeus (Ulvophyceae, Chlorophyta)

To gain further insights into intraspecific evolution of Ulva mitochondrial genomes, the mitogenomes of three morphotypes of the green-tide forming alga, Ulva compressa Linnaeus from China and USA, were sequenced and compared with the available data from Ulvophyceae. The U. compressa mitogenomes displayed significantly variable size at intraspecific level ranging from 61,700 to >66,587 bp, due to difference in acquisition of foreign DNA fragments, gain or loss of both group I and II introns, and non-coding intergenic spacer regions. The U. compressa mitogenomes harbored variable gene content ranging from 69 genes (including orfs) in Uco1 to 76 in Uco5, and contained different intron content from 4 introns in Uco3 and Uco4 to 7 in Uco1. A total of 63 genes and only two group IB introns (intron cox1-1107 and cox1-1125) were shared by these five mitogenomes. The U. compressa mitogenomes accumulated much more inverted repeat (IR) elements, ranging from 45 in Uco1 to 88 in Uco2, than that of the other Ulva species (3 - 34). A locally collinear block of eight genes (rps11-rps19-rps4-rpl16-trnR-trnQ-trnE-trnS) with the size of 3,631 bp has been inverted only in Uco1 indicating that the rearrangement event happened after its divergence from Uco2-5 and might be relate to a specific IR element. A majority of the common genes (≈76%) displayed high identity (>98%) among these five mitogenomes, while some low values were observed in six genes mainly due to duplication and insertion/deletion mutations of small DNA fragments. Our study presented the first case of multiple intraspecific variations in ulvophycean mitogenomes, and indicated that the mitogenome will be a valuable tool for understanding the native or nonindigenous nature of the cosmopolitan Ulva species.

Genome compartmentalization predates species divergence in the plant pathogen genus Zymoseptoria

BackgroundAntagonistic co-evolution can drive rapid adaptation in pathogens and shape genome architecture. Comparative genome analyses of several fungal pathogens revealed highly variable genomes, for many species characterized by specific repeat-rich genome compartments with exceptionally high sequence variability. Dynamic genome structure may enable fast adaptation to host genetics. The wheat pathogen Zymoseptoria tritici with its highly variable genome, has emerged as a model organism to study genome evolution of plant pathogens. Here, we compared genomes of Z. tritici isolates and of sister species infecting wild grasses to address the evolution of genome composition and structure.ResultsUsing long-read technology, we sequenced and assembled genomes of Z. ardabiliae, Z. brevis, Z. pseudotritici and Z. passerinii, together with two isolates of Z. tritici. We report a high extent of genome collinearity among Zymoseptoria species and high conservation of genomic, transcriptomic and epigenomic signatures of compartmentalization. We identify high gene content variability both within and between species. In addition, such variability is mainly limited to the accessory chromosomes and accessory compartments. Despite strong host specificity and non-overlapping host-range between species, predicted effectors are mainly shared among Zymoseptoria species, yet exhibiting a high level of presence-absence polymorphism within Z. tritici. Using in planta transcriptomic data from Z. tritici, we suggest different roles for the shared orthologs and for the accessory genes during infection of their hosts.ConclusionDespite previous reports of high genomic plasticity in Z. tritici, we describe here a high level of conservation in genomic, epigenomic and transcriptomic composition and structure across the genus Zymoseptoria. The compartmentalized genome allows the maintenance of a functional core genome co-occurring with a highly variable accessory genome.

FLAVi: An Enhanced Annotator for Viral Genomes of Flaviviridae.

Responding to the ongoing and severe public health threat of viruses of the family Flaviviridae, including dengue, hepatitis C, West Nile, yellow fever, and Zika, demands a greater understanding of how these viruses emerge and spread. Updated phylogenies are central to this understanding. Most cladograms of Flaviviridae focus on specific lineages and ignore outgroups, hampering the efficacy of the analysis to test ingroup monophyly and relationships. This is due to the lack of annotated Flaviviridae genomes, which has gene content variation among genera. This variation makes analysis without partitioning difficult. Therefore, we developed an annotation pipeline for the genera of Flaviviridae (Flavirirus, Hepacivirus, Pegivirus, and Pestivirus, named “Fast Loci Annotation of Viruses” (FLAVi; http://flavi-web.com/), that combines ab initio and homology-based strategies. FLAVi recovered 100% of the genes in Flavivirus and Hepacivirus genomes. In Pegivirus and Pestivirus, annotation efficiency was 100% except for one partition each. There were no false positives. The combined phylogenetic analysis of multiple genes made possible by annotation has clear impacts over the tree topology compared to phylogenies that we inferred without outgroups or data partitioning. The final tree is largely congruent with previous hypotheses and adds evidence supporting the close phylogenetic relationship between dengue and Zika.

Alternative strategies of nutrient acquisition and energy conservation map to the biogeography of marine ammonia-oxidizing archaea

Ammonia-oxidizing archaea (AOA) are among the most abundant and ubiquitous microorganisms in the ocean, exerting primary control on nitrification and nitrogen oxides emission. Although united by a common physiology of chemoautotrophic growth on ammonia, a corresponding high genomic and habitat variability suggests tremendous adaptive capacity. Here, we compared 44 diverse AOA genomes, 37 from species cultivated from samples collected across diverse geographic locations and seven assembled from metagenomic sequences from the mesopelagic to hadopelagic zones of the deep ocean. Comparative analysis identified seven major marine AOA genotypic groups having gene content correlated with their distinctive biogeographies. Phosphorus and ammonia availabilities as well as hydrostatic pressure were identified as selective forces driving marine AOA genotypic and gene content variability in different oceanic regions. Notably, AOA methylphosphonate biosynthetic genes span diverse oceanic provinces, reinforcing their importance for methane production in the ocean. Together, our combined comparative physiological, genomic, and metagenomic analyses provide a comprehensive view of the biogeography of globally abundant AOA and their adaptive radiation into a vast range of marine and terrestrial habitats.

A Hopeful Sea-Monster: A Very Large Homologous Recombination Event Impacting the Core Genome of the Marine Pathogen Vibrio anguillarum.

Vibrio anguillarum is the causative agent of vibriosis in many species important to aquaculture. We generated whole genome sequence (WGS) data on a diverse collection of 64 V. anguillarum strains, which we supplemented with 41 publicly available genomes to produce a combined dataset of 105 strains. These WGS data resolved six major lineages (L1-L6), and the additional use of multilocus sequence analysis (MLSA) clarified the association of L1 with serotype O1 and Salmonidae hosts (salmon/trout), and L2 with serotypes O2a/O2b/O2c and Gadidae hosts (cod). Our analysis also revealed a large-scale homologous replacement of 526-kb of core genome in an L2 strain from a con-specific donor. Although the strains affected by this recombination event are exclusively associated with Gadidae, we find no clear genetic evidence that it has played a causal role in host specialism. Whilst it is established that Vibrio species freely recombine, to our knowledge this is the first report of a contiguous recombinational replacement of this magnitude in any Vibrio genome. We also note a smaller accessory region of high single nucleotide polymorphism (SNP) density and gene content variation that contains lipopolysaccharide biosynthesis genes which may play a role in determining serotype.

Vast Differences in Strain-Level Diversity in the Gut Microbiota of Two Closely Related Honey Bee Species

SummaryMost bacterial species encompass strains with vastly different gene content. Strain diversity in microbial communities is therefore considered to be of functional importance. Yet little is known about the extent to which related microbial communities differ in diversity at this level and which underlying mechanisms may constrain and maintain strain-level diversity. Here, we used shotgun metagenomics to characterize and compare the gut microbiota of two honey bee species, Apis mellifera and Apis cerana, which diverged about 6 mya. Although the host species are colonized largely by the same bacterial 16S rRNA phylotypes, we find that their communities are host specific when analyzed with genomic resolution. Moreover, despite their similar ecology, A. mellifera displayed a much higher diversity of strains and functional gene content in the microbiota compared to A. cerana, both per colony and per individual bee. In particular, the gene repertoire for polysaccharide degradation was massively expanded in the microbiota of A. mellifera relative to A. cerana. Bee management practices, divergent ecological adaptation, or habitat size may have contributed to the observed differences in microbiota genomic diversity of these key pollinator species. Our results illustrate that the gut microbiota of closely related animal hosts can differ vastly in genomic diversity while displaying similar levels of diversity based on the 16S rRNA gene. Such differences are likely to have consequences for gut microbiota functioning and host-symbiont interactions, highlighting the need for metagenomic studies to understand the ecology and evolution of microbial communities.

Natural killer cell receptor variants and chronic hepatitis B virus infection in the Vietnamese population

ObjectivesGenes of host immunity play an important role in disease pathogenesis and are determinants of clinical courses of infections, including hepatitis B virus (HBV). Killer-cell immunoglobulin-like receptor (KIR), expressed on the surface of natural killer cells (NK), regulate NK cell cytotoxicity by interacting with human leukocyte antigen (HLA) class I molecules and are candidates for influencing the course of HBV. This study evaluated whether variations in KIR gene content and HLA-C ligands are associated with HBV and with the development of liver cirrhosis and hepatocellular carcinoma. MethodsA Vietnamese study cohort (HBV n = 511; controls n = 140) was genotyped using multiplex sequence-specific polymerase chain reaction (PCR-SSP) followed by melting curve analysis. ResultsThe presence of the functional allelic group of KIR2DS4 was associated with an increased risk of chronic HBV (OR = 1.86, pcorr = 0.02), while KIR2DL2+HLA-C1 (OR = 0.62, pcorr = 0.04) and KIR2DL3+HLA-C1 (OR = 0.48, pcorr = 0.04) were associated with a decreased risk. The pair KIR2DL3+HLA-C1 was associated with liver cirrhosis (OR = 0.40, pcorr = 0.01). The presence of five or more activating KIR variants was associated with hepatocellular carcinoma (OR = 0.53, pcorr = 0.04). ConclusionsKIR gene content variation and combinations KIR-HLA influence the outcome of HBV infection.

Selection Is a Significant Driver of Gene Gain and Loss in the Pangenome of the Bacterial Genus Sulfurovum in Geographically Distinct Deep-Sea Hydrothermal Vents.

Microbial genomes have highly variable gene content, and the evolutionary history of microbial populations is shaped by gene gain and loss mediated by horizontal gene transfer and selection. To evaluate the influence of selection on gene content variation in hydrothermal vent microbial populations, we examined 22 metagenome-assembled genomes (MAGs) (70 to 97% complete) from the ubiquitous vent Epsilonbacteraeota genus Sulfurovum that were recovered from two deep-sea hydrothermal vent regions, Axial Seamount in the northeastern Pacific Ocean (13 MAGs) and the Mid-Cayman Rise in the Caribbean Sea (9 MAGs). Genes involved in housekeeping functions were highly conserved across Sulfurovum lineages. However, genes involved in environment-specific functions, and in particular phosphate regulation, were found mostly in Sulfurovum genomes from the Mid-Cayman Rise in the low-phosphate Atlantic Ocean environment, suggesting that nutrient limitation is an important selective pressure for these bacteria. Furthermore, genes that were rare within the pangenome were more likely to undergo positive selection than genes that were highly conserved in the pangenome, and they also appeared to have experienced gene-specific sweeps. Our results suggest that selection is a significant driver of gene gain and loss for dominant microbial lineages in hydrothermal vents and highlight the importance of factors like nutrient limitation in driving microbial adaptation and evolution.IMPORTANCE Microbes can alter their gene content through the gain and loss of genes. However, there is some debate as to whether natural selection or neutral processes play a stronger role in molding the gene content of microbial genomes. In this study, we examined variation in gene content for the Epsilonbacteraeota genus Sulfurovum from deep-sea hydrothermal vents, which are dynamic habitats known for extensive horizontal gene transfer within microbial populations. Our results show that natural selection is a strong driver of Sulfurovum gene content and that nutrient limitation in particular has shaped the Sulfurovum genome, leading to differences in gene content between ocean basins. Our results also suggest that recently acquired genes undergo stronger selection than genes that were acquired in the more distant past. Overall, our results highlight the importance of natural selection in driving the evolution of microbial populations in these dynamic habitats.

A 19-isolate reference-quality global pangenome for the fungal wheat pathogen Zymoseptoria tritici

BackgroundThe gene content of a species largely governs its ecological interactions and adaptive potential. A species is therefore defined by both core genes shared between all individuals and accessory genes segregating presence-absence variation. There is growing evidence that eukaryotes, similar to bacteria, show intra-specific variability in gene content. However, it remains largely unknown how functionally relevant such a pangenome structure is for eukaryotes and what mechanisms underlie the emergence of highly polymorphic genome structures.ResultsHere, we establish a reference-quality pangenome of a fungal pathogen of wheat based on 19 complete genomes from isolates sampled across six continents. Zymoseptoria tritici causes substantial worldwide losses to wheat production due to rapidly evolved tolerance to fungicides and evasion of host resistance. We performed transcriptome-assisted annotations of each genome to construct a global pangenome. Major chromosomal rearrangements are segregating within the species and underlie extensive gene presence-absence variation. Conserved orthogroups account for only ~ 60% of the species pangenome. Investigating gene functions, we find that the accessory genome is enriched for pathogenesis-related functions and encodes genes involved in metabolite production, host tissue degradation and manipulation of the immune system. De novo transposon annotation of the 19 complete genomes shows that the highly diverse chromosomal structure is tightly associated with transposable element content. Furthermore, transposable element expansions likely underlie recent genome expansions within the species.ConclusionsTaken together, our work establishes a highly complex eukaryotic pangenome providing an unprecedented toolbox to study how pangenome structure impacts crop-pathogen interactions.

Transcriptional regulation of organohalide pollutant utilisation in bacteria

Organohalides are organic molecules formed biotically and abiotically, both naturally and through industrial production. They are usually toxic and represent a health risk for living organisms, including humans. Bacteria capable of degrading organohalides for growth express dehalogenase genes encoding enzymes that cleave carbon-halogen bonds. Such bacteria are of potential high interest for bioremediation of contaminated sites. Dehalogenase genes are often part of gene clusters that may include regulators, accessory genes and genes for transporters and other enzymes of organohalide degradation pathways. Organohalides and their degradation products affect the activity of regulatory factors, and extensive genome-wide modulation of gene expression helps dehalogenating bacteria to cope with stresses associated with dehalogenation, such as intracellular increase of halides, dehalogenase-dependent acid production, organohalide toxicity and misrouting and bottlenecks in metabolic fluxes. This review focuses on transcriptional regulation of gene clusters for dehalogenation in bacteria, as studied in laboratory experiments and in situ. The diversity in gene content, organization and regulation of such gene clusters is highlighted for representative organohalide-degrading bacteria. Selected examples illustrate a key, overlooked role of regulatory processes, often strain-specific, for efficient dehalogenation and productive growth in presence of organohalides.

A snapshot of diversity: Intraclonal variation of Escherichia coli clones as commensals and pathogens

Whole-genome sequencing has enabled detailed studies on bacterial evolution during infection, but there is limited knowledge on intraclonal variation. In this study, we sought to provide a snapshot of the intraclonal diversity of Escherichia coli as both commensal in the faecal environment and pathogen during urinary tract infection, respectively. This was performed by whole-genome sequencing and analyses of single nucleotide polymorphisms (SNPs) and gene-content variation in ten isolates belonging to the same clone and isolated from rectal swabs or urine samples. We identified only one clone in eight of the nine urines sampled (89 %). In both the commensal and pathogenic state, the within-host diversity was limited with intraclonal SNP diversity of 0–2 non-synonymous SNPs for each clone. The genetic diversity showed variation in gene content in a range of 2–15 genes in total for all clones, including genes positioned on plasmids, and in the K- and O-antigen cluster. The observed SNP- and gene variation shows that sampling of one colony would be enough for surveillance, outbreak investigations and clonal evolution. However, for studies of adaptation during or between colonization and infection, this variation is relevant to consider.