Explore the vast range of titles available.

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. Yang et al. generated a genomic atlas of protein levels in brain, cerebrospinal fluid and plasma and used human genetics approaches to identify proteins implicated in neurological diseases as well as druggable targets.

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.

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.

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 beta 1–42 , total tau, and phosphorylated tau 181 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 beta 1–42 levels (effect = − 0.5, p = 9.2 × 10 −19 ). No genome-wide loci associated with CSF alpha-synuclein, total tau, or phosphorylated tau 181 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 beta 1–42 (R 2  = 2.29%; p = 2.5 × 10 −11 ). Individuals with higher polygenic risk scores for PD risk presented with lower CSF amyloid beta 1–42 levels (p = 7.3 × 10 −04 ). Two-sample Mendelian Randomization revealed that CSF amyloid beta 1–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 beta 1–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 beta 1–42 , and APOE .

Background: Fanconi syndrome is a disorder of renal proximal tubule transport characterized by metabolic acidosis, amino aciduria, glucosuria, and phosphaturia. There are acquired and hereditary forms of this disorder. A late-onset form of Fanconi syndrome in Basenjis was first described in 1976 and is now recognized as an inherited disease in these dogs. In part because of the late onset of disease signs, the disorder has not been eradicated from the breed by selective mating. A study was therefore undertaken to identify the molecular genetic basis of the disease so that dogs could be screened prior to breeding in order to avoid generating affected offspring. Methods: Linkage analysis within a large family of Basenjis that included both affected and unaffected individuals was performed to localize the causative variant within the genome. Significant linkage was identified between chromosome 3 (CFA3) makers and the disease phenotype. Fine mapping restricted the region to a 2.7 Mb section of CFA3. A whole genome sequence of a Basenji affected with Fanconi syndrome was generated, and the sequence data were examined for the presence of potentially deleterious homozygous variants within the mapped region. Results: A homozygous 317 bp deletion was identified in the last exon of FAN1 of the proband. 78 Basenjis of known disease status were genotyped for the deletion variant. Among these dogs, there was almost complete concordance between genotype and phenotype. The only exception was one dog that was homozygous for the deletion variant but did not exhibit signs of Fanconi syndrome. Conclusions: These data indicate that the disorder is very likely the result of FAN1 deficiency. The mechanism by which this deficiency causes the disease signs remains to be elucidated. FAN1 has endonuclease and exonuclease activity that catalyzes incisions in regions of double-stranded DNA containing interstrand crosslinks. FAN1 inactivation may cause Fanconi syndrome in Basenjis by sensitization of kidney proximal tubule cells to toxin-mediated DNA crosslinking, resulting in the accumulation of genomic and mitochondrial DNA damage in the kidney. Differential exposure to environmental toxins that promote DNA crosslink formation may explain the wide age-at-onset variability for the disorder in Basenjis.

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.

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.

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.

Background L-2-hydroxyglutaric aciduria is a metabolic repair deficiency characterized by elevated levels of L-2-hydroxyglutaric acid in urine, blood and cerebrospinal fluid. Neurological signs associated with the disease in humans and dogs include seizures, ataxia and dementia. Case presentation Here we describe an 8 month old Yorkshire terrier that presented with episodes of hyperactivity and aggressive behavior. Between episodes, the dog’s behavior and neurologic examinations were normal. A T2 weighted MRI of the brain showed diffuse grey matter hyperintensity and a urine metabolite screen showed elevated 2-hydroxyglutaric acid. We sequenced all 10 exons and intron-exon borders of L2HGDH from the affected dog and identified a homozygous A to G transition in the initiator methionine codon. The first inframe methionine is at p.M183 which is past the mitochondrial targeting domain of the protein. Initiation of translation at p.M183 would encode an N-terminal truncated protein unlikely to be functional. Conclusions We have identified a mutation in the initiation codon of L2HGDH that is likely to result in a non-functional gene. The Yorkshire terrier could serve as an animal model to understand the pathogenesis of L-2-hydroxyglutaric aciduria and to evaluate potential therapies.