The inference of population genetic structures is essential in many research areas in population genetics, conservation biology and evolutionary biology. Recently, unsupervised Bayesian clustering algorithms have been developed to detect a hidden population structure from genotypic data, assuming among others that individuals taken from the population are unrelated. Under this assumption, markers in a sample taken from a subpopulation can be considered to be in Hardy-Weinberg and linkage equilibrium. However, close relatives might be sampled from the same subpopulation, and consequently, might cause Hardy-Weinberg and linkage disequilibrium and thus bias a population genetic structure analysis. In this study, we used simulated and real data to investigate the impact of close relatives in a sample on Bayesian population structure analysis. We also showed that, when close relatives were identified by a pedigree reconstruction approach and removed, the accuracy of a population genetic structure analysis can be greatly improved. The results indicate that unsupervised Bayesian clustering algorithms cannot be used blindly to detect genetic structure in a sample with closely related individuals. Rather, when closely related individuals are suspected to be frequent in a sample, these individuals should be first identified and removed before conducting a population structure analysis.

Hybrid zones of ecologically divergent populations are ideal systems to study the interaction between natural selection and gene flow during the initial stages of speciation. Here, we perform an amplified fragment length polymorphism (AFLP) genome scan in parallel hybrid zones between divergent ecotypes of the marine snail Littorina saxatilis, which is considered a model case for the study of ecological speciation. Ridged-Banded (RB) and Smooth-Unbanded (SU) ecotypes are adapted to different shore levels and microhabitats, although they present a sympatric distribution at the mid-shore where they meet and mate (partially assortatively). We used shell morphology, outlier and nonoutlier AFLP loci from RB, SU and hybrid specimens captured in sympatry to determine the level of phenotypic and genetic introgression. We found different levels of introgression at parallel hybrid zones and nonoutlier loci showed more gene flow with greater phenotypic introgression. These results were independent from the phylogeography of the studied populations, but not from the local ecological conditions. Genetic variation at outlier loci was highly correlated with phenotypic variation. In addition, we used the relationship between genetic and phenotypic variation to estimate the heritability of morphological traits and to identify potential Quantitative Trait Loci to be confirmed in future crosses. These results suggest that ecology (exogenous selection) plays an important role in this hybrid zone. Thus, ecologically based divergent natural selection is responsible, simultaneously, for both ecotype divergence and hybridization. On the other hand, genetic introgression occurs only at neutral loci (nonoutliers). In the future, genome-wide studies and controlled crosses would give more valuable information about this process of speciation in the face of gene flow.

Population subdivision must be explicitly considered in the management of conservation programmes, as most populations of wild species at risk of extinction and those kept in captivity are spatially structured. The partition of gene and allelic diversity in within-and between-subpopulation components allows for the integral management of populations. We summarise the main aspects of this partition and some of its applications in terms of priorisation of populations for conservation and establishment of synthetic populations. The procedures for the maintenance of diversity in subdivided populations making use of molecular markers and its implementation by the software METAPOP are illustrated with empirical data.

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.

Before availability of dense SNP data, genetic diversity was characterized and managed with pedigreebased information. Besides this classical approach, 2 methodologies have been proposed in recent years to characterize and manage diversity from dense SNP data: the SNP-by-SNP approach and the alternative based on runs of homozygosity (ROH). The establishment of criteria to identify ROH is a current constraint in the literature dealing with ROH. The objective of this study was, using a medium-density SNP chip, to quantify by 3 methods (pedigree, SNP-by-SNP, and ROH) the genetic diversity on 5 selected French dairy sheep subpopulations and breeds and to assess the effect of the definition of ROH on these estimates. The data set available included individuals from the breeds Basco-Béarnaise, Manech Tête Noire, Manech Tête Rousse, and 2 subpopulations of Lacaune: Lacaune Confederation and Lacaune Ovitest. Animals were genotyped with the Illumina OvineSNP50 BeadChip (Illumina Inc., San Diego, CA). After filtering, the genomic data included 38,287 autosomal SNP and 8,700 individuals, which comprised 72,803 animals in the pedigree. The results indicated that no significant differences were observed in effective population size estimates obtained from pedigree or genomic (SNPby-SNP or ROH) information. In general, estimates of effective population size were above 200 in Lacaune Confederation and Lacaune Ovitest subpopulations and below 200 in Basco-Béarnaise, Manech Tête Noire, and Manech Tête Rousse breeds. The minimum length that constituted a ROH, the minimum number of SNP that constituted a ROH, as well as the minimum density and the maximum distance allowed between 2 homozygous SNP are ROH-defining factors with important implications in the estimation of the rate of inbreeding. The ROH-based rates of inbreeding in concordance with those obtained from pedigree information require a specific set of values. This particular set of values is different from that identified to obtain ROH-based rates of inbreeding similar to those obtained on a SNPby-SNP basis. Factors to define ROH do not change the results much unless extreme values are considered, although further research on ROH-based inbreeding is still required.

Copy number variants (CNV) are structural variants consisting of duplications or deletions of genomic fragments longer than 1 kb that present variability in the population and are heritable. The objective of this study was to identify CNV regions (CNVR) associated with 7 economically important traits (production, functional, and type traits) in Holstein cattle: fat yield, protein yield, somatic cell count, days open, stature, foot angle, and udder depth. Copy number variants were detected by using deep-sequencing data from 10 sequenced bulls and the Bovine SNP chip array hybridization signals. To reduce the number of false-positive calls, only CNV identified by both sequencing and Bovine SNP chip assays were kept in the final data set. This resulted in 823 CNVR. After filtering by minor allele frequency >0.01, a total of 90 CNVR appeared segregating in the bulls that had phenotypic data. Linear and quadratic CNVR effects were estimated using Bayesian approaches. A total of 15 CNVR were associated with the traits included in the analysis. One CNVR was associated with fat and protein yield, another 1 with fat yield, 3 with stature, 1 with foot angle, 7 with udder depth, and only 1 with days open. Among the genes located within these regions, highlighted were the MTHFSD gene that belongs to the folate metabolism genes, which play critical roles in regulating milk protein synthesis; the SNRPE gene that is related to several morphological pathologies; and the NF1 gene, which is associated with potential effects on fertility traits. The results obtained in the current study revealed that these CNVR segregate in the Holstein population, and therefore some potential exists to increase the frequencies of the favorable alleles in the population after independent validation of results in this study. However, genetic variance explained by the variants reported in this study was small.

Equalization of parental contributions is one of the most simple and widely recognized methods to maintain genetic diversity in conservation programs, as it halves the rate of increase in inbreeding and genetic drift. It has, however, the negative side effect of implying a reduced intensity of natural selection so that deleterious genes are less efficiently removed from the population with possible negative consequences on the reproductive capacity of the individuals. Theoretical results suggest that the lower fitness resulting from equalization of family sizes relative to that for free contribution schemes is expected to be substantial only for relatively large population sizes and after many generations. We present a long-term experiment with Drosophila melanogaster, comparing the fitness performance of lines maintained with equalization of contributions (EC) and others maintained with no management (NM), allowing for free matings and contributions from parents. Two (five) replicates of size N ¼ 100 (20) individuals of each type of line were maintained for 38 generations. As expected, EC lines retained higher gene diversity and allelic richness for four microsatellite markers and a higher heritability for sternopleural bristle number. Measures of life-history traits, such as egg-to-adult viability, mating success, and global fitness declined with generations, but no significant differences were observed between EC and NM lines. Our results, therefore, provide no evidence to suggest that equalization of family sizes entails a disadvantage on the reproductive capacity of conserved populations in comparison with no management procedures, even after long periods of captivity.

One of the main objectives of genetic breeding programs in turbot industry is to reduce disease-related mortality. In the present study, a genome scan to detect quantitative trait loci (QTL) affecting resistance and survival to viral haemorrhagic septicaemia (VHS) was carried out. Three full-sib families with approximately 90 individuals each were genotyped and evaluated by linear regression and maximum likelihood approaches. In addition, a comparison between QTL detected for resistance and survival time to other important bacterial and parasite diseases affecting turbot (furunculosis and scuticociliatosis) was also carried out. Finally, the relationship between QTL affecting resistance/survival time to the virus and growth-related QTL was also evaluated. Several genomic regions controlling resistance and survival time to VHS were detected. Also significant associations between the evaluated traits and genotypes at particular markers were identified, explaining up to 14 % of the phenotypic variance. Several genomic regions controlling general and specific resistance to different diseases in turbot were detected. A preliminary gene mining approach identified candidate genes related to general or specific immunity. This information will be valuable to develop marker-assisted selection programs and to discover candidate genes related to disease resistance to improve turbot production.