Identification of horizontally transferred genes in the genus Colletotrichum reveals a steady tempo of bacterial to fungal gene transfer

Background

Horizontal gene transfer (HGT) is the stable transmission of genetic material between organisms by means other than vertical inheritance. HGT has an important role in the evolution of prokaryotes but is relatively rare in eukaryotes. HGT has been shown to contribute to virulence in eukaryotic pathogens. We studied the importance of HGT in plant pathogenic fungi by identifying horizontally transferred genes in the genomes of three members of the genus Colletotrichum.

Results

We identified eleven HGT events from bacteria into members of the genus Colletotrichum or their ancestors. The HGT events include genes involved in amino acid, lipid and sugar metabolism as well as lytic enzymes. Additionally, the putative minimal dates of transference were calculated using a time calibrated phylogenetic tree. This analysis reveals a constant flux of genes from bacteria to fungi throughout the evolution of subphylum Pezizomycotina.

Conclusions

Genes that are typically transferred by HGT are those that are constantly subject to gene duplication and gene loss. The functions of some of these genes suggest roles in niche adaptation and virulence. We found no evidence of a burst of HGT events coinciding with major geological events. In contrast, HGT appears to be a constant, albeit rare phenomenon in the Pezizomycotina, occurring at a steady rate during their evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-16-2) contains supplementary material, which is available to authorized users.     

Nitrile hydratase genes are present in multiple eukaryotic supergroups

Background

Nitrile hydratases are enzymes involved in the conversion of nitrile-containing compounds into ammonia and organic acids. Although they are widespread in prokaryotes, nitrile hydratases have only been reported in two eukaryotes: the choanoflagellate Monosiga brevicollis and the stramenopile Aureococcus anophagefferens. The nitrile hydratase gene in M. brevicollis was believed to have arisen by lateral gene transfer from a prokaryote, and is a fusion of beta and alpha nitrile hydratase subunits. Only the alpha subunit has been reported in A. anophagefferens.

Methodology/Principal Findings

Here we report the detection of nitrile hydratase genes in five eukaryotic supergroups: opisthokonts, amoebozoa, archaeplastids, CCTH and SAR. Beta-alpha subunit fusion genes are found in the choanoflagellates, ichthyosporeans, apusozoans, haptophytes, rhizarians and stramenopiles, and potentially also in the amoebozoans. An individual alpha subunit is found in a dinoflagellate and an individual beta subunit is found in a haptophyte. Phylogenetic analyses recover a clade of eukaryotic-type nitrile hydratases in the Opisthokonta, Amoebozoa, SAR and CCTH; this is supported by analyses of introns and gene architecture. Two nitrile hydratase sequences from an animal and a plant resolve in the prokaryotic nitrile hydratase clade.

Conclusions/Significance

The evidence presented here demonstrates that nitrile hydratase genes are present in multiple eukaryotic supergroups, suggesting that a subunit fusion gene was present in the last common ancestor of all eukaryotes. The absence of nitrile hydratase from several sequenced species indicates that subunits were lost in multiple eukaryotic taxa. The presence of nitrile hydratases in many other eukaryotic groups is unresolved due to insufficient data and taxon sampling. The retention and expression of the gene in distantly related eukaryotic species suggests that it plays an important metabolic role. The novel family of eukaryotic nitrile hydratases presented in this paper represents a promising candidate for research into their molecular biology and possible biotechnological applications.     

Horizontal gene transfer in an acid mine drainage microbial community

Background

Horizontal gene transfer (HGT) has been widely identified in complete prokaryotic genomes. However, the roles of HGT among members of a microbial community and in evolution remain largely unknown. With the emergence of metagenomics, it is nontrivial to investigate such horizontal flow of genetic materials among members in a microbial community from the natural environment. Because of the lack of suitable methods for metagenomics gene transfer detection, microorganisms from a low-complexity community acid mine drainage (AMD) with near-complete genomes were used to detect possible gene transfer events and suggest the biological significance.

Results

Using the annotation of coding regions by the current tools, a phylogenetic approach, and an approximately unbiased test, we found that HGTs in AMD organisms are not rare, and we predicted 119 putative transferred genes. Among them, 14 HGT events were determined to be transfer events among the AMD members. Further analysis of the 14 transferred genes revealed that the HGT events affected the functional evolution of archaea or bacteria in AMD, and it probably shaped the community structure, such as the dominance of G-plasma in archaea in AMD through HGT.

Conclusions

Our study provides a novel insight into HGT events among microorganisms in natural communities. The interconnectedness between HGT and community evolution is essential to understand microbial community formation and development.

Electronic supplementary material

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

Reference genome of wild goat (capra aegagrus) and sequencing of goat breeds provide insight into genic basis of goat domestication

Background

Domestic goats (Capra hircus) have been selected to play an essential role in agricultural production systems, since being domesticated from their wild progenitor, bezoar (Capra aegagrus). A detailed understanding of the genetic consequences imparted by the domestication process remains a key goal of evolutionary genomics.

Results

We constructed the reference genome of bezoar and sequenced representative breeds of domestic goats to search for genomic changes that likely have accompanied goat domestication and breed formation. Thirteen copy number variation genes associated with coat color were identified in domestic goats, among which ASIP gene duplication contributes to the generation of light coat-color phenotype in domestic goats. Analysis of rapidly evolving genes identified genic changes underlying behavior-related traits, immune response and production-related traits.

Conclusion

Based on the comparison studies of copy number variation genes and rapidly evolving genes between wild and domestic goat, our findings and methodology shed light on the genetic mechanism of animal domestication and will facilitate future goat breeding.

Electronic supplementary material

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

Horizontal gene transfers link a human MRSA pathogen to contagious bovine mastitis bacteria

Background

Acquisition of virulence factors and antibiotic resistance by many clinically important bacteria can be traced to horizontal gene transfer (HGT) between related or evolutionarily distant microflora. Comparative genomic analysis has become an important tool for identifying HGT DNA in emerging pathogens. We have adapted the multi-genome alignment tool EvoPrinter to facilitate discovery of HGT DNA sequences within bacterial genomes and within their mobile genetic elements.

Principal Findings

EvoPrinter analysis of 13 different Staphylococcus aureus genomes revealed that one of the human isolates, the hospital epidemic methicillin-resistant MRSA252 strain, uniquely shares multiple putative HGT DNA sequences with different causative agents of bovine mastitis that are not found in the other human S. aureus isolates. MRSA252 shares over 14 different DNA sequence blocks with the bovine mastitis ET3 S. aureus strain RF122, and many of the HGT DNAs encode virulence factors. EvoPrinter analysis of the MRSA252 chromosome also uncovered virulence-factor encoding HGT events with the genome of Listeria monocytogenes and a Staphylococcus saprophyticus associated plasmid. Both bacteria are also causal agents of contagious bovine mastitis.

Conclusions

EvoPrinter analysis reveals that the human MRSA252 strain uniquely shares multiple DNA sequence blocks with different causative agents of bovine mastitis, suggesting that HGT events may be occurring between these pathogens. These findings have important implications with regard to animal husbandry practices that inadvertently enhance the contact of human and livestock bacterial pathogens.     

Ancient origin of animal U-box ubiquitin ligases

Background

The patterns of emergence and diversification of the families of ubiquitin ligases provide insights about the evolution of the eukaryotic ubiquitination system. U-box ubiquitin ligases (UULs) are proteins characterized by containing a peculiar protein domain known as U box. In this study, the origin of the animal UUL genes is described.

Results

Phylogenetic and structural data indicate that six of the seven main UUL-encoding genes found in humans (UBE4A, UBE4B, UIP5, PRP19, CHIP and CYC4) were already present in the ancestor of all current metazoans and the seventh (WDSUB1) is found in placozoans, cnidarians and bilaterians. The fact that only 4 - 5 genes orthologous to the human ones are present in the choanoflagellate Monosiga brevicollis suggests that several animal-specific cooptions of the U box to generate new genes occurred. Significantly, Monosiga contains five additional UUL genes that are not present in animals. One of them is also present in distantly-related protozoans. Along animal evolution, losses of UUL-encoding genes are rare, except in nematodes, which lack three of them. These general patterns are highly congruent with those found for other two families (RBR, HECT) of ubiquitin ligases.

Conclusions

Finding that the patterns of emergence, diversification and loss of three unrelated families of ubiquitin ligases (RBR, HECT and U-box) are parallel indicates that there are underlying, linage-specific evolutionary forces shaping the complexity of the animal ubiquitin system.

     

Origin of Saxitoxin Biosynthetic Genes in Cyanobacteria

Background

Paralytic shellfish poisoning (PSP) is a potentially fatal syndrome associated with the consumption of shellfish that have accumulated saxitoxin (STX). STX is produced by microscopic marine dinoflagellate algae. Little is known about the origin and spread of saxitoxin genes in these under-studied eukaryotes. Fortuitously, some freshwater cyanobacteria also produce STX, providing an ideal model for studying its biosynthesis. Here we focus on saxitoxin-producing cyanobacteria and their non-toxic sisters to elucidate the origin of genes involved in the putative STX biosynthetic pathway.

Methodology/Principal Findings

We generated a draft genome assembly of the saxitoxin-producing (STX+) cyanobacterium Anabaena circinalis ACBU02 and searched for 26 candidate saxitoxin­genes (named sxtA to sxtZ) that were recently identified in the toxic strain Cylindrospermopsis raciborskii T3. We also generated a draft assembly of the non-toxic (STX−) sister Anabaena circinalis ACFR02 to aid the identification of saxitoxin-specific genes. Comparative phylogenomic analyses revealed that nine putative STX genes were horizontally transferred from non-cyanobacterial sources, whereas one key gene (sxtA) originated in STX+ cyanobacteria via two independent horizontal transfers followed by fusion. In total, of the 26 candidate saxitoxin-genes, 13 are of cyanobacterial provenance and are monophyletic among the STX+ taxa, four are shared amongst STX+ and STX-cyanobacteria, and the remaining nine genes are specific to STX+ cyanobacteria.

Conclusions/Significance

Our results provide evidence that the assembly of STX genes in ACBU02 involved multiple HGT events from different sources followed presumably by coordination of the expression of foreign and native genes in the common ancestor of STX+ cyanobacteria. The ability to produce saxitoxin was subsequently lost multiple independent times resulting in a nested relationship of STX+ and STX− strains among Anabaena circinalis strains.     

HGTector: an automated method facilitating genome-wide discovery of putative horizontal gene transfers

Background

First pass methods based on BLAST match are commonly used as an initial step to separate the different phylogenetic histories of genes in microbial genomes, and target putative horizontal gene transfer (HGT) events. This will continue to be necessary given the rapid growth of genomic data and the technical difficulties in conducting large-scale explicit phylogenetic analyses. However, these methods often produce misleading results due to their inability to resolve indirect phylogenetic links and their vulnerability to stochastic events.

Results

A new computational method of rapid, exhaustive and genome-wide detection of HGT was developed, featuring the systematic analysis of BLAST hit distribution patterns in the context of a priori defined hierarchical evolutionary categories. Genes that fall beyond a series of statistically determined thresholds are identified as not adhering to the typical vertical history of the organisms in question, but instead having a putative horizontal origin. Tests on simulated genomic data suggest that this approach effectively targets atypically distributed genes that are highly likely to be HGT-derived, and exhibits robust performance compared to conventional BLAST-based approaches. This method was further tested on real genomic datasets, including Rickettsia genomes, and was compared to previous studies. Results show consistency with currently employed categories of HGT prediction methods. In-depth analysis of both simulated and real genomic data suggests that the method is notably insensitive to stochastic events such as gene loss, rate variation and database error, which are common challenges to the current methodology. An automated pipeline was created to implement this approach and was made publicly available at: https://github.com/DittmarLab/HGTector. The program is versatile, easily deployed, has a low requirement for computational resources.

Conclusions

HGTector is an effective tool for initial or standalone large-scale discovery of candidate HGT-derived genes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-717) contains supplementary material, which is available to authorized users.     

Genome diversification within a clonal population of pandemic Vibrio parahaemolyticus seems to depend on the life circumstances of each individual bacteria

Background

New strains of Vibrio parahaemolyticus that cause diarrhea in humans by seafood ingestion periodically emerge through continuous evolution in the ocean. Influx and expansion in the Southern Chilean ocean of a highly clonal V. parahaemolyticus (serotype O3:K6) population from South East Asia caused one of the largest seafood-related diarrhea outbreaks in the world. Here, genomics analyses of isolates from this rapidly expanding clonal population offered an opportunity to observe the molecular evolutionary changes often obscured in more diverse populations.

Results

Whole genome sequence comparison of eight independent isolates of this population from mussels or clinical cases (from different years) was performed. Differences of 1366 to 217,729 bp genome length and 13 to 164 bp single nucleotide variants (SNVs) were found. Most genomic differences corresponded to the presence of regions unique to only one or two isolates, and were probably acquired by horizontal gene transfer (HGT). Some DNA gain was chromosomal but most was in plasmids. One isolate had a large region (8,644 bp) missing, which was probably caused by excision of a prophage. Genome innovation by the presence of unique DNA, attributable to HGT from related bacteria, varied greatly among the isolates, with values of 1,366 (ten times the number of highest number of SNVs) to 217,729 (a thousand times more than the number of highest number of SNVs).

Conclusions

The evolutionary forces (SNVs, HGT) acting on each isolate of the same population were found to differ to an extent that probably depended on the ecological scenario and life circumstances of each bacterium.

Electronic supplementary material

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

Extended Spectrum Beta-Lactamase-Producing Enterobacteriaceae in International Travelers and Non-Travelers in New York City

Background

We performed this study 1) to determine the prevalence of community-associated extended spectrum beta-lactamase producing Enterobacteriaceae (ESBLPE) colonization and infection in New York City (NYC); 2) to determine the prevalence of newly-acquired ESBLPE during travel; 3) to look for similarilties in contemporaneous hospital-associated bloodstream ESBLPE and travel-associated ESBLPE.

Methods

Subjects were recruited from a travel medicine practice and consented to submit pre- and post-travel stools, which were assessed for the presence of ESBLPE. Pre-travel stools and stools submitted for culture were used to estimate the prevalence of community-associated ESBLPE. The prevalence of ESBLPE-associated urinary tract infections was calculated from available retrospective data. Hospital-associated ESBLPE were acquired from saved bloodstream isolates. All ESBLPE underwent multilocus sequence typing (MLST) and ESBL characterization.

Results

One of 60 (1.7%) pre- or non-travel associated stool was colonized with ESBLPE. Among community-associated urine specimens, 1.3% of Escherichia coli and 1.4% of Klebsiella pneumoniae were identified as ESBLPE. Seven of 28 travelers (25.0%) acquired a new ESBLPE during travel. No similarities were found between travel-associated ESBLPE and hospital-associated ESBLPE. A range of imported ESBL genes were found, including CTX-M-14 and CTX-15.

Conclusion

ESBL colonization and infection were relatively low during the study period in NYC. A signficant minority of travelers acquired new ESBLPE during travel.     

Towards a better detection of horizontally transferred genes by combining unusual properties effectively

Background

Horizontal gene transfer (HGT) is one of the major mechanisms contributing to microbial genome diversification. A number of computational methods for finding horizontally transferred genes have been proposed in the past decades; however none of them has provided a reliable detector yet. In existing parametric approaches, only one single compositional property can participate in the detection process, or the results obtained through each single property are just simply combined. It’s known that different properties may mean different information, so the single property can’t sufficiently contain the information encoded by gene sequences. In addition, the class imbalance problem in the datasets, which also results in great errors for the gene detection, hasn’t been considered by the published methods. Here we developed an effective classifier system (Hgtident) that used support vector machine (SVM) by combining unusual properties effectively for HGT detection.

Results

Our approach Hgtident includes the introduction of more representative datasets, optimization of SVM model, feature selection, handling of imbalance problem in the datasets and extensive performance evaluation via systematic cross-validation methods. Through feature selection, we found that JS-DN and JS-CB have higher discriminating power for HGT detection, while GC1–GC3 and k-mer (k = 1, 2, …, 7) make the least contribution. Extensive experiments indicated the new classifier could reduce Mean error dramatically, and also improve Recall by a certain level. For the testing genomes, compared with the existing popular multiple-threshold approach, on average, our Recall and Mean error was respectively improved by 2.81% and reduced by 26.32%, which means that numerous false positives were identified correctly.

Conclusions

Hgtident introduced here is an effective approach for better detecting HGT. Combining multiple features of HGT is also essential for a wider range of HGT events detection.     

Functional biogeography as evidence of gene transfer in hypersaline microbial communities

Background

Horizontal gene transfer (HGT) plays a major role in speciation and evolution of bacteria and archaea by controlling gene distribution within an environment. However, information that links HGT to a natural community using relevant population-genetics parameters and spatial considerations is scarce. The Great Salt Lake (Utah, USA) provides an excellent model for studying HGT in the context of biogeography because it is a contiguous system with dispersal limitations due to a strong selective salinity gradient. We hypothesize that in spite of the barrier to phylogenetic dispersal, functional characteristics—in the form of HGT—expand beyond phylogenetic limitations due to selective pressure.

Methodology and Results

To assay the functional genes and microorganisms throughout the GSL, we used a 16S rRNA oligonucleotide microarray (Phylochip) and a functional gene array (GeoChip) to measure biogeographic patterns of nine microbial communities. We found a significant difference in biogeography based on microarray analyses when comparing Sørensen similarity values for presence/absence of function and phylogeny (Student''s t-test; p = 0.005).

Conclusion and Significance

Biogeographic patterns exhibit behavior associated with horizontal gene transfer in that informational genes (16S rRNA) have a lower similarity than functional genes, and functional similarity is positively correlated with lake-wide selective pressure. Specifically, high concentrations of chromium throughout GSL correspond to an average similarity of chromium resistance genes that is 22% higher than taxonomic similarity. This suggests active HGT may be measured at the population level in microbial communities and these biogeographic patterns may serve as a model to study bacteria adaptation and speciation.     

Genetic diversity of selected genes that are potentially economically important in feral sheep of New Zealand

Background

Feral sheep are considered to be a source of genetic variation that has been lost from their domestic counterparts through selection.

Methods

This study investigates variation in the genes KRTAP1-1, KRT33, ADRB3 and DQA2 in Merino-like feral sheep populations from New Zealand and its offshore islands. These genes have previously been shown to influence wool, lamb survival and animal health.

Results

All the genes were polymorphic, but no new allele was identified in the feral populations. In some of these populations, allele frequencies differed from those observed in commercial Merino sheep and other breeds found in New Zealand. Heterozygosity levels were comparable to those observed in other studies on feral sheep. Our results suggest that some of the feral populations may have been either inbred or outbred over the duration of their apparent isolation.

Conclusion

The variation described here allows us to draw some conclusions about the likely genetic origin of the populations and selective pressures that may have acted upon them, but they do not appear to be a source of new genetic material, at least for these four genes.     

Expression analysis of asthma candidate genes during human and murine lung development

Background

Little is known about the role of most asthma susceptibility genes during human lung development. Genetic determinants for normal lung development are not only important early in life, but also for later lung function.

Objective

To investigate the role of expression patterns of well-defined asthma susceptibility genes during human and murine lung development. We hypothesized that genes influencing normal airways development would be over-represented by genes associated with asthma.

Methods

Asthma genes were first identified via comprehensive search of the current literature. Next, we analyzed their expression patterns in the developing human lung during the pseudoglandular (gestational age, 7-16 weeks) and canalicular (17-26 weeks) stages of development, and in the complete developing lung time series of 3 mouse strains: A/J, SW, C57BL6.

Results

In total, 96 genes with association to asthma in at least two human populations were identified in the literature. Overall, there was no significant over-representation of the asthma genes among genes differentially expressed during lung development, although trends were seen in the human (Odds ratio, OR 1.22, confidence interval, CI 0.90-1.62) and C57BL6 mouse (OR 1.41, CI 0.92-2.11) data. However, differential expression of some asthma genes was consistent in both developing human and murine lung, e.g. NOD1, EDN1, CCL5, RORA and HLA-G. Among the asthma genes identified in genome wide association studies, ROBO1, RORA, HLA-DQB1, IL2RB and PDE10A were differentially expressed during human lung development.

Conclusions

Our data provide insight about the role of asthma susceptibility genes during lung development and suggest common mechanisms underlying lung morphogenesis and pathogenesis of respiratory diseases.     

Agrobacterium-mediated transformation of tomato elicits unexpected flower phenotypes with similar gene expression profiles

Background

Genetic transformation mediated by Agrobacterium tumefaciens is known to cause unexpected phenotypes. Mutations of a specific set of homeotic genes can result in alterred floral structure.

Methodology/Principal Findings

Previously we identified two genes (LeTGA1 and SOLly GLB1) induced by nutrient availability in tomato. To further elucidate their function, we sought to knock out the genes using antisense RNAi. When antisense constructs for the two different tomato genes were each transformed into Micro-Tina tomato plants, one primary transformant with similar mutant flower phenotypes was identified from transformation of each construct. Microarray analysis shows that a similar set of genes were up- or downregulated in both mutants. Sequencing of insertion sites indicates that each is inserted into a repetitive region which could impact expression of affected genes but direct alteration of floral homeotic gene sequences was not detected.

Conclusion

This is the first report that dominant flower mutations could be caused by genetic transformation designed to knock out two nutrient stress related genes.     

Detection of Prokaryotic Genes in the Amphimedon queenslandica Genome

Horizontal gene transfer (HGT) is common between prokaryotes and phagotrophic eukaryotes. In metazoans, the scale and significance of HGT remains largely unexplored but is usually linked to a close association with parasites and endosymbionts. Marine sponges (Porifera), which host many microorganisms in their tissues and lack an isolated germ line, are potential carriers of genes transferred from prokaryotes. In this study, we identified a number of potential horizontally transferred genes within the genome of the sponge, Amphimedon queenslandica. We further identified homologs of some of these genes in other sponges. The transferred genes, most of which possess catalytic activity for carbohydrate or protein metabolism, have assimilated host genome characteristics and are actively expressed. The diversity of functions contributed by the horizontally transferred genes is likely an important factor in the adaptation and evolution of A. queenslandica. These findings highlight the potential importance of HGT on the success of sponges in diverse ecological niches.     

Genome sequence of Erinnyis ello granulovirus (ErelGV), a natural cassava hornworm pesticide and the first sequenced sphingid-infecting betabaculovirus

Background

Cassava (Manihot esculenta) is the basic source for dietary energy of 500 million people in the world. In Brazil, Erinnyis ello ello (Lepidoptera: Sphingidae) is a major pest of cassava crops and a bottleneck for its production. In the 1980s, a naturally occurring baculovirus was isolated from E. ello larva and successfully applied as a bio-pesticide in the field. Here, we described the structure, the complete genome sequence, and the phylogenetic relationships of the first sphingid-infecting betabaculovirus.

Results

The baculovirus isolated from the cassava hornworm cadavers is a betabaculovirus designated Erinnyis ello granulovirus (ErelGV). The 102,759 bp long genome has a G + C content of 38.7%. We found 130 putative ORFs coding for polypeptides of at least 50 amino acid residues. Only eight genes were found to be unique. ErelGV is closely related to ChocGV and PiraGV isolates. We did not find typical homologous regions and cathepsin and chitinase homologous genes are lacked. The presence of he65 and p43 homologous genes suggests horizontal gene transfer from Alphabaculovirus. Moreover, we found a nucleotide metabolism-related gene and two genes that could be acquired probably from Densovirus.

Conclusions

The ErelGV represents a new virus species from the genus Betabaculovirus and is the closest relative of ChocGV. It contains a dUTPase-like, a he65-like, p43-like genes, which are also found in several other alpha- and betabaculovirus genomes, and two Densovirus-related genes. Importantly, recombination events between insect viruses from unrelated families and genera might drive baculovirus genomic evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-856) contains supplementary material, which is available to authorized users.