Sequencing of candidate genes for obesity in Labrador retriever dogs identified a 14 bp deletion in pro-opiomelanocortin (POMC) with an allele frequency of 12%. The deletion disrupts the β-MSH and β-endorphin coding sequences and is associated with body weight (per allele effect of 0.33 SD), adiposity, and greater food motivation. Among other dog breeds, the deletion was only found in the closely related flat-coat retriever (FCR), where it is similarly associated with body weight and food motivation. The mutation is significantly more common in Labrador retrievers selected to become assistance dogs than pets. In conclusion, the deletion in POMC is a significant modifier of weight and appetite in Labrador retrievers and FCRs and may influence other behavioral traits.

Muscle development and lipid accumulation in muscle critically affect meat quality of livestock. However, the genetic factors underlying myofiber-type specification and intramuscular fat (IMF) accumulation remain to be elucidated. Using two independent intercrosses between Western commercial breeds and Korean native pigs (KNPs) and a joint linkage-linkage disequilibrium analysis, we identified a 488.1-kb region on porcine chromosome 12 that affects both reddish meat color (a*) and IMF. In this critical region, only the MYH3 gene, encoding myosin heavy chain 3, was found to be preferentially overexpressed in the skeletal muscle of KNPs. Subsequently, MYH3-transgenic mice demonstrated that this gene controls both myofiber-type specification and adipogenesis in skeletal muscle. We discovered a structural variant in the promotor/regulatory region of MYH3 for which Q allele carriers exhibited significantly higher values of a* and IMF than q allele carriers. Furthermore, chromatin immunoprecipitation and cotransfection assays showed that the structural variant in the 5′-flanking region of MYH3 abrogated the binding of the myogenic regulatory factors (MYF5, MYOD, MYOG, and MRF4). The allele distribution of MYH3 among pig populations worldwide indicated that the MYH3 Q allele is of Asian origin and likely predates domestication. In conclusion, we identified a functional regulatory sequence variant in porcine MYH3 that provides novel insights into the genetic basis of the regulation of myofiber type ratios and associated changes in IMF in pigs. The MYH3 variant can play an important role in improving pork quality in current breeding programs. Local native livestock breeds, such as pigs and cattle, are important reservoirs of genetic variation, because they have adapted to local conditions and have not been as strongly selected as commercial livestock breeds. Therefore, the study of local native breeds is of utmost interest to identify causative mutations for traits in which large differences exist between local native breeds and breeds used internationally for meat production. Herein, we apply these principles to meat quality related traits in pigs. A combined linkage-linkage disequilibrium analysis identified the MYH3 gene on pig chromosome 12, which encodes myosin heavy chain 3 protein. Subsequently, we discovered a regulatory mutation in MYH3, which occurs at high frequency in Korean Native Pigs compared to the frequency at which it occurs in Western pig breeds. We confirmed the effect of MYH3 on muscle fiber-type composition and intramuscular fat content using transgenic mice and transfection assays in porcine fibroblast cells. This is the first causal mutation identified in livestock for these traits and the result is immediately applicable for the genetic improvement of pork quality. These results also enhance our understanding of the genetic basis regarding the regulation of myofiber type ratios and associated changes in intramuscular fat contents in pigs and mammals in general.

Gliomas are the most common form of malignant primary brain tumors in humans and second most common in dogs, occurring with similar frequencies in both species. Dogs are valuable spontaneous models of human complex diseases including cancers and may provide insight into disease susceptibility and oncogenesis. Several brachycephalic breeds such as Boxer, Bulldog and Boston Terrier have an elevated risk of developing glioma, but others, including Pug and Pekingese, are not at higher risk. To identify glioma-associated genetic susceptibility factors, an across-breed genome-wide association study (GWAS) was performed on 39 dog glioma cases and 141 controls from 25 dog breeds, identifying a genome-wide significant locus on canine chromosome (CFA) 26 (p = 2.8 x 10⁻⁸). Targeted re-sequencing of the 3.4 Mb candidate region was performed, followed by genotyping of the 56 SNVs that best fit the association pattern between the re-sequenced cases and controls. We identified three candidate genes that were highly associated with glioma susceptibility: CAMKK2, P2RX7 and DENR. CAMKK2 showed reduced expression in both canine and human brain tumors, and a non-synonymous variant in P2RX7, previously demonstrated to have a 50% decrease in receptor function, was also associated with disease. Thus, one or more of these genes appear to affect glioma susceptibility. Gliomas are devastating malignant brain tumors that are very rarely curable. Despite extensive research to define pathways and genes involved in the development of disease, there is still an urgent need to improve therapy. Some dog breeds have a considerable elevated risk of glioma, making the dog a suitable model for locating genes potentially of importance also for development of human glioma. In this study we defined a genomic region strongly associated with glioma in dogs. We also showed that this genomic region had likely been under selection in the dog breeds with the highest risk of developing glioma. Sometimes selection for breed specific traits results in amplification of disease causing mutations together with the variant selected for. We located three candidate genes in the identified region: CAMKK2, P2RX7 and DENR. We performed further functional studies to evaluate the potential role of these genes in both canine and human glioma. By comparing normal and tumor tissue we could show that two of the genes—CAMKK2 and P2RX7, were affected at the level of gene expression and protein structure, respectively. We propose that further investigation of all three genes could be of interest with potential benefit to both dog and human.

Background: Long terminal repeats (LTRs, consisting of U3-R-U5 portions) are important elements of retroviruses and related retrotransposons. They are difficult to analyse due to their variability. The aim was to obtain a more comprehensive view of structure, diversity and phylogeny of LTRs than hitherto possible. Results: Hidden Markov models (HMM) were created for 11 clades of LTRs belonging to Retroviridae (class III retroviruses), animal Metaviridae (Gypsy/Ty3) elements and plant Pseudoviridae (Copia/Ty1) elements, complementing our work with Orthoretrovirus HMMs. The great variation in LTR length of plant Metaviridae and the few divergent animal Pseudoviridae prevented building HMMs from both of these groups. Animal Metaviridae LTRs had the same conserved motifs as retroviral LTRs, confirming that the two groups are closely related. The conserved motifs were the short inverted repeats (SIRs), integrase recognition signals (5´TGTTRNR…YNYAACA 3´); the polyadenylation signal or AATAAA motif; a GT-rich stretch downstream of the polyadenylation signal; and a less conserved AT-rich stretch corresponding to the core promoter element, the TATA box. Plant Pseudoviridae LTRs differed slightly in having a conserved TATA-box, TATATA, but no conserved polyadenylation signal, plus a much shorter R region. The sensitivity of the HMMs for detection in genomic sequences was around 50% for most models, at a relatively high specificity, suitable for genome screening. The HMMs yielded consensus sequences, which were aligned by creating an HMM model (a ‘Superviterbi’ alignment). This yielded a phylogenetic tree that was compared with a Pol-based tree. Both LTR and Pol trees supported monophyly of retroviruses. In both, Pseudoviridae was ancestral to all other LTR retrotransposons. However, the LTR trees showed the chromovirus portion of Metaviridae clustering together with Pseudoviridae, dividing Metaviridae into two portions with distinct phylogeny. Conclusion: The HMMs clearly demonstrated a unitary conserved structure of LTRs, supporting that they arose once during evolution. We attempted to follow the evolution of LTRs by tracing their functional foundations, that is, acquisition of RNAse H, a combined promoter/ polyadenylation site, integrase, hairpin priming and the primer binding site (PBS). Available information did not support a simple evolutionary chain of events.

Objectives Systemic lupus erythematosus (SLE) is an autoimmune disease with extensive heterogeneity in disease presentation between patients, which is likely due to an underlying molecular diversity. Here, we aimed at elucidating the genetic aetiology of SLE from the immunity pathway level to the single variant level, and stratify patients with SLE into distinguishable molecular subgroups, which could inform treatment choices in SLE. Methods We undertook a pathway-centred approach, using sequencing of immunological pathway genes. Altogether 1832 candidate genes were analysed in 958 Swedish patients with SLE and 1026 healthy individuals. Aggregate and single variant association testing was performed, and we generated pathway polygenic risk scores (PRS). Results We identified two main independent pathways involved in SLE susceptibility: T lymphocyte differentiation and innate immunity, characterised by HLA and interferon, respectively. Pathway PRS defined pathways in individual patients, who on average were positive for seven pathways. We found that SLE organ damage was more pronounced in patients positive for the T or B cell receptor signalling pathways. Further, pathway PRS-based clustering allowed stratification of patients into four groups with different risk score profiles. Studying sets of genes with priors for involvement in SLE, we observed an aggregate common variant contribution to SLE at genes previously reported for monogenic SLE as well as at interferonopathy genes. Conclusions Our results show that pathway risk scores have the potential to stratify patients with SLE beyond clinical manifestations into molecular subsets, which may have implications for clinical follow-up and therapy selection.