Standard Codon Substitution Models Overestimate Purifying Selection for Nonstationary Data

Nonsynonymous Polymorphism Counts in Bacterial Genomes: a Comparative Examination

Genomic data reveal single-nucleotide polymorphisms (SNPs) that may carry information about the evolutionary history of bacteria. However, it remains unclear what inferences about selection can be made from genomic SNP data. Bacterial species are often sampled during epidemic outbreaks or within hosts during the course of chronic infections. SNPs obtained from genomic analysis of these data are not necessarily fixed. Treating them as fixed during analysis by using measures such as the ratio of nonsynonymous to synonymous evolutionary changes (dN/dS) may lead to incorrect inferences about the strength and direction of selection. In this study, we consider data from a range of whole-genome sequencing studies of bacterial pathogens and explore patterns of nonsynonymous variation to assess whether evidence of selection can be identified by investigating SNP counts alone across multiple WGS studies. We visualize these SNP data in ways that highlight their relationship to neutral baseline expectations. These neutral expectations are based on a simple model of mutation, from which we simulate SNP accumulation to investigate how SNP counts are distributed under alternative assumptions about positive and negative selection. We compare these patterns with empirical SNP data and illustrate the general difficulty of detecting positive selection from SNP data. Finally, we consider whether SNP counts observed at the between-host population level differ from those observed at the within-host level and find some evidence that suggests that dynamics across these two scales are driven by different underlying processes.IMPORTANCE Identifying selection from SNP data obtained from whole-genome sequencing studies is challenging. Some current measures used to identify and quantify selection acting on genomes rely on fixed differences; thus, these are inappropriate for SNP data where variants are not fixed. With the increase in whole-genome sequencing studies, it is important to consider SNP data in the context of evolutionary processes. How SNPs are counted and analyzed can help in understanding mutation accumulation and trajectories of strains. We developed a tool for identifying possible evidence of selection and for comparative analysis with other SNP data. We propose a model that provides a rule-of-thumb guideline and two new visualization techniques that can be used to interpret and compare SNP data. We quantify the expected proportion of nonsynonymous SNPs in coding regions under neutrality and demonstrate its use in identifying evidence of positive and negative selection from simulations and empirical data.

Genome-wide identification and characterization of the soybean SOD family during alkaline stress

Background

Superoxide dismutase (SOD) proteins, as one kind of the antioxidant enzymes, play critical roles in plant response to various environment stresses. Even though its functions in the oxidative stress were very well characterized, the roles of SOD family genes in regulating alkaline stress response are not fully reported.

Methods

We identified the potential family members by using Hidden Markov model and soybean genome database. The neighbor-joining phylogenetic tree and exon-intron structures were generated by using software MEGA 5.0 and GSDS online server, respectively. Furthermore, the conserved motifs were analyzed by MEME online server. The syntenic analysis was conducted using Circos-0.69. Additionally, the expression levels of soybean SOD genes under alkaline stress were identified by qRT-PCR.

Results

In this study, we identified 13 potential SOD genes in soybean genome. Phylogenetic analysis suggested that SOD genes could be classified into three subfamilies, including MnSODs (GmMSD1-2), FeSODs (GmFSD1-5) and Cu/ZnSODs (GmCSD1-6). We further investigated the gene structure, chromosomal locations and gene-duplication, conserved domains and promoter cis-elements of the soybean SOD genes. We also explored the expression profiles of soybean SOD genes in different tissues and alkaline, salt and cold stresses, based on the transcriptome data. In addition, we detected their expression patterns in roots and leaves by qRT-PCR under alkaline stress, and found that different SOD subfamily genes may play different roles in response to alkaline stress. These results also confirmed the hypothesis that the great evolutionary divergence may contribute to the potential functional diversity in soybean SOD genes. Taken together, we established a foundation for further functional characterization of soybean SOD genes in response to alkaline stress in the future.

Repertoire-wide phylogenetic models of B cell molecular evolution reveal evolutionary signatures of aging and vaccination

In order to produce effective antibodies, B cells undergo rapid somatic hypermutation (SHM) and selection for binding affinity to antigen via a process called affinity maturation. The similarities between this process and evolution by natural selection have led many groups to use phylogenetic methods to characterize the development of immunological memory, vaccination, and other processes that depend on affinity maturation. However, these applications are limited by the fact that most phylogenetic models are designed to be applied to individual lineages comprising genetically diverse sequences, while B cell repertoires often consist of hundreds to thousands of separate low-diversity lineages. Further, several features of affinity maturation violate important assumptions in standard phylogenetic models. Here, we introduce a hierarchical phylogenetic framework that integrates information from all lineages in a repertoire to more precisely estimate model parameters while simultaneously incorporating the unique features of SHM. We demonstrate the power of this repertoire-wide approach by characterizing previously undescribed phenomena in affinity maturation. First, we find evidence consistent with age-related changes in SHM hot-spot targeting. Second, we identify a consistent relationship between increased tree length and signs of increased negative selection, apparent in the repertoires of recently vaccinated subjects and those without any known recent infections or vaccinations. This suggests that B cell lineages shift toward negative selection over time as a general feature of affinity maturation. Our study provides a framework for undertaking repertoire-wide phylogenetic testing of SHM hypotheses and provides a means of characterizing dynamics of mutation and selection during affinity maturation.

What is an archaeon and are the Archaea really unique?

The recognition of the group Archaea as a major branch of the tree of life (ToL) prompted a new view of the evolution of biodiversity. The genomic representation of archaeal biodiversity has since significantly increased. In addition, advances in phylogenetic modeling of multi-locus datasets have resolved many recalcitrant branches of the ToL. Despite the technical advances and an expanded taxonomic representation, two important aspects of the origins and evolution of the Archaea remain controversial, even as we celebrate the 40th anniversary of the monumental discovery. These issues concern (i) the uniqueness (monophyly) of the Archaea, and (ii) the evolutionary relationships of the Archaea to the Bacteria and the Eukarya; both of these are relevant to the deep structure of the ToL. To explore the causes for this persistent ambiguity, I examine multiple datasets and different phylogenetic approaches that support contradicting conclusions. I find that the uncertainty is primarily due to a scarcity of information in standard datasets-universal core-genes datasets-to reliably resolve the conflicts. These conflicts can be resolved efficiently by comparing patterns of variation in the distribution of functional genomic signatures, which are less diffused unlike patterns of primary sequence variation. Relatively lower heterogeneity in distribution patterns minimizes uncertainties and supports statistically robust phylogenetic inferences, especially of the earliest divergences of life. This case study further highlights the limitations of primary sequence data in resolving difficult phylogenetic problems, and raises questions about evolutionary inferences drawn from the analyses of sequence alignments of a small set of core genes. In particular, the findings of this study corroborate the growing consensus that reversible substitution mutations may not be optimal phylogenetic markers for resolving early divergences in the ToL, nor for determining the polarity of evolutionary transitions across the ToL.