Identification of genomic regions associated with phenotypic variation between dog breeds using selection mapping

Characterizing differences in the muscle transcriptome between cattle with alternative LCORL-NCAPG haplotypes

BACKGROUND

The LCORL-NCAPG locus is a major quantitative trait locus (QTL) on bovine chromosome 6 (BTA6) that influences growth and carcass composition in cattle. To further understand the molecular mechanism responsible for the phenotypic changes associated with this locus, twenty-four Charolais-sired calves were selected for muscle transcriptome analysis based on alternative homozygous LCORL-NCAPG haplotypes (i.e., 12 "QQ" and 12 "qq", where "Q" is a haplotype harboring variation associated with increased growth). At 300 days of age, a biopsy of the longissimus dorsi muscle was collected from each animal for RNA sequencing.

RESULTS

Gene expression analysis identified 733 genes as differentially expressed between QQ and qq animals (q-value < 0.05). Notably, LCORL and genes known to be important regulators of growth such as IGF2 were upregulated in QQ individuals, while genes associated with adiposity such as FASN and LEP were downregulated, reflecting the increase in lean growth associated with this locus. Gene set enrichment analysis demonstrated QQ individuals had downregulation of pathways associated with adipogenesis, alongside upregulation of transcripts for cellular machinery essential for protein synthesis and energy metabolism, particularly ribosomal and mitochondrial components.

CONCLUSIONS

The differences in the muscle transcriptome between QQ and qq animals imply that muscle hypertrophy may be metabolically favored over accumulation of fat in animals with the QQ haplotype. Our findings also suggest this haplotype could be linked to a difference in LCORL expression that potentially influences the downstream transcriptional effects observed, though further research will be needed to confirm the molecular mechanisms underlying the associated changes in phenotype.

Analysis of canine gene constraint identifies new variants for orofacial clefts and stature

Dog breeding promotes within-group homogeneity through conformation to strict breed standards, while simultaneously driving between-group heterogeneity. There are over 350 recognized dog breeds that provide the foundation for investigating the genetic basis of phenotypic diversity. Typically, breed standard phenotypes such as stature, pelage, and craniofacial structure are analyzed through genetic association studies. However, such analyses are limited to assayed phenotypes only, leaving difficult-to-measure phenotypic subtleties easily overlooked. We investigated coding variation from over 2000 dogs, leading to discoveries of variants related to craniofacial morphology and stature. Breed-enriched variants were prioritized according to gene constraint, which was calculated using a mutation model derived from trinucleotide substitution probabilities. Among the newly found variants is a splice-acceptor variant in PDGFRA associated with bifid nose, a characteristic trait of Çatalburun dogs, implicating the gene's role in midline closure. Two additional LCORL variants, both associated with canine body size are also discovered: a frameshift that causes a premature stop in large breeds (>25 kg) and an intronic substitution found in small breeds (<10 kg), thus highlighting the importance of allelic heterogeneity in selection for breed traits. Most variants prioritized in this analysis are not associated with genomic signatures for breed differentiation, as these regions are enriched for constrained genes intolerant to nonsynonymous variation. This indicates trait selection in dogs is likely a balancing act between preserving essential gene functions and maximizing regulatory variation to drive phenotypic extremes.

Deep sequencing of Korean Jindo dog reveals evolutionary trajectory of coat color variations

The inherent diversity of canines is closely intertwined with the unique color patterns of each dog population. These variations in color patterns are believed to have originated through mutations and selective breeding practices that occurred during and after the domestication of dogs from wolves. To address the significant gaps that persist in comprehending the evolutionary processes that underlie the development of these patterns, we generated and analyzed deep-sequenced genomes of 113 Korean Indigenous Jindo dogs that represent five distinct color patterns to identify the associated mutations in CBD103, ASIP, and MC1R. The degree of linkage disequilibrium and estimated allelic ages consistently indicate that the black-and-tan dogs descend from the first major founding population on Jindo island, compatible with the documented literature. We additionally demonstrate that black-and-tan dogs, in contrast to other color variations within the breed, exhibit a closer genetic affinity to ancient wolves from western Eurasia than those from eastern Eurasia. Lastly, the alleles unique to black-and-tan dogs indicated a pleiotropic effect by significantly decreasing the body size phenotypes compared to non-carriers. Overall, comparisons of whole genome sequences of each coat color population diverged from the same breed provided an unprecedented glimpse into the properties of evolutionary processes maintaining variation in Korean Jindo dog populations.

Measures of Homozygosity and Relationship to Genetic Diversity in the Bearded Collie Breed

Background: Genetic diversity in closed populations, such as pedigree dogs, is of concern for maintaining the health and vitality of the population in the face of evolving challenges. Measures of genetic diversity rely upon estimates of homozygosity without consideration of whether the homozygosity is desirable or undesirable or if heterozygosity has a functional impact. Pedigree coefficients of inbreeding have been the classical approach yet they are inadequate unless based upon the entire population. Methods: Homozygosity measures based upon pedigree analyses (n = 11,898), SNP array data (n = 244), and whole genome sequencing (n = 23) were compared in the Bearded Collie, as well as a comparison of SNP array data to a pedigree cohort (n = 5042) and a mixed-breed cohort (n = 1171). Results: Molecular measures based upon DNA are more informative on an individual's homozygosity levels than pedigree analyses, although SNP coefficients of inbreeding overestimate the level of inbreeding based on the nature of SNP array methodology. Whole genome sequence (WGS) analyses revealed that the heterozygosity observed is generally in variants having neutral or low impact, which would indicate that the variability may not contribute substantially to functional diversity in the population. The majority of high-impact variants were observed in the shortest runs of homozygosity (ROH) reflecting ancestral breeding and domestication practices. As expected, mixed-breed dogs displayed higher measures of genomic diversity than either Bearded Collies or other pedigree dogs as a whole using the current paradigm algorithm models to calculate homozygosity. Conclusions: Using typical DNA-based measures reflect only a single individual and not the population thereby failing to account for regions of homozygosity that reflect ancestral breeding, domestication history, breed-defining regions, or regions positively selected for health traits. Incorporating measures of genetic diversity into dog breeding schemes is meritorious. However, until measures of diversity can distinguish between breed-defining homozygosity and homozygosity associated with positive health alleles, the measures to use as selection tools need refinement before their widespread implementation.

Ancient dog introgression into the Iberian wolf genome may have facilitated adaptation to human-dominated landscapes

Understanding how large carnivores respond to increasingly human-dominated landscapes will determine their future adaptive potential. The Iberian wolf (Canis lupus signatus), a gray wolf subspecies endemic to the Iberian Peninsula (Portugal and Spain), has uniquely persisted in human-dominated landscapes, unlike many other wolf populations that faced widespread extinction across Europe during the twentieth century. In this study, we conducted a comprehensive genome-wide analysis of 145 historical and contemporary Iberian wolf samples to investigate whether hybridization with domestic dogs resulted in genetic introgression. We identified a dog-derived block on Chromosome 2 in Iberian wolves, displaying signatures consistent with introgression and high nucleotide similarity among introgressed individuals. Additionally, our estimates place the average timing of introgression between 6100 and 3000 years ago, with low sequence divergence to dogs from the Iberian Peninsula suggesting a single local origin for the hybridization event. Using forward genetic simulations, we show that the introgressed haplotype is most likely being maintained in Iberian wolves by selection. The introgressed dog variants are located within the MAST4 gene, which has been linked to neurological disorders, including cognitive and motor developmental delays, hinting at a potential role in cognitive behavior in Iberian wolves. This study uncovers a case of putative adaptive introgression from domestic dogs into wolves, offering new insights into wild canids' adaptation to human-dominated landscapes.

Behavioural differences and similarities between dog breeds: proposing an ecologically valid approach for canine behavioural research

The behaviour of dogs holds great relevance for not only scientists from fundamental and applied research areas, but also due to the widespread roles of dogs in our societies as companions and working animals; their behaviour is also an important factor in animal and human welfare. A large proportion of dogs currently under human supervision belong to one of roughly 400 recognised breeds. Dog breeds can be characterised by distinctive, predictable and reproducible features, including some of their behavioural traits. To the scientist, the comparative analysis of the behaviour of dog breeds provides an opportunity for investigating an array of intriguing phenomena within an easily accessible model organism created from natural and human-driven evolutionary processes. There are many ways to design and conduct breed-related behavioural investigations, but such endeavours should always be based around biologically relevant research questions and lead to ecologically valid conclusions. In this review, we surveyed recent research efforts that included dog behaviour-related comparisons and applied a critical evaluation according to their methods of breed choice and the subsequent research design. Our aim was to assess whether these two fundamentally important components of experimental design provide a solid basis to reach valid conclusions. Based on 97 publications that fulfilled our selection criteria, we identified three primary methods used by researchers to select breeds for their investigations: (i) convenience sampling; (ii) hypothesis-driven, ancestry-based sampling; and (iii) hypothesis-driven, functional sampling. By using the SWOT (Strengths, Weaknesses, Opportunities, Threats) evaluation system, we highlight each of these techniques' merits and shortcomings. We identify when particular methods may be inherently unable to produce biologically meaningful results due to a mismatch between breed choice and the initial research goals. We hope that our evaluation will help researchers adopt best practices in experimental design regarding future dog breed comparisons.

Genomic Rewilding of Domestic Animals: The Role of Hybridization and Selection in Wolfdog Breeds

Background/Objectives: The domestication of the grey wolf (Canis lupus) and subsequent creation of modern dog breeds have significantly shaped the genetic landscape of domestic canines. This study investigates the genomic effects of hybridization and breeding management practices in two hybrid wolfdog breeds: the Czechoslovakian Wolfdog (CSW) and the Saarloos Wolfdog (SAW). Methods: We analyzed the genomes of 46 CSWs and 20 SAWs, comparing them to 12 German Shepherds (GSHs) and 20 wolves (WLFs), which served as their ancestral populations approximately 70-90 years ago. Results: Our findings highlight that hybridization can increase genetic variability and mitigate the effects of inbreeding, as evidenced by the observed heterozygosity levels in both wolfdog breeds. However, the SAW genome revealed a higher coefficient of inbreeding and longer runs of homozygosity compared to the CSW, reflecting significant inbreeding during its development. Discriminant Analysis of Principal Components and fixation index analyses demonstrate that the CSW exhibits closer genetic proximity to the GSH than the SAW, likely due to differences in the numbers of GSHs used during their creation. Maximum likelihood clustering further confirmed clear genetic differentiation between these hybrid breeds and their respective ancestors, while shared ancestral polymorphism was detectable in all populations. Conclusions: These results highlight the role of controlled hybridization with captive-bred wolves and peculiar breeding strategies in shaping the genetic structure of wolfdog breeds. To ensure the long-term genetic health of these breeds, it is recommended to promote adequate and sustainable breeding practices to maintain genetic diversity, minimize inbreeding, and incorporate the careful selection of unrelated individuals from diverse lineages, while avoiding additional, uncontrolled crossings with wild wolves.

Comparative Whole-Genome Analysis of Production Traits and Genetic Structure in Baiyu and Chuanzhong Black Goats

Natural selection and artificial breeding are crucial methods for developing new animal groups. The Baiyu black goats and Chuanzhong black goats are indigenous goat breeds from distinct ecological regions in Sichuan Province, with dramatically different growth and reproductivity. This study aimed to systematically elucidate the differences in production performance and genetic traits between Baiyu black goats and Chuanzhong black goats. We quantified growth and reproductive attributes for both breeds. Furthermore, we conducted a comprehensive analysis of genetic diversity, population structure, and selection signatures using whole-genome resequencing data. This dataset included 30 individuals from the Baiyu black goat breed, 41 from the Chuanzhong black goat breed, and an additional 59 individuals representing Chengdu grey goats, Tibetan cashmere goats, and Jianchang black goats, totaling 130 individuals across five goat breeds. The comparative analysis of production performance revealed that the weight and body size of Chuanzhong black goats were significantly higher than those of Baiyu black goats (p < 0.01). At the same time, the average kidding rate and kid-weaning survival rate of Chuanzhong black goats were also notably superior to those of Baiyu black goats (p < 0.01). The Baiyu black goats exhibited a more abundant genetic diversity and distinct genetic differences compared to the Chuanzhong black goat, according to an analysis grounded on genomic variation. The Baiyu black goats are more closely related to Tibetan cashmere goats, whereas Chuanzhong black goats share a closer genetic relationship with Chengdu grey goats. Additionally, we employed the π, Fst, and XP-EHH methodologies to identify genes related to immunity (TRIM10, TRIM15, TRIM26, and TRIM5), neurodevelopment (FOXD4L1, PCDHB14, PCDHB4, PCDHB5, PCDHB6, and PCDHB7), reproduction (BTNL2 and GABBR1), body size (NCAPG, IBSP, and MKNK1), and meat quality traits (SUCLG2 and PGM5). These results provide a theoretical basis for further resource conservation and breeding improvement of the Baiyu black goat and Chuanzhong black goat.

Deciphering the genetic basis of behavioral traits in dogs: Observed-trait GWAS and latent-trait GWAS analysis reveal key genes and variants

Dogs exhibit remarkable phenotypic diversity, particularly in behavioral traits, making them an excellent model for studying the genetic basis of complex behaviors. Behavioral traits such as aggression and fear are highly heritable among different dog breeds, but their genetic basis is largely unknown. We used the genome-wide association study (GWAS) to identify candidate genes associated with nine behavioral traits including; stranger-directed aggression (SDA), owner-directed aggression (ODA), dog-directed aggression (DDA), stranger-directed fear (SDF), nonsocial fear (NF), dog-directed fear (DDF), touch sensitivity (TS), separation-related behavior (SRB) and attachment attention-seeking (AAS). The observed behavioral traits were collected from 38,714 to 40,460 individuals across 108 modern dog breeds. We performed a GWAS based on a latent trait extracted using the confirmatory factor analysis (CFA) method with nine observable behavioral traits and compared the results with those from the GWAS of the observed traits. Using both observed-trait and latent-trait GWAS, we identified 41 significant SNPs that were common between both GWAS methods, of which 26 were pleiotropic, as well as 10 SNPs unique to the latent-trait GWAS, and 5 SNPs unique to the observed-trait GWAS discovered. These SNPs were associated with 21 genes in latent-trait GWAS and 22 genes in the observed-trait GWAS, with 19 genes shared by both. According to previous studies, some of the genes from this study have been reported to be related to behavioral and neurological functions in dogs. In the human population, these identified genes play a role in either the formation of the nervous system or are linked to various mental health conditions. Taken together, our findings suggest that latent-trait GWAS for behavioral traits in dogs identifies significant latent genes that are neurologically prioritized.

MtDNA genetic diversity and phylogeographic insights into giant domestic pigeon (Columba livia domestica) breeds: connections between Central Europe and the Middle East

Humans have selectively bred domestic pigeons (Columba livia domestica) to create breeds with a diversity of shapes, colors and other attributes. Since Darwin, the domestic pigeon has always been a popular model species for scientific research because of its richness of form, colouration and behaviour. It is believed that the world's squab pigeon industry uses breeds and hybrids from the Mediterranean region. An exception is the indigenous giant pigeon breeds of the Carpathian Basin, whose origin is not known. Therefore, our aims were 1) to understand the phylogenetic relationships of giant pigeons, which sheds light on the origin of Hungarian breeds and their relationship to the Mediterranean giant pigeon breed group; 2) to contribute molecular genetic data to the genealogy of 2 Iraqi pigeon breeds close to the pigeon domestication center, including the culturally important Iraqi Red Pigeon, and 3) to compare the genetic diversity of European and Middle Eastern domestic pigeon populations and to draw conclusions on the phylogenetic relationships between pigeon breeds and molecular clues to their different breeding practices of both regions. A 655-bp-long sequence of the cytochrome oxidase 1 (COI) region of the mitochondrial DNA was studied in a total of 276 pigeons (19 breeds). A total of 27 haplotypes were found, of which 22 were unique. The highest genetic diversity was found in the Carpathian Basin, and the lowest in the Iraqi region. STRUCTURE analysis revealed low structurality, K=3 was the most likely. The majority of the samples belong to the most ancient haplotype H_2=219, however the Jacobin pigeon is on a very separate evolutionary branch with a large number of mutations. None of the 19 breeds investigated in this study have been previously studied in phylogenetics, and most of these breeds have potential as squab pigeons, and have good meat forms for utilization, therefore the results of this study may also be of help to the squab pigeon industry.

A map of canine sequence variation relative to a Greenland wolf outgroup

For over 15 years, canine genetics research relied on a reference assembly from a Boxer breed dog named Tasha (i.e., canFam3.1). Recent advances in long-read sequencing and genome assembly have led to the development of numerous high-quality assemblies from diverse canines. These assemblies represent notable improvements in completeness, contiguity, and the representation of gene promoters and gene models. Although genome graph and pan-genome approaches have promise, most genetic analyses in canines rely upon the mapping of Illumina sequencing reads to a single reference. The Dog10K consortium, and others, have generated deep catalogs of genetic variation through an alignment of Illumina sequencing reads to a reference genome obtained from a German Shepherd Dog named Mischka (i.e., canFam4, UU_Cfam_GSD_1.0). However, alignment to a breed-derived genome may introduce bias in genotype calling across samples. Since the use of an outgroup reference genome may remove this effect, we have reprocessed 1929 samples analyzed by the Dog10K consortium using a Greenland wolf (mCanLor1.2) as the reference. We efficiently performed remapping and variant calling using a GPU-implementation of common analysis tools. The resulting call set removes the variability in genetic differences seen across samples and breed relationships revealed by principal component analysis are not affected by the choice of reference genome. Using this sequence data, we inferred the history of population sizes and found that village dog populations experienced a 9-13 fold reduction in historic effective population size relative to wolves.

Genotypic and phenotypic consequences of domestication in dogs

Runs of homozygosity (ROH) are genomic regions that arise when two copies of identical haplotypes are inherited from a shared common ancestor. In this study, we leverage ROH to identify associations between genetic diversity and non-disease phenotypes in Canis lupus familiaris (dogs). We find significant association between the ROH inbreeding coefficient (FROH) and several phenotypic traits. These traits include height, weight, lifespan, muscled, white coloring of the head and chest, furnishings, and fur length. After correcting for population structure, we identified more than 45 genes across the examined quantitative traits that exceed the threshold for suggestive significance. We observe distinct distributions of inbreeding and elevated levels of long ROH in modern breed dogs compared to more ancient breeds, which aligns with breeding practices during Victorian era breed establishment. Our results highlight the impact of non-additive variation and of polygenicity on complex quantitative phenotypes in dogs due to domestication and the breed formation bottleneck.

Evolutionary consequences of domestication on the selective effects of new amino acid changing mutations in canids

The domestication of wild canids led to dogs no longer living in the wild but instead residing alongside humans. Extreme changes in behavior and diet associated with domestication may have led to the relaxation of the selective pressure on traits that may be less important in the domesticated context. Thus, here we hypothesize that strongly deleterious mutations may have become less deleterious in domesticated populations. We test this hypothesis by estimating the distribution of fitness effects (DFE) for new amino acid changing mutations using whole-genome sequence data from 24 gray wolves and 61 breed dogs. We find that the DFE is strikingly similar across canids, with 26-28% of new amino acid changing mutations being neutral/nearly neutral (|s| < 1e-5), and 41-48% under strong purifying selection (|s| > 1e-2). Our results are robust to different model assumptions suggesting that the DFE is stable across short evolutionary timescales, even in the face of putative drastic changes in the selective pressure caused by artificial selection during domestication and breed formation. On par with previous works describing DFE evolution, our data indicate that the DFE of amino acid changing mutations depends more strongly on genome structure and organismal characteristics, and less so on shifting selective pressures or environmental factors. Given the constant DFE and previous data showing that genetic variants that differentiate wolf and dog populations are enriched in regulatory elements, we speculate that domestication may have had a larger impact on regulatory variation than on amino acid changing mutations.

Genome-wide association study for morphological and hunting-behavior traits in Braque Français Type Pyrénées dogs: A preliminary study

High-throughput genotyping offers great potential to increase our understanding of the genomic basis of canid variation. Braque Français Type Pyrénées (BRA) are smart, agile, and friendly dogs originally developed for tracking, hunting, and retrieving feathered game. On a population of 44 unrelated BRA dogs, single nucleotide polymorphism (SNP) genotype data from the CanineHD Whole-Genome Genotyping BeadChip and evaluation scores for 12 traits related to morphology and hunting performance were available. After quality filtering, 95,859 SNPs on the 38 dog autosomes (CFA) were retained. Phenotypic scores were expressed on a scale from 1 (worst) to 6 (best) and were mostly poorly to moderately correlated except for some morphological traits (e.g. r = 0.81 between the conformation of the head and that of the eye). From GWAS, a total of 378 SNP-phenotype associations with posterior odds of association > 1 have been detected. The strongest associations were found for the eye conformation, for the skull/muzzle ratio, and for connection to the hunter. These included both new and previously identified markers and genes potentially involved with type and behavior traits in BRA. Six of the significant markers mapped within SETDB2, a gene known to be related to pointing behavior in dogs. These results advance our understanding of the genetic basis for morphology and hunting behavior in dogs and identify new variants which are potential targets for further research.

Dogs and their genes: what ever will they think of next?

The Edward Novitski Prize recognizes creativity and intellectual ingenuity in the solution of problems in genetics research. The prize honors scientific experimental work-either a single experimental accomplishment or a body of work. Ostrander is recognized for work developing the domestic dog as an experimental system for solving fundamental biological problems and identifying genetic sequences of relevance to human health and disease. Including work on disease and behavioral health, Ostrander has shown a dedication to creative methods for understanding canine genetics and the value of translating research organisms to human genetics.

Historic dog Furs Unravel the Origin and Artificial Selection of Modern Nordic Lapphund and Elkhound dog Breeds

The origins and extreme morphological evolution of the modern dog breeds are poorly studied because the founder populations are extinct. Here, we analyse eight 100 to 200 years old dog fur samples obtained from traditional North Swedish clothing, to explore the origin and artificial selection of the modern Nordic Lapphund and Elkhound dog breeds. Population genomic analysis confirmed the Lapphund and Elkhound breeds to originate from the local dog population, and showed a distinct decrease in genetic diversity in agreement with intense breeding. We identified eleven genes under positive selection during the breed development. In particular, the MSRB3 gene, associated with breed-related ear morphology, was selected in all Lapphund and Elkhound breeds, and functional assays showed that a SNP mutation in the 3'UTR region suppresses its expression through miRNA regulation. Our findings demonstrate analysis of near-modern dog artifacts as an effective tool for interpreting the origin and artificial selection of the modern dog breeds.

Improving laboratory animal genetic reporting: LAG-R guidelines

The biomedical research community addresses reproducibility challenges in animal studies through standardized nomenclature, improved experimental design, transparent reporting, data sharing, and centralized repositories. The ARRIVE guidelines outline documentation standards for laboratory animals in experiments, but genetic information is often incomplete. To remedy this, we propose the Laboratory Animal Genetic Reporting (LAG-R) framework. LAG-R aims to document animals' genetic makeup in scientific publications, providing essential details for replication and appropriate model use. While verifying complete genetic compositions may be impractical, better reporting and validation efforts enhance reliability of research. LAG-R standardization will bolster reproducibility, peer review, and overall scientific rigor.

Genomic analyses of withers height and linear conformation traits in German Warmblood horses using imputed sequence-level genotypes

BACKGROUND

Body conformation, including withers height, is a major selection criterion in horse breeding and is associated with other important traits, such as health and performance. However, little is known about the genomic background of equine conformation. Therefore, the aim of this study was to use imputed sequence-level genotypes from up to 4891 German Warmblood horses to identify genomic regions associated with withers height and linear conformation traits. Furthermore, the traits were genetically characterised and putative causal variants for withers height were detected.

RESULTS

A genome-wide association study (GWAS) for withers height confirmed the presence of a previously known quantitative trait locus (QTL) on Equus caballus (ECA) chromosome 3 close to the LCORL/NCAPG locus, which explained 16% of the phenotypic variance for withers height. An additional significant association signal was detected on ECA1. Further investigations of the region on ECA3 identified a few promising candidate causal variants for withers height, including a nonsense mutation in the coding sequence of the LCORL gene. The estimated heritability for withers height was 0.53 and ranged from 0 to 0.34 for the conformation traits. GWAS identified significantly associated variants for more than half of the investigated conformation traits, among which 13 showed a peak on ECA3 in the same region as withers height. Genetic parameter estimation revealed high genetic correlations between these traits and withers height for the QTL on ECA3.

CONCLUSIONS

The use of imputed sequence-level genotypes from a large study cohort led to the discovery of novel QTL associated with conformation traits in German Warmblood horses. The results indicate the high relevance of the QTL on ECA3 for various conformation traits, including withers height, and contribute to deciphering causal mutations for body size in horses.

Identification of Copy Number Variations and Selection Signatures in Wannan Spotted Pigs by Whole Genome Sequencing Data: A Preliminary Study

Copy number variation (CNV) is an important structural variation used to elucidate complex economic traits. In this study, we sequenced 25 Wannan spotted pigs (WSPs) to detect their CNVs and identify their selection signatures compared with those of 10 Asian wild boars. A total of 14,161 CNVs were detected in the WSPs, accounting for 0.72% of the porcine genome. The fixation index (Fst) was used to identify the selection signatures, and 195 CNVs with the top 1% of the Fst value were selected. Eighty genes were identified in the selected CNV regions. Functional GO and KEGG analyses revealed that the genes within these selected CNVs are associated with key traits such as reproduction (GAL3ST1 and SETD2), fatty acid composition (PRKG1, ACACA, ACSL3, UGT8), immune system (LYZ), ear size (WIF1), and feed efficiency (VIPR2). The findings of this study contribute novel insights into the genetic CNVs underlying WSP characteristics and provide essential information for the protection and utilization of WSP populations.

Increased susceptibility to Mycobacterium avium complex infection in miniature Schnauzer dogs caused by a codon deletion in CARD9

Mammals are generally resistant to Mycobacterium avium complex (MAC) infections. We report here on a primary immunodeficiency disorder causing increased susceptibility to MAC infections in a canine breed. Adult Miniature Schnauzers developing progressive systemic MAC infections were related to a common founder, and pedigree analysis was consistent with an autosomal recessive trait. A genome-wide association study and homozygosity mapping using 8 infected, 9 non-infected relatives, and 160 control Miniature Schnauzers detected an associated region on chromosome 9. Whole genome sequencing of 2 MAC-infected dogs identified a codon deletion in the CARD9 gene (c.493_495del; p.Lys165del). Genotyping of Miniature Schnauzers revealed the presence of this mutant CARD9 allele worldwide, and all tested MAC-infected dogs were homozygous mutants. Peripheral blood mononuclear cells from a dog homozygous for the CARD9 variant exhibited a dysfunctional CARD9 protein with impaired TNF-α production upon stimulation with the fungal polysaccharide β-glucan that activates the CARD9-coupled C-type lectin receptor, Dectin-1. While CARD9-deficient knockout mice are susceptible to experimental challenges by fungi and mycobacteria, Miniature Schnauzer dogs with systemic MAC susceptibility represent the first spontaneous animal model of CARD9 deficiency, which will help to further elucidate host defense mechanisms against mycobacteria and fungi and assess potential therapies for animals and humans.

Defining a Haplotype Encompassing the LCORL-NCAPG Locus Associated with Increased Lean Growth in Beef Cattle

Numerous studies have shown genetic variation at the LCORL-NCAPG locus is strongly associated with growth traits in beef cattle. However, a causative molecular variant has yet to be identified. To define all possible candidate variants, 34 Charolais-sired calves were whole-genome sequenced, including 17 homozygous for a long-range haplotype associated with increased growth (QQ) and 17 homozygous for potential ancestral haplotypes for this region (qq). The Q haplotype was refined to an 814 kb region between chr6:37,199,897-38,014,080 and contained 218 variants not found in qq individuals. These variants include an insertion in an intron of NCAPG, a previously documented mutation in NCAPG (rs109570900), two coding sequence mutations in LCORL (rs109696064 and rs384548488), and 15 variants located within ATAC peaks that were predicted to affect transcription factor binding. Notably, rs384548488 is a frameshift variant likely resulting in loss of function for long isoforms of LCORL. To test the association of the coding sequence variants of LCORL with phenotype, 405 cattle from five populations were genotyped. The two variants were in complete linkage disequilibrium. Statistical analysis of the three populations that contained QQ animals revealed significant (p < 0.05) associations with genotype and birth weight, live weight, carcass weight, hip height, and average daily gain. These findings affirm the link between this locus and growth in beef cattle and describe DNA variants that define the haplotype. However, further studies will be required to define the true causative mutation.

Transposons in the Williams-Beuren Syndrome Critical Region are Associated with Social Behavior in Assistance Dogs

A strong signature of selection in the domestic dog genome is found in a five-megabase region of chromosome six in which four structural variants derived from transposons have previously been associated with human-oriented social behavior, such as attentional bias to social stimuli and social interest in strangers. To explore these genetic associations in more phenotypic detail-as well as their role in training success in a specialized assistance dog program-we genotyped 1001 assistance dogs from Canine Companions for Independence®, including both successful graduates and dogs released from the training program for behaviors incompatible with their working role. We collected phenotypes on each dog using puppy-raiser questionnaires, trainer questionnaires, and both cognitive and behavioral tests. Using Bayesian mixed models, we found strong associations (95% credibility intervals excluding zero) between genotypes and certain behavioral measures, including separation-related problems, aggression when challenged or corrected, and reactivity to other dogs. Furthermore, we found moderate differences in the genotypes of dogs who graduated versus those who did not; insertions in GTF2I showed the strongest association with training success (β = 0.23, CI95% = - 0.04, 0.49), translating to an odds-ratio of 1.25 for one insertion. Our results provide insight into the role of each of these four transposons in canine sociability and may inform breeding and training practices for working dog organizations. Furthermore, the observed importance of the gene GTF2I supports the emerging consensus that variation in GTF2I genotypes and expression have important consequences for social behavior broadly.

Whole-Genome Scanning for Selection Signatures Reveals Candidate Genes Associated with Growth and Tail Length in Sheep

Compared to Chinese indigenous sheep, Western sheep have rapid growth rate, larger physique, and higher meat yield. These excellent Western sheep were introduced into China for crossbreeding to expedite the enhancement of production performance and mutton quality in local breeds. Here, we investigated population genetic structure and genome-wide selection signatures among the Chinese indigenous sheep and the introduced sheep based on whole-genome resequencing data. The PCA, N-J tree and ADMIXTURE results showed significant genetic difference between Chinese indigenous sheep and introduced sheep. The nucleotide diversity (π) and linkage disequilibrium (LD) decay results indicated that the genomic diversity of introduced breeds were lower. Then, Fst & π ratio, XP-EHH, and de-correlated composite of multiple signals (DCMS) methods were used to detect the selection signals. The results showed that we identified important candidate genes related to growth rate and body size in the introduced breeds. Selected genes with stronger selection signatures are associated with growth rate (CRADD), embryonic development (BVES, LIN28B, and WNT11), body size (HMGA2, MSRB3, and PTCH1), muscle development and fat metabolism (MSTN, PDE3A, LGALS12, GGPS1, and SAR1B), wool color (ASIP), and hair development (KRT71, KRT74, and IRF2BP2). Thus, these genes have the potential to serve as candidate genes for enhancing the growth traits of Chinese indigenous sheep. We also identified tail-length trait-related candidate genes (HOXB13, LIN28A, PAX3, and VEGFA) in Chinese long-tailed breeds. Among these genes, HOXB13 is the main candidate gene for sheep tail length phenotype. LIN28A, PAX3, and VEGFA are related to embryonic development and angiogenesis, so these genes may be candidate genes for sheep tail type traits. This study will serve as a foundation for further genetic improvement of Chinese indigenous sheep and as a reference for studies related to growth and development of sheep.

Genome-wide association study of growth traits and validation of key mutations (MSTN c.C861T) associated with the muscle mass of meat pigeons

Selective breeding of meat pigeons is primarily based on growth traits, especially muscle mass (MM). Identification of functional genes and molecular markers of growth and slaughter traits through a genome-wide association study (GWAS) will help to elucidate the underlying molecular mechanisms and provide a theoretical basis for the selective breeding of meat pigeons. The phenotypic data of body weight (BW) and body size (BS) of 556 meat pigeons at 52 and 80 weeks of age were collected. In total, 160 434 high-quality single nucleotide polymorphism sites were obtained by restriction site-associated DNA sequencing. The GWAS analysis revealed that MSTN, IGF2BP3 and NCAPG/LCORL were important candidate genes affecting the growth traits of meat pigeons. IGF2BP3 and NCAPG/LCORL were highly correlated to BW and BS, which are related to overall growth and development, while MSTN was associated with pectoral thickness and BW. Phenotypic association validation with the use of two meat pigeon populations found that the MSTN mutation c.C861T determines the MM. These results provide new insights into the genetic mechanisms underlying phenotypic variations of growth traits and MM in commercial meat pigeons. The identified markers and genes provide a theoretical basis for the selective breeding of meat pigeons.

A purebred South American breed showing high effective population size and independent breed ancestry: The Chilean Terrier

The Chilean Terrier is a known breed in Chile that has not been genetically assessed despite its distinctive color patterns, agility, and hardiness across the diversity of climates encountered within the Chilean landscape. The population structure and its relatedness with other breeds, as well as the actual origin of the breed, remain unknown. We estimated several population parameters using samples from individuals representing the distribution of the Chilean Terrier across the country. By utilizing the Illumina HD canine genotyping array, we computed the effective population size (Ne ), individual inbreeding, and relatedness to evaluate the genetic diversity of the breed. The results show that linkage disequilibrium was relatively low and decayed rapidly; in fact, Ne was very high when compared to other breeds, and similar to other American indigenous breeds (such as the Chihuahua with values of Ne near 500). These results are in line with the low estimates of genomic inbreeding and relatedness and the relatively large number of effective chromosome segments (Me  = 2467) obtained using the properties of the genomic relationship matrix. Between population analysis (cross-population extended haplotype homozygosity, di ) with other breeds such as the Jack Russell Terrier, the Peruvian-Inca Orchid, and the Chihuahua suggested that candidate regions harboring FGF5, PAX3, and ASIP, probably explained some morphological traits, such as the distinctive color pattern characteristic of the breed. When considering Admixture estimates and phylogenetic analysis, together with other breeds of American and European origin, the Chilean Terrier does not have a recent European ancestry. Overall, the results suggest that the breed has evolved independently in Chile from other terrier breeds, from an unknown European terrier ancestor.

Sequenced-based GWAS for linear classification traits in Belgian Blue beef cattle reveals new coding variants in genes regulating body size in mammals

BACKGROUND

Cohorts of individuals that have been genotyped and phenotyped for genomic selection programs offer the opportunity to better understand genetic variation associated with complex traits. Here, we performed an association study for traits related to body size and muscular development in intensively selected beef cattle. We leveraged multiple trait information to refine and interpret the significant associations.

RESULTS

After a multiple-step genotype imputation to the sequence-level for 14,762 Belgian Blue beef (BBB) cows, we performed a genome-wide association study (GWAS) for 11 traits related to muscular development and body size. The 37 identified genome-wide significant quantitative trait loci (QTL) could be condensed in 11 unique QTL regions based on their position. Evidence for pleiotropic effects was found in most of these regions (e.g., correlated association signals, overlap between credible sets (CS) of candidate variants). Thus, we applied a multiple-trait approach to combine information from different traits to refine the CS. In several QTL regions, we identified strong candidate genes known to be related to growth and height in other species such as LCORL-NCAPG or CCND2. For some of these genes, relevant candidate variants were identified in the CS, including three new missense variants in EZH2, PAPPA2 and ADAM12, possibly two additional coding variants in LCORL, and candidate regulatory variants linked to CCND2 and ARMC12. Strikingly, four other QTL regions associated with dimension or muscular development traits were related to five (recessive) deleterious coding variants previously identified.

CONCLUSIONS

Our study further supports that a set of common genes controls body size across mammalian species. In particular, we added new genes to the list of those associated with height in both humans and cattle. We also identified new strong candidate causal variants in some of these genes, strengthening the evidence of their causality. Several breed-specific recessive deleterious variants were identified in our QTL regions, probably as a result of the extreme selection for muscular development in BBB cattle.

Significant association of mutations close to LCORL gene with growth performance in Zhedong white geese (Anser cygnoides)

This study aimed to investigate the role of the LCORL gene in regulating the growth performance of Zhedong white (ZDW) geese, belonging to the swan geese (Anser cygnoides), and identify possible selective signatures in diverse goose breeds. Single nucleotide polymorphisms around LCORL were genotyped, and their associations with body-size-related (BSR) traits were estimated. The results showed that the genotyped loci upstream of LCORL were significantly related to the body weight and breast width of ZDW geese aged 10 weeks (p < 0.05). A genome scan comparing expected heterozygosity among different breeds identified a ~150 kb long genomic region with extremely low heterozygosity downstream of LCORL among swan geese. Further, significant associations of variants within the low heterozygosity region among ZDW geese with BSR traits, including body weight, body length and breast width (p < 0.05) were also detected. Overall, mutations adjacent to LCORL were related to the growth performance of swan geese, and the significant effects of variants in a low-heterozygosity region on BSR traits provided valuable insights into the molecular mechanism of artificial selection reshaping body stature in swan geese.

Lymphoma in Border Collies: Genome-Wide Association and Pedigree Analysis

There has been considerable interest in studying cancer in dogs and its potential as a model system for humans. One area of research has been the search for genetic risk variants in canine lymphoma, which is amongst the most common canine cancers. Previous studies have focused on a limited number of breeds, but none have included Border Collies. The aims of this study were to identify relationships between Border Collie lymphoma cases through an extensive pedigree investigation and to utilise relationship information to conduct genome-wide association study (GWAS) analyses to identify risk regions associated with lymphoma. The expanded pedigree analysis included 83,000 Border Collies, with 71 identified lymphoma cases. The analysis identified affected close relatives, and a common ancestor was identified for 54 cases. For the genomic study, a GWAS was designed to incorporate lymphoma cases, putative "carriers", and controls. A case-control GWAS was also conducted as a comparison. Both analyses showed significant SNPs in regions on chromosomes 18 and 27. Putative top candidate genes from these regions included DLA-79, WNT10B, LMBR1L, KMT2D, and CCNT1.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture

BACKGROUND

The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function.

RESULTS

We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection.

CONCLUSIONS

We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

The Shepherd and the Hunter: A Genomic Comparison of Italian Dog Breeds

Shepherd and hunting dogs have undergone divergent selection for specific tasks, resulting in distinct phenotypic and behavioural differences. Italy is home to numerous recognized and unrecognized breeds of both types, providing an opportunity to compare them genomically. In this study, we analysed SNP data obtained from the CanineHD BeadChip, encompassing 116 hunting dogs (representing 6 breeds) and 158 shepherd dogs (representing 9 breeds). We explored the population structure, genomic background, and phylogenetic relationships among the breeds. To compare the two groups, we employed three complementary methods for selection signature detection: FST, XP-EHH, and ROH. Our results reveal a clear differentiation between shepherd and hunting dogs as well as between gun dogs vs. hounds and guardian vs. herding shepherd dogs. The genomic regions distinguishing these groups harbour several genes associated with domestication and behavioural traits, including gregariousness (WBSRC17) and aggressiveness (CDH12 and HTT). Additionally, genes related to morphology, such as size and coat colour (ASIP and TYRP1) and texture (RSPO2), were identified. This comparative genomic analysis sheds light on the genetic underpinnings of the phenotypic and behavioural variations observed in Italian hunting and shepherd dogs.

Homozygosity mapping in the Kazakh national dog breed Tazy

The Tazy is a breed of sighthound common in Kazakhstan. The identification of runs of homozygosity (ROH) is an informative approach to assessing the history and possible patterns of directional selection pressure. To our knowledge, the present study is the first to provide an overview of the ROH pattern in the Tazy dogs from a genome-wide perspective. The ROH of the Tazy was found to be mainly composed of shorter segments (1-2 Mb), accounting for approximately 67% of the total ROH. The estimated ROH-based inbreeding coefficients (F_ROH) ranged from 0.028 to 0.058 with a mean of 0.057. Five genomic regions under positive selection were identified on chromosomes 18, 22, and 25. The regions on chromosomes 18 and 22 may be breed specific, while the region on chromosome 22 overlaps with regions of hunting traits in other hunting dog breeds. Among the 12 candidate genes located in these regions, the gene CAB39L may be a candidate that affects running speed and endurance of the Tazy dog. Eight genes could belong to an evolutionarily conserved complex as they were clustered in a large protein network with strong linkages. The results may enable effective interventions when incorporated into conservation planning and selection of the Tazy breed.

Systematically identifying genetic signatures including novel SNP-clusters, nonsense variants, frame-shift INDELs, and long STR expansions that potentially link to unknown phenotypes existing in dog breeds

BACKGROUND

In light of previous studies that profiled breed-specific traits or used genome-wide association studies to refine loci associated with characteristic morphological features in dogs, the field has gained tremendous genetic insights for known dog traits observed among breeds. Here we aim to address the question from a reserve perspective: whether there are breed-specific genotypes that may underlie currently unknown phenotypes. This study provides a complete set of breed-specific genetic signatures (BSGS). Several novel BSGS with significant protein-altering effects were highlighted and validated.

RESULTS

Using the next generation whole-genome sequencing technology coupled with unsupervised machine learning for pattern recognitions, we constructed and analyzed a high-resolution sequence map for 76 breeds of 412 dogs. Genomic structures including novel single nucleotide polymorphisms (SNPs), SNP clusters, insertions, deletions (INDELs) and short tandem repeats (STRs) were uncovered mutually exclusively among breeds. We also partially validated some novel nonsense variants by Sanger sequencing with additional dogs. Four novel nonsense BSGS were found in the Bernese Mountain Dog, Samoyed, Bull Terrier, and Basset Hound, respectively. Four INDELs resulting in either frame-shift or codon disruptions were found in the Norwich Terrier, Airedale Terrier, Chow Chow and Bernese Mountain Dog, respectively. A total of 15 genomic regions containing three types of BSGS (SNP-clusters, INDELs and STRs) were identified in the Akita, Alaskan Malamute, Chow Chow, Field Spaniel, Keeshond, Shetland Sheepdog and Sussex Spaniel, in which Keeshond and Sussex Spaniel each carried one amino-acid changing BSGS in such regions.

CONCLUSION

Given the strong relationship between human and dog breed-specific traits, this study might be of considerable interest to researchers and all. Novel genetic signatures that can differentiate dog breeds were uncovered. Several functional genetic signatures might indicate potentially breed-specific unknown phenotypic traits or disease predispositions. These results open the door for further investigations. Importantly, the computational tools we developed can be applied to any dog breeds as well as other species. This study will stimulate new thinking, as the results of breed-specific genetic signatures may offer an overarching relevance of the animal models to human health and disease.

Predicting health outcomes in dogs using insurance claims data

In this paper we propose a machine learning-based approach to predict a multitude of insurance claim categories related to canine diseases. We introduce several machine learning approaches that are evaluated on a pet insurance dataset consisting of 785,565 dogs from the US and Canada whose insurance claims have been recorded over 17 years. 270,203 dogs with a long insurance tenure were used to train a model while the inference is applicable to all dogs in the dataset. Through this analysis we demonstrate that with this richness of data, supported by the right feature engineering, and machine learning approaches, 45 disease categories can be predicted with high accuracy.

Psychological Stress Is Associated with Increased Cancer Risk in Dogs

Although there is evidence that psychological stress may be associated with increased cancer risk, the effect of stress on cancer risk is difficult to study, both in humans, due to socioeconomic factors, and in animal models, due to questionable biological relevance. Here, we test whether heritable canine temperament that increases psychological stress is associated with cancer risk. The study data are breed-specific averages of incidences of multiple cancer types and of temperament classes. The latter are derived from a latent class analysis of behavioral questionnaires completed by owners (C-BARQ). We thus classified the dogs according to whether they are calm vs. reactive within and across breeds. Using meta-analysis approaches, we modeled the risk of multiple cancer types in calm vs. reactive dogs. We adjusted for breed averages of body mass and lifespan, which are common confounders that impact cancer. Our study confirms that body size has a significant effect of on risk of multiple types of cancers in dogs and shows for the first time that temperament also has a moderate effect. These findings suggest dog models of heritable psychological stress are suitable for molecular epidemiological and translational studies on its effects on cancer risk.

Genome-wide association study reveals that the IBSP locus affects ear size in cattle

Ear size is a classical model for hot climate adaptation following the evolution, but the genetic basis of the traits associated with ear size remains to be elucidated. Here, we performed a genome-wide association study on 158 cattle to explain the genetic mechanism of ear size. One region on BTA6 between 36.79 and 38.80 Mb included 50 suggestive SNPs and 4 significant SNPs that were significantly associated with ear size. The most significant locus (P = 1.30 × 10-8) was a missense mutation (T250I) on the seventh exon of integrin-binding sialoprotein (IBSP), which had an allele substitution effect of 23.46 cm2 for ear size. Furthermore, this mutation will cause changes in the three-dimensional structure of the protein. To further identify genes underlying this typical feature, we performed a genome scan among nine cattle breeds with different ear sizes by using SweeD. Results suggested that IBSP was under positive selection among four breeds with relatively large ear sizes. The expression levels of IBSP in ear tissues of large- and small-ear cattle were significantly different. A haplotype diversity survey of this missense mutation in worldwide cattle breeds strongly implied that the origin of this missense mutation event was Bos taurus. These findings have important theoretical importance for the exploration of major genes associated with ear size and provide important molecular markers for the identification of cattle germplasm resources.

Social determinants of health and disease in companion dogs: a cohort study from the Dog Aging Project

Exposure to social environmental adversity is associated with health and survival across many social species, including humans. However, little is known about how these health and mortality effects vary across the lifespan and may be differentially impacted by various components of the environment. Here, we leveraged a relatively new and powerful model for human aging, the companion dog, to investigate which components of the social environment are associated with dog health and how these associations vary across the lifespan. We drew on comprehensive survey data collected on 21,410 dogs from the Dog Aging Project and identified five factors that together explained 33.7% of the variation in a dog's social environment. Factors capturing financial and household adversity were associated with poorer health and lower physical mobility in companion dogs, while factors that captured social support, such as living with other dogs, were associated with better health when controlling for dog age and weight. Notably, the effects of each environmental component were not equal: the effect of social support was 5× stronger than financial factors. The strength of these associations depended on the age of the dog, including a stronger relationship between the owner's age and the dog's health in younger as compared to older dogs. Taken together, these findings suggest the importance of income, stability and owner's age on owner-reported health outcomes in companion dogs and point to potential behavioral and/or environmental modifiers that can be used to promote healthy aging across species.

Genetic dissection of behavioral traits related to successful training of drug detection dogs

Drug detection dogs play integral roles in society. However, the interplay between their behaviors and genetic characteristics underlying their performance remains uninvestigated. Herein, more than 120,000 genetic variants were evaluated in 326 German Shepherd or Labrador Retriever dogs to profile the genetic traits associated with various behavioral traits related to the successful training of drug detection dogs. Behavioral breed differences were observed in 'friendliness to humans' and 'tolerance to dogs.' A genome-wide association study within both breeds identified 11 regions potentially associated with drug detection dog characteristics as well as 'interest in the target' and 'friendliness to humans,' which are related to drug detection abilities. Among them, 63 protein-coding genes, including Atat1 and Pfn2 known to be associated with anxiety-related or exploration behavior in mice, respectively, were located surrounding the identified candidate polymorphisms. This study highlights genetic characteristics associated with behavioral traits that are important for the successful training of drug detection dogs. Thus, these findings may facilitate improved breeding and training of these dogs.

Identification of Genomic Signatures in Bullmastiff Dogs Using Composite Selection Signals Analysis of 23 Purebred Clades

Dog breeds represent canine sub-populations with distinctive phenotypic features and limited genetic diversity. We have established a resource to study breed-specific genetic diversity. Utilising genetic resources within our laboratory biobank, public domain genotype data and the phylogenetic framework of 23 breed clades, the primary objective for this study was to identify genomic regions that differentiate the Bullmastiff breed. Through application of a composite index analysis (CSS), genomic signatures were identified in Bullmastiffs when compared to the formative breeds, Mastiffs and Bulldogs, and to 22 other breed groups. Significant regions were identified on 15 chromosomes, with the most differentiated regions found on CFA1, CFA9, and CFA18. These regions may reflect genetic drift following establishment of the breed or the effects of selective breeding during development of the modern Bullmastiff. This was supported by analysis of genes from the identified genomic regions, including 458 genes from the multi-clade analysis, which revealed enriched pathways that may be related to characteristic traits and distinct morphology of the breed. The study demonstrates the utility of the CSS method in breed-specific genome analysis and advances our understanding of genetic diversity in Bullmastiff dogs.

Identification of Genetic Risk Factors for Monogenic and Complex Canine Diseases

Advances in DNA sequencing and other technologies have greatly facilitated the identification of genetic risk factors for inherited diseases in dogs. We review recent technological developments based on selected examples from canine disease genetics. The identification of disease-causing variants in dogs with monogenic diseases may become a widely employed diagnostic approach in clinical veterinary medicine in the not-too-distant future. Diseases with complex modes of inheritance continue to pose challenges to researchers but have also become much more tangible than in the past. In addition to strategies for identifying genetic risk factors, we provide some thoughts on the interpretation of sequence variants that are largely inspired by developments in human clinical genetics.

Analysis of dog breed diversity using a composite selection index

During breed development, domestic dogs have undergone genetic bottlenecks and sustained selective pressures, as a result distinctive genomic diversity occurs to varying degrees within and between breed groups. This diversity can be identified using standard methods or combinations of these methods. This study explored the application of a combined selection index, composite selection signals (CSS), derived from multiple methods to an existing genotype dataset from three breed groups developed in distinct regions of Asia: Qinghai-Tibet plateau dogs (adapted to living at altitude), Xi dogs (with superior running ability) and Mountain hounds (used for hunting ability). The CSS analysis confirmed top ranked genomic regions on CFA10 and CFA21 in Qinghai-Tibet plateau dogs, CFA1 in Xi dogs and CFA5 in Mountain hounds. CSS analysis identified additional significant genomic regions in each group, defined by a total of 1,397, 1,475 and 1,675 significant SNPs in the Qinghai-Tibetan Plateau dogs, Xi dogs and Mountain hounds, respectively. Chitinase 3 Like 1 (CHI3L1) and Leucine Rich Repeat Containing G Protein-Coupled Receptor 6 (LGR6) genes were located in the top ranked region on CFA7 (0.02-1 Mb) in the Qinghai-Tibetan Plateau dogs. Both genes have been associated with hypoxia responses or altitude adaptation in humans. For the Xi dogs, the top ranked region on CFA25 contained the Transient Receptor Potential Cation Channel Subfamily C Member 4 (TRPC4) gene. This calcium channel is important for optimal muscle performance during exercise. The outstanding signals in the Mountain dogs were on CFA5 with 213 significant SNPs that spanned genes involved in cardiac development, sight and generation of biochemical energy. These findings support the use of the combined index approach for identifying novel regions of genome diversity in dogs. As with other methods, the results do not prove causal links between these regions and phenotypes, but they may assist in focusing future studies that seek to identify functional pathways that contribute to breed diversity.

Exploration of body weight in 115 000 young adult dogs of 72 breeds

High body weight (BW), due to large size or excess body fat, has been associated with developmental and metabolic alterations, and degenerative diseases in dogs. Study objectives were to determine mean BW in young adult dogs of different breeds, including changes over a 10-year period. Body weight data from the official Swedish hip dysplasia screening program were used, including data from dogs screened at 1-2.5 years of age, in breeds with ≥ 15 individual observations/year during 2007-2016. Mean BW per breed and sex was established from 114 568 dogs representing 72 breeds. Estimates of breed BW showed significant change in 33 (45%) breeds over the 10-year period. Body weight increased in five breeds (2-14% change) and decreased in 26 breeds (1-8% change). In two breeds, BW increased in male and decreased in female dogs. This observational study provides extensive breed BW data on young adult dogs. The change in breed BW, noted in almost half of the breeds, could be due to changes either in size or in body fat mass. In certain breeds, the change in BW over time might have an impact on overall health. Studies with simultaneous evaluation of BW and body condition over time are warranted.

A comprehensive review of genomic perspectives of canine diseases as a model to study human disorders

The domestic dog has been given considerable attention as a system for investigating the genetics of human diseases. Population diversity and breed structure are unique features that make dogs particularly amenable to genetic studies. Dogs show distinguished features of breed-specific homogeneity, which is associated with striking interbreed heterogeneity. This review discusses the significance of studying the genetic maps, genome-wide association studies (GWAS), and usefulness of this species as an animal model. Most canine genetic disorders are similar to those of humans, including inherited, psychiatric, and genetic disorders. In addition to revealing new candidate genes, canine models allow access to experimental resources, such as cells, tissues, and even live animals, for research and intervention purposes.

Selection Signatures in Italian Livestock Guardian and Herding Shepherd Dogs

Livestock guardian (LGD) and herding shepherd (HSD) dogs have distinct morphological and behavioural characteristics, long selected by farmers and breeders, to accomplish different tasks. This study aimed to find the genomic regions that best differentiate and characterise Italian LGD and HSD. Genomic data of 158 dogs of four LGD and five HSD breeds, obtained with the 170K canine SNPchip, were collected. The two groups were compared using FST and XP-EHH analyses, identifying regions containing 29 genes. Moreover, 16 islands of runs of homozygosity were found in LGD, and 15 in HSD; 4 of them were partially shared. Among the genes found that better differentiated HSD and LGD, several were associated with dog domestication and behavioural aspects; particularly, MSRB3 and LLPH were linked to herding behaviour in previous studies. Others, DYSK, MAP2K5, and RYR, were related to body size and muscle development. Prick ears prevailed in sampled HSD, and drop ears in LGD; this explains the identification of WIF1 and MSRB3 genes. Unexpectedly, a number of genes were also associated with eye development and functionality. These results shed further light on the differences that human selection introduced in dogs aimed at different duties, even in a limited geographic area such as Italy.

Genomic adaptation of Ethiopian indigenous cattle to high altitude

The mountainous areas of Ethiopia represent one of the most extreme environmental challenges in Africa faced by humans and other inhabitants. Selection for high-altitude adaptation is expected to have imprinted the genomes of livestock living in these areas. Here we assess the genomic signatures of positive selection for high altitude adaptation in three cattle populations from the Ethiopian mountainous areas (Semien, Choke, and Bale mountains) compared to three Ethiopian lowland cattle populations (Afar, Ogaden, and Boran), using whole-genome resequencing and three genome scan approaches for signature of selection (iHS, XP-CLR, and PBS). We identified several candidate selection signature regions and several high-altitude adaptation genes. These include genes such as ITPR2, MB, and ARNT previously reported in the human population inhabiting the Ethiopian highlands. Furthermore, we present evidence of strong selection and high divergence between Ethiopian high- and low-altitude cattle populations at three new candidate genes (CLCA2, SLC26A2, and CBFA2T3), putatively linked to high-altitude adaptation in cattle. Our findings provide possible examples of convergent selection between cattle and humans as well as unique African cattle signature to the challenges of living in the Ethiopian mountainous regions.

Domestic dog lineages reveal genetic drivers of behavioral diversification

Selective breeding of domestic dogs has generated diverse breeds often optimized for performing specialized tasks. Despite the heritability of breed-typical behavioral traits, identification of causal loci has proven challenging due to the complexity of canine population structure. We overcome longstanding difficulties in identifying genetic drivers of canine behavior by developing a framework for understanding relationships between breeds and the behaviors that define them, utilizing genetic data for over 4,000 domestic, semi-feral, and wild canids and behavioral survey data for over 46,000 dogs. We identify ten major canine genetic lineages and their behavioral correlates and show that breed diversification is predominantly driven by non-coding regulatory variation. We determine that lineage-associated genes converge in neurodevelopmental co-expression networks, identifying a sheepdog-associated enrichment for interrelated axon guidance functions. This work presents a scaffold for canine diversification that positions the domestic dog as an unparalleled system for revealing the genetic origins of behavioral diversity.

Across-breed genetic investigation of canine hip dysplasia, elbow dysplasia, and anterior cruciate ligament rupture using whole-genome sequencing

Here, we report the use of genome-wide association study (GWAS) for the analysis of canine whole-genome sequencing (WGS) repository data using breed phenotypes. Single-nucleotide polymorphisms (SNPs) were called from WGS data from 648 dogs that included 119 breeds from the Dog10K Genomes Project. Next, we assigned breed phenotypes for hip dysplasia (Orthopedic Foundation for Animals (OFA) HD, n = 230 dogs from 27 breeds; hospital HD, n = 279 dogs from 38 breeds), elbow dysplasia (ED, n = 230 dogs from 27 breeds), and anterior cruciate ligament rupture (ACL rupture, n = 279 dogs from 38 breeds), the three most important canine spontaneous complex orthopedic diseases. Substantial morbidity is common with these diseases. Previous within- and between-breed GWAS for HD, ED, and ACL rupture using array SNPs have identified disease-associated loci. Individual disease phenotypes are lacking in repository data. There is a critical knowledge gap regarding the optimal approach to undertake categorical GWAS without individual phenotypes. We considered four GWAS approaches: a classical linear mixed model, a haplotype-based model, a binary case-control model, and a weighted least squares model using SNP average allelic frequency. We found that categorical GWAS was able to validate HD candidate loci. Additionally, we discovered novel candidate loci and genes for all three diseases, including FBX025, IL1A, IL1B, COL27A1, SPRED2 (HD), UGDH, FAF1 (ED), TGIF2 (ED & ACL rupture), and IL22, IL26, CSMD1, LDHA, and TNS1 (ACL rupture). Therefore, categorical GWAS of ancestral dog populations may contribute to the understanding of any disease for which breed epidemiological risk data are available, including diseases for which GWAS has not been performed and candidate loci remain elusive.

Population Structure and Selection Signatures Underlying Domestication Inferred from Genome-Wide Copy Number Variations in Chinese Indigenous Pigs

Single nucleotide polymorphism was widely used to perform genetic and evolution research in pigs. However, little is known about the effect of copy number variation (CNV) on characteristics in pigs. This study performed a genome-wide comparison of CNVs between Wannan black pigs (WBP) and Asian wild boars (AWB), using whole genome resequencing data. By using Manta, we detected in total 28,720 CNVs that covered approximately 1.98% of the pig genome length. We identified 288 selected CNVs (top 1%) by performing Fst statistics. Functional enrichment analyses for genes located in selected CNVs were found to be muscle related (NDN, TMOD4, SFRP1, and SMYD3), reproduction related (GJA1, CYP26B1, WNT5A, SRD5A2, PTPN11, SPEF2, and CCNB1), residual feed intake (RFI) related (MAP3K5), and ear size related (WIF1). This study provides essential information on selected CNVs in Wannan black pigs for further research on the genetic basis of the complex phenotypic and provides essential information for direction in the protection and utilization of Wannan black pig.

Balancing at the Borderline of a Breed: A Case Study of the Hungarian Short-Haired Vizsla Dog Breed, Definition of the Breed Profile Using Simple SNP-Based Methods

The aim of this study was to determine the breed boundary of the Hungarian Short-haired Vizsla (HSV) dog breed. Seventy registered purebred HSV dogs were genotyped on approximately 145,000 SNPs. Principal Component Analysis (PCA) and Admixture analysis certified that they belong to the same population. The outer point of the breed demarcation was a single Hungarian Wire-haired Vizsla (HWV) individual, which was the closest animal genetically to the HSV population in the PCA analysis. Three programs were used for the breed assignment calculations, including the widely used GeneClass2.0 software and two additional approaches developed here: the 'PCA-distance' and 'IBS-central' methods. Both new methods calculate a single number that represents how closely a dog fits into the actual reference population. The former approach calculates this number based on the PCA distances from the median of HSV animals. The latter calculates it from identity by state (IBS) data, measuring the distance from a central animal that is the best representative of the breed. Having no mixed-breed dogs with known HSV genome proportion, admixture animals were simulated by using data of HSV and HWV individuals to calibrate the inclusion/exclusion probabilities for the assignment. The numbers generated from these relatively simple calculations can be used by breeders and clubs to keep their populations under genetic supervision.

Comparison of behavioural tendencies between "dangerous dogs" and other domestic dog breeds - Evolutionary context and practical implications

Aggressive behaviour by dogs is a considerable social problem, but the ability to predict which individuals may have increased aggressive tendencies is very limited, restricting the development of efficient preventive measures. There is a common perception that certain breeds are more likely to exhibit aggressive behaviour, which has contributed to the introduction of breed-specific legislation. The rationale for such legislation explicitly assumes high heritability of this trait while also implying relatively little variation within breeds; these assumptions are largely untested. We compared behavioural tendencies between 8 breeds that are subject to legislation in at least one country and 17 breeds that are not subject to legislation using two validated psychometric tools: the Dog Impulsivity Assessment Scale (DIAS), which scores elements of impulsivity, including a tendency for aggressive behaviour, and Positive and Negative Activation Scale (PANAS), which scores sensitivity to positive and negative stimuli (which may trigger aggressive responses). We found that the two groups of breeds do not differ significantly in the specific DIAS factor relating to aggressive behaviour, "Aggression Threshold and Response to Novelty", or any other DIAS and PANAS factors. We found large variations in all behavioural tendencies measured by both psychometric scales within both groups and within each breed studied. Taken together, our findings indicate that breed alone is not a reliable predictor of individual behavioural tendencies, including those related to aggression, and therefore breed-specific legislation is unlikely to be an effective instrument for reducing risk.

Current Analytical Methods and Research Trends Are Used to Identify Domestic Pig and Wild Boar DNA in Meat and Meat Products

The pig, one of the most important livestock species, is a meaningful source of global meat production. It is necessary, however, to prove whether a food product that a discerning customer selects in a store is actually made from pork or venison, or does not contain it at all. The problem of food authenticity is widespread worldwide, and cases of meat adulteration have accelerated the development of food and the identification methods of feed species. It is worth noting that several different molecular biology techniques can identify a porcine component. However, the precise differentiation between wild boar and a domestic pig in meat products is still challenging. This paper presents the current state of knowledge concerning the species identification of the domestic pig and wild boar DNA in meat and its products.