Horizontal Gene Transfer Events Research Articles

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.

The ESKAPE mobilome contributes to the spread of antimicrobial resistance and CRISPR-mediated conflict between mobile genetic elements

Mobile genetic elements (MGEs) mediate the shuffling of genes among organisms. They contribute to the spread of virulence and antibiotic resistance (AMR) genes in human pathogens, such as the particularly problematic group of ESKAPE pathogens. Here, we performed the first systematic analysis of MGEs, including plasmids, prophages, and integrative and conjugative/mobilizable elements (ICEs/IMEs), across all ESKAPE pathogens. We found that different MGE types are asymmetrically distributed across these pathogens, and that most horizontal gene transfer (HGT) events are restricted by phylum or genus. We show that the MGEs proteome is involved in diverse functional processes and distinguish widespread proteins within the ESKAPE context. Moreover, anti-CRISPRs and AMR genes are overrepresented in the ESKAPE mobilome. Our results also underscore species-specific trends shaping the number of MGEs, AMR, and virulence genes across pairs of conspecific ESKAPE genomes with and without CRISPR-Cas systems. Finally, we observed that CRISPR spacers found on prophages, ICEs/IMEs, and plasmids have different targeting biases: while plasmid and prophage CRISPRs almost exclusively target other plasmids and prophages, respectively, ICEs/IMEs CRISPRs preferentially target prophages. Overall, our study highlights the general importance of the ESKAPE mobilome in contributing to the spread of AMR and mediating conflict among MGEs.

Genome-wide signatures of adaptation to extreme environments in red algae

The high temperature, acidity, and heavy metal-rich environments associated with hot springs have a major impact on biological processes in resident cells. One group of photosynthetic eukaryotes, the Cyanidiophyceae (Rhodophyta), has successfully thrived in hot springs and associated sites worldwide for more than 1 billion years. Here, we analyze chromosome-level assemblies from three representative Cyanidiophyceae species to study environmental adaptation at the genomic level. We find that subtelomeric gene duplication of functional genes and loss of canonical eukaryotic traits played a major role in environmental adaptation, in addition to horizontal gene transfer events. Shared responses to environmental stress exist in Cyanidiales and Galdieriales, however, most of the adaptive genes (e.g., for arsenic detoxification) evolved independently in these lineages. Our results underline the power of local selection to shape eukaryotic genomes that may face vastly different stresses in adjacent, extreme microhabitats.

Horizontal gene transfer: Driving the evolution and adaptation of plants.

Horizontal gene transfer greatly contributes to the diversification and long-term evolution of green plants. Recent studies suggest that horizontal gene transfer events drove the evolution and adaptation of charophyte green algae and land plants.

Horizontal gene transfer: Implications in human health and diseases

Horizontal gene transfer (HGT) between the different organisms of the same or different species can cause both beneficial and adverse biological consequences. HGT has been demonstrated among all the phylogenetic groups including human beings. There is rising evidence that prokaryotic HGT is a dominant evolutionary force. Integration of microbes with human beings and their interaction influences genetic, physiological, and biotic sustainability of host. In the present era, this mechanism in one end is used for constructing desired genome for gene therapy, and in other hand, it is decisive in progression and surfacing many human diseases because of rapid evolution of microorganism and resistant of these microorganism to the treated antibiotics. Transformation, transduction, and conjugation are the principal HGT mechanisms in prokaryote, whereas exosomes, apoptotic bodies, cell free DNA acts as mediators in eukaryotes. Conducive environment and diversified microbiome in human gut and their interaction with host facilitate HGT in gastrointestinal microbiome, which can modulate human health and diseases. Integration of viral and bacterial genome by HGT into human genomic DNA can lead to cancers such as human papillomavirus, acute myeloid leukemia, adenocarcinoma of stomach, and many more. Many metabolic, immunological, and neurodegenerative diseases are also emerging because of HGT. The development of resistance to antibiotics or multidrug in human beings are the consequences of HGT among microbiome and simultaneously clinicians therapeutic concern. Hence, the way of understanding different mediators of HGT and their integration may facilitate better diagnostic and therapeutic modalities in control of diseases.

Robust analysis of prokaryotic pangenome gene gain and loss rates with Panstripe.

Horizontal gene transfer (HGT) plays a critical role in the evolution and diversification of many microbial species. The resulting dynamics of gene gain and loss can have important implications for the development of antibiotic resistance and the design of vaccine and drug interventions. Methods for the analysis of gene presence/absence patterns typically do not account for errors introduced in the automated annotation and clustering of gene sequences. In particular, methods adapted from ecological studies, including the pangenome gene accumulation curve, can be misleading as they may reflect the underlying diversity in the temporal sampling of genomes rather than a difference in the dynamics of HGT. Here, we introduce Panstripe, a method based on generalized linear regression that is robust to population structure, sampling bias, and errors in the predicted presence/absence of genes. We show using simulations that Panstripe can effectively identify differences in the rate and number of genes involved in HGT events, and illustrate its capability by analyzing several diverse bacterial genome data sets representing major human pathogens.