Computing linkage disequilibrium aware genome embeddings using autoencoders.

MOTIVATION: The completion of the genome has paved the way for genome-wide association studies (GWAS), which explained certain proportions of heritability. GWAS are not optimally suited to detect non-linear effects in disease risk, possibly hidden in non-additive interactions (epistasis). Alternative methods for epistasis detection using, e.g. deep neural networks (DNNs) are currently under active development. However, DNNs are constrained by finite computational resources, which can be rapidly depleted due to increasing complexity with the sheer size of the genome. Besides, the curse of dimensionality complicates the task of capturing meaningful genetic patterns for DNNs; therefore necessitates dimensionality reduction. RESULTS: We propose a method to compress single nucleotide polymorphism (SNP) data, while leveraging the linkage disequilibrium (LD) structure and preserving potential epistasis. This method involves clustering correlated SNPs into haplotype blocks and training per-block autoencoders to learn a compressed representation of the block's genetic content. We provide an adjustable autoencoder design to accommodate diverse blocks and bypass extensive hyperparameter tuning. We applied this method to genotyping data from Project MinE, and achieved 99% average test reconstruction accuracy-i.e. minimal information loss-while compressing the input to nearly 10% of the original size. We demonstrate that haplotype-block based autoencoders outperform linear Principal Component Analysis (PCA) by approximately 3% chromosome-wide accuracy of reconstructed variants. To the extent of our knowledge, our approach is the first to simultaneously leverage haplotype structure and DNNs for dimensionality reduction of genetic data. AVAILABILITY AND IMPLEMENTATION: Data are available for academic use through Project MinE at https://www.projectmine.com/research/data-sharing/, contingent upon terms and requirements specified by the source studies. Code is available at https://github.com/gizem-tas/haploblock-autoencoders.

Longitudinal Effects of Asymptomatic C9orf72 Carriership on Brain Morphology.

OBJECTIVE: We investigated effects of C9orf72 repeat expansion and gene expression on longitudinal cerebral changes before symptom onset. METHODS: We enrolled 79 asymptomatic family members (AFMs) from 9 families with C9orf72 repeat expansion. Twenty-eight AFMs carried the mutation (C9+). Participants had up to 3 magnetic resonance imaging (MRI) scans, after which we compared motor cortex and motor tracts between C9+ and C9- AFMs using mixed effects models, incorporating kinship to correct for familial relations and lessen effects of other genetic factors. We also compared cortical, subcortical, cerebellar, and connectome structural measurements in a hypothesis-free analysis. We correlated regional C9orf72 expression in donor brains with the pattern of cortical thinning in C9+ AFMs using meta-regression. For comparison, we included 42 C9+ and 439 C9- patients with amyotrophic lateral sclerosis (ALS) in this analysis. RESULTS: C9+ AFM motor cortex had less gyrification and was thinner than in C9- AFMs, without differences in motor tracts. Whole brain analysis revealed thinner cortex and less gyrification in parietal, occipital, and temporal regions, smaller thalami and right hippocampus, and affected frontotemporal connections. Thinning of bilateral precentral, precuneus, and left superior parietal cortex was faster in C9+ than in C9- AFMs. Higher C9orf72 expression correlated with thinner cortex in both C9+ AFMs and C9+ ALS patients. INTERPRETATION: In asymptomatic C9orf72 repeat expansion carriers, brain MRI reveals widespread features suggestive of impaired neurodevelopment, along with faster decline of motor and parietal cortex than found in normal aging. C9orf72 expression might play a role in cortical development, and consequently explain the specific brain abnormalities of mutation carriers. ANN NEUROL 2023;93:668-680.

Common Genetic Variants Contribute to Risk of Transposition of the Great Arteries.

RATIONALE: Dextro-transposition of the great arteries (D-TGA) is a severe congenital heart defect which affects approximately 1 in 4,000 live births. While there are several reports of D-TGA patients with rare variants in individual genes, the majority of D-TGA cases remain genetically elusive. Familial recurrence patterns and the observation that most cases with D-TGA are sporadic suggest a polygenic inheritance for the disorder, yet this remains unexplored. OBJECTIVE: We sought to study the role of common single nucleotide polymorphisms (SNPs) in risk for D-TGA. METHODS AND RESULTS: We conducted a genome-wide association study in an international set of 1,237 patients with D-TGA and identified a genome-wide significant susceptibility locus on chromosome 3p14.3, which was subsequently replicated in an independent case-control set (rs56219800, meta-analysis P=8.6x10-10, OR=0.69 per C allele). SNP-based heritability analysis showed that 25% of variance in susceptibility to D-TGA may be explained by common variants. A genome-wide polygenic risk score derived from the discovery set was significantly associated to D-TGA in the replication set (P=4x10-5). The genome-wide significant locus (3p14.3) co-localizes with a putative regulatory element that interacts with the promoter of WNT5A, which encodes the Wnt Family Member 5A protein known for its role in cardiac development in mice. We show that this element drives reporter gene activity in the developing heart of mice and zebrafish and is bound by the developmental transcription factor TBX20. We further demonstrate that TBX20 attenuates Wnt5a expression levels in the developing mouse heart. CONCLUSIONS: This work provides support for a polygenic architecture in D-TGA and identifies a susceptibility locus on chromosome 3p14.3 near WNT5A. Genomic and functional data support a causal role of WNT5A at the locus.

Analysis of aneurysmal subarachnoid hemorrhage as a multistep process.

BACKGROUND AND PURPOSE: Aneurysmal subarachnoid hemorrhage (ASAH) is a complex disease with higher incidence in women compared to men and in Japan compared to other countries. It was hypothesized that ASAH is consistent with a multistep model of disease. The following assessments were made: (1) the number of steps needed for the disease to occur and (2) whether this number may be different in female versus male and in Japanese versus non-Japanese patients. METHODS: Incidence data were generated from a meta-analysis on ASAH incidence until 2017, which was supplemented with a literature search from 2017 to April 2023. Age- and sex-adjusted incidences per 10-year age groups were calculated and the logarithm of age-specific incidence against the logarithm of age was regressed with least-squares regression. RESULTS: In 2317 ASAH patients a linear relationship between logarithm of incidence and logarithm of age was found with a slope estimate of 3.13 (95% confidence interval 2.60-3.65), consistent with a four-step process. Similar estimates were found for female, male, Japanese and non-Japanese patients. CONCLUSIONS: Our results suggest that ASAH is a four-step process, also in subgroups with higher ASAH incidence. Elucidation of the exact nature of these steps can provide important clues for identification of disease mechanisms underlying ASAH.

Unexpected frequency of the pathogenic AR CAG repeat expansion in the general population.

CAG repeat expansions in exon 1 of the AR gene on the X chromosome cause spinal and bulbar muscular atrophy, a male-specific progressive neuromuscular disorder associated with a variety of extra-neurological symptoms. The disease has a reported male prevalence of approximately 1:30 000 or less, but the AR repeat expansion frequency is unknown. We established a pipeline, which combines the use of the ExpansionHunter tool and visual validation, to detect AR CAG expansion on whole-genome sequencing data, benchmarked it to fragment PCR sizing, and applied it to 74 277 unrelated individuals from four large cohorts. Our pipeline showed sensitivity of 100% [95% confidence interval (CI) 90.8-100%], specificity of 99% (95% CI 94.2-99.7%), and a positive predictive value of 97.4% (95% CI 84.4-99.6%). We found the mutation frequency to be 1:3182 (95% CI 1:2309-1:4386, n = 117 734) X chromosomes-10 times more frequent than the reported disease prevalence. Modelling using the novel mutation frequency led to estimate disease prevalence of 1:6887 males, more than four times more frequent than the reported disease prevalence. This discrepancy is possibly due to underdiagnosis of this neuromuscular condition, reduced penetrance, and/or pleomorphic clinical manifestations.

Genetic Risk Score for Intracranial Aneurysms: Prediction of Subarachnoid Hemorrhage and Role in Clinical Heterogeneity.

BACKGROUND: Recently, common genetic risk factors for intracranial aneurysm (IA) and aneurysmal subarachnoid hemorrhage (ASAH) were found to explain a large amount of disease heritability and therefore have potential to be used for genetic risk prediction. We constructed a genetic risk score to (1) predict ASAH incidence and IA presence (combined set of unruptured IA and ASAH) and (2) assess its association with patient characteristics. METHODS: A genetic risk score incorporating genetic association data for IA and 17 traits related to IA (so-called metaGRS) was created using 1161 IA cases and 407 392 controls from the UK Biobank population study. The metaGRS was validated in combination with risk factors blood pressure, sex, and smoking in 828 IA cases and 68 568 controls from the Nordic HUNT population study. Furthermore, we assessed association between the metaGRS and patient characteristics in a cohort of 5560 IA patients. RESULTS: Per SD increase of metaGRS, the hazard ratio for ASAH incidence was 1.34 (95% CI, 1.20-1.51) and the odds ratio for IA presence 1.09 (95% CI, 1.01-1.18). Upon including the metaGRS on top of clinical risk factors, the concordance index to predict ASAH hazard increased from 0.63 (95% CI, 0.59-0.67) to 0.65 (95% CI, 0.62-0.69), while prediction of IA presence did not improve. The metaGRS was statistically significantly associated with age at ASAH (ß=-4.82×10-3 per year [95% CI, -6.49×10-3 to -3.14×10-3]; P=1.82×10-8), and location of IA at the internal carotid artery (odds ratio=0.92 [95% CI, 0.86-0.98]; P=0.0041). CONCLUSIONS: The metaGRS was predictive of ASAH incidence, although with limited added value over clinical risk factors. The metaGRS was not predictive of IA presence. Therefore, we do not recommend using this metaGRS in daily clinical care. Genetic risk does partly explain the clinical heterogeneity of IA warranting prioritization of clinical heterogeneity in future genetic prediction studies of IA and ASAH.

The contribution of Neanderthal introgression and natural selection to neurodegenerative diseases.

Humans are thought to be more susceptible to neurodegeneration than equivalently-aged primates. It is not known whether this vulnerability is specific to anatomically-modern humans or shared with other hominids. The contribution of introgressed Neanderthal DNA to neurodegenerative disorders remains uncertain. It is also unclear how common variants associated with neurodegenerative disease risk are maintained by natural selection in the population despite their deleterious effects. In this study, we aimed to quantify the genome-wide contribution of Neanderthal introgression and positive selection to the heritability of complex neurodegenerative disorders to address these questions. We used stratified-linkage disequilibrium score regression to investigate the relationship between five SNP-based signatures of natural selection, reflecting different timepoints of evolution, and genome-wide associated variants of the three most prevalent neurodegenerative disorders: Alzheimer's disease, amyotrophic lateral sclerosis and Parkinson's disease. We found no evidence for enrichment of positively-selected SNPs in the heritability of Alzheimer's disease, amyotrophic lateral sclerosis and Parkinson's disease, suggesting that common deleterious disease variants are unlikely to be maintained by positive selection. There was no enrichment of Neanderthal introgression in the SNP-heritability of these disorders, suggesting that Neanderthal admixture is unlikely to have contributed to disease risk. These findings provide insight into the origins of neurodegenerative disorders within the evolution of Homo sapiens and addresses a long-standing debate, showing that Neanderthal admixture is unlikely to have contributed to common genetic risk of neurodegeneration in anatomically-modern humans.

MRI Clustering Reveals Three ALS Subtypes With Unique Neurodegeneration Patterns.

OBJECTIVE: The purpose of this study was to identify subtypes of amyotrophic lateral sclerosis (ALS) by comparing patterns of neurodegeneration using brain magnetic resonance imaging (MRI) and explore their phenotypes. METHODS: We performed T1-weighted and diffusion tensor imaging in 488 clinically well-characterized patients with ALS and 338 control subjects. Measurements of whole-brain cortical thickness and white matter connectome fractional anisotropy were adjusted for disease-unrelated variation. A probabilistic network-based clustering algorithm was used to divide patients into subgroups of similar neurodegeneration patterns. Clinical characteristics and cognitive profiles were assessed for each subgroup. In total, 512 follow-up scans were used to validate clustering results longitudinally. RESULTS: The clustering algorithm divided patients with ALS into 3 subgroups of 187, 163, and 138 patients. All subgroups displayed involvement of the precentral gyrus and are characterized, respectively, by (1) pure motor involvement (pure motor cluster [PM]), (2) orbitofrontal and temporal involvement (frontotemporal cluster [FT]), and (3) involvement of the posterior cingulate cortex, parietal white matter, temporal operculum, and cerebellum (cingulate-parietal-temporal cluster [CPT]). These subgroups had significantly distinct clinical profiles regarding male-to-female ratio, age at symptom onset, and frequency of bulbar symptom onset. FT and CPT revealed higher rates of cognitive impairment on the Edinburgh cognitive and behavioral ALS screen (ECAS). Longitudinally, clustering remained stable: at 90.4% of their follow-up visits, patients clustered in the same subgroup as their baseline visit. INTERPRETATION: ALS can manifest itself in 3 main patterns of cerebral neurodegeneration, each associated with distinct clinical characteristics and cognitive profiles. Besides the pure motor and frontotemporal dementia (FTD)-like variants of ALS, a new neuroimaging phenotype has emerged, characterized by posterior cingulate, parietal, temporal, and cerebellar involvement. ANN NEUROL 2022;92:1030-1045.

UNC13A in amyotrophic lateral sclerosis: from genetic association to therapeutic target.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease with limited treatment options and an incompletely understood pathophysiology. Although genomewide association studies (GWAS) have advanced our understanding of the disease, the precise manner in which risk polymorphisms contribute to disease pathogenesis remains unclear. Of relevance, GWAS have shown that a polymorphism (rs12608932) in the UNC13A gene is associated with risk for both ALS and frontotemporal dementia (FTD). Homozygosity for the C-allele at rs12608932 modifies the ALS phenotype, as these patients are more likely to have bulbar-onset disease, cognitive impairment and FTD at baseline as well as shorter survival. UNC13A is expressed in neuronal tissue and is involved in maintaining synaptic active zones, by enabling the priming and docking of synaptic vesicles. In the absence of functional TDP-43, risk variants in UNC13A lead to the inclusion of a cryptic exon in UNC13A messenger RNA, subsequently leading to nonsense mediated decay, with loss of functional protein. Depletion of UNC13A leads to impaired neurotransmission. Recent discoveries have identified UNC13A as a potential target for therapy development in ALS, with a confirmatory trial with lithium carbonate in UNC13A cases now underway and future approaches with antisense oligonucleotides currently under consideration. Considering UNC13A is a potent phenotypic modifier, it may also impact clinical trial outcomes. This present review describes the path from the initial discovery of UNC13A as a risk gene in ALS to the current therapeutic options being explored and how knowledge of its distinct phenotype needs to be taken into account in future trials.

Glutaminase Deficiency Caused by Short Tandem Repeat Expansion in GLS.

We report an inborn error of metabolism caused by an expansion of a GCA-repeat tract in the 5' untranslated region of the gene encoding glutaminase (GLS) that was identified through detailed clinical and biochemical phenotyping, combined with whole-genome sequencing. The expansion was observed in three unrelated patients who presented with an early-onset delay in overall development, progressive ataxia, and elevated levels of glutamine. In addition to ataxia, one patient also showed cerebellar atrophy. The expansion was associated with a relative deficiency of GLS messenger RNA transcribed from the expanded allele, which probably resulted from repeat-mediated chromatin changes upstream of the GLS repeat. Our discovery underscores the importance of careful examination of regions of the genome that are typically excluded from or poorly captured by exome sequencing.

Characterising ALS disease progression according to El Escorial and Gold Coast criteria.

BACKGROUND: The Gold Coast criteria (GCC) have been proposed as a means of selecting patients for amyotrophic lateral sclerosis (ALS) clinical trials. We aimed to characterise disease progression according to the GCC. METHODS: Data from population-based ALS registries from the Netherlands and Belgium were analysed. The GCC additionally define ALS as lower motor neuron (LMN) dysfunction in ≥2 body regions without upper motor neuron dysfunction. Therefore, the revised El Escorial criteria (rEEC) were supplemented with a 'Gold Coast ALS' category for patients with only LMN dysfunction in ≥2 body regions. We assessed survival time, ALS Functional Rating Scale (ALSFRS-R) progression rates and between-patient variability per diagnostic category. RESULTS: We included 5957 ALS patients, of whom 600 (10.1%) fulfilled the GCC but not the rEEC, and 95 (1.6%) fulfilled only the rEEC. ALSFRS-R progression rates were similar for the rEEC (0.84 points/month) and GCC (0.81 points/month) with similar variability (standard deviation of 0.59 vs. 0.60) and median survival time (17.8 vs.18.7 months). Survival time and average progression rates varied (p<0.001) between categories. Per category, however, there was considerable between-patient variability with progression rates ranging from: -2.10 to -0.14 (definite), -1.94 to -0.06 (probable), -2.10 to -0.02 (probable laboratory supported), -1.79 to -0.02 (possible) and -1.31 to 0.08 (Gold Coast). CONCLUSIONS: The GCC broaden the definition of ALS, allowing more patients to participate in trials, while minimally impacting population heterogeneity. Given the large variability per diagnostic category, selecting only specific categories for trials may not result in a more homogeneous study population.

Sensitivity of brain MRI and neurological examination for detection of upper motor neurone degeneration in amyotrophic lateral sclerosis.

OBJECTIVES: To investigate sensitivity of brain MRI and neurological examination for detection of upper motor neuron (UMN) degeneration in patients with amyotrophic lateral sclerosis (ALS). METHODS: We studied 192 patients with ALS and 314 controls longitudinally. All patients visited our centre twice and underwent full neurological examination and brain MRI. At each visit, we assessed UMN degeneration by measuring motor cortex thickness (CT) and pyramidal tract fibre density (FD) corresponding to five body regions (bulbar region and limbs). For each body region, we measured degree of clinical UMN and lower motor neuron (LMN) symptom burden using a validated scoring system. RESULTS: We found deterioration over time of CT of motor regions (p≤0.0081) and progression of UMN signs of bulbar region and left arm (p≤0.04). FD was discriminative between controls and patients with moderate/severe UMN signs (all regions, p≤0.034), but did not change longitudinally. Higher clinical UMN burden correlated with reduced CT, but not lower FD, for the bulbar region (p=2.2×10-10) and legs (p≤0.025). In the arms, we found that severe LMN signs may reduce the detectability of UMN signs (p≤0.043). With MRI, UMN degeneration was detectable before UMN signs became clinically evident (CT: p=1.1×10-10, FD: p=6.3×10-4). Motor CT, but not FD, deteriorated more than UMN signs during the study period. CONCLUSIONS: Motor CT is a more sensitive measure of UMN degeneration than UMN signs. Motor CT and pyramidal tract FD are discriminative between patients and controls. Brain MRI can monitor UMN degeneration before signs become clinically evident. These findings promote MRI as a potential biomarker for UMN progression in clinical trials in ALS.

Transethnic Genome-Wide Association Study Provides Insights in the Genetic Architecture and Heritability of Long QT Syndrome.

BACKGROUND: Long QT syndrome (LQTS) is a rare genetic disorder and a major preventable cause of sudden cardiac death in the young. A causal rare genetic variant with large effect size is identified in up to 80% of probands (genotype positive) and cascade family screening shows incomplete penetrance of genetic variants. Furthermore, a proportion of cases meeting diagnostic criteria for LQTS remain genetically elusive despite genetic testing of established genes (genotype negative). These observations raise the possibility that common genetic variants with small effect size contribute to the clinical picture of LQTS. This study aimed to characterize and quantify the contribution of common genetic variation to LQTS disease susceptibility. METHODS: We conducted genome-wide association studies followed by transethnic meta-analysis in 1656 unrelated patients with LQTS of European or Japanese ancestry and 9890 controls to identify susceptibility single nucleotide polymorphisms. We estimated the common variant heritability of LQTS and tested the genetic correlation between LQTS susceptibility and other cardiac traits. Furthermore, we tested the aggregate effect of the 68 single nucleotide polymorphisms previously associated with the QT-interval in the general population using a polygenic risk score. RESULTS: Genome-wide association analysis identified 3 loci associated with LQTS at genome-wide statistical significance (P<5×10-8) near NOS1AP, KCNQ1, and KLF12, and 1 missense variant in KCNE1(p.Asp85Asn) at the suggestive threshold (P<10-6). Heritability analyses showed that ≈15% of variance in overall LQTS susceptibility was attributable to common genetic variation (h2SNP 0.148; standard error 0.019). LQTS susceptibility showed a strong genome-wide genetic correlation with the QT-interval in the general population (rg=0.40; P=3.2×10-3). The polygenic risk score comprising common variants previously associated with the QT-interval in the general population was greater in LQTS cases compared with controls (P<10-13), and it is notable that, among patients with LQTS, this polygenic risk score was greater in patients who were genotype negative compared with those who were genotype positive (P<0.005). CONCLUSIONS: This work establishes an important role for common genetic variation in susceptibility to LQTS. We demonstrate overlap between genetic control of the QT-interval in the general population and genetic factors contributing to LQTS susceptibility. Using polygenic risk score analyses aggregating common genetic variants that modulate the QT-interval in the general population, we provide evidence for a polygenic architecture in genotype negative LQTS.

Role of Rare Genetic Variants Found in Families With Intracranial Aneurysms in the General Dutch and UK Population.

[Figure: see text].

Advances in the genetic classification of amyotrophic lateral sclerosis.

PURPOSE OF REVIEW: Amyotrophic lateral sclerosis (ALS) is an archetypal complex disease wherein disease risk and severity are, for the majority of patients, the product of interaction between multiple genetic and environmental factors. We are in a period of unprecedented discovery with new large-scale genome-wide association study (GWAS) and accelerating discovery of risk genes. However, much of the observed heritability of ALS is undiscovered and we are not yet approaching elucidation of the total genetic architecture, which will be necessary for comprehensive disease subclassification. RECENT FINDINGS: We summarize recent developments and discuss the future. New machine learning models will help to address nonlinear genetic interactions. Statistical power for genetic discovery may be boosted by reducing the search-space using cell-specific epigenetic profiles and expanding our scope to include genetically correlated phenotypes. Structural variation, somatic heterogeneity and consideration of environmental modifiers represent significant challenges which will require integration of multiple technologies and a multidisciplinary approach, including clinicians, geneticists and pathologists. SUMMARY: The move away from fully penetrant Mendelian risk genes necessitates new experimental designs and new standards for validation. The challenges are significant, but the potential reward for successful disease subclassification is large-scale and effective personalized medicine.

The Distinct Traits of the UNC13A Polymorphism in Amyotrophic Lateral Sclerosis.

OBJECTIVE: The rs12608932 single nucleotide polymorphism in UNC13A is associated with amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) susceptibility, and may underlie differences in treatment response. We aimed to characterize the clinical, cognitive, behavioral, and neuroimaging phenotype of UNC13A in patients with ALS. METHODS: We included 2,216 patients with ALS without a C9orf72 mutation to identify clinical characteristics associated with the UNC13A polymorphism. A subcohort of 428 patients with ALS was used to study cognitive and behavioral profiles, and 375 patients to study neuroimaging characteristics. Associations were analyzed under an additive genetic model. RESULTS: Genotyping rs12608932 resulted in 854 A/A, 988 A/C, and 374 C/C genotypes. The C allele was associated with a higher age at symptom onset (median years A/A 63.5, A/C 65.6, and C/C 65.5; p < 0.001), more frequent bulbar onset (A/A 29.6%, A/C 31.8%, and C/C 43.1%; p < 0.001), higher incidences of ALS-FTD (A/A 4.3%, A/C 5.2%, and C/C 9.5%; p = 0.003), lower forced vital capacity at diagnosis (median percentage A/A 92.0, A/C 90.0, and C/C 86.5; p < 0.001), and a shorter survival (median in months A/A 33.3, A.C 30.7, and C/C 26.6; p < 0.001). UNC13A was associated with lower scores on ALS-specific cognition tests (means A/A 79.5, A/C 78.1, and C/C 76.6; p = 0.037), and more frequent behavioral disturbances (A/A 16.7%, A/C 24.4%, and C/C 27.7%; p = 0.045). Thinner left inferior temporal and right fusiform cortex were associated with the UNC13A single nucleotide polymorphism (SNP; p = 0.045 and p = 0.036). INTERPRETATION: Phenotypical distinctions associated with UNC13A make it an important factor to take into account in clinical trial design, studies on cognition and behavior, and prognostic counseling. ANN NEUROL 2020;88:796-806.

RNA-Sequencing Highlights Inflammation and Impaired Integrity of the Vascular Wall in Brain Arteriovenous Malformations.

Background and Purpose- Interventional treatment of unruptured brain arteriovenous malformations (BAVMs) has become increasingly controversial. Because medical therapy is still lacking, we aimed to obtain insight into the disease mechanisms implicated in BAVMs and to identify potential targets for medical treatment to prevent rupture of a BAVM. Methods- We used next-generation RNA sequencing to identify differential expression on a transcriptome-wide level comparing tissue samples of 12 BAVMs to 16 intracranial control arteries. We identified differentially expressed genes by negative binominal generalized log-linear regression (false discovery rate corrected P<0.05). We selected 10 genes for validation using droplet digital polymerase chain reaction. We performed functional pathway analysis accounting for potential gene-length bias, to establish enhancement of biological pathways involved in BAVMs. We further assessed which Gene Ontology terms were enriched. Results- We found 736 upregulated genes in BAVMs including genes implicated in the cytoskeletal machinery and cell-migration and genes encoding for inflammatory cytokines and secretory products of neutrophils and macrophages. Furthermore, we found 498 genes downregulated including genes implicated in extracellular matrix composition, the binary angiopoietin-TIE system, and TGF (transforming growth factor)-ß signaling. We confirmed the differential expression of top 10 ranked genes. Functional pathway analysis showed enrichment of the protein digestion and absorption pathway (false discovery rate-adjusted P=1.70×10-2). We identified 47 enriched Gene Ontology terms (false discovery rate-adjusted P<0.05) implicated in cytoskeleton network, cell-migration, endoplasmic reticulum, transmembrane transport, and extracellular matrix composition. Conclusions- Our genome-wide RNA-sequencing study points to involvement of inflammatory mediators, loss of cerebrovascular quiescence, and impaired integrity of the vascular wall in the pathophysiology of BAVMs. Our study may lend support to potential receptivity of BAVMs to medical therapeutics, including those promoting vessel maturation, and anti-inflammatory and immune-modifying drugs.

Detection of long repeat expansions from PCR-free whole-genome sequence data.

Identifying large expansions of short tandem repeats (STRs), such as those that cause amyotrophic lateral sclerosis (ALS) and fragile X syndrome, is challenging for short-read whole-genome sequencing (WGS) data. A solution to this problem is an important step toward integrating WGS into precision medicine. We developed a software tool called ExpansionHunter that, using PCR-free WGS short-read data, can genotype repeats at the locus of interest, even if the expanded repeat is larger than the read length. We applied our algorithm to WGS data from 3001 ALS patients who have been tested for the presence of the C9orf72 repeat expansion with repeat-primed PCR (RP-PCR). Compared against this truth data, ExpansionHunter correctly classified all (212/212, 95% CI [0.98, 1.00]) of the expanded samples as either expansions (208) or potential expansions (4). Additionally, 99.9% (2786/2789, 95% CI [0.997, 1.00]) of the wild-type samples were correctly classified as wild type by this method with the remaining three samples identified as possible expansions. We further applied our algorithm to a set of 152 samples in which every sample had one of eight different pathogenic repeat expansions, including those associated with fragile X syndrome, Friedreich's ataxia, and Huntington's disease, and correctly flagged all but one of the known repeat expansions. Thus, ExpansionHunter can be used to accurately detect known pathogenic repeat expansions and provides researchers with a tool that can be used to identify new pathogenic repeat expansions.

ExpansionHunter: a sequence-graph-based tool to analyze variation in short tandem repeat regions.

SUMMARY: We describe a novel computational method for genotyping repeats using sequence graphs. This method addresses the long-standing need to accurately genotype medically important loci containing repeats adjacent to other variants or imperfect DNA repeats such as polyalanine repeats. Here we introduce a new version of our repeat genotyping software, ExpansionHunter, that uses this method to perform targeted genotyping of a broad class of such loci. AVAILABILITY AND IMPLEMENTATION: ExpansionHunter is implemented in C++ and is available under the Apache License Version 2.0. The source code, documentation, and Linux/macOS binaries are available at https://github.com/Illumina/ExpansionHunter/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Pharmacogenetic interactions in amyotrophic lateral sclerosis: a step closer to a cure?

Genetic mutations related to amyotrophic lateral sclerosis (ALS) act through distinct pathophysiological pathways, which may lead to varying treatment responses. Here we assess the genetic interaction between C9orf72, UNC13A, and MOBP with creatine and valproic acid treatment in two clinical trials. Genotypic data was available for 309 of the 338 participants (91.4%). The UNC13A genotype affected mortality (p = 0.012), whereas C9orf72 repeat-expansion carriers exhibited a faster rate of decline in overall (p = 0.051) and bulbar functioning (p = 0.005). A dose-response pharmacogenetic interaction was identified between creatine and the A allele of the MOBP genotype (p = 0.027), suggesting a qualitative interaction in a recessive model (HR 3.96, p = 0.015). Not taking genetic information into account may mask evidence of response to treatment or be an unrecognized source of bias. Incorporating genetic data could help investigators to identify critical treatment clues in patients with ALS.