Mycobacterial phylogenomics: an enhanced method for gene turnover analysis reveals uneven levels of gene gain and loss among species and gene families

Nonsynonymous Synonymous Variants Demand for a Paradigm Shift in Genetics

Synonymous (also known as silent) variations are by definition not considered to change the coded protein. Still many variations in this category affect either protein abundance or properties. As this situation is confusing, we have recently introduced systematics for synonymous variations and those that may on the surface look like synonymous, but these may affect the coded protein in various ways. A new category, unsense variation, was introduced to describe variants that do not introduce a stop codon into the variation site, but which lead to different types of changes in the coded protein. Many of these variations lead to mRNA degradation and missing protein. Here, consequences of the systematics are discussed from the perspectives of variation annotation and interpretation, evolutionary calculations, nonsynonymous-to-synonymous substitution rates, phylogenetics and other evolutionary inferences that are based on the principle of (nearly) neutral synonymous variations. It may be necessary to reassess published results. Further, databases for synonymous variations and prediction methods for such variations should consider unsense variations. Thus, there is a need to evaluate and reflect principles of numerous aspects in genetics, ranging from variation naming and classification to evolutionary calculations.

Evolution and Expression of the Immune System of a Facultatively Anadromous Salmonid

Vertebrates have evolved a complex immune system required for the identification of and coordinated response to harmful pathogens. Migratory species spend periods of their life-cycle in more than one environment, and their immune system consequently faces a greater diversity of pathogens residing in different environments. In facultatively anadromous salmonids, individuals may spend parts of their life-cycle in freshwater and marine environments. For species such as the brown trout Salmo trutta, sexes differ in their life-histories with females more likely to migrate to sea while males are more likely to stay and complete their life-cycle in their natal river. Salmonids have also undergone a lineage-specific whole genome duplication event, which may provide novel immune innovations but our current understanding of the differences in salmonid immune expression between the sexes is limited. We characterized the brown trout immune gene repertoire, identifying a number of canonical immune genes in non-salmonid teleosts to be duplicated in S. trutta, with genes involved in innate and adaptive immunity. Through genome-wide transcriptional profiling (“RNA-seq”) of male and female livers to investigate sex differences in gene expression amplitude and alternative splicing, we identified immune genes as being generally male-biased in expression. Our study provides important insights into the evolutionary consequences of whole genome duplication events on the salmonid immune gene repertoire and how the sexes differ in constitutive immune expression.

The many lives of nontuberculous mycobacteria

Nontuberculous mycobacteria (NTM) include species that colonize human epithelia, as well as species that are ubiquitous in soil and aquatic environments. NTM that primarily inhabit soil and aquatic environments include the Mycobacterium avium complex (MAC, M. avium and Mycobacterium intracellulare) and the Mycobacterium abscessus complex (MABSC, M. abscessus subspecies abscessus, massiliense, and bolletii), and can be free-living, biofilm-associated, or amoeba-associated. Although NTM are rarely pathogenic in immunocompetent individuals, those who are immunocompromised - due to either an inherited or acquired immunodeficiency - are highly susceptible to NTM infection (NTMI). Several characteristics such as biofilm formation and the ability of select NTM species to form distinct colony morphotypes all may play a role in pathogenesis not observed in the related, well-characterized pathogen Mycobacterium tuberculosis The recognition of different morphotypes of NTM has been established and characterized since the 1950s, but the mechanisms that underlie colony phenotype change and subsequent differences in pathogenicity are just beginning to be explored. Advances in genomic analysis have led to progress in identifying genes important to the pathogenesis and persistence of MAC disease as well as illuminating genetic aspects of different colony morphotypes. Here we review recent literature regarding NTM ecology and transmission, as well as the factors which regulate colony morphotype and pathogenicity.

Genomic Comparisons Reveal Microevolutionary Differences in Mycobacterium abscessus Subspecies

Mycobacterium abscessus, a rapid-growing non-tuberculous mycobacterium, has been the cause of sporadic and outbreak infections world-wide. The subspecies in M. abscessus complex (M. abscessus, M. massiliense, and M. bolletii) are associated with different biologic and pathogenic characteristics and are known to be among the most frequently isolated opportunistic pathogens from clinical material. To date, the evolutionary forces that could have contributed to these biological and clinical differences are still unclear. We compared genome data from 243 M. abscessus strains downloaded from the NCBI ftp Refseq database to understand how the microevolutionary processes of homologous recombination and positive selection influenced the diversification of the M. abscessus complex at the subspecies level. The three subspecies are clearly separated in the Minimum Spanning Tree. Their MUMi-based genomic distances support the separation of M. massiliense and M. bolletii into two subspecies. Maximum Likelihood analysis through dN/dS (the ratio of number of non-synonymous substitutions per non-synonymous site, to the number of synonymous substitutions per synonymous site) identified distinct genes in each subspecies that could have been affected by positive selection during evolution. The results of genome-wide alignment based on concatenated locally-collinear blocks suggest that (a) recombination has affected the M. abscessus complex more than mutation and positive selection; (b) recombination occurred more frequently in M. massiliense than in the other two subspecies; and (c) the recombined segments in the three subspecies have come from different intra-species and inter-species origins. The results lead to the identification of possible gene sets that could have been responsible for the subspecies-specific features and suggest independent evolution among the three subspecies, with recombination playing a more significant role than positive selection in the diversification among members in this complex.

Rapid Gene Turnover as a Significant Source of Genetic Variation in a Recently Seeded Population of a Healthcare-Associated Pathogen

Genome sequencing has been useful to gain an understanding of bacterial evolution. It has been used for studying the phylogeography and/or the impact of mutation and recombination on bacterial populations. However, it has rarely been used to study gene turnover at microevolutionary scales. Here, we sequenced Mexican strains of the human pathogen Acinetobacter baumannii sampled from the same locale over a 3 year period to obtain insights into the microevolutionary dynamics of gene content variability. We found that the Mexican A. baumannii population was recently founded and has been emerging due to a rapid clonal expansion. Furthermore, we noticed that on average the Mexican strains differed from each other by over 300 genes and, notably, this gene content variation has accrued more frequently and faster than the accumulation of mutations. Moreover, due to its rapid pace, gene content variation reflects the phylogeny only at very short periods of time. Additionally, we found that the external branches of the phylogeny had almost 100 more genes than the internal branches. All in all, these results show that rapid gene turnover has been of paramount importance in producing genetic variation within this population and demonstrate the utility of genome sequencing to study alternative forms of genetic variation.

Adaptive Evolution of Extreme Acidophile Sulfobacillus thermosulfidooxidans Potentially Driven by Horizontal Gene Transfer and Gene Loss

Recent phylogenomic analysis has suggested that three strains isolated from different copper mine tailings around the world were taxonomically affiliated with Sulfobacillusthermosulfidooxidans Here, we present a detailed investigation of their genomic features, particularly with respect to metabolic potentials and stress tolerance mechanisms. Comprehensive analysis of the Sulfobacillus genomes identified a core set of essential genes with specialized biological functions in the survival of acidophiles in their habitats, despite differences in their metabolic pathways. The Sulfobacillus strains also showed evidence for stress management, thereby enabling them to efficiently respond to harsh environments. Further analysis of metabolic profiles provided novel insights into the presence of genomic streamlining, highlighting the importance of gene loss as a main mechanism that potentially contributes to cellular economization. Another important evolutionary force, especially in larger genomes, is gene acquisition via horizontal gene transfer (HGT), which might play a crucial role in the recruitment of novel functionalities. Also, a successful integration of genes acquired from archaeal donors appears to be an effective way of enhancing the adaptive capacity to cope with environmental changes. Taken together, the findings of this study significantly expand the spectrum of HGT and genome reduction in shaping the evolutionary history of Sulfobacillus strains.IMPORTANCE Horizontal gene transfer (HGT) and gene loss are recognized as major driving forces that contribute to the adaptive evolution of microbial genomes, although their relative importance remains elusive. The findings of this study suggest that highly frequent gene turnovers within microorganisms via HGT were necessary to incur additional novel functionalities to increase the capacity of acidophiles to adapt to changing environments. Evidence also reveals a fascinating phenomenon of potential cross-kingdom HGT. Furthermore, genome streamlining may be a critical force in driving the evolution of microbial genomes. Taken together, this study provides insights into the importance of both HGT and gene loss in the evolution and diversification of bacterial genomes.

Comparative genome analyses of mycobacteria give better insights into their evolution

Mycobacteria a genus of Actinobacteria are widespread in nature ranging from soil-dwelling saprophytes to human and animal pathogens. The rate of growth has been a classifying factor for the Mycobacterium spp., dividing them into the rapid growers and the slow growers. Here we have performed a comparative genome study of mycobacterial species in order to get better understanding of their evolution, particularly to understand the distinction between the rapid and slow growers. Our study shows that the slow growers had generally gained and lost more genes compared to the rapid growers. The slow growers might haved eventually lost genes (LivFGMH operon, shaACDEFG genes and MspA porin) that could contribute to the slow growth rate of the slow growers. The genes gained and lost in mycobacteria had eventually helped these bacteria to adapt to different environments and have led to the evolution of the present day rapid and slow growers. Our results also show high number of Mycobacterium abscessus specific genes (811 genes) and some of them are associated with the known bacterial quorum sensing genes that might be important for Mycobacterium abscessus to adapt and survive in variety of unfavorable environments. Mycobacterium abscessus also does not contains genes involved in the bacterial defense system and together with the quorum sensing genes may have contributed to the high gene gain rate of Mycobacterium abscessus.

Gene Turnover Contributes to the Evolutionary Adaptation of Acidithiobacillus caldus: Insights from Comparative Genomics

Acidithiobacillus caldus is an extremely acidophilic sulfur-oxidizer with specialized characteristics, such as tolerance to low pH and heavy metal resistance. To gain novel insights into its genetic complexity, we chosen six A. caldus strains for comparative survey. All strains analyzed in this study differ in geographic origins as well as in ecological preferences. Based on phylogenomic analysis, we clustered the six A. caldus strains isolated from various ecological niches into two groups: group 1 strains with smaller genomes and group 2 strains with larger genomes. We found no obvious intraspecific divergence with respect to predicted genes that are related to central metabolism and stress management strategies between these two groups. Although numerous highly homogeneous genes were observed, high genetic diversity was also detected. Preliminary inspection provided a first glimpse of the potential correlation between intraspecific diversity at the genome level and environmental variation, especially geochemical conditions. Evolutionary genetic analyses further showed evidence that the difference in environmental conditions might be a crucial factor to drive the divergent evolution of A. caldus species. We identified a diverse pool of mobile genetic elements including insertion sequences and genomic islands, which suggests a high frequency of genetic exchange in these harsh habitats. Comprehensive analysis revealed that gene gains and losses were both dominant evolutionary forces that directed the genomic diversification of A. caldus species. For instance, horizontal gene transfer and gene duplication events in group 2 strains might contribute to an increase in microbial DNA content and novel functions. Moreover, genomes undergo extensive changes in group 1 strains such as removal of potential non-functional DNA, which results in the formation of compact and streamlined genomes. Taken together, the findings presented herein show highly frequent gene turnover of A. caldus species that inhabit extremely acidic environments, and shed new light on the contribution of gene turnover to the evolutionary adaptation of acidophiles.

Evolution by gene loss

The recent increase in genomic data is revealing an unexpected perspective of gene loss as a pervasive source of genetic variation that can cause adaptive phenotypic diversity. This novel perspective of gene loss is raising new fundamental questions. How relevant has gene loss been in the divergence of phyla? How do genes change from being essential to dispensable and finally to being lost? Is gene loss mostly neutral, or can it be an effective way of adaptation? These questions are addressed, and insights are discussed from genomic studies of gene loss in populations and their relevance in evolutionary biology and biomedicine.

Trends in substitution models of molecular evolution

Substitution models of evolution describe the process of genetic variation through fixed mutations and constitute the basis of the evolutionary analysis at the molecular level. Almost 40 years after the development of first substitution models, highly sophisticated, and data-specific substitution models continue emerging with the aim of better mimicking real evolutionary processes. Here I describe current trends in substitution models of DNA, codon and amino acid sequence evolution, including advantages and pitfalls of the most popular models. The perspective concludes that despite the large number of currently available substitution models, further research is required for more realistic modeling, especially for DNA coding and amino acid data. Additionally, the development of more accurate complex models should be coupled with new implementations and improvements of methods and frameworks for substitution model selection and downstream evolutionary analysis.

High Gene Family Turnover Rates and Gene Space Adaptation in the Compact Genome of the Carnivorous Plant Utricularia gibba

Utricularia gibba is an aquatic carnivorous plant with highly specialized morphology, featuring fibrous floating networks of branches and leaf-like organs, no recognizable roots, and bladder traps that capture and digest prey. We recently described the compressed genome of U. gibba as sufficient to control the development and reproduction of a complex organism. We hypothesized intense deletion pressure as a mechanism whereby most noncoding DNA was deleted, despite evidence for three independent whole-genome duplications (WGDs). Here, we explore the impact of intense genome fractionation in the evolutionary dynamics of U. gibba's functional gene space. We analyze U. gibba gene family turnover by modeling gene gain/death rates under a maximum-likelihood statistical framework. In accord with our deletion pressure hypothesis, we show that the U. gibba gene death rate is significantly higher than those of four other eudicot species. Interestingly, the gene gain rate is also significantly higher, likely reflecting the occurrence of multiple WGDs and possibly also small-scale genome duplications. Gene ontology enrichment analyses of U. gibba-specific two-gene orthogroups, multigene orthogroups, and singletons highlight functions that may represent adaptations in an aquatic carnivorous plant. We further discuss two homeodomain transcription factor gene families (WOX and HDG/HDZIP-IV) showing conspicuous differential expansions and contractions in U. gibba. Our results 1) reconcile the compactness of the U. gibba genome with its accommodation of a typical number of genes for a plant genome, and 2) highlight the role of high gene family turnover in the evolutionary diversification of U. gibba's functional gene space and adaptations to its unique lifestyle and highly specialized body plan.