Genetic drift and within-host metapopulation dynamics of HIV-1 infection

Most HIV replication occurs in solid lymphoid tissue, which has prominent architecture at the histological level, which separates groups of productively infected CD4+ cells. Nevertheless, current population models of HIV assume panmixis within lymphoid tissue. We present a simple “metapopulation” model of HIV replication, where the population of infected cells is comprised of a large number of small populations, each of which is established by a few founder viruses and undergoes turnover. To test this model, we analyzed viral genetic variation of infected cell subpopulations within the spleen and demonstrated the action of founder effects as well as significant variation in the extent of genetic differentiation between subpopulations among patients. The combination of founder effects and subpopulation turnover can result in an effective population size much lower than the actual population size and may contribute to the importance of genetic drift in HIV evolution despite a large number of infected cells.

A pre-Columbian Y chromosome-specific transition and its implications for human evolutionary history.

A polymorphic C-->T transition located on the human Y chromosome was found by the systematic comparative sequencing of Y-specific sequence-tagged sites by denaturing high-performance liquid chromatography. The results of genotyping representative global indigenous populations indicate that the locus is polymorphic exclusively within the Western Hemisphere. The pre-Columbian T allele occurs at > 90% frequency within the native South and Central American populations examined, while its occurrence in North America is approximately 50%. Concomitant genotyping at the polymorphic tetranucleotide microsatellite DYS19 locus revealed that the C-->T mutation displayed significant linkage disequilibrium with the 186-bp allele. The data suggest a single origin of linguistically diverse native Americans with subsequent haplotype differentiation within radiating indigenous populations as well as post-Columbian European and African gene flow. The mutation may have originated either in North America at a very early time during the expansion or before it, in the ancestral population(s) from which all Americans may have originated. The analysis of linkage of the DYS199 and the DYS19 tetranucleotide loci suggests that the C-->T mutation may have occurred around 30,000 years ago. We estimate the nucleotide diversity over 4.2 kb of the nonrecombining portion of the Y chromosome to be 0.00014. compared to autosomes, the majority of variation is due to the smaller effective population size of the Y chromosome rather than selective sweeps. There begins to emerge a pattern of pronounced geographical localization of Y-specific nucleotide substitution polymorphisms.

Deriva genética de caracteres quantitativos em milho

A obtenção de genótipos superiores no melhoramento de plantas depende da existência de variabilidade genética. A existência de coleções de germoplasma representativas e a utilização de um tamanho adequado de amostra são fundamentais para a preservação das frequências alélicas e genotípicas, diminuindo a perda de variabilidade genética e postergando o aparecimento dos efeitos da deriva genética. Assim, teve-se como objetivo avaliar os efeitos da deriva genética em caracteres quantitativos em subpopulações de milho. Este estudo foi realizado a partir das populações originais BR-105 e BR-106, das quais 10 subpopulações foram obtidas em cada um dos cinco ciclos sucessivos de amostragem com tamanho efetivo reduzido, totalizando 50 subpopulações para cada população original, as quais foram posteriormente autofecundadas, gerando um nível a mais de endogamia. Os tratamentos foram constituídos de 10 amostras da população original sem autofecundação, 10 amostras com autofecundação, 50 subpopulações obtidas da população original e 50 subpopulações autofecundadas, totalizando 120 tratamentos para cada população, avaliados separadamente. Utilizou-se o delineamento em blocos casualizados no esquema de parcelas subdivididas em faixas hierárquico, em quatro ambientes com duas repetições por ambiente. Os caracteres avaliados foram produção de grãos (PG), prolificidade (PROL), comprimento e diâmetro de espigas (CE e DE), número de fileiras por espiga (NFE), número de grãos por fileira (NGF), altura de planta e espiga (AP e AE), florescimento masculino e feminino (FM e FF) e número de ramificações do pendão (NRP). Foram estimados os efeitos da deriva genética entre as médias das subpopulações nos dois níveis de endogamia e os efeitos da depressão por endogamia nas subpopulações dentro dos ciclos. Posteriormente, realizaram-se análises de regressão linear para as subpopulações nos dois níveis de endogamia, separadamente, e em conjunto. Foi verificada uma grande variação nas médias das subpopulações ao longo dos ciclos, indicando que a deriva genética causou diferenciação entre as mesmas e que estas se diferenciaram das populações originais. Detectaram-se efeitos significativos da deriva genética nas populações não autofecundadas para todos os caracteres avaliados, em maior número para PG, já que este caráter é mais sensível à deriva genética por possuir maior grau de dominância que os demais. Houve diminuição no número de estimativas de deriva significativas para as populações autofecundadas, incluindo mudanças na magnitude e no sinal das mesmas em relação às populações não autofecundadas. Para as estimativas de depressão por endogamia, os caracteres PG, NGF, FM e FF apresentaram maior quantidade de estimativas significativas que os demais. Para a maioria dos caracteres, a regressão linear explicou a maior parte da variação encontrada com o aumento dos coeficientes de endogamia. As populações BR-105 e BR-106, por terem estruturas genéticas distintas, apresentaram performances diferentes quanto aos efeitos da deriva genética. Enfim, como a deriva genética interfere na integridade genética das populações, torna-se importante considerar seus efeitos na coleta e manutenção dos bancos de germoplasma e nas populações utilizadas no melhoramento genético de plantas.

Genetic variation and inference of demographic histories in non-model species

Both long-term environmental changes such as those driven by the glacial cycles and more recent anthropogenic impacts have had major effects on the past demography in wild organisms. Within species, these changes are reflected in the amount and distribution of neutral genetic variation. In this thesis, mitochondrial and microsatellite DNA was analysed to investigate how environmental and anthropogenic factors have affected genetic diversity and structure in four ecologically different animal species. Paper I describes the post-glacial recolonisation history of the speckled-wood butterfly (Pararge aegeria) in Northern Europe. A decrease in genetic diversity with latitude and a marked population structure were uncovered, consistent with a hypothesis of repeated founder events during the postglacial recolonisation. Moreover, Approximate Bayesian Computation analyses indicate that the univoltine populations in Scandinavia and Finland originate from recolonisations along two routes, one on each side of the Baltic. Paper II aimed to investigate how past sea-level rises affected the population history of the convict surgeonfish (Acanthurus triostegus) in the Indo-Pacific. Assessment of the species’ demographic history suggested a population expansion that occurred approximately at the end of the last glaciation. Moreover, the results demonstrated an overall lack of phylogeographic structure, probably due to the high dispersal rates associated with the species’ pelagic larval stage. Populations at the species’ eastern range margin were significantly differentiated from other populations, which likely is a consequence of their geographic isolation. In Paper III, we assessed the effect of human impact on the genetic variation of European moose (Alces alces) in Sweden. Genetic analyses revealed a spatial structure with two genetic clusters, one in northern and one in southern Sweden, which were separated by a narrow transition zone. Moreover, demographic inference suggested a recent population bottleneck. The inferred timing of this bottleneck coincided with a known reduction in population size in the 19th and early 20th century due to high hunting pressure. In Paper IV, we examined the effect of an indirect but well-described human impact, via environmental toxic chemicals (PCBs), on the genetic variation of Eurasian otters (Lutra lutra) in Sweden. Genetic clustering assignment revealed differentiation between otters in northern and southern Sweden, but also in the Stockholm region. ABC analyses indicated a decrease in effective population size in both northern and southern Sweden. Moreover, comparative analyses of historical and contemporary samples demonstrated a more severe decline in genetic diversity in southern Sweden compared to northern Sweden, in agreement with the levels of PCBs found.

Bayesian inferences on the recent island colonization history by the bird Zosterops lateralis lateralis

The founding of new populations by small numbers of colonists has been considered a potentially important mechanism promoting evolutionary change in island populations. Colonizing species, such as members of the avian species complex Zosterops lateralis, have been used to support this idea. A large amount of background information on recent colonization history is available for one Zosterops subspecies, Z. lateralis lateralis, providing the opportunity to reconstruct the population dynamics of its colonization sequence. We used a Bayesian approach to combine historical and demographic information available on Z. l. lateralis with genotypic data from six microsatellite loci, and a rejection algorithm to make simultaneous inferences on the demographic parameters describing the recent colonization history of this subspecies in four southwest Pacific islands. Demographic models assuming mutation–drift equilibrium or a large number of founders were better supported than models assuming founder events for three of four recently colonized island populations. Posterior distributions of demographic parameters supported (i) a large stable effective population size of several thousands individuals with point estimates around 4000–5000; (ii) a founder event of very low intensity with a large effective number of founders around 150–200 individuals for each island in three of four islands, suggesting the colonization of those islands by one flock of large size or several flocks of average size; and (iii) a founder event of higher intensity on Norfolk Island with an effective number of founders around 20 individuals, suggesting colonization by a single flock of moderate size. Our inferences on demographic parameters, especially those on the number of founders, were relatively insensitive to the precise choice of prior distributions for microsatellite mutation processes and demographic parameters, suggesting that our analysis provides a robust description of the recent colonization history of the subspecies.

Limited effect of anthropogenic habitat fragmentation on molecular diversity in a rain forest skink, Gnypetoscincus queenslandiae

To examine the effects of recent habitat fragmentation, we assayed genetic diversity in a rain forest endemic lizard, the prickly forest skink (Gnypetoscincus queenslandiae), from seven forest fragments and five sites in continuous forest on the Atherton tableland of northeastern Queensland, Australia. The rain forest in this region was fragmented by logging and clearing for dairy farms in the early 1900s and most forest fragments studied have been isolated for 50-80 years or nine to 12 skink generations. We genotyped 411 individuals at nine microsatellite DNA loci and found fewer alleles per locus in prickly forest skinks from small rain forest fragments and a lower ratio of allele number to allele size range in forest fragments than in continuous forest, indicative of a decrease in effective population size. In contrast, and as expected for populations with small neighbourhood sizes, neither heterozygosity nor variance in allele size differed between fragments and sites in continuous forests. Considering measures of among population differentiation, there was no increase in F-ST among fragments and a significant isolation by distance pattern was identified across all 12 sites. However, the relationship between genetic (F-ST) and geographical distance was significantly stronger for continuous forest sites than for fragments, consistent with disruption of gene flow among the latter. The observed changes in genetic diversity within and among populations are small, but in the direction predicted by the theory of genetic erosion in recently fragmented populations. The results also illustrate the inherent difficulty in detecting genetic consequences of recent habitat fragmentation, even in genetically variable species, and especially when effective population size and dispersal rates are low.

Genetic analysis of complex demographic scenarios: Spatially expanding populations of the cane toad, Bufo marinus

Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection-regression algorithm. to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration-foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker. demographic, and historical features. The method presented here is the only robust, model-based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

High nuclear genetic diversity, high levels of outcrossing and low differentiation among remnant populations of Quercus petraea at the margin of its range in Ireland

Background and Aims Quercus petraea colonized Ireland after the last glaciation from refugia on mainland Europe. Deforestation. however. beginning in Neolithic times, has resulted in small, scattered forest fragments, now covering less than 12 000 ha. Methods Plastid (three fragments) and microsatellite variation (13 loci) were characterized in seven Irish populations sampled along a north-south gradient. Using Bayesian approaches and Wright's F-statistics, the effects of colonization and fragmentation on the genetic structure and mating patterns of extant oak populations were investigated. Key-Results All Populations possessed cytotypes common to the Iberian Peninsula. Despite the distance from the refugial core and the extensive deforestation in Ireland, nuclear genetic variation was high and comparable to mainland Europe. Low population differentiation was observed within Ireland and populations showed no evidence for isolation by distance. As expected of a marker with an effective Population size of one-quarter relative to the nuclear genome, plastid variation indicated higher differentiation. Individual inbreeding coefficients indicated high levels of outcrossing. Conclusions Consistent with a large effective Population size in the historical migrant gene pool and/or with high gene flow among populations, high within-population diversity and low population differentiation was observred within Ireland. It is proposed that native Q. petraea populations in Ireland share a common phylogeographic history and that the present genetic structure does not reflect founder effects. (C) 2004 Annals of Botany Company.

Sociality and the rate of molecular evolution

The molecular clock does not tick at a uniform rate in all taxa but maybe influenced by species characteristics. Eusocial species (those with reproductive division of labor) have been predicted to have faster rates of molecular evolution than their nonsocial relatives because of greatly reduced effective population size; if most individuals in a population are nonreproductive and only one or few queens produce all the offspring, then eusocial animals could have much lower effective population sizes than their solitary relatives, which should increase the rate of substitution of nearly neutral mutations. An earlier study reported faster rates in eusocial honeybees and vespid wasps but failed to correct for phylogenetic nonindependence or to distinguish between potential causes of rate variation. Because sociality has evolved independently in many different lineages, it is possible to conduct a more wide-ranging study to test the generality of the relationship. We have conducted a comparative analysis of 25 phylogenetically independent pairs of social lineages and their nonsocial relatives, including bees, wasps, ants, termites, shrimps, and mole rats, using a range of available DNA sequences (mitochondrial and nuclear DNA coding for proteins and RNAs, and nontranslated sequences). By including a wide range of social taxa, we were able to test whether there is a general influence of sociality on rates of molecular evolution and to test specific predictions of the hypothesis: (1) that social species have faster rates because they have reduced effective population sizes; (2) that mitochondrial genes would show a greater effect of sociality than nuclear genes; and (3) that rates of molecular evolution should be correlated with the degree of sociality. We find no consistent pattern in rates of molecular evolution between social and nonsocial lineages and no evidence that mitochondrial genes show faster rates in social taxa. However, we show that the most highly eusocial Hymenoptera do have faster rates than their nonsocial relatives. We also find that social parasites (that utilize the workers from related species to produce their own offspring) have faster rates than their social relatives, which is consistent with an effect of lower effective population size on rate of molecular evolution. Our results illustrate the importance of allowing for phylogenetic nonindependence when conducting investigations of determinants of variation in rate of molecular evolution.

A multilocus perspective on refugial isolation and divergence in rainforest skinks (Carlia)

To explore the evolutionary consequences of climate-induced fluctuations in distribution of rainforest habitat we contrasted demographic histories of divergence among three lineages of Australian rainforest endemic skinks. The red-throated rainbow skink, Carlia rubrigularis, consists of morphologically indistinguishable northern and southern mitochondrial DNA (mtDNA) lineages that are partially reproductively isolated at their parapatric boundary. The third lineage (C. rhomboidalis) inhabits rainforests just to the south of C. rubrigularis, has blue, rather than red-throated males, and for mtDNA is more closely related to southern C. rubrigularis than is northern C. rubrigularis. Multigene coalescent analyses supported more recent divergence between morphologically distinct lineages than between morphologically conservative lineages. There was effectively no migration and therefore stronger isolation between southern C. rubrigularis and C. rhomboidalis, and low unidirectional migration between morphologically conservative lineages of C. rubrigularis. We found little or no evidence for strong differences in effective population size, and hence different contributions of genetic drift in the demographic history of the three lineages. Overall the results suggest contrasting responses to long-term fluctuations in rainforest habitats, leading to varying opportunities for speciation.

The genetic structure of Australasian green turtles (Chelonia mydas): exploring the geographical scale of genetic exchange

Ecological and genetic studies of marine turtles generally support the hypothesis of natal homing, but leave open the question of the geographical scale of genetic exchange and the capacity of turtles to shift breeding sites. Here we combine analyses of mitochondrial DNA (mtDNA) variation and recapture data to assess the geographical scale of individual breeding populations and the distribution of such populations through Australasia. We conducted multiscale assessments of mtDNA variation among 714 samples from 27 green turtle rookeries and of adult female dispersal among nesting sites in eastern Australia. Many of these rookeries are on shelves that were flooded by rising sea levels less than 10 000 years (c. 450 generations) ago. Analyses of sequence variation among the mtDNA control region revealed 25 haplotypes, and their frequency distributions indicated 17 genetically distinct breeding stocks (Management Units) consisting either of individual rookeries or groups of rookeries in general that are separated by more than 500 km. The population structure inferred from mtDNA was consistent with the scale of movements observed in long-term mark-recapture studies of east Australian rookeries. Phylogenetic analysis of the haplotypes revealed five clades with significant partitioning of sequence diversity (Phi = 68.4) between Pacific Ocean and Southeast Asian/Indian Ocean rookeries. Isolation by distance was indicated for rookeries separated by up to 2000 km but explained only 12% of the genetic structure. The emerging general picture is one of dynamic population structure influenced by the capacity of females to relocate among proximal breeding sites, although this may be conditional on large population sizes as existed historically across this region.

Genetic structure and male-mediated gene flow in the ghost bat (Macroderma gigas)

The Australian ghost bat is a large, opportunistic carnivorous species that has undergone a marked range contraction toward more mesic, tropical sites over the past century. Comparison of mitochondrial DNA (mtDNA) control region sequences and six nuclear microsatellite loci in 217 ghost bats from nine populations across subtropical and tropical Australia revealed strong population subdivision (mtDNA phi(ST) = 0.80; microsatellites URST = 0.337). Low-latitude (tropical) populations had higher heterozygosity and less marked phylogeographic structure and lower subdivision among sites within regions (within Northern Territory [NT] and within North Queensland [NQ]) than did populations at higher latitudes (subtropical sites; central Queensland [CQ]), although sampling of geographically proximal breeding sites is unavoidably restricted for the latter. Gene flow among populations within each of the northern regions appears to be male biased in that the difference in population subdivision for mtDNA and microsatellites (NT phi(ST) = 0.39, URST = 0.02; NQ phi(ST) = 0.60, URST = -0.03) is greater than expected from differences in the effective population size of haploid versus diploid loci. The high level of population subdivision across the range of the ghost bat contrasts with evidence for high gene flow in other chiropteran species and may be due to narrow physiological tolerances and consequent limited availability of roosts for ghost bats, particularly across the subtropical and relatively arid regions. This observation is consistent with the hypothesis that the contraction of the species' range is associated with late Holocene climate change. The extreme isolation among higher-latitude populations may predispose them to additional local extinctions if the processes responsible for the range contraction continue to operate.

Biogeographical distribution of Rimicaris exoculata resident gut epibiont communities along the Mid-Atlantic Ridge hydrothermal vent sites

Rimicaris exoculata is a deep-sea hydrothermal vent shrimp which enlarged gill chamber houses a complex trophic epibiotic community. Its gut harbours an autochthonous and distinct microbial community. This species dominates hydrothermal ecosystems megafauna along the Mid-Atlantic Ridge, regardless of contrasted geochemical conditions prevailing in them. Here, the resident gut epibiont community at four contrasted hydrothermal vent sites (Rainbow/TAG/Logatchev/Ashadze) was analysed and compiled with previous data to evaluate the possible influence of site location, using 16S rRNA surveys and microscopic observations (TEM, SEM and FISH analyses). Filamentous epibionts inserted between the epithelial cells microvilli were observed on all examined samples. Results confirmed resident gut community affiliation to Deferribacteres, Mollicutes, Epsilonproteobacteria and to a lesser extent Gammaproteobacteria lineages. Still a single Deferribacteres phylotype was retrieved at all sites. Four Mollicutes-related OTUs were distinguished, one being only identified on Rainbow specimens. The topology of ribotypes median-joining networks illustrated a community diversification possibly following demographic expansions, suggesting a more ancient evolutionary history and/or a larger effective population size at Rainbow. Finally, the gill chamber community distribution was also analysed through ribotypes networks based on sequences from R. exoculata collected at Rainbow/Snake Pit/TAG/Logatchev/Ashadze sites. Results allow refining hypotheses on the epibiont role and transmission pathways.

Phylodynamic modelling of foot-and-mouth disease virus sequence data

The under-reporting of cases of infectious diseases is a substantial impediment to the control and management of infectious diseases in both epidemic and endemic contexts. Information about infectious disease dynamics can be recovered from sequence data using time-varying coalescent approaches, and phylodynamic models have been developed in order to reconstruct demographic changes of the numbers of infected hosts through time. In this study I have demonstrated the general concordance between empirically observed epidemiological incidence data and viral demography inferred through analysis of foot-and-mouth disease virus VP1 coding sequences belonging to the CATHAY topotype over large temporal and spatial scales. However a more precise and robust relationship between the effective population size (

Aplicação de microssatélites na crioconservação da diversidade genética do cavalo lusitano

O género Equus teve origem na América do Norte e alguns exemplares migraram para a Eurásia pelo Estreito de Bering, durante a última glaciação. No fim da glaciação, todos os cavalos do continente americano extinguiram-se, mas sobreviveram nas estepes da Eurásia, na Peninsula lbérica e nas florestas da Europa Ocidental e Central. O cavalo Lusitano teve a sua origem em cavalos selvagens e domesticados da Peninsula lbérica, ocorrendo uma mistura com outros animais trazidos por eventos migratórios ocorridos no passado. Os cavalos deste gene pool contribuiram para o desenvolvimento de outras raças modernas na Europa e foram mais tarde introduzidos e dispersos pelo continente Americano, tornando-se fundadores de numerosas raças do novo mundo. A raça Lusitana é uma raça equina autóctone portuguesa, com especial relevancia económica no panorama nacional e internacional. Apesar de não ser uma raça ameaçada, alguns autores defendem que a informação genealógica disponivel (pedigrees) indica que uma utilização excessiva de um reduzido número de reprodutores machos esta a diminuir a diversidade genética da raça, tendo como consequência o aumento da consanguinidade e a diminuição do tamanho efetivo da população para cerca de metade dos valores recomendados pela FAO. No entanto, a anàlise da diversidade genética com base em 16 microssatélites (Marcadores de DNA) a um grupo de 2699 machos da raça Lusitana, nascidos entre 1985 e 2010 e inscritos como reprodutores no Livro Genealógico da raça, revelou um elevado nível de diversidade, idêntico ao encontrado na maioria das raças equinas. Dada a crescente relevância da Crioconservação, omo estratégia complementar para a conservação da diversidade genética in situ, e tendo em conta que não existe criopreservação de oocitos, embriões ou sémen, do cavalo de raça Lusitana em Banco de Genes, selecionaram-se 62 machos reprodutores (garanhões) com interesse genético para a criopreservação de sémen, quer no sentido de preservar a diversidade da raça quer no da salvaguardar em caso de calamidade; ABSTRACT: The genus Equus originated in North America and some exemplary migrated to Eurasia through the Bering Strait during the last glaciation. By the end of the last glaciation, all horses on the American continent became extinct but the genus survived in the steppes of Eurasia, in the Iberian Peninsula and on the Central and West Europe forests. The Lusitano horse breed has its origins in wild and domesticated horses of the Iberian Peninsula, where a mixture with other animals brought by migratory events in the past occurred. The horses of this gene pool contributed to the development of other modern breeds in Europe and were later introduced and spread throughout the American continent, becoming founders of numerous breeds of the New World. The Lusitano horse breed, is a Portuguese native equine breed, with special economic relevance in the national and international scene. Although not being an endangered breed, some authors argue that the available genealogical information (pedigrees) indicates that an excessive use of a limited number of stallions is decreasing the genetic diversity of the breed, resulting in the increase of inbreeding and on the decrease of the effective population size to about half of the values recommended by FAO. However, the analysis of genetic diversity based on 16 microsatellites (DNA markers) in a group of 2699 males of the Lusitano horse breed, born between 1985 and 2010 and registered as Stallions in the Studbook, revealed a high level of diversity similar to that found in the majority of equine breeds. Given the growing relevance of Cryopreservation as a complementary strategy for the conservation of genetic diversity in situ and, taking into consideration the inexistence of criopreservation for oocytes, embryos and semen, in a Gene Bank, for the Lusitano horse breed, 62 breeding males (stallions) with genetic interest for semen cryopreservation were selected in order either to preserve the diversity of the breed or as safeguard in case of calamity.