Reliable selfing rate estimates from imperfect population genetic data.

Genotypic frequencies at codominant marker loci in population samples convey information on mating systems. A classical way to extract this information is to measure heterozygote deficiencies (FIS) and obtain the selfing rate s from FIS = s/(2 - s), assuming inbreeding equilibrium. A major drawback is that heterozygote deficiencies are often present without selfing, owing largely to technical artefacts such as null alleles or partial dominance. We show here that, in the absence of gametic disequilibrium, the multilocus structure can be used to derive estimates of s independent of FIS and free of technical biases. Their statistical power and precision are comparable to those of FIS, although they are sensitive to certain types of gametic disequilibria, a bias shared with progeny-array methods but not FIS. We analyse four real data sets spanning a range of mating systems. In two examples, we obtain s = 0 despite positive FIS, strongly suggesting that the latter are artefactual. In the remaining examples, all estimates are consistent. All the computations have been implemented in a open-access and user-friendly software called rmes (robust multilocus estimate of selfing) available at http://ftp.cefe.cnrs.fr, and can be used on any multilocus data. Being able to extract the reliable information from imperfect data, our method opens the way to make use of the ever-growing number of published population genetic studies, in addition to the more demanding progeny-array approaches, to investigate selfing rates.

Genetic improvement and population structure of the Nelore breed in Northern Brazil

The objective of this work was to evaluate the population structure and the genetic and phenotypic progress of Nelore cattle in Northern Brazil. Pedigree information concerning animals born between 1942 and 2006 were analyzed. Population structure was performed using the Endog program. Out of the 140,628 animals studied, 67.7, 14.52 and 3.18% had complete pedigree record of the first, second and third parental generation, respectively. Inbreeding and average relatedness coefficients were low: 0.2 and 0.13%, respectively. However, these parameters may have been underestimated, since information on pedigree was incomplete. The effective number of founders was 370 and the genetic contribution of 10, 50 and 448 most influent ancestors explained 13.2, 28 and 50% of the genetic variability in the population, respectively. The genetic variability for growth traits and population structure demonstrates high probability of increasing productivity through selective breeding. Moreover, management strategies to reduce the currently observed age at first calving and generation intervals are important for Nelore cattle genetic improvement.

Characterization of thirty two microsatellite loci for three Atlanto-Mediterranean echinoderm species

Thirty two microsatellites were optimized from 454 pyrosequencing libraries for three Atlanto-Mediterranean echinoderms: Coscinasterias tenuispina, Echinaster sepositus and Arbacia lixula. We observed different frequency of microsatellite types (di-, tri-, tetra- and pentanucleotide) throughout the genome of the species, but no significant differences were observed in allele richness among different microsatellite repeats. No loci showed linkage disequilibrium. Heterozygosity deficit and departure from Hardy Weinberg equilibrium were observed for some loci, in two species, probably due to high levels of inbreeding. Heterozygosity excess observed in C. tenuispina could be explained by selection against homozygotes and/or outcrossing.

Mate availability and male dispersal in the Argentine ant Linepithema humile (Mayr) (=Iridomyrmex humilis)

In the Argentine ant Linepithema humile (=Iridomyrmex humilis) only males disperse whereas female sexuals (unmated winged queens) stay in their mother nest where they mate. This study investigated (1) whether dispersing males are accepted into foreign colonies, (2) whether they can mate with resident female sexuals, and (3) whether the propensity of males to disperse is affected by the expectation of mating in their mother nest. Field experiments demonstrated that males were accepted into foreign colonies only when these colonies contained female sexuals or queen pupae. Before and after the time of (sic) female sexuals, workers attacked and killed most of the foreign males. Laboratory experiments snowed that males that successfully enter foreign colonies can mate with resident female sexuals. The propensity of males to disperse was significantly influenced by the presence of female sexuals in their nest. Males were more likely to fly out from colonies containing no female sexuals than from those with them. These results are consistent with males preferentially dispersing when there is little or no opportunity to mate in their mother nest. Thus there are two mating strategies available for males: staying in their mother nest when an opportunity to mate arises or dispersing and attempting to mate in a foreign nest when there are no female sexuals in their mother nest. This latter behaviour could mediate gene flow between colonies and account for the lack of significant inbreeding previously documented in this species.

Mating system parameters in a high density population of andirobas in the Amazon forest

The objective of this work was to estimate the mating system parameters of a andiroba (Carapa guianensis) population using microsatellite markers and the mixed and correlated mating models. Twelve open‑pollinated progeny arrays of 15 individuals were sampled in an area with C. guianensis estimated density of 25.7 trees per hectare. Overall, the species has a mixed reproductive system, with a predominance of outcrossing. The multilocus outcrossing rate (t m = 0.862) was significantly lower than the unity, indicating that self‑pollination occurred. The rate of biparental inbreeding was substantial (t m ‑ t s = 0.134) and significantly different from zero. The correlation of selfing within progenies was high (r s = 0.635), indicating variation in the individual outcrossing rate. Consistent with this result, the estimate of the individual outcrossing rate ranged from 0.598 to 0.978. The multilocus correlation of paternity was low (r p(m) = 0.081), but significantly different from zero, suggesting that the progenies contain full‑sibs. The coancestry within progenies (Θ = 0.185) was higher and the variance effective size (Ne(v) = 2.7) was lower than expected for true half‑sib progenies (Θ = 0.125; Ne(v) = 4). These results suggest that, in order to maintain a minimum effective size of 150 individuals for breeding, genetic conservation, and environmental reforestation programs, seeds from at least 56 trees must be collected.

Comparative population genetic structure in a plant-pollinator/seed predator system.

Comparative analyses of spatial genetic structure of populations of plants and the insects they interact with provide an indication of how gene flow, natural selection and genetic drift may jointly influence the distribution of genetic variation and potential for local co-adaptation for interacting species. Here, we analysed the spatial scale of genetic structure within and among nine populations of an interacting species pair, the white campion Silene latifolia and the moth Hadena bicruris, along a latitudinal gradient across Northern/Central Europe. This dioecious, short-lived perennial plant inhabits patchy, often disturbed environments. The moth H. bicruris acts both as its pollinator and specialist seed predator that reproduces by laying eggs in S. latifolia flowers. We used nine microsatellite markers for S. latifolia and eight newly developed markers for H. bicruris. We found high levels of inbreeding in most populations of both plant and pollinator/seed predator. Among populations, significant genetic structure was observed for S. latifolia but not for its pollinator/seed predator, suggesting that despite migration among populations of H. bicruris, pollen is not, or only rarely, carried over between populations, thus maintaining genetic structure among plant populations. There was a weak positive correlation between genetic distances of S. latifolia and H. bicruris. These results indicate that while significant structure of S. latifolia populations creates the potential for differentiation at traits relevant for the interaction with the pollinator/seed predator, substantial gene flow in H. bicruris may counteract this process in at least some populations.

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.

Evolutionary aspects of population structure for molecular and quantitative traits in the freshwater snail Radix balthica.

Detecting the action of selection in natural populations can be achieved using the QST-FST comparison that relies on the estimation of FST with neutral markers, and QST using quantitative traits potentially under selection. QST higher than FST suggests the action of directional selection and thus potential local adaptation. In this article, we apply the QST-FST comparison to four populations of the hermaphroditic freshwater snail Radix balthica located in a floodplain habitat. In contrast to most studies published so far, we did not detect evidence of directional selection for local optima for any of the traits we measured: QST calculated using three different methods was never higher than FST. A strong inbreeding depression was also detected, indicating that outcrossing is probably predominant over selfing in the studied populations. Our results suggest that in this floodplain habitat, local adaptation of R. balthica populations may be hindered by genetic drift, and possibly altered by uneven gene flow linked to flood frequency.

Population genetics of Capercaillie (Tetrao urogallus) in the Jura and the Pyrenees: a non-invasive approach to avian conservation genetics

RÉSUMÉ Le Grand tétras est un galliforme de montagne apparenté au faisan et au tétras lyre. Il est distribué de manière continue à travers la toundra et les montagnes de moyenne altitude en Europe de l'ouest. Toutefois, les populations d'Europe de l'ouest ont subi un déclin constant au cours des derniers siècles. Les causes de ce déclin sont probablement liées à l'activité humaine, telle .que l'élevage ou le tourisme, qui ont engendré une modification et une fragmentation de l'habitat de l'espèce. Malheureusement, les populations soumises à de forts déclins démographiques peuvent subir des effets génétiques (augmentation de la consanguinité et perte de diversité génétique) pouvant diminuer leur potentiel de reproduction et conduire irrémédiablement à l'extinction. Cette thèse présente les analyses conduites dans le but d'estimer l'impact du déclin démographique des populations de Grand tétras sur l'étendue et la distribution de leur variabilité génétique dans le Jura et dans les Pyrénées. Du fait de la législation locale protégeant les tétraonidés en général, mais également en raison de la biologie très cryptique du Grand tétras, l'ensemble des analyses de cette étude a été réalisé à partir de matériel génétique extrait des fientes (ou échantillonnage génétique non invasif). Dans la première partie de l'étude, je détaille les protocoles d'extraction. d'ADN et d'amplification par PCR modifiés à partir des protocoles classiques utilisant des échantillons conventionnels, riches en ADN. L'utilisation d'ADN fécal impose des contraintes dues à la mauvaise qualité et à la faible quantité du matériel génétique à disposition dans les fientes. Ces contraintes ont pu être partiellement contournées en réalisant des répétitions multiples du génotypage afin d'obtenir un degré de fiabilité suffisante. J'ai également analysé les causes de la dégradation de l'ADN dans les excréments. Parmi les causes les plus communes, telles que l'activité bactérienne, l'hydrolyse spontanée et la dégradation enzymatique par les DNases libres, c'est ce dernier facteur qui apparaît comme étant la cause majeure et la plus rapide responsable de la dégradation de la qualité des échantillons. La rapidité de l'action enzymatique suggère que les plans d'échantillonnages de excréments sur le terrain pourraient être optimisés en les réalisant dans des conditions climatiques froides et sèches, favorisant ainsi l'inhibition des DNases. La seconde partie de la thèse est une étude par simulation visant à déterminer la capacité du logiciel Structure à identifier les structures génétiques complexes et hiérarchiques fréquemment rencontrées dans les populations naturelles, et ce en utilisant différents types de marqueurs génétiques. Les troisième et quatrième parties de cette thèse décrivent le statut génétique des populations résiduelles du Jura et des Pyrénées à partir de l'analyse de 11 loci microsatellites. Nous n'avons pas pu mettre en évidence dans les deux populations des effets liés à la consanguinité ou à la réduction de la diversité génétique. De plus, la différenciation génétique entre les patches d'habitats favorables reste modérée et corrélée à la distance géographique, ce qui suggère que la dispersion d'individus entre les patches a été importante au moins pendant ces dernières générations. La comparaison des paramètres de la diversité génétique avec ceux d'autres populations de Grand tétras, ou d'autres espèces proches, indique que la population du Jura a retenu une proportion importante de sa diversité originelle. Ces résultats suggèrent que le déclin récent des populations a jusqu'ici eu un impact modéré sur les facteurs génétiques et que ces populations semblent avoir conservé le potentiel génétique nécessaire à leur survie à long terme. Finalement, en cinquième partie, l'analyse de l'apparentement entre les mâles qui participent à la parade sur les places de chant (leks) indique que ces derniers sont distribués en agrégats de manière non aléatoire, préférentiellement entre individus apparentés. De plus, la corrélation entre les distances génétique et géographique entre les leks est en accord avec les motifs d'isolement par la distance mis en évidence à d'autres niveaux hiérarchiques (entre patches d'habitat et populations), ainsi qu'avec les études menées sur d'autres espèces ayant choisi ce même système de reproduction. En conclusion, cette première étude basée uniquement sur de l'ADN nucléaire aviaire extrait à partir de fèces a fourni des informations nouvelles qui n'auraient pas pu être obtenues par une méthode d'observation sur le terrain ou d'échantillonnage génétique classique. Aucun oiseau n'a été dérangé ou capturé, et les résultats sont comparables à d'autres études concernant des espèces proches. Néanmoins, la taille de ces populations approche des niveaux au-dessous desquels la survie à long terme est fortement incertaine. La persistance de la diversité génétique pour les prochaines générations reste en conséquence liée à la survie des adultes et à une reprise du succès de la reproduction. ABSTRACT Capercaillie (Tetrao urogallus) is a large grouse that is continuously distributed across the tundra and the mid-high mountains of Western Europe. However, the populations in Western Europe have been showing a constant decline during the last decades. The causes for this decline are possibly related to human activities, such as cattle breeding and tourism that have both led to habitat modification and fragmentation. Unfortunately, populations that have undergone drastic demographic bottlenecks often go through genetic processes of inbreeding and loss of diversity that decrease their fitness and eventually lead to extinction. This thesis presents the investigations conducted to estimate the impact of the demographic decline of capercaillie populations on the extent and distribution of their genetic variability in the Jura and in the Pyrenees mountains. Because grouse are protected by wildlife legislation, and also because of the cryptic behaviour of capercaillie, all DNA material used in this study was extracted from faeces (non-invasive genetic sampling). In the first part of my thesis, I detail the protocols of DNA extraction and PCR amplification adapted from classical methods using conventional DNA-rich samples. The use of faecal DNA imposes specific constraints due to the low quantity and the highly degraded genetic material available. These constraints are partially overcome by performing multiple genotyping repetitions to obtain sufficient reliability. I also investigate the causes of DNA degradation in faeces. Among the main degraders, namely bacterial activity, spontaneous hydrolysis, and free-¬DNase activities, the latter was pointed out as the most important according to our experiments. These enzymes degrade DNA very rapidly, and, as a consequence, faeces sampling schemes must be planned preferably in cold and dry weather conditions, allowing for enzyme activity inhibition. The second part of the thesis is a simulation study aiming to assess the capacity of the software Structure to detect population structure in hierarchical models relevant to situations encountered in wild populations, using several genetic markers. The methods implemented in Structure appear efficient in detecting the highest hierarchical structure. The third and fourth parts of the thesis describe the population genetics status of the remaining Jura and Pyrenees populations using 11 microsatellite loci. In either of these populations, no inbreeding nor reduced genetic diversity was detected. Furthermore, the genetic differentiation between patches defined by habitat suitability remains moderate and correlated with geographical distance, suggesting that significant dispersion between patches was at work at least until the last generations. The comparison of diversity indicators with other species or other populations of capercaillie indicate that population in the Jura has retained a large part of its original genetic diversity. These results suggest that the recent decline has had so forth a moderate impact on• genetic factors and that these populations might have retained the potential for long term survival, if the decline is stopped. Finally, in the fifth part, the analysis of relatedness between males participating in the reproduction parade, or lek, indicate that capercaillie males, like has been shown for some other grouse species, gather on leks• among individuals that are more related than the average of the population. This pattern appears to be due to both population structure and kin-association. As a conclusion, this first study relying exclusively on nuclear DNA extracted from faeces has provided novel information that was not available through field observation or classical genetic sampling. No bird has been captured or disturbed, and the results are consistent with other studies of closely related species. However, the size of these populations is approaching thresholds below which long-term survival is unlikely. The persistence of genetic diversity for the forthcoming generations remains therefore bond to adult survival and to the increase of reproduction success.

Altruism, spite and choosy females

Natural selection favours the genes which are able to introduce replicates of themselves in the next generation with higher certainty than do rival genes (Hamilton 1963). The fitness of an individual, it?s ability to produce future parents, depends on it?s own behaviour as well as on the behaviour of other individuals in the population. For instance, the intensity of competition an individual experience depends on the exploitation of resources by neighbours. The fitness is thus frequency dependent on what neighbours do. Behaviours can be classified according to the costs and benefits they have on the fitness of the behaver and it?s neighbours (Hamilton 1964, Hamilton 1975). According to this classification there exist four distinct social behaviours. (1) A gene confering the ability to use a new ressource is called selfish because it has a positive e_ect on the bearer of the gene but a negative e_ect on neighbours by the concomitant increase in competition. (2) An altruistic behaviour is defined as an action where an individual increases the fitness of a neighbour at the expense of it?s own. The e_ect is deleterious for the actor but positive for the receptor. (3) More surprinsingly, an individual might sacrifice a fraction of it?s ressources to harm another at no direct benefits. This spitefull behaviour incurs a cost for the actor but is also deleterious for the receptor. (4) Finally a cooperative behaviour breeds benefits for both actors and neighbours. In this thesis I will continue on the path traced by numerous evolutionnary biologist which attempt to fine tune our understanding of the evolution of social behaviours since Hamilton?s foundation (1963, 1964). A critical development over the last 40 years has been the realisation that competition between kin can partly or completely cancel out the role of relatedness as an agent favouring altruism (Wilson et al., 1992; Taylor, 1992a,b). Of importance is thus to determine the scale at which competition and altruism occur. One mechanism avoiding the complete dilution of relatedness by competition is the conditionnal expression of the social behaviors. Focus will be given in this thesis at the role played by di_erent recognition mechanism in paving the way to altruism (Komdeur and Hatchwell, 1999) when the population has a spatial structure. Further, the evolution of spite will also be considered in these settings. The thesis is fractionated into two parts. First, di_erent models promoting altruism cooperation and spite will be compared under the same theoretical umbrella. This is a rather informal and more personnal part of my thesis. It also serve as a justification and basis to "Altruism among kin and non-kin individuals" which is an article attempting to clas- sify the mechanisms leading to altruism and cooperation. Second, in the annexe, there are three research papers about kin selection, altruism and dispersal: "Is sociality driven by the costs of dispersal or the benefits of philopatry?: A role for kin-discrimination mechanism", "Altruism, dispersal and phenotype kin recognition" and "Inbreeding avoidance through kin recognition: choosy female boost male dispersal" this last paper incorporates kin recognition as an agent favoring sex-biased dispersal.

Global Diversity Lines-A Five-Continent Reference Panel of Sequenced Drosophila melanogaster Strains.

Reference collections of multiple Drosophila lines with accumulating collections of "omics" data have proven especially valuable for the study of population genetics and complex trait genetics. Here we present a description of a resource collection of 84 strains of Drosophila melanogaster whose genome sequences were obtained after 12 generations of full-sib inbreeding. The initial rationale for this resource was to foster development of a systems biology platform for modeling metabolic regulation by the use of natural polymorphisms as perturbations. As reference lines, they are amenable to repeated phenotypic measurements, and already a large collection of metabolic traits have been assayed. Another key feature of these strains is their widespread geographic origin, coming from Beijing, Ithaca, Netherlands, Tasmania, and Zimbabwe. After obtaining 12.5× coverage of paired-end Illumina sequence reads, SNP and indel calls were made with the GATK platform. Thorough quality control was enabled by deep sequencing one line to >100×, and single-nucleotide polymorphisms and indels were validated using ddRAD-sequencing as an orthogonal platform. In addition, a series of preliminary population genetic tests were performed with these single-nucleotide polymorphism data for assessment of data quality. We found 83 segregating inversions among the lines, and as expected these were especially abundant in the African sample. We anticipate that this will make a useful addition to the set of reference D. melanogaster strains, thanks to its geographic structuring and unusually high level of genetic diversity.

Diversity and genetic struture in cajá tree (spondias mombin l.) populations in northeastern Brazil

Spondias mombin L. is a fruit tree from the American continent from the Anacardiaceae family. In Brazil it is common in different vegetation types but is more frequent in the Atlantic and Amazonian rainforests. It is economically important because of its fruits, which are widely consumed raw or processed as fruit jellies, juices and ice creams. The leaves have great importance in the pharmaceutical industry because of their antibacterial properties. In the state of Pernambuco, cajá tree is widely distributed in the Zona da Mata region and less frequently in the Agreste and Sertão areas. In this work diversity and genetic structure were studied in four populations of cajá tree from Pernambuco's Zona da Mata, Northeast Brazil, using isozymes polymorphism analyses from electrophoreses. The result showed 100% of polymorphism (P) for nine alleles and the average of alleles per locus s was 2.4. The expected heterozygosity ranged from 0.530 to 0.574 and the observed heterozygosity , from 0.572 to 0.735. It was not observed inbreeding and the average F IT was -0.175, whereas within population inbreeding (f) varied from -0.08 to- 0.37. The genetic divergence among the populations (F ST) ranged from 0.006 to 0.028 and the average was 0.026. The average of estimated gene flow (Nm) was high (5.27). The CG-IPA population, corresponding to the germplasm collection of IPA, showed more than 96% of genetic similarity with other populations; therefore, it is a good representative of the existent genetic diversity in the Zona da Mata region.

Diversity and genetic structure in natural populations of Hancornia speciosa var. speciosa Gomes in northeastern Brazil

Hancornia speciosa Gomes is a fruit tree native from Brazil that belongs to Apocinaceae family, and is popularly known as Mangabeira. Its fruits are widely consumed raw or processed as fruit jam, juices and ice creams, which have made it a target of intense exploitation. The extractive activities and intense human activity on the environment of natural occurrence of H. speciosa has caused genetic erosion in the species and little is known about the ecology or genetic structure of natural populations. The objective of this research was the evaluation of the genetic diversity and genetic structure of H. speciosa var. speciosa. The genetic variability was assessed using 11 allozyme loci with a sample of 164 individuals distributed in six natural populations located in the States of Pernambuco and Alagoas, Northeastern Brazil. The results showed a high level of genetic diversity within the species (e= 0.36) seeing that the most of the genetic variability of H. speciosa var. speciosa is within its natural populations with low difference among populations ( or = 0.081). The inbreeding values within ( = -0.555) and among populations ( =-0.428) were low showing lacking of endogamy and a surplus of heterozygotes. The estimated gene flow ( m ) was high, ranging from 2.20 to 13.18, indicating to be enough to prevent the effects of genetic drift and genetic differentiation among populations. The multivariate analyses indicated that there is a relationship between genetic and geographical distances, which was confirmed by a spatial pattern analysis using Mantel test (r = 0.3598; p = 0.0920) with 1000 random permutations. The high genetic diversity index in these populations indicates potential for in situ genetic conservation.

Directional dominance on stature and cognition in diverse human populations.

Homozygosity has long been associated with rare, often devastating, Mendelian disorders, and Darwin was one of the first to recognize that inbreeding reduces evolutionary fitness. However, the effect of the more distant parental relatedness that is common in modern human populations is less well understood. Genomic data now allow us to investigate the effects of homozygosity on traits of public health importance by observing contiguous homozygous segments (runs of homozygosity), which are inferred to be homozygous along their complete length. Given the low levels of genome-wide homozygosity prevalent in most human populations, information is required on very large numbers of people to provide sufficient power. Here we use runs of homozygosity to study 16 health-related quantitative traits in 354,224 individuals from 102 cohorts, and find statistically significant associations between summed runs of homozygosity and four complex traits: height, forced expiratory lung volume in one second, general cognitive ability and educational attainment (P < 1 × 10(-300), 2.1 × 10(-6), 2.5 × 10(-10) and 1.8 × 10(-10), respectively). In each case, increased homozygosity was associated with decreased trait value, equivalent to the offspring of first cousins being 1.2 cm shorter and having 10 months' less education. Similar effect sizes were found across four continental groups and populations with different degrees of genome-wide homozygosity, providing evidence that homozygosity, rather than confounding, directly contributes to phenotypic variance. Contrary to earlier reports in substantially smaller samples, no evidence was seen of an influence of genome-wide homozygosity on blood pressure and low density lipoprotein cholesterol, or ten other cardio-metabolic traits. Since directional dominance is predicted for traits under directional evolutionary selection, this study provides evidence that increased stature and cognitive function have been positively selected in human evolution, whereas many important risk factors for late-onset complex diseases may not have been.