Cross-species gene-family fluctuations reveal the dynamics of horizontal transfers

Whole-genome sequencing of Fusarium spp. causing sugarcane root rot on both chewing cane and sugar-making cane

Previously we isolated three Fusarium strains (a F. sacchari strain namely GXUF-1, and another two F. commune strains namely GXUF-2 and GXUF-3), and we verified that GXUF-3 was able to cause sugarcane root rot to the chewing cane cultivar Badila. Considering that Fusarium spp. are a group of widely distributed fungal pathogens, we tested whether these three Fusarium isolates were able to cause root rot to Badila as well as sugar-making cane cultivar (Guitang42), using a suitable inoculation method established based on infection assays using Badila. We found that the three Fusarium strains were able to cause root rot symptoms to both Badila and Guitang42, to different extents. To better investigate the potential pathogenicity mechanisms, we performed Illumina high-throughput sequencing and analyzed the whole genomic sequence data of these three Fusarium strains. The results reveal that the assembly sizes of the three Fusarium strains were in a range of 44.7-48.2 Mb, with G + C contents of 48.0-48.5%, and 14,154-15,175 coding genes. The coding genes were annotated by multiple public databases, and potential pathogenic genes were predicted using proprietary databases (such as PHI, DFVF, CAZy, etc.). Furthermore, based on evolutionary analysis of the coding sequence, we found that contraction and expansion of gene families occurred in the three Fusarium strains. Overall, our results suggest a potential risk that the root rot disease may occur to the sugar-making canes although it was initially spotted from fruit cane, and provide clues to understand the pathogenic mechanisms of Fusarium spp. causing sugarcane root rot.

Metacommunity structure preserves genome diversity in the presence of gene-specific selective sweeps under moderate rates of horizontal gene transfer

The horizontal transfer of genes is fundamental for the eco-evolutionary dynamics of microbial communities, such as oceanic plankton, soil, and the human microbiome. In the case of an acquired beneficial gene, classic population genetics would predict a genome-wide selective sweep, whereby the genome spreads clonally within the community and together with the beneficial gene, removing genome diversity. Instead, several sources of metagenomic data show the existence of "gene-specific sweeps", whereby a beneficial gene spreads across a bacterial community, maintaining genome diversity. Several hypotheses have been proposed to explain this process, including the decreasing gene flow between ecologically distant populations, frequency-dependent selection from linked deleterious allelles, and very high rates of horizontal gene transfer. Here, we propose an additional possible scenario grounded in eco-evolutionary principles. Specifically, we show by a mathematical model and simulations that a metacommunity where species can occupy multiple patches, acting together with a realistic (moderate) HGT rate, helps maintain genome diversity. Assuming a scenario of patches dominated by single species, our model predicts that diversity only decreases moderately upon the arrival of a new beneficial gene, and that losses in diversity can be quickly restored. We explore the generic behaviour of diversity as a function of three key parameters, frequency of insertion of new beneficial genes, migration rates and horizontal transfer rates.Our results provides a testable explanation for how diversity can be maintained by gene-specific sweeps even in the absence of high horizontal gene transfer rates.

Taxonomic classification method for metagenomics based on core protein families with Core-Kaiju

Characterizing species diversity and composition of bacteria hosted by biota is revolutionizing our understanding of the role of symbiotic interactions in ecosystems. Determining microbiomes diversity implies the assignment of individual reads to taxa by comparison to reference databases. Although computational methods aimed at identifying the microbe(s) taxa are available, it is well known that inferences using different methods can vary widely depending on various biases. In this study, we first apply and compare different bioinformatics methods based on 16S ribosomal RNA gene and shotgun sequencing to three mock communities of bacteria, of which the compositions are known. We show that none of these methods can infer both the true number of taxa and their abundances. We thus propose a novel approach, named Core-Kaiju, which combines the power of shotgun metagenomics data with a more focused marker gene classification method similar to 16S, but based on emergent statistics of core protein domain families. We thus test the proposed method on various mock communities and we show that Core-Kaiju reliably predicts both number of taxa and abundances. Finally, we apply our method on human gut samples, showing how Core-Kaiju may give more accurate ecological characterization and a fresh view on real microbiomes.

Living Organisms Author Their Read-Write Genomes in Evolution

Evolutionary variations generating phenotypic adaptations and novel taxa resulted from complex cellular activities altering genome content and expression: (i) Symbiogenetic cell mergers producing the mitochondrion-bearing ancestor of eukaryotes and chloroplast-bearing ancestors of photosynthetic eukaryotes; (ii) interspecific hybridizations and genome doublings generating new species and adaptive radiations of higher plants and animals; and, (iii) interspecific horizontal DNA transfer encoding virtually all of the cellular functions between organisms and their viruses in all domains of life. Consequently, assuming that evolutionary processes occur in isolated genomes of individual species has become an unrealistic abstraction. Adaptive variations also involved natural genetic engineering of mobile DNA elements to rewire regulatory networks. In the most highly evolved organisms, biological complexity scales with "non-coding" DNA content more closely than with protein-coding capacity. Coincidentally, we have learned how so-called "non-coding" RNAs that are rich in repetitive mobile DNA sequences are key regulators of complex phenotypes. Both biotic and abiotic ecological challenges serve as triggers for episodes of elevated genome change. The intersections of cell activities, biosphere interactions, horizontal DNA transfers, and non-random Read-Write genome modifications by natural genetic engineering provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

Family-specific scaling laws in bacterial genomes

Among several quantitative invariants found in evolutionary genomics, one of the most striking is the scaling of the overall abundance of proteins, or protein domains, sharing a specific functional annotation across genomes of given size. The size of these functional categories change, on average, as power-laws in the total number of protein-coding genes. Here, we show that such regularities are not restricted to the overall behavior of high-level functional categories, but also exist systematically at the level of single evolutionary families of protein domains. Specifically, the number of proteins within each family follows family-specific scaling laws with genome size. Functionally similar sets of families tend to follow similar scaling laws, but this is not always the case. To understand this systematically, we provide a comprehensive classification of families based on their scaling properties. Additionally, we develop a quantitative score for the heterogeneity of the scaling of families belonging to a given category or predefined group. Under the common reasonable assumption that selection is driven solely or mainly by biological function, these findings point to fine-tuned and interdependent functional roles of specific protein domains, beyond our current functional annotations. This analysis provides a deeper view on the links between evolutionary expansion of protein families and the functional constraints shaping the gene repertoire of bacterial genomes.

Ammonia-Oligotrophic and Diazotrophic Heavy Metal-Resistant Serratia liquefaciens Strains from Pioneer Plants and Mine Tailings

Mine tailings are man-made environments characterized by low levels of organic carbon and assimilable nitrogen, as well as moderate concentrations of heavy metals. For the introduction of nitrogen into these environments, a key role is played by ammonia-oligotrophic/diazotrophic heavy metal-resistant guilds. In mine tailings from Zacatecas, Mexico, Serratia liquefaciens was the dominant heterotrophic culturable species isolated in N-free media from bulk mine tailings as well as the rhizosphere, roots, and aerial parts of pioneer plants. S. liquefaciens strains proved to be a meta-population with high intraspecific genetic diversity and a potential to respond to these extreme conditions. The phenotypic and genotypic features of these strains reveal the potential adaptation of S. liquefaciens to oligotrophic and nitrogen-limited mine tailings with high concentrations of heavy metals. These features include ammonia-oligotrophic growth, nitrogen fixation, siderophore and indoleacetic acid production, phosphate solubilization, biofilm formation, moderate tolerance to heavy metals under conditions of diverse nitrogen availability, and the presence of zntA, amtB, and nifH genes. The acetylene reduction assay suggests low nitrogen-fixing activity. The nifH gene was harbored in a plasmid of ∼60 kb and probably was acquired by a horizontal gene transfer event from Klebsiella variicola.