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.

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.

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 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.

Stepwise evaluation of unexplained syncope in a large ambulatory population.

BACKGROUND: Up to 60% of syncopal episodes remain unexplained. We report the results of a standardized, stepwise evaluation of patients referred to an ambulatory clinic for unexplained syncope. METHODS AND RESULTS: We studied 939 consecutive patients referred for unexplained syncope, who underwent a standardized evaluation, including history, physical examination, electrocardiogram, head-up tilt testing (HUTT), carotid sinus massage (CSM) and hyperventilation testing (HYV). Echocardiogram and stress test were performed when underlying heart disease was initially suspected. Electrophysiological study (EPS) and implantable loop recorder (ILR) were used only in patients with underlying structural heart disease or major unexplained syncope. We identified a cause of syncope in 66% of patients, including 27% vasovagal, 14% psychogenic, 6% arrhythmias, and 6% hypotension. Noninvasive testing identified 92% and invasive testing an additional 8% of the causes. HUTT yielded 38%, CSM 28%, HYV 49%, EPS 22%, and ILR 56% of diagnoses. On average, patients with arrhythmic causes were older, had a lower functional capacity, longer P-wave duration, and presented with fewer prodromes than patients with vasovagal or psychogenic syncope. CONCLUSIONS: A standardized stepwise evaluation emphasizing noninvasive tests yielded 2/3 of causes in patients referred to an ambulatory clinic for unexplained syncope. Neurally mediated and psychogenic mechanisms were behind >50% of episodes, while cardiac arrhythmias were uncommon. Sudden syncope, particularly in older patients with functional limitations or a prolonged P-wave, suggests an arrhythmic cause.

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é large contexte phylogénétique et en utilisant des régions géniques adéquates.Ce travail de thèse a contribué à la compréhension de l'évolution à long-terme du CMH aviaire classe II Β (CMH II Β). Nous avons isolé et caractérisé les gènes CMH I IB de la chouette effraie (Tyto alba; Strigiformes, Tytonidae), une espèce appartenant à une lignée aviaire dont le CMH n'avait pas été étudié jusqu'ici. Avec uniquement deux gènes CMHIIB fonctionnels, l'organisation CMH de la chouette effraie semble être comparable à celle "minimale" du poulet (Gallus gallus), ce qui suggère qu'une organisation simple du CMH pourrait être ancestrale chez les oiseaux. En se basant sur l'information obtenue pour cette espèce, nous avons ensuite étudié l'évolution des gènes CMHIIB chez 13 espèces de chouettes "typiques" (Strigiformes, Strigidae). Les analyses phylogénétiques ont révélé que chez les chouettes, selon leur fonction, l'exon 2 codant pour la 'peptide-binding region' (PBR) et l'exon 3 ne codant pas pour la région PBR évoluent de manière distincte. L'exon 2 montre une histoire évolutive de sélection positive et recombinaison, tandis que l'exon 3 retrace l'histoire de duplication et révèle deux paralogues évoluant de façon divergente chez les chouettes, ainsi que chez un limicole, la bécassine double (Gallinago media). Les résultats découlant de l'analyse de l'exon 3 ont montré pour la première fois chez les oiseaux une orthologie de gènes chez des espèces ayant divergé il y a des dizaines de millions d'années. Dans une deuxième étude, nous avons voulu expliquer le pattern surprenant révélé par les arbres phylogénétiques lorsque les séquences CMHIIB des chouettes sont analysées avec les séquences d'autres oiseaux: un paralogue de chouette (appelé DABI) forme un groupe avec les séquences de passereaux et faucons, tandis que l'autre (DAB2) se regroupe avec les gallo-ansériformes, les manchots et les rapaces diurnes. Ceci pourrait être expliqué soit par une duplication qui a précédé l'évolution de ces ordres aviaires soit par évolution convergente dans un certain nombre d'ordres. Grâce à des analyses phylogénétiques approfondies, nous avons pu montrer qu'un événement de duplication a effectivement précédé la radiation aviaire principale d'il y a environ 100 millions d'années, et qu'après cette duplication les paralogues ont évolué sous sélection positive. De plus, les résidus d'acides aminés évoluant de façon divergente dans la chaîne β codée par le CMHIIB interagissent potentiellement avec la chaîne α codée par le CMHII, et que la coévolution moléculaire de ces résidus pourrait être à la base de l'évolution divergente des paralogues CMHIIB.Dans cette thèse sont discutés entre autres l'importance d'inclure la structure des protéines dans l'étude phylogénétique des familles multigéniques, ainsi que les rôles respectifs des pressions sélectives écologique et moléculaire. La conclusion comporte une perspective sur la génomique des populations du CMH aviaire, qui pourrait servir de base pour de futures recherches sur les rôles relatifs joués par Îes processus neutres comprenant les effets de la taille efficace des populations et l'adaptation dans l'évolution de la diversité et l'organisation du MHC aviaire.

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.

Molecular changes underlying mammalian phenotypic innovation

Mammals are characterized by specific phenotypic traits that include lactation, hair, and relatively large brains with unique structures. Individual mammalian lineages have, in turn, evolved characteristic traits that distinguish them from others. These include obvious anatom¬ical differences but also differences related to reproduction, life span, cognitive abilities, be¬havior. and disease susceptibility. However, the molecular basis of the diverse mammalian phenotypes and the selective pressures that shaped their evolution remain largely unknown. In the first part of my thesis, I analyzed the genetic factors associated with the origin of a unique mammalian phenotype lactation and I studied the selective pressures that forged the transition from oviparity to viviparity. Using a comparative genomics approach and evolutionary simulations, I showed that the emergence of lactation, as well as the appear¬ance of the casein gene family, significantly reduced selective pressure on the major egg-yolk proteins (the vitellogenin family). This led to a progressive loss of vitellogenins, which - in oviparous species - act as storage proteins for lipids, amino acids, phosphorous and calcium in the isolated egg. The passage to internal fertilization and placentation in therian mam¬mals rendered vitellogenins completely dispensable, which ended in the loss of the whole gene family in this lineage. As illustrated by the vitellogenin study, changes in gene content are one possible underlying factor for the evolution of mammalian-specific phenotypes. However, more subtle genomic changes, such as mutations in protein-coding sequences, can also greatly affect the phenotypes. In particular, it was proposed that changes at the level of gene reg¬ulation could underlie many (or even most) phenotypic differences between species. In the second part of my thesis, I participated in a major comparative study of mammalian tissue transcriptomes, with the goal of understanding how evolutionary forces affected expression patterns in the past 200 million years of mammalian evolution. I showed that, while com¬parisons of gene expressions are in agreement with the known species phylogeny, the rate of expression evolution varies greatly among lineages. Species with low effective population size, such as monotremes and hominoids, showed significantly accelerated rates of gene expression evolution. The most likely explanation for the high rate of gene expression evolution in these lineages is the accumulation of mildly deleterious mutations in regulatory regions, due to the low efficiency of purifying selection. Thus, our observations are in agreement with the nearly neutral theory of molecular evolution. I also describe substantial differences in evolutionary rates between tissues, with brain being the most constrained (especially in primates) and testis significantly accelerated. The rate of gene expression evolution also varies significantly between chromosomes. In particular, I observed an acceleration of gene expression changes on the X chromosome, probably as a result of adaptive processes associated with the origin of therian sex chromosomes. Lastly, I identified several individual genes as well as co-regulated expression modules that have undergone lineage specific expression changes and likely under¬lie various phenotypic innovations in mammals. The methods developed during my thesis, as well as the comprehensive gene content analyses and transcriptomics datasets made available by our group, will likely prove to be useful for further exploratory analyses of the diverse mammalian phenotypes.

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).

Population demography of an endangered lizard, the Blue Mountains Water Skink.

BACKGROUND: Information on the age structure within populations of an endangered species can facilitate effective management. The Blue Mountains Water Skink (Eulamprus leuraensis) is a viviparous scincid lizard that is restricted to < 40 isolated montane swamps in south-eastern Australia. We used skeletochronology of phalanges (corroborated by mark-recapture data) to estimate ages of 222 individuals from 13 populations. RESULTS: These lizards grow rapidly, from neonatal size (30 mm snout-vent length) to adult size (about 70 mm SVL) within two to three years. Fecundity is low (mean 2.9 offspring per litter) and is affected by maternal body length and age. Offspring quality may decline with maternal age, based upon captive-born neonates (older females gave birth to slower offspring). In contrast to its broadly sympatric (and abundant) congener E. tympanum, E. leuraensis is short-lived (maximum 6 years, vs 15 years for E. tympanum). Litter size and offspring size are similar in the two species, but female E. leuraensis reproduce annually whereas many E. tympanum produce litters biennially. Thus, a low survival rate (rather than delayed maturation or low annual fecundity) is the key reason why E. leuraensis is endangered. Our 13 populations exhibited similar growth rates and population age structures despite substantial variation in elevation, geographic location and swamp size. However, larger populations (based on a genetic estimate of effective population size) contained older lizards, and thus a wider variance in ages. CONCLUSION: Our study suggests that low adult survival rates, as well as specialisation on a rare and fragmented habitat type (montane swamps) contribute to the endangered status of the Blue Mountains Water Skink.

Molecular substitution rate increases in myrmecophilous lycaenid butterflies (Lepidoptera)

Is species diversification driven by neutral- or niche-based processes? Butterflies of the Lycaenidae family have developed mutualistic interactions with ants. This biotic requirement increased the spatial fragmentation of populations of lower effective population size (Ne) compared with autonomous species. The nearly neutral theory predicts that species with smaller Ne should fix more mutations because of the increased strength of drift. Taking into account the phylogenetic relatedness among species, this study shows that species with a stronger dependence on ants displayed more intra-specific Single Nucleotide Polymorphisms compared with species with low or no myrmecophily. This phenomenon can cause more pronounced genetic differentiation between populations and could ultimately promote speciation in a similar manner as on physical islands. The large species diversity observed in this family could be the consequence of this neutral process enhancing the diversification of lineages.

Stepwise evaluation of unexplained syncope in a large ambulatory population

La syncope est un symptôme clinique fréquent mais son origine demeure indéterminée jusque dans 60% des cas de patients admis dans un centre d'urgences. Le développement de consultations spécialisées de la syncope a considérablement modifié l'évaluation des patients avec une syncope inexpliquée en les orientant vers des stratégies d'investigations non-invasives, tels que le tilt-test, le massage du sinus carotidien et le test ^hyperventilation. Cependant, il existe peu de données dans 10 la littérature concernant dans la performance diagnostique réelle de ces tests fonctionnels.Notre travail de recherche porte sur l'analyse des données des 939 premiers patients adressés à la consultation ambulatoire de la syncope du CHUV pour l'investigation d'une syncope d'origine indéterminée. L'objectif de notre travail de thèse est 1) d'évaluer la performance diagnostique de l'algorithme de prise en charge standardisé et de ses différents tests pratiqués dans le cadre de notre 15 consultation et 2) de déterminer les caractéristiques cliniques communes des patients avec un diagnostic final de syncope d'origine rythmique ou vaso-vagale.Notre travail de thèse démontre qu'un algorithme de prise en charge standardisé basé sur des tests non-invasifs permet de déterminer 2/3 des causes de syncope initialement d'origine indéterminée. Par ailleurs, notre travail montre que des étiologies bénignes, telles que la syncope d'origine vaso- 20 vagale ou psychogène, représentent la moitié des causes syncopales alors que les arythmies cardiaques demeurent peu fréquentes. Finalement, notre travail démontre que l'absence de symptomatologie prodromique, en particulier chez les patients âgés avec une limitation fonctionnelle ou un allongement de la durée de l'onde Ρ à l'électrocardiogramme, suggère une syncope d'origine rythmique. Ce travail de thèse contribuera à optimaliser notre algorithme de prise 25 en charge standardisée de la syncope d'origine indéterminée et ouvre de nouvelles perspectives de recherche dans le développement de modèles basés sur des facteurs cliniques permettant de prédire les principales causes syncopales.

Genetic bottlenecks driven by population disconnection.

Connectivity among populations plays a crucial role in maintaining genetic variation at a local scale, especially in small populations affected strongly by genetic drift. The negative consequences of population disconnection on allelic richness and gene diversity (heterozygosity) are well recognized and empirically established. It is not well recognized, however, that a sudden drop in local effective population size induced by such disconnection produces a temporary disequilibrium in allelic frequency distributions that is akin to the genetic signature of a demographic bottleneck. To document this effect, we used individual-based simulations and empirical data on allelic richness and gene diversity in six pairs of isolated versus well-connected (core) populations of European tree frogs. In our simulations, population disconnection depressed allelic richness more than heterozygosity and thus resulted in a temporary excess in gene diversity relative to mutation drift equilibrium (i.e., signature of a genetic bottleneck). We observed a similar excess in gene diversity in isolated populations of tree frogs. Our results show that population disconnection can create a genetic bottleneck in the absence of demographic collapse.