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.
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.
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
, is a widely ...studied cerebrospinal fluid biomarker that increases with aging and early in AD. However, the function of
/YKL-40 in AD is unknown. In a cohort of patients with AD, we observed that a variant in the human
gene, which results in decreased CSF YKL-40 expression, was associated with slower AD progression. At baseline,
deletion in mice had no effect on astrocyte activation while modestly promoting microglial activation. In a mouse APP/PS1 model of AD,
deletion decreased amyloid plaque burden and increased periplaque expression of the microglial lysosomal marker CD68, suggesting that
may suppress glial phagocytic activation and promote amyloid accumulation. Accordingly,
knockdown increased phagocytosis of zymosan particles and of β-amyloid peptide in both astrocytes and microglia in vitro. We further observed that expression of
is regulated by the circadian clock, as deletion of the core clock proteins BMAL1 or CLOCK/NPAS2 strongly suppresses basal
expression, whereas deletion of the negative clock regulators PER1/PER2 increased
expression. Basal
mRNA was nonrhythmic because of a long mRNA half-life in astrocytes. However, inflammatory induction of
was gated by the clock. Our findings reveal
/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
induction.
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;
r
2
= 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.
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.
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.
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).
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.
Abstract A recessive, adult-onset neuronal ceroid-lipofuscinosis (NCL) occurs in Tibetan terriers. A genome-wide association study restricted this NCL locus to a 1.3 Mb region of canine chromosome 2 ...which contains canine ATP13A2 . NCL-affected dogs were homozygous for a single-base deletion in ATP13A2 , predicted to produce a frameshift and premature termination codon. Homozygous truncating mutations in human ATP13A2 have been shown by others to cause Kufor-Rakeb syndrome (KRS), a rare neurodegenerative disease. These findings suggest that KRS is also an NCL, although analysis of KRS brain tissue will be needed to confirm this prediction. Generalized brain atrophy, behavioral changes, and cognitive decline occur in both people and dogs with ATP13A2 mutations; however, other clinical features differ between the species. For example, Tibetan terriers with NCL develop cerebellar ataxia not reported in KRS patients and KRS patients exhibit parkinsonism and pyramidal dysfunction not observed in affected Tibetan terriers. To see if ATP13A2 mutations could be responsible for some cases of human adult-onset NCL (Kufs disease), we resequenced ATP13A2 from 28 Kufs disease patients. None of these patients had ATP13A2 sequence variants likely to be causal for their disease, suggesting that mutations in this gene are not common causes of Kufs disease.
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.
Humans and dogs are both affected by the allergic skin disease atopic dermatitis (AD), caused by an interaction between genetic and environmental factors. The German shepherd dog (GSD) is a high-risk ...breed for canine AD (CAD). In this study, we used a Swedish cohort of GSDs as a model for human AD. Serum IgA levels are known to be lower in GSDs compared to other breeds. We detected significantly lower IgA levels in the CAD cases compared to controls (p = 1.1 × 10(-5)) in our study population. We also detected a separation within the GSD cohort, where dogs could be grouped into two different subpopulations. Disease prevalence differed significantly between the subpopulations contributing to population stratification (λ = 1.3), which was successfully corrected for using a mixed model approach. A genome-wide association analysis of CAD was performed (n cases = 91, n controls = 88). IgA levels were included in the model, due to the high correlation between CAD and low IgA levels. In addition, we detected a correlation between IgA levels and the age at the time of sampling (corr = 0.42, p = 3.0 × 10(-9)), thus age was included in the model. A genome-wide significant association was detected on chromosome 27 (praw = 3.1 × 10(-7), pgenome = 0.03). The total associated region was defined as a ~1.5-Mb-long haplotype including eight genes. Through targeted re-sequencing and additional genotyping of a subset of identified SNPs, we defined 11 smaller haplotype blocks within the associated region. Two blocks showed the strongest association to CAD. The ~209-kb region, defined by the two blocks, harbors only the PKP2 gene, encoding Plakophilin 2 expressed in the desmosomes and important for skin structure. Our results may yield further insight into the genetics behind both canine and human AD.
ObjectivesSystemic 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.MethodsWe 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).ResultsWe 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.ConclusionsOur 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.
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 factorstype 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 strokebaseline (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), hypertension (P=3.93×10), 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.