Horizontal Gene Transfer Events Research Articles

A novel aminotransferase gene and its regulator acquired in Saccharomyces by a horizontal gene transfer event

BackgroundHorizontal gene transfer (HGT) is an evolutionary mechanism of adaptive importance, which has been deeply studied in wine S. cerevisiae strains, where those acquired genes conferred improved traits related to both transport and metabolism of the nutrients present in the grape must. However, little is known about HGT events that occurred in wild Saccharomyces yeasts and how they determine their phenotypes.ResultsThrough a comparative genomic approach among Saccharomyces species, we detected a subtelomeric segment present in the S. uvarum, S. kudriavzevii, and S. eubayanus species, belonging to the first species to diverge in the Saccharomyces genus, but absent in the other Saccharomyces species. The segment contains three genes, two of which were characterized, named DGD1 and DGD2. DGD1 encodes dialkylglicine decarboxylase, whose specific substrate is the non-proteinogenic amino acid 2-aminoisobutyric acid (AIB), a rare amino acid present in some antimicrobial peptides of fungal origin. DGD2 encodes putative zinc finger transcription factor, which is essential to induce the AIB-dependent expression of DGD1. Phylogenetic analysis showed that DGD1 and DGD2 are closely related to two adjacent genes present in Zygosaccharomyces.ConclusionsThe presented results show evidence of an early HGT event conferring new traits to the ancestor of the Saccharomyces genus that could be lost in the evolutionary more recent Saccharomyces species, perhaps due to loss of function during the colonization of new habitats.

GH20 and GH84 β-N-acetylglucosaminidases with different linkage specificities underpin mucin O-glycan breakdown capability of Bifidobacterium bifidum

Intestinal mucous layers mediate symbiosis and dysbiosis of host-microbe interactions. These interactions are influenced by the mucin O-glycan degrading ability of several gut microbes. The identities and prevalence of many glycoside hydrolases (GHs) involved in microbial mucin O-glycan breakdown have been previously reported; however, the exact mechanisms and extent to which these GHs are dedicated to mucin O-glycan degradation pathways warrant further research. Here, using Bifidobacterium bifidum as a model mucinolytic bacterium, we revealed that two β-N-acetylglucosaminidases belonging to the GH20 (BbhI) and GH84 (BbhIV) families play important roles in mucin O-glycan degradation. Using substrate specificity analysis of natural oligosaccharides and O-glycomic analysis of porcine gastric mucin (PGM) incubated with purified enzymes or B.bifidum carrying bbhI and/or bbhIV mutations, we showed that BbhI and BbhIV are highly specific for β-(1→3)- and β-(1→6)-GlcNAc linkages of mucin core structures, respectively. Interestingly, we found that efficient hydrolysis of the β-(1→3)-linkage by BbhI of the mucin core 4 structure [GlcNAcβ1-3(GlcNAcβ1-6)GalNAcα-O-Thr] required prior removal of the β-(1→6)-GlcNAc linkage by BbhIV. Consistent with this, inactivation of bbhIV markedly decreased the ability of B.bifidum to release GlcNAc from PGM. When combined with a bbhI mutation, we observed that the growth of the strain on PGM was reduced. Finally, phylogenetic analysis suggests that GH84 members may have gained diversified functions through microbe-microbe and host-microbe horizontal gene transfer events. Taken together, these data strongly suggest the involvement of GH84 family members in host glycan breakdown.

Investigating Possible Interspecies Communication of Plasmids Associated with Transfer of Third-Generation Cephalosporin, Quinolone, and Colistin Resistance Between Simultaneously Isolated Escherichia Coli and Klebsiella Pneumoniae.

The coinfection process producing multiple species of pathogens provides a specific ecological niche for the exchange of genetic materials between pathogens, in which plasmids play a vital role in horizontal gene transfer, especially for drug resistance, but the underlying transfer pathway remains unclear. Interspecies communication of the plasmids associated with the transfer of third-generation cephalosporins, quinolones, and colistin resistance has been observed in simultaneously isolated Escherichia coli and Klebsiella pneumoniae from abdominal drainage following surgery. The MICs of antimicrobial agents were determined by the broth microdilution method. The complete chromosome and plasmid sequences were obtained by combining Illumina paired-end short reads and MinION long reads. S1-PFGE, southern blot analysis and conjugation assay confirmed the transferability of the mcr-1-harboring plasmid. Both the E. coli isolate EC15255 and K. pneumoniae isolate KP15255 from the same specimen presented multidrug resistance. Each of them harbored one chromosome and three plasmids, and two plasmids and their mediated resistance could be transferred to the recipient by conjugation. Comparison of their genome sequences suggested that several genetic communication events occurred between species, especially among their plasmids, such as whole-plasmid transfer, insertion, deletion, amplification, or inversion. Exchange of plasmids or the genetic elements they harbor plays a critical role in antimicrobial resistance gene transmission and poses a substantial threat to nosocomial infection control, necessitating the continued surveillance of multidrug resistant pathogens, especially during coinfection. IMPORTANCE The genome sequence of bacterial pathogens commonly provides a detailed clue of genetic communication among clones or even distinct species. The intestinal microecological environment is a representative ecological niche for genetic communication. However, it is still difficult to describe the details of horizontal gene transfer or other genetic events within them because the evidence in the genome sequence is incomplete and limited. In this study, the simultaneously isolated Escherichia coli and Klebsiella pneumoniae from a coinfection process provided an excellent example for observation of interspecies communication between the two genomes and the plasmids they harbor. A complete genome sequence acquired by combining the Illumina and MinION sequencing platforms facilitated the understanding of genetic communication events, such as whole-plasmid transfer, insertion, deletion, amplification, or inversion, which contribute to antimicrobial resistance gene transmission and are a substantial threat to nosocomial infection control.

Metagenome-based metabolic modelling predicts unique microbial interactions in deep-sea hydrothermal plume microbiomes

Deep-sea hydrothermal vents are abundant on the ocean floor and play important roles in ocean biogeochemistry. In vent ecosystems such as hydrothermal plumes, microorganisms rely on reduced chemicals and gases in hydrothermal fluids to fuel primary production and form diverse and complex microbial communities. However, microbial interactions that drive these complex microbiomes remain poorly understood. Here, we use microbiomes from the Guaymas Basin hydrothermal system in the Pacific Ocean to shed more light on the key species in these communities and their interactions. We built metabolic models from metagenomically assembled genomes (MAGs) and infer possible metabolic exchanges and horizontal gene transfer (HGT) events within the community. We highlight possible archaea–archaea and archaea–bacteria interactions and their contributions to the robustness of the community. Cellobiose, D-Mannose 1-phosphate, O2, CO2, and H2S were among the most exchanged metabolites. These interactions enhanced the metabolic capabilities of the community by exchange of metabolites that cannot be produced by any other community member. Archaea from the DPANN group stood out as key microbes, benefiting significantly as acceptors in the community. Overall, our study provides key insights into the microbial interactions that drive community structure and organisation in complex hydrothermal plume microbiomes.

Metagenomic analysis reveals the short-term influences on conjugation of blaNDM-1 and microbiome in hospital wastewater by silver nanoparticles at environmental-related concentration

Hospital wastewater contains large amounts of antibiotic-resistant bacteria and serves as an important reservoir for horizontal gene transfer (HGT). However, the response of the microbiome in hospital wastewater to silver remains unclear. In this study, the short-term impacts of silver on the microbiome in hospital wastewater were investigated by metagenome next-generation sequencing. The influence of silver on the conjugation of plasmid carrying blaNDM-1 was further examined. Our results showed that in hospital wastewater, high abundances of antibiotic resistance genes (ARGs) were detected. The distribution tendencies of certain ARG types on chromosomes or plasmids were different. Clinically important ARGs were identified in phage-like contigs, indicating potential transmission via transduction. Pseudomonadales, Enterobacterales, and Bacteroidales were the major ARG hosts. Mobile genetic elements were mainly detected in plasmids and associated with various types of ARGs. The binning approach identified 29 bins that were assigned to three phyla. Various ARGs and virulence factors were identified in 14 and 11 bins, respectively. MetaCHIP identified 49 HGT events. The transferred genes were annotated as ARGs, mobile genetic elements, and functional genes, and they mainly originated from donors belonging to Bacteroides and Pseudomonadales. In addition, 20 nm AgNPs reduced microbial diversity and enhanced the relative abundance of Acinetobacter. The changes induced by 20 nm AgNPs included increases in the abundances of ARGs and genes involved lipid metabolism pathway. Conjugation experiments showed that Ag+ and 20 nm AgNPs caused 2.38-, 3.31-, 4.72-, and 4.57-fold and 1.46-, 1.61-, 3.86-, and 2.16-fold increases in conjugation frequencies of plasmid with blaNDM-1 at 0.1, 1, 10, and 100 μg/L, respectively. Our findings provide insight into the response of the microbiome in hospital wastewater to silver, emphasize the adaptation capability of Acinetobacter inhabiting hospitals against adverse environments, and highlight the promotion of silver for antibiotic resistance.

Vibrio type III secretion system 2 is not restricted to the Vibrionaceae and encodes differentially distributed repertoires of effector proteins.

Vibrio parahaemolyticus is the leading cause of seafood-borne gastroenteritis worldwide. A distinctive feature of the O3:K6 pandemic clone, and its derivatives, is the presence of a second, phylogenetically distinct, type III secretion system (T3SS2) encoded within the genomic island VPaI-7. The T3SS2 allows the delivery of effector proteins directly into the cytosol of infected eukaryotic cells to subvert key host-cell processes, critical for V. parahaemolyticus to colonize and cause disease. Furthermore, the T3SS2 also increases the environmental fitness of V. parahaemolyticus in its interaction with bacterivorous protists; hence, it has been proposed that it contributed to the global oceanic spread of the pandemic clone. Several reports have identified T3SS2-related genes in Vibrio and non-Vibrio species, suggesting that the T3SS2 gene cluster is not restricted to the Vibrionaceae and can mobilize through horizontal gene transfer events. In this work, we performed a large-scale genomic analysis to determine the phylogenetic distribution of the T3SS2 gene cluster and its repertoire of effector proteins. We identified putative T3SS2 gene clusters in 1130 bacterial genomes from 8 bacterial genera, 5 bacterial families and 47 bacterial species. A hierarchical clustering analysis allowed us to define six T3SS2 subgroups (I-VI) with different repertoires of effector proteins, redefining the concepts of T3SS2 core and accessory effector proteins. Finally, we identified a subset of the T3SS2 gene clusters (subgroup VI) that lacks most T3SS2 effector proteins described to date and provided a list of 10 novel effector candidates for this subgroup through bioinformatic analysis. Collectively, our findings indicate that the T3SS2 extends beyond the family Vibrionaceae and suggest that different effector protein repertories could have a differential impact on the pathogenic potential and environmental fitness of each bacterium that has acquired the Vibrio T3SS2 gene cluster.

Spontaneously Produced Lysogenic Phages Are an Important Component of the Soybean Bradyrhizobium Mobilome.

The ability of Bradyrhizobium spp. to nodulate and fix atmospheric nitrogen in soybean root nodules is critical to meeting humanity's nutritional needs. The intricacies of soybean bradyrhizobia-plant interactions have been studied extensively; however, bradyrhizobial ecology as influenced by phages has received somewhat less attention, even though these interactions may significantly impact soybean yield. In batch culture, four soybean bradyrhizobia strains, Bradyrhizobium japonicum S06B (S06B-Bj), B. japonicum S10J (S10J-Bj), Bradyrhizobium diazoefficiens USDA 122 (USDA 122-Bd), and Bradyrhizobium elkanii USDA 76T (USDA 76-Be), spontaneously (without apparent exogenous chemical or physical induction) produced tailed phages throughout the growth cycle; for three strains, phage concentrations exceeded cell numbers by ~3-fold after 48 h of incubation. Phage terminase large-subunit protein phylogeny revealed possible differences in phage packaging and replication mechanisms. Bioinformatic analyses predicted multiple prophage regions within each soybean bradyrhizobia genome, preventing accurate identification of spontaneously produced prophage (SPP) genomes. A DNA sequencing and mapping approach accurately delineated the boundaries of four SPP genomes within three of the soybean bradyrhizobia chromosomes and suggested that the SPPs were capable of transduction. In addition to the phages, S06B-Bj and USDA 76-Be contained three to four times more insertion sequences (IS) and large, conjugable, broad host range plasmids, both of which are known drivers of horizontal gene transfer (HGT) in soybean bradyrhizobia. These factors indicate that SPP along with IS and plasmids participate in HGT, drive bradyrhizobia evolution, and play an outsized role in bradyrhizobia ecology. IMPORTANCE Previous studies have shown that IS and plasmids mediate HGT of symbiotic nodulation (nod) genes in soybean bradyrhizobia; however, these events require close cell-to-cell contact, which could be limited in soil environments. Bacteriophage-assisted gene transduction through spontaneously produced prophages provides a stable means of HGT not limited by the constraints of proximal cell-to-cell contact. These phage-mediated HGT events may shape soybean bradyrhizobia population ecology, with concomitant impacts on soybean agriculture.

Transfer of the Integrative and Conjugative Element ICESt3 of Streptococcus thermophilus in Physiological Conditions Mimicking the Human Digestive Ecosystem.

Metagenome analyses of the human microbiome suggest that horizontal gene transfer (HGT) is frequent in these rich and complex microbial communities. However, so far, only a few HGT studies have been conducted in vivo. In this work, three different systems mimicking the physiological conditions encountered in the human digestive tract were tested, including (i) the TNO gastro-Intestinal tract Model 1 (TIM-1) system (for the upper part of the intestine), (ii) the ARtificial COLon (ARCOL) system (to mimic the colon), and (iii) a mouse model. To increase the likelihood of transfer by conjugation of the integrative and conjugative element studied in the artificial digestive systems, bacteria were entrapped in alginate, agar, and chitosan beads before being placed in the different gut compartments. The number of transconjugants detected decreased, while the complexity of the ecosystem increased (many clones in TIM-1 but only one clone in ARCOL). No clone was obtained in a natural digestive environment (germfree mouse model). In the human gut, the richness and diversity of the bacterial community would offer more opportunities for HGT events to occur. In addition, several factors (SOS-inducing agents, microbiota-derived factors) that potentially increase in vivo HGT efficiency were not tested here. Even if HGT events are rare, expansion of the transconjugant clones can happen if ecological success is fostered by selecting conditions or by events that destabilize the microbial community. IMPORTANCE The human gut microbiota plays a key role in maintaining normal host physiology and health, but its homeostasis is fragile. During their transit in the gastrointestinal tract, bacteria conveyed by food can exchange genes with resident bacteria. New traits acquired by HGT (e.g., new catabolic properties, bacteriocins, antibiotic resistance) can impact the gut microbial composition and metabolic potential. We showed here that TIM-1, a system mimicking the upper digestive tract, is a useful tool to evaluate HGT events in conditions closer to the physiological ones. Another important fact pointed out in this work is that Enterococcus faecalis is a good candidate for foreign gene acquisition. Due to its high ability to colonize the gut and acquire mobile genetic elements, this commensal bacterium could serve as an intermediate for HGT in the human gut.

Genomic characterization and identification of virulence-related genes in Vibrio nigripulchritudo isolated from white leg shrimp Penaeus vannamei.

Vibrio nigripulchritudo causes vibriosis in penaeid shrimps. Here, we used Illumina and Nanopore sequencing technologies to sequence the genomes of three of its strains (TUMSAT-V. nig1, TUMSAT-V. nig2, and TUMSAT-V. nig3) to explore opportunities for disease management. Putative virulence factors and mobile genetic elements were detected while evaluating the phylogenetic relationship of each isolated strain. The genomes consisted of two circular chromosomes (I and II) plus one or two plasmids. The size of chromosome I ranged from 4.02 to 4.07 Mb with an average GC content of 46%, while the number of predicted CDSs ranged from 3563 to 3644. The size of chromosome II ranged from 2.16 to 2.18 Mb, with an average GC content of 45.5%, and the number of predicted CDSs ranged from 1970 to 1987. Numerous virulence genes were identified related to adherence, antiphagocytosis, chemotaxis, motility, iron uptake, quorum sensing, secretion systems, and toxins in all three genomes. Higher numbers of prophages and genomic islands found in TUMSAT-V. nig1 suggest that the strain has experienced numerous horizontal gene transfer events. The presence of antimicrobial resistance genes suggests that the strains have multidrug resistance. Comparative genomic analysis showed that all three strains belonged to the same clade.

Euglenozoan kleptoplasty illuminates the early evolution of photoendosymbiosis

Kleptoplasts (kP) are distinct among photosynthetic organelles in eukaryotes (i.e., plastids) because they are routinely sequestered from prey algal cells and function only temporarily in the new host cell. Therefore, the hosts of kleptoplasts benefit from photosynthesis without constitutive photoendosymbiosis. Here, we report that the euglenozoan Rapaza viridis has only kleptoplasts derived from a specific strain of green alga, Tetraselmis sp., but no canonical plastids like those found in its sister group, the Euglenophyceae. R. viridis showed a dynamic change in the accumulation of cytosolic polysaccharides in response to light-dark cycles, and 13C isotopic labeling of ambient bicarbonate demonstrated that these polysaccharides originate insitu via photosynthesis; these data indicate that the kleptoplasts of R. viridis are functionally active. We also identified 276 sequences encoding putative plastid-targeting proteins and 35 sequences of presumed kleptoplast transporters in the transcriptome of R. viridis. These genes originated in a wide range of algae other than Tetraselmis sp., the source of the kleptoplasts, suggesting a long history of repeated horizontal gene transfer events from different algal prey cells. Many of the kleptoplast proteins, as well as the protein-targeting system, in R. viridis were shared with members of the Euglenophyceae, providing evidence that the early evolutionary stages in the green alga-derived secondary plastids of euglenophytes also involved kleptoplasty.

Palidis: fast discovery of novel insertion sequences.

The diversity of microbial insertion sequences, crucial mobile genetic elements in generating diversity in microbial genomes, needs to be better represented in current microbial databases. Identification of these sequences in microbiome communities presents some significant problems that have led to their underrepresentation. Here, we present a bioinformatics pipeline called Palidis that recognizes insertion sequences in metagenomic sequence data rapidly by identifying inverted terminal repeat regions from mixed microbial community genomes. Applying Palidis to 264 human metagenomes identifies 879 unique insertion sequences, with 519 being novel and not previously characterized. Querying this catalogue against a large database of isolate genomes reveals evidence of horizontal gene transfer events across bacterial classes. We will continue to apply this tool more widely, building the Insertion Sequence Catalogue, a valuable resource for researchers wishing to query their microbial genomes for insertion sequences.

De novo Assembly and Comparative Analyses of Mitochondrial Genomes in Piperales.

The mitochondrial genome of Liriodendron tulipifera exhibits many ancestral angiosperm features and a remarkably slow evolutionary rate, while mitochondrial genomes of other magnoliids remain yet to be characterized. We assembled nine new mitochondrial genomes, representing all genera of perianth-bearing Piperales, as well as for a member of the sister clade: three complete or nearly complete mitochondrial genomes from Aristolochiaceae and six additional draft assemblies including Thottea, Asaraceae, Lactoridaceae, and Hydnoraceae. For comparative purpose, a complete mitochondrial genome was assembled for Saururus, a member of the perianth-less Piperales. The average number of short repeats (50-99 bp) was much larger in genus Aristolochia than in other angiosperm mitochondrial genomes, and approximately 30% of repeats (<350 bp) were found to have the capacity to mediate recombination. We found mitochondrial genomes in perianth-bearing Piperales comprising conserved repertories of protein-coding genes and rRNAs but variable copy numbers of tRNA genes. We identified several shifts from cis- to trans-splicing of the Group II introns of nad1i728, cox2i373, and nad7i209. Two short regions of the cox1 and atp8 genes were likely derived from independent horizontal gene transfer events in perianth-bearing Piperales. We found biased enrichment of specific substitution types in different lineages of magnoliids and the Aristolochiaceae family showed the highest ratio of A:T > T:A substitutions of all other investigated angiosperm groups. Our study reports the first mitochondrial genomes for Piperales and uses this new information for a better understanding of the evolutionary patterns of magnoliids and angiosperms in general.

Ancient Rapid Radiation Explains Most Conflicts Among Gene Trees and Well-Supported Phylogenomic Trees of Nostocalean Cyanobacteria.

Prokaryotic genomes are often considered to be mosaics of genes that do not necessarily share the same evolutionary history due to widespread horizontal gene transfers (HGTs). Consequently, representing evolutionary relationships of prokaryotes as bifurcating trees has long been controversial. However, studies reporting conflicts among gene trees derived from phylogenomic data sets have shown that these conflicts can be the result of artifacts or evolutionary processes other than HGT, such as incomplete lineage sorting, low phylogenetic signal, and systematic errors due to substitution model misspecification. Here, we present the results of an extensive exploration of phylogenetic conflicts in the cyanobacterial order Nostocales, for which previous studies have inferred strongly supported conflicting relationships when using different concatenated phylogenomic data sets. We found that most of these conflicts are concentrated in deep clusters of short internodes of the Nostocales phylogeny, where the great majority of individual genes have low resolving power. We then inferred phylogenetic networks to detect HGT events while also accounting for incomplete lineage sorting. Our results indicate that most conflicts among gene trees are likely due to incomplete lineage sorting linked to an ancient rapid radiation, rather than to HGTs. Moreover, the short internodes of this radiation fit the expectations of the anomaly zone, i.e., a region of the tree parameter space where a species tree is discordant with its most likely gene tree. We demonstrated that concatenation of different sets of loci can recover up to 17 distinct and well-supported relationships within the putative anomaly zone of Nostocales, corresponding to the observed conflicts among well-supported trees based on concatenated data sets from previous studies. Our findings highlight the important role of rapid radiations as a potential cause of strongly conflicting phylogenetic relationships when using phylogenomic data sets of bacteria. We propose that polytomies may be the most appropriate phylogenetic representation of these rapid radiations that are part of anomaly zones, especially when all possible genomic markers have been considered to infer these phylogenies. [Anomaly zone; bacteria; horizontal gene transfer; incomplete lineage sorting; Nostocales; phylogenomic conflict; rapid radiation; Rhizonema.].

Coal-source acid mine drainage reduced the soil multidrug-dominated antibiotic resistome but increased the heavy metal(loid) resistome and energy production-related metabolism

A recent global scale study found that mining-impacted environments have multi-antibiotic resistance gene (ARG)-dominated resistomes with an abundance similar to urban sewage but much higher than freshwater sediment. These findings raised concern that mining may increase the risk of ARG environmental proliferation. The current study assessed how typical multimetal(loid)-enriched coal-source acid mine drainage (AMD) contamination affects soil resistomes by comparing with background soils unaffected by AMD. Both contaminated and background soils have multidrug-dominated antibiotic resistomes attributed to the acidic environment. AMD-contaminated soils had a lower relative abundance of ARGs (47.45 ± 23.34 ×/Gb) than background soils (85.47 ± 19.71 ×/Gb) but held high-level heavy metal(loid) resistance genes (MRGs, 133.29 ± 29.36 ×/Gb) and transposase- and insertion sequence-dominated mobile genetic elements (MGEs, 188.51 ± 21.81 ×/Gb), which was 56.26 % and 412.12 % higher than background soils, respectively.Procrustes analysis showed that the microbial community and MGEs exerted more influence on driving heavy metal(loid) resistome variation than antibiotic resistome. The microbial community increased energy production-related metabolism to fulfill the increasing energy needs required by acid and heavy metal(loid) resistance. Horizontal gene transfer (HGT) events primarily exchanged energy- and information-related genes to adapt to the harsh AMD environment. These findings provide new insight into the risk of ARG proliferation in mining environments.

Genomics discovery of giant fungal viruses from subsurface oceanic crustal fluids

The oceanic igneous crust is a vast reservoir for microbial life, dominated by diverse and active bacteria, archaea, and fungi. Archaeal and bacterial viruses were previously detected in oceanic crustal fluids at the Juan de Fuca Ridge (JdFR). Here we report the discovery of two eukaryotic Nucleocytoviricota genomes from the same crustal fluids by sorting and sequencing single virions. Both genomes have a tRNATyr gene with an intron (20 bps) at the canonical position between nucleotide 37 and 38, a common feature in eukaryotic and archaeal tRNA genes with short introns (<100 bps), and fungal genes acquired through horizontal gene transfer (HGT) events. The dominance of Ascomycota fungi as the main eukaryotes in crustal fluids and the evidence for HGT point to these fungi as the putative hosts, making these the first putative fungi-Nucleocytoviricota specific association. Our study suggests active host-viral dynamics for the only eukaryotic group found in the subsurface oceanic crust and raises important questions about the impact of viral infection on the productivity and biogeochemical cycling in this ecosystem.

The Evolutionary History of a DNA Methylase Reveals Frequent Horizontal Transfer and Within-Gene Recombination.

Inteins, often referred to as protein introns, are highly mobile genetic elements that invade conserved genes throughout the tree of life. Inteins have been found to invade a wide variety of key genes within actinophages. While in the process of conducting a survey of these inteins in actinophages, we discovered that one protein family of methylases contained a putative intein, and two other unique insertion elements. These methylases are known to occur commonly in phages as orphan methylases (possibly as a form of resistance to restriction-modification systems). We found that the methylase family is not conserved within phage clusters and has a disparate distribution across divergent phage groups. We determined that two of the three insertion elements have a patchy distribution within the methylase protein family. Additionally, we found that the third insertion element is likely a second homing endonuclease, and that all three elements (the intein, the homing endonuclease, and what we refer to as the ShiLan domain) have different insertion sites that are conserved in the methylase gene family. Furthermore, we find strong evidence that both the intein and ShiLan domain are partaking in long-distance horizontal gene transfer events between divergent methylases in disparate phage hosts within the already dispersed methylase distribution. The reticulate evolutionary history of methylases and their insertion elements reveals high rates of gene transfer and within-gene recombination in actinophages.

The Aphelenchoides genomes reveal substantial horizontal gene transfers in the last common ancestor of free-living and major plant-parasitic nematodes.

Aphelenchoides besseyi is a plant-parasitic nematode (PPN) in the family Aphelenchoididae capable of infecting more than 200 plant species. A.besseyi is also a species complex with strains exhibiting varying pathogenicity to plants. We present the genome and annotations of six Aphelenchoides species, four of which belonged to the A.besseyi species complex. Most Aphelenchoides genomes have a size of 44.7-47.4Mb and are among the smallest in clade IV, with the exception of A.fujianensis, which has a size of 143.8Mb and is one of the largest. Phylogenomic analysis successfully delimited the species complex into A.oryzae and A.pseudobesseyi and revealed a reduction of transposon elements in the last common ancestor of Aphelenchoides. Synteny analyses between reference genomes indicated that three chromosomes in A.besseyi were derived from fission and fusion events. A systematic identification of horizontal gene transfer (HGT) genes across 27 representative nematodes allowed us to identify two major episodes of acquisition corresponding to the last common ancestor of clade IV or major PPNs, respectively. These genes were mostly lost and differentially retained between clades or strains. Most HGT events were acquired from bacteria, followed by fungi, and also from plants; plant HGT was especially prevalent in Bursaphelenchus mucronatus. Our results comprehensively improve the understanding of HGT in nematodes.

DNA uptake from a laboratory environment drives unexpected adaptation of a thermophile to a minor medium component

DNA uptake is widespread among microorganisms and considered a strategy for rapid adaptation to new conditions. While both DNA uptake and adaptation are referred to in the context of natural environments, they are often studied in laboratories under defined conditions. For example, a strain of the thermophile Thermoanaerobacter kivui had been adapted to growth on high concentrations of carbon monoxide (CO). Unusual phenotypes of the CO-adapted strain prompted us to examine it more closely, revealing a horizontal gene transfer (HGT) event from another thermophile, Thermoanaerobacter sp. strain X514, being cultured in the same laboratory. The transferred genes conferred on T. kivui the ability to utilize trehalose, a trace component of the yeast-extract added to the media during CO-adaptation. This same HGT event simultaneously deleted a native operon for thiamine biosynthesis, which likely explains why the CO-adapted strain grows poorly without added vitamins. Attempts to replicate this HGT by providing T. kivui with genomic DNA from Thermoanaerobacter sp. strain X514 revealed that it is easily reproducible in the lab. This subtle form of “genome contamination” is difficult to detect, since the genome remains predominantly T. kivui, and no living cells from the original contamination remain. Unexpected HGT between two microorganisms as well as simultaneous adaptation to several conditions may occur often and unrecognized in laboratory environments, requiring caution and careful monitoring of phenotype and genotype of microorganisms that are naturally-competent for DNA uptake.

Evolution and function analysis of glycerol kinase GlpK in Pseudomonasaeruginosa

Pseudomonas aeruginosa is a Gram-negative bacterium capable of widespread niches, which is also one of the main bacteria that cause patient infection. The metabolic diversity of Pseudomonas aeruginosa is an essential factor in adapting to a variety of environments. Based on the previous studies, adaptive genetic variation in the glycerol kinase GlpK, the glycerol 3-phosphotransferase, contributes to the fitness of bacteria in human bodies, such as Mycobacterium tuberculosis and Escherichia coli. Thus, this study aimed to explore the molecular evolution and function of glpK in P. aeruginosa. Using extensive population genomic data, we have identified the prevalence of two glpK copies in P. aeruginosa that clustered into distinct branches, which were later known as Clade 1 and 2. The evolution analysis revealed that glpK in Clade 1 derived from an ancestral P. aeruginosa species and the other from an ancient horizontal gene transfer event. In addition, we confirmed that the GlpK in Clade 2 still retained glycerol kinase activity but was much weaker than that of GlpK in Clade 1. We demonstrated the importance of the critical amino acid Q70 in GlpK glycerol kinase activity by point mutation. Furthermore, Co-expression network analysis implied that the two glpK copies of P. aeruginosa regulate separate networks and may be a strategy to improve fitness in P. aeruginosa.

Horizontally acquired fungal killer protein genes affect cell development in mosses.

The moss Physcomitrium patens is crucial for studying plant development and evolution. Although the P. patens genome includes genes acquired from bacteria, fungi and viruses, the functions and evolutionary significance of these acquired genes remain largely unclear. Killer protein 4 (KP4) is a toxin secreted by the phytopathogenic fungus Ustilago maydis that inhibits the growth of sensitive target strains by blocking their calcium uptake. Here, we show that KP4 genes in mosses were acquired from fungi through at least three independent events of horizontal gene transfer. Two paralogous copies of KP4 (PpKP4-1 and PpKP4-2) exist in P. patens. Knockout mutants ppkp4-1 and ppkp4-2 showed cell death at the protonemal stage, and ppkp4-2 also exhibited defects in tip growth. We provide experimental evidence indicating that PpKP4-1/2 affects P. patens protonemal cell development by mediating cytoplasmic calcium and that KP4 genes are functionally conserved between P. patens and fungi. The present study provides additional insights into the role of horizontal gene transfer in land plant development and evolution.