Multi-ancestry GWAS reveals excitotoxicity associated with outcome after ischaemic stroke.

During the first hours after stroke onset, neurological deficits can be highly unstable: some patients rapidly improve, while others deteriorate. This early neurological instability has a major impact on long-term outcome. Here, we aimed to determine the genetic architecture of early neurological instability measured by the difference between the National Institutes of Health Stroke Scale (NIHSS) within 6 h of stroke onset and NIHSS at 24 h. A total of 5876 individuals from seven countries (Spain, Finland, Poland, USA, Costa Rica, Mexico and Korea) were studied using a multi-ancestry meta-analyses. We found that 8.7% of NIHSS at 24 h of variance was explained by common genetic variations, and also that early neurological instability has a different genetic architecture from that of stroke risk. Eight loci (1p21.1, 1q42.2, 2p25.1, 2q31.2, 2q33.3, 5q33.2, 7p21.2 and 13q31.1) were genome-wide significant and explained 1.8% of the variability suggesting that additional variants influence early change in neurological deficits. We used functional genomics and bioinformatic annotation to identify the genes driving the association from each locus. Expression quantitative trait loci mapping and summary data-based Mendelian randomization indicate that ADAM23 (log Bayes factor = 5.41) was driving the association for 2q33.3. Gene-based analyses suggested that GRIA1 (log Bayes factor = 5.19), which is predominantly expressed in the brain, is the gene driving the association for the 5q33.2 locus. These analyses also nominated GNPAT (log Bayes factor = 7.64) ABCB5 (log Bayes factor = 5.97) for the 1p21.1 and 7p21.1 loci. Human brain single-nuclei RNA-sequencing indicates that the gene expression of ADAM23 and GRIA1 is enriched in neurons. ADAM23, a presynaptic protein and GRIA1, a protein subunit of the AMPA receptor, are part of a synaptic protein complex that modulates neuronal excitability. These data provide the first genetic evidence in humans that excitotoxicity may contribute to early neurological instability after acute ischaemic stroke.

Genomic atlas of the proteome from brain, CSF and plasma prioritizes proteins implicated in neurological disorders.

Understanding the tissue-specific genetic controls of protein levels is essential to uncover mechanisms of post-transcriptional gene regulation. In this study, we generated a genomic atlas of protein levels in three tissues relevant to neurological disorders (brain, cerebrospinal fluid and plasma) by profiling thousands of proteins from participants with and without Alzheimer's disease. We identified 274, 127 and 32 protein quantitative trait loci (pQTLs) for cerebrospinal fluid, plasma and brain, respectively. cis-pQTLs were more likely to be tissue shared, but trans-pQTLs tended to be tissue specific. Between 48.0% and 76.6% of pQTLs did not co-localize with expression, splicing, DNA methylation or histone acetylation QTLs. Using Mendelian randomization, we nominated proteins implicated in neurological diseases, including Alzheimer's disease, Parkinson's disease and stroke. This first multi-tissue study will be instrumental to map signals from genome-wide association studies onto functional genes, to discover pathways and to identify drug targets for neurological diseases.

African Americans Have Differences in CSF Soluble TREM2 and Associated Genetic Variants.

OBJECTIVE: To evaluate for racial differences in triggering receptor expressed on myeloid cells 2 (TREM2), a key immune mediator in Alzheimer disease, the levels of CSF soluble TREM2 (sTREM2), and the frequency of associated genetic variants were compared in groups of individuals who self-reported their race as African American (AA) or non-Hispanic White (NHW). METHODS: Community-dwelling older research participants underwent measurement of CSF sTREM2 concentrations and genetic analyses. RESULTS: The primary cohort included 91 AAs and 868 NHWs. CSF sTREM2 levels were lower in the AA compared with the NHW group (1,336 ± 470 vs 1,856 ± 624 pg/mL, p < 0.0001). AAs were more likely to carry TREM2 coding variants (15% vs 3%, p < 0.0001), which were associated with lower CSF sTREM2. AAs were less likely to carry the rs1582763 minor allele (8% vs 37%, p < 0.0001), located near MS4A4A, which was associated with higher CSF sTREM2. These findings were replicated in an independent cohort of 23 AAs and 917 NHWs: CSF sTREM2 levels were lower in the AA group (p = 0.03), AAs were more likely to carry coding TREM2 variants (22% vs 4%, p = 0.002), and AAs were less likely to carry the rs1582763 minor allele (16% vs 37%, p = 0.003). CONCLUSIONS: On average, AAs had lower CSF sTREM2 levels compared with NHWs, potentially because AAs are more likely to carry genetic variants associated with lower CSF sTREM2 levels. Importantly, CSF sTREM2 reflects TREM2-mediated microglial activity, a critical step in the immune response to amyloid plaques. These findings suggest that race may be associated with risk for genetic variants that influence Alzheimer disease-related inflammation.

Quantitative endophenotypes as an alternative approach to understanding genetic risk in neurodegenerative diseases.

Endophenotypes, as measurable intermediate features of human diseases, reflect underlying molecular mechanisms. The use of quantitative endophenotypes in genetic studies has improved our understanding of pathophysiological changes associated with diseases. The main advantage of the quantitative endophenotypes approach to study human diseases over a classic case-control study design is the inferred biological context that can enable the development of effective disease-modifying treatments. Here, we summarize recent progress on biomarkers for neurodegenerative diseases, including cerebrospinal fluid and blood-based, neuroimaging, neuropathological, and clinical studies. This review focuses on how endophenotypic studies have successfully linked genetic modifiers to disease risk, disease onset, or progression rate and provided biological context to genes identified in genome-wide association studies. Finally, we review critical methodological considerations for implementing this approach and future directions.

Molecular pathways in patients with systemic lupus erythematosus revealed by gene-centred DNA sequencing.

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.

Chi3l1/YKL-40 is controlled by the astrocyte circadian clock and regulates neuroinflammation and Alzheimer's disease pathogenesis.

Regulation of glial activation and neuroinflammation are critical factors in the pathogenesis of Alzheimer's disease (AD). YKL-40, a primarily astrocytic protein encoded by the gene Chi3l1, is a widely studied cerebrospinal fluid biomarker that increases with aging and early in AD. However, the function of Chi3l1/YKL-40 in AD is unknown. In a cohort of patients with AD, we observed that a variant in the human CHI3L1 gene, which results in decreased CSF YKL-40 expression, was associated with slower AD progression. At baseline, Chi3l1 deletion in mice had no effect on astrocyte activation while modestly promoting microglial activation. In a mouse APP/PS1 model of AD, Chi3l1 deletion decreased amyloid plaque burden and increased periplaque expression of the microglial lysosomal marker CD68, suggesting that Chi3l1 may suppress glial phagocytic activation and promote amyloid accumulation. Accordingly, Chi3l1 knockdown increased phagocytosis of zymosan particles and of ß-amyloid peptide in both astrocytes and microglia in vitro. We further observed that expression of Chi3l1 is regulated by the circadian clock, as deletion of the core clock proteins BMAL1 or CLOCK/NPAS2 strongly suppresses basal Chi3l1 expression, whereas deletion of the negative clock regulators PER1/PER2 increased Chi3l1 expression. Basal Chi3l1 mRNA was nonrhythmic because of a long mRNA half-life in astrocytes. However, inflammatory induction of Chi3l1 was gated by the clock. Our findings reveal Chi3l1/YKL-40 as a modulator of glial phagocytic activation and AD pathogenesis in both mice and humans and suggest that the astrocyte circadian clock regulates inflammatory Chi3l1 induction.

Multi-ancestry genetic study in 5,876 patients identifies an association between excitotoxic genes and early outcomes after acute ischemic stroke.

During the first hours after stroke onset neurological deficits can be highly unstable: some patients rapidly improve, while others deteriorate. This early neurological instability has a major impact on long-term outcome. Here, we aimed to determine the genetic architecture of early neurological instability measured by the difference between NIH stroke scale (NIHSS) within six hours of stroke onset and NIHSS at 24h (ΔNIHSS). A total of 5,876 individuals from seven countries (Spain, Finland, Poland, United States, Costa Rica, Mexico and Korea) were studied using a multi-ancestry meta-analyses. We found that 8.7% of ΔNIHSS variance was explained by common genetic variations, and also that early neurological instability has a different genetic architecture than that of stroke risk. Seven loci (2p25.1, 2q31.2, 2q33.3, 4q34.3, 5q33.2, 6q26 and 7p21.1) were genome-wide significant and explained 2.1% of the variability suggesting that additional variants influence early change in neurological deficits. We used functional genomics and bioinformatic annotation to identify the genes driving the association from each loci. eQTL mapping and SMR indicate that ADAM23 (log Bayes Factor (LBF)=6.34) was driving the association for 2q33.3. Gene based analyses suggested that GRIA1 (LBF=5.26), which is predominantly expressed in brain, is the gene driving the association for the 5q33.2 locus. These analyses also nominated PARK2 (LBF=5.30) and ABCB5 (LBF=5.70) for the 6q26 and 7p21.1 loci. Human brain single nuclei RNA-seq indicates that the gene expression of ADAM23 and GRIA1 is enriched in neurons. ADAM23 , a pre-synaptic protein, and GRIA1 , a protein subunit of the AMPA receptor, are part of a synaptic protein complex that modulates neuronal excitability. These data provides the first evidence in humans that excitotoxicity may contribute to early neurological instability after acute ischemic stroke. RESEARCH INTO CONTEXT: Evidence before this study: No previous genome-wide association studies have investigated the genetic architecture of early outcomes after ischemic stroke.Added Value of this study: This is the first study that investigated genetic influences on early outcomes after ischemic stroke using a genome-wide approach, revealing seven genome-wide significant loci. A unique aspect of this genetic study is the inclusion of all of the major ethnicities by recruiting from participants throughout the world. Most genetic studies to date have been limited to populations of European ancestry.Implications of all available evidence: The findings provide the first evidence that genes implicating excitotoxicity contribute to human acute ischemic stroke, and demonstrates proof of principle that GWAS of acute ischemic stroke patients can reveal mechanisms involved in ischemic brain injury.

Functional genomic analyses uncover APOE-mediated regulation of brain and cerebrospinal fluid beta-amyloid levels in Parkinson disease.

Alpha-synuclein is the main protein component of Lewy bodies, the pathological hallmark of Parkinson's disease. However, genetic modifiers of cerebrospinal fluid (CSF) alpha-synuclein levels remain unknown. The use of CSF levels of amyloid beta1-42, total tau, and phosphorylated tau181 as quantitative traits in genetic studies have provided novel insights into Alzheimer's disease pathophysiology. A systematic study of the genomic architecture of CSF biomarkers in Parkinson's disease has not yet been conducted. Here, genome-wide association studies of CSF biomarker levels in a cohort of individuals with Parkinson's disease and controls (N = 1960) were performed. PD cases exhibited significantly lower CSF biomarker levels compared to controls. A SNP, proxy for APOE ε4, was associated with CSF amyloid beta1-42 levels (effect = - 0.5, p = 9.2 × 10-19). No genome-wide loci associated with CSF alpha-synuclein, total tau, or phosphorylated tau181 levels were identified in PD cohorts. Polygenic risk score constructed using the latest Parkinson's disease risk meta-analysis were associated with Parkinson's disease status (p = 0.035) and the genomic architecture of CSF amyloid beta1-42 (R2 = 2.29%; p = 2.5 × 10-11). Individuals with higher polygenic risk scores for PD risk presented with lower CSF amyloid beta1-42 levels (p = 7.3 × 10-04). Two-sample Mendelian Randomization revealed that CSF amyloid beta1-42 plays a role in Parkinson's disease (p = 1.4 × 10-05) and age at onset (p = 7.6 × 10-06), an effect mainly mediated by variants in the APOE locus. In a subset of PD samples, the APOE ε4 allele was associated with significantly lower levels of CSF amyloid beta1-42 (p = 3.8 × 10-06), higher mean cortical binding potentials (p = 5.8 × 10-08), and higher Braak amyloid beta score (p = 4.4 × 10-04). Together these results from high-throughput and hypothesis-free approaches converge on a genetic link between Parkinson's disease, CSF amyloid beta1-42, and APOE.

A new domestic cat genome assembly based on long sequence reads empowers feline genomic medicine and identifies a novel gene for dwarfism.

The domestic cat (Felis catus) numbers over 94 million in the USA alone, occupies households as a companion animal, and, like humans, suffers from cancer and common and rare diseases. However, genome-wide sequence variant information is limited for this species. To empower trait analyses, a new cat genome reference assembly was developed from PacBio long sequence reads that significantly improve sequence representation and assembly contiguity. The whole genome sequences of 54 domestic cats were aligned to the reference to identify single nucleotide variants (SNVs) and structural variants (SVs). Across all cats, 16 SNVs predicted to have deleterious impacts and in a singleton state were identified as high priority candidates for causative mutations. One candidate was a stop gain in the tumor suppressor FBXW7. The SNV is found in cats segregating for feline mediastinal lymphoma and is a candidate for inherited cancer susceptibility. SV analysis revealed a complex deletion coupled with a nearby potential duplication event that was shared privately across three unrelated cats with dwarfism and is found within a known dwarfism associated region on cat chromosome B1. This SV interrupted UDP-glucose 6-dehydrogenase (UGDH), a gene involved in the biosynthesis of glycosaminoglycans. Importantly, UGDH has not yet been associated with human dwarfism and should be screened in undiagnosed patients. The new high-quality cat genome reference and the compilation of sequence variation demonstrate the importance of these resources when searching for disease causative alleles in the domestic cat and for identification of feline biomedical models.

The TMEM106B FTLD-protective variant, rs1990621, is also associated with increased neuronal proportion.

Apart from amyloid ß deposition and tau neurofibrillary tangles, Alzheimer's disease (AD) is a neurodegenerative disorder characterized by neuronal loss and astrocytosis in the cerebral cortex. The goal of this study is to investigate genetic factors associated with the neuronal proportion in health and disease. To identify cell-autonomous genetic variants associated with neuronal proportion in cortical tissues, we inferred cellular population structure from bulk RNA-Seq derived from 1536 individuals. We identified the variant rs1990621 located in the TMEM106B gene region as significantly associated with neuronal proportion (p value = 6.40 × 10-07) and replicated this finding in an independent dataset (p value = 7.41 × 10-04) surpassing the genome-wide threshold in the meta-analysis (p value = 9.42 × 10-09). This variant is in high LD with the TMEM106B non-synonymous variant p.T185S (rs3173615; r2 = 0.98) which was previously identified as a protective variant for frontotemporal lobar degeneration (FTLD). We stratified the samples by disease status, and discovered that this variant modulates neuronal proportion not only in AD cases, but also several neurodegenerative diseases and in elderly cognitively healthy controls. Furthermore, we did not find a significant association in younger controls or schizophrenia patients, suggesting that this variant might increase neuronal survival or confer resilience to the neurodegenerative process. The single variant and gene-based analyses also identified an overall genetic association between neuronal proportion, AD and FTLD risk. These results suggest that common pathways are implicated in these neurodegenerative diseases, that implicate neuronal survival. In summary, we identified a protective variant in the TMEM106B gene that may have a neuronal protection effect against general aging, independent of disease status, which could help elucidate the relationship between aging and neuronal survival in the presence or absence of neurodegenerative disorders. Our findings suggest that TMEM106B could be a potential target for neuronal protection therapies to ameliorate cognitive and functional deficits.

Overlap in the Genetic Architecture of Stroke Risk, Early Neurological Changes, and Cardiovascular Risk Factors.

Background and Purpose- The genetic relationships between stroke risk, stroke severity, and early neurological changes are complex and not completely understood. Genetic studies have identified 32 all stroke risk loci. Polygenic risk scores can be used to compare the genetic architecture of related traits. In this study, we compare the genetic architecture of stroke risk, stroke severity, and early neurological changes with that of 2 stroke risk factors: type 2 diabetes mellitus (T2DM) and hypertension. Methods- We assessed the degree of overlap in the genetic architecture of stroke risk, T2DM, hypertension, and 2 acute stroke phenotypes based on the National Institutes of Health Stroke Scale (NIHSS), which ranges from 0 for no stroke symptoms to 21 to 42 for a severe stroke: baseline (within 6 hours after onset) and change in NIHSS (ΔNIHSS=NIHSS at baseline-NIHSS at 24 hours). This was done by (1) single-nucleotide polymorphism by single-nucleotide polymorphism comparison, (2) weighted polygenic risk scores with sentinel variants, and (3) whole-genome polygenic risk scores using multiple P thresholds. Results- We found evidence of genetic architecture overlap between stroke risk and T2DM ( P=2.53×10-169), hypertension ( P=3.93×10-04), and baseline NIHSS ( P=0.03). However, there was no evidence of overlap between ΔNIHSS and stroke risk, T2DM, or hypertension. Conclusions- The genetic architecture of stroke risk is correlated with that of T2DM, hypertension, and initial stroke severity (NIHSS within 6 hours of stroke onset). However, the genetic architecture of early neurological change after stroke (ΔNIHSS) is not correlated with that of ischemic stroke risk, T2DM, or hypertension. Thus, stroke risk and early neurological change after stroke have distinct genetic architectures.

TREM2 brain transcript-specific studies in AD and TREM2 mutation carriers.

BACKGROUND: Low frequency coding variants in TREM2 are associated with Alzheimer disease (AD) risk and cerebrospinal fluid (CSF) TREM2 protein levels are different between AD cases and controls. Similarly, TREM2 risk variant carriers also exhibit differential CSF TREM2 levels. TREM2 has three different alternative transcripts, but most of the functional studies only model the longest transcript. No studies have analyzed TREM2 expression levels or alternative splicing in brains from AD and cognitively normal individuals. We wanted to determine whether there was differential expression of TREM2 in sporadic-AD cases versus AD-TREM2 carriers vs sex- and aged-matched normal controls; and if this differential expression was due to a particular TREM2 transcript. METHODS: We analyzed RNA-Seq data from parietal lobe brain tissue from AD cases with TREM2 variants (n = 33), AD cases (n = 195) and healthy controls (n = 118), from three independent datasets using Kallisto and the R package tximport to determine the read count for each transcript and quantified transcript abundance as transcripts per million. RESULTS: The three TREM2 transcripts were expressed in brain cortex in the three datasets. We demonstrate for the first time that the transcript that lacks the transmembrane domain and encodes a soluble form of TREM2 (sTREM2) has an expression level around 60% of the canonical transcript, suggesting that around 25% of the sTREM2 protein levels could be explained by this transcript. We did not observe a difference in the overall TREM2 expression level between cases and controls. However, the isoform which lacks the 5' exon, but includes the transmembrane domain, was significantly lower in TREM2- p.R62H carriers than in AD cases (p = 0.007). CONCLUSION: Using bulk RNA-Seq data from three different cohorts, we were able to quantify the expression level of the three TREM2 transcripts, demonstrating: (1) all three transcripts of them are highly expressed in the human cortex, (2) that up to 25% of the sTREM2 may be due to the expression of a specific isoform and not TREM2 cleavage; and (3) that TREM2 risk variants do not affect expression levels, suggesting that the effect of the TREM2 variants on CSF levels occurs at post-transcriptional level.

A rare regulatory variant in the MEF2D gene affects gene regulation and splicing and is associated with a SLE sub-phenotype in Swedish cohorts.

Systemic lupus erythematosus (SLE) is an autoimmune disorder with heterogeneous clinical presentation and complex etiology involving the interplay between genetic, epigenetic, environmental and hormonal factors. Many common SNPs identified by genome wide-association studies (GWAS) explain only a small part of the disease heritability suggesting the contribution from rare genetic variants, undetectable in GWAS, and complex epistatic interactions. Using targeted re-sequencing of coding and conserved regulatory regions within and around 215 candidate genes selected on the basis of their known role in autoimmunity and genes associated with canine immune-mediated diseases, we identified a rare regulatory variant rs200395694:G > T located in intron 4 of the MEF2D gene encoding the myocyte-specific enhancer factor 2D transcription factor and associated with SLE in Swedish cohorts (504 SLE patients and 839 healthy controls, p = 0.014, CI = 1.1-10). Fisher's exact test revealed an association between the genetic variant and a triad of disease manifestations including Raynaud, anti-U1-ribonucleoprotein (anti-RNP), and anti-Smith (anti-Sm) antibodies (p = 0.00037) among the patients. The DNA-binding activity of the allele was further studied by EMSA, reporter assays, and minigenes. The region has properties of an active cell-specific enhancer, differentially affected by the alleles of rs200395694:G > T. In addition, the risk allele exerts an inhibitory effect on the splicing of the alternative tissue-specific isoform, and thus may modify the target gene set regulated by this isoform. These findings emphasize the potential of dissecting traits of complex diseases and correlating them with rare risk alleles with strong biological effects.

Evaluation of ADAMTS17 in Chinese Shar-Pei with primary open-angle glaucoma, primary lens luxation, or both.

OBJECTIVE To evaluate the coding regions of ADAMTS17 for potential mutations in Chinese Shar-Pei with a diagnosis of primary open-angle glaucoma (POAG), primary lens luxation (PLL), or both. ANIMALS 63 Shar-Pei and 96 dogs of other breeds. PROCEDURES ADAMTS17 exon resequencing was performed on buccal mucosal DNA from 10 Shar-Pei with a diagnosis of POAG, PLL, or both (affected dogs). A candidate causal variant sequence was identified, and additional dogs (53 Shar-Pei [11 affected and 42 unaffected] and 95 dogs of other breeds) were genotyped for the variant sequence by amplified fragment length polymorphism analysis. Total RNA was extracted from ocular tissues of 1 affected Shar-Pei and 1 ophthalmologically normal Golden Retriever; ADAMTS17 cDNA was reverse transcribed and sequenced, and ADAMTS17 expression was evaluated by quantitative reverse-transcription PCR assay. RESULTS All affected Shar-Pei were homozygous for a 6-bp deletion in exon 22 of ADAMTS17 predicted to affect the resultant protein. All unaffected Shar-Pei were heterozygous or homozygous for the wild-type allele. The variant sequence was significantly associated with affected status (diagnosis of POAG, PLL, or both). All dogs of other breeds were homozygous for the wild-type allele. The cDNA sequencing confirmed presence of the expected variant mRNA sequence in ocular tissue from the affected dog only. Gene expression analysis revealed a 4.24-fold decrease in the expression of ADAMTS17 in ocular tissue from the affected dog. CONCLUSIONS AND CLINICAL RELEVANCE Results supported that the phenotype (diagnosis of POAG, PLL, or both) is an autosomal recessive trait in Shar-Pei significantly associated with the identified mutation in ADAMTS17.

The Novel Evolution of the Sperm Whale Genome.

The sperm whale, made famous by Moby Dick, is one of the most fascinating of all ocean-dwelling species given their unique life history, novel physiological adaptations to hunting squid at extreme ocean depths, and their position as one of the earliest branching toothed whales (Odontoceti). We assembled the sperm whale (Physeter macrocephalus) genome and resequenced individuals from multiple ocean basins to identify new candidate genes for adaptation to an aquatic environment and infer demographic history. Genes crucial for skin integrity appeared to be particularly important in both the sperm whale and other cetaceans. We also find sperm whales experienced a steep population decline during the early Pleistocene epoch. These genomic data add new comparative insight into the evolution of whales.

The practical use of genome sequencing data in the management of a feline colony pedigree.

BACKGROUND: A higher prevalence of inherited disorders among companion animals are often rooted in their historical restricted artificial selection for a variety of observed phenotypes that eventually decreased genetic diversity. Cats have been afflicted with many inherited diseases due to domestication and intense breed selection. Advances in sequencing technology have generated a more comprehensive way to access genetic information from an individual, allowing identification of putative disease-causing variants and in practice a means to avoid their spread and thus better pedigree management. We examine variants in three domestic shorthair cats and then calculated overall genetic diversity to extrapolate the benefits of this data for breeding programs within a feline colony. RESULTS: We generated whole genome sequence (WGS) data for three related cats that belong to a large feline pedigree colony. Genome-wide coverage ranged from 27-32X, from which we identified 18 million variants in total. Previously known disease-causing variants were screened in our cats, but none carry any of these known disease alleles. Loss of function (LoF) variants, that are in genes associated with a detrimental phenotype in human or mice were chosen for further evaluation on the comparative impact inferred. A set of LoF variants were observed in four genes, each with predicted detrimental phenotypes as a result. However, none of our cats displayed the expected disease phenotypes. Inbreeding coefficients and runs of homozygosity were also evaluated as a measure of genetic diversity. We find low inbreeding coefficients and total runs of homozygosity, thus suggesting pedigree management of genetic relatedness is acceptable. CONCLUSIONS: The use of WGS of a small sampling among a large feline colony has enabled us to identify possible disease-causing variants, their genotype state and measure pedigree management of genetic diversity. We contend a limited but strategic sampling of feline colony individuals using WGS can inform veterinarians of future health anomalies and guide breeding practices to ensure healthy genetic diversity.

Comparison of cellular location and expression of Plakophilin-2 in epidermal cells from nonlesional atopic skin and healthy skin in German shepherd dogs.

BACKGROUND: Canine atopic dermatitis (CAD) is an inflammatory and pruritic allergic skin disease caused by interactions between genetic and environmental factors. Previously, a genome-wide significant risk locus on canine chromosome 27 for CAD was identified in German shepherd dogs (GSDs) and Plakophilin-2 (PKP2) was defined as the top candidate gene. PKP2 constitutes a crucial component of desmosomes and also is important in signalling, metabolic and transcriptional activities. OBJECTIVES: The main objective was to evaluate the role of PKP2 in CAD by investigating PKP2 expression and desmosome structure in nonlesional skin from CAD-affected (carrying the top GWAS SNP risk allele) and healthy GSDs. We also aimed at defining the cell types in the skin that express PKP2 and its intracellular location. ANIMALS/METHODS: Skin biopsies were collected from nine CAD-affected and five control GSDs. The biopsies were frozen for immunofluorescence and fixed for electron microscopy immunolabelling and morphology. RESULTS: We observed the novel finding of PKP2 expression in dendritic cells and T cells in dog skin. Moreover, we detected that PKP2 was more evenly expressed within keratinocytes compared to its desmosomal binding-partner plakoglobin. PKP2 protein was located in the nucleus and on keratin filaments attached to desmosomes. No difference in PKP2 abundance between CAD cases and controls was observed. CONCLUSION: Plakophilin-2 protein in dog skin is expressed in both epithelial and immune cells; based on its subcellular location its functional role is implicated in both nuclear and structural processes.

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.

Multiple Changes of Gene Expression and Function Reveal Genomic and Phenotypic Complexity in SLE-like Disease.

The complexity of clinical manifestations commonly observed in autoimmune disorders poses a major challenge to genetic studies of such diseases. Systemic lupus erythematosus (SLE) affects humans as well as other mammals, and is characterized by the presence of antinuclear antibodies (ANA) in patients' sera and multiple disparate clinical features. Here we present evidence that particular sub-phenotypes of canine SLE-related disease, based on homogenous (ANA(H)) and speckled ANA (ANA(S)) staining pattern, and also steroid-responsive meningitis-arteritis (SRMA) are associated with different but overlapping sets of genes. In addition to association to certain MHC alleles and haplotypes, we identified 11 genes (WFDC3, HOMER2, VRK1, PTPN3, WHAMM, BANK1, AP3B2, DAPP1, LAMTOR3, DDIT4L and PPP3CA) located on five chromosomes that contain multiple risk haplotypes correlated with gene expression and disease sub-phenotypes in an intricate manner. Intriguingly, the association of BANK1 with both human and canine SLE appears to lead to similar changes in gene expression levels in both species. Our results suggest that molecular definition may help unravel the mechanisms of different clinical features common between and specific to various autoimmune disease phenotypes in dogs and humans.