Inbreeding and genetic structure in the endangered Sorraia horse breed: implications for its conservation and management

The Sorraia horse is a closed breed with reduced effective population size and considered in critical maintained risk status. The breed exists in 2 main breeding populations, one in Portugal and one in Germany, with a smaller population size. A set of 22 microsatellite loci was used to examine genetic diversity and structure of the Sorraia horse breed and to compare individual inbreeding coefficient F, estimated from pedigree data, with individual heterozygosity and mean d(2). The Sorraia horse shows lower levels of microsatellite diversity when compared with other horse breeds. Due to management strategies, there are clear differences in the genetic structure of the 2 main Sorraia horse populations. Individual heterozygosity was shown to be a good estimator, used together with or as an alternative to inbreeding coefficient, in predicting fitness and evaluating the inbreeding level of the Sorraia horse. The information gathered in this study, combined with information available from previous studies, offers an important and wide information base for the future development of an effective breeding management of the Sorraia horse in order to preserve this endangered breed.     

Microsatellite diversity, population subdivision and gene flow in the Lipizzan horse

Blood samples of 561 Lipizzan horses from subpopulations (studs) of seven European countries representing a large fraction of the breed's population were used to examine the genetic diversity, population subdivision and gene flow in the breed. DNA analysis based on 18 microsatellite loci revealed that genetic diversity (observed heterozygosity = 0.663, gene diversity = 0.675 and the mean number of alleles = 7.056) in the Lipizzan horse is similar to other horse breeds as well as to other domestic animal species. The genetic differentiation between Lipizzan horses from different studs, although moderate, was apparent (pairwise F(ST) coefficients ranged from 0.021 to 0.080). Complementary findings explaining the genetic relationship among studs were revealed by genetic distance and principal component analysis. One genetic cluster consisted of the subpopulations of Austria, Italy and Slovenia, which represent the classical pool of Lipizzan horse breeding. A second cluster was formed by the Croatian, Hungarian and Slovakian subpopulations. The Romanian subpopulation formed a separate unit. The largest genetic differentiation was found between the Romanian and Italian subpopulation. Genetic results are consistent with the known breeding history of the Lipizzan horse. Correct stud assignment was obtained for 80.9% and 92.1% of Lipizzan horses depending on the inclusion or exclusion of migrant horses, respectively. The results of the present study will be useful for the development of breeding strategies, which consider classical horse breeding as well as recent achievements of population and conservation genetics.     

Genetic analysis, breed assignment and conservation priorities of three native Danish horse breeds

A genetic analysis was performed on three indigenous Danish horse breeds using 12 microsatellite markers from a standard kit for parental testing. These three breeds are all considered endangered based on their small population sizes. Genetic variation in these three breeds was comparable to other horse breeds in Europe, and they do not seem to be at immediate danger of extinction caused by genetic deterioration. The Knabstrupper breed had more genetic variation, as measured by expected heterozygosity and allelic richness, than the other two breeds (Frederiksborg and Jutland). F(ST) statistics and population assignments confirmed population differentiation into three distinct breeds. The Frederiksborg and Knabstrupper breeds were closer to each other than to the Jutland breed. When establishing conservation priorities for the breeds, the priorities will depend on the conservation goals. Different methods for establishing conservation priorities are also discussed.     

Semen parameters and level of microsatellite heterozygosity in Noriker draught horse stallions

It was the aim of the present study to determine physiological values for different semen parameters in an endangered draught horse breed, the Austrian Noriker. Because small population size is often believed to cause a decrease in fertility and/or semen quality through inbreeding and a reduction in genetic variation, the general genomic heterogeneity of the breed was estimated on the basis of microsatellite variation and correlated to semen parameters. Semen could be collected from 104 of 139 stallions with semen collection being more often successful in younger stallions. Mean volume of ejaculates was 90.8+/-55.1 ml, density 243+/-114 x 10(6)ml(-1), total sperm count 21.0+/-23.7 x 10(9), percentage of morphologically normal spermatozoa 38+/-18% and total motility 50+/-23%. Total sperm count and semen motility were significantly affected by age. Blood samples of 134 stallions were analysed for 12 microsatellite DNA markers. Genotypes of 110 stallions with at least 11 successfully typed markers were used for calculation of heterozygosity. A total of 82 alleles was identified with a mean of 6.8 alleles per marker. Heterozygosity varied between 35 and 76% for the different markers, mean heterozygosity was calculated to 63%. No correlation between heterozygosity and semen parameters was found.     

Population structure of an endangered species living in contrasted habitats: Parnassia palustris (Saxifragaceae)

In endangered species, it is critical to analyse the level at which populations interact (i.e. dispersal) as well as the levels of inbreeding and local adaptation to set up conservation policies. These parameters were investigated in the endangered species Parnassia palustris living in contrasted habitats. We analysed population structure in 14 populations of northern France for isozymes, cpDNA markers and phenotypic traits related to fitness. Within population genetic diversity and inbreeding coefficients were not correlated to population size. Populations seem not to have undergone severe recent bottleneck. Conversely to pollen migration, seed migration seems limited at a regional scale, which could prevent colonization of new sites even if suitable habitats appear. Finally, the habitat type affects neither within-population genetic diversity nor genetic and phenotypic differentiation among populations. Thus, even if unnoticed local adaptation to habitats exists, it does not influence gene flow between populations.     

Inbreeding,microsatellite heterozygosity,and morphological traits in Lipizzan horses

While the negative effects of inbreeding and reduced heterozygosity on fecundity and survival are well established, only a few investigations have been carried out concerning their influence on morphological traits. This topic is of particular interest for a small and closed population such as the Lipizzan horse. Thus, 27 morphological traits were measured in 360 Lipizzan mares and were regressed on the individual inbreeding coefficients, as well as on the individual heterozygosity and mean squared distances (mean d(2)) between microsatellite alleles within an individual. Both individual heterozygosity and mean d(2) were based on 17 microsatellite loci dispersed over 14 chromosomes. The results obtained by multivariate analysis reveal significant effects of stud (P <.0001), age at measurement (P <.0001), and mean d(2) (P =.0143). In univariate analyses, significant associations were obtained between length of pastern-hindlimbs and inbreeding coefficient (P <.01), length of cannons-hindlimb and mean d(2) (P <.01), and length of neck and mean d(2) (P <.001). After adjustment of single-test P values for multiple tests (Hochberg's step-up Bonferroni method), only the association of the length of neck and mean d(2) remained significant (P =.0213). Thus, no overall large effects of inbreeding, microsatellite heterozygosity, and mean d(2) on morphological traits were observed in the Lipizzan horse.     

Bayesian prediction of breeding values for multivariate binary and continuous traits in simulated horse populations using threshold-linear models with Gibbs sampling

Simulated data were used to determine the properties of multivariate prediction of breeding values for categorical and continuous traits using phenotypic, molecular genetic and pedigree information by mixed linear-threshold animal models via Gibbs sampling. Simulation parameters were chosen such that the data resembled situations encountered in Warmblood horse populations. Genetic evaluation was performed in the context of the radiographic findings in the equine limbs. The simulated pedigree comprised seven generations and 40 000 animals per generation. The simulated data included additive genetic values, residuals and fixed effects for one continuous trait and liabilities of four binary traits. For one of the binary traits, quantitative trait locus (QTL) effects and genetic markers were simulated, with three different scenarios with respect to recombination rate (r) between genetic markers and QTL and polymorphism information content (PIC) of genetic markers being studied: r = 0.00 and PIC = 0.90 (r0p9), r = 0.01 and PIC = 0.90 (r1p9), and r = 0.00 and PIC = 0.70 (r0p7). For each scenario, 10 replicates were sampled from the simulated horse population, and six different data sets were generated per replicate. Data sets differed in number and distribution of animals with trait records and the availability of genetic marker information. Breeding values were predicted via Gibbs sampling using a Bayesian mixed linear-threshold animal model with residual covariances fixed to zero and a proper prior for the genetic covariance matrix. Relative breeding values were used to investigate expected response to multi- and single-trait selection. In the sires with 10 or more offspring with trait information, correlations between true and predicted breeding values ranged between 0.89 and 0.94 for the continuous traits and between 0.39 and 0.77 for the binary traits. Proportions of successful identification of sires of average, favourable and unfavourable genetic value were 81% to 86% for the continuous trait and 57% to 74% for the binary traits in these sires. Expected decrease of prevalence of the QTL trait was 3% to 12% after multi-trait selection for all binary traits and 9% to 17% after single-trait selection for the QTL trait. The combined use of phenotype and genotype data was superior to the use of phenotype data alone. It was concluded that information on phenotypes and highly informative genetic markers should be used for prediction of breeding values in mixed linear-threshold animal models via Gibbs sampling to achieve maximum reduction in prevalences of binary traits.     

Genetic variation in two conserved local Romanian pig breeds using type 1 DNA markers

Analysis of the genetic variation of an endangered population is an important component for the success of conservation. Animals from two local Romanian pig breeds, the Mangalitsa and Bazna, were analysed for variation at a number of genetic loci using PCR-based DNA tests. Polymorphism was assessed at loci which 1) are known to cause phenotypic variation, 2) are potentially involved in trait differences or 3) are putative candidate genes. The traits considered are disease resistance, growth, coat colour, meat quality and prolificacy. Even though the populations are small and the markers are limited to specific genes, we found significant differences in five of the ten characterised loci. In some cases the observed allele frequencies were interesting in relation to gene function and the phenotype of the breed. These breeds are part of a conservation programme in Romania and marker information may be useful in preserving a representative gene pool in the populations. The use of polymorphisms in type 1 (gene) markers may be a useful complement to analysis based on anonymous markers.     

Patterns of genetic and phenotypic variation in Iris haynei and I. atrofusca (Iris sect. Oncocyclus = the royal irises) along an ecogeographical gradient in Israel and the West Bank

Iris haynei and I. atrofusca are two closely related narrow endemics distributed vicariously along an ecogeographical north-south gradient in Israel and the West Bank. To obtain baseline information of the taxonomic status, conservation and population history of these taxa, we investigated patterns of phenotypic variation and the partitioning of genetic variation within and among populations using dominant random amplified polymorphic DNA (RAPD) markers. Multivariate (principal components analysis) and taxonomic distance analyses based on morphometric traits from eight populations revealed no unambiguous separation into two distinct groups. Results of genetic analyses for nine populations differed only slightly when either allele- or marker-based approaches were employed. Mean within-population diversity was high (0.258 for Nei's expected heterozygosity), but there was no significant relationship between genetic diversity and either population size or latitude. Although the range-wide estimate of GST ( approximately 0.20) revealed relatively high differentiation among populations this value was inflated because of a small, but significant, component of molecular variance among regions viz. taxa ( approximately 5%). Limited long-distance dispersal capabilities in conjunction with a linearized habitat distribution are proposed to contribute to the approximate isolation by distance pattern observed. It also appears that extant populations are currently deviating from equilibrium conditions because of primary divergence of a formerly more widespread ancestral population. Given the absence of deep genetic and phenotypic subdivision among northern (I. haynei) vs. central/southern (I. atrofusca) populations, we argue for a revision of their species status. Nonetheless, we recommend conservation attention to these geographically differentiated segments as separate management units, which can be seen as an instructive example of incipient species formation.     

Relationships between the endangered Pustertaler-Sprinzen and three related European cattle breeds as analysed with 20 microsatellite loci

We estimated the genetic relationships between the endangered German Pustertaler-Sprinzen cattle breed and the Pinzgauer, Vosges and Simmental breeds--decided upon after consultation of the available historical literature. Within-breed diversity of the four breeds was also assessed. Twenty microsatellite markers were amplified in 27-50 unrelated individuals from populations of each breed. Within-breed variation was estimated from average heterozygosity values and mean number of alleles. Breed relationships were evaluated by genetic distance and a neighbour-joining tree was calculated from these estimates. Bootstrap resampling of loci tested the robustness of the tree topology obtained. A tree was also constructed from distance matrices using individual animals as operational taxonomic units. From both the average heterozygosity values and mean number of alleles calculated, the Pustertaler breed appears to be no more genetically impoverished than the other breeds analysed. The breed tree showed an 85% support for the Pustertaler-Pinzgauer grouping, and this result is echoed in the genetic distance values and allele-sharing individual tree.     

Development and validation of microsatellite markers for the tiny dragonfly, <Emphasis Type="Italic">Nannophya pygmaea</Emphasis> (Odonata: Libellulidae), which is endangered in South Korea

The tiny dragonfly, Nannophya pygmaea (Odonata: Libellulidae), is listed as a second-degree endangered wild animal in South Korea. The application of molecular markers to assess genetic diversity and population relationships can provide information necessary to establish an effective conservation strategy. In this study, we developed 12 microsatellite markers specific to N. pygmaea using the NextSeq 500 platform. Forty individuals of N. pygmaea collected from three currently known localities in South Korea were genotyped to validate these markers and to preliminarily assess population genetic characteristics. The observed number of alleles, observed heterozygosity, and expected heterozygosity at a locus ranged from 2 to 9, 0.421–1.0, and 0.508–0.766 in a population with the largest sample size (20 individuals), respectively, thereby validating the suitability of the markers for population analysis. Five of 12 loci showed significant deviation from the Hardy–Weinberg equilibrium in the population. Our preliminary data indicate an absence of inbreeding in all populations and an absence of obvious genetic difference. The microsatellite markers developed in this study will be useful for studying the population genetics of N. pygmaea collected from other regions worldwide, including additional sites in South Korea.     

Geographical variation of genetic and phenotypic traits in the Mexican sailfin mollies, Poecilia velifera and P. petenensis

Comparing the patterns of population divergence using both neutral genetic and phenotypic traits provides an opportunity to examine the relative importance of evolutionary mechanisms in shaping population differences. We used microsatellite markers to examine population genetic structure in the Mexican sailfin mollies Poecilia velifera and P. petenensis. We compared patterns of genetic structure and divergence to that in two types of phenotypic traits: morphological characters and mating behaviours. Populations within each species were genetically distinct, and conformed to a model of isolation by distance, with populations within different geographical regions being more genetically similar to one another than were populations from different regions. Bayesian clustering and barrier analyses provided additional support for population separation, especially between geographical regions. In contrast, none of the phenotypic traits showed any type of geographical pattern, and population divergence in these traits was uncorrelated with that found in neutral markers. There was also a weaker pattern of regional differences among geographical regions compared to neutral genetic divergence. These results suggest that while divergence in neutral traits is likely a product of population history and genetic drift, phenotypic divergence is governed by different mechanisms, such as natural and sexual selection, and arises at spatial scales independent from those of neutral markers.     

长蓝猪QTL定位资源群建立及其遗传分析

利用高度选育的8头瘦肉型长白猪为父本,16头中国优良地方品种蓝塘猪为母本,采用远交系杂交的F2代设计方法,建立包括8头F1公猪、40头F1母猪和232头F2个体的猪资源群体。资源群体的遗传分析表明,所测定的32个性状都有一定的变异,主要经济性状的变异系数都在10％以上。进一步的方差组分估计表明,主要经济性状的加性遗传方差较大。在所测定6号染色体连锁群的22个微卫星DNA标记中,12个表现有多态性,全群平均杂合度为0.53,多态信息含量为0.46,能较好地为QTL检测提供信息。因此,F2代有较好的分离状态,有可能利用该资源群体检测到较多的QTL。     

Genetic diversity in the Churra tensina and Churra lebrijana endangered Spanish sheep breeds and relationship with other Churra group breeds and Spanish mouflon

The aim of the present study was to estimate the genetic intra-breed variability of Churra tensina and Churra lebrijana endangered breeds and to establish genetic relationships with Churra, Latxa and Merino breeds, as well as Spanish mouflon, by using 28 microsatellite markers, to provide useful information for their conservation. Allele frequencies and heterozygosity revealed high genetic variation in the two endangered breeds despite their small population size. Estimates of inbreeding coefficient (F_IS) were significant for all breeds studied, except for Churra lebrijana breed. The highest inbreeding coefficient (F_IS = 0.143) was found in the Spanish mouflon. Genetic differentiation tests (F_ST = 0.121) and assignment of individuals to populations indicated the existence of defined breed populations, and low genetic flow between these breeds. The highest pairwise Reynolds distance (D_R) values were observed between Mouflon and the domestic sheep breeds. Considering only domestic sheep breeds, the Churra lebrijana breed showed the highest pairwise D_R values_. The lowest values were found between Latxa and the other domestic sheep, except for Churra lebrijana. Results of pairwise D_R values, as well as phylogenetic tree and bottleneck analysis showed an important genetic isolation of the Churra lebrijana breed from the other Churra types, and genetic signatures of a demographic bottleneck. Finally, structure analysis of populations detected a population subdivision in the Latxa sheep breed. In conclusion, this study presents valuable insight into the existing genetic variability of two Spanish endangered breeds, as well as the first study in Spanish mouflon based on microsatellite analysis. The high degree of variability demonstrated in Churra tensina and Churra lebrijana implies that these populations are rich reservoirs of genetic diversity.     

Indigenous domestic breeds as reservoirs of genetic diversity: the Argentinean Creole cattle

Contrary to highly selected commercial breeds, indigenous domestic breeds are composed of semi-wild or feral populations subjected to reduced levels of artificial selection. As a consequence, many of these breeds have become locally adapted to a wide range of environments, showing high levels of phenotypic variability and increased fitness under natural conditions. Genetic analyses of three loci associated with milk production (alpha(S1)-casein, kappa-casein and prolactin) and the locus BoLA-DRB3 of the major histocompatibility complex indicated that the Argentinean Creole cattle (ACC), an indigenous breed from South America, maintains high levels of genetic diversity and population structure. In contrast to the commercial Holstein breed, the ACC showed considerable variation in heterozygosity (H(e)) and allelic diversity (A) across populations. As expected, bi-allelic markers showed extensive variation in He whereas the highly polymorphic BoLA-DRB3 showed substantial variation in A, with individual populations having 39-74% of the total number of alleles characterized for the breed. An analysis of molecular variance (AMOVA) of nine populations throughout the distribution range of the ACC revealed that 91.9-94.7% of the total observed variance was explained by differences within populations whereas 5.3-8.1% was the result of differences among populations. In addition, the ACC breed consistently showed higher levels of genetic differentiation among populations than Holstein. Results from this study emphasize the importance of population genetic structure within domestic breeds as an essential component of genetic diversity and suggest that indigenous breeds may be considered important reservoirs of genetic diversity for commercial domestic species.     

How to use molecular marker data to measure evolutionary parameters in wild populations

Estimating the genetic basis of phenotypic traits and the selection pressures acting on them are central to our understanding of the evolution and conservation of wild populations. However, obtaining such evolutionary-related parameters is not an easy task as it requires accurate information on both relatedness among individuals and their breeding success. Polymorphic molecular markers are very useful in estimating relatedness between individuals and parentage analyses are now extensively used in most taxa. The next step in the application of molecular data to wild populations is to use them to derive estimates of evolutionary-related parameters for quantitative traits, such as quantitative genetic parameters (e.g. heritability, genetic correlations) and measures of selection (e.g. selection gradients). Despite their great appeal and potential, the optimal use of molecular tools is still debated and it remains unclear how they should best be used to obtain reliable estimates of evolutionary parameters in the wild. Here, we review the methods available for estimating quantitative genetic and selection parameters and discuss their merits and shortcomings, to provide a tool that summarizes the potential uses of molecular data to obtain such parameters in wild populations.     

Conservation genetics of freshwater fish

Genetic markers have helped to resolve many difficult taxonomic problems and map patterns of diversity within and among remnant populations of threatened and endangered species. Knowledge of historical patterns of gene flow can help to manage dispersal among anthropogenically fragmented populations. Genetic considerations are used in the design of captive breeding programmes that avoid inbreeding depression and artificial selection that may impact on Darwinian fitness. Case studies from endangered populations of topminnows from North American deserts are used to illustrate a variety of methods used in conservation genetic studies. Several merits of studying putatively neutral, molecular markers v. adaptive phenotypic traits are discussed.     

闽江水系绒螯蟹的表型性状差异与分子遗传差异的比较与分析

表型性状差异(differentiation in phenotypictraits,PST)和分子遗传差异(differentiation at neutral molecular markers,FST)是近期进化生物学的研究热点之一。闽江水系是我国中华绒螯蟹与合浦绒螯蟹的主要混杂地域,是研究绒螯蟹遗传与进化的理想地之一。为探讨闽江水系绒螯蟹的PST和FST,以2009和2010年度闽江水系的133个绒螯蟹样本为材料,进行了14个表型数量性状差异分析和6个微卫星标记的遗传差异分析。结果发现:除3个表型性状不存在显著差异外,其他表型性状在不同年份间均存在极显著差异(P<0.01);2009年绒螯蟹的平均期望杂合度极显著高于2010年绒螯蟹(P=0.008),而平均等位基因丰富度、观测杂合度和近交系数均不存在显著差异(P=0.136～0.675);年份间的平均FST为0.1429;通过对PST与FST的比较发现,除第二步足掌节长度(F2)性状外,其他表型数量性状的PST值均高于FST值,表明这些性状均受到了较明显的选择压力。该文研究结果为绒螯蟹的分子进化研究积累了资料,也为其他水产生物的PST和FST比较研究提供了参考。     

Population genetic characterization of the endangered dung beetle Copris tripartitus (Coleoptera: Scarabaeidae) using novel microsatellite markers

The dung beetle Copris tripartitus (Coleoptera: Scarabaeidae) has long been considered an endangered insect in South Korea; the detection of recent population increases leaves its endangered status uncertain. Population genetic analysis subsequent to development of molecular markers is essential for establishing proper conservation strategies. In this study, we developed ten microsatellite markers specifically for C. tripartitus. Sixty-eight individuals of C. tripartitus collected from six South Korean localities were genotyped to validate these markers and preliminarily assess population genetic characteristics. Per-locus observed number of alleles, observed heterozygosity (H_O), and expected heterozygosity (H_E) ranged from 5–12, 0.499–0.958, and 0.54–0.743, respectively. All populations showed higher H_O than H_E, negative values of inbreeding coefficient, and, overall, no sign of recent population bottlenecks (excluding one population, Seosan). This suggests that C. tripartitus did not suffer from genetic drift and inbreeding, which are typically severe in small, isolated populations. Nevertheless, detection of only one of the two gene pools in some populations and resultant genetic subdivision into two population groups may suggest that the population size is not enough to cover both gene pools. Thus, a more extended period of protection may be required to ensure higher genetic diversity of widespread populations and achieve the long-term conservation goal.     

Genome-wide analyses of the Jeju,Thoroughbred, and Jeju crossbred horse populations using the high density SNP array

The Jeju horse is an indigenous Korean horse breed that is currently registered with the Food and Agriculture Organization of the United Nations. However, there is severe lack of genomic studies on Jeju horse. This study was conducted to investigate genetic characteristics of horses including Jeju horse, Thoroughbred and Jeju crossbred (Jeju?×?Thoroughbred) populations. We compared the genomes of three horse populations using the Equine SNP70 Beadchip array. Short-range Linkage disequilibrium was the highest in Thoroughbred, whereas r² values were lowest in Jeju horse. Expected heterozygosity was the highest in Jeju crossbred (0.351), followed by the Thoroughbred (0.337) and Jeju horse (0.311). The level of inbreeding was slightly higher in Thoroughbred (??0.009) than in Jeju crossbred (??0.035) and Jeju horse (??0.038). F_ST value was the highest between Jeju horse and Thoroughbred (0.113), whereas Jeju crossbred and Thoroughbred showed the lowest value (0.031). The genetic relationship was further assessed by principal component analysis, suggesting that Jeju crossbred is more genetically similar to Thoroughbred than Jeju horse population. Additionally, we detected potential selection signatures, for example, in loci located on LCORL/NCAPG and PROP1 genes that are known to influence body. Genome-wide analyses of the three horse populations showed that all the breeds had somewhat a low level of inbreeding within each population. In the population structure analysis, we found that Jeju crossbred was genetically closer to Thoroughbred than Jeju horse. Furthermore, we identified several signatures of selection which might be associated with traits of interest. To our current knowledge, this study is the first genomic research, analyzing genetic relationships of Jeju horse, Thoroughbred and Jeju crossbred.