Copy number variations in the amylase gene (AMY2B) in Japanese native dog breeds

A recent study suggested that increased copy numbers of the AMY2B gene might be a crucial genetic change that occurred during the domestication of dogs. To investigate AMY2B expansion in ancient breeds, which are highly divergent from modern breeds of presumed European origins, we analysed copy numbers in native Japanese dog breeds. Copy numbers in the Akita and Shiba, two ancient breeds in Japan, were higher than those in wolves. However, compared to a group of various modern breeds, Akitas had fewer copy numbers, whereas Shibas exhibited the same level of expansion as modern breeds. Interestingly, average AMY2B copy numbers in the Jomon‐Shiba, a unique line of the Shiba that has been bred to maintain their appearance resembling ancestors of native Japanese dogs and that originated in the same region as the Akita, were lower than those in the Shiba. These differences may have arisen from the earlier introduction of rice farming to the region in which the Shiba originated compared to the region in which the Akita and the Jomon‐Shiba originated. Thus, our data provide insights into the relationship between the introduction of agriculture and AMY2B expansion in dogs.     

Dietary Variation and Evolution of Gene Copy Number among Dog Breeds

Prolonged human interactions and artificial selection have influenced the genotypic and phenotypic diversity among dog breeds. Because humans and dogs occupy diverse habitats, ecological contexts have likely contributed to breed-specific positive selection. Prior to the advent of modern dog-feeding practices, there was likely substantial variation in dietary landscapes among disparate dog breeds. As such, we investigated one type of genetic variant, copy number variation, in three metabolic genes: glucokinase regulatory protein (GCKR), phytanol-CoA 2-hydroxylase (PHYH), and pancreatic α-amylase 2B (AMY2B). These genes code for proteins that are responsible for metabolizing dietary products that originate from distinctly different food types: sugar, meat, and starch, respectively. After surveying copy number variation among dogs with diverse dietary histories, we found no correlation between diet and positive selection in either GCKR or PHYH. Although it has been previously demonstrated that dogs experienced a copy number increase in AMY2B relative to wolves during or after the dog domestication process, we demonstrate that positive selection continued to act on amylase copy number in dog breeds that consumed starch-rich diets in time periods after domestication. Furthermore, we found that introgression with wolves is not responsible for deterioration of positive selection on AMY2B among diverse dog breeds. Together, this supports the hypothesis that the amylase copy number expansion is found universally in dogs.     

Individual differences in AMY1 gene copy number, salivary α-amylase levels, and the perception of oral starch

Background

The digestion of dietary starch in humans is initiated by salivary α-amylase, an endo-enzyme that hydrolyzes starch into maltose, maltotriose and larger oligosaccharides. Salivary amylase accounts for 40 to 50% of protein in human saliva and rapidly alters the physical properties of starch. Importantly, the quantity and enzymatic activity of salivary amylase show significant individual variation. However, linking variation in salivary amylase levels with the oral perception of starch has proven difficult. Furthermore, the relationship between copy number variations (CNVs) in the AMY1 gene, which influence salivary amylase levels, and starch viscosity perception has not been explored.

Principal Findings

Here we demonstrate that saliva containing high levels of amylase has sufficient activity to rapidly hydrolyze a viscous starch solution in vitro. Furthermore, we show with time-intensity ratings, which track the digestion of starch during oral manipulation, that individuals with high amylase levels report faster and more significant decreases in perceived starch viscosity than people with low salivary amylase levels. Finally, we demonstrate that AMY1 CNVs predict an individual''s amount and activity of salivary amylase and thereby, ultimately determine their perceived rate of oral starch viscosity thinning.

Conclusions

By linking genetic variation and its consequent salivary enzymatic differences to the perceptual sequellae of these variations, we show that AMY1 copy number relates to salivary amylase concentration and enzymatic activity level, which, in turn, account for individual variation in the oral perception of starch viscosity. The profound individual differences in salivary amylase levels and salivary activity may contribute significantly to individual differences in dietary starch intake and, consequently, to overall nutritional status.     

Rethinking the starch digestion hypothesis for AMY1 copy number variation in humans

Alpha‐amylase exists across taxonomic kingdoms with a deep evolutionary history of gene duplications that resulted in several α‐amylase paralogs. Copy number variation (CNV) in the salivary α‐amylase gene (AMY1) exists in many taxa, but among primates, humans appear to have higher average AMY1 copies than nonhuman primates. Additionally, AMY1 CNV in humans has been associated with starch content of diets, and one known function of α‐amylase is its involvement in starch digestion. Thus high AMY1 CNV is considered to result from selection favoring more efficient starch digestion in the Homo lineage. Here, we present several lines of evidence that challenge the hypothesis that increased AMY1 CNV is an adaptation to starch consumption. We observe that α‐ amylase plays a very limited role in starch digestion, with additional steps required for starch digestion and glucose metabolism. Specifically, we note that α‐amylase hydrolysis only produces a minute amount of free glucose with further enzymatic digestion and glucose absorption being rate‐limiting steps for glucose availability. Indeed α‐amylase is nonessential for starch digestion since sucrase‐isomaltase and maltase‐glucoamylase can hydrolyze whole starch granules while releasing glucose. While higher AMY1 CN and CNV among human populations may result from natural selection, existing evidence does not support starch digestion as the major selective force. We report that in humans α‐amylase is expressed in several other tissues where it may have potential roles of evolutionary significance.     

Genome Sequencing Highlights the Dynamic Early History of Dogs

To identify genetic changes underlying dog domestication and reconstruct their early evolutionary history, we generated high-quality genome sequences from three gray wolves, one from each of the three putative centers of dog domestication, two basal dog lineages (Basenji and Dingo) and a golden jackal as an outgroup. Analysis of these sequences supports a demographic model in which dogs and wolves diverged through a dynamic process involving population bottlenecks in both lineages and post-divergence gene flow. In dogs, the domestication bottleneck involved at least a 16-fold reduction in population size, a much more severe bottleneck than estimated previously. A sharp bottleneck in wolves occurred soon after their divergence from dogs, implying that the pool of diversity from which dogs arose was substantially larger than represented by modern wolf populations. We narrow the plausible range for the date of initial dog domestication to an interval spanning 11–16 thousand years ago, predating the rise of agriculture. In light of this finding, we expand upon previous work regarding the increase in copy number of the amylase gene (AMY2B) in dogs, which is believed to have aided digestion of starch in agricultural refuse. We find standing variation for amylase copy number variation in wolves and little or no copy number increase in the Dingo and Husky lineages. In conjunction with the estimated timing of dog origins, these results provide additional support to archaeological finds, suggesting the earliest dogs arose alongside hunter-gathers rather than agriculturists. Regarding the geographic origin of dogs, we find that, surprisingly, none of the extant wolf lineages from putative domestication centers is more closely related to dogs, and, instead, the sampled wolves form a sister monophyletic clade. This result, in combination with dog-wolf admixture during the process of domestication, suggests that a re-evaluation of past hypotheses regarding dog origins is necessary.     

The Shepherds’ Tale: A Genome-Wide Study across 9 Dog Breeds Implicates Two Loci in the Regulation of Fructosamine Serum Concentration in Belgian Shepherds

Diabetes mellitus is a serious health problem in both dogs and humans. Certain dog breeds show high prevalence of the disease, whereas other breeds are at low risk. Fructosamine and glycated haemoglobin (HbA1c) are two major biomarkers of glycaemia, where serum concentrations reflect glucose turnover over the past few weeks to months. In this study, we searched for genetic factors influencing variation in serum fructosamine concentration in healthy dogs using data from nine dog breeds. Considering all breeds together, we did not find any genome-wide significant associations to fructosamine serum concentration. However, by performing breed-specific analyses we revealed an association on chromosome 3 (p_corrected ≈ 1:68 × 10^-6) in Belgian shepherd dogs of the Malinois subtype. The associated region and its close neighbourhood harbours interesting candidate genes such as LETM1 and GAPDH that are important in glucose metabolism and have previously been implicated in the aetiology of diabetes mellitus. To further explore the genetics of this breed specificity, we screened the genome for reduced heterozygosity stretches private to the Belgian shepherd breed. This revealed a region with reduced heterozygosity that shows a statistically significant interaction (p = 0.025) with the association region on chromosome 3. This region also harbours some interesting candidate genes and regulatory regions but the exact mechanisms underlying the interaction are still unknown. Nevertheless, this finding provides a plausible explanation for breed-specific genetic effects for complex traits in dogs. Shepherd breeds are at low risk of developing diabetes mellitus. The findings in Belgian shepherds could be connected to a protective mechanism against the disease. Further insight into the regulation of glucose metabolism could improve diagnostic and therapeutic methods for diabetes mellitus.     

Studies of modern Italian dog populations reveal multiple patterns for domestic breed evolution

Through thousands of years of breeding and strong human selection, the dog (Canis lupus familiaris) exists today within hundreds of closed populations throughout the world, each with defined phenotypes. A singular geographic region with broad diversity in dog breeds presents an interesting opportunity to observe potential mechanisms of breed formation. Italy claims 14 internationally recognized dog breeds, with numerous additional local varieties. To determine the relationship among Italian dog populations, we integrated genetic data from 263 dogs representing 23 closed dog populations from Italy, seven Apennine gray wolves, and an established dataset of 161 globally recognized dog breeds, applying multiple genetic methods to characterize the modes by which breeds are formed within a single geographic region. Our consideration of each of five genetic analyses reveals a series of development events that mirror historical modes of breed formation, but with variations unique to the codevelopment of early dog and human populations. Using 142,840 genome‐wide SNPs and a dataset of 1,609 canines, representing 182 breeds and 16 wild canids, we identified breed development routes for the Italian breeds that included divergence from common populations for a specific purpose, admixture of regional stock with that from other regions, and isolated selection of local stock with specific attributes.     

Genetic characterization of four native Italian shepherd dog breeds and analysis of their relationship to cosmopolitan dog breeds using microsatellite markers

Very little research into genetic diversity of Italian native dog breeds has been carried out so far. In this study we aimed to estimate and compare the genetic diversity of four native Italian shepherd dog breeds: the Maremma, Bergamasco, Lupino del Gigante and Oropa shepherds. Therefore, some cosmopolitan dog breeds, which have been widely raised in Italy for a long time past, have also been considered to check possible influence of these dog populations on the Italian autochthonous breeds considered here. A total of 212 individuals, belonging to 10 different dog breeds, were sampled and genotyped using 18 autosomal microsatellite loci. We analyzed the genetic diversity of these breeds, within breed diversity, breed relationship and population structure. The 10 breeds considered in this study were clearly genetically differentiated from each other, regardless of current population sizes and the onset of separate breeding history. The level of genetic diversity explained 20% of the total genetic variation. The level of H_E found here is in agreement with that found by other studies. The native Italian breeds showed generally higher genetic diversity compared with the long established, well-defined cosmopolitan dog breeds. As the Border Collie seems closer to the Italian breeds than the other cosmopolitan shepherd dogs considered here, a possible utilization of this breed to improve working performance in Italian traditional working shepherd dogs cannot be ignored. The data and information found here can be utilized in the organization of conservation programs planned to reduce inbreeding and to minimize loss of genetic variability.     

Genetic variants of glucose phosphate isomerase (E.C. 5.3.1.9) in canine erythrocytes*

Genetic polymorphism of glucose phosphate isomerase (GPI) was found in the erythrocytes of dogs of six Japanese breeds by using starch gel electrophoresis. Analysis of parentage records of dogs revealed that the phenotypic variation of erythrocyte glucose phosphate isomerase was controlled by one autosomal locus with two codominant alleles, GPI^Aand GPI^B. The allele GPI^Bwas observed in the following breeds: San'in-Shiba, Shinshu-Shiba, Shikoku, Kai and Kishu, but not in Hokkaidoes and Akitas. All the dogs belonging to 25 European breeds, 5 oriental or China-origin (except Japan) breeds examined in this experiments had the genotype constitution GPI^A/GPI^A, whereas one Dalmation dog was heterozygous GPI^A/GPI^B.     

Measurements of Salivary Alpha Amylase and Salivary Cortisol in Hominoid Primates Reveal Within-Species Consistency and Between-Species Differences

Salivary alpha amylase (sAA) is the most abundant enzyme in saliva. Studies in humans found variation in enzymatic activity of sAA across populations that could be linked to the copy number of loci for salivary amylase (AMY1), which was seen as an adaptive response to the intake of dietary starch. In addition to diet dependent variation, differences in sAA activity have been related to social stress. In a previous study, we found evidence for stress-induced variation in sAA activity in the bonobos, a hominoid primate that is closely related to humans. In this study, we explored patterns of variation in sAA activity in bonobos and three other hominoid primates, chimpanzee, gorilla, and orangutan to (a) examine if within-species differences in sAA activity found in bonobos are characteristic for hominoids and (b) assess the extent of variation in sAA activity between different species. The results revealed species-differences in sAA activity with gorillas and orangutans having higher basal sAA activity when compared to Pan. To assess the impact of stress, sAA values were related to cortisol levels measured in the same saliva samples. Gorillas and orangutans had low salivary cortisol concentrations and the highest cortisol concentration was found in samples from male bonobos, the group that also showed the highest sAA activity. Considering published information, the differences in sAA activity correspond with differences in AMY1 copy numbers and match with general features of natural diet. Studies on sAA activity have the potential to complement molecular studies and may contribute to research on feeding ecology and nutrition.     

Isolation,characterization and genetic analysis of canine <Emphasis Type="Italic">GATA4</Emphasis> gene in a family of Doberman Pinschers with an atrial septal defect

GATA4 is expressed early in the developing heart where it plays a key role in regulating the expression of genes encoding myocardial contractile proteins. Gene mutations in the human GATA4 have been implicated in various congenital heart defects (CHD), including atrial septal defect (ASD). Although ASD is the third most common CHD in humans, it is generally rare in dogs and cats. There is also no obvious predilection for ASD in dogs and cats, based on sex or breed. However, among dogs, the incidence rate of ASD is relatively high in Samoyeds and Doberman Pinschers, where its inheritance and genetic aetiology are not well understood. In this study, we identified and investigated the genetic aetiology of an ASD affected family in a pure breed dog population. Although the GATA4 gene was screened, we did not find any mutations that would result in the alteration of the coding sequence and hence, the predicted GATA4 structure and function. Although the aetiology of ASD is multifactorial, our findings indicate that GATA4 may not be responsible for the ASD in the dogs used in this study. However, this does not eliminate GATA4 as a candidate for ASD in other dog breeds.     

Validation of a urine metabolome fingerprint in dog for phenotypic classification

Selective breeding of dogs over hundreds of years has inadvertently resulted in breed-specific propensities to particular diseases. Furthermore, it has likely induced more subtle affects on the physiology of certain breeds and moved them from their evolutionary optima. In the absence of obvious disease phenotypes such subtle changes could have yet unrecognised breed-specific implications for health and well-being. Here we have applied NMR metabolomics as a discovery-driven approach to identify the impact of breed on the urinary profile of dog and to determine if non-disease-related breed differences can be identified. Multiple urines were collected non-invasively over a two-week period from seven neutered male Labrador retrievers and miniature Schnauzers. Following NMR analyses by 1-dimensional ¹H and 2-dimensional ¹H J-resolved (JRES) spectroscopy, principal component analysis revealed that the metabolic variability within each individual is relatively small compared to inter-individual variability, and that some separation between breeds was evident. A supervised model, using partial least squares discriminant analysis (PLS-DA) with class based upon breed, was trained using the JRES data. The model predicted correctly the breed of seven additional urines, yielding a model sensitivity and specificity of 100%. Several significant metabolic differences between the breeds were identified. A second model was developed using PLS-DA with class based upon individual dogs, which again achieved high classification accuracy for the test set. Overall, this confirms that canine urine is information-rich and that breed is a major determinant of urinary metabolic fingerprints. In the future this may enable a more accurate development of specific nutritional care for an individual or breed.     

Genetic characterisation of a domestic dog<Emphasis Type="Italic">Canis familiaris</Emphasis> breed endemic to South African rural areas

Allozyme electrophoresis (horizontal starch gel and PAGE) and histochemical staining techniques were used to study the genetic composition of an endemic southern African domestic dogCanis familiaris Linnaeus, 1758, the Africanis breed. Genetic differentiation was analysed at 21 protein-coding loci. The results were compared to those for three other populations/breeds: representatives of established Western breeds, crossbred dogs of Western descent from rural areas in South Africa, and indigenous Saluki dogs from the Middle East. Nine polymorphic loci were found (Ak-1,-2, Ck, Per, Hb, Po-A-1 to-3 andPo-Tf). Two unique alleles at theCk andPo-A-2 loci separated the Africanis breed from the other groups. There were also significant differences between Africanis and the other breeds in pair-wise comparisons of allelic frequencies at polymorphic loci. An assignment test, fixation index values, gene flow and genetic distance values indicated a closer genetic association between the Africanis and Saluki breeds than with dogs of Western origin. This finding supports archaeological evidence that the endemic Africanis breed was introduced from the Middle East into Africa thousands of years ago, and not through later western influences. The average heterozygosity ranged from 0.106–0.15, with least heterozygosity in the Africanis and most in the rural crossbred group. The percentage of polymorphic loci, the mean number of alleles per locus (biologically more significant than heterozygosity), and conformation of genotypes to Hardy-Weinberg proportions showed no evidence of recent loss of genetic diversity in Africanis. Genetic differentiation and support of archaeological evidence by genetics indicate that the endemic southern African domestic dog breed is unique.     

Phylogenetic distinctiveness of Middle Eastern and Southeast Asian village dog Y chromosomes illuminates dog origins

Modern genetic samples are commonly used to trace dog origins, which entails untested assumptions that village dogs reflect indigenous ancestry or that breed origins can be reliably traced to particular regions. We used high-resolution Y chromosome markers (SNP and STR) and mitochondrial DNA to analyze 495 village dogs/dingoes from the Middle East and Southeast Asia, along with 138 dogs from >35 modern breeds to 1) assess genetic divergence between Middle Eastern and Southeast Asian village dogs and their phylogenetic affinities to Australian dingoes and gray wolves (Canis lupus) and 2) compare the genetic affinities of modern breeds to regional indigenous village dog populations. The Y chromosome markers indicated that village dogs in the two regions corresponded to reciprocally monophyletic clades, reflecting several to many thousand years divergence, predating the Neolithic ages, and indicating long-indigenous roots to those regions. As expected, breeds of the Middle East and East Asia clustered within the respective regional village dog clade. Australian dingoes also clustered in the Southeast Asian clade. However, the European and American breeds clustered almost entirely within the Southeast Asian clade, even sharing many haplotypes, suggesting a substantial and recent influence of East Asian dogs in the creation of European breeds. Comparison to 818 published breed dog Y STR haplotypes confirmed this conclusion and indicated that some African breeds reflect another distinct patrilineal origin. The lower-resolution mtDNA marker consistently supported Y-chromosome results. Both marker types confirmed previous findings of higher genetic diversity in dogs from Southeast Asia than the Middle East. Our findings demonstrate the importance of village dogs as windows into the past and provide a reference against which ancient DNA can be used to further elucidate origins and spread of the domestic dog.     

Molecular structure in peripheral dog breeds: Portuguese native breeds as a case study

Genetic variability in purebred dogs is known to be highly structured, with differences among breeds accounting for ∼30% of the genetic variation. However, analysis of the genetic structure in non-cosmopolitan breeds and local populations is still limited. Nine Portuguese native dog breeds, and other peripheral dog populations (five) with regional affinities, were characterized using 16 microsatellites and 225 amplified fragment length polymorphism (AFLP) markers, and the pattern of genetic differentiation was investigated. Although the level of breed differentiation detected is below that of other dog breeds, there is in most cases a correlation between breed affiliation and molecular structure. AFLP markers and Bayesian clustering methods allowed an average of 73.1% of individuals to be correctly assigned to source populations, providing robust genotypic assessment of breed affiliation. A geographical genetic structure was also detected, which suggests a limited influence of African dogs on the Iberian breeds. The sampling effect on the estimation of population structure was evaluated and there was a 2.2% decrease in genetic differentiation among breeds when working animals were included. Genetic diversity of stray dogs was also assessed and there is no evidence that they pose a threat to the preservation of the gene pool of native dog breeds.     

Similar recent selection criteria associated with different behavioural effects in two dog breeds

Selection during the last decades has split some established dog breeds into morphologically and behaviourally divergent types. These breed splits are interesting models for behaviour genetics since selection has often been for few and well‐defined behavioural traits. The aim of this study was to explore behavioural differences between selection lines in golden and Labrador retriever, in both of which a split between a common type (pet and conformation) and a field type (hunting) has occurred. We hypothesized that the behavioural profiles of the types would be similar in both breeds. Pedigree data and results from a standardized behavioural test from 902 goldens (698 common and 204 field) and 1672 Labradors (1023 and 649) were analysed. Principal component analysis revealed six behavioural components: curiosity, play interest, chase proneness, social curiosity, social greeting and threat display. Breed and type affected all components, but interestingly there was an interaction between breed and type for most components. For example, in Labradors the common type had higher curiosity than the field type (F_1,1668 = 18.359; P < 0.001), while the opposite was found in goldens (F_1,897 = 65.201; P < 0.001). Heritability estimates showed considerable genetic contributions to the behavioural variations in both breeds, but different heritabilities between the types within breeds was also found, suggesting different selection pressures. In conclusion, in spite of similar genetic origin and similar recent selection criteria, types behave differently in the breeds. This suggests that the genetic architecture related to behaviour differs between the breeds.     

The Kintamani dog: genetic profile of an emerging breed from Bali, Indonesia

The Kintamani dog is an evolving breed indigenous to the Kintamani region of Bali. Kintamani dogs cohabitate with feral Bali street dogs, although folklore has the breed originating 600 years ago from a Chinese Chow Chow. The physical and personality characteristics of the Kintamani dog make it a popular pet for the Balinese, and efforts are currently under way to have the dog accepted by the Federation Cynologique Internationale as a recognized breed. To study the genetic background of the Kintamani dog, 31 highly polymorphic short tandem repeat markers were analyzed in Kintamani dogs, Bali street dogs, Australian dingoes, and nine American Kennel Club (AKC) recognized breeds of Asian or European origin. The Kintamani dog was identical to the Bali street dog at all but three loci. The Bali street dog and Kintamani dog were most closely aligned with the Australian dingo and distantly related to AKC recognized breeds of Asian but not European origin. Therefore, the Kintamani dog has evolved from Balinese feral dogs with little loss of genetic diversity.     

Cloning and characterization of a third type of human α-amylase gene, AMY2B

We have previously reported concerning the existence of a third type of human α-amylase gene, AMY3 [Emi et al., Gene 62 (1988) 229–235; Tomita et al., Gene 76 (1989) 11–18], which is expressed in a lung carcinoid tissue, and differs in nucleotide sequence from the two previously characterized human α-amylase genes coding for salivary and pancreatic isozymes, termed AMY1 and AMY2, respectively.Here, we rename this gene AMY2B to coincide with the designation by Gumucio et al. [Mol. Cell Biol. 8 (1988) 1197–1205] and describe its genetic properties as revealed by sequencing studies. It consists of ten major exons whose sequences are highly homologous to those of AMY1 and AMY2. Not only the exons, but also most of the introns seem to be highly conserved, as judged from physical mapping data. The AMY2B gene identified from mRNA in a lung carcinoid tissue has at least two additional untranslated exons in its 5′ region; hence the promoter lies far upstream relative to the other two AMY genes.     

Genotype imputation in the domestic dog

Application of imputation methods to accurately predict a dense array of SNP genotypes in the dog could provide an important supplement to current analyses of array-based genotyping data. Here, we developed a reference panel of 4,885,283 SNPs in 83 dogs across 15 breeds using whole genome sequencing. We used this panel to predict the genotypes of 268 dogs across three breeds with 84,193 SNP array-derived genotypes as inputs. We then (1) performed breed clustering of the actual and imputed data; (2) evaluated several reference panel breed combinations to determine an optimal reference panel composition; and (3) compared the accuracy of two commonly used software algorithms (Beagle and IMPUTE2). Breed clustering was well preserved in the imputation process across eigenvalues representing 75 % of the variation in the imputed data. Using Beagle with a target panel from a single breed, genotype concordance was highest using a multi-breed reference panel (92.4 %) compared to a breed-specific reference panel (87.0 %) or a reference panel containing no breeds overlapping with the target panel (74.9 %). This finding was confirmed using target panels derived from two other breeds. Additionally, using the multi-breed reference panel, genotype concordance was slightly higher with IMPUTE2 (94.1 %) compared to Beagle; Pearson correlation coefficients were slightly higher for both software packages (0.946 for Beagle, 0.961 for IMPUTE2). Our findings demonstrate that genotype imputation from SNP array-derived data to whole genome-level genotypes is both feasible and accurate in the dog with appropriate breed overlap between the target and reference panels.