Paul Glyn Jones

Genome-wide SNP and haplotype analyses reveal a rich history underlying dog domestication

Advances in genome technology have facilitated a new understanding of the historical and genetic processes crucial to rapid phenotypic evolution under domestication. To understand the process of dog diversification better, we conducted an extensive genome-wide survey of more than 48,000 single nucleotide polymorphisms in dogs and their wild progenitor, the grey wolf. Here we show that dog breeds share a higher proportion of multi-locus haplotypes unique to grey wolves from the Middle East, indicating that they are a dominant source of genetic diversity for dogs rather than wolves from east Asia, as suggested by mitochondrial DNA sequence data. Furthermore, we find a surprising correspondence between genetic and phenotypic/functional breed groupings but there are exceptions that suggest phenotypic diversification depended in part on the repeated crossing of individuals with novel phenotypes. Our results show that Middle Eastern wolves were a critical source of genome diversity, although interbreeding with local wolf populations clearly occurred elsewhere in the early history of specific lineages. More recently, the evolution of modern dog breeds seems to have been an iterative process that drew on a limited genetic toolkit to create remarkable phenotypic diversity.

An Expressed Fgf4 Retrogene Is Associated with Breed-Defining Chondrodysplasia in Domestic Dogs

Going Retro In a year celebrating Darwin, the question of how new functional genes arise during evolution is of particular interest. Through a multibreed genetic analysis of the domestic dog, Parker et al. (p. 995, published online 16 July; see the Perspective by Kaessmann) find that the short-legged phenotype that characterizes at least 19 common dog breeds, including the corgi, dachshund, and basset hound, is specifically associated with the expression in developing bone of a gene encoding fibroblast growth factor 4 (fgf4), a member of a gene family previously implicated in dwarfism in humans. Interestingly, the culprit fgf4 gene in dogs has the hallmarks of a “retrogene,” a gene that arises when a parental gene is duplicated through an RNA-based copying mechanism. The short legs that characterize certain dog breeds are associated with a gene that arose recently by RNA-based gene duplication. Retrotransposition of processed mRNAs is a common source of novel sequence acquired during the evolution of genomes. Although the vast majority of retroposed gene copies, or retrogenes, rapidly accumulate debilitating mutations that disrupt the reading frame, a small percentage become new genes that encode functional proteins. By using a multibreed association analysis in the domestic dog, we demonstrate that expression of a recently acquired retrogene encoding fibroblast growth factor 4 (fgf4) is strongly associated with chondrodysplasia, a short-legged phenotype that defines at least 19 dog breeds including dachshund, corgi, and basset hound. These results illustrate the important role of a single evolutionary event in constraining and directing phenotypic diversity in the domestic dog.

Genetic Mapping of Fixed Phenotypes: Disease Frequency as a Breed Characteristic

Traits that have been stringently selected to conform to specific criteria in a closed population are phenotypic stereotypes. In dogs, Canis familiaris, such stereotypes have been produced by breeding for conformation, performance (behaviors), etc. We measured phenotypes on a representative sample to establish breed stereotypes. DNA samples from 147 dog breeds were used to characterize single nucleotide polymorphism allele frequencies for association mapping of breed stereotypes. We identified significant size loci (quantitative trait loci [QTLs]), implicating candidate genes appropriate to regulation of size (e.g., IGF1, IGF2BP2 SMAD2, etc.). Analysis of other morphological stereotypes, also under extreme selection, identified many additional significant loci. Behavioral loci for herding, pointing, and boldness implicated candidate genes appropriate to behavior (e.g., MC2R, DRD1, and PCDH9). Significant loci for longevity, a breed characteristic inversely correlated with breed size, were identified. The power of this approach to identify loci regulating the incidence of specific polygenic diseases is demonstrated by the association of a specific IGF1 haplotype with hip dysplasia, patella luxation, and pancreatitis.

Gene (mRNA) expression in canine atopic dermatitis: microarray analysis

Genes potentially involved in the pathology of canine atopic dermatitis (AD) were identified using gene expression microarrays. Total RNA extracted from skin biopsies was hybridized to an Agilent Technologies custom-designed 22K canine array. The arrays were analysed using Genedata Analyst software. Data were corrected for multiple hypothesis testing and tested for significance using the National Institute on Aging array analysis tool. For comparison, data were divided into separate groups: lesional atopic (n = 16), nonlesional atopic (n = 17) and healthy controls (n = 9). Fifty-four genes were differentially expressed at a significance level of 0.05 in canine AD compared to healthy controls. Sixteen genes were differentially expressed in both nonlesional and lesional atopic skin, 26 genes only in nonlesional skin and 12 only in lesional skin. These genes were associated with innate immune and inflammatory responses, cell cycle, apoptosis, barrier formation and transcriptional regulation. The most dysregulated gene in lesional skin was S100A8, which showed an almost 23-fold increase in expression. This is a pro-inflammatory cytokine located in the epidermal differentiation complex. Microarray analysis is a novel technique in canine AD. Significant changes in gene expression were identified in atopic skin. These were relevant to skin barrier formation and the immune response, suggesting that they both participate in AD. Gene expression restricted to lesional skin may be involved in inflammatory changes, whereas those shared or restricted to nonlesional skin may reflect the atopic phenotype. Investigating gene polymorphisms in the targets identified in this study will help improve our understanding of the genetic basis of this disease.

Canine hip dysplasia is predictable by genotyping

OBJECTIVE To establish a predictive method using whole genome genotyping for early intervention in canine hip dysplasia (CHD) risk management, for the prevention of the progression of secondary osteoarthritis (OA), and for selective breeding. DESIGN Two sets of dogs (six breeds) were genotyped with dense SNPs covering the entire canine genome. The first set contained 359 dogs upon which a predictive formula for genomic breeding value (GBV) was derived by using their estimated breeding value (EBV) of the Norberg angle (a measure of CHD) and their genotypes. To investigate how well the formula would work for an individual dog with genotype only (without using EBV), a cross validation was performed by masking the EBV of one dog at a time. The genomic data and the EBV of the remaining dogs were used to predict the GBV for the single dog that was left out. The second set of dogs included 38 new Labrador retriever dogs, which had no pedigree relationship to the dogs in the first set. RESULTS The cross validation showed a strong correlation (R>0.7) between the EBV and the GBV. The independent validation showed a moderate correlation (R=0.5) between GBV for the Norberg angle and the observed Norberg angle (no EBV was available for the new 38 dogs). Sensitivity, specificity, positive and negative predictive values of the genomic data were all above 70%. CONCLUSIONS Prediction of CHD from genomic data is feasible, and can be applied for risk management of CHD and early selection for genetic improvement to reduce the prevalence of CHD in breeding programs. The prediction can be implemented before maturity, at which age current radiographic screening programs are traditionally applied, and as soon as DNA is available.

Deep Illumina-Based Shotgun Sequencing Reveals Dietary Effects on the Structure and Function of the Fecal Microbiome of Growing Kittens

Background Previously, we demonstrated that dietary protein:carbohydrate ratio dramatically affects the fecal microbial taxonomic structure of kittens using targeted 16S gene sequencing. The present study, using the same fecal samples, applied deep Illumina shotgun sequencing to identify the diet-associated functional potential and analyze taxonomic changes of the feline fecal microbiome. Methodology & Principal Findings Fecal samples from kittens fed one of two diets differing in protein and carbohydrate content (high–protein, low–carbohydrate, HPLC; and moderate-protein, moderate-carbohydrate, MPMC) were collected at 8, 12 and 16 weeks of age (n = 6 per group). A total of 345.3 gigabases of sequence were generated from 36 samples, with 99.75% of annotated sequences identified as bacterial. At the genus level, 26% and 39% of reads were annotated for HPLC- and MPMC-fed kittens, with HPLC-fed cats showing greater species richness and microbial diversity. Two phyla, ten families and fifteen genera were responsible for more than 80% of the sequences at each taxonomic level for both diet groups, consistent with the previous taxonomic study. Significantly different abundances between diet groups were observed for 324 genera (56% of all genera identified) demonstrating widespread diet-induced changes in microbial taxonomic structure. Diversity was not affected over time. Functional analysis identified 2,013 putative enzyme function groups were different (p<0.000007) between the two dietary groups and were associated to 194 pathways, which formed five discrete clusters based on average relative abundance. Of those, ten contained more (p<0.022) enzyme functions with significant diet effects than expected by chance. Six pathways were related to amino acid biosynthesis and metabolism linking changes in dietary protein with functional differences of the gut microbiome. Conclusions These data indicate that feline feces-derived microbiomes have large structural and functional differences relating to the dietary protein:carbohydrate ratio and highlight the impact of diet early in life.

Molecular recognition of metabotropic glutamate receptor type 1 (mGluR1): synergistic understanding with free energy perturbation and linear response modeling.

Metabotropic glutamate receptors (mGluRs) constitute an important family of the G-protein coupled receptors. Due to their widespread distribution in the central nervous system (CNS), these receptors are attractive candidates for understanding the molecular basis of various cognitive processes as well as for designing inhibitors for relevant psychiatric and neurological disorders. Despite many studies on drugs targeting the mGluR receptors to date, the molecular level details on the ligand binding dynamics still remain unclear. In this study, we performed in silico experiments for mGluR1 with 29 different ligands including known synthetic agonists and antagonists as well as natural amino acids. The ligand-receptor binding affinities were estimated by the use of atomistic simulations combined with the mathematically rigorous, Free Energy Perturbation (FEP) method, which successfully recognized the native agonist l-glutamate among the highly favorable binders, and also accurately distinguished antagonists from agonists. Comparative contact analysis also revealed the binding mode differences between natural and non-natural amino acid-based ligands. Several factors potentially affecting the ligand binding affinity and specificity were identified including net charges, dipole moments, and the presence of aromatic rings. On the basis of these findings, linear response models (LRMs) were built for different sets of ligands that showed high correlations (R(2) > 0.95) to the corresponding FEP binding affinities. These results identify some key factors that determine ligand-mGluR1 binding and could be used for future inhibitor designs and support a role for in silico modeling for understanding receptor ligand interactions.