Dog colour patterns explained by modular promoters of ancient canid origin

Taurine pangenome uncovers a segmental duplication upstream of KIT associated with depigmentation in white-headed cattle

Cattle have been selectively bred for coat color, spotting, and depigmentation patterns. The assumed autosomal dominant inherited genetic variants underlying the characteristic white head of Fleckvieh, Simmental, and Hereford cattle have not been identified yet, although the contribution of structural variation upstream of the KIT gene has been proposed. Here, we construct a graph pangenome from 24 haplotype assemblies representing seven taurine cattle breeds to identify and characterize the white-head-associated locus for the first time based on long-read sequencing data and pangenome analyses. We introduce a pangenome-wide association mapping approach that examines assembly path similarities within the graph to reveal an association between two most likely serial alleles of a complex structural variant (SV) 66 kb upstream of KIT and facial depigmentation. The complex SV contains a variable number of tandemly duplicated 14.3 kb repeats, consisting of LTRs, LINEs, and other repetitive elements, leading to misleading alignments of short and long reads when using a linear reference. We align 250 short-read sequencing samples spanning 15 cattle breeds to the pangenome graph, further validating that the alleles of the SV segregate with head depigmentation. We estimate an increased count of repeats in Hereford relative to Simmental and other white-headed cattle breeds from the graph alignment coverage, suggesting a large under-assembly in the current Hereford-based cattle reference genome, which had fewer copies. Our work shows that exploiting assembly path similarities within graph pangenomes can reveal trait-associated complex SVs.

Two novel SNP variants at ASIP and SNAI2 genes are associated with yellow coat color in rabbits

Rabbits display a wide range of coat colors, with yellow being a particular phenotype that aids in exploring the molecular mechanisms of coat pigmentation. The Fujian yellow (FJY) rabbit, as China's only indigenous breed with a yellow coat, serves as a valuable genetic resource. Fujian yellow rabbits have predominantly yellow fur, with a diluted white hue on the distal limbs and tail. However, the genetic mechanism underlying yellow coat color remains unclear. To address this, we conducted selection signature analysis to identify candidate genes and potential casual mutations underlying the yellow phenotype in rabbits. Utilizing whole-genome resequencing, a total of 22 486 177 high-quality SNPs were identified from 30 individuals belonging to three Chinese indigenous rabbit breeds featured with yellow or non-yellow phenotype. The results revealed that the ASIP gene on chromosome 4 and the SNAI2 gene on chromosome 3 were under strong selection pressure, both of which play pivotal roles in determining coat color phenotypes. The ASIP gene is involved in melanogenesis across various livestock species, while the SNAI2 gene is linked to hypopigmentation in the distal regions such as the limbs and tail. We further identified two SNP variants, g.23870943C>T in the fourth intron of the ASIP gene, which is closely associated with the yellow phenotype, and g.73725380A>G downstream of the SNAI2 gene, probably contributing to the white shading in Fujian yellow rabbits' limb and tail regions. These variants are key determinants in the development of the yellow coat color in rabbits. These findings advance the understanding of coat color pigmentation in domestic animals.

Selection Signatures in Italian Goat Populations Sharing the "facciuto" Phenotype

Background: The presence of light-pigmented facial stripes, parallel on both sides of the cranial region, is a widespread characteristic in various goat breeds of European origin and beyond. In Italy, this phenotype is relatively common from the north to the south of the peninsula. The availability of genotypic data at single-nucleotide polymorphism (SNP) loci for breeds and populations characterized by such a pigmentation pattern enabled us to study the genomic regions potentially correlated with this phenotype, for simplicity referred to as "facciuto". Methods: We adopted an FST-outlier approach to detect signals of differential selection in 18 pairwise comparisons, each involving 6 genetic goat types with the "facciuto" phenotype (Facciuta Lucana, Facciuta della Valnerina, Valfortorina, Teramana, Capestrina, and Roccaverano) contrasted with each of 3 "non-facciuto" goat breeds selected as reference populations (Red Mediterranean, light brown; Saanen, white; Malagueña, mahogany solid). Results: The analysis of the region ±200 kbps upstream and downstream of the two significant signals on chromosome 13 and 15 allowed us to identify, among the annotated genes, ASIP, AHCY, ITCH, DYNLRB1, MAP1LC3A, PIGU, LOC102177263, and DTX4, whose functions could be related to several mechanisms underlying the phenotype under investigation. Conclusions: This study confirmed the fundamental role of ASIP in pigmentation, although additional pathways may concurrently contribute to the determinism of the considered "facciuto" phenotype in Italian goats.

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.

An Evolutionary Perspective on Dog Behavioral Genetics

Dogs have played an outsized role in the field of behavioral genetics since its earliest days. Their unique evolutionary history and ubiquity in the modern world make them a potentially powerful model system for discovering how genetic changes lead to changes in behavior. Genomic technology has supercharged this potential by enabling scientists to sequence the DNA of thousands of dogs and test for correlations with behavioral traits. However, fractures in the early history of animal behavior between biological and psychological subfields may be impeding progress. In addition, canine behavioral genetics has included almost exclusively dogs from modern breeds, who represent just a small fraction of all dog diversity. By expanding the scope of dog behavior studies, and incorporating an evolutionary perspective on canine behavioral genetics, we can move beyond associations to understanding the complex interactions between genes and environment that lead to dog behavior.

The DoGA consortium expression atlas of promoters and genes in 100 canine tissues

The dog, Canis lupus familiaris, is an important model for studying human diseases. Unlike many model organisms, the dog genome has a comparatively poor functional annotation, which hampers gene discovery for development, morphology, disease, and behavior. To fill this gap, we established a comprehensive tissue biobank for both the dog and wolf samples. The biobank consists of 5485 samples representing 132 tissues from 13 dogs, 12 dog embryos, and 24 wolves. In a subset of 100 tissues from nine dogs and 12 embryos, we characterized gene expression activity for each promoter, including alternative and novel, i.e., previously not annotated, promoter regions, using the 5' targeting RNA sequencing technology STRT2-seq. We identified over 100,000 promoter region candidates in the recent canine genome assembly, CanFam4, including over 45,000 highly reproducible sites with gene expression and respective tissue enrichment levels. We provide a promoter and gene expression atlas with interactive, open data resources, including a data coordination center and genome browser track hubs. We demonstrated the applicability of Dog Genome Annotation (DoGA) data and resources using multiple examples spanning canine embryonic development, morphology and behavior, and diseases across species.

Comparative Population Genomics of Arctic Sled Dogs Reveals a Deep and Complex History

Recent evidence demonstrates genomic and morphological continuity in the Arctic ancestral lineage of dogs. Here, we use the Siberian Husky to investigate the genomic legacy of the northeast Eurasian Arctic lineage and model the deep population history using genome-wide single nucleotide polymorphisms. Utilizing ancient dog-calibrated molecular clocks, we found that at least two distinct lineages of Arctic dogs existed in ancient Eurasia at the end of the Pleistocene. This pushes back the origin of sled dogs in the northeast Siberian Arctic with humans likely intentionally selecting dogs to perform different functions and keeping breeding populations that overlap in time and space relatively reproductively isolated. In modern Siberian Huskies, we found significant population structure based on how they are used by humans, recent European breed introgression in about half of the dogs that participate in races, moderate levels of inbreeding, and fewer potentially harmful variants in populations under strong selection for form and function (show, sled show, and racing populations of Siberian Huskies). As the struggle to preserve unique evolutionary lineages while maintaining genetic health intensifies across pedigreed dogs, understanding the genomic history to guide policies and best practices for breed management is crucial to sustain these ancient lineages and their unique evolutionary identity.

A sequence of SVA retrotransposon insertions in ASIP shaped human pigmentation

Retrotransposons comprise about 45% of the human genome1, but their contributions to human trait variation and evolution are only beginning to be explored2,3. Here, we find that a sequence of SVA retrotransposon insertions in an early intron of the ASIP (agouti signaling protein) gene has probably shaped human pigmentation several times. In the UK Biobank (n = 169,641), a recent 3.3-kb SVA insertion polymorphism associated strongly with lighter skin pigmentation (0.22 [0.21-0.23] s.d.; P = 2.8 × 10-351) and increased skin cancer risk (odds ratio = 1.23 [1.18-1.27]; P = 1.3 × 10-28), appearing to underlie one of the strongest common genetic influences on these phenotypes within European populations4-6. ASIP expression in skin displayed the same association pattern, with the SVA insertion allele exhibiting 2.2-fold (1.9-2.6) increased expression. This effect had an unusual apparent mechanism: an earlier, nonpolymorphic, human-specific SVA retrotransposon 3.9 kb upstream appeared to have caused ASIP hypofunction by nonproductive splicing, which the new (polymorphic) SVA insertion largely eliminated. Extended haplotype homozygosity indicated that the insertion allele has risen to allele frequencies up to 11% in European populations over the past several thousand years. These results indicate that a sequence of retrotransposon insertions contributed to a species-wide increase, then a local decrease, of human pigmentation.

Exploring the Maintaining Period and the Differentially Expressed Genes between the Yellow and Black Stripes of the Juvenile Stripe in the Offspring of Wild Boar and Duroc

Coloration is a crucial trait that allows species to adapt and survive in different environments. Wild boars exhibit alternating black (dark) and yellow (light) longitudinal stripes on their back during their infancy (juvenile stripes), and as adults, they transform into uniform wild-type coat color. Aiming to record the procedure of juvenile stripes disappearing, piglets (WD) with juvenile stripes were produced by crossing a wild boar with Duroc sows, and photos of their coat color were taken from 20 d to 220 d. The pigments in the hairs from the black and yellow stripes were determined. Furthermore, the differentially expressed genes between the black and yellow stripes were investigated in 5 WD with the age of 30 d using whole-transcriptome sequencing to explore the genetic mechanism of the juvenile stripes. The juvenile stripes started to disappear at about 70 d, and stripes were not distinguished with the naked eye at about 160 d; that is, the juvenile stripe completely disappeared. A hotspot of a differentially expressing (DE) region was found on chromosome 13, containing/covering 2 of 13 DE genes and 8 of 10 DE lncRNAs in this region. A network among ZIC4, ssc-miR-532-3p, and ENSSSCG00000056225 might regulate the formation of juvenile stripes. Altogether, this study provides new insights into spatiotemporal coat color pattern.

Coat colour in marsupials: genetic variants at the ASIP locus determine grey and black fur of the brushtail possum

The possession of fur or hair is a defining characteristic of mammals and can occur in a variety of colours and patterns. While genetic determinants of coat colour are well described in eutherian 'placental' mammals, the other major mammalian infraclass, marsupials, is grossly understudied. The fur of the common brushtail possum (Trichosurus vulpecula), an iconic native mammal found throughout Australia and introduced into Aotearoa New Zealand, possesses two main colour morphs: grey and black. To identify genetic variants associated with coat colour, we performed a genome-wide association study (GWAS) with genotype by sequencing (GBS) data. Single nucleotide variants (SNVs) on chromosome 3, close to the agouti signalling protein (ASIP) gene that controls the temporal and spatial distribution of pigments in eutherian mammals, were identified. Fine-mapping identified a C>T variant at chr3:100483705 that results in a ASIP:p.Arg115Cys missense substitution, and animals homozygous for this variant have black fur. In addition to uncovering the first genetic determinant of coat colour in a natural marsupial population, comparative analysis of ASIP in divergent marsupial species identified the dasyurids as having accelerated evolution, reflecting their well described diversity of coat colour and pattern.

Ancestry dynamics and trait selection in a designer cat breed

The Bengal cat breed was developed from intercrosses between the Asian leopard cat, Prionailurus bengalensis, and the domestic cat, Felis catus, with a last common ancestor approximately 6 million years ago. Predicted to derive ∼94% of their genome from domestic cats, regions of the leopard cat genome are thought to account for the unique pelage traits and ornate color patterns of the Bengal breed, which are similar to those of ocelots and jaguars. We explore ancestry distribution and selection signatures in the Bengal breed by using reduced representation and whole-genome sequencing from 947 cats. The mean proportion of leopard cat DNA in the Bengal breed is 3.48%, lower than predicted from breed history, and is broadly distributed, covering 93% of the Bengal genome. Overall, leopard cat introgressions do not show strong signatures of selection across the Bengal breed. However, two popular color traits in Bengal cats, charcoal and pheomelanin intensity, are explained by selection of leopard cat genes whose expression is reduced in a domestic cat background, consistent with genetic incompatibility resulting from hybridization. We characterize several selective sweeps in the Bengal genome that harbor candidate genes for pelage and color pattern and that are associated with domestic, rather than leopard, cat haplotypes. We identify the molecular and phenotypic basis of one selective sweep as reduced expression of the Fgfr2 gene, which underlies glitter, a trait desired by breeders that affects hair texture and light reflectivity.

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.

An evolutionarily distinct ringed seal in the Ilulissat Icefjord

The Earth's polar regions are low rates of inter- and intraspecific diversification. An extreme mammalian example is the Arctic ringed seal (Pusa hispida hispida), which is assumed to be panmictic across its circumpolar Arctic range. Yet, local Inuit communities in Greenland and Canada recognize several regional variants; a finding supported by scientific studies of body size variation. It is however unclear whether this phenotypic variation reflects plasticity, morphs or distinct ecotypes. Here, we combine genomic, biologging and survey data, to document the existence of a unique ringed seal ecotype in the Ilulissat Icefjord (locally 'Kangia'), Greenland; a UNESCO World Heritage site, which is home to the most productive marine-terminating glacier in the Arctic. Genomic analyses reveal a divergence of Kangia ringed seals from other Arctic ringed seals about 240 kya, followed by secondary contact since the Last Glacial Maximum. Despite ongoing gene flow, multiple genomic regions appear under strong selection in Kangia ringed seals, including candidate genes associated with pelage coloration, growth and osmoregulation, potentially explaining the Kangia seal's phenotypic and behavioural uniqueness. The description of 'hidden' diversity and adaptations in yet another Arctic species merits a reassessment of the evolutionary processes that have shaped Arctic diversity and the traditional view of this region as an evolutionary freezer. Our study highlights the value of indigenous knowledge in guiding science and calls for efforts to identify distinct populations or ecotypes to understand how these might respond differently to environmental change.

Exploring the Domestication Syndrome Hypothesis in Dogs: Pigmentation Does Not Predict Cortisol Levels

Previous research has found connections between pigmentation, behavior, and the physiological stress response in both wild and domestic animals; however, to date, no extensive research has been devoted to answering these questions in domestic dogs. Modern dogs are exposed to a variety of stressors; one well-studied stressor is residing in an animal shelter. To explore the possible relationships between dogs' responses to stress and their pigmentation, we conducted statistical analyses of the cortisol:creatinine ratios of 208 American shelter dogs as a function of their coat color/pattern, eumelanin pigmentation, or white spotting. These dogs had been enrolled in previous welfare studies investigating the effect of interventions during which they left the animal shelter and spent time with humans. In the current investigation, we visually phenotype dogs based on photographs in order to classify their pigmentation and then conduct post hoc analyses to examine whether they differentially experience stress as a function of pigmentation. We found that the dogs did not differ significantly in their urinary cortisol:creatinine ratios based on coat color/pattern, eumelanin pigmentation, or white spotting, either while they were residing in the animal shelter or during the human interaction intervention. These preliminary data suggest that pigmentation alone does not predict the stress responses of shelter dogs; however, due to the small sample size and retrospective nature of the study, more research is needed.

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.

A Chinese indicine pangenome reveals a wealth of novel structural variants introgressed from other Bos species

Chinese indicine cattle harbor a much higher genetic diversity compared with other domestic cattle, but their genome architecture remains uninvestigated. Using PacBio HiFi sequencing data from 10 Chinese indicine cattle across southern China, we assembled 20 high-quality partially phased genomes and integrated them into a multiassembly graph containing 148.5 Mb (5.6%) of novel sequence. We identified 156,009 high-confidence nonredundant structural variants (SVs) and 206 SV hotspots spanning ∼195 Mb of gene-rich sequence. We detected 34,249 archaic introgressed fragments in Chinese indicine cattle covering 1.93 Gb (73.3%) of the genome. We inferred an average of 3.8%, 3.2%, 1.4%, and 0.5% of introgressed sequence originating, respectively, from banteng-like, kouprey-like, gayal-like, and gaur-like Bos species, as well as 0.6% of unknown origin. Introgression from multiple donors might have contributed to the genetic diversity of Chinese indicine cattle. Altogether, this study highlights the contribution of interspecies introgression to the genomic architecture of an important livestock population and shows how exotic genomic elements can contribute to the genetic variation available for selection.

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.

The LASSIE MPS panel: Predicting externally visible traits in dogs for forensic purposes

Predicting the outward appearance of dogs via their DNA, also known as Canine DNA Phenotyping, is a young, emerging field of research in forensic genetics. The few previous studies published in this respect were restricted to the consecutive analysis of single DNA markers, a process that is time- and sample-consuming and therefore not a viable option for limited forensic specimens. Here, we report on the development and evaluation of a Massively Parallel Sequencing (MPS) based molecular genetic assay, the LASSIE MPS Panel. This panel aims to predict externally visible as well as skeletal traits, which include coat color, coat pattern, coat structure, tail morphology, skull shape, ear shape, eye color and body size from DNA using 44 genetic markers in a single molecular genetic assay. A biostatistical naïve Bayes classification approach was applied to identify the most informative marker combinations for predicting phenotypes. Overall, the predictive performance was characterized by a very high classification success for some of the trait categories, and high to moderate success for others. The performance of the developed predictive framework was further evaluated using blind samples from three randomly selected dog individuals, whose appearance was well predicted.

The influence of transposable elements on animal colouration

Transposable elements (TEs) are mobile genetic sequences present within host genomes. TEs can contribute to the evolution of host traits, since transposition is mutagenic and TEs often contain host regulatory and protein coding sequences. We review cases where TEs influence animal colouration, reporting major patterns and outstanding questions. TE-induced colouration phenotypes typically arise via introduction of novel regulatory sequences and splice sites, affecting pigment cell development or pigment synthesis. We discuss if particular TE types may be more frequently involved in the evolution of colour variation in animals, given that examples involving long terminal repeat (LTR) elements appear to dominate. Currently, examples of TE-induced colouration phenotypes in animals mainly concern model and domesticated insect and mammal species. However, several influential recent examples, coupled with increases in genome sequencing, suggest cases reported from wild species will increase considerably.

Possible origins and implications of atypical morphologies and domestication-like traits in wild golden jackals (Canis aureus)

Deciphering the origins of phenotypic variations in natural animal populations is a challenging topic for evolutionary and conservation biologists. Atypical morphologies in mammals are usually attributed to interspecific hybridisation or de-novo mutations. Here we report the case of four golden jackals (Canis aureus), that were observed during a camera-trapping wildlife survey in Northern Israel, displaying anomalous morphological traits, such as white patches, an upturned tail, and long thick fur which resemble features of domesticated mammals. Another individual was culled under permit and was genetically and morphologically examined. Paternal and nuclear genetic profiles, as well as geometric morphometric data, identified this individual as a golden jackal rather than a recent dog/wolf-jackal hybrid. Its maternal haplotype suggested past introgression of African wolf (Canis lupaster) mitochondrial DNA, as previously documented in other jackals from Israel. When viewed in the context of the jackal as an overabundant species in Israel, the rural nature of the surveyed area, the abundance of anthropogenic waste, and molecular and morphological findings, the possibility of an individual presenting incipient stages of domestication should also be considered.

Genomic Diversity and Runs of Homozygosity in Bernese Mountain Dogs

Bernese mountain dogs are a large dog breed formed in the early 1900s in Switzerland. While originally farm dogs that were used for pulling carts, guarding, and driving cattle, today they are considered multi-purpose companion and family dogs. The breed is predisposed to several complex diseases, such as histiocytic sarcoma, degenerative myelopathy, or hip dysplasia. Using whole-genome sequencing (WGS) data, we assessed the genomic architecture of 33 unrelated dogs from four countries: France, Sweden, Switzerland, and the United States. Analysis of runs of homozygosity (ROH) identified 12,643 ROH with an average length of 2.29 Mb and an average inbreeding coefficient of 0.395. Multidimensional scaling analysis of the genetic relatedness revealed limited clustering of European versus USA dogs, suggesting exchanges of breeding stock between continents. Furthermore, only two mtDNA haplotypes were detected in the 33 studied dogs, both of which are widespread throughout multiple dog breeds. WGS-based ROH analyses revealed several fixed or nearly fixed regions harboring discreet morphological trait-associated as well as disease-associated genetic variants. Several genes involved in the regulation of immune cells were found in the ROH shared by all dogs, which is notable in the context of the breed’s strong predisposition to hematopoietic cancers. High levels of inbreeding and relatedness, strongly exaggerated in the last 30 years, have likely led to the high prevalence of specific genetic disorders in this breed.

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.

Deciphering the Genetic Basis of Silkworm Cocoon Colors Provides New Insights into Biological Coloration and Phenotypic Diversification

The genetic basis of phenotypic variation is a long-standing concern of evolutionary biology. Coloration has proven to be a visual, easily quantifiable, and highly tractable system for genetic analysis and is an ever-evolving focus of biological research. Compared with the homogenized brown-yellow cocoons of wild silkworms, the cocoons of domestic silkworms are spectacularly diverse in color, such as white, green, and yellow-red; this provides an outstanding model for exploring the phenotypic diversification and biological coloration. Herein, the molecular mechanism underlying silkworm green cocoon formation was investigated, which was not fully understood. We demonstrated that five of the seven members of a sugar transporter gene cluster were specifically duplicated in the Bombycidae and evolved new spatial expression patterns predominantly expressed in silk glands, accompanying complementary temporal expression; they synergistically facilitate the uptake of flavonoids, thus determining the green cocoon. Subsequently, polymorphic cocoon coloring landscape involving multiple loci and the evolution of cocoon color from wild to domestic silkworms were analyzed based on the pan-genome sequencing data. It was found that cocoon coloration involved epistatic interaction between loci; all the identified cocoon color-related loci existed in wild silkworms; the genetic segregation, recombination, and variation of these loci shaped the multicolored cocoons of domestic silkworms. This study revealed a new mechanism for flavonoids-based biological coloration that highlights the crucial role of gene duplication followed by functional diversification in acquiring new genetic functions; furthermore, the results in this work provide insight into phenotypic innovation during domestication.

Evaluation of Genetic Diversity in Dog Breeds Using Pedigree and Molecular Analysis: A Review

Domestic dogs are important for many economic and social reasons, and they have become a well-known model species for human disease. According to research, dog breeds exhibit significant levels of inbreeding and genetic diversity loss, decreasing the population’s ability to adapt in certain conditions, and indicating the need of conservation strategies. Before the development of molecular markers, pedigree information was used for genetic diversity management. In recent years, genomic tools are frequently applied for accurate estimation of genetic diversity and improved genetic conservation due to incomplete pedigrees and pedigree errors. The most frequently used molecular markers include PCR-based microsatellite markers (STRs) and DNA sequencing-based single-nucleotide polymorphism markers (SNP). The aim of this review was to highlight genetic diversity studies on dog breeds conducted using pedigree and molecular markers, as well as the importance of genetic diversity conservation in increasing the adaptability and survival of dog breed populations.

Ruminant-specific retrotransposons shape regulatory evolution of bovine immunity

Cattle are an important livestock species, and mapping the genomic architecture of agriculturally relevant traits such as disease susceptibility is a major challenge in the bovine research community. Lineage-specific transposable elements (TEs) are increasingly recognized to contribute to gene regulatory evolution and variation, but this possibility has been largely unexplored in ruminant genomes. We conducted epigenomic profiling of the type II interferon (IFN) response in bovine cells and found thousands of ruminant-specific TEs including MER41_BT and Bov-A2 elements predicted to act as IFN-inducible enhancer elements. CRISPR knockout experiments in bovine cells established that critical immune factors including IFNAR2 and IL2RB are transcriptionally regulated by TE-derived enhancers. Finally, population genomic analysis of 38 individuals revealed that a subset of polymorphic TE insertions may function as enhancers in modern cattle. Our study reveals that lineage-specific TEs have shaped the evolution of ruminant IFN responses and potentially continue to contribute to immune gene regulatory differences across modern breeds and individuals. Together with previous work in human cells, our findings demonstrate that lineage-specific TEs have been independently co-opted to regulate IFN-inducible gene expression in multiple species, supporting TE co-option as a recurrent mechanism driving the evolution of IFN-inducible transcriptional networks.

An enhancer of Agouti contributes to parallel evolution of cryptically colored beach mice

Identifying the genetic basis of repeatedly evolved traits provides a way to reconstruct their evolutionary history and ultimately investigate the predictability of evolution. Here, we focus on the oldfield mouse (Peromyscus polionotus), which occurs in the southeastern United States, where it exhibits considerable color variation. Dorsal coats range from dark brown in mainland mice to near white in mice inhabiting sandy beaches; this light pelage has evolved independently on Florida's Gulf and Atlantic coasts as camouflage from predators. To facilitate genomic analyses, we first generated a chromosome-level genome assembly of Peromyscus polionotus subgriseus. Next, in a uniquely variable mainland population (Peromyscus polionotus albifrons), we scored 23 pigment traits and performed targeted resequencing in 168 mice. We find that pigment variation is strongly associated with an ∼2-kb region ∼5 kb upstream of the Agouti signaling protein coding region. Using a reporter-gene assay, we demonstrate that this regulatory region contains an enhancer that drives expression in the dermis of mouse embryos during the establishment of pigment prepatterns. Moreover, extended tracts of homozygosity in this Agouti region indicate that the light allele experienced recent and strong positive selection. Notably, this same light allele appears fixed in both Gulf and Atlantic coast beach mice, despite these populations being separated by >1,000 km. Together, our results suggest that this identified Agouti enhancer allele has been maintained in mainland populations as standing genetic variation and from there, has spread to and been selected in two independent beach mouse lineages, thereby facilitating their rapid and parallel evolution.

A 13.42-kb tandem duplication at the ASIP locus is strongly associated with the depigmentation phenotype of non-classic Swiss markings in goats

BACKGROUND

The pigmentation phenotype diversity is rich in domestic goats, and identification of the genetic loci affecting coat color in goats has long been of interest. Via the detections of selection signatures, a duplication upstream ASIP was previously reported to be a variant affecting the Swiss markings depigmentation phenotype in goats.

RESULTS

We conducted a genome-wide association study using whole-genome sequencing (WGS) data to identify the genetic loci and causal variants affecting the pigmentation phenotype in 65 Jintang black (JT) goats (i.e., 48 solid black vs. 17 non-classic Swiss markings). Although a single association peak harboring the ASIP gene at 52,619,845-72,176,538 bp on chromosome 13 was obtained using a linear mixed model approach, all the SNPs and indels in this region were excluded as causal variants for the pigmentation phenotype. We then found that all 17 individuals with non-classic Swiss markings carried a 13,420-bp duplication (CHI13:63,129,198-63,142,617 bp) nearly 101 kb upstream of ASIP, and this variant was strongly associated (P = 1.48 × 10^- 12) with the coat color in the 65 JT goats. The copy numbers obtained from the WGS data also showed that the duplication was present in all 53 goats from three European breeds with Swiss markings and absent in 45 of 51 non-Swiss markings goats from four other breeds and 21 Bezoars, which was further validated in 314 samples from seven populations based on PCR amplification. The copy numbers of the duplication vary in different goat breeds with Swiss markings, indicating a threshold effect instead of a dose-response effect at the molecular level. Furthermore, breakpoint flanking repeat analysis revealed that the duplication was likely to be a result of the Bov-B-mediated nonallelic homologous recombination.

CONCLUSION

We confirmed that a genomic region harboring the ASIP gene is a major locus affecting the coat color phenotype of Swiss markings in goats. Although the molecular genetic mechanisms remain unsolved, the 13,420-bp duplication upstream of ASIP is a necessary but not sufficient condition for this phenotype in goats. Moreover, the variations in the copy number of the duplication across different goat breeds do not lead to phenotypic heterogeneity.

ASIP Promoter Variants Predict the Sesame Coat Color in Shiba Inu Dogs

Animals exhibit a wide variety of genetically determined coat colors and pigmentation patterns that serve important roles in adaptation and communication. Although the genetics of the main coat colors in dogs have been studied extensively, there are types of coat pigmentation that have not been explained yet. Recently, an association between the variants in the ASIP gene Ventral (VP) and Hair Cycle (HCP) promoters with different coat colors in dogs has been established. Here, we used the new findings as a basis to investigate the genetics of the red sesame coat color in Shiba Inu dogs. Our study revealed that red sesame dogs carry a specific heterozygous ASIP promoter diplotype, VP2-HCP1/VP2-HCP3, where VP2-HCP1 is responsible for the red coat with a dark overlay, and VP2-HCP3 for a tan point-like pattern. This finding explains the inheritance of this coat color pattern and can be used by breeders to produce dogs with this rare phenotype. A comparison of sesame dogs (VP2-HCP1/VP2-HCP3) to a dog homozygous for the VP2-HCP1 promoter haplotype suggests that the incomplete dominance between the ASIP alleles may be involved in the sesame coat formation. These results are in good agreement with the new model explaining how different levels of ASIP gene expression affect the regulation of pigment synthesis in melanocytes.

Natural and human-driven selection of a single non-coding body size variant in ancient and modern canids

Domestic dogs (Canis lupus familiaris) are the most variable-sized mammalian species on Earth, displaying a 40-fold size difference between breeds.¹ Although dogs of variable size are found in the archeological record,^2-4 the most dramatic shifts in body size are the result of selection over the last two centuries, as dog breeders selected and propagated phenotypic extremes within closed breeding populations.⁵ Analyses of over 200 domestic breeds have identified approximately 20 body size genes regulating insulin processing, fatty acid metabolism, TGFβ signaling, and skeletal formation.^6-10 Of these, insulin-like growth factor 1 (IGF1) predominates, controlling approximately 15% of body size variation between breeds.⁸ The identification of a functional mutation associated with IGF1 has thus far proven elusive.^6,10,11 Here, to identify and elucidate the role of an ancestral IGF1 allele in the propagation of modern canids, we analyzed 1,431 genome sequences from 13 species, including both ancient and modern canids, thus allowing us to define the evolutionary history of both ancestral and derived alleles at this locus. We identified a single variant in an antisense long non-coding RNA (IGF1-AS) that interacts with the IGF1 gene, creating a duplex. While the derived mutation predominates in both modern gray wolves and large domestic breeds, the ancestral allele, which predisposes to small size, was common in small-sized breeds and smaller wild canids. Our analyses demonstrate that this major regulator of canid body size nearly vanished in Pleistocene wolves, before its recent resurgence resulting from human-imposed selection for small-sized breed dogs.