Population genomics of eusocial insects: the costs of a vertebrate-like effective population size

The evolution of reproductive division of labour and social life in social insects has lead to the emergence of several life-history traits and adaptations typical of larger organisms: social insect colonies can reach masses of several kilograms, they start reproducing only when they are several years old, and can live for decades. These features and the monopolization of reproduction by only one or few individuals in a colony should affect molecular evolution by reducing the effective population size. We tested this prediction by analysing genome-wide patterns of coding sequence polymorphism and divergence in eusocial vs. noneusocial insects based on newly generated RNA-seq data. We report very low amounts of genetic polymorphism and an elevated ratio of nonsynonymous to synonymous changes - a marker of the effective population size - in four distinct species of eusocial insects, which were more similar to vertebrates than to solitary insects regarding molecular evolutionary processes. Moreover, the ratio of nonsynonymous to synonymous substitutions was positively correlated with the level of social complexity across ant species. These results are fully consistent with the hypothesis of a reduced effective population size and an increased genetic load in eusocial insects, indicating that the evolution of social life has important consequences at both the genomic and population levels.

Substitution rates at neutral genes depend on population size under fluctuating demography and overlapping generations.

It is widely accepted that the rate of evolution (substitution rate) at neutral genes is unaffected by population size fluctuations. This result has implications for the analysis of genetic data in population genetics and phylogenetics, and provides, in particular, a justification for the concept of the molecular clock. Here, we show that the substitution rate at neutral genes does depend on population size fluctuations in the presence of overlapping generations. As both population size fluctuations and overlapping generations are expected to be the norm rather than the exception in natural populations, this observation may be relevant for understanding variation in substitution rates within and between lineages.

Coevolution of adaptive technology, maladaptive culture and population size in a producer-scrounger game.

Technology (i.e. tools, methods of cultivation and domestication, systems of construction and appropriation, machines) has increased the vital rates of humans, and is one of the defining features of the transition from Malthusian ecological stagnation to a potentially perpetual rising population growth. Maladaptations, on the other hand, encompass behaviours, customs and practices that decrease the vital rates of individuals. Technology and maladaptations are part of the total stock of culture carried by the individuals in a population. Here, we develop a quantitative model for the coevolution of cumulative adaptive technology and maladaptive culture in a 'producer-scrounger' game, which can also usefully be interpreted as an 'individual-social' learner interaction. Producers (individual learners) are assumed to invent new adaptations and maladaptations by trial-and-error learning, insight or deduction, and they pay the cost of innovation. Scroungers (social learners) are assumed to copy or imitate (cultural transmission) both the adaptations and maladaptations generated by producers. We show that the coevolutionary dynamics of producers and scroungers in the presence of cultural transmission can have a variety of effects on population carrying capacity. From stable polymorphism, where scroungers bring an advantage to the population (increase in carrying capacity), to periodic cycling, where scroungers decrease carrying capacity, we find that selection-driven cultural innovation and transmission may send a population on the path of indefinite growth or to extinction.

Individual reproductive success and effective population size in the greater white-toothed shrew Crocidura russula.

In order to investigate the determinants of effective population size in the socially monogamous Crocidura russula, the reproductive output of 44 individuals was estimated through genetic assignment methods. The individual variance in breeding success turned out to be surprisingly high, mostly because the males were markedly less monogamous than expected from previous behavioural data. Males paired simultaneously with up to four females and polygynous males had significantly more offspring than monogamous ones. The variance in female reproductive success also exceeded that of a Poisson distribution (though to a lesser extent), partly because females paired with multiply mated males weaned significantly more offspring. Polyandry also occurred occasionally, but only sequentially (i.e. without multiple paternity of litters). Estimates of the effective to census size ratio were ca. 0.60, which excluded the mating system as a potential explanation for the high genetic variance found in this shrew's populations. Our data suggest that gene flow from the neighbourhood (up to one-third of the total recruitment) is the most likely cause of the high levels of genetic diversity observed in this shrew's subpopulations.

Size-dependent predation by Dugesia lugubris (Turbellaria) on Physa acuta (Gastropoda): experiments and model

1. We investigated experimentally predation by the flatworm Dugesia lugubris on the snail Physa acuta in relation to predator body length and to prey morphology [shell length (SL) and aperture width (AW)]. 2. SL and AW correlate strongly in the field, but display significant and independent variance among populations. In the laboratory, predation by Dugesia resulted in large and significant selection differentials on both SL and AW. Analysis of partial effects suggests that selection on AW was indirect, and mediated through its strong correlation with SL. 3. The probability P(ij) for a snail of size category i (SL) to be preyed upon by a flatworm of size category j was fitted with a Poisson-probability distribution, the mean of which increased linearly with predator size (i). Despite the low number of parameters, the fit was excellent (r2 = 0.96). We offer brief biological interpretations of this relationship with reference to optimal foraging theory. 4. The largest size class of Dugesia (>2 cm) did not prey on snails larger than 7 mm shell length. This size threshold might offer Physa a refuge against flatworm predation and thereby allow coexistence in the field. 5. Our results are further discussed with respect to previous field and laboratory observations on P acuta life-history patterns, in particular its phenotypic variance in adult body size.

Testing demographic models of effective population size.

Due to practical difficulties in obtaining direct genetic estimates of effective sizes, conservation biologists have to rely on so-called 'demographic models' which combine life-history and mating-system parameters with F-statistics in order to produce indirect estimates of effective sizes. However, for the same practical reasons that prevent direct genetic estimates, the accuracy of demographic models is difficult to evaluate. Here we use individual-based, genetically explicit computer simulations in order to investigate the accuracy of two such demographic models aimed at investigating the hierarchical structure of populations. We show that, by and large, these models provide good estimates under a wide range of mating systems and dispersal patterns. However, one of the models should be avoided whenever the focal species' breeding system approaches monogamy with no sex bias in dispersal or when a substructure within social groups is suspected because effective sizes may then be strongly overestimated. The timing during the life cycle at which F-statistics are evaluated is also of crucial importance and attention should be paid to it when designing field sampling since different demographic models assume different timings. Our study shows that individual-based, genetically explicit models provide a promising way of evaluating the accuracy of demographic models of effective size and delineate their field of applicability.

Optimal growth strategies when mortality and production rates are size-dependent

Pontryagin's maximum principle from optimal control theory is used to find the optimal allocation of energy between growth and reproduction when lifespan may be finite and the trade-off between growth and reproduction is linear. Analyses of the optimal allocation problem to date have generally yielded bang-bang solutions, i.e. determinate growth: life-histories in which growth is followed by reproduction, with no intermediate phase of simultaneous reproduction and growth. Here we show that an intermediate strategy (indeterminate growth) can be selected for if the rates of production and mortality either both increase or both decrease with increasing body size, this arises as a singular solution to the problem. Our conclusion is that indeterminate growth is optimal in more cases than was previously realized. The relevance of our results to natural situations is discussed.

Opportunity for sexual selection and effective population size in the lek-breeding European treefrog (Hyla arborea).

Sexual selection in lek-breeding species might drastically lower male effective population size, with potentially important consequences for evolutionary and conservation biology. Using field-monitoring and parental-assignment methods, we analyzed sex-specific variances in breeding success in a population of European treefrogs, to (1) help understanding the dynamics of genetic variance at sex-specific loci, and (2) better quantify the risk posed by genetic drift in this species locally endangered by habitat fragmentation. The variance in male mating success turned out to be markedly lower than values obtained from other amphibian species with polygamous mating systems. The ratio of effective breeding size to census breeding size was only slightly lower in males (0.44) than in females (0.57), in line with the patterns of genetic diversity previously reported from H. arborea sex chromosomes. Combining our results with data on age at maturity and adult survival, we show that the negative effect of the mating system is furthermore compensated by the effect of delayed maturity, so that the estimated instantaneous effective size broadly corresponded to census breeding size. We conclude that the lek-breeding system of treefrogs impacts only weakly the patterns of genetic diversity on sex-linked genes and the ability of natural populations to resist genetic drift.

Population demography and the evolution of helping behaviors.

Limited dispersal may favor the evolution of helping behaviors between relatives as it increases their relatedness, and it may inhibit such evolution as it increases local competition between these relatives. Here, we explore one way out of this dilemma: if the helping behavior allows groups to expand in size, then the kin-competition pressure opposing its evolution can be greatly reduced. We explore the effects of two kinds of stochasticity allowing for such deme expansion. First, we study the evolution of helping under environmental stochasticity that may induce complete patch extinction. Helping evolves if it results in a decrease in the probability of extinction or if it enhances the rate of patch recolonization through propagules formed by fission of nonextinct groups. This mode of dispersal is indeed commonly found in social species. Second, we consider the evolution of helping in the presence of demographic stochasticity. When fecundity is below its value maximizing deme size (undersaturation), helping evolves, but under stringent conditions unless positive density dependence (Allee effect) interferes with demographic stochasticity. When fecundity is above its value maximizing deme size (oversaturation), helping may also evolve, but only if it reduces negative density-dependent competition.

The population genetics of clonal and partially clonal diploids.

The consequences of variable rates of clonal reproduction on the population genetics of neutral markers are explored in diploid organisms within a subdivided population (island model). We use both analytical and stochastic simulation approaches. High rates of clonal reproduction will positively affect heterozygosity. As a consequence, nearly twice as many alleles per locus can be maintained and population differentiation estimated as F(ST) value is strongly decreased in purely clonal populations as compared to purely sexual ones. With increasing clonal reproduction, effective population size first slowly increases and then points toward extreme values when the reproductive system tends toward strict clonality. This reflects the fact that polymorphism is protected within individuals due to fixed heterozygosity. Contrarily, genotypic diversity smoothly decreases with increasing rates of clonal reproduction. Asexual populations thus maintain higher genetic diversity at each single locus but a lower number of different genotypes. Mixed clonal/sexual reproduction is nearly indistinguishable from strict sexual reproduction as long as the proportion of clonal reproduction is not strongly predominant for all quantities investigated, except for genotypic diversities (both at individual loci and over multiple loci).

Within- and among-population impact of genetic erosion on adult fitness-related traits in the European tree frog Hyla arborea.

Assessing in wild populations how fitness is impacted by inbreeding and genetic drift is a major goal for conservation biology. An approach to measure the detrimental effects of inbreeding on fitness is to estimate correlations between molecular variation and phenotypic performances within and among populations. Our study investigated the effect of individual multilocus heterozygosity on body size, body condition and reproductive investment of males (that is, chorus attendance) and females (that is, clutch mass and egg size) in both small fragmented and large non-fragmented populations of European tree frog (Hyla arborea). Because adult size and/or condition and reproductive investment are usually related, genetic erosion may have detrimental effects directly on reproductive investment, and also on individual body size and condition that in turn may affect reproductive investment. We confirmed that the reproductive investment was highly size-dependent for both sexes. Larger females invested more in offspring production, and larger males attended the chorus in the pond more often. Our results did not provide evidence for a decline in body size, condition and reproductive effort with decreased multilocus heterozygosity both within and among populations. We showed that the lack of heterozygosity-fitness correlations within populations probably resulted from low inbreeding levels (inferior to ca. 20% full-sib mating rate), even in the small fragmented populations. The detrimental effects of fixation load were either low in adults or hidden by environmental variation among populations. These findings will be useful to design specific management actions to improve population persistence.

Population Normalization with Ammonium in Wastewater-Based Epidemiology: Application to Illicit Drug Monitoring

Fluctuations in ammonium (NH4+), measured as NH4-N loads using an ion-selective electrode installed at the inlet of a sewage treatment plant, showed a distinctive pattern which was associated to weekly (i.e., commuters) and seasonal (i.e., holidays) fluctuations of the population. Moreover, population size estimates based on NH4-N loads were lower compared to census data. Diurnal profiles of benzoylecgonine (BE) and 11-nor-9-carboxy-Δ9-tetrahydrocannabinol (THC-COOH) were shown to be strongly correlated to NH4-N. Characteristic patterns, which reflect the prolonged nocturnal activity of people during the weekend, could be observed for BE, cocaine, and a major metabolite of MDMA (i.e., 4-hydroxy-3-methoxymethamphetamine). Additional 24 h composite samples were collected between February and September 2013. Per-capita loads (i.e., grams per day per 1000 inhabitants) were computed using census data and NH4-N measurements. Normalization with NH4-N did not modify the overall pattern, suggesting that the magnitude of fluctuations in the size of the population is negligible compared to those of illicit drug loads. Results show that fluctuations in the size of the population over longer periods of time or during major events can be monitored using NH4-N loads: either using raw NH4-N loads or population size estimates based on NH4-N loads, if information about site-specific NH4-N population equivalents is available.

The tumor suppressor p16Ink4a regulates T lymphocyte survival.

In contrast to other cell cycle inhibitors, the tumor suppressor p16Ink4a is not detectable or expressed at very low levels in embryonic and adult mouse tissues, and therefore it has often been considered as a specialized checkpoint protein that does not participate in the control of normal cell cycle progression. However, Ink4a-/- mice possess increased thymus size and cellularity, thus suggesting the involvement of p16(Ink4a) in the control of thymocyte proliferation. In this study, we found increased numbers of CD8 and CD4 T lymphocytes in thymus and spleen from Ink4a-/- mice. Unexpectedly, this was not related to an increase in T-cell division rates, which were similar in lymphoid organs of Ink4a-/- and wild-type mice. In contrast, T-cell apoptosis rates were significantly decreased in thymus and spleen from Ink4a-/- mice. Moreover, whereas p16Ink4a-deficient and wild-type T cells were equally sensitive to Fas or TCR-mediated apoptosis, the former were clearly more resistant to apoptosis induced by oxidative stress or gamma irradiation. Our results indicate that p16Ink4a function is associated with T-cell apoptosis, and subsequently contributes to the control of T-cell population size in lymphoid organs.

New insights into the long-term evolutionary history of the avian MHC

SummaryGene duplication and neofunctidnalization are important processes in the evolution of phenotypic complexity. They account for important evolutionary novelties that confer ecological adaptation, such as the major histocompatibility complex (MHC), a multigene family with a central role in vertebrates' adaptive immune system. Multigene families, which evolved in large part through duplication, represent promising systems to study the still strongly depbated relative roles of neutral and adaptive processes in the evolution of phenotypic complexity. Detailed knowledge on ecological function and a well-characterized evolutionary history place the mammals' MHC amongst ideal study systems. However mammalian MHCs usually encompass several million base pairs and hold a large number of functional and non-functional duplicate genes, which makes their study complex. Avian MHCs on the other hand are usually way more compact, but the reconstruction of. their evolutionary history has proven notoriously difficult. However, no focused attempt has been undertaken so far to study the avian MHC evolutionary history in a broad phylogenetic context and using adequate gene regions.In the present PhD, we were able to make important contributions to the understanding of the long-term evolution of the avian MHC class II Β (MHCI1B). First, we isolated and characterized MHCIIB genes in barn owl (Tyto alba?, Strigiformes, Tytonidae), a species from an avian lineage in which MHC has not been studied so far. Our results revealed that with only two functional MHCIIB genes the MHC organization of barn owl may be similar to the 'minimal essential' MHC of chicken (Gallus gallus), indicating that simple MHC organization may be ancestral to birds. Taking advantage of the sequence information from barn owl, we studied the evolution of MHCIIB genes in 13 additional species of 'typical' owls (Strigiformes, Strigidae). Phylogenetic analyses revealed that according to their function, in owls the peptide-binding region (PBR) encoding exon 2 and the non-PBR encoding exon 3 evolve by different patterns. Exon 2 exhibited an evolutionary history of positive selection and recombination, while exon 3 traced duplication history and revealed two paralogs evolving divergently from each other in owls, and in a shorebird, the great snipe {Gallinago media). The results from exon 3 were the first ever from birds to demonstrate gene orthology in species that diverged tens of millions of years ago, and strongly questioned whether the taxa studied before provided an adequate picture of avian MHC evolution. In a follow-up study, we aimed at explaining a striking pattern revealed by phylogenetic trees analyzing the owl sequences along with MHCIIB sequences from other birds: One owl paralog (termed DAB1) grouped with sequences of passerines and falcons, while the other (DAB2) grouped with wildfowl, penguins and birds of prey. This could be explained by either a duplication event preceding the evolution of these bird orders, or by convergent evolution of similar sequences in a number of orders. With extensive phylogenetic analyses we were able to show, that indeed a duplication event preceeded the major avian radiation -100 my ago, and that following this duplication, the paralogs evolved under positive selection. Furthermore, we showed that the divergently evolving amino acid residues in the MHCIIB-encoded β-chain potentially interact with the MHCI I α-chain, and that molecular coevolution of the interacting residues may have been involved in the divergent evolution of the MHCIIB paralogs.The findings of this PhD are of particular interest to the understanding of the evolutionary history of the avian MHC and, by providing essential information on long-term gene history in the avian MHC, open promising perspectives for advances in the understanding of the evolution of multigene families in general, and for avian MHC organization in particular. Amongst others I discuss the importance of including protein structure in the phylogenetic study of multigene families, and the roles of ecological versus molecular selection pressures. I conclude by providing a population genomic perspective on avian MHC, which may serve as a basis for future research to investigate the relative roles of neutral processes involving effective population size effects and of adaptation in the evolution of avian MHC diversity and organization.RésuméLa duplication de gènes et leur néo-fonctionnalisation sont des processus importants dans l'évolution de la complexité phénotypique. Ils sont impliqués dans l'apparition d'importantes nouveautés évolutives favorisant l'adaptation écologique, comme c'est le cas pour le complexe majeur d'histocompatibilité

Introduction of Trojan sex chromosomes to boost population growth.

Conservation programs that deal with small or declining populations often aim at a rapid increase of population size to above-critical levels in order to avoid the negative effects of demographic stochasticity and genetic problems like inbreeding depression, fixation of deleterious alleles, or a general loss of genetic variability and hence of evolutionary potential. In some situations, population growth is determined by the number of females available for reproduction, and manipulation of family sex ratios towards more daughters has beneficial effects. If sex determination is predominantly genetic but environmentally reversible, as is the case in many amphibia, reptiles, and fish, Trojan sex chromosomes could be introduced into populations in order to change sex ratios towards more females. We analyse the possible consequences for the introduction of XX-males (XX individuals that have been changed to phenotypic males in a XY/XX sex determination system) and ZW males, WW males, or WW females (in a ZZ/ZW sex determination system). We find that the introduction of WW individuals can be most effective for an increase of population growth, especially if the induced sex change has little or no effect on viability.