Testing for Allele-specific Expression from Human Brain Samples.

Many single nucleotide polymorphisms (SNPs) identified by genome-wide association studies exert their effects on disease risk as expression quantitative trait loci (eQTL) via allele-specific expression (ASE). While databases for probing eQTLs in tissues from normal individuals exist, one may wish to ascertain eQTLs or ASE in specific tissues or disease-states not characterized in these databases. Here, we present a protocol to assess ASE of two possible target genes (GPNMB and KLHL7) of a known genome-wide association study (GWAS) Parkinson's disease (PD) risk locus in postmortem human brain tissue from PD and neurologically normal individuals. This was done using a sequence of RNA isolation, cDNA library generation, enrichment for transcripts of interest using customizable cDNA capture probes, paired-end RNA sequencing, and subsequent analysis. This method provides increased sensitivity relative to traditional bulk RNAseq-based and a blueprint that can be extended to the study of other genes, tissues, and disease states. Key features • Analysis of GPNMB allele-specific expression (ASE) in brain lysates from cognitively normal controls (NC) and Parkinson's disease (PD) individuals. • Builds on the ASE protocol of Mayba et al. (2014) and extends application from cells to human tissue. • Increased sensitivity by enrichment for desired transcript via RNA CaptureSeq (Mercer et al., 2014). • Optimized for human brain lysates from cingulate gyrus, caudate nucleus, and cerebellum.

ATTCT and ATTCC repeat expansions in the ATXN10 gene affect disease penetrance of spinocerebellar ataxia type 10.

Spinocerebellar ataxia type 10 (SCA10) is an autosomal-dominant disorder caused by an expanded pentanucleotide repeat in the ATXN10 gene. This repeat expansion, when fully penetrant, has a size of 850-4,500 repeats. It has been shown that the repeat composition can be a modifier of disease, e.g., seizures. Here, we describe a Mexican kindred in which we identified both pure (ATTCT)n and mixed (ATTCT)n-(ATTCC)n expansions in the same family. We used amplification-free targeted sequencing and optical genome mapping to decipher the composition of these repeat expansions. We found a considerable degree of mosaicism of the repeat expansion. This mosaicism was confirmed in skin fibroblasts from individuals with ATXN10 expansions with RNAScope in situ hybridization. All affected family members with the mixed ATXN10 repeat expansion showed typical clinical signs of spinocerebellar ataxia and epilepsy. In contrast, individuals with the pure ATXN10 expansion present with Parkinson's disease or are unaffected, even in individuals more than 20 years older than the average age at onset for SCA10. Our findings suggest that the pure (ATTCT)n expansion is non-pathogenic, while repeat interruptions, e.g., (ATTCC)n, are necessary to cause SCA10. This mechanism has been recently described for several other repeat expansions including SCA31 (BEAN1), SCA37 (DAB1), and three loci for benign adult familial myoclonic epilepsy BAFME (SAMD12, TNRC6A, RAPGEF2). Therefore, long-read sequencing and optical genome mapping of the entire genomic structure of repeat expansions are critical for clinical practice and genetic counseling, as variations in the repeat can affect disease penetrance, symptoms, and disease trajectory.

GPNMB confers risk for Parkinson's disease through interaction with α-synuclein.

Many risk loci for Parkinson's disease (PD) have been identified by genome-wide association studies (GWASs), but target genes and mechanisms remain largely unknown. We linked the GWAS-derived chromosome 7 locus (sentinel single-nucleotide polymorphism rs199347) to GPNMB through colocalization analyses of expression quantitative trait locus and PD risk signals, confirmed by allele-specific expression studies in the human brain. In cells, glycoprotein nonmetastatic melanoma protein B (GPNMB) coimmunoprecipitated and colocalized with α-synuclein (aSyn). In induced pluripotent stem cell-derived neurons, loss of GPNMB resulted in loss of ability to internalize aSyn fibrils and develop aSyn pathology. In 731 PD and 59 control biosamples, GPNMB was elevated in PD plasma, associating with disease severity. Thus, GPNMB represents a PD risk gene with potential for biomarker development and therapeutic targeting.

No Evidence of the Vertical Transmission of Non-Virulent Infectious Salmon Anaemia Virus (ISAV-HPR0) in Farmed Atlantic Salmon.

The nonvirulent infectious salmon anaemia virus (ISAV-HPR0) is the putative progenitor for virulent-ISAV, and a potential risk factor for the development of infectious salmon anaemia (ISA). Understanding the transmission dynamics of ISAV-HPR0 is fundamental to proper management and mitigation strategies. Here, we demonstrate that ISAV-HPR0 causes prevalent and transient infections in all three production stages of Atlantic salmon in the Faroe Islands. Phylogenetic analysis of the haemagglutinin-esterase gene from 247 salmon showed a clear geographical structuring into two significantly distinct HPR0-subgroups, which were designated G2 and G4. Whereas G2 and G4 co-circulated in marine farms, Faroese broodfish were predominantly infected by G2, and smolt were predominantly infected by G4. This infection pattern was confirmed by our G2- and G4-specific RT-qPCR assays. Moreover, the HPR0 variants detected in Icelandic and Norwegian broodfish were never detected in the Faroe Islands, despite the extensive import of ova from both countries. Accordingly, the vertical transmission of HPR0 from broodfish to progeny is uncommon. Phylogenetic and statistical analysis suggest that HPR0 persists in the smolt farms as "house-strains", and that new HPR0 variants are occasionally introduced from the marine environment, probably by HPR0-contaminated sea-spray. Thus, high biosecurity-including water and air intake-is required to avoid the introduction of pathogens to the smolt farms.

Genomic Epidemiology of Salmonid Alphavirus in Norwegian Aquaculture Reveals Recent Subtype-2 Transmission Dynamics and Novel Subtype-3 Lineages.

Viral disease poses a major barrier to sustainable aquaculture, with outbreaks causing large economic losses and growing concerns for fish welfare. Genomic epidemiology can support disease control by providing rapid inferences on viral evolution and disease transmission. In this study, genomic epidemiology was used to investigate salmonid alphavirus (SAV), the causative agent of pancreas disease (PD) in Atlantic salmon. Our aim was to reconstruct SAV subtype-2 (SAV2) diversity and transmission dynamics in recent Norwegian aquaculture, including the origin of SAV2 in regions where this subtype is not tolerated under current legislation. Using nanopore sequencing, we captured ~90% of the SAV2 genome for n = 68 field isolates from 10 aquaculture production regions sampled between 2018 and 2020. Using time-calibrated phylogenetics, we infer that, following its introduction to Norway around 2010, SAV2 split into two clades (SAV2a and 2b) around 2013. While co-present at the same sites near the boundary of Møre og Romsdal and Trøndelag, SAV2a and 2b were generally detected in non-overlapping locations at more Southern and Northern latitudes, respectively. We provide evidence for recent SAV2 transmission over large distances, revealing a strong connection between Møre og Romsdal and SAV2 detected in 2019/20 in Rogaland. We also demonstrate separate introductions of SAV2a and 2b outside the SAV2 zone in Sognefjorden (Vestland), connected to samples from Møre og Romsdal and Trøndelag, respectively, and a likely 100 km Northward transmission of SAV2b within Trøndelag. Finally, we recovered genomes of SAV2a and SAV3 co-infecting single fish in Rogaland, involving novel SAV3 lineages that diverged from previously characterized strains >25 years ago. Overall, this study demonstrates useful applications of genomic epidemiology for tracking viral disease spread in aquaculture.

Genomic epidemiology reveals multiple introductions of SARS-CoV-2 from mainland Europe into Scotland.

Coronavirus disease 2019 (COVID-19) was first diagnosed in Scotland on 1 March 2020. During the first month of the outbreak, 2,641 cases of COVID-19 led to 1,832 hospital admissions, 207 intensive care admissions and 126 deaths. We aimed to identify the source and number of introductions of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) into Scotland using a combined phylogenetic and epidemiological approach. Sequencing of 1,314 SARS-CoV-2 viral genomes from available patient samples enabled us to estimate that SARS-CoV-2 was introduced to Scotland on at least 283 occasions during February and March 2020. Epidemiological analysis confirmed that early introductions of SARS-CoV-2 originated from mainland Europe (the majority from Italy and Spain). We identified subsequent early outbreaks in the community, within healthcare facilities and at an international conference. Community transmission occurred after 2 March, 3 weeks before control measures were introduced. Earlier travel restrictions or quarantine measures, both locally and internationally, would have reduced the number of COVID-19 cases in Scotland. The risk of multiple reintroduction events in future waves of infection remains high in the absence of population immunity.

Nanopore whole genome sequencing and partitioned phylogenetic analysis supports a new salmonid alphavirus genotype (SAV7).

Salmon pancreas disease virus, more commonly known as salmonid alphavirus (SAV), is a single-stranded positive sense RNA virus and the causative agent of pancreas disease and sleeping disease in salmonids. In this study, a unique strain of SAV previously isolated from ballan wrasse was subjected to whole genome sequencing using nanopore sequencing. In order to accurately examine the evolutionary history of this strain in comparison to other SAV strains, a partitioned phylogenetic analysis was performed to account for variation in the rate of evolution for both individual genes and codon positions. Partitioning the genome alignments almost doubled the observed branch lengths in the phylogenetic tree when compared to the more common approach of applying one model of substitution across the genome and significantly increased the statistical fit of the best-fitting models of nucleotide substitution. Based on the genomic data, a valid case can be made for the viral strain examined in this study to be considered a new SAV genotype. In addition, this study adds to a growing number of studies in which SAV has been found to infect non-salmonid fish, and as such we have suggested that the viral species name be amended to the more inclusive 'piscine alphavirus'.

The structural variation landscape in 492 Atlantic salmon genomes.

Structural variants (SVs) are a major source of genetic and phenotypic variation, but remain challenging to accurately type and are hence poorly characterized in most species. We present an approach for reliable SV discovery in non-model species using whole genome sequencing and report 15,483 high-confidence SVs in 492 Atlantic salmon (Salmo salar L.) sampled from a broad phylogeographic distribution. These SVs recover population genetic structure with high resolution, include an active DNA transposon, widely affect functional features, and overlap more duplicated genes retained from an ancestral salmonid autotetraploidization event than expected. Changes in SV allele frequency between wild and farmed fish indicate polygenic selection on behavioural traits during domestication, targeting brain-expressed synaptic networks linked to neurological disorders in humans. This study offers novel insights into the role of SVs in genome evolution and the genetic architecture of domestication traits, along with resources supporting reliable SV discovery in non-model species.

Genome Sequencing of SAV3 Reveals Repeated Seeding Events of Viral Strains in Norwegian Aquaculture.

Understanding the dynamics of pathogen transfer in aquaculture systems is essential to manage and mitigate disease outbreaks. The goal of this study was to understand recent transmission dynamics of salmonid alphavirus (SAV) in Norway. SAV causes significant economic impacts on farmed salmonids in European aquaculture. SAV is classified into six subtypes, with Norway having ongoing epidemics of SAV subtypes 2 and 3. These two viral subtypes are present in largely distinct geographic regions of Norway, with SAV2 present in Trondelag, SAV3 in Rogaland, Sogn og Fjordane, and Hordaland, and Møre og Romsdal having outbreaks of both subtypes. To determine likely transmission routes of Norwegian SAV an established Nanopore amplicon sequencing approach was used in the current study. After confirming the accuracy of this approach for distinguishing subtype level co-infections of SAV2 and SAV3, a hypothetical possibility in regions of neighboring epidemics, twenty-four SAV3 genomes were sequenced to characterize the current genetic diversity of SAV3 in Norwegian aquaculture. Sequencing was performed on naturally infected heart tissues originating from a range of geographic locations sampled between 2016 and 2019. Phylogenetic analyses revealed that the currently active SAV3 strains sampled comprise several distinct lineages sharing an ancestor that existed ∼15 years ago (95% HPD, 12.51-17.7 years) and likely in Hordaland. At least five of these lineages have not shared a common ancestor for 7.85 years (95% HPD, 5.39-10.96 years) or more. Furthermore, the ancestor of the strains that were sampled outside of Hordaland (Sogn of Fjordane and Rogaland) existed less than 8 years ago, indicating a lack of long-term viral reservoirs in these counties. This evident lack of geographically distinct subclades is compatible with a source-sink transmission dynamic explaining the long-term movements of SAV around Norway. Such anthropogenic transport of the virus indicates that at least for sink counties, biosecurity strategies might be effective in mitigating the ongoing SAV epidemic. Finally, genomic analyses of SAV sequences were performed, offering novel insights into the prevalence of SAV genomes containing defective deletions. Overall, this study improves our understanding of the recent transmission dynamics and biology of the SAV epidemic affecting Norwegian aquaculture.

Cell Type-Specific Intralocus Interactions Reveal Oligodendrocyte Mechanisms in MS.

Multiple sclerosis (MS) is an autoimmune disease characterized by attack on oligodendrocytes within the central nervous system (CNS). Despite widespread use of immunomodulatory therapies, patients may still face progressive disability because of failure of myelin regeneration and loss of neurons, suggesting additional cellular pathologies. Here, we describe a general approach for identifying specific cell types in which a disease allele exerts a pathogenic effect. Applying this approach to MS risk loci, we pinpoint likely pathogenic cell types for 70%. In addition to T cell loci, we unexpectedly identified myeloid- and CNS-specific risk loci, including two sites that dysregulate transcriptional pause release in oligodendrocytes. Functional studies demonstrated inhibition of transcriptional elongation is a dominant pathway blocking oligodendrocyte maturation. Furthermore, pause release factors are frequently dysregulated in MS brain tissue. These data implicate cell-intrinsic aberrations outside of the immune system and suggest new avenues for therapeutic development. VIDEO ABSTRACT.

Genetic diversity among geographically distant isolates of Neoparamoeba perurans.

The present study explored the use of 2 common genetic markers, the mitochondrial cytochrome oxidase I (COI) and the ribosomal internal transcribed spacer (ITS) to infer the relationship between geographically distant isolates of the protozoan gill parasite Neoparamoeba perurans, the agent responsible for amoebic gill disease in farmed Atlantic salmon worldwide. Present data confirmed that the ITS marker is suitable for Neoparamoeba species discrimination; however, it is not recommended as a population marker due to the presence of multiple copies of ITS within both N. perurans clonal and polycultures. On the other hand, in the partial COI gene analysed here, a low variability was observed, with 8 haplotypes recovered from N. perurans samples collected from Europe (Ireland, Norway, Scotland) and Tasmania (Australia). In Europe, the COI haplotypes which have more recently been detected in aquaculture are different to the haplotypes associated with the original gill disease emergence in Ireland in 1997 and Norway in 2006. The presence of unique COI haplotypes in different continents suggests the presence of multiple distinct reservoirs of the pathogen in both Europe and Tasmania. Isolates from additional geographical locations are required to fully understand the origins and routes for the spread of N. perurans worldwide.

Nanopore sequencing for rapid diagnostics of salmonid RNA viruses.

Analysis of pathogen genome variation is essential for informing disease management and control measures in farmed animals. For farmed fish, the standard approach is to use PCR and Sanger sequencing to study partial regions of pathogen genomes, with second and third-generation sequencing tools yet to be widely applied. Here we demonstrate rapid and accurate sequencing of two disease-causing viruses affecting global salmonid aquaculture, salmonid alphavirus (SAV) and infectious salmon anaemia virus (ISAV), using third-generation nanopore sequencing on the MinION platform (Oxford Nanopore Technologies). Our approach complements PCR from infected material with MinION sequencing to recover genomic information that matches near perfectly to Sanger-verified references. We use this method to present the first SAV subtype-6 genome, which branches as the sister to all other SAV lineages in a genome-wide phylogenetic reconstruction. MinION sequencing offers an effective strategy for fast, genome-wide analysis of fish viruses, with major potential applications for diagnostics and robust investigations into the origins and spread of disease outbreaks.

The Post-GWAS Era: From Association to Function.

During the past 12 years, genome-wide association studies (GWASs) have uncovered thousands of genetic variants that influence risk for complex human traits and diseases. Yet functional studies aimed at delineating the causal genetic variants and biological mechanisms underlying the observed statistical associations with disease risk have lagged. In this review, we highlight key advances in the field of functional genomics that may facilitate the derivation of biological meaning post-GWAS. We highlight the evidence suggesting that causal variants underlying disease risk often function through regulatory effects on the expression of target genes and that these expression effects might be modest and cell-type specific. We moreover discuss specific studies as proof-of-principle examples for current statistical, bioinformatic, and empirical bench-based approaches to downstream elucidation of GWAS-identified disease risk loci.

A Dementia-Associated Risk Variant near TMEM106B Alters Chromatin Architecture and Gene Expression.

Neurodegenerative diseases pose an extraordinary threat to the world's aging population, yet no disease-modifying therapies are available. Although genome-wide association studies (GWASs) have identified hundreds of risk loci for neurodegeneration, the mechanisms by which these loci influence disease risk are largely unknown. Here, we investigated the association between common genetic variants at the 7p21 locus and risk of the neurodegenerative disease frontotemporal lobar degeneration. We showed that variants associated with disease risk correlate with increased expression of the 7p21 gene TMEM106B and no other genes; co-localization analyses implicated a common causal variant underlying both association with disease and association with TMEM106B expression in lymphoblastoid cell lines and human brain. Furthermore, increases in the amount of TMEM106B resulted in increases in abnormal lysosomal phenotypes and cell toxicity in both immortalized cell lines and neurons. We then combined fine-mapping, bioinformatics, and bench-based approaches to functionally characterize all candidate causal variants at this locus. This approach identified a noncoding variant, rs1990620, that differentially recruits CTCF in lymphoblastoid cell lines and human brain to influence CTCF-mediated long-range chromatin-looping interactions between multiple cis-regulatory elements, including the TMEM106B promoter. Our findings thus provide an in-depth analysis of the 7p21 locus linked by GWASs to frontotemporal lobar degeneration, nominating a causal variant and causal mechanism for allele-specific expression and disease association at this locus. Finally, we show that genetic variants associated with risk of neurodegenerative diseases beyond frontotemporal lobar degeneration are enriched in CTCF-binding sites found in brain-relevant tissues, implicating CTCF-mediated gene regulation in risk of neurodegeneration more generally.

Evolution and Expression of Tissue Globins in Ray-Finned Fishes.

The globin gene family encodes oxygen-binding hemeproteins conserved across the major branches of multicellular life. The origins and evolutionary histories of complete globin repertoires have been established for many vertebrates, but there remain major knowledge gaps for ray-finned fish. Therefore, we used phylogenetic, comparative genomic and gene expression analyses to discover and characterize canonical "non-blood" globin family members (i.e., myoglobin, cytoglobin, neuroglobin, globin-X, and globin-Y) across multiple ray-finned fish lineages, revealing novel gene duplicates (paralogs) conserved from whole genome duplication (WGD) and small-scale duplication events. Our key findings were that: (1) globin-X paralogs in teleosts have been retained from the teleost-specific WGD, (2) functional paralogs of cytoglobin, neuroglobin, and globin-X, but not myoglobin, have been conserved from the salmonid-specific WGD, (3) triplicate lineage-specific myoglobin paralogs are conserved in arowanas (Osteoglossiformes), which arose by tandem duplication and diverged under positive selection, (4) globin-Y is retained in multiple early branching fish lineages that diverged before teleosts, and (5) marked variation in tissue-specific expression of globin gene repertoires exists across ray-finned fish evolution, including several previously uncharacterized sites of expression. In this respect, our data provide an interesting link between myoglobin expression and the evolution of air breathing in teleosts. Together, our findings demonstrate great-unrecognized diversity in the repertoire and expression of nonblood globins that has arisen during ray-finned fish evolution.

Common variant rs356182 near SNCA defines a Parkinson's disease endophenotype.

OBJECTIVE: Parkinson's disease (PD) presents clinically with several motor subtypes that exhibit variable treatment response and prognosis. Here, we investigated genetic variants for their potential association with PD motor phenotype and progression. METHODS: We screened 10 SNPs, previously associated with PD risk, for association with tremor-dominant (TD) versus postural-instability gait disorder (PIGD) motor subtypes. SNPs that correlated with the TD/PIGD ratio in a discovery cohort of 251 PD patients were then evaluated in a multi-site replication cohort of 559 PD patients. SNPs associated with motor phenotype in both cross-sectional cohorts were next evaluated for association with (1) rates of motor progression in a longitudinal subgroup of 230 PD patients and (2) brain alpha-synuclein (SNCA) expression in the GTEx (Genotype-Tissue Expression project) consortium database. RESULTS: Genotype at rs356182, near SNCA, correlated with the TD/PIGD ratio in both the discovery (Bonferroni-corrected P = 0.04) and replication cohorts (P = 0.02). The rs356182 GG genotype was associated with a more tremor-predominant phenotype and predicted a slower rate of motor progression (1-point difference in annual rate of UPDRS-III motor score change, P = 0.01). The rs356182 genotype was associated with SNCA expression in the cerebellum (P = 0.005). INTERPRETATION: Our study demonstrates that the GG genotype at rs356182 provides molecular definition for a clinically important endophenotype associated with (1) more tremor-predominant motor phenomenology, (2) slower rates of motor progression, and (3) decreased brain expression of SNCA. Such molecularly defined endophenotyping in PD may benefit both clinical trial design and tailoring of clinical care as we enter the era of precision medicine.

Lower plasma apolipoprotein A1 levels are found in Parkinson's disease and associate with apolipoprotein A1 genotype.

The discovery of novel plasma-based biomarkers could lead to new approaches in the treatment of Parkinson's disease (PD). Here, we explore the role of plasma apolipoprotein A1 (ApoA1) as a risk marker for PD and evaluate the influence of APOA1 promoter variation on plasma ApoA1 levels. Plasma ApoA1 and the single-nucleotide polymorphism, rs670, were assayed in a discovery cohort (cohort 1) of 301 PD patients, 80 normal controls (NCs), and 165 subjects with other neurodegenerative diseases, as well as a cohort (cohort 2) of 158 PD patients from a second clinical site. Additionally, rs670 was genotyped in a third cohort of 1,494 PD and 925 NC subjects from both clinical sites. Compared to both normal and disease controls, PD patients have lower plasma ApoA1 (P < 0.001 for both comparisons). Moreover, in PD patients, plasma ApoA1 levels are correlated with genotype at the APOA1 promoter polymorphism, rs670. Specifically, lower plasma ApoA1 levels were found in rs670 major allele (G) homozygotes in both cohort 1 (P = 0.009) and in a replication cohort (cohort 2; n = 158 PD patients; P = 0.024). Finally, evaluating rs670 genotype frequencies in 1,930 PD cases versus 997 NCs, the rs670 GG genotype shows a trend toward association (odds ratio: 1.1; P = 0.10) with PD. Our results are compatible with a model whereby circulating ApoA1 levels may be useful in risk-stratifying subjects for the development of PD, with higher ApoA1 levels suggesting relative protection. Future studies evaluating modulation of ApoA1 as a novel therapeutic strategy in PD are warranted. © 2014 International Parkinson and Movement Disorder Society.

TMEM106B is a genetic modifier of frontotemporal lobar degeneration with C9orf72 hexanucleotide repeat expansions.

Hexanucleotide repeat expansions in chromosome 9 open reading frame 72 (C9orf72) have recently been linked to frontotemporal lobar degeneration (FTLD) and amyotrophic lateral sclerosis, and may be the most common genetic cause of both neurodegenerative diseases. Genetic variants at TMEM106B influence risk for the most common neuropathological subtype of FTLD, characterized by inclusions of TAR DNA-binding protein of 43 kDa (FTLD-TDP). Previous reports have shown that TMEM106B is a genetic modifier of FTLD-TDP caused by progranulin (GRN) mutations, with the major (risk) allele of rs1990622 associating with earlier age at onset of disease. Here, we report that rs1990622 genotype affects age at death in a single-site discovery cohort of FTLD patients with C9orf72 expansions (n = 14), with the major allele correlated with later age at death (p = 0.024). We replicate this modifier effect in a 30-site international neuropathological cohort of FTLD-TDP patients with C9orf72 expansions (n = 75), again finding that the major allele associates with later age at death (p = 0.016), as well as later age at onset (p = 0.019). In contrast, TMEM106B genotype does not affect age at onset or death in 241 FTLD-TDP cases negative for GRN mutations or C9orf72 expansions. Thus, TMEM106B is a genetic modifier of FTLD with C9orf72 expansions. Intriguingly, the genotype that confers increased risk for developing FTLD-TDP (major, or T, allele of rs1990622) is associated with later age at onset and death in C9orf72 expansion carriers, providing an example of sign epistasis in human neurodegenerative disease.