Analysis of large versus small dogs reveals three genes on the canine X chromosome associated with body weight, muscling and back fat thickness

The Demography, Longevity and Mortality of Bullmastiffs Attending Veterinary Practices in Australia

The Bullmastiff is a giant breed of dog, and there is evidence the breed is predisposed to cancer and musculoskeletal disorders, though the disease investigation of the breed is limited. This study reports on the demography, longevity and mortality of Bullmastiffs attending veterinary practices in Australia over a ten-year period. VetCompass Australia collects patient data from veterinary practices across Australia for epidemiological analysis. All patient records of Bullmastiffs available in the VetCompass Australia database during this decade period were reviewed, with demographic information on the breed inclusive of coat colour, sex, neuter status, weight and location collated. Standardised veterinary diagnostic (VeNom) codes for the most appropriate cause of death were assigned to deceased dogs. The population comprised 2771 Bullmastiffs with an overall median age of 2.8 years. Within the group, 1259 were female (45.4%), 1491 were male (53.8%), and 21 dogs (0.8%) had no recorded sex or neuter status. Dogs grew rapidly in their first year, with an average gain of approximately 1 kg every 10 days. A slower growth rate continued in their second year, and growth plateaued as adulthood was reached, with the mean body weight of adult male dogs (46.6 kg) being heavier than that of females (40.5 kg). The age at death for the group was 8.5 years. The most common causes of death in the breed were mass lesions (28.2%), old age (9.9%), musculoskeletal-related disease (9.9%) and neurological (5.3%) and behavioural disorders (4.8%). Neutering was protective against mortality from urogenital causes (OR: 0.14; CI: 0.02-0.52; p = 0.003) and had a positive effect on longevity. This study provides demographic and health information on a population of Bullmastiffs attending veterinary practices in Australia, which will benefit evidence-based veterinary decisions for this breed. Additionally, the results may assist owners and breeders in making informed decisions on health risks and breeding programmes in the population.

Elucidation of the radius and ulna fracture mechanisms in toy poodle dogs using finite element analysis

Fractures occurring in the distal radius and ulna of toy breed dogs pose distinctive challenges for veterinary practitioners, requiring specialized treatment approaches primarily based on anatomical features. Finite Element Analysis (FEA) was applied to conduct numerical experiments to determine stress distribution across the bone. This methodology offers an alternative substitute for directly investigating these phenomena in living dog experiments, which could present ethical obstacles. A three-dimensional bone model of the metacarpal, carpal, radius, ulna, and humerus was reconstructed from Computed Tomography (CT) images of the toy poodle and dachshund forelimb. The model was designed to simulate the jumping and landing conditions from a vertical distance of 40 cm to the ground within a limited timeframe. The investigation revealed considerable variations in stress distribution patterns between the radius and ulna of toy poodles and dachshunds, indicating notably elevated stress levels in toy poodles compared to dachshunds. In static and dynamic stress analysis, toy poodles exhibit peak stress levels at the distal radius and ulna. The Von Mises stresses for toy poodles reach 90.07 MPa (static) and 1,090.75 MPa (dynamic) at the radius and 1,677.97 MPa (static) and 1,047.98 MPa (dynamic) at the ulna. Conversely, dachshunds demonstrate lower stress levels for 5.39 MPa (static) and 231.79 MPa (dynamic) at the radius and 390.56 MPa (static) and 513.28 MPa (dynamic) at the ulna. The findings offer valuable insights for modified treatment approaches in managing fractures in toy breed dogs, optimizing care and outcomes.

Genome wide association study in Swedish Labrador retrievers identifies genetic loci associated with hip dysplasia and body weight

Genome wide association studies (GWAS) have been utilized to identify genetic risk loci associated with both simple and complex inherited disorders. Here, we performed a GWAS in Labrador retrievers to identify genetic loci associated with hip dysplasia and body weight. Hip dysplasia scores were available for 209 genotyped dogs. We identified a significantly associated locus for hip dysplasia on chromosome 24, with three equally associated SNPs (p = 4.3 × 10-7) in complete linkage disequilibrium located within NDRG3, a gene which in humans has been shown to be differentially expressed in osteoarthritic joint cartilage. Body weight, available for 85 female dogs, was used as phenotype for a second analysis. We identified two significantly associated loci on chromosome 10 (p = 4.5 × 10-7) and chromosome 31 (p = 2.5 × 10-6). The most associated SNPs within these loci were located within the introns of the PRKCE and CADM2 genes, respectively. PRKCE has been shown to play a role in regulation of adipogenesis whilst CADM2 has been associated with body weight in multiple human GWAS. In summary, we identified credible candidate loci explaining part of the genetic inheritance for hip dysplasia and body weight in Labrador retrievers with strong candidate genes in each locus previously implicated in the phenotypes investigated.

Uncovering structural variants associated with body weight and obesity risk in labrador retrievers: a genome-wide study

Although obesity in the domestic dog (Canis lupus familiaris) is known to decrease well-being and shorten lifespan, the genetic risk variants associated with canine obesity remain largely unknown. In our study, which focused on the obesity-prone Labrador Retriever breed, we conducted a genome-wide analysis to identify structural variants linked to body weight and obesity. Obesity status was based on a 5-point body condition score (BCS) and the obese dog group included all dogs with a BCS of 5, along with dogs with the highest body weight within the BCS 4 group. Data from whole-gene sequencing of fifty dogs, including 28 obese dogs, were bioinformatically analyzed to identify potential structural variants that varied in frequency between obese and healthy dogs. The seven most promising variants were further analyzed by droplet digital PCR in a group of 110 dogs, including 63 obese. Our statistical evidence suggests that common structural mutations in or near six genes, specifically ALPL, KCTD8, SGSM1, SLC12A6, RYR3, and VPS26C, may contribute to the variability observed in body weight and body condition scores among Labrador Retriever dogs. These findings emphasize the need for additional research to validate the associations and explore the specific functions of these genes in relation to canine obesity.

DNA methylation networks underlying mammalian traits

Using DNA methylation profiles (n = 15,456) from 348 mammalian species, we constructed phyloepigenetic trees that bear marked similarities to traditional phylogenetic ones. Using unsupervised clustering across all samples, we identified 55 distinct cytosine modules, of which 30 are related to traits such as maximum life span, adult weight, age, sex, and human mortality risk. Maximum life span is associated with methylation levels in HOXL subclass homeobox genes and developmental processes and is potentially regulated by pluripotency transcription factors. The methylation state of some modules responds to perturbations such as caloric restriction, ablation of growth hormone receptors, consumption of high-fat diets, and expression of Yamanaka factors. This study reveals an intertwined evolution of the genome and epigenome that mediates the biological characteristics and traits of different mammalian species.

Evaluating instruments for assessing healthspan: a multi-center cross-sectional study on health-related quality of life (HRQL) and frailty in the companion dog

Developing valid tools that assess key determinants of canine healthspan such as frailty and health-related quality of life (HRQL) is essential to characterizing and understanding aging in dogs. Additionally, because the companion dog is an excellent translational model for humans, such tools can be applied to evaluate gerotherapeutics and investigate mechanisms underlying longevity in both dogs and humans. In this multi-center, cross-sectional study, we investigated the use of a clinical questionnaire (Canine Frailty Index; CFI; Banzato et al., 2019) to assess frailty and an owner assessment tool (VetMetrica HRQL) to evaluate HRQL in 451 adult companion dogs. Results demonstrated validity of the tools by confirming expectations that frailty score increases and HRQL scores deteriorate with age. CFI scores were significantly higher (higher frailty) and HRQL scores significantly lower (worse HRQL) in old dogs (≥ 7 years of age) compared to young dogs (≥ 2 and < 6 years of age). Body size (small < 11.3 kg (25 lbs) or large > 22.7 kg (50 lbs)) was not associated with CFI or total HRQL score. However, older, larger dogs showed faster age-related decline in HRQL scores specific to owner-reported activity and comfort. Findings suggest that the clinician-assessed CFI and owner-reported VetMetrica HRQL are useful tools to evaluate two determinants of healthspan in dogs: the accumulation of frailty and the progressive decline in quality of life. Establishing tools that operationalize the assessment of canine healthspan is critical for the advancement of geroscience and the development of gerotherapeutics that benefit both human and veterinary medicine. Graphical summary of the design, results, and conclusions of the study.

Chromosome-Level Analysis of the Pelochelys cantorii Genome Provides Insights to Its Immunity, Growth and Longevity

The Asian giant soft-shelled turtle, Pelochelys cantorii (Trionychidae), is one of the largest aquatic turtles in China and was designated as a First-Grade Protected Animal in China in 1989. Previous investigation based on a combination of Illumina short-read, PacBio long-read and Hi-C scaffolding technologies acquired a high-quality chromosome-level genome of Pc. cantorii. In this study, comparative genomic analysis between Pc. cantorii and 16 other vertebrate genomes indicated that turtles separated from the ancestor of archosaurians approximately 256.6 (95% highest posterior density interval, 263.6-251.9) million years ago (Mya) (Upper Permian to Triassic) and that Pc. cantorii separated from the ancestor of Pd. sinensis and R. swinhoei approximately 59.3 (95% highest posterior density interval, 64.3-54.3) Mya. Moreover, several candidate genes, such as VWA5A, ABCG2, A2M and IGSF1, associated with tumor suppression, growth and age were expanded, implicating their potential roles in the exceptional longevity of turtles. This new chromosome-level assembly has important scientific value in the study of conservation of Pc. cantorii and also enriches the evolutionary investigation of turtle species.

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.

The Australasian dingo archetype: de novo chromosome-length genome assembly, DNA methylome, and cranial morphology

BACKGROUND

One difficulty in testing the hypothesis that the Australasian dingo is a functional intermediate between wild wolves and domesticated breed dogs is that there is no reference specimen. Here we link a high-quality de novo long-read chromosomal assembly with epigenetic footprints and morphology to describe the Alpine dingo female named Cooinda. It was critical to establish an Alpine dingo reference because this ecotype occurs throughout coastal eastern Australia where the first drawings and descriptions were completed.

FINDINGS

We generated a high-quality chromosome-level reference genome assembly (Canfam_ADS) using a combination of Pacific Bioscience, Oxford Nanopore, 10X Genomics, Bionano, and Hi-C technologies. Compared to the previously published Desert dingo assembly, there are large structural rearrangements on chromosomes 11, 16, 25, and 26. Phylogenetic analyses of chromosomal data from Cooinda the Alpine dingo and 9 previously published de novo canine assemblies show dingoes are monophyletic and basal to domestic dogs. Network analyses show that the mitochondrial DNA genome clusters within the southeastern lineage, as expected for an Alpine dingo. Comparison of regulatory regions identified 2 differentially methylated regions within glucagon receptor GCGR and histone deacetylase HDAC4 genes that are unmethylated in the Alpine dingo genome but hypermethylated in the Desert dingo. Morphologic data, comprising geometric morphometric assessment of cranial morphology, place dingo Cooinda within population-level variation for Alpine dingoes. Magnetic resonance imaging of brain tissue shows she had a larger cranial capacity than a similar-sized domestic dog.

CONCLUSIONS

These combined data support the hypothesis that the dingo Cooinda fits the spectrum of genetic and morphologic characteristics typical of the Alpine ecotype. We propose that she be considered the archetype specimen for future research investigating the evolutionary history, morphology, physiology, and ecology of dingoes. The female has been taxidermically prepared and is now at the Australian Museum, Sydney.

The Australasian dingo archetype: De novo chromosome-length genome assembly, DNA methylome, and cranial morphology

Background

One difficulty in testing the hypothesis that the Australasian dingo is a functional intermediate between wild wolves and domesticated breed dogs is that there is no reference specimen. Here we link a high-quality de novo long read chromosomal assembly with epigenetic footprints and morphology to describe the Alpine dingo female named Cooinda. It was critical to establish an Alpine dingo reference because this ecotype occurs throughout coastal eastern Australia where the first drawings and descriptions were completed.

Findings

We generated a high-quality chromosome-level reference genome assembly (Canfam_ADS) using a combination of Pacific Bioscience, Oxford Nanopore, 10X Genomics, Bionano, and Hi-C technologies. Compared to the previously published Desert dingo assembly, there are large structural rearrangements on Chromosomes 11, 16, 25 and 26. Phylogenetic analyses of chromosomal data from Cooinda the Alpine dingo and nine previously published de novo canine assemblies show dingoes are monophyletic and basal to domestic dogs. Network analyses show that the mtDNA genome clusters within the southeastern lineage, as expected for an Alpine dingo. Comparison of regulatory regions identified two differentially methylated regions within glucagon receptor GCGR and histone deacetylase HDAC4 genes that are unmethylated in the Alpine dingo genome but hypermethylated in the Desert dingo. Morphological data, comprising geometric morphometric assessment of cranial morphology place dingo Cooinda within population-level variation for Alpine dingoes. Magnetic resonance imaging of brain tissue show she had a larger cranial capacity than a similar-sized domestic dog.

Conclusions

These combined data support the hypothesis that the dingo Cooinda fits the spectrum of genetic and morphological characteristics typical of the Alpine ecotype. We propose that she be considered the archetype specimen for future research investigating the evolutionary history, morphology, physiology, and ecology of dingoes. The female has been taxidermically prepared and is now at the Australian Museum, Sydney.

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.

DNA methylation clocks for dogs and humans

DNA methylation profiles have been used to develop biomarkers of aging known as epigenetic clocks, which predict chronological age with remarkable accuracy and show promise for inferring health status as an indicator of biological age. Epigenetic clocks were first built to monitor human aging, but their underlying principles appear to be evolutionarily conserved, as they have now been successfully developed for many mammalian species. Here, we describe reliable and highly accurate epigenetic clocks shown to apply to 93 domestic dog breeds. The methylation profiles were generated using the mammalian methylation array, which utilizes DNA sequences that are conserved across all mammalian species. Canine epigenetic clocks were constructed to estimate age and also average time to death. We also present two highly accurate human–dog dual species epigenetic clocks (R = 0.97), which may facilitate the ready translation from canine to human use (or vice versa) of antiaging treatments being developed for longevity and preventive medicine. Finally, epigenome-wide association studies here reveal individual methylation sites that may underlie the inverse relationship between breed weight and lifespan. Overall, we describe robust biomarkers to measure aging and, potentially, health status in canines.

Ancestry-inclusive dog genomics challenges popular breed stereotypes

Behavioral genetics in dogs has focused on modern breeds, which are isolated subgroups with distinctive physical and, purportedly, behavioral characteristics. We interrogated breed stereotypes by surveying owners of 18,385 purebred and mixed-breed dogs and genotyping 2155 dogs. Most behavioral traits are heritable [heritability (h²) > 25%], and admixture patterns in mixed-breed dogs reveal breed propensities. Breed explains just 9% of behavioral variation in individuals. Genome-wide association analyses identify 11 loci that are significantly associated with behavior, and characteristic breed behaviors exhibit genetic complexity. Behavioral loci are not unusually differentiated in breeds, but breed propensities align, albeit weakly, with ancestral function. We propose that behaviors perceived as characteristic of modern breeds derive from thousands of years of polygenic adaptation that predates breed formation, with modern breeds distinguished primarily by aesthetic traits.

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.

Forensic DNA phenotyping: Canis familiaris breed classification and skeletal phenotype prediction using functionally significant skeletal SNPs and indels

This review highlights a novel application of breed identification and prediction of skeletal traits in forensic investigations using canine DNA evidence. Currently, genotyping methods used for canine breed classification involve the application of highly polymorphic short tandem repeats in addition to larger commercially available SNP arrays. Both applications face technical challenges. An additional approach to breed identification could be through genotyping SNPs and indels that characterise the array of skeletal differences displayed across domestic dog populations. Research has shown that a small number of genetic variants of large effect drive differences in skeletal phenotypes among domestic dog breeds. This feature makes functionally significant canine skeletal variants a cost-effective target for forensic investigators to classify individuals according to their breed. Further analysis of these skeletal variants would enable the prediction of external appearance. To date, functionally significant genes with genetic variants associated with differences in size, bulk, skull shape, ear shape, limb length, digit type, and tail morphology have been uncovered. Recommendations of a cost-effective genotyping method that can be readily designed and applied by forensic investigators have been given. Further advances to improve the field of canine skeletal forensic DNA phenotyping include the refinement of phenotyping methods, further biological validation of the skeletal genetic variants and establishing a publicly available database for storage of allele frequencies of the skeletal genetic variants in the wider domestic dog population.

The effect of inbreeding, body size and morphology on health in dog breeds

Background

Dog breeds are known for their distinctive body shape, size, coat color, head type and behaviors, features that are relatively similar across members of a breed. Unfortunately, dog breeds are also characterized by distinct predispositions to disease. We explored the relationships between inbreeding, morphology and health using genotype based inbreeding estimates, body weight and insurance data for morbidity.

Results

The average inbreeding based on genotype across 227 breeds was F_adj = 0.249 (95% CI 0.235–0.263). There were significant differences in morbidity between breeds with low and high inbreeding (H = 16.49, P = 0.0004). There was also a significant difference in morbidity between brachycephalic breeds and non-brachycephalic breeds (P = 0.0048) and between functionally distinct groups of breeds (H = 14.95 P < 0.0001). Morbidity was modeled using robust regression analysis and both body weight (P < 0.0001) and inbreeding (P = 0.013) were significant (r² = 0.77). Smaller less inbred breeds were healthier than larger more inbred breeds.

Conclusions

In this study, body size and inbreeding along with deleterious morphologies contributed to increases in necessary health care in dogs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40575-021-00111-4.

Insight into the Candidate Genes and Enriched Pathways Associated with Height, Length, Length to Height Ratio and Body-Weight of Korean Indigenous Breed, Jindo Dog Using Gene Set Enrichment-Based GWAS Analysis

As a companion and hunting dog, height, length, length to height ratio (LHR) and body-weight are the vital economic traits for Jindo dog. Human selection and targeted breeding have produced an extraordinary diversity in these traits. Therefore, the identification of causative markers, genes and pathways that help us to understand the genetic basis of this variability is essential for their selection purposes. Here, we performed a genome-wide association study (GWAS) combined with enrichment analysis on 757 dogs using 118,879 SNPs. The genomic heritability (h²) was 0.33 for height and 0.28 for weight trait in Jindo. At p-value < 5 × 10^-5, ten, six, thirteen and eleven SNPs on different chromosomes were significantly associated with height, length, LHR and body-weight traits, respectively. Based on our results, HHIP, LCORL and NCAPG for height, IGFI and FGFR3 for length, DLK1 and EFEMP1 for LHR and PTPN2, IGFI and RASAL2 for weight can be the potential candidate genes because of the significant SNPs located in their intronic or upstream regions. The gene-set enrichment analysis highlighted here nine and seven overlapping significant (p < 0.05) gene ontology (GO) terms and pathways among traits. Interestingly, the highlighted pathways were related to hormone synthesis, secretion and signalling were generally involved in the metabolism, growth and development process. Our data provide an insight into the significant genes and pathways if verified further, which will have a significant effect on the breeding of the Jindo dog's population.

Whole Genome Sequencing Reveals Signatures for Artificial Selection for Different Sizes in Japanese Primitive Dog Breeds

Body size is an important trait in companion animals. Recently, a primitive Japanese dog breed, the Shiba Inu, has experienced artificial selection for smaller body size, resulting in the "Mame Shiba Inu" breed. To identify loci and genes that might explain the difference in the body size of these Shiba Inu dogs, we applied whole genome sequencing of pooled samples (pool-seq) on both Shiba Inu and Mame Shiba Inu. We identified a total of 13,618,261 unique SNPs in the genomes of these two breeds of dog. Using selective sweep approaches, including F _ST, H _p and XP-CLR with sliding windows, we identified a total of 12 genomic windows that show signatures of selection that overlap with nine genes (PRDM16, ZNF382, ZNF461, ERGIC2, ENSCAFG00000033351, CCDC61, ALDH3A2, ENSCAFG00000011141, and ENSCAFG00000018533). These results provide candidate genes and specific sites that might be associated with body size in dogs. Some of these genes are associated with body size in other mammals, but 8 of the 9 genes are novel candidate genes that need further study.

Introgression contributes to distribution of structural variations in cattle

Structural variations (SVs) are an important source of phenotypic diversity in cattle. Here, 72 whole genome sequences representing taurine and zebu cattle were used to identify SVs. Applying multiple approaches, 16,738 SVs were identified. A comparison against the Database of Genomic Variants archives revealed that 1575 SVs were novel in our data. A novel duplication covering the entire GALNT15 gene, was observed only in N'Dama. A duplication, which was previously reported only in zebu and associated with navel length, was also observed in N'Dama. Investigation of a novel deletion located upstream of CAST13 gene and identified only in Italian cattle and zebu, revealed its introgressed origin in the former. Overall, our data highlights how the SVs distribution in cattle is also shaped by forces such as demographical differences and gene flow. The cattle SVs of this study and its meta-data can be visualized on an interactive genome browser at https://tinyurl.com/svCowArs.

Variants That Differentiate Wolf and Dog Populations Are Enriched in Regulatory Elements

Research on the genetics of domestication most often focuses on the protein-coding exons. However, exons cover only a minor part (1–2%) of the canine genome, whereas functional mutations may be located also in regions beyond the exome, in regulatory regions. Therefore, a large proportion of phenotypical differences between dogs and wolves may remain genetically unexplained. In this study, we identified variants that have high allelic frequency differences (i.e., highly differentiated variants) between wolves and dogs across the canine genome and investigated the potential functionality. We found that the enrichment of highly differentiated variants was substantially higher in promoters than in exons and that such variants were enriched also in enhancers. Several enriched pathways were identified including oxytocin signaling, carbohydrate digestion and absorption, cancer risk, and facial and body features, many of which reflect phenotypes of potential importance during domestication, including phenotypes of the domestication syndrome. The results highlight the importance of regulatory mutations during dog domestication and motivate the functional annotation of the noncoding part of the canine genome.

Identification of Major Loci and Candidate Genes for Meat Production-Related Traits in Broilers

Background

Carcass traits are crucial characteristics of broilers. However, the underlying genetic mechanisms are not well understood. In the current study, significant loci and major-effect candidate genes affecting nine carcass traits related to meat production were analyzed in 873 purebred broilers using an imputation-based genome-wide association study.

Results

The heritability estimates of nine carcass traits, including carcass weight, thigh muscle weight, and thigh muscle percentage, were moderate to high and ranged from 0.21 to 0.39. Twelve genome-wide significant SNPs and 118 suggestively significant SNPs of 546,656 autosomal variants were associated with carcass traits. All SNPs for six weight traits (body weight at 42 days of age, carcass weight, eviscerated weight, whole thigh weight, thigh weight, and thigh muscle weight) were clustered around the 24.08 Kb region (GGA24: 5.73–5.75 Mb) and contained only one candidate gene (DRD2). The most significant SNP, rs15226023, accounted for 4.85–7.71% of the estimated genetic variance of the six weight traits. The remaining SNPs for carcass composition traits (whole thigh percentage and thigh percentage) were clustered around the 42.52 Kb region (GGA3: 53.03–53.08 Mb) and contained only one candidate gene (ADGRG6). The most significant SNP in this region, rs13571431, accounted for 11.89–13.56% of the estimated genetic variance of two carcass composition traits. Some degree of genetic differentiation in ADGRG6 between large and small breeds was observed.

Conclusion

We identified one 24.08 Kb region for weight traits and one 42.52 Kb region for thigh-related carcass traits. DRD2 was the major-effect candidate gene for weight traits, and ADGRG6 was the major-effect candidate gene for carcass composition traits. Our results supply essential information for causative mutation identification of carcass traits in broilers.

Whole genome variants across 57 pig breeds enable comprehensive identification of genetic signatures that underlie breed features

Background

A large number of pig breeds are distributed around the world, their features and characteristics vary among breeds, and they are valuable resources. Understanding the underlying genetic mechanisms that explain across-breed variation can help breeders develop improved pig breeds.

Results

In this study, we performed GWAS using a standard mixed linear model with three types of genome variants (SNP, InDel, and CNV) that were identified from public, whole-genome, sequencing data sets. We used 469 pigs of 57 breeds, and we identified and analyzed approximately 19 million SNPs, 1.8 million InDels, and 18,016 CNVs. We defined six biological phenotypes by the characteristics of breed features to identify the associated genome variants and candidate genes, which included coat color, ear shape, gradient zone, body weight, body length, and body height. A total of 37 candidate genes was identified, which included 27 that were reported previously (e.g., PLAG1 for body weight), but the other 10 were newly detected candidate genes (e.g., ADAMTS9 for coat color).

Conclusion

Our study indicated that using GWAS across a modest number of breeds with high density genome variants provided efficient mapping of complex traits.

Supplementary Information

Supplementary information accompanies this paper at 10.1186/s40104-020-00520-8.

The Genetic Basis of Obesity and Related Metabolic Diseases in Humans and Companion Animals

Obesity is one of the most prevalent health conditions in humans and companion animals globally. It is associated with premature mortality, metabolic dysfunction, and multiple health conditions across species. Obesity is, therefore, of importance in the fields of medicine and veterinary medicine. The regulation of adiposity is a homeostatic process vulnerable to disruption by a multitude of genetic and environmental factors. It is well established that the heritability of obesity is high in humans and laboratory animals, with ample evidence that the same is true in companion animals. In this review, we provide an overview of how genes link to obesity in humans, drawing on a wealth of information from laboratory animal models, and summarise the mechanisms by which obesity causes related disease. Throughout, we focus on how large-scale human studies and niche investigations of rare mutations in severely affected patients have improved our understanding of obesity biology and can inform our ability to interpret results of animal studies. For dogs, cats, and horses, we compare the similarities in obesity pathophysiology to humans and review the genetic studies that have been previously reported in those species. Finally, we discuss how veterinary genetics may learn from humans about studying precise, nuanced phenotypes and implementing large-scale studies, but also how veterinary studies may be able to look past clinical findings to mechanistic ones and demonstrate translational benefits to human research.

Identification of Candidate Genes and Pathways Associated with Obesity-Related Traits in Canines via Gene-Set Enrichment and Pathway-Based GWAS Analysis

The present study aimed to identify causative loci and genes enriched in pathways associated with canine obesity using a genome-wide association study (GWAS). The GWAS was first performed to identify candidate single-nucleotide polymorphisms (SNPs) associated with obesity and obesity-related traits including body weight and blood sugar in 18 different breeds of 153 dogs. A total of 10 and 2 SNPs were found to be significantly (p < 3.74 × 10^-7) associated with body weight and blood sugar, respectively. None of the SNPs were identified to be significantly associated with obesity trait. We subsequently followed up the GWAS analysis with gene-set enrichment and pathway analyses. A gene-set with 1057, 1409, and 1243 SNPs annotated to 449, 933 and 820 genes for obesity, body weight, and blood sugar, respectively was created by sub-setting the GWAS result at a threshold of p < 0.01 for the gene-set enrichment analysis. In total, 84 GO and 21 KEGG pathways for obesity, 114 GO and 44 KEGG pathways for blood sugar, 120 GO and 24 KEGG pathways for body weight were found to be enriched. Among the pathways and GO terms, we highlighted five enriched pathways (Wnt signaling pathway, adherens junction, pathways in cancer, axon guidance, and insulin secretion) and seven GO terms (fat cell differentiation, calcium ion binding, cytoplasm, nucleus, phospholipid transport, central nervous system development, and cell surface) that were found to be shared among all the traits. Our data provide insights into the genes and pathways associated with obesity and obesity-related traits.

A large deletion on CFA28 omitting ACSL5 gene is associated with intestinal lipid malabsorption in the Australian Kelpie dog breed

Inborn errors of metabolism are genetic conditions that can disrupt intermediary metabolic pathways and cause defective absorption and metabolism of dietary nutrients. In an Australian Kelpie breeding population, 17 puppies presented with intestinal lipid malabsorption. Juvenile dogs exhibited stunted postnatal growth, steatorrhea, abdominal distension and a wiry coat. Using genome-wide association analysis, an associated locus on CFA28 (P_raw = 2.87E^-06) was discovered and validated in a closely related population (P_raw = 1.75E^-45). A 103.3 kb deletion NC_006610.3CFA28:g.23380074_23483377del, containing genes Acyl-CoA Synthetase Long Chain Family Member 5 (ACSL5) and Zinc Finger DHHC-Type Containing 6 (ZDHHC6), was characterised using whole transcriptomic data. Whole transcriptomic sequencing revealed no expression of ACSL5 and disrupted splicing of ZDHHC6 in jejunal tissue of affected Kelpies. The ACSL5 gene plays a key role in long chain fatty acid absorption, a phenotype similar to that of our affected Kelpies has been observed in a knockout mouse model. A PCR-based diagnostic test was developed and confirmed fully penetrant autosomal recessive mode of inheritance. We conclude the structural variant causing a deletion of the ACSL5 gene is the most likely cause for intestinal lipid malabsorption in the Australian Kelpie.

Deciphering the puzzles of dog domestication

The domestic dog, as a highly successful domestication model, is well known as a favored human companion. Exploring its domestication history should provide great insight into our understanding of the prehistoric development of human culture and productivity. Furthermore, investigation on the mechanisms underpinning the morphological and behavioral traits associated with canid domestication syndrome is of significance not only for scientific study but also for human medical research. Current development of a multidisciplinary canine genome database, which includes enormous omics data, has substantially improved our understanding of the genetic makeup of dogs. Here, we reviewed recent advances associated with the original history and genetic basis underlying environmental adaptations and phenotypic diversities in domestic dogs, which should provide perspectives on improving the communicative relationship between dogs and humans.

Whole genome sequencing for the investigation of canine mammary tumor inheritance - an initial assessment of high-risk breast cancer genes reveal BRCA2 and STK11 variants potentially associated with risk in purebred dogs

Background

Although, in general, cancer is considered a multifactorial disease, clustering of particular cancers in pedigrees suggests a genetic predisposition and could explain why some dog breeds appear to have an increased risk of certain cancers. To our knowledge, there have been no published reports of whole genome sequencing to investigate inherited canine mammary tumor (CMT) risk, and with little known about CMT genetic susceptibility, we carried out whole genome sequencing on 14 purebred dogs diagnosed with mammary tumors from four breed-specific pedigrees. Following sequencing, each dog’s data was processed through a bioinformatics pipeline. This initial report highlights variants in orthologs of human breast cancer susceptibility genes.

Results

The overall whole genome and exome coverage averages were 26.0X and 25.6X, respectively, with 96.1% of the genome and 96.7% of the exome covered at least 10X. Of the average 7.9 million variants per dog, initial analyses involved surveying variants in orthologs of human breast cancer susceptibility genes, BRCA1, BRCA2, CDH1, PTEN, STK11, and TP53, and identified 19 unique coding variants that were validated through PCR and Sanger sequencing. Statistical analyses identified variants in BRCA2 and STK11 that appear to be associated with CMT, and breed-specific analyses revealed the breeds at the highest risk. Several additional BRCA2 variants showed trends toward significance, but have conflicting interpretations of pathogenicity, and correspond to variants of unknown significance in humans, which require further investigation. Variants in other genes were noted but did not appear to be associated with disease.

Conclusions

Whole genome sequencing proves to be an effective method to elucidate risk of CMT. Risk variants in orthologs of human breast cancer susceptibility genes have been identified. Ultimately, these whole genome sequencing efforts have provided a plethora of data that can also be assessed for novel discovery and have the potential to lead to breakthroughs in canine and human research through comparative analyses.

Genetic Variants Affecting Skeletal Morphology in Domestic Dogs

Purebred dog breeds provide a powerful resource for the discovery of genetic variants affecting skeletal morphology. Domesticated and subsequently purebred dogs have undergone strong artificial selection for a broad range of skeletal variation, which include both the size and shapes of their bones. While the phenotypic variation between breeds is high, within-breed morphological variation is typically low. Approaches for defining genetic variants associated with canine morphology include quantitative within-breed analyses, as well as across-breed analyses, using breed standards as proxies for individual measurements. The ability to identify variants across the genomes of individual dogs can now be paired with precise measures of morphological variation to define the genetic interactions and the phenotypic effect of variants on skeletal morphology.

X-linked duplication copy number variation in a familial overgrowth condition

We describe an overgrowth condition associated with X-linked copy number variation. Three brothers displayed an overgrowth pattern at birth that continued postnatally. Clinical findings included macrocephaly, distinctive facial features, developmental delay and variable clubfoot. Normal fetal growth was noted until the third trimester by Hadlock standards, revealing a late gestational overgrowth pattern. Microarray analysis in the family showed a maternally inherited 680 kb copy number duplication at Xq26.1-q26.2 in all three brothers. Molecular sequencing for known overgrowth conditions including GPC3, Sotos 1 (NSD1), Malan (NFIX), Perlman (DIS3L2), Weaver (EZH2), Opitz-Kaveggia (MED12) loci were negative. BWS IC1 and IC2 methylation and CDKN1C testing was also negative. Normal IGF1 levels excluded X-linked acrogiantism. The duplicated region Xq26.1-q26.2 contained IGSF1 and at least part of the lncRNA FIRRE. IGSF1, a highly expressed pituitary immunoglobulin superfamily gene, was recently implicated in a genome-wide association study of canine size. IGSF1 variants were associated with large canine breeds compared to smaller breeds. Our findings support the hypothesis that an X-linked variant encompassing the IGSF1 region may be associated with body size. Although IGSF1 loss has been noted in human hypothyroidism, this is the first reported phenotype in a family with copy number duplication in the region. Our findings suggest that prenatal evaluation, cross-species evaluation, Mendelian, and GWAS studies may describe a distinctive familial condition and its corresponding phenotypic features.

Genetic Pathways of Aging and Their Relevance in the Dog as a Natural Model of Human Aging

Aging research has experienced a burst of scientific efforts in the last decades as the growing ratio of elderly people has begun to pose an increased burden on the healthcare and pension systems of developed countries. Although many breakthroughs have been reported in understanding the cellular mechanisms of aging, the intrinsic and extrinsic factors that contribute to senescence on higher biological levels are still barely understood. The dog, Canis familiaris, has already served as a valuable model of human physiology and disease. The possible role the dog could play in aging research is still an open question, although utilization of dogs may hold great promises as they naturally develop age-related cognitive decline, with behavioral and histological characteristics very similar to those of humans. In this regard, family dogs may possess unmatched potentials as models for investigations on the complex interactions between environmental, behavioral, and genetic factors that determine the course of aging. In this review, we summarize the known genetic pathways in aging and their relevance in dogs, putting emphasis on the yet barely described nature of certain aging pathways in canines. Reasons for highlighting the dog as a future aging and gerontology model are also discussed, ranging from its unique evolutionary path shared with humans, its social skills, and the fact that family dogs live together with their owners, and are being exposed to the same environmental effects.

Genetic mapping of distal femoral, stifle, and tibial radiographic morphology in dogs with cranial cruciate ligament disease

Cranial cruciate ligament disease (CCLD) is a complex trait. Ten measurements were made on orthogonal distal pelvic limb radiographs of 161 pure and mixed breed dogs with, and 55 without, cranial cruciate partial or complete ligament rupture. Dogs with CCLD had significantly smaller infrapatellar fat pad width, higher average tibial plateau angle, and were heavier than control dogs. The first PC weightings captured the overall size of the dog’s stifle and PC2 weightings reflected an increasing tibial plateau angle coupled with a smaller fat pad width. Of these dogs, 175 were genotyped, and 144,509 polymorphisms were used in a genome-wide association study with both a mixed linear and a multi-locus model. For both models, significant (p_genome <3.46×10⁻⁷ for the mixed and< 6.9x10^-8 for the multilocus model) associations were found for PC1, tibial diaphyseal length and width, fat pad base length, and femoral and tibial condyle width at LCORL, a known body size-regulating locus. Other body size loci with significant associations were growth hormone 1 (GH1), which was associated with the length of the fat pad base and the width of the tibial diaphysis, and a region on CFAX near IRS4 and ACSL4 in the multilocus model. The tibial plateau angle was associated significantly with a locus on CFA10 in the linear mixed model with nearest candidate genes BET1 and MYH9 and on CFA08 near candidate genes WDHD1 and GCH1. MYH9 has a major role in osteoclastogenesis. Our study indicated that tibial plateau slope is associated with CCLD and a compressed infrapatellar fat pad, a surrogate for stifle osteoarthritis. Because of the association between tibial plateau slope and CCLD, and pending independent validation, these candidate genes for tibial plateau slope may be tested in breeds susceptible to CCLD before they develop disease or are bred.

Imputation of canine genotype array data using 365 whole-genome sequences improves power of genome-wide association studies

Genomic resources for the domestic dog have improved with the widespread adoption of a 173k SNP array platform and updated reference genome. SNP arrays of this density are sufficient for detecting genetic associations within breeds but are underpowered for finding associations across multiple breeds or in mixed-breed dogs, where linkage disequilibrium rapidly decays between markers, even though such studies would hold particular promise for mapping complex diseases and traits. Here we introduce an imputation reference panel, consisting of 365 diverse, whole-genome sequenced dogs and wolves, which increases the number of markers that can be queried in genome-wide association studies approximately 130-fold. Using previously genotyped dogs, we show the utility of this reference panel in identifying potentially novel associations, including a locus on CFA20 significantly associated with cranial cruciate ligament disease, and fine-mapping for canine body size and blood phenotypes, even when causal loci are not in strong linkage disequilibrium with any single array marker. This reference panel resource will improve future genome-wide association studies for canine complex diseases and other phenotypes.

Complex traits are controlled by more than one gene and as such are difficult to map. For complex trait mapping in the domestic dog, researchers use the current array of 173,000 variants, with only minimal success. Here, we use a method called imputation to increase the number of variants–from 173,000 to 24 million–that can be queried in canine association studies. We use sequence data from the whole genomes of 365 dogs and wolves to accurately predict variants, in a separate cohort of dogs, that are not present on the array. Using dog body size, blood phenotypes, and a common orthopedic disease that involves rupture of the cranial cruciate ligament, we show that the increase in variants results in an increase in mapping power, through the identification of new associations and the narrowing of regions of interest. This imputation panel is particularly important because of its usefulness in improving complex trait mapping in the dog, which has significant implications for discovery of variants in humans with similar diseases.

Network analysis of canine brain morphometry links tumour risk to oestrogen deficiency and accelerated brain ageing

Structural ‘brain age’ is a valuable but complex biomarker for several brain disorders. The dog is an unrivalled comparator for neurological disease modeling, however canine brain morphometric diversity creates computational and statistical challenges. Using a data-driven approach, we explored complex interactions between patient metadata, brain morphometry, and neurological disease. Twenty-four morphometric parameters measured from 286 canine brain magnetic resonance imaging scans were combined with clinical parameters to generate 9,438 data points. Network analysis was used to cluster patients according to their brain morphometry profiles. An ‘aged-brain’ profile, defined by a small brain width and volume combined with ventriculomegaly, was revealed in the Boxer breed. Key features of this profile were paralleled in neutered female dogs which, relative to un-neutered females, had an 11-fold greater risk of developing brain tumours. Boxer dog and geriatric dog groups were both enriched for brain tumour diagnoses, despite a lack of geriatric Boxers within the cohort. Our findings suggest that advanced brain ageing enhances brain tumour risk in dogs and may be influenced by oestrogen deficiency—a risk factor for dementia and brain tumours in humans. Morphometric features of brain ageing in dogs, like humans, might better predict neurological disease risk than patient chronological age.

Dog10K: an international sequencing effort to advance studies of canine domestication, phenotypes and health

Dogs are the most phenotypically diverse mammalian species, and they possess more known heritable disorders than any other non-human mammal. Efforts to catalog and characterize genetic variation across well-chosen populations of canines are necessary to advance our understanding of their evolutionary history and genetic architecture. To date, no organized effort has been undertaken to sequence the world's canid populations. The Dog10K Consortium (http://www.dog10kgenomes.org) is an international collaboration of researchers from across the globe who will generate 20× whole genomes from 10 000 canids in 5 years. This effort will capture the genetic diversity that underlies the phenotypic and geographical variability of modern canids worldwide. Breeds, village dogs, niche populations and extended pedigrees are currently being sequenced, and de novo assemblies of multiple canids are being constructed. This unprecedented dataset will address the genetic underpinnings of domestication, breed formation, aging, behavior and morphological variation. More generally, this effort will advance our understanding of human and canine health.

Genome Sequence Analysis Reveals Selection Signatures in Endangered Trypanotolerant West African Muturu Cattle

Like most West African Bos taurus, the shorthorn Muturu is under threat of replacement or crossbreeding with zebu. Their populations are now reduced to a few hundred breeding individuals and they are considered endangered. So far, the genetic variation and genetic basis of the trypanotolerant Muturu environmental adaptation have not been assessed. Here, we present genome-wide candidate positive selection signatures in Muturu following within-population iHS and between population Rsb signatures of selection analysis. We compared the results in Muturu with the ones obtained in N’Dama, a West African longhorn trypanotolerant taurine, and in two European taurine (Holstein and Jersey). The results reveal candidate signatures of selection regions in Muturu including genes linked to the innate (e.g., TRIM10, TRIM15, TRIM40, and TRIM26) and the adaptive (e.g., JSP.1, BOLA-DQA2, BOLA-DQA5, BOLA-DRB3, and BLA-DQB) immune responses. The most significant regions are identified on BTA 23 at the bovine major histocompatibility complex (MHC) (iHS analysis) and on BTA 12 at genes including a heat tolerance gene (INTS6) (Rsb analysis). Other candidate selected regions include genes related to growth traits/stature (e.g., GHR and GHRHR), coat color (e.g., MITF and KIT), feed efficiency (e.g., ZRANB3 and MAP3K5) and reproduction (e.g., RFX2, SRY, LAP3, and GPX5). Genes under common signatures of selection regions with N’Dama, including for adaptive immunity and heat tolerance, suggest shared mechanisms of adaptation to environmental challenges for these two West African taurine cattle. Interestingly, out of the 242,910 SNPs identified within the candidate selected regions in Muturu, 917 are missense SNPs (0.4%), with an unequal distribution across 273 genes. Fifteen genes including RBBP8, NID1, TEX15, LAMA3, TRIM40, and OR12D3 comprise 220 missense variants, each between 11 and 32. Our results, while providing insights into the candidate genes under selection in Muturu, are paving the way to the identification of genes and their polymorphisms linked to the unique tropical adaptive traits of the West Africa taurine, including trypanotolerance.

Whole genome sequencing of canids reveals genomic regions under selection and variants influencing morphology

Domestic dog breeds are characterized by an unrivaled diversity of morphologic traits and breed-associated behaviors resulting from human selective pressures. To identify the genetic underpinnings of such traits, we analyze 722 canine whole genome sequences (WGS), documenting over 91 million single nucleotide and small indels, creating a large catalog of genomic variation for a companion animal species. We undertake both selective sweep analyses and genome wide association studies (GWAS) inclusive of over 144 modern breeds, 54 wild canids and a hundred village dogs. Our results identify variants of strong impact associated with 16 phenotypes, including body weight variation which, when combined with existing data, explain greater than 90% of body size variation in dogs. We thus demonstrate that GWAS and selection scans performed with WGS are powerful complementary methods for expanding the utility of companion animal systems for the study of mammalian growth and biology.

Mutations in IRS4 are associated with central hypothyroidism

BACKGROUND

Four genetic causes of isolated congenital central hypothyroidism (CeH) have been identified, but many cases remain unexplained. We hypothesised the existence of other genetic causes of CeH with a Mendelian inheritance pattern.

METHODS

We performed exome sequencing in two families with unexplained isolated CeH and subsequently Sanger sequenced unrelated idiopathic CeH cases. We performed clinical and biochemical characterisation of the probands and carriers identified by family screening. We investigated IRS4 mRNA expression in human hypothalamus and pituitary tissue, and measured serum thyroid hormones and Trh and Tshb mRNA expression in hypothalamus and pituitary tissue of Irs4 knockout mice.

RESULTS

We found mutations in the insulin receptor substrate 4 (IRS4) gene in two pairs of brothers with CeH (one nonsense, one frameshift). Sequencing of IRS4 in 12 unrelated CeH cases negative for variants in known genes yielded three frameshift mutations (two novel) in three patients and one male sibling. All male carriers (n=8) had CeH with plasma free thyroxine concentrations below the reference interval. MRI of the hypothalamus and pituitary showed no structural abnormalities (n=12). 24-hour thyroid-stimulating hormone (TSH) secretion profiles in two adult male patients showed decreased basal, pulsatile and total TSH secretion. IRS4 mRNA was expressed in human hypothalamic nuclei, including the paraventricular nucleus, and in the pituitary gland. Female knockout mice showed decreased pituitary Tshb mRNA levels but had unchanged serum thyroid hormone concentrations.

CONCLUSIONS

Mutations in IRS4 are associated with isolated CeH in male carriers. As IRS4 is involved in leptin signalling, the phenotype may be related to disrupted leptin signalling.

Genes Whose Gain or Loss-Of-Function Increases Skeletal Muscle Mass in Mice: A Systematic Literature Review

Skeletal muscle mass differs greatly in mice and humans and this is partially inherited. To identify muscle hypertrophy candidate genes we conducted a systematic review to identify genes whose experimental loss or gain-of-function results in significant skeletal muscle hypertrophy in mice. We found 47 genes that meet our search criteria and cause muscle hypertrophy after gene manipulation. They are from high to small effect size: Ski, Fst, Acvr2b, Akt1, Mstn, Klf10, Rheb, Igf1, Pappa, Ppard, Ikbkb, Fstl3, Atgr1a, Ucn3, Mcu, Junb, Ncor1, Gprasp1, Grb10, Mmp9, Dgkz, Ppargc1a (specifically the Ppargc1a4 isoform), Smad4, Ltbp4, Bmpr1a, Crtc2, Xiap, Dgat1, Thra, Adrb2, Asb15, Cast, Eif2b5, Bdkrb2, Tpt1, Nr3c1, Nr4a1, Gnas, Pld1, Crym, Camkk1, Yap1, Inhba, Tp53inp2, Inhbb, Nol3, Esr1. Knock out, knock down, overexpression or a higher activity of these genes causes overall muscle hypertrophy as measured by an increased muscle weight or cross sectional area. The mean effect sizes range from 5 to 345% depending on the manipulated gene as well as the muscle size variable and muscle investigated. Bioinformatical analyses reveal that Asb15, Klf10, Tpt1 are most highly expressed hypertrophy genes in human skeletal muscle when compared to other tissues. Many of the muscle hypertrophy-regulating genes are involved in transcription and ubiquitination. Especially genes belonging to three signaling pathways are able to induce hypertrophy: (a) Igf1-Akt-mTOR pathway, (b) myostatin-Smad signaling, and (c) the angiotensin-bradykinin signaling pathway. The expression of several muscle hypertrophy-inducing genes and the phosphorylation of their protein products changes after human resistance and high intensity exercise, in maximally stimulated mouse muscle or in overloaded mouse plantaris.