Deciphering the evolution of microbial interactions: in silico studies of two-member microbial communities

Microbes thrive in communities, embedded in a complex web of interactions. These interactions, particularly metabolic interactions, play a crucial role in maintaining the community structure and function. As the organisms thrive and evolve, a variety of evolutionary processes alter the interactions among the organisms in the community, although the community function remains intact. In this work, we simulate the evolution of two-member microbial communities in silico to study how evolutionary forces can shape the interactions between organisms. We employ genomescale metabolic models of organisms from the human gut, which exhibit a range of interaction patterns, from mutualism to parasitism. We observe that the evolution of microbial interactions varies depending upon the starting interaction and also on the metabolic capabilities of the organisms in the community. We find that evolutionary constraints play a significant role in shaping the dependencies of organisms in the community. Evolution of microbial communities yields fitness benefits in only a small fraction of the communities, and is also dependent on the interaction type of the wild-type communities. The metabolites cross-fed in the wild-type communities appear in only less than 50% of the evolved communities. A wide range of new metabolites are cross-fed as the communities evolve. Further, the dynamics of microbial interactions are not specific to the interaction of the wild-type community but vary depending on the organisms present in the community. Our approach of evolving microbial communities in silico provides an exciting glimpse of the dynamics of microbial interactions and offers several avenues for future investigations.

Systematic investigation of the link between enzyme catalysis and cold adaptation

Cold temperature is prevalent across the biosphere and slows the rates of chemical reactions. Increased catalysis has been predicted to be a dominant adaptive trait of enzymes to reduced temperature, and this expectation has informed physical models for enzyme catalysis and influenced bioprospecting strategies. To systematically test rate enhancement as an adaptive trait to cold, we paired kinetic constants of 2223 enzyme reactions with their organism’s optimal growth temperature (TGrowth) and analyzed trends of rate constants as a function of TGrowth. These data do not support a general increase in rate enhancement in cold adaptation. In the model enzyme ketosteroid isomerase (KSI), there is prior evidence for temperature adaptation from a change in an active site residue that results in a tradeoff between activity and stability. Nevertheless, we found that little of the rate constant variation for 20 KSI variants was accounted for by TGrowth. In contrast, and consistent with prior expectations, we observed a correlation between stability and TGrowth across 433 proteins. These results suggest that temperature exerts a weaker selection pressure on enzyme rate constants than stability and that evolutionary forces other than temperature are responsible for the majority of enzymatic rate constant variation.

Re-Thinking the Roles of Pathogens and Mutualists: Exploring the Continuum of Symbiosis in the Context of Microbial Ecology and Evolution

Systems of classification are important for guiding research activities and providing a common platform for discussion and investigation. One such system is assigning microbial taxa to the roles of mutualists and pathogens. Yet, there are often challenges and even inconsistencies in reports of research findings when microbial taxa display behaviors outside of these two static conditions (e.g. commensal). Over the last two decades, there has been some effort to highlight a continuum of symbiosis, wherein certain microbial taxa may exhibit mutualistic or pathogenic traits depending on environmental contexts, life stages, and plant host associations. However, gaps remain in understanding how to apply the continuum approach to host-microbe pairs across a range of environmental and ecological factors. This commentary presents an alternative framework for evaluating the continuum of symbiosis using dominant archetypes that define symbiotic ranges. We focus particularly on fungi and bacteria, though we recognize that archaea and other microeukaryotes play important roles in host-microbe interactions that may be described by this approach. This framework is centered in eco-evolutionary theory and aims to enhance communication among researchers, as well as prioritize holistic consideration of the factors shaping microbial life strategies. We discuss the influence of plant-mediated factors, habitat constraints, co-evolutionary forces, and the genetic contributions which shape different microbial lifestyles. Looking to the future, using a continuum of symbiosis paradigm will enable greater flexibility in defining the roles of target microbes and facilitate a more holistic view of the complex and dynamic relationship between microbes and plants.

Linkage disequilibrium between rare mutations

The statistical associations between mutations, collectively known as linkage disequilibrium (LD), encode important information about the evolutionary forces acting within a population. Yet in contrast to single-site analogues like the site frequency spectrum, our theoretical understanding of linkage disequilibrium remains limited. In particular, little is currently known about how mutations with different ages and fitness costs contribute to expected patterns of LD, even in simple settings where recombination and genetic drift are the major evolutionary forces. Here, I introduce a forward-time framework for predicting linkage disequilibrium between pairs of neutral and deleterious mutations as a function of their present-day frequencies. I show that the dynamics of linkage disequilibrium become much simpler in the limit that mutations are rare, where they admit a simple heuristic picture based on the trajectories of the underlying lineages. I use this approach to derive analytical expressions for a family of frequency-weighted LD statistics as a function of the recombination rate, the frequency scale, and the additive and epistatic fitness costs of the mutations. I find that the frequency scale can have a dramatic impact on the shapes of the resulting LD curves, reflecting the broad range of time scales over which these correlations arise. I also show that the differences between neutral and deleterious LD are not purely driven by differences in their mutation frequencies, and can instead display qualitative features that are reminiscent of epistasis. I conclude by discussing the implications of these results for recent LD measurements in bacteria. This forward-time approach may provide a useful framework for predicting linkage disequilibrium across a range of evolutionary scenarios.

Comparative studies of codon usage profile of Anisakis simplex (Nematoda) and Carassius gibelio (Prussian carp)

Aim: The aim of the present study was to understand the molecular relationship between nematode (parasite) and fish (host) through codon usage bias (CUB) analysis. Methodology: The Codon usage bias analysis has been performed in fish Carassius gibelio (Prussian carp) and nematode fish parasite Anisakis simplex. The complete coding sequences (CDS) of C. gibelio (Prussian carp) and A. simplex (Nematode) were retrieved from National Center for Biotechnology Information and followed to that we have performed bioinformatics analysis to understand the codon usage pattern between host and parasite. Results: Different CUB indices like Relative synonymous codon usage (RSCU), Effective number of codons (ENC), Codon adaptation index (CAI) and Codon bias index (CBI) revealed a similar pattern in the codon usage in C. gibelio and A. simplex. In addition, inclusive analysis using different plots (ENC, parity, neutrality) had shown the influence of both the evolutionary forces i.e mutational and translational selection on codon usage pattern. This describes the role of evolutionary forces in determining the conserved genome to establish species-specific function-level differences for efficient survival. Interpretation: The present study elucidated the association between Carassiusgibelio (host) and Anisakis simplex (parasite) based on the similar pattern of codon usage bias between both the species.

The evolution of insect visual opsin genes with specific consideration of the influence of ocelli and life history traits

Background Visual opsins are expressed in the compound eyes and ocelli of insects and enable light detection. Three distinct phylogenetic groups of visual opsins are found in insects, named long (LW), short (SW) and ultraviolet (UV) wavelength sensitive opsins. Recently, the LW group was found to be duplicated into the LW2b and the LW2a opsins. The expression of LW2b opsins is ocelli specific in some insects (e.g., bees, cricket, scorpion flies), but the gene was not found in other orders possessing three or less ocelli (e.g., dragonflies, beetles, moths, bugs). In flies, two LW2b homologs have been characterised, with one expressed in the ocelli and the other in the compound eyes. To date, it remains unclear which evolutionary forces have driven gains and losses of LW opsins in insects. Here we take advantage of the recent rapid increase in available sequence data (i.e., from insect genomes, targeted PCR amplification, RNAseq) to characterize the phylogenetic relationships of 1000 opsin sequences in 18 orders of Insects. The resulting phylogeny discriminates between four main groups of opsins, and onto this phylogeny we mapped relevant morphological and life history traits. Results Our results demonstrate a conserved LW2b opsin only present in insects with three ocelli. Only two groups (Brachycera and Odonata) possess more than one LW2b opsin, likely linked to their life history. In flies, we hypothesize that the duplication of the LW2b opsin occurred after the transition from aquatic to terrestrial larvae. During this transition, higher flies (Brachycera) lost a copy of the LW2a opsin, still expressed and duplicated in the compound eyes of lower flies (Nematocera). In higher flies, the LW2b opsin has been duplicated and expressed in the compound eyes while the ocelli and the LW2b opsin were lost in lower flies. In dragonflies, specialisation of flight capabilities likely drove the diversification of the LW2b visual opsins. Conclusion The presence of the LW2b opsin in insects possessing three ocelli suggests a role in specific flight capabilities (e.g., stationary flight). This study provides the most complete view of the evolution of visual opsin genes in insects yet, and provides new insight into the influence of ocelli and life history traits on opsin evolution in insects.

The evolution of sexual dichromatism in a large radiation of landfowls: re-examining the female-biased selection in Wallace's model

Sexual dichromatism, the colour difference between males and females, has been particularly important for studying the interplay between sexual and natural selection. However, previous studies on the evolutionary forces of sexual dichromatism examing the Darwin's and Wallace's model have produced mixed results. Phasianidae is a species-rich family with worldwide distribution, occupancy in nearly all terrestrial habitats, and a wide diversity of plumage patterns and colourations. Here, we use phylogenetic comparative methods to test the relationship between sexual dichromatism and colour complexity of males and females on both evolutionary direction and tempo including all species in Phasianidae. We show that the evolutionary direction of sexual dichromatism is negatively correlated with colour complexity in females but not males, and the evolutionary rates of sexual dichromatism are positively correlated with the evolutionary rates of colour complexity in both sexes. These results highlight the important role of female colour evolution in shaping sexual dichromatism in the pheasant family, and provide strong empirical supports for Wallace's hypothesis via a mosaic of sexual and natural selection in both sexes.

Genomic and Meiotic Changes Accompanying Polyploidization

Hybridization and polyploidy have been considered as significant evolutionary forces in adaptation and speciation, especially among plants. Interspecific gene flow generates novel genetic variants adaptable to different environments, but it is also a gene introgression mechanism in crops to increase their agronomical yield. An estimate of 9% of interspecific hybridization has been reported although the frequency varies among taxa. Homoploid hybrid speciation is rare compared to allopolyploidy. Chromosome doubling after hybridization is the result of cellular defects produced mainly during meiosis. Unreduced gametes, which are formed at an average frequency of 2.52% across species, are the result of altered spindle organization or orientation, disturbed kinetochore functioning, abnormal cytokinesis, or loss of any meiotic division. Meiotic changes and their genetic basis, leading to the cytological diploidization of allopolyploids, are just beginning to be understood especially in wheat. However, the nature and mode of action of homoeologous recombination suppressor genes are poorly understood in other allopolyploids. The merger of two independent genomes causes a deep modification of their architecture, gene expression, and molecular interactions leading to the phenotype. We provide an overview of genomic changes and transcriptomic modifications that particularly occur at the early stages of allopolyploid formation.

Counteracting Forces of Introgressive Hybridization and Interspecific Competition Shape the Morphological Traits of Cryptic Lberian Eptesicus Bats

Cryptic species that coexist in sympatry are likely to simultaneously experience strong competition and hybridization. The first phenomenon would lead to character displacement, whereas the second can potentially promote morphological similarity through adaptive introgression. The main goal of this work was to investigate the effect of introgressive hybridization on the morphology of cryptic Iberian Eptesicus bats when facing counteracting evolutionary forces from interspecific competition. We found substantial overlap both in dentition and in wing morphology traits, though mainly in individuals in sympatry. The presence of hybrids contributes to a fifth of this overlap, with hybrids showing traits with intermediate morphometry. Thus, introgressive hybridization may contribute to species adaptation to trophic and ecological space responding directly to the macro-habitats characteristics of the sympatric zone and to local prey availability. On the other hand, fur shade tended to be browner and brighter in hybrids than parental species. Colour differences could result from partitioning of resources as an adaptation to environmental factors such as roost and microhabitats. We argue that a balance between adaptive introgression and niche partitioning shapes species interactions with the environment through affecting morphological traits under selection.

Effects of homoeologous exchange on gene expression and alternative splicing in a newly formed allotetraploid wheat

Homoeologous exchange (HE) is a major mechanism generating post-polyploidization genetic variation with important evolutionary consequences. However, the direct impacts of HE without entangling with additional evolutionary forces on gene expression remains to be fully understood. Here, we analyzed high-throughput RNA-seq data of young leaves from individuals of a synthetic allotetraploid wheat (AADD), which contain variable numbers of HEs. We aimed to investigate if and to which extent HE directly impacts gene expression and alternative splicing (AS). We found that HE impacts expression of genes located within HE regions primarily via cis-acting dosage effect, which led to significant changes in the total expression of homoeolog pairs, especially for homoeologs whose original expression was biased. In parallel, HE influences expression of a large amount of genes residing in non-HE regions by trans-regulation leading to convergent expression of homoeolog pairs. Intriguingly, when taking into account of the original relative homoeolog expression states, homoeolog pairs under trans-effect are more prone to showing convergent response to HE whereas those under cis-effect trended to show subgenome-specific expression. Moreover, HE induced quantitative, largely individual-specific, changes of alternative splicing (AS) events. Like homoeologs expression, homoeo-AS events which related to trans effect were more responsive to HE. HE therefore exerts multifaceted immediate effects on gene expression and, to a less extent, also transcript diversity in nascent allopolyploidy.