Natural selection on fecundity variance in subdivided populations: kin selection meets bet hedging.

In a series of seminal articles in 1974, 1975, and 1977, J. H. Gillespie challenged the notion that the "fittest" individuals are those that produce on average the highest number of offspring. He showed that in small populations, the variance in fecundity can determine fitness as much as mean fecundity. One likely reason why Gillespie's concept of within-generation bet hedging has been largely ignored is the general consensus that natural populations are of large size. As a consequence, essentially no work has investigated the role of the fecundity variance on the evolutionary stable state of life-history strategies. While typically large, natural populations also tend to be subdivided in local demes connected by migration. Here, we integrate Gillespie's measure of selection for within-generation bet hedging into the inclusive fitness and game theoretic measure of selection for structured populations. The resulting framework demonstrates that selection against high variance in offspring number is a potent force in large, but structured populations. More generally, the results highlight that variance in offspring number will directly affect various life-history strategies, especially those involving kin interaction. The selective pressures on three key traits are directly investigated here, namely within-generation bet hedging, helping behaviors, and the evolutionary stable dispersal rate. The evolutionary dynamics of all three traits are markedly affected by variance in offspring number, although to a different extent and under different demographic conditions.

FTO genotype is associated with phenotypic variability of body mass index.

There is evidence across several species for genetic control of phenotypic variation of complex traits, such that the variance among phenotypes is genotype dependent. Understanding genetic control of variability is important in evolutionary biology, agricultural selection programmes and human medicine, yet for complex traits, no individual genetic variants associated with variance, as opposed to the mean, have been identified. Here we perform a meta-analysis of genome-wide association studies of phenotypic variation using ∼170,000 samples on height and body mass index (BMI) in human populations. We report evidence that the single nucleotide polymorphism (SNP) rs7202116 at the FTO gene locus, which is known to be associated with obesity (as measured by mean BMI for each rs7202116 genotype), is also associated with phenotypic variability. We show that the results are not due to scale effects or other artefacts, and find no other experiment-wise significant evidence for effects on variability, either at loci other than FTO for BMI or at any locus for height. The difference in variance for BMI among individuals with opposite homozygous genotypes at the FTO locus is approximately 7%, corresponding to a difference of ∼0.5 kilograms in the standard deviation of weight. Our results indicate that genetic variants can be discovered that are associated with variability, and that between-person variability in obesity can partly be explained by the genotype at the FTO locus. The results are consistent with reported FTO by environment interactions for BMI, possibly mediated by DNA methylation. Our BMI results for other SNPs and our height results for all SNPs suggest that most genetic variants, including those that influence mean height or mean BMI, are not associated with phenotypic variance, or that their effects on variability are too small to detect even with samples sizes greater than 100,000.

On stand by: host genetics of HIV control.

The impact of host genetic variation on determining the differential outcomes after HIV infection has been studied by two approaches: targeting of candidate genes and genome-wide association studies (GWASs). The overlap in genetic variants that has been identified by these two means has essentially been restricted to variants near to the human leukocyte antigen (HLA) class I genes, although variation in the CCR5 locus, which was first shown to have an effect on HIV outcomes using the candidate gene approach, does reach significance genome-wide when very large samples sizes (i.e. thousands) are used in GWAS. Overall, many of the variants identified by the candidate gene approach are likely to be spurious, as no additional variants apart from a novel variant near the HLA-C gene have been consistently identified by GWAS. Variants with low frequency and/or low impact on HIV outcomes are likely to exist in the genome and there could be many of them, but these are not identifiable, given current GWAS sample sizes. Several loci centrally involved in the immune response, including the immunoglobulin genes, T-cell receptor loci, or leukocyte receptor complex, are either poorly covered on the GWAS chips or difficult to interpret due to their repetitive nature and/or the presence of insertion/deletion polymorphisms in the region. These loci warrant further interrogation, but genetic characterization of these regions across a range of individuals will first be required. Finally, synergistic interactions between loci may affect outcome after infection, as suggested by associations of specific, functionally relevant HLA and killer cell immunoglobulin-like receptor variants with HIV disease outcomes, and these require further consideration as well.

Common variants near MC4R are associated with fat mass, weight and risk of obesity.

To identify common variants influencing body mass index (BMI), we analyzed genome-wide association data from 16,876 individuals of European descent. After previously reported variants in FTO, the strongest association signal (rs17782313, P = 2.9 x 10(-6)) mapped 188 kb downstream of MC4R (melanocortin-4 receptor), mutations of which are the leading cause of monogenic severe childhood-onset obesity. We confirmed the BMI association in 60,352 adults (per-allele effect = 0.05 Z-score units; P = 2.8 x 10(-15)) and 5,988 children aged 7-11 (0.13 Z-score units; P = 1.5 x 10(-8)). In case-control analyses (n = 10,583), the odds for severe childhood obesity reached 1.30 (P = 8.0 x 10(-11)). Furthermore, we observed overtransmission of the risk allele to obese offspring in 660 families (P (pedigree disequilibrium test average; PDT-avg) = 2.4 x 10(-4)). The SNP location and patterns of phenotypic associations are consistent with effects mediated through altered MC4R function. Our findings establish that common variants near MC4R influence fat mass, weight and obesity risk at the population level and reinforce the need for large-scale data integration to identify variants influencing continuous biomedical traits.

Phylogeography and Pleistocene refugia of the adder (Vipera berus) as inferred from mitochondrial DNA sequence data.

In order to contribute to the debate about southern glacial refugia used by temperate species and more northern refugia used by boreal or cold-temperate species, we examined the phylogeography of a widespread snake species (Vipera berus) inhabiting Europe up to the Arctic Circle. The analysis of the mitochondrial DNA (mtDNA) sequence variation in 1043 bp of the cytochrome b gene and in 918 bp of the noncoding control region was performed with phylogenetic approaches. Our results suggest that both the duplicated control region and cytochrome b evolve at a similar rate in this species. Phylogenetic analysis showed that V. berus is divided into three major mitochondrial lineages, probably resulting from an Italian, a Balkan and a Northern (from France to Russia) refugial area in Eastern Europe, near the Carpathian Mountains. In addition, the Northern clade presents an important substructure, suggesting two sequential colonization events in Europe. First, the continent was colonized from the three main refugial areas mentioned above during the Lower-Mid Pleistocene. Second, recolonization of most of Europe most likely originated from several refugia located outside of the Mediterranean peninsulas (Carpathian region, east of the Carpathians, France and possibly Hungary) during the Mid-Late Pleistocene, while populations within the Italian and Balkan Peninsulas fluctuated only slightly in distribution range, with larger lowland populations during glacial times and with refugial mountain populations during interglacials, as in the present time. The phylogeographical structure revealed in our study suggests complex recolonization dynamics of the European continent by V. berus, characterized by latitudinal as well as altitudinal range shifts, driven by both climatic changes and competition with related species.

Comparative phylogeography and postglacial colonization routes in Europe.

The Quaternary cold periods in Europe are thought to have heavily influenced the amount and distribution of intraspecific genetic variation in both animals and plants. The phylogeographies of 10 taxa, including mammals (Ursus arctos, Sorex spp., Crocidura suaveolens, Arvicola spp.), amphibians (Triturus spp.), arthropods (Chorthippus parallelus), and plants (Abies alba, Picea abies, Fagus sylvatica, Quercus spp.), were analysed to elucidate general trends across Europe. Only a small degree of congruence was found amongst the phylogeographies of the 10 taxa, but the likely postglacial colonization routes exhibit some similarities. A Brooks parsimony analysis produced an unrooted area phylogram, showing that: (i) the northern regions were colonized generally from the Iberic and Balkanic refugia; and (ii) the Italian lineages were often isolated due to the presence of the Alpine barrier. The comparison of colonization routes highlighted four main suture-zones where lineages from the different refugia meet. Some of the intraspecific genetic distances among lineages indicated a prequaternary divergence that cannot be connected to any particular cold period, but are probably related mainly to the date of arrival of each taxon in the European continent. As a consequence, molecular genetics so far appears to be of limited use in dating Quaternary events.

The vitamin D receptor gene bAt (CCA) haplotype impairs the response to pegylated-interferon/ribavirin-based therapy in chronic hepatitis C patients.

BACKGROUND: Chronic hepatitis C infection is a major cause of end-stage liver disease. Therapy outcome is influenced by 25-OH vitamin D deficiency. To further address this observation, our study investigates the impact of the vitamin D receptor (NR1I1) haplotype and combined effects of plasma vitamin D levels in a well-described cohort of hepatitis C patients. METHODS: A total of 155 chronic hepatitis C patients were recruited from the Swiss Hepatitis C Cohort Study for NR1I1 genotyping and plasma 25-OH vitamin D level measurement. NR1I1 genotype data and combined effects of plasma 25-OH vitamin D level were analysed regarding therapy response (sustained virological response). RESULTS: A strong association was observed between therapy non-response and the NR1I1 CCA (bAt) haplotype consisting of rs1544410 (BsmI) C, rs7975232 (ApaI) C and rs731236 (TaqI) A alleles. Of the HCV patients carrying the CCA haplotype, 50.3% were non-responders (odds ratio [OR] 1.69, 95% CI 1.07, 2.67; P=0.028). A similar association was observed for the combinational CCCCAA genotype (OR 2.94, 95% CI 1.36, 6.37; P=0.007). The combinational CCCCAA genotype was confirmed as an independent risk factor for non-response in multivariate analysis (OR 2.50, 95% CI 1.07, 5.87; P=0.034). Analysing combined effects, a significant impact of low 25-OH vitamin D levels on sustained virological response were only seen in patients with the unfavourable NR1I1 CCA (bAt) haplotype (OR for non-SVR 3.55; 95% CI 1.005, 12.57; P=0.049). CONCLUSIONS: NR1I1 vitamin D receptor polymorphisms influence response to pegylated-interferon/ribavirin-based therapy in chronic hepatitis C and exert an additive genetic predisposition to previously described low 25-OH vitamin D serum levels.

Host genome influences on susceptibility to HIV-1

In vitro and in vivo analyses identified a significant component of heritability in cellular or host susceptibility to HIV-1. The bases for susceptibility can be traced to genetic differences (inter-species) resulting from evolutionary adaptation to exogenous (and endogenous) retroviral infections, and to intra-species and inter-individual (human) differences associated with genetic variation. We have completed large scale evolutionary analysis of genes involved in HIV life cycle and pathogenesis, as well as participating and conducting genome-wide association studies, linkage analysis, and transcriptome analysis. These studies allowed a better understanding of the influence of common human variants in HIV-1 susceptibility and define a number of experimental challenges in the filed: understanding of the role of rare and private mutations in susceptibility, and the development of better tools for the integration of data from large-scale studies.

Melanin-based colorations signal strategies to cope with poor and rich environments

One hypothesis for the maintenance of genetic variation states that alternative genotypes are adapted to different environmental conditions (i.e., genotype-by-environment interaction GxE) that vary in space and time. Although GxE has been demonstrated for morphological traits, little evidence has been given whether these GxE are associated with traits used as signal in mate choice. In three wild bird species, we investigated whether the degree of melanin-based coloration, a heritable trait, covaries with nestling growth rate in rich and poor environments. Variation in the degree of reddish-brown phaeomelanism is pronounced in the barn owl (Tyto alba) and tawny owl (Strix aluco), and variation in black eumelanism in the barn owl and Alpine swift (Apus melba). Melanin-based coloration has been shown to be a criterion in mate choice in the barn owl. We cross-fostered hatchlings to test whether nestlings sired by parents displaying melanin-based colorations to different extent exhibit alternative growth trajectories when raised by foster parents in poor (experimentally enlarged broods) and rich (experimentally reduced broods) environments. With respect to phaeomelanism, barn owl and tawny owl offspring sired by redder parents grew more rapidly in body mass only in experimentally reduced broods. With respect to eumelanism, Alpine swift offspring of darker fathers grew their wings more rapidly only in experimentally enlarged broods, a difference that was not detected in reduced broods. These interactions between parental melanism and offspring growth rate indicate that individuals display substantial plasticity in response to the rearing environment which is associated with the degree of melanism: at least with respect to nestling growth, phaeomelanic and eumelanic individuals are best adapted to rich and poor environments, respectively. It now remains to be investigated why eumelanism and phaeomelanism have a different signaling function and what the lifelong consequences of these melanism-dependent allocation strategies are. This is important to fully appraise the role played by environmental heterogeneity in maintaining variation in the degree of melanin-based coloration.

The co-evolution of culturally inherited altruistic helping and cultural transmission under random group formation.

Limited migration results in kin selective pressure on helping behaviors under a wide range of ecological, demographic and life-history situations. However, such genetically determined altruistic helping can evolve only when migration is not too strong and group size is not too large. Cultural inheritance of helping behaviors may allow altruistic helping to evolve in groups of larger size because cultural transmission has the potential to markedly decrease the variance within groups and augment the variance between groups. Here, we study the co-evolution of culturally inherited altruistic helping behaviors and two alternative cultural transmission rules for such behaviors. We find that conformist transmission, where individuals within groups tend to copy prevalent cultural variants (e.g., beliefs or values), has a strong adverse effect on the evolution of culturally inherited helping traits. This finding is at variance with the commonly held view that conformist transmission is a crucial factor favoring the evolution of altruistic helping in humans. By contrast, we find that under one-to-many transmission, where individuals within groups tend to copy a "leader" (or teacher), altruistic helping can evolve in groups of any size, although the cultural transmission rule itself hitchhikes rather weakly with a selected helping trait. Our results suggest that culturally determined helping behaviors are more likely to be driven by "leaders" than by popularity, but the emergence and stability of the cultural transmission rules themselves should be driven by some extrinsic factors.

HIV rebounds from latently infected cells, rather than from continuing low-level replication.

Rapid rebound of plasma viremia in patients after interruption of long-term combination antiretroviral therapy (cART) suggests persistence of low-level replicating cells or rapid reactivation of latently infected cells. To further characterize rebounding virus, we performed extensive longitudinal clonal evolutionary studies of HIV env C2-V3-C3 regions and exploited the temporal relationships of rebounding plasma viruses with regard to pretreatment sequences in 20 chronically HIV-1-infected patients having undergone multiple 2-week structured treatment interruptions (STI). Rebounding virus during the short STI was homogeneous, suggesting mono- or oligoclonal origin during reactivation. No evidence for a temporal structure of rebounding virus in regard to pretreatment sequences was found. Furthermore, expansion of distinct lineages at different STI cycles emerged. Together, these findings imply stochastic reactivation of different clones from long-lived latently infected cells rather than expansion of viral populations replicating at low levels. After treatment was stopped, diversity increased steadily, but pretreatment diversity was, on average, achieved only >2.5 years after the start of STI when marked divergence from preexisting quasispecies also emerged. In summary, our results argue against persistence of ongoing low-level replication in patients on suppressive cART. Furthermore, a prolonged delay in restoration of pretreatment viral diversity after treatment interruption demonstrates a surprisingly sustained evolutionary bottleneck induced by punctuated antiretroviral therapy.

Identity and combinations of arbuscular mycorrhizal fungal isolates influence plant resistance and insect preference

1. Accumulating evidence indicates that plant resistance against above-ground herbivores can be affected by the presence of arbuscular mycorrhizal fungi (AMF) in association with the host plant. Little is known, however, about how AMF composition can influence herbivore choice to feed on a particular plant. 2. Unravelling the preference-performance hypothesis in a multitrophic context is needed to expand our knowledge of complex multitrophic interactions in natural systems. If given mycorrhizal fungal genotypes increase attractiveness for a herbivore (reduced plant resistance), then the benefits of increased unpalatability provided by the mycorrhizal fungi (increased plant resistance) might be outweighed by the increased herbivore recruitment. 3. This was addressed by designing three experiments to test the effects of different AMF genotypes, inoculated either alone or in combination, to measure intraspecific AMF effects on plant resistance and insect herbivore preference. Using strawberry (Fragaria vesca L.) plants that were colonised by eight different combinations of Rhizophagus irregularis isolates, we measured effects on plant growth, insect growth and survival, as well as feeding preferences of a generalist herbivore caterpillar (Spodoptera littoralis Boisduval). 4. Overall, it was found that: (i) AMF influenced plant resistance in an AMF genotype-specific manner; (ii) some AMF inoculations decreased insect performance; (iii) insects preferentially chose to feed more on leaves originating from non-mycorrhizal plants; but also that (iv) in a whole plant bioassay, insects preferentially chose the biggest plant, regardless of their mycorrhizal status. 5. Therefore, AMF-mediated trade-offs between growth and resistance against herbivores have been shown. Such trade-offs, particularly driven by plant attractiveness to herbivores, buffer the positive effects of the mycorrhizal symbiosis on enhanced plant growth.

High quantitative and no molecular differentiation of a freshwater snail (Galba truncatula) between temporary and permanent water habitats.

We investigate the variation in quantitative and molecular traits in the freshwater snail Galba truncatula, from permanent and temporary water habitats. Using a common garden experiment, we measured 20 quantitative traits and molecular variation using seven microsatellites in 17 populations belonging to these two habitats. We estimated trait means in each habitat. We also estimated the distributions of overall genetic quantitative variation (QST), and of molecular variation (FST), within and between habitats. Overall, we observed a lack of association between molecular and quantitative variance. Among habitats, we found QST>FST, an indication of selection for different optima. Individuals from temporary water habitat matured older, at a larger size and were less fecund than individuals from permanent water habitat. We discuss these findings in the light of several theories for life-history traits evolution.