Partitioning of the genetic trends of French dairy sheep in Mendelian samplings and long-term contributions.

The genetic trend of milk yield for 4 French dairy sheep breeds (Lacaune, Basco-Béarnaise, Manech Tête Noire, and Manech Tête Rousse) was partitioned in Mendelian sampling trends by categories of animals defined by sex and by selection pathways. Five categories were defined, as follows: (1) artificial insemination (AI) males (after progeny testing), (2) males discarded after progeny testing, (3) natural mating males, (4) dams of males, and (5) dams of females. Dams of males and AI males were the most important sources of genetic progress, as observed in the decomposition in Mendelian sampling trends. The yearly contributions were more erratic for AI males than for dams of males, as AI males are averaged across a smaller number of individuals. Natural mating males and discarded males did not contribute to the trend in terms of Mendelian sampling, as their estimated Mendelian sampling term is either null (natural mating males) or negative (discarded males). Overall, in terms of Mendelian sampling, females contributed more than males to the total genetic gain, and we interpret that this is because females constitute a larger pool of genetic diversity. In addition, we computed long-term contributions from each individual to the following pseudo-generations (one pseudo-generation spanning 4 years). With this information, we studied the selection decisions (selected or not selected) for females, and the contributions to the following generations. Mendelian sampling was more important than parent average to determine the selection of individuals and their long-term contributions. Long-term contributions were greater for AI males (with larger progeny sizes than females) and in Basco-Béarnaise than in Lacaune (with the latter being a larger population).

Genomic and pedigree estimation of inbreeding depression for semen traits in the Basco-Béarnaise dairy sheep breed.

Inbreeding depression is associated with a decrease in performance and fitness of the animals. The goal of this study was to evaluate pedigree-based and genomic methods to estimate the level of inbreeding and inbreeding depression for 3 semen traits (volume, concentration, and motility score) in the Basco-Béarnaise sheep breed. Data comprised 16,196 (or 15,071) phenotypic records from 620 rams (of which 533 rams had genotypes of 36,464 SNPs). The pedigree included 8,266 animals, composed of the 620 rams and their ancestors. The number of equivalent complete generations for the 620 rams was 7.04. Inbreeding coefficients were estimated using genomic and pedigree-based information. Genomic inbreeding coefficients were estimated from individual SNP and using segments of homozygous SNP (runs of homozygosity, ROH). Short ROH are of old origin, whereas long ROH are due to recent inbreeding. Considering that the equivalent number of generations in Basco-Béarnaise was 6, inbreeding coefficients for ROH with a length >4 Mb refer to all (recent + old) inbreeding, those with a length >17 Mb correspond to recent inbreeding, and the difference between them indicates old inbreeding. Pedigree-based inbreeding coefficients were also estimated classically, or accounting for nonzero relationships for unknown parents, or including metafounder relationships (estimated using markers) to account for missing pedigree information. Finally, inbreeding coefficients combining genotyped and nongenotyped animal information were computed from matrix H of the single-step approach, also including metafounders. Inbreeding depression was estimated differently depending on the approach used to compute inbreeding coefficients. These 8 estimators of inbreeding coefficients were included as covariates in different animal models. No inbreeding depression was detected for sperm volume or sperm concentration. Inbreeding depression was significant for the motility of spermatozoa. The effect of old and recent inbreeding on motility was null and negative, respectively, demonstrating the existence of purging by selection of deleterious recessive alleles affecting motility. A 10% increase in inbreeding would result in a reduction in mean motility ranging between 0.09 and 0.22 points in the score (from 0 to 5). Motility is unfavorably affected by increasing recent inbreeding but the impact is very small. Runs of homozygosity and metafounders allow us to accurately estimate inbreeding depression and detect recent inbreeding.

Islands of runs of homozygosity indicate selection signatures in Ovis aries 6 (OAR6) of French dairy sheep.

Runs of homozygosity (ROH) are contiguous homozygous segments of the genome where the haplotypes inherited from each parent are identical. The occurrence of ROH is not randomly distributed across the genome, and ROH islands across many animals may be the result of selective pressure. The objective of this study was to demonstrate that the presence of ROH islands may be indicative of selection signatures in French dairy sheep breeds and subpopulations. The data set available included animals (artificial insemination males) from various breeds and subpopulations: Basco-Béarnaise breed (321 individuals), Manech Tête Noire breed (329 individuals), Manech Tête Rousse breed (1,906 individuals), Lacaune Confederation subpopulation (3,030 individuals), and Lacaune Ovitest subpopulation (3,114 individuals). Animals were genotyped with the Illumina OvineSNP50 BeadChip. After applying filtering criteria, the genomic data included 38,287 autosomal SNP distributed across 26 chromosomes and 8,700 individuals. One island of ROH was detected on OAR6 in the same genomic position across animals (between 30 and 40 Mb). Global Wright's differentiation coefficients for 2 SNP within this ROH island were high (0.67-0.68). The linkage disequilibrium between both SNP was also elevated (0.98). The divergence in allele frequencies in those SNP grouped Basco-Béarnaise, Manech Tête Noire, and Manech Tête Rousse breeds in one cluster and Lacaune Confederation and Lacaune Ovitest subpopulations in another cluster. The closest candidate genes are NCAPG and LCORL, which have been reported to be under positive selection and suggested to control weight and height in sheep. The preliminary identification of ROH suggests the presence of selection. However, for the identification of potential candidate genes, ROH detection should be combined with other approaches to improve mapping accuracy.

Inbreeding, effective population size, and coancestry in the Latxa dairy sheep breed.

Traditionally, breeding programs have estimated and managed inbreeding based on pedigree information. The availability of genomic marker panels has made possible new alternatives to achieve more precise estimates, for example in case of missing pedigree. The objective of the present study was to assess and compare, different estimation methods (pedigree-based methodologies, single SNP-based approach (homozygosity) and runs of homozygosity-based method) to analyze the evolution of genetic diversity measured as inbreeding or as coancestry of 3 selected populations of Latxa dairy sheep (Latxa Cara Rubia and Latxa Cara Negra from Euskadi and Navarre). Genomic data came from 972 artificial insemination rams genotyped with the Illumina OvineSNP50 BeadChip (Illumina Inc., San Diego, CA) whose genealogy consisted of 4,484 animals. Inbreeding estimates based on molecular data were more similar between them than compared with those based on pedigree information. However, the SNP-based approach estimations of effective population size differed more, reflecting the sensitivity of effective population size to small changes in the evolution of inbreeding. The 2 Latxa Cara Negra populations showed increases of inbreeding rates with time and effective population sizes between 64 and 103 animals, depending on breed and methodology used. The Latxa Cara Rubia population did not show an increase in inbreeding rate, mainly due to semen importation from the related French population of Manech Tête Rousse. The effective size estimates based on coancestry increase show a higher variability and they are more sensitive to the source of information and the data structure considered. Realized effective population size based on individual increase in inbreeding were in agreement with the previous estimates. Coancestry evolution analysis based on DNA information showed an increase on coancestry during the last 10 yr in all breeds, as a consequence of the selection process. Moreover, the increase on coancestry between Latxa Cara Rubia and Manech Tête Rousse was more noticeable between than within each of those breeds.

Inbreeding and effective population size in French dairy sheep: Comparison between genomic and pedigree estimates.

Before availability of dense SNP data, genetic diversity was characterized and managed with pedigree-based 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 (SNP-by-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 SNP-by-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.

The role of local ecology during hybridization at the initial stages of ecological speciation in a marine snail.

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.

Uncovering QTL for resistance and survival time to Philasterides dicentrarchi in turbot (Scophthalmus maximus).

Disease resistance-related traits have received increasing importance in aquaculture breeding programs worldwide. Currently, genomic information offers new possibilities in breeding to address the improvement of this kind of traits. The turbot is one of the most promising European aquaculture species, and Philasterides dicentrarchi is a scuticociliate parasite causing fatal disease in farmed turbot. An appealing approach to fight against disease is to achieve a more robust broodstock, which could prevent or diminish the devastating effects of scuticociliatosis on farmed individuals. In the present study, a genome scan for quantitative trait loci (QTL) affecting resistance and survival time to P. dicentrarchi in four turbot families was carried out. The objectives were to identify QTL using different statistical approaches [linear regression (LR) and maximum likelihood (ML)] and to locate significantly associated markers for their application in genetic breeding strategies. Several genomic regions controlling resistance and survival time to P. dicentrarchi were detected. When analyzing each family separately, significant QTL for resistance were identified by the LR method in two linkage groups (LG1 and LG9) and for survival time in LG1, while the ML methodology identified QTL for resistance in LG9 and LG23 and for survival time in LG6 and LG23. The analysis of the total data set identified an additional significant QTL for resistance and survival time in LG3 with the LR method. Significant association between disease resistance-related traits and genotypes was detected for several markers, a single one explaining up to 22% of the phenotypic variance. Obtained results will be essential to identify candidate genes for resistance and to apply them in marker-assisted selection programs to improve turbot production.

Relaxation of selection with equalization of parental contributions in conservation programs: an experimental test with Drosophila melanogaster.

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.

Mutation-selection balance accounting for genetic variation for viability in Drosophila melanogaster as deduced from an inbreeding and artificial selection experiment.

We carried out an experiment of inbreeding and upward artificial selection for egg-to-adult viability in a recently captured population of Drosophila melanogaster, as well as computer simulations of the experimental design, in order to obtain information on the nature of genetic variation for this important fitness component. The inbreeding depression was linear with a rate of 0.70 +/- 0.11% of the initial mean per 1% increase in inbreeding coefficient, and the realized heritability was 0.06 +/- 0.07. We compared the empirical observations of inbreeding depression and selection response with computer simulations assuming a balance between the occurrence of partially recessive deleterious mutations and their elimination by selection. Our results suggest that a model assuming mutation-selection balance with realistic mutational parameters can explain the genetic variation for viability in the natural population studied. Several mutational models are incompatible with some observations and can be discarded. Mutational models assuming a low rate of mutations of large average effect and highly recessive gene action, and others assuming a high rate of mutations of small average effect and close to additive gene action, are compatible with all the observations.

Lack of nonadditive genetic effects on early fecundity in Drosophila melanogaster.

Fecundity is usually considered as a trait closely connected to fitness and is expected to exhibit substantial nonadditive genetic variation and inbreeding depression. However, two independent experiments, using populations of different geographical origin, indicate that early fecundity in Drosophila melanogaster behaves as a typical additive trait of low heritability. The first experiment involved artificial selection in inbred and non-inbred lines, all of them started from a common base population previously maintained in the laboratory for about 35 generations. The realized heritability estimate was 0.151 +/- 0.075 and the inbreeding depression was very small and nonsignificant (0.09 +/- 0.09% of the non-inbred mean per 1% increase in inbreeding coefficient). With inbreeding, the observed decrease in the within-line additive genetic variance and the corresponding increase of the between-line variance were very close to their expected values for pure additive gene action. This result is at odds with previous studies showing inbreeding depression and, therefore, directional dominance for the same trait and species. All experiments, however, used laboratory populations, and it is possible that the original genetic architecture of the trait in nature was subsequently altered by the joint action of random drift and adaptation to captivity. Thus, we carried out a second experiment, involving inbreeding without artificial selection in a population recently collected from the wild. In this case we obtained, again, a maximum-likelihood heritability estimate of 0.210 +/- 0.027 and very little nonsignificant inbreeding depression (0.06 +/- 0.12%). The results suggest that, for fitness-component traits, low levels of additive genetic variance are not necessarily associated with large inbreeding depression or high levels of nonadditive genetic variance.