Correction: The genetic consequences of dog breed formation-Accumulation of deleterious genetic variation and fixation of mutations associated with myxomatous mitral valve disease in cavalier King Charles spaniels.

[This corrects the article DOI: 10.1371/journal.pgen.1009726.].

Identification and functional characterization of a novel susceptibility locus for small vessel vasculitis with MPO-ANCA.

OBJECTIVE: To identify and characterize genetic loci associated with the risk of developing ANCA-associated vasculitides (AAV). METHODS: Genetic association analyses were performed after Illumina sequencing of 1853 genes and subsequent replication with genotyping of selected single nucleotide polymorphisms in a total cohort of 1110 Scandinavian cases with granulomatosis with polyangiitis or microscopic polyangiitis, and 1589 controls. A novel AAV-associated single nucleotide polymorphism was analysed for allele-specific effects on gene expression using luciferase reporter assay. RESULTS: PR3-ANCA+ AAV was significantly associated with two independent loci in the HLA-DPB1/HLA-DPA1 region [rs1042335, P = 6.3 × 10-61, odds ratio (OR) 0.10; rs9277341, P = 1.5 × 10-44, OR 0.22] and with rs28929474 in the SERPINA1 gene (P = 2.7 × 10-10, OR 2.9). MPO-ANCA+ AAV was significantly associated with the HLA-DQB1/HLA-DQA2 locus (rs9274619, P = 5.4 × 10-25, OR 3.7) and with a rare variant in the BACH2 gene (rs78275221, P = 7.9 × 10-7, OR 3.0), the latter a novel susceptibility locus for MPO-ANCA+ granulomatosis with polyangiitis/microscopic polyangiitis. The rs78275221-A risk allele reduced luciferase gene expression in endothelial cells, specifically, as compared with the non-risk allele. CONCLUSION: We identified a novel susceptibility locus for MPO-ANCA+ AAV and propose that the associated variant is of mechanistic importance, exerting a regulatory function on gene expression in specific cell types.

The genetic consequences of dog breed formation-Accumulation of deleterious genetic variation and fixation of mutations associated with myxomatous mitral valve disease in cavalier King Charles spaniels.

Selective breeding for desirable traits in strictly controlled populations has generated an extraordinary diversity in canine morphology and behaviour, but has also led to loss of genetic variation and random entrapment of disease alleles. As a consequence, specific diseases are now prevalent in certain breeds, but whether the recent breeding practice led to an overall increase in genetic load remains unclear. Here we generate whole genome sequencing (WGS) data from 20 dogs per breed from eight breeds and document a ~10% rise in the number of derived alleles per genome at evolutionarily conserved sites in the heavily bottlenecked cavalier King Charles spaniel breed (cKCs) relative to in most breeds studied here. Our finding represents the first clear indication of a relative increase in levels of deleterious genetic variation in a specific breed, arguing that recent breeding practices probably were associated with an accumulation of genetic load in dogs. We then use the WGS data to identify candidate risk alleles for the most common cause for veterinary care in cKCs-the heart disease myxomatous mitral valve disease (MMVD). We verify a potential link to MMVD for candidate variants near the heart specific NEBL gene in a dachshund population and show that two of the NEBL candidate variants have regulatory potential in heart-derived cell lines and are associated with reduced NEBL isoform nebulette expression in papillary muscle (but not in mitral valve, nor in left ventricular wall). Alleles linked to reduced nebulette expression may hence predispose cKCs and other breeds to MMVD via loss of papillary muscle integrity.

Association of Protective HLA-A With HLA-B∗27 Positive Ankylosing Spondylitis.

OBJECTIVES: To further elucidate the role of the MHC in ankylosing spondylitis by typing 17 genes, searching for HLA-B∗27 independent associations and assessing the impact of sex on this male biased disease. METHODS: High-confidence two-field resolution genotyping was performed on 310 cases and 2196 controls using an n-1 concordance method. Protein-coding variants were called from next-generation sequencing reads using up to four software programs and the consensus result recorded. Logistic regression tests were applied to the dataset as a whole, and also in stratified sets based on sex or HLA-B∗27 status. The amino acids driving association were also examined. RESULTS: Twenty-five HLA protein-coding variants were significantly associated to disease in the population. Three novel protective associations were found in a HLA-B∗27 positive population, HLA-A∗24:02 (OR = 0.4, CI = 0.2-0.7), and HLA-A amino acids Leu95 and Gln156. We identified a key set of seven loci that were common to both sexes, and robust to change in sample size. Stratifying by sex uncovered three novel risk variants restricted to the female population (HLA-DQA1∗04.01, -DQB1∗04:02, -DRB1∗08:01; OR = 2.4-3.1). We also uncovered a set of neutral variants in the female population, which in turn conferred strong effects in the male set, highlighting how population composition can lead to the masking of true associations. CONCLUSION: Population stratification allowed for a nuanced investigation into the tightly linked MHC region, revealing novel HLA-B∗27 signals as well as replicating previous HLA-B∗27 dependent results. This dissection of signals may help to elucidate sex biased disease predisposition and clinical progression.

Targeted sequencing reveals the somatic mutation landscape in a Swedish breast cancer cohort.

Breast cancer (BC) is a genetically heterogeneous disease with high prevalence in Northern Europe. However, there has been no detailed investigation into the Scandinavian somatic landscape. Here, in a homogeneous Swedish cohort, we describe the somatic events underlying BC, leveraging a targeted next-generation sequencing approach. We designed a 20.5 Mb array targeting coding and regulatory regions of genes with a known role in BC (n = 765). The selected genes were either from human BC studies (n = 294) or from within canine mammary tumor associated regions (n = 471). A set of predominantly estrogen receptor positive tumors (ER + 85%) and their normal tissue counterparts (n = 61) were sequenced to ~ 140 × and 85 × mean target coverage, respectively. MuTect2 and VarScan2 were employed to detect single nucleotide variants (SNVs) and copy number aberrations (CNAs), while MutSigCV (SNVs) and GISTIC (CNAs) algorithms estimated the significance of recurrent somatic events. The significantly mutated genes (q ≤ 0.01) were PIK3CA (28% of patients), TP53 (21%) and CDH1 (11%). However, histone modifying genes contained the largest number of variants (KMT2C and ARID1A, together 28%). Mutations in KMT2C were mutually exclusive with PI3KCA mutations (p ≤ 0. 001) and half of these affect the formation of a functional PHD domain. The tumor suppressor CDK10 was deleted in 80% of the cohort while the oncogene MDM4 was amplified. Mutational signature analyses pointed towards APOBEC deaminase activity (COSMIC signature 2) and DNA mismatch repair (COSMIC signature 6). We noticed two significantly distinct patterns related to patient age; TP53 being more mutated in the younger group (29% vs 9% of patients) and CDH23 mutations were absent from the older group. The increased somatic mutation prevalence in the histone modifying genes KMT2C and ARID1A distinguishes the Swedish cohort from previous studies. KMT2C regulates enhancer activation and assists tumor proliferation in a hormone-rich environment, possibly pointing to a role in ER + BC, especially in older cases. Finally, age of onset appears to affect the mutational landscape suggesting that a larger age-diverse population incorporating more molecular subtypes should be studied to elucidate the underlying mechanisms.

Whole-genome genotyping and resequencing reveal the association of a deletion in the complex interferon alpha gene cluster with hypothyroidism in dogs.

BACKGROUND: Hypothyroidism is a common complex endocrinopathy that typically has an autoimmune etiology, and it affects both humans and dogs. Genetic and environmental factors are both known to play important roles in the disease development. In this study, we sought to identify the genetic risk factors potentially involved in the susceptibility to the disease in the high-risk Giant Schnauzer dog breed. RESULTS: By employing genome-wide association followed by fine-mapping (top variant p-value = 5.7 × 10- 6), integrated with whole-genome resequencing and copy number variation analysis, we detected a ~ 8.9 kbp deletion strongly associated (p-value = 0.0001) with protection against development of hypothyroidism. The deletion is located between two predicted Interferon alpha (IFNA) genes and it may eliminate functional elements potentially involved in the transcriptional regulation of these genes. Remarkably, type I IFNs have been extensively associated to human autoimmune hypothyroidism and general autoimmunity. Nonetheless, the extreme genomic complexity of the associated region on CFA11 warrants further long-read sequencing and annotation efforts in order to ascribe functions to the identified deletion and to characterize the canine IFNA gene cluster in more detail. CONCLUSIONS: Our results expand the current knowledge on genetic determinants of canine hypothyroidism by revealing a significant link with the human counterpart disease, potentially translating into better diagnostic tools across species, and may contribute to improved canine breeding strategies.

Multiple regulatory variants located in cell type-specific enhancers within the PKP2 locus form major risk and protective haplotypes for canine atopic dermatitis in German shepherd dogs.

BACKGROUND: Canine atopic dermatitis (CAD) is a chronic inflammatory skin disease triggered by allergic reactions involving IgE antibodies directed towards environmental allergens. We previously identified a ~1.5 Mb locus on canine chromosome 27 associated with CAD in German shepherd dogs (GSDs). Fine-mapping indicated association closest to the PKP2 gene encoding plakophilin 2. RESULTS: Additional genotyping and association analyses in GSDs combined with control dogs from five breeds with low-risk for CAD revealed the top SNP 27:19,086,778 (p = 1.4 × 10(-7)) and a rare ~48 kb risk haplotype overlapping the PKP2 gene and shared only with other high-risk CAD breeds. We selected altogether nine SNPs (four top-associated in GSDs and five within the ~48 kb risk haplotype) that spanned ~280 kb forming one risk haplotype carried by 35 % of the GSD cases and 10 % of the GSD controls (OR = 5.1, p = 5.9 × 10(-5)), and another haplotype present in 85 % of the GSD cases and 98 % of the GSD controls and conferring a protective effect against CAD in GSDs (OR = 0.14, p = 0.0032). Eight of these SNPs were analyzed for transcriptional regulation using reporter assays where all tested regions exerted regulatory effects on transcription in epithelial and/or immune cell lines, and seven SNPs showed allelic differences. The DNA fragment with the top-associated SNP 27:19,086,778 displayed the highest activity in keratinocytes with 11-fold induction of transcription by the risk allele versus 8-fold by the control allele (pdifference = 0.003), and also mapped close (~3 kb) to an ENCODE skin-specific enhancer region. CONCLUSIONS: Our experiments indicate that multiple CAD-associated genetic variants located in cell type-specific enhancers are involved in gene regulation in different cells and tissues. No single causative variant alone, but rather multiple variants combined in a risk haplotype likely contribute to an altered expression of the PKP2 gene, and possibly nearby genes, in immune and epithelial cells, and predispose GSDs to CAD.

Variants within the SP110 nuclear body protein modify risk of canine degenerative myelopathy.

Canine degenerative myelopathy (DM) is a naturally occurring neurodegenerative disease with similarities to some forms of amyotrophic lateral sclerosis (ALS). Most dogs that develop DM are homozygous for a common superoxide dismutase 1 gene (SOD1) mutation. However, not all dogs homozygous for this mutation develop disease. We performed a genome-wide association analysis in the Pembroke Welsh Corgi (PWC) breed comparing DM-affected and -unaffected dogs homozygous for the SOD1 mutation. The analysis revealed a modifier locus on canine chromosome 25. A haplotype within the SP110 nuclear body protein (SP110) was present in 40% of affected compared with 4% of unaffected dogs (P = 1.5 × 10(-5)), and was associated with increased probability of developing DM (P = 4.8 × 10(-6)) and earlier onset of disease (P = 1.7 × 10(-5)). SP110 is a nuclear body protein involved in the regulation of gene transcription. Our findings suggest that variations in SP110-mediated gene transcription may underlie, at least in part, the variability in risk for developing DM among PWCs that are homozygous for the disease-related SOD1 mutation. Further studies are warranted to clarify the effect of this modifier across dog breeds.

Genome-Wide Analysis Identifies Germ-Line Risk Factors Associated with Canine Mammary Tumours.

Canine mammary tumours (CMT) are the most common neoplasia in unspayed female dogs. CMTs are suitable naturally occurring models for human breast cancer and share many characteristics, indicating that the genetic causes could also be shared. We have performed a genome-wide association study (GWAS) in English Springer Spaniel dogs and identified a genome-wide significant locus on chromosome 11 (praw = 5.6x10-7, pperm = 0.019). The most associated haplotype spans a 446 kb region overlapping the CDK5RAP2 gene. The CDK5RAP2 protein has a function in cell cycle regulation and could potentially have an impact on response to chemotherapy treatment. Two additional loci, both on chromosome 27, were nominally associated (praw = 1.97x10-5 and praw = 8.30x10-6). The three loci explain 28.1±10.0% of the phenotypic variation seen in the cohort, whereas the top ten associated regions account for 38.2±10.8% of the risk. Furthermore, the ten GWAS loci and regions with reduced genetic variability are significantly enriched for snoRNAs and tumour-associated antigen genes, suggesting a role for these genes in CMT development. We have identified several candidate genes associated with canine mammary tumours, including CDK5RAP2. Our findings enable further comparative studies to investigate the genes and pathways in human breast cancer patients.

Utilizing the Dog Genome in the Search for Novel Candidate Genes Involved in Glioma Development-Genome Wide Association Mapping followed by Targeted Massive Parallel Sequencing Identifies a Strongly Associated Locus.

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-8). 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.

A Ham1p-dependent mechanism and modulation of the pyrimidine biosynthetic pathway can both confer resistance to 5-fluorouracil in yeast.

5-Fluorouracil (5-FU) is an anticancer drug and pyrimidine analogue. A problem in 5-FU therapy is acquired resistance to the drug. To find out more about the mechanisms of resistance, we screened a plasmid library in yeast for genes that confer 5-FU resistance when overexpressed. We cloned five genes: CPA1, CPA2, HMS1, HAM1 and YJL055W. CPA1 and CPA2 encode a carbamoyl phosphate synthase involved in arginine biosynthesis and HMS1 a helix-loop-helix transcription factor. Our results suggest that CPA1, CPA2, and HMS1 confer 5-FU resistance by stimulating pyrimidine biosynthesis. Thus, they are unable to confer 5-FU resistance in a ura2 mutant, and inhibit the uptake and incorporation into RNA of both uracil and 5-FU. In contrast, HAM1 and YJL055W confer 5-FU resistance in a ura2 mutant, and selectively inhibit incorporation into RNA of 5-FU but not uracil. HAM1 is the strongest resistance gene, but it partially depends on YJL055W for its function. This suggests that HAM1 and YJL055W function together in mediating resistance to 5-FU. Ham1p encodes an inosine triphosphate pyrophosphatase that has been implicated in resistance to purine analogues. Our results suggest that Ham1p could have a broader specificity that includes 5-FUTP and other pyrimidine analogoue triphosphates.

Thorough investigation of a canine autoinflammatory disease (AID) confirms one main risk locus and suggests a modifier locus for amyloidosis.

Autoinflammatory disease (AID) manifests from the dysregulation of the innate immune system and is characterised by systemic and persistent inflammation. Clinical heterogeneity leads to patients presenting with one or a spectrum of phenotypic signs, leading to difficult diagnoses in the absence of a clear genetic cause. We used separate genome-wide SNP analyses to investigate five signs of AID (recurrent fever, arthritis, breed specific secondary dermatitis, otitis and systemic reactive amyloidosis) in a canine comparative model, the pure bred Chinese Shar-Pei. Analysis of 255 DNA samples revealed a shared locus on chromosome 13 spanning two peaks of association. A three-marker haplotype based on the most significant SNP (p<2.6×10(-8)) from each analysis showed that one haplotypic pair (H13-11) was present in the majority of AID individuals, implicating this as a shared risk factor for all phenotypes. We also noted that a genetic signature (F ST) distinguishing the phenotypic extremes of the breed specific Chinese Shar-Pei thick and wrinkled skin, flanked the chromosome 13 AID locus; suggesting that breed development and differentiation has played a parallel role in the genetics of breed fitness. Intriguingly, a potential modifier locus for amyloidosis was revealed on chromosome 14, and an investigation of candidate genes from both this and the chromosome 13 regions revealed significant (p<0.05) renal differential expression in four genes previously implicated in kidney or immune health (AOAH, ELMO1, HAS2 and IL6). These results illustrate that phenotypic heterogeneity need not be a reflection of genetic heterogeneity, and that genetic modifiers of disease could be masked if syndromes were not first considered as individual clinical signs and then as a sum of their component parts.

Rescue and characterization of episomally replicating DNA from the moss Physcomitrella.

The moss Physcomitrella is unique among plants in that it permits efficient gene targeting by homologous recombination. Furthermore, transformed DNA can replicate episomally in Physcomitrella. Here we show that episomally replicating DNA can be rescued back into Escherichia coli, and we use such rescue to study the fate of the transformed DNA. Significantly, plasmids rescued from moss transformed with circular DNA are identical to the original plasmid, whereas plasmids rescued from moss transformed with linearized DNA frequently have deletions created by direct repeat recombination. These events are highly predictable in that they target the longest direct repeat on the plasmid if this repeat is at least 12 bp. Episomal transformants obtained with linearized DNA show a more than 1,000-fold amplification of the DNA whereas transformants obtained with circular DNA have much lower copy numbers. Most episomal transformants quickly lose the plasmid in the absence of selection, but a semistable type of transformant that loses the plasmid at a much lower frequency was also observed. The consistent rescue of the original plasmid, or of predictable derivatives thereof, suggests that molecular genetics methods which rely on shuttle plasmids are feasible in Physcomitrella.

Combinatorial control of gene expression by the three yeast repressors Mig1, Mig2 and Mig3.

BACKGROUND: Expression of a large number of yeast genes is repressed by glucose. The zinc finger protein Mig1 is the main effector in glucose repression, but yeast also has two related proteins: Mig2 and Mig3. We have used microarrays to study global gene expression in all possible combinations of mig1, mig2 and mig3 deletion mutants. RESULTS: Mig1 and Mig2 repress a largely overlapping set of genes on 2% glucose. Genes that are upregulated in a mig1 mig2 double mutant were grouped according to the contribution of Mig2. Most of them show partially redundant repression, with Mig1 being the major repressor, but some genes show complete redundancy, and some are repressed only by Mig1. Several redundantly repressed genes are involved in phosphate metabolism. The promoters of these genes are enriched for Pho4 sites, a novel GGGAGG motif, and a variant Mig1 site which is absent from genes repressed only by Mig1. Genes repressed only by Mig1 on 2% glucose include the hexose transporter gene HXT4, but Mig2 contributes to HXT4 repression on 10% glucose. HXT6 is one of the few genes that are more strongly repressed by Mig2. Mig3 does not seem to overlap in function with Mig1 and Mig2. Instead, Mig3 downregulates the SIR2 gene encoding a histone deacetylase involved in gene silencing and the control of aging. CONCLUSION: Mig2 fine-tunes glucose repression by targeting a subset of the Mig1-repressed genes, and by responding to higher glucose concentrations. Mig3 does not target the same genes as Mig1 and Mig2, but instead downregulates the SIR2 gene.

Functional genomics of monensin sensitivity in yeast: implications for post-Golgi traffic and vacuolar H+-ATPase function.

We have screened a complete collection of yeast knockout mutants for sensitivity to monensin, an ionophore that interferes with intracellular transport. A total of 63 sensitive strains were found. Most of the strains were deleted for genes involved in post-Golgi traffic, with an emphasis on vacuolar biogenesis. A high correlation was thus seen with VPS and VAM genes, but there were also significant differences between the three sets of genes. A weaker correlation was seen with sensitivity to NaCl, in particular rate of growth effects. Interestingly, all 14 genes encoding subunits of the vacuolar H(+)-ATPase (V-ATPase) were absent in our screen, even though they appeared in the VPS or VAM screens. All monensin-sensitive mutants that could be tested interact synthetically with a deletion of the A subunit of the V-ATPase, Vma1. Synthetic lethality was limited to mutations affecting endocytosis or retrograde transport to Golgi. In addition, vma1 was epistatic over the monensin sensitivity of vacuolar transport mutants, but not endocytosis mutants. Deletions of the two isoforms of the V-ATPase a subunit, Vph1 and Stv1 had opposite effects on the monensin sensitivity of a ypt7 mutant. These findings are consistent with a model where monensin inhibits growth by interfering with the maintenance of an acidic pH in the late secretory pathway. The synthetic lethality of vma1 with mutations affecting retrograde transport to the Golgi further suggests that it is in the late Golgi that a low pH must be maintained.

The jmjN and jmjC domains of the yeast zinc finger protein Gis1 interact with 19 proteins involved in transcription, sumoylation and DNA repair.

The jumonji domain is a highly conserved bipartite domain made up of two subdomains, jmjN and jmjC, which is found in many eukaryotic transcription factors. The jmjC domain was recently shown to possess the histone demethylase activity. Here we show that the jmjN and jmjC domains of the yeast zinc finger protein Gis1 interact in a two-hybrid system with 19 yeast proteins that include the RecQ helicase Sgs1, the silencing factors Esc1 and Sir4, the URI-type prefoldin Bud27 and the PIAS type SUMO ligase Nfi1/Siz2. Extensive interaction cross dependencies further suggest that the proteins form a larger complex. Consistent with this, 16 of the proteins also interact with a Bud27 two-hybrid bait, and three of them co-precipitate with TAP-tagged Gis1. The Gis1 jumonji domain can repress transcription when recruited to a promoter as a lexA fusion. This effect is dependent on both the jmjN and jmjC subdomains, as were all 19 two-hybrid interactions, indicating that the two subdomains form a single functional unit. The human Sgs1 homolog WRN also interacts with the Gis1 jumonji domain. Finally, we note that several jumonji domain interactors are related to proteins that are found in mammalian PML nuclear bodies.

Yeast Mon1p/Aut12p functions in vacuolar fusion of autophagosomes and cvt-vesicles.

Here we identify Mon1p as being essential for the cvt-pathway and autophagy. Thus, mon1Delta cells are impaired in proaminopeptidase I maturation and homozygous diploid mon1Delta cells do not sporulate. Quantitative autophagy measurements suggest a complete autophagy block. The autophagosomal marker protein GFP-Aut7p accumulates in mon1Delta cells at punctate structures outside the vacuole. Furthermore, proaminopeptidase I accumulates in mon1Delta cells in a proteinase-protected form. Our data demonstrate that mon1Delta cells are defective in the fusion of cvt-vesicles and autophagosomes with the vacuole. Consistent with this, GFP-Mon1p localizes to the cytosol and to punctate structures within the cytosol.