A burden of rare copy number variants in obsessive-compulsive disorder.

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-3). The excess rate of these CNVs in cases versus controls was around 0.07 (95% CI 0.02-0.11, P = 2.58 × 10-3). This signal was largely driven by CNVs overlapping protein-coding regions (OR = 1.19, P = 3.08 × 10-4), particularly deletions impacting loss-of-function intolerant genes (pLI >0.995, OR = 4.12, P = 2.54 × 10-5). 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-3). 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.

Unraveling the Genetics of Shared Clinical and Serological Manifestations in Patients With Systemic Inflammatory Autoimmune Diseases.

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.

In sickness and health - a questionnaire based study regarding immune mediated diseases and neoplasia in Swedish Nova Scotia Duck Tolling Retrievers.

BACKGROUND: The Nova Scotia Duck Tolling Retriever (NSDTR) has previously been highlighted as a breed at risk for developing immune mediated diseases and cancer. The immune response is of great importance for the development of neoplastic disease and a dysregulated immune response may predispose to cancer. Two of the commonly seen immune mediated diseases in NSDTRs are immune mediated rheumatic disease (IMRD), which bears similarities to systemic lupus erythematosus (SLE) affecting humans, and steroid-responsive meningitis-arteritis (SRMA), which is a non-infectious inflammation of the meninges and the leptomeningeal vessels. The aim of this survey study was to investigate the lifetime prevalence of immune mediated diseases and tumors among Swedish NSDTRs based on owners' information. The study design was cross-sectional. A questionnaire was sent to 4102 persons who owned or had previously owned a NSDTR. The questions concerned information about the dog and its overall health status as well as specific diseases. RESULTS: The response rate was 30%, including 935 live NSDTRs, corresponding to approximately 20% of the current population registered in Sweden (n = 4564), and 177 dead dogs. The surveyed dogs were spread over different ages and sex and corresponded to the typical demographic profile of the general dog population. Of the 935 individuals that were alive, 28 dogs (3%) were reported as previously diagnosed with IMRD and 33 dogs (3.5%) were reported as previously diagnosed with SRMA, one dog was reported to have been diagnosed with both SRMA and IMRD. There were 129 dogs (14%) reported to have or have had a neoplasia of some kind. For the dead dogs (n = 177), almost 40% of the owners reported neoplasia as the main reason for death/euthanasia. CONCLUSION: This study reports an estimated lifetime prevalence of IMRD and SRMA, in the studied population of Swedish NSDTRs, of 3.0 and 3.5% respectively. In this study, 14% of the living dogs (n = 935) were reported to have a neoplasia of some kind and almost 40% of the deceased dogs (n = 177) were euthanized due to neoplasia or suspicion of it.

Elevated levels of IL-12/IL-23p40 in Nova Scotia Duck Tolling Retrievers with autoimmune disease and lymphoma.

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.

Hemangiosarcoma Cells Promote Conserved Host-derived Hematopoietic Expansion.

Hemangiosarcoma and angiosarcoma are soft-tissue sarcomas of blood vessel-forming cells in dogs and humans, respectively. These vasoformative sarcomas are aggressive and highly metastatic, with disorganized, irregular blood-filled vascular spaces. Our objective was to define molecular programs which support the niche that enables progression of canine hemangiosarcoma and human angiosarcoma. Dog-in-mouse hemangiosarcoma xenografts recapitulated the vasoformative and highly angiogenic morphology and molecular characteristics of primary tumors. Blood vessels in the tumors were complex and disorganized, and they were lined by both donor and host cells. In a series of xenografts, we observed that the transplanted hemangiosarcoma cells created exuberant myeloid hyperplasia and gave rise to lymphoproliferative tumors of mouse origin. Our functional analyses indicate that hemangiosarcoma cells generate a microenvironment that supports expansion and differentiation of hematopoietic progenitor populations. Furthermore, gene expression profiling data revealed hemangiosarcoma cells expressed a repertoire of hematopoietic cytokines capable of regulating the surrounding stromal cells. We conclude that canine hemangiosarcomas, and possibly human angiosarcomas, maintain molecular properties that provide hematopoietic support and facilitate stromal reactions, suggesting their potential involvement in promoting the growth of hematopoietic tumors. SIGNIFICANCE: We demonstrate that hemangiosarcomas regulate molecular programs supporting hematopoietic expansion and differentiation, providing insights into their potential roles in creating a permissive stromal-immune environment for tumor progression.

Assessment of glial fibrillary acidic protein and anti-glial fibrillary acidic protein autoantibody concentrations and necrotising meningoencephalitis risk genotype in dogs with pug dog myelopathy.

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 HLA region in ANCA-associated vasculitis: characterisation of genetic associations in a Scandinavian patient population.

OBJECTIVE: The antineutrophil cytoplasmic antibody (ANCA)-associated vasculitides (AAV) are inflammatory disorders with ANCA autoantibodies recognising either proteinase 3 (PR3-AAV) or myeloperoxidase (MPO-AAV). PR3-AAV and MPO-AAV have been associated with distinct loci in the human leucocyte antigen (HLA) region. While the association between MPO-AAV and HLA has been well characterised in East Asian populations where MPO-AAV is more common, studies in populations of European descent are limited. The aim of this study was to thoroughly characterise associations to the HLA region in Scandinavian patients with PR3-AAV as well as MPO-AAV. METHODS: Genotypes of single-nucleotide polymorphisms (SNPs) located in the HLA region were extracted from a targeted exome-sequencing dataset comprising Scandinavian AAV cases and controls. Classical HLA alleles were called using xHLA. After quality control, association analyses were performed of a joint SNP/classical HLA allele dataset for cases with PR3-AAV (n=411) and MPO-AAV (n=162) versus controls (n=1595). Disease-associated genetic variants were analysed for association with organ involvement, age at diagnosis and relapse, respectively. RESULTS: PR3-AAV was significantly associated with both HLA-DPB1*04:01 and rs1042335 at the HLA-DPB1 locus, also after stepwise conditional analysis. MPO-AAV was significantly associated with HLA-DRB1*04:04. Neither carriage of HLA-DPB1*04:01 alleles in PR3-AAV nor of HLA-DRB1*04:04 alleles in MPO-AAV were associated with organ involvement, age at diagnosis or relapse. CONCLUSIONS: The association to the HLA region was distinct in Scandinavian cases with MPO-AAV compared with cases of East Asian descent. In PR3-AAV, the two separate signals of association to the HLD-DPB1 region mediate potentially different functional effects.

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.

Genome-wide association study identifies new loci associated with OCD.

To date, four genome-wide association studies (GWAS) of obsessive-compulsive disorder (OCD) have been published, reporting a high single-nucleotide polymorphism (SNP)-heritability of 28% but finding only one significant SNP. A substantial increase in sample size will likely lead to further identification of SNPs, genes, and biological pathways mediating the susceptibility to OCD. We conducted a GWAS meta-analysis with a 2-3-fold increase in case sample size (OCD cases: N = 37,015, controls: N = 948,616) compared to the last OCD GWAS, including six previously published cohorts (OCGAS, IOCDF-GC, IOCDF-GC-trio, NORDiC-nor, NORDiC-swe, and iPSYCH) and unpublished self-report data from 23andMe Inc. We explored the genetic architecture of OCD by conducting gene-based tests, tissue and celltype enrichment analyses, and estimating heritability and genetic correlations with 74 phenotypes. To examine a potential heterogeneity in our data, we conducted multivariable GWASs with MTAG. We found support for 15 independent genome-wide significant loci (14 new) and 79 protein-coding genes. Tissue enrichment analyses implicate multiple cortical regions, the amygdala, and hypothalamus, while cell type analyses yielded 12 cell types linked to OCD (all neurons). The SNP-based heritability of OCD was estimated to be 0.08. Using MTAG we found evidence for specific genetic underpinnings characteristic of different cohort-ascertainment and identified additional significant SNPs. OCD was genetically correlated with 40 disorders or traits-positively with all psychiatric disorders and negatively with BMI, age at first birth and multiple autoimmune diseases. The GWAS meta-analysis identified several biologically informative genes as important contributors to the aetiology of OCD. Overall, we have begun laying the groundwork through which the biology of OCD will be understood and described.

A Burden of Rare Copy Number Variants in Obsessive-Compulsive Disorder.

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 large, 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 2,248 deeply phenotyped OCD cases and 3,608 unaffected controls from Sweden and Norway. We found that in general cases carry an elevated burden of large (>30kb, at least 15 probes) CNVs (OR=1.12, P=1.77×10-3). The excess rate of these CNVs in cases versus controls was around 0.07 (95% CI 0.02-0.11, P=2.58×10-3). This signal was largely driven by CNVs overlapping protein-coding regions (OR=1.19, P=3.08×10-4), particularly deletions impacting loss-of-function intolerant genes (pLI>0.995, OR=4.12, P=2.54×10-5). 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-3). 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 large, 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.

A survey of ficolin-3 activity in Systemic Lupus Erythematosus reveals a link to hematological disease manifestations and autoantibody profile.

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.

Distinct HLA associations with autoantibody-defined subgroups in idiopathic inflammatory myopathies.

BACKGROUND: In patients with idiopathic inflammatory myopathies (IIM), autoantibodies are associated with specific clinical phenotypes suggesting a pathogenic role of adaptive immunity. We explored if autoantibody profiles are associated with specific HLA genetic variants and clinical manifestations in IIM. METHODS: We included 1348 IIM patients and determined the occurrence of 14 myositis-specific or -associated autoantibodies. We used unsupervised cluster analysis to identify autoantibody-defined subgroups and logistic regression to estimate associations with clinical manifestations, HLA-DRB1, HLA-DQA1, HLA-DQB1 alleles, and amino acids imputed from genetic information of HLA class II and I molecules. FINDINGS: We identified eight subgroups with the following dominant autoantibodies: anti-Ro52, -U1RNP, -PM/Scl, -Mi2, -Jo1, -Jo1/Ro52, -TIF1γ or negative for all analysed autoantibodies. Associations with HLA-DRB1∗11, HLA-DRB1∗15, HLA-DQA1∗03, and HLA-DQB1∗03 were present in the anti-U1RNP-dominated subgroup. HLA-DRB1∗03, HLA-DQA1∗05, and HLA-DQB1∗02 alleles were overrepresented in the anti-PM/Scl and anti-Jo1/Ro52-dominated subgroups. HLA-DRB1∗16, HLA-DRB1∗07 alleles were most frequent in anti-Mi2 and HLA-DRB1∗01 and HLA-DRB1∗07 alleles in the anti-TIF1γ subgroup. The HLA-DRB1∗13, HLA-DQA1∗01 and HLA-DQB1∗06 alleles were overrepresented in the negative subgroup. Significant signals from variations in class I molecules were detected in the subgroups dominated by anti-Mi2, anti-Jo1/Ro52, anti-TIF1γ, and the negative subgroup. INTERPRETATION: Distinct HLA class II and I associations were observed for almost all autoantibody-defined subgroups. The associations support autoantibody profiles use for classifying IIM which would likely reflect underlying pathogenic mechanisms better than classifications based on clinical symptoms and/or histopathological features. FUNDING: See a detailed list of funding bodies in the Acknowledgements section at the end of the manuscript.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture.

BACKGROUND: The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function. RESULTS: We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection. CONCLUSIONS: We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Relation between HLA and copy number variation of steroid 21-hydroxylase in a Swedish cohort of patients with autoimmune Addison's disease.

OBJECTIVE: Autoantibodies against the adrenal enzyme 21-hydroxylase is a hallmark manifestation in autoimmune Addison's disease (AAD). Steroid 21-hydroxylase is encoded by CYP21A2, which is located in the human leucocyte antigen (HLA) region together with the highly similar pseudogene CYP21A1P. A high level of copy number variation is seen for the 2 genes, and therefore, we asked whether genetic variation of the CYP21 genes is associated with AAD. DESIGN: Case-control study on patients with AAD and healthy controls. METHODS: Using next-generation DNA sequencing, we estimated the copy number of CYP21A2 and CYP21A1P, together with HLA alleles, in 479 Swedish patients with AAD and autoantibodies against 21-hydroxylase and in 1393 healthy controls. RESULTS: With 95% of individuals carrying 2 functional 21-hydroxylase genes, no difference in CYP21A2 copy number was found when comparing patients and controls. In contrast, we discovered a lower copy number of the pseudogene CYP21A1P among AAD patients (P = 5 × 10-44), together with associations of additional nucleotide variants, in the CYP21 region. However, the strongest association was found for HLA-DQB1*02:01 (P = 9 × 10-63), which, in combination with the DRB1*04:04-DQB1*03:02 haplotype, imposed the greatest risk of AAD. CONCLUSIONS: We identified strong associations between copy number variants in the CYP21 region and risk of AAD, although these associations most likely are due to linkage disequilibrium with disease-associated HLA class II alleles.

Using evolutionary constraint to define novel candidate driver genes in medulloblastoma.

Current knowledge of cancer genomics remains biased against noncoding mutations. To systematically search for regulatory noncoding mutations, we assessed mutations in conserved positions in the genome under the assumption that these are more likely to be functional than mutations in positions with low conservation. To this end, we use whole-genome sequencing data from the International Cancer Genome Consortium and combined it with evolutionary constraint inferred from 240 mammals, to identify genes enriched in noncoding constraint mutations (NCCMs), mutations likely to be regulatory in nature. We compare medulloblastoma (MB), which is malignant, to pilocytic astrocytoma (PA), a primarily benign tumor, and find highly different NCCM frequencies between the two, in agreement with the fact that malignant cancers tend to have more mutations. In PA, a high NCCM frequency only affects the BRAF locus, which is the most commonly mutated gene in PA. In contrast, in MB, >500 genes have high levels of NCCMs. Intriguingly, several loci with NCCMs in MB are associated with different ages of onset, such as the HOXB cluster in young MB patients. In adult patients, NCCMs occurred in, e.g., the WASF-2/AHDC1/FGR locus. One of these NCCMs led to increased expression of the SRC kinase FGR and augmented responsiveness of MB cells to dasatinib, a SRC kinase inhibitor. Our analysis thus points to different molecular pathways in different patient groups. These newly identified putative candidate driver mutations may aid in patient stratification in MB and could be valuable for future selection of personalized treatment options.

Mutations in the CYP27B1 gene cause vitamin D dependent rickets in pugs.

Rickets is a disorder of bone development and can be the result of either dietary or genetic causes. Here, related pugs from 2 litters were included. Three pugs had clinical signs including, lameness, bone deformities, and dyspnea. One other pug was found dead. Radiographs of 2 affected pugs, 5 and 6 months old, showed generalized widening, and irregular margination of the physes of both the appendicular and the axial skeleton with generalized decrease in bone opacity and bulbous swelling of the costochondral junctions. Two pugs had low serum calcium and 1,25 (OH)2 D3 concentrations. Test results further indicated secondary hyperparathyroidism with adequate concentrations of 25-hydroxyvitamin D. Necropsy revealed tongue-like projections of cartilage extending into the metaphysis consistent with rickets, loss of metaphyseal mineralization and lung pathology. Vitamin D-dependent rickets was diagnosed. A truncating mutation in the 1α-hydroxylase gene (CYP27B1) was identified by genome sequence analysis of the pugs with VDDR type 1A. Vitamin D-dependent rickets type 1A can occur in young pugs, and if left untreated is a life-threatening condition. Early medical intervention can reverse clinical signs and should be instituted as soon as possible.

Stratified genetic analysis reveals sex differences in MPO-ANCA-associated vasculitis.

OBJECTIVE: To identify and genetically characterize subgroups of patients with ANCA-associated vasculitides (AAV) based on sex and ANCA subtype. METHODS: A previously established SNP dataset derived from DNA sequencing of 1853 genes and genotyping of 1088 Scandinavian cases with AAV and 1589 controls was stratified for sex and ANCA subtype and analysed for association with five top AAV SNPs. rs9274619, a lead variant at the HLA-DQB1/HLA-DQA2 locus previously associated with AAV positive for myeloperoxidase (MPO)-ANCA, was analysed for association with the cumulative disease involvement of ten different organ systems. RESULTS: rs9274619 showed a significantly stronger association to MPO-ANCA-positive females than males [P = 2.0 × 10-4, OR = 2.3 (95% CI 1.5, 3.5)], whereas proteinase 3 (PR3)-ANCA-associated variants rs1042335, rs9277341 (HLA-DPB1/A1) and rs28929474 (SERPINA1) were equally associated with females and males with PR3-ANCA. In MPO-ANCA-positive cases, carriers of the rs9274619 risk allele were more prone to disease engagement of eyes [P = 0.021, OR = 11 (95% CI 2.2, 205)] but less prone to pulmonary involvement [P = 0.026, OR = 0.52 (95% CI 0.30, 0.92)]. Moreover, AAV with both MPO-ANCA and PR3-ANCA was associated with the PR3-ANCA lead SNP rs1042335 [P = 0.0015, OR = 0.091 (95% CI 0.0022, 0.55)] but not with rs9274619. CONCLUSIONS: Females and males with MPO-ANCA-positive AAV differ in genetic predisposition to disease, suggesting at least partially distinct disease mechanisms between the sexes. Double ANCA-positive AAV cases are genetically similar to PR3-ANCA-positive cases, providing clues to the clinical follow-up and treatment of these patients.

Insights into mammalian TE diversity through the curation of 248 genome assemblies.

We examined transposable element (TE) content of 248 placental mammal genome assemblies, the largest de novo TE curation effort in eukaryotes to date. We found that although mammals resemble one another in total TE content and diversity, they show substantial differences with regard to recent TE accumulation. This includes multiple recent expansion and quiescence events across the mammalian tree. Young TEs, particularly long interspersed elements, drive increases in genome size, whereas DNA transposons are associated with smaller genomes. Mammals tend to accumulate only a few types of TEs at any given time, with one TE type dominating. We also found association between dietary habit and the presence of DNA transposon invasions. These detailed annotations will serve as a benchmark for future comparative TE analyses among placental mammals.

The contribution of historical processes to contemporary extinction risk in placental mammals.

Species persistence can be influenced by the amount, type, and distribution of diversity across the genome, suggesting a potential relationship between historical demography and resilience. In this study, we surveyed genetic variation across single genomes of 240 mammals that compose the Zoonomia alignment to evaluate how historical effective population size (Ne) affects heterozygosity and deleterious genetic load and how these factors may contribute to extinction risk. We find that species with smaller historical Ne carry a proportionally larger burden of deleterious alleles owing to long-term accumulation and fixation of genetic load and have a higher risk of extinction. This suggests that historical demography can inform contemporary resilience. Models that included genomic data were predictive of species' conservation status, suggesting that, in the absence of adequate census or ecological data, genomic information may provide an initial risk assessment.