Logo: to the web site of Uppsala University

Current genetic research on obsessive-compulsive disorder (OCD) supports contributions to risk specifically from common single nucleotide variants (SNVs), along with rare coding SNVs and small insertion-deletions (indels). The contribution to OCD risk from rare copy number variants (CNVs), however, has not been formally assessed at a similar scale. Here we describe an analysis of rare CNVs called from genotype array data in 2248 deeply phenotyped OCD cases and 3608 unaffected controls from Sweden and Norway. Cases carry an elevated burden of CNVs ≥30 kb in size (OR = 1.12, P = 1.77 × 10⁻³). The excess rate of these CNVs in cases versus controls was around 0.07 (95% CI 0.02–0.11, P = 2.58 × 10⁻³). This signal was largely driven by CNVs overlapping protein-coding regions (OR = 1.19, P = 3.08 × 10⁻⁴), particularly deletions impacting loss-of-function intolerant genes (pLI >0.995, OR = 4.12, P = 2.54 × 10⁻⁵). We did not identify any specific locus where CNV burden was associated with OCD case status at genome-wide significance, but we noted non-random recurrence of CNV deletions in cases (permutation P = 2.60 × 10⁻³). In cases where sufficient clinical data were available (n = 1612) we found that carriers of neurodevelopmental duplications were more likely to have comorbid autism (P < 0.001), and that carriers of deletions overlapping neurodevelopmental genes had lower treatment response (P = 0.02). The results demonstrate a contribution of rare CNVs to OCD risk, and suggest that studies of rare coding variation in OCD would have increased power to identify risk genes if this class of variation were incorporated into formal tests.

Diffuse large B-cell lymphoma (DLBCL) is an aggressive hematopoietic neoplasm that affects humans as well as dogs. While previous studies on canine DLBCL (cDLBCL) have significantly advanced our understanding of the disease, the majority of this research has relied on whole-exome sequencing, which is limited in its ability to detect copy number aberrations and other genomic changes beyond coding regions. Furthermore, many of these studies lack sufficient clinical follow-up data, making it difficult to draw meaningful associations between genetic variants and patient outcomes. Our study aimed to characterize the mutational landscape of cDLBCL using whole-genome sequencing of matched tumor-normal samples obtained from a cohort of 43 dogs previously enrolled in a clinical trial for which longitudinal follow-up was available. We focused on identifying genes that were significantly or recurrently mutated with coding point mutations, copy number aberrations, and their associations with patient outcomes. We identified 26 recurrently mutated genes, 18 copy number gains, and 8 copy number losses. Consistent with prior studies, the most commonly mutated genes included TRAF3, FBXW7, POT1, TP53, SETD2, DDX3X and TBL1XR1. The most prominent copy number gain occurred on chromosome 13, overlapping key oncogenes such as MYC and KIT, while the most frequent deletion was a focal loss on chromosome 26, encompassing IGL, PRAME, GNAZ, RAB36, RSPH14, and ZNF280B. Notably, our set of recurrently mutated genes was significantly enriched with genes involved in epigenetic regulation. In particular, we identified hotspot mutations in two histone genes, H3C8, and LOC119877878, resulting in H3K27M alterations predicted to dysregulate gene expression. Finally, a survival analysis revealed that H3K27M mutations in H3C8 were associated with increased hazard ratios for progression-free survival. No copy number aberrations were associated with survival. These findings underscore the critical role of epigenetic dysregulation in cDLBCL and affirm the dog as a relevant large animal model for interrogating the biological activity of novel histone-modifying treatment strategies.

Mammalian genomes are biased towards GC bases at third codon positions, likely due to a GC-biased ancestral genome and the selectively neutral recombination-related process of GC-biased gene conversion. The unwanted transcript hypothesis posits that this high GC content at synonymous sites may be beneficial for protecting against spurious transcripts, particularly in species with low effective population sizes. Utilising a 240 placental mammal genome alignment and single-base resolution conservation scores, we interpret sequence conservation at mammalian four-fold degenerate sites in this context and find evidence in support of the unwanted transcript hypothesis, including a strong GC bias, high conservation at sites relating to exon splicing, less human genetic variation at conserved four-fold degenerate sites, and conservation of sites important for epigenetic regulation of developmental genes. Additionally, we show that high conservation of four-fold degenerate sites in essential developmental genes, including homeobox genes, likely relates to the low mutation rates experienced by these genes.

Background:

Genetic diversity in closed populations, such as pedigree dogs, is of concern for maintaining the health and vitality of the population in the face of evolving challenges. Measures of genetic diversity rely upon estimates of homozygosity without consideration of whether the homozygosity is desirable or undesirable or if heterozygosity has a functional impact. Pedigree coefficients of inbreeding have been the classical approach yet they are inadequate unless based upon the entire population.

Methods:

Homozygosity measures based upon pedigree analyses (n = 11,898), SNP array data (n = 244), and whole genome sequencing (n = 23) were compared in the Bearded Collie, as well as a comparison of SNP array data to a pedigree cohort (n = 5042) and a mixed-breed cohort (n = 1171).

Results:

Molecular measures based upon DNA are more informative on an individual's homozygosity levels than pedigree analyses, although SNP coefficients of inbreeding overestimate the level of inbreeding based on the nature of SNP array methodology. Whole genome sequence (WGS) analyses revealed that the heterozygosity observed is generally in variants having neutral or low impact, which would indicate that the variability may not contribute substantially to functional diversity in the population. The majority of high-impact variants were observed in the shortest runs of homozygosity (ROH) reflecting ancestral breeding and domestication practices. As expected, mixed-breed dogs displayed higher measures of genomic diversity than either Bearded Collies or other pedigree dogs as a whole using the current paradigm algorithm models to calculate homozygosity.

Conclusions:

Using typical DNA-based measures reflect only a single individual and not the population thereby failing to account for regions of homozygosity that reflect ancestral breeding, domestication history, breed-defining regions, or regions positively selected for health traits. Incorporating measures of genetic diversity into dog breeding schemes is meritorious. However, until measures of diversity can distinguish between breed-defining homozygosity and homozygosity associated with positive health alleles, the measures to use as selection tools need refinement before their widespread implementation.

Objective: Systemic inflammatory autoimmune diseases (SIADs) such as systemic lupus erythematosus (SLE), primary Sjögren disease (pSS), and idiopathic inflammatory myopathies (myositis) are complex conditions characterized by shared circulating autoantibodies and clinical manifestations, including skin rashes, among others. This study was aimed at elucidating the genetics underlying these common features.

Methods: We performed targeted DNA sequencing of coding and regulatory regions from approximately 1,900 immune-related genes in a large cohort of 2,292 well-characterized Scandinavian patients with SIADs with SLE, pSS, and myositis as well as 1,252 controls. A gene-based functionally weighted genetic score for aggregate testing of all genetic variants, including rare variants, was complemented by in silico functional analyses and in vitro reporter experiments.

Results: Case–control association analysis detected known and potentially novel genetic loci in agreement with previous genetic and transcriptomics findings linked to the SIAD autoimmune background. Intriguingly, case–case comparisons between patient subgroups with and without specific autoantibodies revealed that the subgroups defined by antinuclear antibodies and anti–double-stranded DNA antibodies have unique genetic profiles reflecting their heterogeneity. When focusing on clinical features, we overall showed that dual-specificity phosphatase 1 (DUSP1) protective genetic variants lead to increased gene expression and potentially to anti-inflammatory effects on the SIAD-associated skin phenotype. This is consistent with recent genetic findings on eczema and with the previously reported down-regulation of the MAPK signaling-related gene DUSP1 in other skin disorders.

Conclusion: Together, this suggests common molecular mechanisms potentially underlying overlapping clinical manifestations shared among different disorders and informs clinical heterogeneity, which could be translated to improve disease diagnostic and treatment, also in more generalized disease frameworks.

Precision medicine has the potential to provide more accurate diagnosis, appropriate treatment and timely prevention strategies by considering patients’ biological makeup. However, this cannot be realized without integrating clinical and omics data in a data-sharing framework that achieves large sample sizes. Systems that integrate clinical and genetic data from multiple sources are scarce due to their distinct data types, interoperability, security and data ownership issues. Here we present a secure framework that allows immutable storage, querying and analysis of clinical and genetic data using blockchain technology. Our platform allows clinical and genetic data to be harmonized by combining them under a unified framework. It supports combined genotype–phenotype queries and analysis, gives institutions control of their data and provides immutable user access logs, improving transparency into how and when health information is used. We demonstrate the value of our framework for precision medicine by creating genotype–phenotype cohorts and examining relationships within them. We show that combining data across institutions using our secure platform increases statistical power for rare disease analysis. By offering an integrated, secure and decentralized framework, we aim to enhance reproducibility and encourage broader participation from communities and patients in data sharing.

The complement system plays a central role in the pathogenesis of Systemic Lupus Erythematosus (SLE), but most studies have focused on the classical pathway. Ficolin-3 is the main initiator of the lectin pathway of complement in humans, but its role in systemic autoimmune disease has not been conclusively determined. Here, we combined biochemical and genetic approaches to assess the contribution of ficolin-3 to SLE risk and disease manifestations. Ficolin-3 activity was measured by a functional assay in serum or plasma samples from Swedish SLE patients (n = 786) and controls matched for age and sex (n = 566). Genetic variants in an extended 300 kb genomic region spanning the FCN3 locus were analyzed for their association with ficolin-3 activity and SLE manifestations in a Swedish multicenter cohort (n = 985). Patients with ficolin-3 activity in the highest tertile showed a strong enrichment in an SLE cluster defined by anti-Sm/DNA/nucleosome antibodies (OR 3.0, p < 0.001) and had increased rates of hematological disease (OR 1.4, p = 0.078) and lymphopenia (OR = 1.6, p = 0.039). Genetic variants associated with low ficolin-3 activity mapped to an extended haplotype in high linkage disequilibrium upstream of the FCN3 gene. Patients carrying the lead genetic variant associated with low ficolin-3 activity had a lower frequency of hematological disease (OR 0.67, p = 0.018) and lymphopenia (OR 0.63, p = 0.031) and fewer autoantibodies (p = 0.0019). Loss-of-function variants in the FCN3 gene were not associated with SLE, but four (0.5 %) SLE patients developed acquired ficolin-3 deficiency where ficolin-3 activity in serum was depleted following diagnosis of SLE.

Taken together, our results provide genetic and biochemical evidence that implicate the lectin pathway in hematological SLE manifestations. We also identify lectin pathway activation through ficolin-3 as a factor that contributes to the autoantibody response in SLE.

Background

Pugs commonly present with thoracolumbar myelopathy, also known as pug dog myelopathy (PDM), which is clinically characterised by progressive signs involving the pelvic limbs, no apparent signs of pain and, often, incontinence. In addition to meningeal fibrosis and focal spinal cord destruction, histopathology has confirmed lymphohistiocytic infiltrates in the central nervous system (CNS) in a considerable number of pugs with PDM. Lymphohistiocytic CNS inflammation also characterises necrotising meningoencephalitis (NME) in pugs. This study aimed to investigate the potential contribution of an immunological aetiology to the development of PDM.

Methods

The concentrations of glial fibrillary acidic protein (GFAP) in serum and CSF and of anti-GFAP autoantibodies in CSF were measured with an ELISA. In addition, a commercial test was used for genetic characterisation of the dog leukocyte antigen class II haplotype, which is associated with NME susceptibility.

Results

This study included 87 dogs: 52 PDM pugs, 14 control pugs, four NME pugs and 17 dogs of breeds other than pugs that were investigated for neurological disease (neuro controls). Anti-GFAP autoantibodies were present in 15 of 19 (79%) of the PDM pugs tested versus six of 16 (38%) of the neuro controls tested (p = 0.018). All 18 PDM pugs evaluated had detectable CSF GFAP. Serum GFAP was detected in two of three (67%) of the NME pugs and in two of 11 (18%) of the control pugs but not in any of the 40 tested PDM pugs. Male pugs heterozygous for the NME risk haplotype had an earlier onset of clinical signs (70 months) compared to male pugs without the risk haplotype (78 months) (p = 0.036).

Limitations

The study was limited by the lack of healthy dogs of breeds other than pugs and the small numbers of control pugs and pugs with NME.

Conclusions

The high proportion of PDM pugs with anti-GFAP autoantibodies and high CSF GFAP concentrations provide support for a potential immunological contribution to the development of PDM.

The Nova Scotia Duck Tolling Retriever (NSDTR) is predisposed to immune mediated rheumatic disease (IMRD), steroid-responsive meningitis-arteritis (SRMA) and certain forms of cancer. Cytokines are the main regulators of the immune system. Interleukin 2 is a cytokine involved in activation of T regulatory cells, playing a role in central tolerance and tumor immunity. Interleukin 12 and interleukin 23 share the same subunit, p40, and are both pro-inflammatory cytokines. The aim of this study was to compare levels of IL-2 in healthy NSDTRs to those with cancer or autoimmune disease and to compare levels of IL-12/IL-23p40 in healthy NSDTRs and beagles versus NSDTRs with cancer or autoimmune disease. 62 dogs were included in the analysis of IL-12/IL-23p40; healthy NSDTRs (n = 16), healthy beagles (n = 16), NSDTRs autoimmune (n = 18) and NDSTRs lymphoma/mastocytoma (n = 12) and 68 dogs for IL-2; healthy (n = 20), autoimmune (n = 36) and lymphoma/mastocytoma/adenocarcinoma (n = 12). NSDTRs with autoimmune disease had higher levels of IL-12/IL-23p40 compared to healthy dogs (p = 0.008). NSDTRs with lymphoma also had higher levels of IL-12/IL-23p40 compared to healthy NSDTRs (p = 0.002). There was no difference in levels of IL-2 between healthy and diseased NSDTR. Statistical analysis was performed using Bonferroni corrections for multiple testing. These findings can contribute to the knowledge of autoimmune disease and cancer in dogs.

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 x 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 x 10^-7) and chromosome 31 (p = 2.5 x 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.