Comparative evaluation of genomic inbreeding parameters in seven commercial and autochthonous pig breeds

Single nucleotide polymorphism (SNP) genotyping tools, which can analyse thousands of SNPs covering the whole genome, have opened new opportunities to estimate the inbreeding level of animals directly using genome information. One of the most commonly used genomic inbreeding measures considers the proportion of the autosomal genome covered by runs of homozygosity (ROH), which are defined as continuous and uninterrupted chromosome portions showing homozygosity at all loci. In this study, we analysed the distribution of ROH in three commercial pig breeds (Italian Large White, n = 1968; Italian Duroc, n = 573; and Italian Landrace, n = 46) and four autochthonous breeds (Apulo-Calabrese, n = 90; Casertana, n = 90; Cinta Senese, n = 38; and Nero Siciliano, n = 48) raised in Italy, using SNP data generated from Illumina PorcineSNP60 BeadChip. We calculated ROH-based inbreeding coefficients (F_ROH) using ROH of different minimum length (1, 2, 4, 8, 16 Mbp) and compared them with several other genomic inbreeding coefficients (including the difference between observed and expected number of homozygous genotypes (F_HOM)) and correlated all these genomic-based measures with the pedigree inbreeding coefficient (F_PED) calculated for the pigs of some of these breeds. Autochthonous breeds had larger mean size of ROH than all three commercial breeds. F_HOM was highly correlated (0.671 to 0.985) with F_ROH measures in all breeds. Apulo-Calabrese and Casertana had the highest F_ROH values considering all ROH minimum lengths (ranging from 0.273 to 0.189 and from 0.226 to 0.152, moving from ROH of minimum size of 1 Mbp (F_ROH1) to 16 Mbp (F_ROH16)), whereas the lowest F_ROH values were for Nero Siciliano (from 0.072 to 0.051) and Italian Large White (from 0.117 to 0.042). F_ROH decreased as the minimum length of ROH increased for all breeds. Italian Duroc had the highest correlations between all F_ROH measures and F_PED (from 0.514 to 0.523) and between F_HOM and F_PED (0.485). Among all analysed breeds, Cinta Senese had the lowest correlation between F_ROH and F_PED. This might be due to the imperfect measure of F_PED, which, mainly in local breeds raised in extensive production systems, cannot consider a higher level of pedigree errors and a potential higher relatedness of the founder population. It appeared that ROH better captured inbreeding information in the analysed breeds and could complement pedigree-based inbreeding coefficients for the management of these genetic resources.     

An evaluation of inbreeding measures using a whole-genome sequenced cattle pedigree

The estimation of the inbreeding coefficient (F) is essential for the study of inbreeding depression (ID) or for the management of populations under conservation. Several methods have been proposed to estimate the realized F using genetic markers, but it remains unclear which one should be used. Here we used whole-genome sequence data for 245 individuals from a Holstein cattle pedigree to empirically evaluate which estimators best capture homozygosity at variants causing ID, such as rare deleterious alleles or loci presenting heterozygote advantage and segregating at intermediate frequency. Estimators relying on the correlation between uniting gametes (F_UNI) or on the genomic relationships (F_GRM) presented the highest correlations with these variants. However, homozygosity at rare alleles remained poorly captured. A second group of estimators relying on excess homozygosity (F_HOM), homozygous-by-descent segments (F_HBD), runs-of-homozygosity (F_ROH) or on the known genealogy (F_PED) was better at capturing whole-genome homozygosity, reflecting the consequences of inbreeding on all variants, and for young alleles with low to moderate frequencies (0.10 < . < 0.25). The results indicate that F_UNI and F_GRM might present a stronger association with ID. However, the situation might be different when recessive deleterious alleles reach higher frequencies, such as in populations with a small effective population size. For locus-specific inbreeding measures or at low marker density, the ranking of the methods can also change as F_HBD makes better use of the information from neighboring markers. Finally, we confirmed that genomic measures are in general superior to pedigree-based estimates. In particular, F_PED was uncorrelated with locus-specific homozygosity.Subject terms: Conservation genomics, Animal breeding, Inbreeding     

Analysis of runs of homozygosity and their relationship with inbreeding in five cattle breeds farmed in Italy

Increased inbreeding is an inevitable consequence of selection in livestock populations. The analysis of high‐density single nucleotide polymorphisms (SNPs) facilitates the identification of long and uninterrupted runs of homozygosity (ROH) that can be used to identify chromosomal regions that are identical by descent. In this work, the distribution of ROH of different lengths in five Italian cattle breeds is described. A total of 4095 bulls from five cattle breeds (2093 Italian Holstein, 749 Italian Brown, 364 Piedmontese, 410 Marchigiana and 479 Italian Simmental) were genotyped at 54K SNP loci. ROH were identified and used to estimate molecular inbreeding coefficients (F_ROH), which were compared with inbreeding coefficients estimated from pedigree information (F_PED) and using the genomic relationship matrix (F_GRM). The average number of ROH per animal ranged from 54 ± 7.2 in Piedmontese to 94.6 ± 11.6 in Italian Brown. The highest number of short ROH (related to ancient consanguinity) was found in Piedmontese, followed by Simmental. The Italian Brown and Holstein had a higher proportion of longer ROH distributed across the whole genome, revealing recent inbreeding. The F_PED were moderately correlated with F_ROH > 1 Mb (0.662, 0.700 and 0.669 in Italian Brown, Italian Holstein and Italian Simmental respectively) but poorly correlated with F_GRM (0.134, 0.128 and 0.448 for Italian Brown, Italian Holstein and Italian Simmental respectively). The inclusion of ROH > 8 Mb in the inbreeding calculation improved the correlation of F_ROH with F_PED and F_GRM. ROH are a direct measure of autozygosity at the DNA level and can overcome approximations and errors resulting from incomplete pedigree data. In populations with high linkage disequilibrium (LD) and recent inbreeding (e.g. Italian Holstein and Italian Brown), a medium‐density marker panel, such as the one used here, may provide a good estimate of inbreeding. However, in populations with low LD and ancient inbreeding, marker density would have to be increased to identify short ROH that are identical by descent more precisely.     

A genome-wide scan study identifies a single nucleotide substitution in ASIP associated with white versus non-white coat-colour variation in sheep (Ovis aries)

In sheep, coat colour (and pattern) is one of the important traits of great biological, economic and social importance. However, the genetics of sheep coat colour has not yet been fully clarified. We conducted a genome-wide association study of sheep coat colours by genotyping 47 303 single-nucleotide polymorphisms (SNPs) in the Finnsheep population in Finland. We identified 35 SNPs associated with all the coat colours studied, which cover genomic regions encompassing three known pigmentation genes (TYRP1, ASIP and MITF) in sheep. Eighteen of these associations were confirmed in further tests between white versus non-white individuals, but none of the 35 associations were significant in the analysis of only non-white colours. Across the tests, the s66432.1 in ASIP showed significant association (P=4.2 × 10⁻¹¹ for all the colours; P=2.3 × 10⁻¹¹ for white versus non-white colours) with the variation in coat colours and strong linkage disequilibrium with other significant variants surrounding the ASIP gene. The signals detected around the ASIP gene were explained by differences in white versus non-white alleles. Further, a genome scan for selection for white coat pigmentation identified a strong and striking selection signal spanning ASIP. Our study identified the main candidate gene for the coat colour variation between white and non-white as ASIP, an autosomal gene that has been directly implicated in the pathway regulating melanogenesis. Together with ASIP, the two other newly identified genes (TYRP1 and MITF) in the Finnsheep, bordering associated SNPs, represent a new resource for enriching sheep coat-colour genetics and breeding.     

Genetic diversity analyses reveal first insights into breed‐specific selection signatures within Swiss goat breeds

We used genotype data from the caprine 50k Illumina BeadChip for the assessment of genetic diversity within and between 10 local Swiss goat breeds. Three different cluster methods allowed the goat samples to be assigned to the respective breed groups, whilst the samples of Nera Verzasca and Tessin Grey goats could not be differentiated from each other. The results of the different genetic diversity measures show that Appenzell, Toggenburg, Valais and Booted goats should be prioritized in future conservation activities. Furthermore, we examined runs of homozygosity (ROH) and compared genomic inbreeding coefficients based on ROH (F_ROH) with pedigree‐based inbreeding coefficients (F_PED). The linear relationship between F_ROH and F_PED was confirmed for goats by including samples from the three main breeds (Saanen, Chamois and Toggenburg goats). F_ROH appears to be a suitable measure for describing levels of inbreeding in goat breeds with missing pedigree information. Finally, we derived selection signatures between the breeds. We report a total of 384 putative selection signals. The 25 most significant windows contained genes known for traits such as: coat color variation (MITF, KIT, ASIP), growth (IGF2, IGF2R, HRAS, FGFR3) and milk composition (PITX2). Several other putative genes involved in the formation of populations, which might have been selected for adaptation to the alpine environment, are highlighted. The results provide a contemporary background for the management of genetic diversity in local Swiss goat breeds.     

Runs of homozygosity: current knowledge and applications in livestock

This review presents a broader approach to the implementation and study of runs of homozygosity (ROH) in animal populations, focusing on identifying and characterizing ROH and their practical implications. ROH are continuous homozygous segments that are common in individuals and populations. The ability of these homozygous segments to give insight into a population's genetic events makes them a useful tool that can provide information about the demographic evolution of a population over time. Furthermore, ROH provide useful information about the genetic relatedness among individuals, helping to minimize the inbreeding rate and also helping to expose deleterious variants in the genome. The frequency, size and distribution of ROH in the genome are influenced by factors such as natural and artificial selection, recombination, linkage disequilibrium, population structure, mutation rate and inbreeding level. Calculating the inbreeding coefficient from molecular information from ROH (F_ROH) is more accurate for estimating autozygosity and for detecting both past and more recent inbreeding effects than are estimates from pedigree data (F_PED). The better results of F_ROH suggest that F_ROH can be used to infer information about the history and inbreeding levels of a population in the absence of genealogical information. The selection of superior animals has produced large phenotypic changes and has reshaped the ROH patterns in various regions of the genome. Additionally, selection increases homozygosity around the target locus, and deleterious variants are seen to occur more frequently in ROH regions. Studies involving ROH are increasingly common and provide valuable information about how the genome's architecture can disclose a population's genetic background. By revealing the molecular changes in populations over time, genome‐wide information is crucial to understanding antecedent genome architecture and, therefore, to maintaining diversity and fitness in endangered livestock breeds.     

Individual increase in inbreeding allows estimating effective sizes from pedigrees

We present here a simple approach to obtain reliable estimates of the effective population size in real world populations via the computation of the increase in inbreeding for each individual (delta F_i) in a given population. The values of delta F_i are computed as t-root of 1 - (1 - F_i) where F_i is the inbreeding coefficient and t is the equivalent complete generations for each individual. The values of delta F computed for a pre-defined reference subset can be averaged and used to estimate effective size. A standard error of this estimate of N_e can be further computed from the standard deviation of the individual increase in inbreeding. The methodology is demonstrated by applying it to several simulated examples and to a real pedigree in which other methodologies fail when considering reference subpopulations. The main characteristics of the approach and its possible use are discussed both for predictive purposes and for analyzing genealogies.     

Genetic consequences of fragmentation in “arbor vitae,” eastern white cedar (Thuja occidentalis L.), toward the northern limit of its distribution range

We tested the hypothesis that marginal fragmented populations of eastern white cedar (EWC) are genetically isolated due to reduced pollen and gene flow. In accordance with the central-marginal model, we predicted a decrease in population genetic diversity and an increase in differentiation along the latitudinal gradient from the boreal mixed-wood to northern coniferous forest. A total of 24 eastern white cedar populations were sampled along the north-south latitudinal gradient for microsatellite genotyping analysis. Positive F_is values and heterozygote deficiency were observed in populations from the marginal (F_is = 0.244; P_HW = 0.0042) and discontinuous zones (F_is = 0.166; P_HW = 0.0042). However, populations from the continuous zone were in HW equilibrium (F_is = −0.007; P_HW = 0.3625). There were no significant latitudinal effects on gene diversity (H_s), allelic richness (AR), or population differentiation (F_st). Bayesian and NJT (neighbor-joining tree) analyses demonstrated the presence of a population structure that was partly consistent with the geographic origins of the populations. The impact of population fragmentation on the genetic structure of EWC is to create a positive inbreeding coefficient, which was two to three times higher on average than that of a population from the continuous zone. This result indicated a higher occurrence of selfing within fragmented EWC populations coupled with a higher degree of gene exchange among near-neighbor relatives, thereby leading to significant inbreeding. Increased population isolation was apparently not correlated with a detectable effect on genetic diversity. Overall, the fragmented populations of EWC appear well-buffered against effects of inbreeding on genetic erosion.     

Comparative Analysis of Genome Diversity in Bullmastiff Dogs

Management and preservation of genomic diversity in dog breeds is a major objective for maintaining health. The present study was undertaken to characterise genomic diversity in Bullmastiff dogs using both genealogical and molecular analysis. Genealogical analysis of diversity was conducted using a database consisting of 16,378 Bullmastiff pedigrees from year 1980 to 2013. Additionally, a total of 188 Bullmastiff dogs were genotyped using the 170,000 SNP Illumina CanineHD Beadchip. Genealogical parameters revealed a mean inbreeding coefficient of 0.047; 142 total founders (f); an effective number of founders (f_e) of 79; an effective number of ancestors (f_a) of 62; and an effective population size of the reference population of 41. Genetic diversity and the degree of genome-wide homogeneity within the breed were also investigated using molecular data. Multiple-locus heterozygosity (MLH) was equal to 0.206; runs of homozygosity (ROH) as proportion of the genome, averaged 16.44%; effective population size was 29.1, with an average inbreeding coefficient of 0.035, all estimated using SNP Data. Fine-scale population structure was analysed using NETVIEW, a population analysis pipeline. Visualisation of the high definition network captured relationships among individuals within and between subpopulations. Effects of unequal founder use, and ancestral inbreeding and selection, were evident. While current levels of Bullmastiff heterozygosity, inbreeding and homozygosity are not unusual, a relatively small effective population size indicates that a breeding strategy to reduce the inbreeding rate may be beneficial.     

Genomic inbreeding estimation in small populations: evaluation of runs of homozygosity in three local dairy cattle breeds

In the local breeds with small population size, one of the most important problems is the increase of inbreeding coefficient (F). High levels of inbreeding lead to reduced genetic diversity and inbreeding depression. The availability of high-density single nucleotide polymorphism (SNP) arrays has facilitated the quantification of F by genomic markers in farm animals. Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes and represent an estimate of the degree of autozygosity at genome-wide level. The current study aims to quantify the genomic F derived from ROH (F_ROH) in three local dairy cattle breeds. F_ROH values were compared with F estimated from the genomic relationship matrix (F_GRM), based on the difference between observed v. expected number of homozygous genotypes (F_HOM) and the genomic homozygosity of individual i (F_{MOL i}). The molecular coancestry coefficient (f_{MOL ij}) between individuals i and j was also estimated. Individuals of Cinisara (71), Modicana (72) and Reggiana (168) were genotyped with the 50K v2 Illumina BeadChip. Genotypes from 96 animals of Italian Holstein cattle breed were also included in the analysis. We used a definition of ROH as tracts of homozygous genotypes that were >4 Mb. Among breeds, 3661 ROH were identified. Modicana showed the highest mean number of ROH per individual and the highest value of F_ROH, whereas Reggiana showed the lowest ones. Differences among breeds existed for the ROH lengths. The individuals of Italian Holstein showed high number of short ROH segments, related to ancient consanguinity. Similar results showed the Reggiana with some extreme animals with segments covering 400 Mb and more of genome. Modicana and Cinisara showed similar results between them with the total length of ROH characterized by the presence of large segments. High correlation was found between F_HOM and F_ROH ranged from 0.83 in Reggiana to 0.95 in Cinisara and Modicana. The correlations among F_ROH and other estimated F coefficients were generally lower ranged from 0.45 (F_{MOL i}−F_ROH) in Cinisara to 0.17 (F_GRM−F_ROH) in Modicana. On the basis of our results, recent inbreeding was observed in local breeds, considering that 16 Mb segments are expected to present inbreeding up to three generations ago. Our results showed the necessity of implementing conservation programs to control the rise of inbreeding and coancestry in the three Italian local dairy cattle breeds.     

Microsatellite support for active inbreeding in a cichlid fish

In wild animal populations, the degree of inbreeding differs between species and within species between populations. Because mating with kin often results in inbreeding depression, observed inbreeding is usually regarded to be caused by limited outbreeding opportunities due to demographic factors like small population size or population substructuring. However, theory predicts inclusive benefits from mating with kin, and thus part of the observed variation in inbreeding might be due to active inbreeding preferences. Although some recent studies indeed report kin mating preferences, the evidence is still highly ambiguous. Here, we investigate inbreeding in a natural population of the West African cichlid fish Pelvicachromis taeniatus which showed clear kin mating preferences in standardized laboratory experiments but no inbreeding depression. The presented microsatellite analysis reveals that the natural population has, in comparison to two reference populations, a reduced allelic diversity (A = 3) resulting in a low heterozygosity (H_o = 0.167) pointing to a highly inbred population. Furthermore, we found a significant heterozygote deficit not only at population (F_is = 0.116) but also at subpopulation level (F_is = 0.081) suggesting that inbreeding is not only a by-product of population substructuring but possibly a consequence of behavioral kin preferences.     

Coho Salmon (Oncorhynchus kisutch) Prefer and Are Less Aggressive in Darker Environments

Fish are capable of excellent vision and can be profoundly influenced by the visual properties of their environment. Ambient colours have been found to affect growth, survival, aggression and reproduction, but the effect of background darkness (i.e., the darkness vs. lightness of the background) on preference and aggression has not been evaluated systematically. One-hundred Coho salmon (Oncorhynchus kisutch), a species that is increasing in popularity in aquaculture, were randomly assigned to 10 tanks. Using a Latin-square design, every tank was bisected to allow fish in each tank to choose between all the following colour choices (8 choices in total): black vs. white, light grey, dark grey, and a mixed dark grey/black pattern, as well as industry-standard blue vs. white, light grey, dark grey, and black. Fish showed a strong preference for black backgrounds over all other options (p < 0.01). Across tests, preference strength increased with background darkness (p < 0.0001). Moreover, having darker backgrounds in the environment resulted in less aggressive behaviour throughout the tank (p < 0.0001). These results provide the first evidence that darker tanks are preferred by and decrease aggression in salmonids, which points to the welfare benefits of housing farmed salmon in enclosures containing dark backgrounds.     

A spectral theory for Wright’s inbreeding coefficients and related quantities

Wright’s inbreeding coefficient, F_ST, is a fundamental measure in population genetics. Assuming a predefined population subdivision, this statistic is classically used to evaluate population structure at a given genomic locus. With large numbers of loci, unsupervised approaches such as principal component analysis (PCA) have, however, become prominent in recent analyses of population structure. In this study, we describe the relationships between Wright’s inbreeding coefficients and PCA for a model of K discrete populations. Our theory provides an equivalent definition of F_ST based on the decomposition of the genotype matrix into between and within-population matrices. The average value of Wright’s F_ST over all loci included in the genotype matrix can be obtained from the PCA of the between-population matrix. Assuming that a separation condition is fulfilled and for reasonably large data sets, this value of F_ST approximates the proportion of genetic variation explained by the first (K − 1) principal components accurately. The new definition of F_ST is useful for computing inbreeding coefficients from surrogate genotypes, for example, obtained after correction of experimental artifacts or after removing adaptive genetic variation associated with environmental variables. The relationships between inbreeding coefficients and the spectrum of the genotype matrix not only allow interpretations of PCA results in terms of population genetic concepts but extend those concepts to population genetic analyses accounting for temporal, geographical and environmental contexts.     

Inbreeding depression in intra-provenance crosses driven by founder relatedness in white spruce

Genetic means for height growth differed between intra- and inter-provenance crosses, which we hypothesized was due partly to unidentified relatedness among intra-provenance base parents resulting in mild inbreeding and inbreeding depression among their offspring. A dense array of 5,844 single nucleotide polymorphisms was used to directly construct a genomic relationship matrix (G) that had four elements ranging from 0.17 to 0.24, between five intra-provenance base parents. Adjusting the numerator relationship matrix for this relatedness among base parents produced inbreeding coefficients of F _i ≈ 0.1 in their offspring, which displayed depressed height growth. Accounting for inbreeding level as a covariate in a mixed model decreased grossly overestimated (up to 2 ×) dominance variance in models without the covariate adjustment. Height growth decreased 39 cm (～ 6 %) for every 0.1 increase in F _i .     

Using genomic information for management planning of an endangered perennial,Viola uliginosa

Species occupying habitats subjected to frequent natural and/or anthropogenic changes are a challenge for conservation management. We studied one such species, Viola uliginosa, an endangered perennial wetland species typically inhabiting sporadically flooded meadows alongside rivers/lakes. In order to estimate genomic diversity, population structure, and history, we sampled five sites in Finland, three in Estonia, and one each in Slovenia, Belarus, and Poland using genomic SNP data with double‐digest restriction site‐associated DNA sequencing (ddRAD‐seq). We found monophyletic populations, high levels of inbreeding (mean population F_SNP = 0.407–0.945), low effective population sizes (N_e = 0.8–50.9), indications of past demographic expansion, and rare long‐distance dispersal. Our results are important in implementing conservation strategies for V. uliginosa, which should include founding of seed banks, ex situ cultivations, and reintroductions with individuals of proper origin, combined with continuous population monitoring and habitat management.     

Measuring individual inbreeding in the age of genomics: marker-based measures are better than pedigrees

Inbreeding (mating between relatives) can dramatically reduce the fitness of offspring by causing parts of the genome to be identical by descent. Thus, measuring individual inbreeding is crucial for ecology, evolution and conservation biology. We used computer simulations to test whether the realized proportion of the genome that is identical by descent (IBD_G) is predicted better by the pedigree inbreeding coefficient (F_P) or by genomic (marker-based) measures of inbreeding. Genomic estimators of IBD_G included the increase in individual homozygosity relative to mean Hardy–Weinberg expected homozygosity (F_H), and two measures (F_ROH and F_E) that use mapped genetic markers to estimate IBD_G. IBD_G was more strongly correlated with F_H, F_E and F_ROH than with F_P across a broad range of simulated scenarios when thousands of SNPs were used. For example, IBD_G was more strongly correlated with F_ROH, F_H and F_E (estimated with ⩾10 000 SNPs) than with F_P (estimated with 20 generations of complete pedigree) in populations with a recent reduction in the effective populations size (from N_e=500 to N_e=75). F_ROH, F_H and F_E generally explained >90% of the variance in IBD_G (among individuals) when 35 K or more SNPs were used. F_P explained <80% of the variation in IBD_G on average in all simulated scenarios, even when pedigrees included 20 generations. Our results demonstrate that IBD_G can be more precisely estimated with large numbers of genetic markers than with pedigrees. We encourage researchers to adopt genomic marker-based measures of IBD_G as thousands of loci can now be genotyped in any species.     

Genetically-determined hyperfunction of the S100B/RAGE axis is a risk factor for aspergillosis in stem cell transplant recipients

Invasive aspergillosis (IA) is a major threat to the successful outcome of hematopoietic stem cell transplantation (HSCT), although individual risk varies considerably. Recent evidence has established a pivotal role for a danger sensing mechanism implicating the S100B/receptor for advanced glycation end products (RAGE) axis in antifungal immunity. The association of selected genetic variants in the S100B/RAGE axis with susceptibility to IA was investigated in 223 consecutive patients undergoing HSCT. Furthermore, studies addressing the functional consequences of these variants were performed. Susceptibility to IA was significantly associated with two distinct polymorphisms in RAGE (-374T/A) and S100B (+427C/T) genes, the relative contribution of each depended on their presence in both transplantation counterparts [patient SNP_RAGE, adjusted hazard ratio (HR), 1.97; P = 0.042 and donor SNP_RAGE, HR, 2.03; P = 0.047] or in donors (SNP_S100B, HR, 3.15; P = 7.8e-⁴) only, respectively. Functional assays demonstrated a gain-of-function phenotype of both variants, as shown by the enhanced expression of inflammatory cytokines in RAGE polymorphic cells and increased S100B secretion in vitro and in vivo in the presence of the S100B polymorphism. These findings point to a relevant role of the danger sensing signaling in human antifungal immunity and highlight a possible contribution of a genetically-determined hyperfunction of the S100B/RAGE axis to susceptibility to IA in the HSCT setting.     

Genome-wide population structure and evolutionary history of the Frizarta dairy sheep

In the present study, we used genomic data, generated with a medium density single nucleotide polymorphisms (SNP) array, to acquire more information on the population structure and evolutionary history of the synthetic Frizarta dairy sheep. First, two typical measures of linkage disequilibrium (LD) were estimated at various physical distances that were then used to make inferences on the effective population size at key past time points. Population structure was also assessed by both multidimensional scaling analysis and k-means clustering on the distance matrix obtained from the animals’ genomic relationships. The Wright’s fixation F_ST index was also employed to assess herds’ genetic homogeneity and to indirectly estimate past migration rates. The Wright’s fixation F_IS index and genomic inbreeding coefficients based on the genomic relationship matrix as well as on runs of homozygosity were also estimated. The Frizarta breed displays relatively low LD levels with r² and |Dʹ| equal to 0.18 and 0.50, respectively, at an average inter-marker distance of 31 kb. Linkage disequilibrium decayed rapidly by distance and persisted over just a few thousand base pairs. Rate of LD decay (β) varied widely among the 26 autosomes with larger values estimated for shorter chromosomes (e.g. β=0.057, for OAR6) and smaller values for longer ones (e.g. β=0.022, for OAR2). The inferred effective population size at the beginning of the breed’s formation was as high as 549, was then reduced to 463 in 1981 (end of the breed’s formation) and further declined to 187, one generation ago. Multidimensional scaling analysis and k-means clustering suggested a genetically homogenous population, F_ST estimates indicated relatively low genetic differentiation between herds, whereas a heat map of the animals’ genomic kinship relationships revealed a stratified population, at a herd level. Estimates of genomic inbreeding coefficients suggested that most recent parental relatedness may have been a major determinant of the current effective population size. A denser than the 50k SNP panel may be more beneficial when performing genome wide association studies in the breed.     

Genetic heterogeneity and selection signature at the KIT gene in pigs showing different coat colours and patterns

Mutations in the porcine KIT gene (Dominant white locus) have been shown to affect coat colours and colour distribution in pigs. We analysed this gene in several pig breeds and populations (Sicilian black, completely black or with white patches; Cinta Senese; grey local population; Large White; Duroc; Hampshire; Pietrain; wild boar; Meishan) with different coat colours and patterns, genotyping a few polymorphisms. The 21 exons and parts of the intronic regions were sequenced in these pigs and 69 polymorphisms were identified. The grey-roan coat colour observed in a local grey population was completely associated with a 4-bp deletion of intron 18 in a single copy KIT gene, providing evidence that this mutation characterizes the I^d allele described in the early genetic literature. The white patches observed in black Sicilian pigs were not completely associated with the presence of a duplicated KIT allele (I^p), suggesting that genetic heterogeneity is a possible cause of different coat colours in this breed. Selection signature was evident at the KIT gene in two different belted pig breeds, Hampshire and Cinta Senese. The same mutation(s) may cause the belted phenotype in these breeds that originated in the 18th–19th centuries from English pigs (Hampshire) and in Tuscany (Italy) in the 14th century (Cinta Senese). Phylogenetic relationships of 28 inferred KIT haplotypes indicated two clades: one of Asian origin that included Meishan and a few Sicilian black haplotypes and another of European origin.     

Quantifying realized inbreeding in wild and captive animal populations

Most molecular measures of inbreeding do not measure inbreeding at the scale that is most relevant for understanding inbreeding depression—namely the proportion of the genome that is identical-by-descent (IBD). The inbreeding coefficient F_Ped obtained from pedigrees is a valuable estimator of IBD, but pedigrees are not always available, and cannot capture inbreeding loops that reach back in time further than the pedigree. We here propose a molecular approach to quantify the realized proportion of the genome that is IBD (propIBD), and we apply this method to a wild and a captive population of zebra finches (Taeniopygia guttata). In each of 948 wild and 1057 captive individuals we analyzed available single-nucleotide polymorphism (SNP) data (260 SNPs) spread over four different genomic regions in each population. This allowed us to determine whether any of these four regions was completely homozygous within an individual, which indicates IBD with high confidence. In the highly nomadic wild population, we did not find a single case of IBD, implying that inbreeding must be extremely rare (propIBD=0–0.00094, 95% CI). In the captive population, a five-generation pedigree strongly underestimated the average amount of realized inbreeding (F_Ped=0.013<propIBD=0.064), as expected given that pedigree founders were already related. We suggest that this SNP-based technique is generally useful for quantifying inbreeding at the individual or population level, and we show analytically that it can capture inbreeding loops that reach back up to a few hundred generations.