Estimating coancestry within open-pollinated progenies of a dioecious species: the case study of Myracrodruon urundeuva

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Genome-Wiide estimates of coancestry, in breeding and effective population size in the Spanish Holstein population

Estimates of effective population size in the Holstein cattle breed have usually been low despite the large number of animals that constitute this breed. Effective population size is inversely related to the rates at which coancestry and inbreeding increase and these rates have been high as a consequence of intense and accurate selection. Traditionally, coancestry and inbreeding coefficients have been calculated from pedigree data. However, the development of genome-wide single nucleotide polymorphisms has increased the interest of calculating these coefficients from molecular data in order to improve their accuracy. In this study, genomic estimates of coancestry, inbreeding and effective population size were obtained in the Spanish Holstein population and then compared with pedigree-based estimates. A total of 11,135 animals genotyped with the Illumina BovineSNP50 BeadChip were available for the study. After applying filtering criteria, the final genomic dataset included 36,693 autosomal SNPs and 10,569 animals. Pedigree data from those genotyped animals included 31,203 animals. These individuals represented only the last five generations in order to homogenise the amount of pedigree information across animals. Genomic estimates of coancestry and inbreeding were obtained from identity by descent segments (coancestry) or runs of homozygosity (inbreeding). The results indicate that the percentage of variance of pedigree-based coancestry estimates explained by genomic coancestry estimates was higher than that for inbreeding. Estimates of effective population size obtained from genome-wide and pedigree information were consistent and ranged from about 66 to 79. These low values emphasize the need of controlling the rate of increase of coancestry and inbreeding in Holstein selection programmes.

NEESTIMATOR v2: re-implementation of software for the estimation of contemporary effective population size (N-e) from genetic data

NeEstimator v2 is a completely revised and updated implementation of software that produces estimates of contemporary effective population size, using several different methods and a single input file. NeEstimator v2 includes three single-sample estimators (updated versions of the linkage disequilibrium and heterozygote-excess methods, and a new method based on molecular coancestry), as well as the two-sample (moment-based temporal) method. New features include the following: (i) an improved method for accounting for missing data; (ii) options for screening out rare alleles; (iii) confidence intervals for all methods; (iv) the ability to analyse data sets with large numbers of genetic markers (10000 or more); (v) options for batch processing large numbers of different data sets, which will facilitate cross-method comparisons using simulated data; and (vi) correction for temporal estimates when individuals sampled are not removed from the population (Plan I sampling). The user is given considerable control over input data and composition, and format of output files. The freely available software has a new JAVA interface and runs under MacOS, Linux and Windows.

Population structure and evolution in the ant Plagiolepis pygmaea and its two social parasites Plagiolepis xene and Plagiolepis grassei

Social behaviour affects dispersal of animals and is an important modifier of genetic population structures. The female sex is often philopatric, which maintains coancestry within the breeding groups and promotes cooperative behaviours. This enables also inclusive fitness returns from altruism and explains why some individuals sacrifice personal reproduction for the good of others in social insects such as ants. However, reduced dispersal and population substructuring at the level of colonies may also entail inbreeding, loss of genetic diversity, and vulnerability. In addition, the most vulnerable ants are species that are evolved to parasitize colonies of other ants, and which compromise between abilities to disperse and the efficiency to parasitize the host. On the other hand, certain social organisations of ant colonies may facilitate a species to disperse outside its natural range and become a pest. Altogether, knowledge on genetic structuring of ant populations, as well as the evolution of their life histories can contribute to conservation biology and population management. The aim of this thesis was to investigate population structures and phylogenetic evolution of the ant Plagiolepis pygmaea and its two obligatory, workerless social parasites (inquilines) P. xene and P. grassei with genetic markers and DNA sequence data. The results support the general assumption that populations of inquiline parasites are highly fragmented and genetically vulnerable. Comparison of the two parasites suggests that differences in their relative abundance may follow from their interaction with the host, i.e. how well the species is adapted to reproduce in the host colonies. The results also indicate that the most recent free living ancestor to these two parasite species is their common host. This is considered to provide evidence for the controversial issue of sympatric speciation. Further, given that the level of adaptations to parasitic life history depends on the evolutionary time since the free-living ancestor, the results establish a link between species rarity and its evolutionary age. The populations of the host species P. pygmaea displayed significantly reduced dispersal both among the females (queens) and males, and high levels of inbreeding which may enhance worker altruism. In addition, the queens were found to mate with multiple males. Given the high relatedness between the queens and their mates, this occurs probably for non-genetic reasons, e.g. without benefits associated in genetically more diverse offspring. The results hence caution that the contribution of non-genetic factors to the prevailing mating patterns and genetic population structures should not be underestimated.

Relatedness testing in subdivided populations

Assigning probabilities to alleged relationships, given DNA profiles, requires, among other things, calculation of a likelihood ratio (LR). Such calculations usually assume independence of genes: this assumption is not appropriate when the tested individuals share recent ancestry due to population substructure. Adjusted LR formulae, incorporating the coancestry coefficient F(ST), are presented here for various two-person relationships, and the issue of mutations in parentage testing is also addressed.

Implications for DNA identification arising from an analysis of Australian forensic databases

Previous analyses of Australian samples have suggested that populations of the same broad racial group (Caucasian, Asian, Aboriginal) tend to be genetically similar across states. This suggests that a single national Australian database for each such group may be feasible, which would greatly facilitate casework. We have investigated samples drawn from each of these groups in different Australian states, and have quantified the genetic homogeneity across states within each racial group in terms of the "coancestry coefficient" F(ST). In accord with earlier results, we find that F(ST) values, as estimated from these data, are very small for Caucasians and Asians, usually <0.5%. We find that "declared" Aborigines (which includes many with partly Aboriginal genetic heritage) are also genetically similar across states, although they display some differentiation from a "pure" Aboriginal population (almost entirely of Aboriginal genetic heritage).

Mating system of a population of Myracrodruon urundeuva F.F. & M.F. Allemão using the fAFLP molecular marker

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Optimum contribution selection using differential evolution

Um programa baseado na técnica de evolução diferencial foi desenvolvido para a definição da contribuição genética ótima na seleção de candidatos a reprodução. A função- objetivo a ser otimizada foi composta pelo mérito genético esperado da futura progênie e pela coascendência média dos animais em reprodução. Conjuntos de dados reais e simulados de populações com gerações sobrepostas foram usados para validar e testar o desempenho do programa desenvolvido. O programa se mostrou computacionalmente eficiente e viável para ser aplicado na prática e as consequências esperadas de sua aplicação, em comparação a procedimentos empíricos de controle da endogamia e/ou com a seleção baseada apenas no valor genético esperado, seriam a melhora da resposta genética futura e limitação mais efetiva da taxa de endogamia.

Sistema de reprodução em pequenas populações fragmentadas e em árvores isoladas de Hymenaea stigonocarpa

Six microsatellite loci were used to quantify the mating system of two small fragmented populations (Selviria - SEL and Aparecida do Tabuado APT, Mato Grosso do Sul State) and isolated trees in pastures, of the bat-pollinated tropical tree Hymenaea stignocarpa, growing in the Center-west region of Brazil. In SEL population, seeds were collected from 11 mother-trees; in APT, from three trees and, in the case of isolated trees, from six individuals growing at least 500 m apart in pastures. To investigate if there are differences on mating system between trees in populations and isolated trees, trees from populations were pooled as a group and, likewise, the isolated trees were pooled to another group. The outcrossing rate was higher in the populations ((t) over cap (m)=0.873) than in isolated trees ((t) over cap (m)=0.857), but the difference was not significant. Significant and high differences between multi-locus and single-locus outcrossing rate were detected in populations ((t) over cap (m)-(t) over cap (s)=0.301, P<0.05) and isolated trees (<(t)over cap>(m)-(t) over cap (s) = 0.276, P < 0.05), suggesting mating between relatives. Higher paternity correlation was observed in trees from population (<(r)over cap>(p)=0.636) than in isolated trees ((r) over cap (p)=0.377), indicating the occurrence of some correlated matings and that part of offspring are full-sibs. It was not observed increased in self-fertilization rate in isolated trees in pastures. In general terms, the unique observed difference in mating system between populations and isolate trees was the high rate of correlated matings in trees from populations, due probably to the small distance among coespecifics and the pollinator behavior, visiting near trees.

Diversidade genética e tamanho efetivo de duas populações de Myracrodruon urundeuva Fr. All., sob conservação ex situ

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Mating system in Myracrodruon urundeuva (Anarcadiaceae): implications for conservation genetics

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Contemporary pollen flow, mating patterns and effective population size inferred from paternity analysis in a small fragmented population of the Neotropical tree Copaifera langsdorffii Desf. (Leguminosae-Caesalpinioideae)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Correlated matings in dioecious tropical tree, Myracrodruon urundeuva fr. all.

The correlated matings in two populations (Selvíria - SEL and Paulo de Faria - PFA) of dioecious Myracrodruon urundeuva were studied in the Southwest of Brazil, by allozyme analysis of progeny arrays using the sibling-pair model. Open-pollinated seeds were collected from 25 to 30 trees within populations. Departure from random matings were evident from the differences in pollen and ovule allele frequencies. The high and significant correlation of paternity (SEL r̂p = 0.671 ± 0.074; PFA r̂p = 0.371 ± 0.062) and a low number of effective pollinating trees (ranging from 2 to 3) were detected in the populations, suggesting high proportion of full-sibs progenies. According to these results, the estimate of coancestry within families (θ = 0.209 - SEL; θ = 0.171 - PFA.) exceeded the expectation of the half-sib progenies (θ = 0.125). Result outcomes are discussed from a conservation and breeding point of view.

Sistema de reprodução em pequenas populações fragmentadas e em árvores isoladas de Hymenaea stigonocarpa

Six microsatellite loci were used to quantify the mating system of two small fragmented populations (Selviria - SEL and Aparecida do Tabuado - APT, Mato Grosso do Sul State) and isolated trees in pastures, of the bat-pollinated tropical tree Hymenaea stignocarpa, growing in the Center-west region of Brazil. In SEL population, seeds were collected from 11 mother-trees; in APT, from three trees and, in the case of isolated trees, from six individuals growing at least 500 m apart in pastures. To investigate if there are differences on mating system between trees in populations and isolated trees, trees from populations were pooled as a group and, likewise, the isolated trees were pooled to another group. The outcrossing rate was higher in the populations (t̂ m= 0.873) than in isolated trees (t̂ m=0.857), but the difference was not significant. Significant and high differences between multi-locus and single-locus outcrossing rate were detected in populations (t̂ m- t̂ s=0.301, P<0.05) and isolated trees (t̂ m- t̂ s=0.276, P<0.05), suggesting mating between relatives. Higher paternity correlation was observed in trees from population (r̂ p=0.636) than in isolated trees (r̂ p=0.377), indicating the occurrence of some correlated matings and that part of offspring are full-sibs. It was not observed increased in self-fertilization rate in isolated trees in pastures. In general terms, the unique observed difference in mating system between populations and isolate trees was the high rate of correlated matings in trees from populations, due probably to the small distance among coespecifics and the pollinator behavior, visiting near trees.

Sistema de reprodução, parentesco e tamanho efetivo em sementes de polinização livre de populações fragmentadas de Copaifera langsdorffii Desf. por análise de locos microssatélites

Anthropogenic changes in the reproductive population density can affect the mating system and result in an increase in selfing and correlated matings. This study investigated the mating system in small fragmented populations of the insect pollinated tropical tree species Copaifera langsdorffii, using microsatellite loci and the mixed and correlated mating models. Open-pollinated seeds were collected from 15 seed-trees located in a small forest fragment (4.8 ha), denominated Bosque and from 14 other seed-trees located in other small forest fragments of the north-western region of São Paulo State. No significant differences were observed between the seed-trees from Bosque (tm=0.933±0.028) and other fragments (tm=0.971±0.032), although these estimates were significantly different from 1.0, suggesting that selfing was occurring. Differences between multilocus and unilocus outcrossing rate were significantly high in both seed-trees of Bosque (tm -ts=0.478±0.05) and other forest fragments (tm -ts=0.475±0.018), suggesting a spatial genetic structure in those stands. The results also showed high rates of correlated mating in the samples, indicating that a good part of the offspring were full-sibs. As a consequence of selfing, mating among relatives and correlated matings, the coancestry within families was equally high in the seed-trees of Bosque ((Θ=0.237) and in the seed-trees of the other forest fragments (Θ=0.241) and the effective population size was lower than expected in panmitic populations (Ne<4). The results were discussed, focusing on the sample size of seed-trees to collect seeds for genetic conservation and enviromental reforestation.