Next Article in Journal
Biofortification of Hens Eggs with Polyunsaturated Fatty Acids by New Dietary Formulation: Supercritical Microalgal Extract
Previous Article in Journal
Validation of a New Resource-Efficient Feeding System for Fattening Pigs Using Increased Crude Fiber Concentrations in Diets: Feed Intake and Ammonia Emissions
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 10
Issue 3
10.3390/ani10030498
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The German Shorthair Pointer Dog Breed (Canis lupus familiaris): Genomic Inbreeding and Variability

by
Antonio Boccardo
1,
Stefano Paolo Marelli
1,
Davide Pravettoni
1,
Alessandro Bagnato
1,
Giuseppe Achille Busca
2 and
Maria Giuseppina Strillacci
1,*
1
Department of Veterinary Medicine, Università degli Studi di Milano, Via dell’Università 6, 26900 Lodi, Italy
2
Centro Clinico-Veterinario e Zootecnico-Sperimentale, Università degli Studi di Milano, Via dell’Università 6, 26900 Lodi, Italy
*
Author to whom correspondence should be addressed.
Animals 2020, 10(3), 498; https://doi.org/10.3390/ani10030498
Submission received: 24 February 2020 / Revised: 11 March 2020 / Accepted: 13 March 2020 / Published: 17 March 2020
(This article belongs to the Section Animal Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

In order to protect domestic animals’ biodiversity, a deep knowledge of the genomic makeup is required. The authors describe the genomic architecture of the German Short Hair Pointer breed and analyze the inbreeding levels under a genomic and a genealogic perspective. Twenty-four dogs from Italy were genotyped and analyzed jointly with 10 dogs from USA, whose genotypes were available from a published research. The authors investigated the genomic structural variation of the breed using runs of homozygosity—the direct measurement of the proportion of homozygous DNA, i.e., genomic inbreeding. Some traits clearly revealed the selection objectives addressed in the breed. The results describe a low inbred population with quite good levels of genetic variability.

Abstract

The German Shorthaired Pointer (GSHP) is a breed worldwide known for its hunting versatility. Dogs of this breed are appreciated as valuable companions, effective trackers, field trailers and obedience athletes. The aim of the present work is to describe the genomic architecture of the GSHP breed and to analyze inbreeding levels under a genomic and a genealogic perspective. A total of 34 samples were collected (24 Italian, 10 USA), and the genomic and pedigree coefficients of inbreeding have been calculated. A total of 3183 runs of homozygosity (ROH) across all 34 dogs have been identified. The minimum and maximum number of Single Nucleotide Polymorphisms (SNPs) defining all ROH are 40 and 3060. The mean number of ROH for the sample was 93.6. ROH were found on all chromosomes. A total of 854 SNPs (TOP_SNPs) defined 11 ROH island regions (TOP_ROH), in which some gene already associated with behavioral and morphological canine traits was annotated. The proportion of averaged observed homozygotes estimated on total number of SNPs was 0.70. The genomic inbreeding coefficient based on ROH was 0.17. The mean inbreeding based on genealogical information resulted 0.023. The results describe a low inbred population with quite a good level of genetic variability.

1. Introduction

Dog selection is based on the breed standards that usually define the selection objectives for morphological traits and the behavioral predisposition of each breed [1]. Going back in human history, the domestication of dogs determined a genetic divergence from wild wolves, and the subsequent directional selection for specific characteristics led to a variety of dog breeds, counted to be around 500 [2]. Directional selection leads to a multitude of phenotypic variants for morphologic and behavioral traits [3,4].
Physical and physiological characteristics are determined by the genetic composition of the breeds which may have been affected by population bottlenecks due to specific directional selection practices. In this occurrence, a reduced number of ancestors, jointly to the large use of popular sires (e.g., the top winning male of the year) may result in an increase in homozygosity at population level [3,5,6].
The creation in the mid-19th century of kennel clubs that defined breeds’ standards somehow formalized the reproductive isolation within the multiple canine morphologies present at that time [2]. A “breed-barrier rule” was then created, defining that every puppy may not enter a studbook unless its parents are registered in the same one [2]. The Kennel clubs’ studbooks activity in breeds preservation is based on strict pedigree control and evaluations of closeness to breeds standards. This process may lead to the fast reduction in the genetic variability within breeds, unless specific reproductive plans to control inbreeding are applied [7,8,9,10].
Previous studies have shown a decline in the genetic variability of some modern dog populations, characterized by high inbreeding levels like the Lundehund population where the pedigree-based inbreeding was 0.10 calculated on 11 generations [11,12,13,14,15]. The evaluation of genetic diversity represents a pivotal step for the establishment of a breeding program or to take conservation measures. An increase in homozygosity often leads to a loss in biological fitness due to a concomitant increase in the likelihood of the offspring being affected by deleterious or recessive traits [15]. Published studies report lower genetic diversity in breeds used for shows with respect to working dog breeds. These latter seem to maintain better genetic health due to the presence of many genetic traits still similar to those of their wild ancestors [13,16,17,18]. The characterization of genetic architecture, the study of inbreeding levels and the attention to the genetic health of dog in general should be considered the basis of effective selection plans, based on objective scientific information [16].
The traditional approach to inbreeding coefficient (IC) calculation is based on pedigree information: this value is clearly dependent on entries’ data accuracy and known generations in the pedigree (pedigree depth) [9,19]. The application of molecular genetics, and the construction of a high-density SNP maps [3] provide a powerful tool for the genomic IC calculation based on molecular information. Genomic IC can be measured via runs of homozygosity (ROH) (FROH), representing a direct measurement of the fractions of homozygous stretches more likely to be transmitted identically by descent [9,10,16,18,19,20,21]. Another method is based on Wright’s F statistics, FIS that returns the proportion of observed homozygosity respect the expected one [22]. The same powerful genomic tools can be used to provide molecular information to study the differentiation among breeds [20].
During domestication, hunting ability could have been considered among the most appreciated traits to select for; evidences have been found that 4000 years ago in the Middle East and North Africa, selective breeding for hunting aptitude was applied [23,24]. The typical hunting behavior of pointing breeds is represented by an interruption of the natural wolf predatory sequence, stopping in front of their prey [25,26]. Furthermore, the introduction of firearms in hunting drove a double goal in dogs’ performance ability: both pointing the pray and keeping away from the rifle [24].
The German Shorthaired Pointer (GSHP) breed represents the optimal evolution of pointing ability in German dog breeds. According to FCI standards, the history of the breed traces back to those dogs from Mediterranean countries that were used to hunt feathered game with nets and with raptors. Those dogs arrived in the German courts where they were selected primarily for their hunting ability and their versatility as gun dogs following the evolution of fire arms. The studbook of the breed, the “Zuchtbook Deutch-Kurzhaar” was firstly published in 1897. The standard of the breed—defining morphology and working characteristics—was written by Prince Albrecht zu Solms-Braunfeld. In Italy, the introduction of the breed began in the 1930s. GSHPs were exported to North America in the 1920s, and they entered the AKC studbook in 1930; recognized worldwide for their hunting versatility. Furthermore, they are appreciated as valuable companions, effective trackers, field trailers and obedience athletes. [27].
The aim of this study is to describe the genomic architecture of GSHP breed and to analyze the inbreeding levels of the breed under a genomic and a genealogic perspective, to supply effective data helpful for the design of breeding plans and conservation projects.

2. Materials and Methods

2.1. Sampling and Genotyping

The DNA of twenty-four pedigree GSHP samples (12 males, 12 female) was extracted from residual blood samples collected for routine screening of health status and according to the University Ethical Committee opinion n. 2/16 on re-use of collected samples. The samples were taken from unrelated, second-generation dogs from different Italian regions. This approach was aimed to make the sample here considered as much as possible representative of the Italian population and of its genetic variability. According to storage protocols, the dogs’ identities were certified by a veterinarian after the chip control: chip number, ENCI studbook registration number, name, sex and date of birth are the data commonly filled in the automated format used by the ENCI official labs in dog identification and recording procedures.
All dogs were genotyped with the Illumina CanineHD BeadChip array containing about 220,853 SNPs and the CanFam3.1 genome assembly (GCF_000002285.3) was considered as reference genome.
Ten genotyped German Shorthair Pointer dogs (GSHP_USA) from [8] (173662 SNPs - GSE90441 code project) were used in this study in order to increase the number of samples and to perform a better evaluation of selection targets related to breed morphology and pointing abilities.
A final SNP dataset of 102,045 autosomal markers was obtained merging the two genotypes dataset: only SNP on autosomes were considered (from 1 to 38) and SNP genotypes data were filtered for minor allele frequency (MAF) (value ≤ 0.01) and for call rate (0.99), in order to reduce the bias ascribable to missing genotypes in inbreeding coefficients evaluation and in other statistics.
ADMIXTURE (v. 1.3.0) software [28] has been used to estimate the individual ancestries in the considered populations (the Italian GSHP_ITALY and the American GSHP_USA ones), and then to calculate the number of underlying possible subpopulations. ADMIXTURE was run from K = 1 to K = 4, and the optimal number of clusters (K-value) was determined as the one having the lowest cross-validation error (--cv flag added in command line). In addition, a Principal component analysis (PCA), based on the allele frequencies of SNPs, using SVS 8.4 software (SVS) (Golden Helix Inc., Bozeman, MT, USA) was performed.

2.2. Runs of homozygosity (ROH)

SVS software has been used to perform ROH detection. The ROH were defined setting a minimum of 1000 kb in size and 40 homozygous SNPs, no heterozygote SNPs are permitted, no missing SNPs were allowed in the ROH, and a maximum gap between SNPs of 1000 Kb was predefined in order to assure that the SNP density did not affect the ROH.
ROHs were grouped into 5 classes of length (0-2 Mb, 2-4 Mb, 4-8 Mb, 8-16 Mb and, >16 Mb) and all the statistics were calculated across individuals within country (n = 24 for GSHP_ITALY and n = 10 for GSHP_USA) and for overall_GSHP (all 34 dogs).
The genomic regions with the highest frequency of ROH (TOP_ROH), and thus potentially under selection, have been identified by selecting the SNPs most commonly occurring in ROH, i.e., those occurring in at least 50% of samples plus one. The SNPs in the TOP_ROH are here referred to as TOP_SNPs. The TOP_SNPs and the TOP_ROH were annotated on CanFam3.1 genome (GCA_000002285.2, Release 99) using the Variant Effect Predictor (VEP) and the BioMart tools of Ensembl genome browser [29], respectively. Only genes with an official gene name were considered. Functional analysis of annotated genes has been performed using DAVID 6.8 database [30].
Graphical representations of ROH were obtained using the R package DetectRUNs [31].

2.3. Inbreeding Coefficients

In this study, genomic molecular inbreeding coefficients (FHOM and FROH) have been calculated for the overall_GSHP and the genealogical F for GSHP_ITALY (FPED) from pedigree information available only for Italian GSHP population:
(i) Inbreeding coefficients based on the excess in the observed number of homozygous genotypes (FHOM) with SVS, following [22]
FHOM = (HomOb − HomEx)/(n. of non-missing-data loci − HomEx)
where HomOb and HomEx are the observed and expected numbers of homozygous genotypes in a dog.
(ii) Inbreeding coefficient based on ROH (FROH) with DetectRun, according to the following formula:
FROH = LROH/Laut
where LROH is the total length of all ROHs of an individual, Laut the specified length of the autosomal genome (from 1 to 38) covered by the SNPs used in this study (2,201,412,378 bp).
(iii) Inbreeding coefficient based on pedigree information (FPED) was calculated using Pedigree Viewer software [32]. Genealogical data of each sample dog (4 known generations) have been obtained using ENCI (Ente Nazionale della Cinofilia Italiana, FCI).

3. Results

ADMIXTURE analysis revealed that the GSHP_ITALY and GSHP_USA shared the same unique ancestor—as the lowest CV value has been obtained with K = 1 (Figure 1A)—and are considered then a unique breed. In addition, the PCA supported the results of the ADMIXTURE: as shown in Figure 1B, all individuals belong to a unique cluster. An example of a GSHP dog is shown in Figure 2C.
All the results of this study are presented for overall_GSHP (n = 34 GSHP dogs).

3.1. Runs of Homozygosity (ROH)

The SVS software identified a total of 3183 runs across all 34 dogs (Table S1). Results revealed marked differences in terms of the number and length of ROH across individuals.
The minimum and maximum number of SNPs defining all ROH are 40 and 3060, respectively.
ROH are mainly short in length; in fact, the ROH of 0-2 Mb and 2–4 Mb are the most frequent classes of length identified (i.e., 70%) (Figure 2A). ROH were also found within the >16 Mb length class (Figure 2A).
The number of ROH per individual ranged from 71 to 131, with a mean number of ROH for sample of 93.6 (Table 1). Only three samples showed a very large number of ROH, with counts of 123, 124 and 131. The average size of the ROH of these three individuals is nevertheless similar to the one of the other dogs (5551411 bp, 2,902,950 bp and 3,203,416 bp). Figure 2B shows the relationship between ROH count and the average total length of ROH for each individual (mainly ranged 3 to 5 Mb). The amount of the genome covered by ROH per dog ranged (as mean values) from 2,557,174 bp to 7,312,103 bp.
ROH were found on all chromosomes: less than 50% of chromosomes (n.17) (Figure 2C) have a mean length over 4.072 Mb (mean ROH length, as reported in Table 1), and no evident correlation between chromosomes length and mean ROH length resulted. A graphical representation of ROH frequencies on autosomes is shown in Figure 2D. The number of ROH within each class of length per chromosome was also calculated (Table S2): the longest ROH (>16 Mb) were identified on almost all the chromosomes, except for chr23, chr27, chr32, chr33, chr34, and chr36. A graphical representation of ROH statistics for Italian and USA populations are shown separately in Figure S1 A–E.
The Figure 3A shows the SNP occurrences in ROH segments across the genome, highlighting also that the genomic distribution of ROH is clearly non-uniform across autosomes. A total of 11 TOP_ROH on 9 chromosomes were identified in at least 18 samples (50% of sampled dogs +1) (red line in Figure 3A) (Table 2). The higher chromosomal peaks were identified on CFA 3, 10 and 13. Table 2 also reports the number of samples in which the TOP_ROH have been identified (min = 18 and max = 29), the list of genes (n.116) present in the identified TOP_ROH, and the list of genes (n. 15) for which an already known association with canine traits has been described (Table 2 – in bold, and detailed in Table S6).
We found a total of 854 TOP_SNPs, and Figure 3B is the graphical representation of the annotated position of these SNPs (11 classes of position) according to Ensembl VEP. The major part of SNPs mapped in intergenic (n. 401; 46,96%) and in intronic positions (n. 401; 46,96%) (Figure 3B). Details of all the TOP_SNPs positions and annotation information are reported in Table S4.
The highest number of homozygotes SNPs (> 25) was annotated within the FMN1 (n. 26), RYR3 (n. 28), CSMD3 (n.29), OPCML (n. 33), and NTM (n. 50) genes. In addition, four SNPs resulted annotated in missense positions: BICF2P1347925, BICF2S23147347 and TIGRP2P48933_rs8545710 on chr3 and BICF2P847459 on chr30 of HTT, HTT, ADD1, and RYR3 genes, respectively. SNP genotypes of the dogs for each of these positions were verified and, as reported in Table S3, resulted in: all wildtype genotype for BICF2P1347925, TIGRP2P48933_rs8545710 and BICF2P847459 (n. 22, 22 and 20, respectively) with AA variant coded as AA genotype; all mutated variants were found for BICF2S23147347 (n. 22) with variant CC coded as BB genotype. Table S3 also includes genotypes for these missense positions for the dogs that did not contribute to these ROH.
Functional classification of genes annotated in TOP_ROH (Table 3) provided by DAVID database (112 gene IDs recognized) revealed that these genes were significantly enriched (nominal p-value <0.05) in three GO terms in biological processes, only one GO term in cellular components and seven KEGG pathways, for which the CACNG5, CACNG4, CACNG1 are the main involved genes (Table 3).
Table S5 reports the GO terms and KEGG pathways resulted with nominal p-value >0.05.

3.2. Inbreeding Coefficients

The average observed and expected homozygotes calculated using 102,045 SNPs were 71785.44 and 70513.85, respectively. The proportion of average observed and expected homozygotes estimated on the total number of SNPs were then 0.70 and 0.69.
The inbreeding coefficient estimated from SNP markers (FHOM, equation 1) was 0.04 (averaged values) and individual values ranged from −0.06 to 0.21 for dogs with lower than average homozygosity and vice versa, respectively.
Instead, the inbreeding coefficients based on ROH (FROH, equation 2) were slightly higher: a mean value of 0.17, ranging between 0.09 to 0.32. Differences in FROH were found along all chromosomes (Figure S2). The higher FROH values (≥ 0.25) have been identified for chrs 9, 13 and 25.
The correlation coefficient, and that of determination (R 2) calculated between FHOM and FROH, were 0.988 and 0.975 (as reported in Figure 4), respectively.
The four-generation pedigree analysis resulted in a very low FPED coefficient for all the samples: 14 subjects showed no genealogical inbreeding, only four GSHPs revealed a FPED higher than 0.05, the maximum FPED was 0.161 (Table S6). The mean FPED of the studied population was 0.023 with a standard deviation of 0.043. The regression between FHOM and FPED (Figure S3A) showed a similar relationship to and FPED and FROH, with a regression coefficient of 0.83 vs. 0.84, but with a lower coefficient of determination (R 2 of 0.427 vs. 0.97). Figure S3B shows the regression between FROH and FPED, resulting in a lower coefficient but similar R2 to FHOM and FPED.

4. Discussion

Not many studies have explored the ROH pattern and its possible association with inbreeding depression, to morphological and health traits in dog breeds and populations in comparison to researches performed and available in other species (i.e., cattle). The ROH patterns and their distributions in specific dog breeds remain largely unexplored.
To the best of our knowledge, this study is the first attempt to characterize the ROH distribution in the German Shorthair Point dog breed using canine high-density SNP arrays. The genealogical pedigree information was used to sample the 24 individuals that were as un-related as possible. Ten available online genotyped GSHP samples from the GSPH_USA populations were available from a previous published study [8] and added to the 24 Italian dogs in order to increase the sample size of this study. According to the results of ADMIXTURE and the PCA which hereinbefore presented the detection of ROH are here presented jointly. Details of the ROH for GSPH_USA and GSPH_Italy are shown in Figure S1 A–E. The results obtained allow us to disclose the variability in the GSPH breed and to identify the genomic region harboring ROH, even if enlarging the sample size in future efforts may disclose additional information on genomic variation in this population.
The detection of genomic regions (i.e., TOP_ROH) underlying breed-specific phenotypic characteristics or attitude, can provide indication about a specific genetic structure of a population or breed. ROH, in concordance with their length, inform on a possible artificial or natural selection pressure on specific genomic tracts and on inbreeding levels also providing information on bottlenecks to which the populations have been subjected over the years. Long ROH (~ 10 Mb) occur as a result of recent inbreeding and when recombination events do not interrupt long chromosome segments. Short ROH (~ 1 Mb), instead, are produced by IBD (identical by descent) genomic regions from old ancestors (up to 50 generations ago) [33].
Similar to what has been identified in some hunting dog breeds—the French Pointing dog type the Pyrenee and in Rhodesian Ridgeback dogs, respectively—most of the ROH identified in this study were short in length (0-2 Mb, with a mean length of 1.42 Mb), suggesting that dogs of these populations are involved in more ancient relatedness [18,34]. The occurrence of TOP_ROH hotspots in genomic regions that harbor candidate genes may be involved in directional selection pressure. Out of 854 TOP_SNPs, 452 are annotated within genes according with VEP tool (taking into account different positions, as in Figure 3A). The most representative genes (those in which at least 25 TOP_SNPs map) resulted NTM, OPCML, CSMD3, RYR3, FMN1, TMTC2, and HTT. All these genes (except for NTM), together with others (n. 15) have been already associated with dog traits, comprising behavior and morphological traits, as reported in Table 4.
NTM and OPCML (the most representative) are genes involved in central nervous system functioning. Studies performed in humans proposed that the NTM gene (Neurotrimin) is associated with IQ level and cognitive function performances. In addition, Gurgul et al. [37] supposed in their study on diversifying selection signature between draft and light horses, the potential contributes of MNT and OPCML (Opioid Binding Protein/Cell Adhesion Molecule Like) genes to the differing horses’ temperaments and to the ability to develop different gaits as a function of motor coordination. It could be supposed that these genes’ expression, considering the presence of many TOP_SNPs, could be related to the peculiar hunting abilities and style of GSHPs. These traits have been strongly selected by the breeders since the 19th century, and they are clearly described in the breed working standards [26,35]. RYR3 (Ryanodine Receptor 3), together with GRK4 (G protein-coupled receptor kinase 4) located in TOP_ROH_02, resulted genes positively selected in athletic dog breeds (sports and hunting aptitudes). In fact, these two genes could have an important role in field activity: RYR3 because of its function in the activation of muscle skeletal contraction through the coordinate activation of voltage dependent Ca2+ [38], and GRK4 in increasing cardiac output.
The only gene for which an SNP missense position (BICF2S23147347) with alternative variant G (coded as B allele) has been found (Table S4), according to VEP tool, is the HTT gene. No specific association study results are yet available for this mutation or gene in the dog species. The only exception is for chasing behavior trait, as reported in Table 4.
Another gene that we could hypothesize under selection in this breed, according to breed-specific characteristics described in the standard, is that the encoding factors involved in skeletal muscle contraction processes—i.e., CACNG1 (calcium voltage-gated channel auxiliary subunit gamma 1) —: a polymorphism of CACNG1 gene has been associated with elite strength athlete status [39]. CACNG1, CACNG4 (subunit gamma 4), and CACNG5 (subunit gamma 5) are genes implicated in five KEGG pathways related to cardiac functions (contraction and healthy). These three genes are also part of the MAPK signaling pathway, active on the proliferation of mammalian cells, and play an important role at different levels—e.g., the regulation of various diseases (cardiovascular), taking part in anti-inflammatory effects, responses to stress, protection against injury, and the maintenance of gastrointestinal functions [40].
The loss of genetic variability and the increase in inbreeding levels coupled with the consequent reduction in the effective population size in a high number of pedigree canine populations could be considered the main risk factor in conservation project definitions [12,15,41,42,43,44]. In general, the reported results show low FPED in the studied subjects [44]. The effectiveness of the molecular dissection of inbreeding is shown by the FROH (equation 2) results, which reveal the genomic relationship among the studied individuals as demonstrated by different researchers [8,16,45]. The proportion between observed and expected homozygosity reveals similar results to those reported by [18] about the French Pointing dog Type Pyrenee (1.03 vs. 1.02). The obtained inbreeding coefficient FHOM (equation 1) based on SNPs is lower compared to all the results reported in 11 breeds by [8], which ranged from 0.179 in the Papillon to 0.53 in the Basenji. When considering the inbreeding coefficient based on ROH the obtained results are the same of those presented by [34] with a registered level of 0.17. The same coefficient calculated in the French Pointing dog Type Pyrenee was slightly lower: 0.112 [18]. Considering pedigree data analysis, the obtained results in the studied population are slightly lower when compared to other pointing breeds registered in the Italian studbook (ENCI): in the Italian Pointing Dog, the pedigree average inbreeding coefficient was 0.041 [46], and in the French Pointing Dog type Pyrenees F was calculated to be 0.033 [47]—in both cases, the whole registered population was considered. The results we obtained on a sample population are close to those reported by [8], who compared pedigree average genealogical inbreeding with SNP chip and whole genome sequence (WGS) heterozygosity: the closeness of our results it is not only related to the pedigree (five generations) inbreeding coefficient of very famous breeds known worldwide like the Labrador Retriever (0.026), Golden Retriever (0.027) and Bernese Mountain Dog (0.022), but the same differences reported comparing pedigree inbreeding with SNPs and WGS inbreeding levels has been described in the Italian GSHP sample population. Inbreeding levels considering pedigree data were considerably lower than the average inbreeding levels obtained by molecular genomic data analysis.

5. Conclusions

The investigation of the genomic architecture and the molecular dissection of inbreeding play a pivotal role in the conservation and protection plans for animal populations. Furthermore, genomic analysis supplies effective tools with high quality molecular data in the demographic profiling of canine populations. The presented results represent the first description of the molecular dissection of the inbreeding levels in the Italian and USA pointing GSHP breeds. The genomic molecular investigation allowed the authors to finely describe the genetic assets of the breed, even if enlarging the sample size may disclose other regions under recent inbreeding. The obtained results could be considered a powerful and objective tool in breed conservation management. The genomic data here analysed showed important homozygosity for morphological and behavioural traits, whose expression is contributed by known annotated genes. The findings of the research also underline the importance of integrating genealogical and molecular information on dog breeds, in order to obtain meaningful data for breeds’ protection actions.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2615/10/3/498/s1, Figure S1: Graphical representation of ROH statistics for Italian and USA GSHP. Figure S2: Graphical representation of FROH value for each chromosome. Figure S3A: Regression and coefficient of determination (R2) calculated between FHOM (Y) and FPED (X); Figure S3B: Regression and coefficient of determination (R2) calculated between FROH (Y) and FPED (X). Table S1: List of ROH identified in GSHP dog breed. Table S2: ROH count per chromosomes and per classes of ROH. Table S3: TOP_SNPs annotation. Table S4: Samples genotypes per TOP_SNPs missense positions. Table S5: Functional annotation of genes mapping in TOP_ROH according with DAVID database (nominal p-value >0.05).

Author Contributions

Conceptualization, A.B. (Antonio Boccardo); methodology, S.P.M. and M.G.S.; formal analysis, S.P.M., A.B. (Alessandro Bagnato) and M.G.S.; investigation, M.G.S.; data curation, A.B. (Antonio Boccardo), S.P.M., D.P., A.B. (Alessandro Bagnato) and G.A.B.; writing—original draft preparation, S.P.M. and M.G.S.; writing—review and editing, A.B. (Antonio Boccardo), D.P., A.B. (Alessandro Bagnato) and G.A.B.; supervision, M.G.S.; funding acquisition, A.B. (Antonio Boccardo). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Università degli Studi di Milano (Linea 2, Azione A; grant code: LINEA2_ABOCC_AA_2017).

Acknowledgments

The authors would like to thank Ente Nazionale della Cinofilia Italiana and Kurzhaar Club Italia for data supplying and the Lombardy delegation of the Kurzhaar Club Italia for the support in developing the study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Asher, L.; Diesel, G.; Summers, J.F.; McGreevy, P.D.; Collins, L.M. Inherited defects in pedigree dogs. Part 1: Disorders related to breed standards. Vet. J. 2009, 182, 402–411. [Google Scholar] [CrossRef]
  2. Parker, H.G.; Kim, L.V.; Sutter, N.B.; Carlson, S.; Lorentzen, T.D.; Malek, T.B.; Johnson, G.S.; DeFrance, H.B.; Ostrander, E.A.; Kruglyak, L. Genetic structure of the purebred domestic dog. Science 2004, 304, 1160–1164. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Lindblad-Toh, K.; Wade, C.M.; Mikkelsen, T.S.; Karlsson, E.K.; Jaffe, D.B.; Kamal, M.; Clamp, M.; Chang, J.L.; Kulbokas, E.J., III; Zody, M.C. Genome sequence, comparative analysis and haplotype structure of the domestic dog. Nature 2005, 438, 803–819. [Google Scholar] [CrossRef] [PubMed]
  4. Schoenebeck, J.J.; Ostrander, E.A. Insights into morphology and disease from the dog genome project. Ann. Rev. Cell Dev. Biol. 2014, 30, 535–560. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Pollinger, J.P.; Lohmueller, K.E.; Han, E.; Parker, H.G.; Quignon, P.; Degenhardt, J.D.; Boyko, A.R.; Earl, D.A.; Auton, A.; Reynolds, A. Genome-wide SNP and haplotype analyses reveal a rich history underlying dog domestication. Nature 2010, 464, 898–902. [Google Scholar]
  6. Vilà, C.; Leonard, J.A. Canid phylogeny and origin of the domestic dog. In The Genetics of the Dog; CABI: Wallingford, UK, 2012. [Google Scholar]
  7. Dreger, D.L.; Davis, B.W.; Cocco, R.; Sechi, S.; Di Cerbo, A.; Parker, H.G.; Polli, M.; Marelli, S.P.; Crepaldi, P.; Ostrander, E.A. Commonalities in development of pure breeds and population isolates revealed in the genome of the sardinian Fonni’s Dog. Genetics 2016, 204, 737–755. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Dreger, D.L.; Rimbault, M.; Davis, B.W.; Bhatnagar, A.; Parker, H.G.; Ostrander, E.A. Whole-genome sequence, SNP chips and pedigree structure: Building demographic profiles in domestic dog breeds to optimize genetic-trait mapping. Dis. Model. Mech. 2016, 9, 1445–1460. [Google Scholar] [CrossRef] [Green Version]
  9. Chu, E.T.; Simpson, M.J.; Diehl, K.; Page, R.L.; Sams, A.J.; Boyko, A.R. Inbreeding depression causes reduced fecundity in Golden Retrievers. Mamm. Genome 2019, 30, 166–172. [Google Scholar] [CrossRef] [Green Version]
  10. Mortlock, S.-A.; Khatkar, M.S.; Williamson, P. Comparative analysis of genome diversity in bullmastiff dogs. PLoS ONE 2016, 11, e0147941. [Google Scholar] [CrossRef] [Green Version]
  11. Pertoldi, C.; Kristensen, T.N.; Loeschcke, V.; Berg, P.; Praebel, A.; Stronen, A.V.; Proschowsky, H.F.; Fredholm, M. Characterization of the genetic profile of five Danish dog breeds. J. Anim. Sci. 2013, 91, 5122–5127. [Google Scholar] [CrossRef] [Green Version]
  12. Pfahler, S.; Distl, O. Effective population size, extended linkage disequilibrium and signatures of selection in the rare dog breed Lundehund. PLoS ONE 2015, 10, e0122680. [Google Scholar] [CrossRef] [PubMed]
  13. Pedersen, N.C.; Pooch, A.S.; Liu, H. A genetic assessment of the English bulldog. Canine Genet. Epidemiol. 2016, 3, 6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Pedersen, N.; Liu, H.; Theilen, G.; Sacks, B. The effects of dog breed development on genetic diversity and the relative influences of performance and conformation breeding. J. Anim. Breed. Genet. 2013, 130, 236–248. [Google Scholar] [CrossRef]
  15. Keijser, S.F.A.; Fieten, H.; Vos-Loohuis, M.; Piek, C.J.; Anderson, H.; Donner, J.; Scholten, I.; Nielen, M.; Hesselink, J.W.; van Steenbeek, F.G. Heterozygosity testing and multiplex DNA panel screening as a potential tool to monitor health and inbreeding in a small, closed dog population. Canine Genet. Epidemiol. 2018, 5, 12. [Google Scholar] [CrossRef] [PubMed]
  16. Wiener, P.; Sánchez-Molano, E.; Clements, D.N.; Woolliams, J.A.; Haskell, M.J.; Blott, S.C. Genomic data illuminates demography, genetic structure and selection of a popular dog breed. BMC Genom. 2017, 18, 609. [Google Scholar] [CrossRef] [PubMed]
  17. Bigi, D.; Marelli, S.P.; Liotta, L.; Frattini, S.; Talenti, A.; Pagnacco, G.; Polli, M.; Crepaldi, P. Investigating the population structure and genetic differentiation of livestock guard dog breeds. Animal 2018, 12, 2009–2016. [Google Scholar] [CrossRef] [PubMed]
  18. Mastrangelo, S.; Biscarini, F.; Auzino, B.; Ragatzu, M.; Spaterna, A.; Ciampolini, R. Genome-wide diversity and runs of homozygosity in the “Braque Français, type Pyrénées” dog breed. BMC Res. Notes 2018, 11, 13. [Google Scholar] [CrossRef] [Green Version]
  19. Kardos, M.; Luikart, G.; Allendorf, F.W. Measuring individual inbreeding in the age of genomics: Marker-based measures are better than pedigrees. Heredity (Edinb) 2015, 115, 63. [Google Scholar] [CrossRef]
  20. Mastrangelo, S.; Biscarini, F.; Tolone, M.; Auzino, B.; Ragatzu, M.; Spaterna, A.; Ciampolini, R. Genomic characterization of the Braque Français type Pyrénées dog and relationship with other breeds. PLoS ONE 2018, 13, e0208548. [Google Scholar] [CrossRef]
  21. Mortlock, S.-A.; Booth, R.; Mazrier, H.; Khatkar, M.S.; Williamson, P. Visualization of genome diversity in german shepherd dogs. Bioinform. Biol. Insights 2015, 9, BBI-S30524. [Google Scholar] [CrossRef] [Green Version]
  22. Wright, S. Genetics of populations. Encycl. Br. 1948, 10, 111-AD. [Google Scholar]
  23. Clutton-Brock, J. A Natural History of Domesticated Mammals; Cambridge University Press: Cambridge, UK, 1999; ISBN 0521634954. [Google Scholar]
  24. Parra, D.; Méndez, S.; Canon, J.; Dunner, S. Genetic differentiation in pointing dog breeds inferred from microsatellites and mitochondrial DNA sequence. Anim. Genet. 2008, 39, 1–7. [Google Scholar] [CrossRef] [PubMed]
  25. Udell, M.A.R.; Ewald, M.; Dorey, N.R.; Wynne, C.D.L. Exploring breed differences in dogs (Canis familiaris): Does exaggeration or inhibition of predatory response predict performance on human-guided tasks? Anim. Behav. 2014, 89, 99–105. [Google Scholar] [CrossRef]
  26. Turcsán, B.; Kubinyi, E.; Miklósi, Á. Trainability and boldness traits differ between dog breed clusters based on conventional breed categories and genetic relatedness. Appl. Anim. Behav. Sci. 2011, 132, 61–70. [Google Scholar] [CrossRef]
  27. Bell, J.; Cavanagh, K.; Tilley, L.; Smith, F.W.K. Veterinary Medical Guide to Dog and Cat Breeds; CRC Press: Boca Raton, FL, USA, 2012; ISBN 1482241412. [Google Scholar]
  28. Alexander, D.H.; Novembre, J.; Lange, K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009, 19, 1655–1664. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. ENSEMBL Genome Browser. Available online: https://www.ensembl.org/index.html (accessed on 16 March 2020).
  30. DAVID 6.8 Data Base. Available online: https://david.ncifcrf.gov (accessed on 16 March 2020).
  31. Biscarini, F.; Cozzi, P.; Gaspa, G.; Marras, G. detectRUNS: Detect Runs of Homozygosity and Runs of Heterozygosity in Diploid Genomes. 2018. Available online: http://orca.cf.ac.uk/id/eprint/108906 (accessed on 16 March 2020).
  32. Kinghorn, B.P.; Kinghorn, A.J. Pedigree Viewer 6.5; University of New England: Armidale, Australia, 2010. [Google Scholar]
  33. Mastrangelo, S.; Portolano, B.; Di Gerlando, R.; Ciampolini, R.; Tolone, M.; Sardina, M.T.; Consortium, I.S.G. Genome-wide analysis in endangered populations: A case study in Barbaresca sheep. Animal 2017, 11, 1107–1116. [Google Scholar] [CrossRef]
  34. Lamac, N.; Sölkner, J.; Mészáros, G. Analysis of Excessive Homozygous Regions in Rhodesian Ridgeback Dogs. Agric. Conspec. Sci. 2017, 82, 277–280. [Google Scholar]
  35. MacLean, E.; Snyder-Mackler, N.; Serpell, J. Highly heritable and functionally relevant breed differences in dog behavior. BioRxiv 2019, 286, 509315. [Google Scholar]
  36. Vaysse, A.; Ratnakumar, A.; Derrien, T.; Axelsson, E.; Pielberg, G.R.; Sigurdsson, S.; Fall, T.; Seppälä, E.H.; Hansen, M.S.T.; Lawley, C.T. Identification of genomic regions associated with phenotypic variation between dog breeds using selection mapping. PLoS Genet. 2011, 7, e1002316. [Google Scholar] [CrossRef] [Green Version]
  37. Gurgul, A.; Jasielczuk, I.; Semik-Gurgul, E.; Pawlina-Tyszko, K.; Stefaniuk-Szmukier, M.; Szmatoła, T.; Polak, G.; Tomczyk-Wrona, I.; Bugno-Poniewierska, M. A genome-wide scan for diversifying selection signatures in selected horse breeds. PLoS ONE 2019, 14, e0210751. [Google Scholar] [CrossRef] [Green Version]
  38. Kim, J.; Williams, F.J.; Dreger, D.L.; Plassais, J.; Davis, B.W.; Parker, H.G.; Ostrander, E.A. Genetic selection of athletic success in sport-hunting dogs. Proc. Natl. Acad. Sci. USA 2018, 115, E7212–E7221. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Maciejewska-Skrendo, A.; Cięszczyk, P.; Chycki, J.; Sawczuk, M.; Smółka, W. Genetic Markers Associated with Power Athlete Status. J. Hum. Kinet. 2019, 68, 17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Du, W.; Hu, H.; Zhang, J.; Bao, G.; Chen, R.; Quan, R. The Mechanism of MAPK Signal Transduction Pathway Involved with Electroacupuncture Treatment for Different Diseases. Evidence-Based Complement. Altern. Med. 2019, 2019, 8138017. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Kettunen, A.; Daverdin, M.; Helfjord, T.; Berg, P. Cross-breeding is inevitable to conserve the highly inbred population of puffin hunter: The Norwegian Lundehund. PLoS ONE 2017, 12, e0170039. [Google Scholar] [CrossRef]
  42. Leroy, G. Genetic diversity, inbreeding and breeding practices in dogs: Results from pedigree analyses. Vet. J. 2011, 189, 177–182. [Google Scholar] [CrossRef] [PubMed]
  43. Leroy, G.; Mary-Huard, T.; Verrier, E.; Danvy, S.; Charvolin, E.; Danchin-Burge, C. Methods to estimate effective population size using pedigree data: Examples in dog, sheep, cattle and horse. Genet. Sel. Evol. 2013, 45, 1. [Google Scholar] [CrossRef] [Green Version]
  44. Jansson, M.; Laikre, L. Pedigree data indicate rapid inbreeding and loss of genetic diversity within populations of native, traditional dog breeds of conservation concern. PLoS ONE 2018, 13, e0202849. [Google Scholar] [CrossRef] [Green Version]
  45. Talenti, A.; Dreger, D.L.; Frattini, S.; Polli, M.; Marelli, S.; Harris, A.C.; Liotta, L.; Cocco, R.; Hogan, A.N.; Bigi, D. Studies of modern Italian dog populations reveal multiple patterns for domestic breed evolution. Ecol. Evol. 2018, 8, 2911–2925. [Google Scholar] [CrossRef]
  46. Cecchi, F.; Paci, G.; Spaterna, A.; Ciampolini, R. Genetic variability in Bracco Italiano dog breed assessed by pedigree data. Ital. J. Anim. Sci. 2013, 12, e54. [Google Scholar] [CrossRef] [Green Version]
  47. Cecchi, F.; Paci, G.; Spaterna, A.; Ragatzu, M.; Ciampolini, R. Demographic approach on the study of genetic parameters in the dog Braque Français type Pyrénées italian population. Ital. J. Anim. Sci. 2016, 15, 30–36. [Google Scholar] [CrossRef]
Animals 10 00498 g001 550
Figure 1. A) Cross-validation (CV) distribution calculated using in ADMIXTURE. Plot display CV errors versus K 1-4. B) PCA results and samples distribution (PC_2 = 0.67; PC_3 = 0.65). C) Example of GSHP dog (female).
Figure 1. A) Cross-validation (CV) distribution calculated using in ADMIXTURE. Plot display CV errors versus K 1-4. B) PCA results and samples distribution (PC_2 = 0.67; PC_3 = 0.65). C) Example of GSHP dog (female).
Animals 10 00498 g001
Animals 10 00498 g002 550
Figure 2. Graphical representation of runs of homozygosity (ROH) statistics. A) Frequencies and counts of ROH for each class of length; B) Relationship between number and averaged total length (Mb) of ROH in each dog; C) Mean length (Mb) for each chromosome; D) Frequencies of ROH per chromosomes.
Figure 2. Graphical representation of runs of homozygosity (ROH) statistics. A) Frequencies and counts of ROH for each class of length; B) Relationship between number and averaged total length (Mb) of ROH in each dog; C) Mean length (Mb) for each chromosome; D) Frequencies of ROH per chromosomes.
Animals 10 00498 g002
Animals 10 00498 g003 550
Figure 3. SNPs occurrences in identified ROH and their annotation on canine genome assembly CanFam3.1: A) Manhattan plot: red line indicates the adopted threshold: 50% of samples + 1; B) Annotated position of TOP_SNPs: numbers in legend of graph represent percentages of each identified position respect to a gene.
Figure 3. SNPs occurrences in identified ROH and their annotation on canine genome assembly CanFam3.1: A) Manhattan plot: red line indicates the adopted threshold: 50% of samples + 1; B) Annotated position of TOP_SNPs: numbers in legend of graph represent percentages of each identified position respect to a gene.
Animals 10 00498 g003
Animals 10 00498 g004 550
Figure 4. Regression and coefficient of determination (R2) calculated between FHOM and FROH: y = FROH; x = FHOM.
Figure 4. Regression and coefficient of determination (R2) calculated between FHOM and FROH: y = FROH; x = FHOM.
Animals 10 00498 g004
Table
Table 1. Descriptive statistics of ROHs.
Table 1. Descriptive statistics of ROHs.
POPN. ROHsMin - Max (Mean) of ROH Per Samples Min ROH LengthMax ROH LengthMean ROH Length(%) Total Coverage
overall_GSHP318371-131 (93.6)10007175689302540724075.9
Table
Table 2. TOP_ROHs found in at least 50% of overall_GSHP samples +1 (n = 18) and annotated genes.
Table 2. TOP_ROHs found in at least 50% of overall_GSHP samples +1 (n = 18) and annotated genes.
TOP_ROH_idChrStart PositionEnd Position*N. SamplesGenes
TOP_ROH_1111200964811232418318CNFN, MEGF8, TMEM145, PRR19, PAFAH1B3, CIC, ERF, GSK3A, ZNF526, DEDD2, POU2F2, ZNF574, GRIK5
TOP_ROH_23608929116282414822-24DOK7, HGFAC, RGS12, MSANTD1, HTT, GRK4, NOP14, MFSD10, ADD1, SH3BP2, TNIP2, FAM193A, RNF4, CFAP99, ZFYVE28, MXD4, HAUS3, POLN, NAT8L, C3H4orf48, NELFA, SCARNA22, NSD2, LETM1, FGFR3, TACC3, TMEM129, SLBP, FAM53A, NKX1-1, UVSSA, MAEA, FAM184B
TOP_ROH_35247752348097618-23B3GAT1, GLB1L2 (LOC606786), GLB1L3, ACAD8, THYN1, cfa-mir-8868, VPS26B, NCAPD3, JAM3, IGSF9B, SPATA19, OPCML, NTM
TOP_ROH_454900070497914818PRDM10, NFRKB
TOP_ROH_59125396691466812118KPNA2, C9H17orf58, BPTF, NOL11, PITPNC1, PSMD12, HELZ, CACNG1, CACNG4, CACNG5, PRKCA, APOH, CEP112, AXIN2
TOP_ROH_6107278956841877122-27RASSF3, GNS, TBC1D30, WIF1, LEMD3, MSRB3, HMGA2
TOP_ROH_7132962719450560521-29ZNF706, GRHL2, NCALD, RRM2B, UBR5, ODF1, KLF10, AZIN1
TOP_ROH_813101513571303698918-20SYBU, KCNV1, CSMD3
TOP_ROH_915242048302508930618CCDC59, METTL25, TMTC2
TOP_ROH_10252268042388738618-20UFM1, TRPC4, POSTN, SUPT20H, EXOSC8, ALG5, SMAD9, RFXAP
TOP_ROH_1130927947277566718-20EMC7, CHRM5, AVEN, RYR3, FMN1, GREM1, SCG5, ARHGAP11A, GJD2, ACTC1, AQR, ZNF770
(*) = min and max number of samples with TOP_ROH; in Italics = genes in which map the highest number of TOP_SNPs; in bold = genes already associated with canine traits.
Table
Table 3. Functional classification of genes according with DAVID database: BP = Biological Process, CC = Cellular Component: KEGG = KEGG pathways.
Table 3. Functional classification of genes according with DAVID database: BP = Biological Process, CC = Cellular Component: KEGG = KEGG pathways.
CategoryTerm (Nominal p-Value)Genes
GO:0019226~transmission of nerve impulse (0.005)CHRM5, CACNG5, JAM3
GOTERM (BP) GO:0071407~cellular response to organic cyclic compound (0.02)SMAD9, GSK3A, AXIN2
GO:0043161~proteasome-mediated ubiquitin-dependent protein catabolic process (0.035)MAEA, RNF4, PSMD12, GSK3A
GOTERM (CC)GO:0005654~nucleoplasm (0.01)HAUS3, SMAD9, HTT, NFRKB, RRM2B, GRHL2, BPTF, UBR5, SYBU, POLN, PITPNC1, AXIN2, TNIP2, KPNA2, DEDD2, ADD1
KEGGcfa04260: Cardiac muscle contraction (0.005);
cfa05410: Hypertrophic cardiomyopathy (HCM) (0.007);
cfa05414: Dilated cardiomyopathy (0.008);
cfa04261: Adrenergic signaling in cardiomyocytes (0.028).		ACTC1, CACNG5, CACNG4, CACNG1
cfa04921: Oxytocin signaling pathway (0.005);
cfa04010: MAPK signaling pathway (0.032).		PRKCA, RYR3, CACNG5, CACNG4, CACNG1
cfa05412: Arrhythmogenic right ventricular cardiomyopathy (ARVC) (0.042)CACNG5, CACNG4, CACNG1
Table
Table 4. Genes already associated with canine traits in the scientific literature.
Table 4. Genes already associated with canine traits in the scientific literature.
TOP_ROH_idGene NameBehavior Traits [35] Morphological Traits [36]
TOP_ROH_2FAM193AChasing, Dog Aggression
TOP_ROH_2HTTChasing
TOP_ROH_2RNF4Chasing, Dog Aggression
TOP_ROH_3NCAPD3Separation Problems
TOP_ROH_3SPATA19Chasing
TOP_ROH_5NOL11Energy
TOP_ROH_5OPCMLDog Aggression
TOP_ROH_6HMGA2 Body size
TOP_ROH_6MSRB3Dog FearEar morphology
TOP_ROH_6TBC1D30Attachment/attention-seeking, Dog Fear
TOP_ROH_6WIF1 Ear morphology
TOP_ROH_8CSMD3Excitability
TOP_ROH_9TMTC2Attachment/attention-seeking, Excitability
TOP_ROH_11FMN1Energy, Trainability
TOP_ROH_11RYR3attachment/attention-seeking, Energy

Share and Cite

MDPI and ACS Style

Boccardo, A.; Marelli, S.P.; Pravettoni, D.; Bagnato, A.; Busca, G.A.; Strillacci, M.G. The German Shorthair Pointer Dog Breed (Canis lupus familiaris): Genomic Inbreeding and Variability. Animals 2020, 10, 498. https://doi.org/10.3390/ani10030498

AMA Style

Boccardo A, Marelli SP, Pravettoni D, Bagnato A, Busca GA, Strillacci MG. The German Shorthair Pointer Dog Breed (Canis lupus familiaris): Genomic Inbreeding and Variability. Animals. 2020; 10(3):498. https://doi.org/10.3390/ani10030498

Chicago/Turabian Style

Boccardo, Antonio, Stefano Paolo Marelli, Davide Pravettoni, Alessandro Bagnato, Giuseppe Achille Busca, and Maria Giuseppina Strillacci. 2020. "The German Shorthair Pointer Dog Breed (Canis lupus familiaris): Genomic Inbreeding and Variability" Animals 10, no. 3: 498. https://doi.org/10.3390/ani10030498

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop