Rapid evolution in response to increased temperature maintains population viability despite genetic erosion in a tropical ectotherm

Climate change is predicted to increase the average global air temperature by up to 4.0 °C by the end of the century. This increased temperature could have negative effects on many life history traits that are closely linked to fitness. Many species will therefore have to adapt to the warmer environment, but life history traits often have limited additive genetic variance. Here, we investigated population demographics and the evolutionary response of life history traits, as well as genetic diversity in guppies (Poecilia reticulata), in response to an experimentally increased temperature. There were fewer successful pregnancies, smaller brood sizes, and males matured earlier at a higher temperature as compared to control populations. However, there was no sign of an evolutionary response in these traits after 24 months of exposure to the increased temperature. We also found that population size, brood survivorship, sex ratio, and male length at maturity were unaffected by the increased temperature. Genetic diversity decreased rapidly in the increased temperature populations at a rate equivalent to an effective population size of only one quarter of the controls, revealing a clear signature of selection in response to increased temperature. This genetic erosion, however, could hamper the adaptive potential of the populations to other environmental changes associated with climate change.     

Temporal genetic stability and high effective population size despite fisheries-induced life-history trait evolution in the North Sea sole

Heavy fishing and other anthropogenic influences can have profound impact on a species' resilience to harvesting. Besides the decrease in the census and effective population size, strong declines in mature adults and recruiting individuals may lead to almost irreversible genetic changes in life-history traits. Here, we investigated the evolution of genetic diversity and effective population size in the heavily exploited sole (Solea solea), through the analysis of historical DNA from a collection of 1379 sole otoliths dating back from 1957. Despite documented shifts in life-history traits, neutral genetic diversity inferred from 11 microsatellite markers showed a remarkable stability over a period of 50 years of heavy fishing. Using simulations and corrections for fisheries induced demographic variation, both single-sample estimates and temporal estimates of effective population size (N(e) ) were always higher than 1000, suggesting that despite the severe census size decrease over a 50-year period of harvesting, genetic drift is probably not strong enough to significantly decrease the neutral diversity of this species in the North Sea. However, the inferred ratio of effective population size to the census size (N(e) /N(c) ) appears very small (10(-5) ), suggesting that overall only a low proportion of adults contribute to the next generation. The high N(e) level together with the low N(e) /N(c) ratio is probably caused by a combination of an equalized reproductive output of younger cohorts, a decrease in generation time and a large variance in reproductive success typical for marine species. Because strong evolutionary changes in age and size at first maturation have been observed for sole, changes in adaptive genetic variation should be further monitored to detect the evolutionary consequences of human-induced selection.     

Dissecting the pleiotropic consequences of a quantitative trait nucleotide

The downstream consequences of a single quantitative trait polymorphism can provide important insight into the molecular basis of a trait. However, the molecular consequences of a polymorphism may be complex and only a subset of these may influence the trait of interest. In natural isolates of Saccharomyces cerevisiae , a nonsynonymous polymorphism in cystathione β-synthase ( CYS4 ) causes a deficiency in both cysteine and glutathione that results in rust-colored colonies and drug-dependent growth defects. Using a single-nucleotide allele replacement, we characterized the effects of this polymorphism on gene expression levels across the genome. To determine whether any of the differentially expressed genes are necessary for the production of rust-colored colonies, we screened the yeast deletion collection for genes that enhance or suppress rust coloration. We found that genes in the sulfur assimilation pathway are required for the production of rust color but not the drug-sensitivity phenotype. Our results show that a single quantitative trait polymorphism can generate a complex set of downstream changes, providing a molecular basis for pleiotropy.     

Mature habitats associated with genetic divergence despite strong dispersal ability in an arthropod

Background  

Populations may be bound by contemporary gene flow, selective sweeps, and extinction-recolonization processes. Indeed, existing molecular estimates indicate that species with low levels of gene flow are rare. However, strong priority effects and local selective regimes may hinder gene flow (despite dispersal) sending populations on independent evolutionary trajectories. In this scenario (the monopolization hypothesis), population differentiation will increase with time and genealogical evidence should yield ample private haplotypes. Cyclical parthenogens (e.g. rotifers and cladocerans such as Daphnia) have an increased capacity for rapid local adaptation and priority effects because sexual reproduction is followed by multiple generations of clonal selection and massive egg bank formation. We aimed to better understand the history of population differentiation and ongoing gene flow in Daphnia rosea s.l., by comparing population and regional divergences in mature unglaciated areas and younger previously glaciated areas. We also examined the timing and paths of colonization of previously-glaciated areas to assess the dispersal limitations of D. rosea s.l. We used DNA sequence variation (84 populations and >400 individuals) at the mitochondrial ND2 and nuclear HSP90 loci from Holarctic populations for our genetic analyses.     

Shared quantitative trait loci underlying the genetic correlation between continuous traits

We review genetic correlations among quantitative traits in light of their underlying quantitative trait loci (QTL). We derive an expectation of genetic correlation from the effects of underlying loci and test whether published genetic correlations can be explained by the QTL underlying the traits. While genetically correlated traits shared more QTL (33%) on average than uncorrelated traits (11%), the actual number of shared QTL shared was small. QTL usually predicted the sign of the correlation with good accuracy, but the quantitative prediction was poor. Approximately 25% of trait pairs in the data set had at least one QTL with antagonistic effects. Yet a significant minority (20%) of such trait pairs have net positive genetic correlations due to such antagonistic QTL 'hidden' within positive genetic correlations. We review the evidence on whether shared QTL represent single pleiotropic loci or closely linked monotropic genes, and argue that strict pleiotropy can be viewed as one end of a continuum of recombination rates where r=0. QTL studies of genetic correlation will likely be insufficient to predict evolutionary trajectories over long time spans in large panmictic populations, but will provide important insights into the trade-offs involved in population and species divergence.     

Constrained evolution of a quantitative character by pleiotropic mutation

The long-term response to directional selection and its selection limit are derived for a quantitative character that is controlled by pleiotropic mutations with direct deleterious effect on fitness. Directional selection is assumed to be weaker than the selection acting directly on mutations via deleterious effects (purging selection), which renders all mutations to eventual elimination. The analysis embedding this restrictive assumption indicates that the evolutionary response of the character starting from an equilibrium state, in which mutation and purging selection balance but no directional selection is operating, decreases monotonically with time at an exponential rate. And the fading rate of responses is mostly determined by the direct deleterious effect. Contrary to the expectation by the standard selection limit theory based on fixation of extant genetic variation, the present model predicts that the selection limit depends on the intensity of directional selection, the limit being proportional to the ratio of the directional selection intensity to the direct deleterious effect. A slightly larger genetic variance is maintained at the selection limit than would be without directional selection.     

A bivariate quantitative genetic model for a threshold trait and a survival trait

Many of the functional traits considered in animal breeding can be analyzed as threshold traits or survival traits with examples including disease traits, conformation scores, calving difficulty and longevity. In this paper we derive and implement a bivariate quantitative genetic model for a threshold character and a survival trait that are genetically and environmentally correlated. For the survival trait, we considered the Weibull log-normal animal frailty model. A Bayesian approach using Gibbs sampling was adopted in which model parameters were augmented with unobserved liabilities associated with the threshold trait. The fully conditional posterior distributions associated with parameters of the threshold trait reduced to well known distributions. For the survival trait the two baseline Weibull parameters were updated jointly by a Metropolis-Hastings step. The remaining model parameters with non-normalized fully conditional distributions were updated univariately using adaptive rejection sampling. The Gibbs sampler was tested in a simulation study and illustrated in a joint analysis of calving difficulty and longevity of dairy cattle. The simulation study showed that the estimated marginal posterior distributions covered well and placed high density to the true values used in the simulation of data. The data analysis of calving difficulty and longevity showed that genetic variation exists for both traits. The additive genetic correlation was moderately favorable with marginal posterior mean equal to 0.37 and 95% central posterior credibility interval ranging between 0.11 and 0.61. Therefore, this study suggests that selection for improving one of the two traits will be beneficial for the other trait as well.     

Rapid evolution of body fluid regulation following independent invasions into freshwater habitats

Colonizations from marine to freshwater environments constitute among the most dramatic evolutionary transitions in the history of life. Colonizing dilute environments poses great challenges for acquiring essential ions against steep concentration gradients. This study explored the evolution of body fluid regulation following freshwater invasions by the copepod Eurytemora affinis. The goals of this study were to determine (1) whether invasions from saline to freshwater habitats were accompanied by evolutionary shifts in body fluid regulation (hemolymph osmolality) and (2) whether parallel shifts occurred during independent invasions. We measured hemolymph osmolality for ancestral saline and freshwater invading populations reared across a range of common-garden salinities (0.2-25 PSU). Our results revealed the evolution of increased hemolymph osmolality (by 16-31%) at lower salinities in freshwater populations of E. affinis relative to their saline ancestors. Moreover, we observed the same evolutionary shifts across two independent freshwater invasions. Such increases in hemolymph osmolality are consistent with evidence of increased ion uptake in freshwater populations at low salinity, found in a previous study, and are likely to entail increased energetic costs upon invading freshwater habitats. Our findings are consistent with the evolution of increased physiological regulation accompanying transitions into stressful environments.     

Successful biological invasion despite a severe genetic load

Understanding the factors that influence the success of ecologically and economically damaging biological invasions is of prime importance. Recent studies have shown that invasive populations typically exhibit minimal, if any, reductions in genetic diversity, suggesting that large founding populations and/or multiple introductions are required for the success of biological invasions, consistent with predictions of the propagule pressure hypothesis. Through population genetic analysis of neutral microsatellite markers and a gene experiencing balancing selection, we demonstrate that the solitary bee Lasioglossum leucozonium experienced a single and severe bottleneck during its introduction from Europe. Paradoxically, the success of L. leucozonium in its introduced range occurred despite the severe genetic load caused by single-locus complementary sex-determination that still turns 30% of female-destined eggs into sterile diploid males, thereby substantially limiting the growth potential of the introduced population. Using stochastic modeling, we show that L. leucozonium invaded North America through the introduction of a very small number of propagules, most likely a singly-mated female. Our results suggest that chance events and ecological traits of invaders are more important than propagule pressure in determining invasion success, and that the vigilance required to prevent invasions may be considerably greater than has been previously considered.     

Derivation of a Bayes factor to distinguish between linked or pleiotropic quantitative trait loci

A simple procedure to calculate the Bayes factor between linked and pleiotropic QTL models is presented. The Bayes factor is calculated from the marginal prior and posterior densities of the locations of the QTL under a linkage and a pleiotropy model. The procedure is computed with a Gibbs sampler, and it can be easily applied to any model including the location of the QTL as a variable. The procedure was compared with a multivariate least-squares method. The proposed procedure showed better results in terms of power of detection of linkage when low information is available. As information increases, the performance of both procedures becomes similar. An example using data provided by an Iberian by Landrace pig intercross is presented. The results showed that three different QTL segregate in SSC6: a pleiotropic QTL affects myristic, palmitic, and eicosadienoic fatty acids; another pleiotropic QTL affects palmitoleic, stearic, and vaccenic fatty acids; and a third QTL affects the percentage of linoleic acid. In the example, the Bayes factor approach was more powerful than the multivariate least-squares approach.     

Pathogen-induced rapid evolution in a vertebrate life-history trait

Anthropogenic factors, including climate warming, are increasing the incidence and prevalence of infectious diseases worldwide. Infectious diseases caused by pathogenic parasites can have severe impacts on host survival, thereby altering the selection regime and inducing evolutionary responses in their hosts. Knowledge about such evolutionary consequences in natural populations is critical to mitigate potential ecological and economic effects. However, studies on pathogen-induced trait changes are scarce and the pace of evolutionary change is largely unknown, particularly in vertebrates. Here, we use a time series from long-term monitoring of perch to estimate temporal trends in the maturation schedule before and after a severe pathogen outbreak. We show that the disease induced a phenotypic change from a previously increasing to a decreasing size at maturation, the most important life-history transition in animals. Evolutionary rates imposed by the pathogen were high and comparable to those reported for populations exposed to intense human harvesting. Pathogens thus represent highly potent drivers of adaptive phenotypic evolution in vertebrates.     

Complex genetic interactions in a quantitative trait locus

Whether in natural populations or between two unrelated members of a species, most phenotypic variation is quantitative. To analyze such quantitative traits, one must first map the underlying quantitative trait loci. Next, and far more difficult, one must identify the quantitative trait genes (QTGs), characterize QTG interactions, and identify the phenotypically relevant polymorphisms to determine how QTGs contribute to phenotype. In this work, we analyzed three Saccharomyces cerevisiae high-temperature growth (Htg) QTGs (MKT1, END3, and RHO2). We observed a high level of genetic interactions among QTGs and strain background. Interestingly, while the MKT1 and END3 coding polymorphisms contribute to phenotype, it is the RHO2 3′UTR polymorphisms that are phenotypically relevant. Reciprocal hemizygosity analysis of the Htg QTGs in hybrids between S288c and ten unrelated S. cerevisiae strains reveals that the contributions of the Htg QTGs are not conserved in nine other hybrids, which has implications for QTG identification by marker-trait association. Our findings demonstrate the variety and complexity of QTG contributions to phenotype, the impact of genetic background, and the value of quantitative genetic studies in S. cerevisiae.     

Phenotypic variation and natural selection at catsup, a pleiotropic quantitative trait gene in Drosophila

Quantitative traits are shaped by networks of pleiotropic genes . To understand the mechanisms that maintain genetic variation for quantitative traits in natural populations and to predict responses to artificial and natural selection, we must evaluate pleiotropic effects of underlying quantitative trait genes and define functional allelic variation at the level of quantitative trait nucleotides (QTNs). Catecholamines up (Catsup), which encodes a negative regulator of tyrosine hydroxylase , the rate-limiting step in the synthesis of the neurotransmitter dopamine, is a pleiotropic quantitative trait gene in Drosophila melanogaster. We used association mapping to determine whether the same or different QTNs at Catsup are associated with naturally occurring variation in multiple quantitative traits. We sequenced 169 Catsup alleles from a single population and detected 33 polymorphisms with little linkage disequilibrium (LD). Different molecular polymorphisms in Catsup are independently associated with variation in longevity, locomotor behavior, and sensory bristle number. Most of these polymorphisms are potentially functional variants in protein coding regions, have large effects, and are not common. Thus, Catsup is a pleiotropic quantitative trait gene, but individual QTNs do not have pleiotropic effects. Molecular population genetic analyses of Catsup sequences are consistent with balancing selection maintaining multiple functional polymorphisms.     

A statistical model for testing the pleiotropic control of phenotypic plasticity for a count trait

The differences of a phenotypic trait produced by a genotype in response to changes in the environment are referred to as phenotypic plasticity. Despite its importance in the maintenance of genetic diversity via genotype-by-environment interactions, little is known about the detailed genetic architecture of this phenomenon, thus limiting our ability to predict the pattern and process of microevolutionary responses to changing environments. In this article, we develop a statistical model for mapping quantitative trait loci (QTL) that control the phenotypic plasticity of a complex trait through differentiated expressions of pleiotropic QTL in different environments. In particular, our model focuses on count traits that represent an important aspect of biological systems, controlled by a network of multiple genes and environmental factors. The model was derived within a multivariate mixture model framework in which QTL genotype-specific mixture components are modeled by a multivariate Poisson distribution for a count trait expressed in multiple clonal replicates. A two-stage hierarchic EM algorithm is implemented to obtain the maximum-likelihood estimates of the Poisson parameters that specify environment-specific genetic effects of a QTL and residual errors. By approximating the number of sylleptic branches on the main stems of poplar hybrids by a Poisson distribution, the new model was applied to map QTL that contribute to the phenotypic plasticity of a count trait. The statistical behavior of the model and its utilization were investigated through simulation studies that mimic the poplar example used. This model will provide insights into how genomes and environments interact to determine the phenotypes of complex count traits.     

Rapid genetic diversification and high fitness penalties associated with pathogenicity evolution in a plant virus

Under the gene-for-gene model of host-pathogen coevolution, recognition of pathogen avirulence factors by host resistance factors triggers host defenses and limits infection. Theory predicts that the evolution of higher levels of pathogenicity will be associated with fitness penalties and that the cost of higher pathogenicity must be much smaller than that of not infecting the host. The analysis of pathogenicity costs is of academic and applied relevance, as these are determinants for the success of resistance genes bred into crops for disease control. However, most previous attempts of addressing this issue in plant pathogens yielded conflicting and inconclusive results. We have analyzed the costs of pathogenicity in pepper-infecting tobamoviruses defined by their ability to infect pepper plants with different alleles at the resistance locus L. We provide conclusive evidence of pathogenicity-associated costs by comparison of pathotype frequency with the fraction of the crop carrying the various resistance alleles, by timescaled phylogenies, and by temporal analyses of population dynamics and selection pressures using nucleotide sequences. In addition, experimental estimates of relative fitness under controlled conditions also provided evidence of high pathogenicity costs. These high pathogenicity costs may reflect intrinsic properties of plant virus genomes and should be considered in future models of host-parasite coevolution.     

Global isolation by distance despite strong regional phylogeography in a small metazoan

Background  

Small vagile eukaryotic organisms, which comprise a large proportion of the Earth's biodiversity, have traditionally been thought to lack the extent of population structuring and geographic speciation observed in larger taxa. Here we investigate the patterns of genetic diversity, amongst populations of the salt lake microscopic metazoan Brachionus plicatilis s. s. (sensu stricto) (Rotifera: Monogononta) on a global scale. We examine the phylogenetic relationships of geographic isolates from four continents using a 603 bp fragment of the mitochondrial COI gene to investigate patterns of phylogeographic subdivision in this species. In addition we investigate the relationship between genetic and geographic distances on a global scale to try and reconcile the paradox between the high vagility of this species and the previously reported patterns of restricted gene flow, even over local spatial scales.     

Birth-death symmetry in the evolution of a social trait

Studies of the evolution of a social trait often make ecological assumptions (of population structure, life history), and thus a trait can be studied many different times with different assumptions. Here, I consider a Moran model of continuous reproduction and use an inclusive fitness analysis to investigate the relationships between fecundity or survival selection and birth-death (BD) or death-birth (DB) demography on the evolution of a social trait. A simple symmetry obtains: fecundity (respectively survival) effects under BD behave the same as survival (respectively fecundity) effects under DB. When these results are specialized to a homogeneous population, greatly simplified conditions for a positive inclusive fitness effect are obtained in both a finite and an infinite population. The results are established using the elegant formalism of mathematical group theory and are illustrated with an example of a finite population arranged in a cycle with asymmetric offspring dispersal.     

Successful maintenance of a stingless bee population despite a severe genetic bottleneck

Stingless bees play an important ecological role as pollinators of many wild plant species in the tropics and have significant potential for the pollination of agricultural crops. Nevertheless, conservation efforts as well as commercial breeding programmes require better guidelines on the amount of genetic variation that is needed to maintain viable populations. In this context, we carried out a long-term genetic study on the stingless bee Melipona scutellaris to evaluate the population viability consequences of prolonged breeding from a small number of founder colonies. In particular, it was artificially imposed a genetic bottleneck by setting up a population starting from only two founder colonies, and continued breeding from it for a period of over 10?years in a location outside its natural area of occurrence. We show that despite a great reduction in the number of alleles present at both neutral microsatellite loci and the sex-determining locus relative to its natural source population, and an increased frequency in the production of sterile diploid males, the genetically impoverished population could be successfully bred and maintained for at least 10?years. This shows that in stingless bees, breeding from a small stock of colonies may have less severe consequences than previously suspected. In addition, we provide a simulation model to determine the number of colonies that are needed to maintain a certain number of sex alleles in a population, thereby providing useful guidelines for stingless bee breeding and conservation efforts.     

Rapid molecular evolution in a living fossil

The tuatara of New Zealand is a unique reptile that coexisted with dinosaurs and has changed little morphologically from its Cretaceous relatives. Tuatara have very slow metabolic and growth rates, long generation times and slow rates of reproduction. This suggests that the species is likely to exhibit a very slow rate of molecular evolution. Our analysis of ancient and modern tuatara DNA shows that, surprisingly, tuatara have the highest rate of molecular change recorded in vertebrates. Our work also suggests that rates of neutral molecular and phenotypic evolution are decoupled.