Isolation with asymmetric gene flow during the nonsynchronous divergence of dry forest birds

Dry forest bird communities in South America are often fragmented by intervening mountains and rainforests, generating high local endemism. The historical assembly of dry forest communities likely results from dynamic processes linked to numerous population histories among codistributed species. Nevertheless, species may diversify in the same way through time if landscape and environmental features, or species ecologies, similarly structure populations. Here we tested whether six co‐distributed taxon pairs that occur in the dry forests of the Tumbes and Marañón Valley of northwestern South America show concordant patterns and modes of diversification. We employed a genome reduction technique, double‐digest restriction site‐associated DNA sequencing, and obtained 4407–7186 genomewide SNPs. We estimated demographic history in each taxon pair and inferred that all pairs had the same best‐fit demographic model: isolation with asymmetric gene flow from the Tumbes into the Marañón Valley, suggesting a common diversification mode. Overall, we also observed congruence in effective population size (N_e) patterns where ancestral N_e were 2.9–11.0× larger than present‐day Marañón Valley populations and 0.3–2.0× larger than Tumbesian populations. Present‐day Marañón Valley N_e was smaller than Tumbes. In contrast, we found simultaneous population isolation due to a single event to be unlikely as taxon pairs diverged over an extended period of time (0.1–2.9 Ma) with multiple nonoverlapping divergence periods. Our results show that even when populations of codistributed species asynchronously diverge, the mode of their differentiation can remain conserved over millions of years. Divergence by allopatric isolation due to barrier formation does not explain the mode of differentiation between these two bird assemblages; rather, migration of individuals occurred before and after geographic isolation.     

Tracing the origins of widespread highland species: a case of Neogene diversification across the Mexican sierras in an endemic lizard

The evolutionary history of the Mexican sierras has been shaped by various geological and climatic events over the past several million years. The relative impacts of these historical events on diversification in highland taxa, however, remain largely uncertain owing to a paucity of studies on broadly‐distributed montane species. We investigated the origins of genetic diversification in widely‐distributed endemic alligator lizards in the genus Barisia to help develop a better understanding of the complex processes structuring biological diversity in the Mexican highlands. We estimated lineage divergence dates and the diversification rate from mitochondrial DNA sequences, and combined divergence dates with reconstructions of ancestral geographical ranges to track lineage diversification across geography through time. Based on our results, we inferred ten geographically structured, well supported mitochondrial lineages within Barisia. Diversification of a widely‐distributed ancestor appears tied to the formation of the Trans‐Mexican Volcanic Belt across central Mexico during the Miocene and Pliocene. The formation of filter barriers such as major river drainages may have later subdivided lineages. The results of the present study provide additional support for the increasing number of studies that suggest Neogene events heavily impacted genetic diversification in widespread montane taxa. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 382–394.     

Diversification in the Mexican horned lizard Phrynosoma orbiculare across a dynamic landscape

The widespread montane Mexican horned lizard Phrynosoma orbiculare (Squamata: Phrynosomatidae) represents an ideal species to investigate the relative impacts of Neogene vicariance and Quaternary climate change on lineage diversification across the Mexican highlands. We used mitochondrial DNA to examine the maternal history of P. orbiculare and estimate the timing and tempo of lineage diversification. Based on our results, we inferred 11 geographically structured, well supported mitochondrial lineages within this species, suggesting P. orbiculare represents a species complex. Six divergences between lineages likely occurred during the Late Miocene and Pliocene, and four splits probably happened during the Pleistocene. Diversification rate appeared relatively constant through time. Spatial and temporal divergences between lineages of P. orbiculare and co-distributed taxa suggest that a distinct period of uplifting of the Transvolcanic Belt around 7.5-3 million years ago broadly impacted diversification in taxa associated with this mountain range. To the north, several river drainages acting as filter barriers differentially subdivided co-distributed highland taxa through time. Diversification patterns observed in P. orbiculare provide additional insight into the mechanisms that impacted differentiation of highland taxa across the complex Mexican highlands.     

Divergent timing and patterns of species accumulation in lowland and highland neotropical birds

Late Pliocene and Pleistocene climatic instability has been invoked to explain the buildup of Neotropical biodiversity, although other theories date Neotropical diversification to earlier periods. If these climatic fluctuations drove Neotropical diversification, then a large proportion of species should date to this period and faunas should exhibit accelerated rates of speciation. However, the unique role of recent climatic fluctuations in promoting diversification could be rejected if late Pliocene and Pleistocene rates declined. To test these temporal predictions, dateable molecular phylogenies for 27 avian taxa were used to contrast the timing and rates of diversification in lowland and highland Neotropical faunas. Trends in diversification rates were analyzed in two ways. First, rates within taxa were analyzed for increasing or decreasing speciation rates through time. There was a significant trend within lowland taxa towards decreasing speciation rates, but no significant trend was observed within most highland taxa. Second, fauna wide diversification rates through time were estimated during one-million-year intervals by combining rates across taxa. In the lowlands, rates were highest during the late Miocene and then decreased towards the present. The decline in rates observed both within taxa and for the fauna as a whole probably resulted from density dependent cladogenesis. In the highlands, faunawide rates did not vary greatly before the Pleistocene but did increase significantly during the last one million years of the Pleistocene following the onset of severe glacial cycles in the Andes. These contrasting patterns of species accumulation suggest that lowland and highland regions were affected differently by recent climatic fluctuations. Evidently, habitat alterations associated with global climate change were not enough to promote an increase in the rate of diversification in lowland faunas. In contrast, direct fragmentation of habitats by glaciers and severe altitudinal migration of montane vegetation zones during climatic cycles may have resulted in the late Pleistocene increase in highland diversification rates. This increase resulted in a fauna with one third of its species dating to the last one million years.     

Biases in Demographic Modeling Affect Our Understanding of Recent Divergence

Testing among competing demographic models of divergence has become an important component of evolutionary research in model and non-model organisms. However, the effect of unaccounted demographic events on model choice and parameter estimation remains largely unexplored. Using extensive simulations, we demonstrate that under realistic divergence scenarios, failure to account for population size (N_e) changes in daughter and ancestral populations leads to strong biases in divergence time estimates as well as model choice. We illustrate these issues reconstructing the recent demographic history of North Sea and Baltic Sea turbots (Scophthalmus maximus) by testing 16 isolation with migration (IM) and 16 secondary contact (SC) scenarios, modeling changes in N_e as well as the effects of linked selection and barrier loci. Failure to account for changes in N_e resulted in selecting SC models with long periods of strict isolation and divergence times preceding the formation of the Baltic Sea. In contrast, models accounting for N_e changes suggest recent (<6 kya) divergence with constant gene flow. We further show how interpreting genomic landscapes of differentiation can help discerning among competing models. For example, in the turbot data, islands of differentiation show signatures of recent selective sweeps, rather than old divergence resisting secondary introgression. The results have broad implications for the study of population divergence by highlighting the potential effects of unmodeled changes in N_e on demographic inference. Tested models should aim at representing realistic divergence scenarios for the target taxa, and extreme caution should always be exercised when interpreting results of demographic modeling.     

LIFE-HISTORY VARIATION AND COMPARATIVE PHYLOGEOGRAPHY OF DARTERS (PISCES: PERCIDAE) FROM THE NORTH AMERICAN CENTRAL HIGHLANDS

Phylogeography seeks to evaluate the relationship between genetic variation and geographic distribution of a species to examine the influence of historical events on divergence. Congruent phylogeographic patterns in codistributed species indicate historical association of the taxa being compared, and the uniform action of biogeographic events in shaping genetic variation. We sought to evaluate the congruence of patterns of genetic variation of five closely related fish species across a well-defined biogeographic boundary. We gathered allozymic and mitochondrial DNA (mtDNA) sequence data for five species of darters (Percidae: Etheostomatini) from populations distributed among biogeographic regions in the Ozark and Ouachita highlands of the south-central United States. Comparisons among species revealed noncongruence in the magnitude of genetic divergence in both allozymes and mtDNA sequences. We hypothesized that noncongruence resulted, in part, from differences in life histories of the species in our comparison. To address this hypothesis, we evaluated the association between gene flow (measured by N_em) and variation in body size and fecundity variables because they have been shown to influence gene flow in fishes. Correlation analysis revealed an association between gene flow and fecundity (r = 0.88), but not with body size (r = 0.36) or reproductive investment per individual (r = –0.23). The result was similar when independent contrasts of the original variables were used in correlation analyses. Phylogeographic analysis of mtDNA sequence data indicated the importance of history, evident in gene trees of Percina nasuta and Percina phoxocephala. Divergence rates between these two taxa may differ because of historically persistent differences in population sizes, reflected in present-day abundance and fecundity differences. Conversely, Percina caprodes showed little evidence of divergence in mtDNA sequences and yielded the highest mean N_em values from allozyme data. Comparisons among closely related, codistributed taxa may help discriminate among the possibilities for noncongruence in biogeographic studies by allowing evaluation of the influence of life history on patterns of gene flow.     

A new Pliocene mollusk fauna from Mejillones, northern Chile

A new Pliocene (3.4 Ma) mollusk fauna from Mejillones Peninsula, northern Chile is described and compared with the Pliocene La Cueva fauna of little constrained age from central Chile and some species from the Huenteguapi Sandstone overlying the Ranquil Formation on Arauco Peninsula, south central Chile. Preliminary correlation is based on faunal similarities. A total of 45 taxa were identified, of which Cyclocardia kieli sp. nov. is new to science. New combinations are Macron escalonia (Vermeij and DeVries, 1997), Austrofusus steinmanni (Möricke, 1896) and Leukoma antiqua (King, 1832). For several species, the oldest occurrences and range extensions are reported. Co-occurrence of warm water taxa, previously assigned to MIS 11, and typical Pliocene taxa on Mejillones cannot be confirmed. Pliocene and Pleistocene mollusk faunas from Mejillones are listed for comparison.     

Repetitive genomic regions and the inference of demographic history

Inference of demographic histories using whole-genome datasets has provided insights into diversification, adaptation, hybridization, and plant–pathogen interactions, and stimulated debate on the impact of anthropogenic interventions and past climate on species demography. However, the impact of repetitive genomic regions on these inferences has mostly been ignored by masking of repeats. We use the Populus trichocarpa genome (Pop_tri_v3) to show that masking of repeat regions leads to lower estimates of effective population size (N_e) in the distant past in contrast to an increase in N_e estimates in recent times. However, in human datasets, masking of repeats resulted in lower estimates of N_e at all time points. We demonstrate that repeats affect demographic inferences using diverse methods like PSMC, MSMC, SMC++, and the Stairway plot. Our genomic analysis revealed that the biases in N_e estimates were dependent on the repeat class type and its abundance in each atomic interval. Notably, we observed a weak, yet consistently significant negative correlation between the repeat abundance of an atomic interval and the N_e estimates for that interval, which potentially reflects the recombination rate variation within the genome. The rationale for the masking of repeats has been that variants identified within these regions are erroneous. We find that polymorphisms in some repeat classes occur in callable regions and reflect reliable coalescence histories (e.g., LTR Gypsy, LTR Copia). The current demography inference methods do not handle repeats explicitly, and hence the effect of individual repeat classes needs careful consideration in comparative analysis. Deciphering the repeat demographic histories might provide a clear understanding of the processes involved in repeat accumulation.Subject terms: Genetic variation, Genome evolution     

Simple life-history traits explain key effective population size ratios across diverse taxa

Effective population size (N_e) controls both the rate of random genetic drift and the effectiveness of selection and migration, but it is difficult to estimate in nature. In particular, for species with overlapping generations, it is easier to estimate the effective number of breeders in one reproductive cycle (N_b) than N_e per generation. We empirically evaluated the relationship between life history and ratios of N_e, N_b and adult census size (N) using a recently developed model (agene) and published vital rates for 63 iteroparous animals and plants. N_b/N_e varied a surprising sixfold across species and, contrary to expectations, N_b was larger than N_e in over half the species. Up to two-thirds of the variance in N_b/N_e and up to half the variance in N_e/N was explained by just two life-history traits (age at maturity and adult lifespan) that have long interested both ecologists and evolutionary biologists. These results provide novel insights into, and demonstrate a close general linkage between, demographic and evolutionary processes across diverse taxa. For the first time, our results also make it possible to interpret rapidly accumulating estimates of N_b in the context of the rich body of evolutionary theory based on N_e per generation.     

Contrasting support for alternative models of genomic variation based on microhabitat preference: species‐specific effects of climate change in alpine sedges

Deterministic processes may uniquely affect codistributed species’ phylogeographic patterns such that discordant genetic variation among taxa is predicted. Yet, explicitly testing expectations of genomic discordance in a statistical framework remains challenging. Here, we construct spatially and temporally dynamic models to investigate the hypothesized effect of microhabitat preferences on the permeability of glaciated regions to gene flow in two closely related montane species. Utilizing environmental niche models from the Last Glacial Maximum and the present to inform demographic models of changes in habitat suitability over time, we evaluate the relative probabilities of two alternative models using approximate Bayesian computation (ABC) in which glaciated regions are either (i) permeable or (ii) a barrier to gene flow. Results based on the fit of the empirical data to data sets simulated using a spatially explicit coalescent under alternative models indicate that genomic data are consistent with predictions about the hypothesized role of microhabitat in generating discordant patterns of genetic variation among the taxa. Specifically, a model in which glaciated areas acted as a barrier was much more probable based on patterns of genomic variation in Carex nova, a wet‐adapted species. However, in the dry‐adapted Carex chalciolepis, the permeable model was more probable, although the difference in the support of the models was small. This work highlights how statistical inferences can be used to distinguish deterministic processes that are expected to result in discordant genomic patterns among species, including species‐specific responses to climate change.     

Comparative Phylogeography of Mississippi Embayment Fishes

The Mississippi Embayment is a prominent physiographic feature of eastern North America consisting of primarily lowland aquatic habitats and a fish fauna that is largely distinct from nearby highland regions. Numerous studies have demonstrated that both pre-Pleistocene and Pleistocene events have had a strong influence on the distributions and relationships of highland fishes in eastern North America. However, the extent to which these same events affected Embayment distributed taxa remains largely unexplored. The purpose of this study was to investigate the relative roles of pre-Pleistocene and Pleistocene events in shaping phylogeographic relationships of four stream dwelling fishes in the Mississippi Embayment. Molecular genetic analyses of the mitochondrial gene cytochrome b were performed for three ictalurid catfish species (Noturus miurus, n = 67; Noturus hildebrandi, n = 93, and Noturus phaeus, n = 44) and one minnow species (Cyprinella camura, n = 78), all distributed in tributary streams of the Mississippi Embayment. Phylogenetic relationships and divergence times among haplotypes for each species were estimated using maximum likelihood and Bayesian methods. Phylogenetic analyses recovered 6 major haplotype clades within N. miurus, 5 within N. hildbrandi, 8 within N. phaeus, and 8 within C. camura. All three Noturus species show a high degree of isolation by drainage, which is less evident in C. camura. A clade of haplotypes from tributaries in the southern portion of the Mississippi Embayment was consistently recovered in all four species. Divergence times among clades spanned the Pleistocene, Pliocene, and Miocene. Novel relationships presented here for C. camura and N. phaeus suggest the potential for cryptic species. Pre-Pleistocene and Pleistocene era sea level fluctuations coincide with some divergence events, but no single event explains any common divergence across all taxa. Like their highland relatives, a combination of both pre-Pleistocene and Pleistocene era events have driven divergences among Embayment lineages.     

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.     

Pleistocene Niche Stability and Lineage Diversification in the Subtropical Spider Araneus omnicolor (Araneidae)

The influence of Quaternary climate oscillations on the diversification of the South American fauna is being increasingly explored. However, most of these studies have focused on taxa that are endemic to tropical environments, and relatively few have treated organisms restricted to subtropical biomes. Here we used an integrative phylogeographical framework to investigate the effects of these climate events on the ecological niche and genetic patterns of the subtropical orb-weaver spider Araneus omnicolor (Araneidae). We analyzed the mitochondrial (Cytochrome Oxidase I, COI) and nuclear (Internal Transcribed Subunit II, ITS2) DNA of 130 individuals throughout the species’ range, and generated distribution models in three different climate scenarios [present, Last Glacial Maximum (LGM), and Last Interglacial Maximum (LIG)]. Additionally, we used an Approximate Bayesian Computation (ABC) approach to compare possible demographic scenarios and select the hypothesis that better explains the genetic patterns of A. omnicolor. We obtained high haplotype diversity but low nucleotide variation among sequences. The population structure and demographic analyses showed discrepancies between markers, suggesting male-biased dispersal in the species. The time-calibrated COI phylogenetic inference showed a recent diversification of lineages (Middle/Late Pleistocene), while the paleoclimate modeling indicated niche stability since ~120 Kya. The ABC results agreed with the niche models, supporting a panmictic population as the most likely historical scenario for the species. These results indicate that A. omnicolor experienced no niche or population reductions during the Late Pleistocene, despite the intense landscape modifications that occurred in the subtropical region, and that other factors beside LGM and LIG climate oscillations might have contributed to the demographic history of this species. This pattern may be related to the high dispersal ability and wide environmental tolerance of A. omnicolor, highlighting the need for more phylogeographical studies with invertebrates and other generalist taxa, in order to understand the effects of Quaternary climate changes on Neotropical biodiversity.     

A temporally dynamic approach for cladistic biogeography and the processes underlying the biogeographic patterns of North American deserts

We implemented a temporally dynamic approach to the cladistic biogeographic analysis of 13 areas of North American deserts and several plant and animal taxa. We undertook a parsimony analysis of paralogy‐free subtrees based on 43 phylogenetic hypotheses of arthropod, vertebrate and plant taxa, assigning their nodes to three different time slices based on their estimated minimum ages: Early‐Mid‐Miocene (23?7 Ma), Late Miocene/Pliocene (6.9?2.5 Ma) and Pleistocene (2.4?0.011 Ma). The analyses resulted in three general area cladograms, one for each time slice, showing different area relationships. They allowed us to detect influences of different geological and palaeoclimatological events of the Early‐Mid‐Miocene, Late Miocene/Pliocene and Pleistocene that might have affected the diversification of the desert biota. Several diversification events in the deserts of North America might have been driven by Neogene uplift, marine incursion and the opening of the California Gulf during the Miocene–Pliocene, whereas climatic fluctuations had the highest impact during the Pleistocene.     

Taking multiple infections of cells and recombination into account leads to small within-host effective-population-size estimates of HIV-1

Whether HIV-1 evolution in infected individuals is dominated by deterministic or stochastic effects remains unclear because current estimates of the effective population size of HIV-1 in vivo, N_e, are widely varying. Models assuming HIV-1 evolution to be neutral estimate N_e∼10²–10⁴, smaller than the inverse mutation rate of HIV-1 (∼10⁵), implying the predominance of stochastic forces. In contrast, a model that includes selection estimates N_e>10⁵, suggesting that deterministic forces would hold sway. The consequent uncertainty in the nature of HIV-1 evolution compromises our ability to describe disease progression and outcomes of therapy. We perform detailed bit-string simulations of viral evolution that consider large genome lengths and incorporate the key evolutionary processes underlying the genomic diversification of HIV-1 in infected individuals, namely, mutation, multiple infections of cells, recombination, selection, and epistatic interactions between multiple loci. Our simulations describe quantitatively the evolution of HIV-1 diversity and divergence in patients. From comparisons of our simulations with patient data, we estimate N_e∼10³–10⁴, implying predominantly stochastic evolution. Interestingly, we find that N_e and the viral generation time are correlated with the disease progression time, presenting a route to a priori prediction of disease progression in patients. Further, we show that the previous estimate of N_e>10⁵ reduces as the frequencies of multiple infections of cells and recombination assumed increase. Our simulations with N_e∼10³–10⁴ may be employed to estimate markers of disease progression and outcomes of therapy that depend on the evolution of viral diversity and divergence.     

Effective/census population size ratio estimation: a compendium and appraisal

With an ecological-evolutionary perspective increasingly applied toward the conservation and management of endangered or exploited species, the genetic estimation of effective population size (N_e) has proliferated. Based on a comprehensive analysis of empirical literature from the past two decades, we asked: (i) how often do studies link N_e to the adult census population size (N)? (ii) To what extent is N_e correctly linked to N? (iii) How readily is uncertainty accounted for in both N_e and N when quantifying N_e/N ratios? and (iv) how frequently and to what degree might errors in the estimation of N_e or N affect inferences of N_e/N ratios? We found that only 20% of available N_e estimates (508 of 2617; 233 studies) explicitly attempted to link N_e and N; of these, only 31% (160 of 508) correctly linked N_e and N. Moreover, only 7% (41 of 508) of N_e/N ratios (correctly linked or not) reported confidence intervals for both N_e and N; for those cases where confidence intervals were reported for N_e only, 31% of N_e/N ratios overlapped with 1, of which more than half also reached below N_e/N = 0.01. Uncertainty in N_e/N ratios thus sometimes spanned at least two orders of magnitude. We conclude that the estimation of N_e/N ratios in natural populations could be significantly improved, discuss several options for doing so, and briefly outline some future research directions.     

Genetic divergence of the Mesoamerican azure‐crowned hummingbird (Amazilia cyanocephala,Trochilidae) across the Motagua‐Polochic‐Jocotán fault system

The cloud forests of Mesoamerica are highly endangered habitats and the existence of narrowly distributed cryptic endemics will increase the number of taxa at potential risk of extinction. Here, we investigate genetic divergence between populations of the azure‐crowned hummingbird (Amazilia cyanocephala), a species complex of endemic hummingbirds to the montane forests of Mesoamerica, by analysing DNA sequences of four mitochondrial markers, morphological data and ecological niche modelling. Our results revealed the presence of two mtDNA lineages corresponding to subspecies A. c. cyanocephala distributed from Tamaulipas to Chiapas in Mexico and Amazilia c. guatemalensis distributed from southern Chiapas to Guatemala. The lineage split can be explained as a consequence of relative isolation of the populations in the different mountain ranges separated by the Motagua‐Polochic‐Jocotán fault system and corresponds to differences in morphology and to the lack of overlap in environmental space between subspecies. The divergence time estimates do not support the proposed model of a highly constrained temporal window at the end of the Pliocene as divergence at this barrier between cyanocephala and guatemalensis and splits of other bird taxa occurred during the Pleistocene.     

Plate tectonics and biogeographical patterns of the Pseudophoxinus (Pisces: Cypriniformes) species complex of central Anatolia, Turkey

We investigated the phylogenetic relationships of Pseudophoxinus (Cyprinidae: Leuciscinae) species from central Anatolia, Turkey to test the hypothesis of geographic speciation driven by early Pliocene orogenic events. We analyzed 1141 aligned base pairs of the complete cytochrome b mitochondrial gene. Phylogenetic relationships reconstructed by maximum likelihood, Bayesian likelihood, and maximum parsimony methods are identical, and generally well supported. Species and clades are restricted to geologically well-defined units, and are deeply divergent from each other. The basal diversification of central Anatolian Pseudophoxinus is estimated to have occurred approximately 15 million years ago. Our results are in agreement with a previous study of the Anatolian fish genus Aphanius that also shows a diversification pattern driven by the Pliocene orogenic events. The distribution of clades of Aphanius and Pseudophoxinus overlap, and areas of distribution comprise the same geological units. The geological history of Anatolia is likely to have had a major impact on the diversification history of many taxa occupying central Anatolia; many of these taxa are likely to be still unrecognized as distinct.     

Effective Population Size Predicts Local Rates but Not Local Mitigation of Read-through Errors

In correctly predicting that selection efficiency is positively correlated with the effective population size (N_e), the nearly neutral theory provides a coherent understanding of between-species variation in numerous genomic parameters, including heritable error (germline mutation) rates. Does the same theory also explain variation in phenotypic error rates and in abundance of error mitigation mechanisms? Translational read-through provides a model to investigate both issues as it is common, mostly nonadaptive, and has good proxy for rate (TAA being the least leaky stop codon) and potential error mitigation via “fail-safe” 3′ additional stop codons (ASCs). Prior theory of translational read-through has suggested that when population sizes are high, weak selection for local mitigation can be effective thus predicting a positive correlation between ASC enrichment and N_e. Contra to prediction, we find that ASC enrichment is not correlated with N_e. ASC enrichment, although highly phylogenetically patchy, is, however, more common both in unicellular species and in genes expressed in unicellular modes in multicellular species. By contrast, N_e does positively correlate with TAA enrichment. These results imply that local phenotypic error rates, not local mitigation rates, are consistent with a drift barrier/nearly neutral model.     

Flight loss linked to faster molecular evolution in insects

The loss of flight ability has occurred thousands of times independently during insect evolution. Flight loss may be linked to higher molecular evolutionary rates because of reductions in effective population sizes (N_e) and relaxed selective constraints. Reduced dispersal ability increases population subdivision, may decrease geographical range size and increases (sub)population extinction risk, thus leading to an expected reduction in N_e. Additionally, flight loss in birds has been linked to higher molecular rates of energy-related genes, probably owing to relaxed selective constraints on energy metabolism. We tested for an association between insect flight loss and molecular rates through comparative analysis in 49 phylogenetically independent transitions spanning multiple taxa, including moths, flies, beetles, mayflies, stick insects, stoneflies, scorpionflies and caddisflies, using available nuclear and mitochondrial protein-coding DNA sequences. We estimated the rate of molecular evolution of flightless (FL) and related flight-capable lineages by ratios of non-synonymous-to-synonymous substitutions (dN/dS) and overall substitution rates (OSRs). Across multiple instances of flight loss, we show a significant pattern of higher dN/dS ratios and OSRs in FL lineages in mitochondrial but not nuclear genes. These patterns may be explained by relaxed selective constraints in FL ectotherms relating to energy metabolism, possibly in combination with reduced N_e.