Invertebrate Species with Nonpelagic Larvae Have Elevated Levels of Nonsynonymous Substitutions and Reduced Nucleotide Diversities

Under a nearly neutral model in which most amino acid substitutions are slightly deleterious, variation in demography, population structure, and other ecological factors among closely related species can potentially modify the effective population size or the selective regime, leading to differences in the rate of nonsynonymous substitution. Ratios of nonsynonymous to synonymous substitutions (d_N/d_S) between species were analyzed in a sea star genus (Patiriella) and a molluscan genus (Littorina), each with diverse modes of reproduction, including multiple lineages with pelagic and nonpelagic larvae. In both genera, lineages with nonpelagic larvae had significantly higher d_N/d_S ratios than lineages with pelagic larvae. The hypothesis that the elevated d_N/d_S ratios in species with nonpelagic larvae was due to reduced effective population size was tested by comparing nucleotide diversities in three genera of gastropod mollusks (Littorina, Crepidula, and Hydrobia), each with several modes of reproduction. Overall, there was a significant (p < 0.05) reduction in nucleotide diversity in species with nonpelagic larvae compared to species with pelagic larvae.     

Kr/Kc but not dN/dS correlates positively with body mass in birds,raising implications for inferring lineage-specific selection

Background

The ratio of the rates of non-synonymous and synonymous substitution (d_N/d_S) is commonly used to estimate selection in coding sequences. It is often suggested that, all else being equal, d_N/d_S should be lower in populations with large effective size (N_e) due to increased efficacy of purifying selection. As N_e is difficult to measure directly, life history traits such as body mass, which is typically negatively associated with population size, have commonly been used as proxies in empirical tests of this hypothesis. However, evidence of whether the expected positive correlation between body mass and d_N/d_S is consistently observed is conflicting.

Results

Employing whole genome sequence data from 48 avian species, we assess the relationship between rates of molecular evolution and life history in birds. We find a negative correlation between d_N/d_S and body mass, contrary to nearly neutral expectation. This raises the question whether the correlation might be a method artefact. We therefore in turn consider non-stationary base composition, divergence time and saturation as possible explanations, but find no clear patterns. However, in striking contrast to d_N/d_S, the ratio of radical to conservative amino acid substitutions (K_r/K_c) correlates positively with body mass.

Conclusions

Our results in principle accord with the notion that non-synonymous substitutions causing radical amino acid changes are more efficiently removed by selection in large populations, consistent with nearly neutral theory. These findings have implications for the use of d_N/d_S and suggest that caution is warranted when drawing conclusions about lineage-specific modes of protein evolution using this metric.

Electronic supplementary material

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

Adaptive evolution of insect selective excitatory β-type sodium channel neurotoxins from scorpion venom

Insect selective excitatory β-type sodium channel neurotoxins from scorpion venom (β-NaScTxs) are composed of about 70–76 amino acid residues and share a common scaffold stabilized by four unique disulfide bonds. The phylogenetic analysis of these toxins was hindered by limited sequence data. In our recent study, two new insect selective excitatory β-NaScTxs, LmIT and ImIT, were isolated from Lychas mucronatus and Isometrus maculatus, respectively. With the sequences previously reported, we examined the adaptive molecular evolution of insect selective excitatory β-NaScTxs by estimating the nonsynonymous-to-synonymous rate ratio (ω = d_N/d_S). The results revealed 12 positively selected sites in the genes of insect selective excitatory β-NaScTxs. Moreover, these positively selected sites match well with the sites important for interacting with sodium channels, as demonstrated in previous mutagenesis study. These results reveal that adaptive evolution after gene duplication is one of the most important genetic mechanisms of scorpion neurotoxin diversification.     

Testing single-sample estimators of effective population size in genetically structured populations

The effective population size (N _e ) is a key parameter in evolutionary and population genetics. Single-sample N _e estimation provides an alternative to traditional approaches requiring two or more samples. Single-sample methods assume that the study population has no genetic sub-structure, which is unlikely to be true in wild populations. Here we empirically investigated two single-sample estimators (onesamp and L d N e) in replicated and controlled genetically structured populations of Drosophila melanogaster. Using experimentally controlled population parameters, we calculated the Wright–Fisher expected N _e for the structured population ( ^Total N _e ) and demonstrated that the loss of heterozygosity did not significantly differ from Wright’s model. We found that disregarding the population substructure resulted in ^Total N _e estimates with a low coefficient of variation but these estimates were systematically lower than the expected values, whereas hierarchical estimates accounting for population structure were closer to the expected values but had a higher coefficient of variation. Analysis of simulated populations demonstrated that incomplete sampling, initial allelic diversity and balancing selection may have contributed to deviations from the Wright–Fisher model. Overall the approximate-Bayesian onesamp method performed better than L d N e (with appropriate priors). Both methods performed best when dispersal rates were high and the population structure was approaching panmixia.     

The Genomic Selfing Syndrome Accompanies the Evolutionary Breakdown of Heterostyly

The evolutionary transition from outcrossing to selfing can have important genomic consequences. Decreased effective population size and the reduced efficacy of selection are predicted to play an important role in the molecular evolution of the genomes of selfing species. We investigated evidence for molecular signatures of the genomic selfing syndrome using 66 species of Primula including distylous (outcrossing) and derived homostylous (selfing) taxa. We complemented our comparative analysis with a microevolutionary study of P. chungensis, which is polymorphic for mating system and consists of both distylous and homostylous populations. We generated chloroplast and nuclear genomic data sets for distylous, homostylous, and distylous–homostylous species and identified patterns of nonsynonymous to synonymous divergence (d_N/d_S) and polymorphism (π_N/π_S) in species or lineages with contrasting mating systems. Our analysis of coding sequence divergence and polymorphism detected strongly reduced genetic diversity and heterozygosity, decreased efficacy of purifying selection, purging of large-effect deleterious mutations, and lower rates of adaptive evolution in samples from homostylous compared with distylous populations, consistent with theoretical expectations of the genomic selfing syndrome. Our results demonstrate that self-fertilization is a major driver of molecular evolutionary processes with genomic signatures of selfing evident in both old and relatively young homostylous populations.     

The Trouble with Sliding Windows and the Selective Pressure in BRCA1

Sliding-window analysis has widely been used to uncover synonymous (silent, d_S) and nonsynonymous (replacement, d_N) rate variation along the protein sequence and to detect regions of a protein under selective constraint (indicated by d_N<d_S) or positive selection (indicated by d_N>d_S). The approach compares two or more protein-coding genes and plots estimates dˆ _S and dˆ _N from each sliding window along the sequence. Here we demonstrate that the approach produces artifactual trends of synonymous and nonsynonymous rate variation, with greater variation in dˆ _S than in dˆ _N. Such trends are generated even if the true d_S and d_N are constant along the whole protein and different codons are evolving independently. Many published tests of negative and positive selection using sliding windows that we have examined appear to be invalid because they fail to correct for multiple testing. Instead, likelihood ratio tests provide a more rigorous framework for detecting signals of natural selection affecting protein evolution. We demonstrate that a previous finding that a particular region of the BRCA1 gene experienced a synonymous rate reduction driven by purifying selection is likely an artifact of the sliding window analysis. We evaluate various sliding-window analyses in molecular evolution, population genetics, and comparative genomics, and argue that the approach is not generally valid if it is not known a priori that a trend exists and if no correction for multiple testing is applied.     

Difference in larval type explains patterns of nonsynonymous substitutions in two ancient paralogs of the histone H3 gene in sea stars

SUMMARY Paralogous genes frequently show differences in patterns and rates of substitution that are typically attributed to different selection regimes, mutation rates, or local recombination rates. Here, two anciently diverged paralogous copies of the histone H3 gene in sea stars, the tandem‐repetitive early‐stage gene and a newly isolated gene with lower copy number that was termed the “putative late‐stage histone H3 gene” were analyzed in 69 species with varying mode of larval development. The two genes showed differences in relative copy number, overall substitution rates, nucleotide composition, and codon usage, but similar patterns of relative nonsynonymous substitution rates, when analyzed by the d_N/d_S ratio. Sea stars with a nonpelagic and nonfeeding larval type (i.e., brooding lineages) were observed to have d_N/d_S ratios that were larger than for nonbrooders but equal between the two paralogs. This finding suggested that demographic differences between brooding and nonbrooding lineages were responsible for the elevated d_N/d_S ratios observed for brooders and refuted a suggestion from a previous analysis of the early‐stage gene that the excess nonsynonymous substitutions were due to either (1) gene expression differences at the larval stage between brooders and nonbrooders or (2) the highly repetitive structure of the early‐stage histone H3 gene.     

Molecular clock in neutral protein evolution

Background

A frequent observation in molecular evolution is that amino-acid substitution rates show an index of dispersion (that is, ratio of variance to mean) substantially larger than one. This observation has been termed the overdispersed molecular clock. On the basis of in silico protein-evolution experiments, Bastolla and coworkers recently proposed an explanation for this observation: Proteins drift in neutral space, and can temporarily get trapped in regions of substantially reduced neutrality. In these regions, substitution rates are suppressed, which results in an overall substitution process that is not Poissonian. However, the simulation method of Bastolla et al. is representative only for cases in which the product of mutation rate μ and population size N_e is small. How the substitution process behaves when μN_e is large is not known.

Results

Here, I study the behavior of the molecular clock in in silico protein evolution as a function of mutation rate and population size. I find that the index of dispersion decays with increasing μN_e, and approaches 1 for large μN_e . This observation can be explained with the selective pressure for mutational robustness, which is effective when μN_e is large. This pressure keeps the population out of low-neutrality traps, and thus steadies the ticking of the molecular clock.

Conclusions

The molecular clock in neutral protein evolution can fall into two distinct regimes, a strongly overdispersed one for small μN_e, and a mostly Poissonian one for large μN_e. The former is relevant for the majority of organisms in the plant and animal kingdom, and the latter may be relevant for RNA viruses.

     

A Survey of the Molecular Evolutionary Dynamics of Twenty-Five Multigene Families from Four Grass Taxa

We surveyed the molecular evolutionary characteristics of 25 plant gene families, with the goal of better understanding general processes in plant gene family evolution. The survey was based on 247 GenBank sequences representing four grass species (maize, rice, wheat, and barley). For each gene family, orthology and paralogy relationships were uncertain. Recognizing this uncertainty, we characterized the molecular evolution of each gene family in four ways. First, we calculated the ratio of nonsynonymous to synonymous substitutions (d _N/d _S) both on branches of gene phylogenies and across codons. Our results indicated that the d _N/d _S ratio was statistically heterogeneous across branches in 17 of 25 (68%) gene families. The vast majority of d _N/d _S estimates were <<1.0, suggestive of selective constraint on amino acid replacements, and no estimates were >1.0, either across phylogenetic lineages or across codons. Second, we tested separately for nonsynonymous and synonymous molecular clocks. Sixty-eight percent of gene families rejected a nonsynonymous molecular clock, and 52% of gene families rejected a synonymous molecular clock. Thus, most gene families in this study deviated from clock-like evolution at either synonymous or nonsynonymous sites. Third, we calculated the effective number of codons and the proportion of G+C synonymous sites for each sequence in each gene family. One or both quantities vary significantly within 18 of 25 gene families. Finally, we tested for gene conversion, and only six gene families provided evidence of gene conversion events. Altogether, evolution for these 25 gene families is marked by selective constraint that varies among gene family members, a lack of molecular clock at both synonymous and nonsynonymous sites, and substantial variation in codon usage. Received: 25 May 2000 / Accepted: 16 October 2000     

Significance of Population Size on the Fixation of Nonsynonymous Mutations in Genes Under Varying Levels of Selection Pressure

Previous studies observed a higher ratio of divergences at nonsynonymous and synonymous sites (ω = d_N/d_S) in species with a small population size compared to that estimated for those with a large population size. Here we examined the theoretical relationship between ω, effective population size (N_e), and selection coefficient (s). Our analysis revealed that when purifying selection is high, ω of species with small N_e is much higher than that of species with large N_e. However the difference between the two ω reduces with the decline in selection pressure (s → 0). We examined this relationship using primate and rodent genes and found that the ω estimated for highly constrained genes of primates was up to 2.9 times higher than that obtained for their orthologous rodent genes. Conversely, for genes under weak purifying selection the ω of primates was only 17% higher than that of rodents. When tissue specificity was used as a proxy for selection pressure we found that the ω of broadly expressed genes of primates was up to 2.1-fold higher than that of their rodent counterparts and this difference was only 27% for tissue specific genes. Since most of the nonsynonymous mutations in constrained or broadly expressed genes are deleterious, fixation of these mutations is influenced by N_e. This results in a higher ω of these genes in primates compared to those from rodents. Conversely, the majority of nonsynonymous mutations in less-constrained or tissue-specific genes are neutral or nearly neutral and therefore fixation of them is largely independent of N_e, which leads to the similarity of ω in primates and rodents.     

ASPM and the Evolution of Cerebral Cortical Size in a Community of New World Monkeys

The ASPM (abnormal spindle-like microcephaly associated) gene has been proposed as a major determinant of cerebral cortical size among primates, including humans. Yet the specific functions of ASPM and its connection to human intelligence remain controversial. This debate is limited in part by a taxonomic focus on Old World monkeys and apes. Here we expand the comparative context of ASPM sequence analyses with a study of New World monkeys, a radiation of primates in which enlarged brain size has evolved in parallel in spider monkeys (genus Ateles) and capuchins (genus Cebus). The primate community of Costa Rica is perhaps a model system because it allows for independent pairwise comparisons of smaller- and larger-brained species within two taxonomic families. Accordingly, we analyzed the complete sequence of exon 18 of ASPM in Ateles geoffroyi, Alouatta palliata, Cebus capucinus, and Saimiri oerstedii. As the analysis of multiple species in a genus improves phylogenetic reconstruction, we also analyzed eleven published sequences from other New World monkeys. Our exon-wide, lineage-specific analysis of eleven genera and the ratio of rates of nonsynonymous to synonymous substitutions (d_N/d_S) on ASPM revealed no detectable evidence for positive selection in the lineages leading to Ateles or Cebus, as indicated by d_N/d_S ratios of <1.0 (0.6502 and 0.4268, respectively). Our results suggest that a multitude of interacting genes have driven the evolution of larger brains among primates, with different genes involved in this process in different encephalized lineages, or at least with evidence for positive selection not readily apparent for the same genes in all lineages. The primate community of Costa Rica may serve as a model system for future studies that aim to elucidate the molecular mechanisms underlying cognitive capacity and cortical size.     

Effective population size and the rate and pattern of nucleotide substitutions

Both the overall rate of nucleotide substitution and the relative proportions of synonymous and non-synonymous substitutions are predicted to vary between species that differ in effective population size (N_e). Our understanding of the genetic processes underlying these lineage-specific differences in molecular evolution is still developing. Empirical analyses indicate that variation in substitution rates and patterns caused by differences in N_e is often substantial, however, and must be accounted for in analyses of molecular evolution.     

Examining the Role of Effective Population Size on Mitochondrial and Multilocus Divergence Time Discordance in a Songbird

Estimates of speciation times are subject to a number of potential errors. One source of bias is that effective population size (N_e) has been shown to influence substitution rates. This issue is of particular interest for phylogeographic studies because population sizes can vary dramatically among genetically structured populations across species’ ranges. In this study, we used multilocus data to examine temporal phylogeographic patterns in a widespread North American songbird, the Northern Cardinal (Cardinalis cardinalis). Species tree estimation indicated that the phylogeographic structure of C. cardinalis was comprised of four well-supported mainland lineages with large population sizes (large N_e) and two island lineages comprised of much smaller populations (small N_e). We inferred speciation times from mtDNA and multilocus data and found there was discordance between events that represented island-mainland divergences, whereas both estimates were similar for divergences among mainland lineages. We performed coalescent simulations and found that the difference in speciation times could be attributed to stochasticity for a recently diverged island lineage. However, the magnitude of the change between speciation times estimated from mtDNA and multilocus data of an older island lineage was substantially greater than predicted by coalescent simulations. For this divergence, we found the discordance in time estimates was due to a substantial increase in the mtDNA substitution rate in the small island population. These findings indicate that in phylogeographic studies the relative tempo of evolution between mtDNA and nuclear DNA can become highly discordant in small populations.     

Fast Accumulation of Nonsynonymous Mutations on the Female-Specific W Chromosome in Birds

Following cessation of recombination during sex chromosome evolution, the nonrecombining sex chromosome is affected by a number of degenerative forces, possibly resulting in the fixation of deleterious mutations. This might take place because of weak selection against recessive or partly recessive deleterious mutations due to permanent heterozygosity of nonrecombining chromosomes. Furthermore, population genetic processes, such as selective sweeps, background selection, and Muller’s ratchet, result in a reduction in N_e, which increase the likelihood of fixation of deleterious mutations. Theory thus predicts that nonrecombining genes should show increased levels of nonsynonymous (d_N) to synonymous substitutions (d_S). We tested this in an avian system by estimating the ratio between d_N and d_S in six gametologous gene pairs located on the Z chromosome and the nonrecombining, female-specific W chromosome. In comparisons, we found a significantly higher d_N/d_S ratio for the W-linked than the Z-linked copy in three of the investigated genes. In a concatenated alignment of all six genes, the d_N/d_S ratio was six times higher for W-linked than Z-linked genes. By using human and mouse as outgroup in maximum likelihood analyses, W-linked genes were found to evolve differently compared with their Z-linked gametologues and outgroup sequences. This seems not to be a consequence of functional diversification because d_N/d_S ratios between gametologous gene copies were consistently low. We conclude that deleterious mutations are accumulating at a high rate on the avian W chromosome, probably as a result of the lack of recombination in this female-specific chromosome. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Deborah Charlesworth]     

Molecular Evolution of Three Avian Neurotrophin Genes: Implications for Proregion Functional Constraints

Neurotrophin proteins are essential for the survival, differentiation, and maintenance of neurons in the peripheral and central nervous systems. Recent studies have shown that the unprocessed proforms of the neurotrophins are preferential high-affinity ligands for p75^NTR and potent inducers of p75^NTR-mediated cell death. Here, we explore differences in the selective constraints acting on the proregions of the three avian neurotrophin genes—NT-3, BDNF, and NGF—in an explicit phylogenetic context. We found a 50-fold difference in levels of constraint as estimated by d _N/d _S ratios, with the NGF proregion showing the lowest degree of constraint and BDNF the highest. These patterns suggest that the high conservation exhibited by the BDNF proregion results from intense functional constraints that are relaxed in NGF and somewhat relaxed in NT-3. The proregion of BDNF is likely to have a function that differentiates it from the corresponding regions of the NGF and NT-3 genes, suggesting that BDNF is the avian neurotrophin most likely to be used both in its precursor and mature forms in vivo.     

Recent relaxation of purifying selection on the tandem-repetitive early-stage histone H3 gene in brooding sea stars

Patterns of nucleotide substitution differ between marine species that have a pelagic feeding (planktotrophic) larval stage and related species that lack such a stage, for both adaptive and non-adaptive reasons. Here, patterns of nucleotide and inferred amino acid substitution are analyzed for the tandem-repetitive early-stage histone H3 gene in 36 sea star species of the order Forcipulatida with documented larval habitat. The relative rate of nonsynonymous substitution (expressed as ω = d_N/d_S) was significantly higher in lineages with a brooded non-feeding (lecithotrophic) larval form than in lineages with a planktotrophic larval form. There was also a significant excess of conservative over radical substitutions. The increase in ω for brooders as compared to non-brooders was much greater than for previously analyzed mitochondrial sequences in echinoderms. These data are consistent with the hypothesis that purifying selection on this gene has been relaxed in brooding lineages compared to non-brooding lineages. The hypotheses of adaptive or neutral evolution are less plausible, although recent pseudogenization following a period of relaxed purifying selection could also explain the results.     

Variance in Female Reproductive Success Differentially Impacts Effective Population Size in the Short-Nosed Fruit Bat, <Emphasis Type="Italic">Cynopterus sphinx</Emphasis>

Effective population size (N _e) quantifies the effects of micro-evolutionary processes and the rate of loss of genetic diversity in a population. Several demographic and mating parameters reduce N _e. Theoretical studies elucidate the impacts of various demographic and mating system parameters on N _e, while empirical studies illustrate realized N _e for species with differing life histories and mating systems. However, effect of intra-specific variation in mating system on effective size remains largely unexplored. In this paper we investigated the effect of promiscuous and polygynous mating on N _e in two wild populations of the short-nosed fruit bat, Cynopterus sphinx. N _e/N (ratio of effective population size to census size) was lower than unity in both populations, and much lower for the polygynous population compared to promiscuous population. Elasticity analyses reveal that N _e/N was sensitive to deviations in the sex ratio. Variance in female reproductive success had a higher impact on N _e compared to variance in male reproductive success in the promiscuous population. However, for the polygynous population, impact of variance in male reproductive success on N _e was higher than that of variance in female reproductive success. Our results suggest that depending on mating system, different populations of the same species could have alternate evolutionary trajectories. The rate of loss of genetic diversity would be lower for the promiscuous population compared to the polygynous population. Our study is the first to highlight which parameters would most significantly impact population specific N _e under different mating systems.     

INCREASED PROBABILITY OF EXTINCTION DUE TO DECREASED GENETIC EFFECTIVE POPULATION SIZE: EXPERIMENTAL POPULATIONS OF CLARKIA PULCHELLA

We established replicated experimental populations of the annual plant Clarkia pulchella to evaluate the existence of a causal relationship between loss of genetic variation and population survival probability. Two treatments differing in the relatedness of the founders, and thus in the genetic effective population size (N_e), were maintained as isolated populations in a natural environment. After three generations, the low N_e treatment had significantly lower germination and survival rates than did the high N_e treatment. These lower germination and survival rates led to decreased mean fitness in the low N_e populations: estimated mean fitness in the low N_e populations was only 21% of the estimated mean fitness in the high N_e populations. This inbreeding depression led to a reduction in population survival: at the conclusion of the experiment, 75% of the high N_e populations were still extant, whereas only 31% of the low N_e populations had survived. Decreased genetic effective population size, which leads to both inbreeding and the loss of alleles by genetic drift, increased the probability of population extinction over that expected from demographic and environmental stochasticity alone. This demonstrates that the genetic effective population size can strongly affect the probability of population persistence.     

Linkage disequilibrium in a finite population that is partially selfing

The linkage disequilibrium expected in a finite, partially selfing population is analyzed, assuming the infinite allele model. Formulas for the expected sum of squares of the linkage disequilibria and the squared standard linkage disequilibrium are derived from the equilibrium values of sixteen inbreeding coefficients required to describe the behavior of the system. These formulas are identical to those obtained with random mating if the effective population size N_e = (1-½S)N and the effective recombination value r_e = (1-S)r/(1-½S), where S is the proportion of selfing, are substituted for the population size and the recombination value. Therefore, the effect of partial selfing at equilibrium is to reduce the population size by a factor 1-½S and the recombination value by a factor (1-S)/(1-½S).     

Genetic Variability and Structuring of Arctic Charr (Salvelinus alpinus) Populations in Northern Fennoscandia

Variation in presumably neutral genetic markers can inform us about evolvability, historical effective population sizes and phylogeographic history of contemporary populations. We studied genetic variability in 15 microsatellite loci in six native landlocked Arctic charr (Salvelinus alpinus) populations in northern Fennoscandia, where this species is considered near threatened. We discovered that all populations were genetically highly (mean F _ST ≈ 0.26) differentiated and isolated from each other. Evidence was found for historical, but not for recent population size bottlenecks. Estimates of contemporary effective population size (N _e) ranged from seven to 228 and were significantly correlated with those of historical N _e but not with lake size. A census size (N _C) was estimated to be approximately 300 individuals in a pond (0.14 ha), which exhibited the smallest N _e (i.e. N _e/N _C = 0.02). Genetic variability in this pond and a connected lake is severely reduced, and both genetic and empirical estimates of migration rates indicate a lack of gene flow between them. Hence, albeit currently thriving, some northern Fennoscandian populations appear to be vulnerable to further loss of genetic variability and are likely to have limited capacity to adapt if selection pressures change.