An intronic VNTR affects splicing of ABCA7 and increases risk of Alzheimer's disease

Multisample motif discovery and visualization for tandem repeats

Tandem repeats (TRs) occupy a significant portion of the human genome and are a source of polymorphisms due to variations in sizes and motif compositions. Some of these variations have been associated with various neuropathological disorders, highlighting the clinical importance of assessing the motif structure of TRs. Moreover, assessing the TR motif variation can offer valuable insights into evolutionary dynamics and population structure. Previously, characterizations of TRs were limited by short-read sequencing technology, which lacks the ability to accurately capture the full TR sequences. As long-read sequencing becomes more accessible and can capture the full complexity of TRs, there is now also a need for tools to characterize and analyze TRs using long-read data across multiple samples. In this study, we present MotifScope, a novel algorithm for the characterization and visualization of TRs based on a de novo k-mer approach for motif discovery. Comparative analysis against established tools reveals that MotifScope can identify a greater number of motifs and more accurately represent the underlying repeat sequences. Moreover, MotifScope has been specifically designed to enable motif composition comparisons across assemblies of different individuals, as well as across long-read sequencing reads within an individual, through combined motif discovery and sequence alignment. We showcase potential applications of MotifScope in diverse fields, including population genetics, clinical settings, and forensic analyses.

Analysis of targeted and whole genome sequencing of PacBio HiFi reads for a comprehensive genotyping of gene-proximal and phenotype-associated Variable Number Tandem Repeats

Variable Number Tandem repeats (VNTRs) refer to repeating motifs of size greater than five bp. VNTRs are an important source of genetic variation, and have been associated with multiple Mendelian and complex phenotypes. However, the highly repetitive structures require reads to span the region for accurate genotyping. Pacific Biosciences HiFi sequencing spans large regions and is highly accurate but relatively expensive. Therefore, targeted sequencing approaches coupled with long-read sequencing have been proposed to improve efficiency and throughput. In this paper, we systematically explored the trade-off between targeted and whole genome HiFi sequencing for genotyping VNTRs. We curated a set of 10 , 787 gene-proximal (G-)VNTRs, and 48 phenotype-associated (P-)VNTRs of interest. Illumina reads only spanned 46% of the G-VNTRs and 71% of P-VNTRs, motivating the use of HiFi sequencing. We performed targeted sequencing with hybridization by designing custom probes for 9,999 VNTRs and sequenced 8 samples using HiFi and Illumina sequencing, followed by adVNTR genotyping. We compared these results against HiFi whole genome sequencing (WGS) data from 28 samples in the Human Pangenome Reference Consortium (HPRC). With the targeted approach only 4,091 (41%) G-VNTRs and only 4 (8%) of P-VNTRs were spanned with at least 15 reads. A smaller subset of 3,579 (36%) G-VNTRs had higher median coverage of at least 63 spanning reads. The spanning behavior was consistent across all 8 samples. Among 5,638 VNTRs with low-coverage ( < 15), 67% were located within GC-rich regions ( > 60%). In contrast, the 40X WGS HiFi dataset spanned 98% of all VNTRs and 49 (98%) of P-VNTRs with at least 15 spanning reads, albeit with lower coverage. Spanning reads were sufficient for accurate genotyping in both cases. Our findings demonstrate that targeted sequencing provides consistently high coverage for a small subset of low-GC VNTRs, but WGS is more effective for broad and sufficient sampling of a large number of VNTRs.

Whole Genome Variable Number Tandem Repeat Analysis in Alzheimer Disease

Background and Objectives

Investigation into different allelic variants may yield new associative genes to predict late-onset Alzheimer disease (LOAD). Variable number tandem repeats (VNTRs) are important polymorphic components of the genome; however, they have been previously overlooked because of their complex genotyping. New software can now determine differing lengths of VNTRs; however, this has not been tested in a large case-control population.

Methods

We used VNTRseek to genotype over 200,000 tandem repeats in 9,501 cases and controls from the Alzheimer's Disease Sequencing Project. We first identified limiting factors of this analysis and then examined the association of VNTRs with AD diagnosis in a subset of non-Hispanic White participants.

Results

We found that VNTRs were highly associated with areas of the genome with a high number of previously identified variants. From our case-control analysis, we identified 9 VNTRs with a repeat allele length associated with LOAD including VNTRs on DSC3, NR2E3, CCNY, PKP4, GRAP, and MAP6.

Discussion

We were able to show the feasibility of this new type of analysis in large-scale whole-genome sequencing data and identify promising VNTRs that are associated with LOAD.

Approximating edit distances between complex tandem repeats efficiently

MOTIVATION

Extended tandem repeats (TRs) have been associated with 60 or more diseases over the past 30 years. Although most TRs have single repeat units (or motifs), complex TRs with different units have recently been correlated with some brain disorders. Of note, a population-scale analysis shows that complex TRs at one locus can be divergent, and different units are often expanded between individuals. To understand the evolution of high TR diversity, it is informative to visualize a phylogenetic tree. To do this, we need to measure the edit distance between pairs of complex TRs by considering duplication and contraction of units created by replication slippage. However, traditional rigorous algorithms for this purpose are computationally expensive.

RESULTS

We here propose an efficient heuristic algorithm to estimate the edit distance with duplication and contraction of units (EDDC, for short). We select a set of frequent units that occur in given complex TRs, encode each unit as a single symbol, compress a TR into an optimal series of unit symbols that partially matches the original TR with the minimum Levenshtein distance, and estimate the EDDC between a pair of complex TRs from their compressed forms. Using substantial synthetic benchmark datasets, we demonstrate that the estimated EDDC is highly correlated with the accurate EDDC, with a Pearson correlation coefficient of >0.983, while the heuristic algorithm achieves orders of magnitude performance speedup.

AVAILABILITY AND IMPLEMENTATION

The software program hEDDC that implements the proposed algorithm is available at https://github.com/Ricky-pon/hEDDC (DOI: 10.5281/zenodo.14732958).

Nanopore Long-Read Sequencing Unveils Genomic Disruptions in Alzheimer's Disease

Studies in laboratory models and postmortem human brain tissue from patients with Alzheimer's disease have revealed disruption of basic cellular processes such as DNA repair and epigenetic control as drivers of neurodegeneration. While genomic alterations in regions of the genome that are rich in repetitive sequences, often termed "dark regions," are difficult to resolve using traditional sequencing approaches, long-read technologies offer promising new avenues to explore previously inaccessible regions of the genome. In the current study, we leverage nanopore-based long-read whole-genome sequencing of DNA extracted from postmortem human frontal cortex at early and late stages of Alzheimer's disease, as well as age-matched controls, to analyze retrotransposon insertion events, non-allelic homologous recombination (NAHR), structural variants and DNA methylation within retrotransposon loci and other repetitive/dark regions of the human genome. Interestingly, we find that retrotransposon insertion events and repetitive element-associated NAHR are particularly enriched within centromeric and pericentromeric regions of DNA in the aged human brain, and that ribosomal DNA (rDNA) is subject to a high degree of NAHR compared to other regions of the genome. We detect a trending increase in potential somatic retrotransposition events of the small interfering nuclear element (SINE) AluY in late-stage Alzheimer's disease, and differential changes in methylation within repetitive elements and retrotransposons according to disease stage. Taken together, our analysis provides the first long-read DNA sequencing-based analysis of retrotransposon sequences, NAHR, structural variants, and DNA methylation in the aged brain, and points toward transposable elements, centromeric/pericentromeric regions and rDNA as hotspots for genomic variation.

Toward understanding the role of genomic repeat elements in neurodegenerative diseases

Neurodegenerative diseases cause great medical and economic burdens for both patients and society; however, the complex molecular mechanisms thereof are not yet well understood. With the development of high-coverage sequencing technology, researchers have started to notice that genomic repeat regions, previously neglected in search of disease culprits, are active contributors to multiple neurodegenerative diseases. In this review, we describe the association between repeat element variants and multiple degenerative diseases through genome-wide association studies and targeted sequencing. We discuss the identification of disease-relevant repeat element variants, further powered by the advancement of long-read sequencing technologies and their related tools, and summarize recent findings in the molecular mechanisms of repeat element variants in brain degeneration, such as those causing transcriptional silencing or RNA-mediated gain of toxic function. Furthermore, we describe how in silico predictions using innovative computational models, such as deep learning language models, could enhance and accelerate our understanding of the functional impact of repeat element variants. Finally, we discuss future directions to advance current findings for a better understanding of neurodegenerative diseases and the clinical applications of genomic repeat elements.

Multi-omic quantitative trait loci link tandem repeat size variation to gene regulation in human brain

Tandem repeat (TR) size variation is implicated in ~50 neurological disorders, yet its impact on gene regulation in the human brain remains largely unknown. In the present study, we quantified the impact of TR size variation on brain gene regulation across distinct molecular phenotypes, based on 4,412 multi-omics samples from 1,597 donors, including 1,586 newly sequenced ones. We identified ~2.2 million TR molecular quantitative trait loci (TR-xQTLs), linking ~139,000 unique TRs to nearby molecular phenotypes, including many known disease-risk TRs, such as the G2C4 expansion in C9orf72 associated with amyotrophic lateral sclerosis. Fine-mapping revealed ~18,700 TRs as potential causal variants. Our in vitro experiments further confirmed the causal and independent regulatory effects of three TRs. Additional colocalization analysis indicated the potential causal role of TR variation in brain-related phenotypes, highlighted by a 3'-UTR TR in NUDT14 linked to cortical surface area and a TG repeat in PLEKHA1, associated with Alzheimer's disease.

Characterizing tandem repeat complexities across long-read sequencing platforms with TREAT and otter

Tandem repeats (TRs) play important roles in genomic variation and disease risk in humans. Long-read sequencing allows for the accurate characterization of TRs; however, the underlying bioinformatics perspectives remain challenging. We present otter and TREAT: otter is a fast targeted local assembler, cross-compatible across different sequencing platforms. It is integrated in TREAT, an end-to-end workflow for TR characterization, visualization, and analysis across multiple genomes. In a comparison with existing tools based on long-read sequencing data from both Oxford Nanopore Technology (ONT, Simplex and Duplex) and Pacific Bioscience (PacBio, Sequel II and Revio), otter and TREAT achieve state-of-the-art genotyping and motif characterization accuracy. Applied to clinically relevant TRs, TREAT/otter significantly identify individuals with pathogenic TR expansions. When applied to a case-control setting, we replicate previously reported associations of TRs with Alzheimer's disease, including those near or within APOC1 (P = 2.63 × 10-9), SPI1 (P = 6.5 × 10-3), and ABCA7 (P = 0.04) genes. Finally, we use TREAT/otter to systematically evaluate potential biases when genotyping TRs using diverse ONT and PacBio long-read sequencing data sets. We show that, in rare cases (0.06%), long-read sequencing from coverage drops in TRs, including the disease-associated TRs in ABCA7 and RFC1 genes. Such coverage drops can lead to TR misgenotyping, hampering the accurate characterization of TR alleles. Taken together, our tools can accurately genotype TRs across different sequencing technologies and with minimal requirements, allowing end-to-end analysis and comparisons of TRs in human genomes, with broad applications in research and clinical fields.

Accurate Detection of Tandem Repeats from Error-Prone Sequences with EquiRep

A tandem repeat is a sequence of nucleotides that occurs as multiple contiguous and near-identical copies positioned next to each other. These repeats play critical roles in genetic diversity, gene regulation, and are strongly linked to various neurological and developmental disorders. While several methods exist for detecting tandem repeats, they often exhibit low accuracy when the repeat unit length increases or the number of copies is low. Furthermore, methods capable of handling highly mutated sequences remain scarce, highlighting a significant opportunity for improvement. We introduce EquiRep, a tool for accurate detection of tandem repeats from erroneous sequences. EquiRep estimates the likelihood of positions originating from the same position in the unit by self-alignment followed by a novel approach that refines the estimation. The built equivalent classes and the consecutive position information will be then used to build a weighted graph, and the cycle in this graph with maximum bottleneck weight while covering most nucleotide positions will be identified to reconstruct the repeat unit. We test EquiRep on simulated and real HOR and RCA datasets where it consistently outperforms or is comparable to state-of-the-art methods. EquiRep is robust to sequencing errors, and is able to make better predictions for long units and low frequencies which underscores its broad usability for studying tandem repeats.

African origin haplotype protective for Alzheimer's disease in APOEε4 carriers: exploring potential mechanisms

APOEε4 is the strongest genetic risk factor for Alzheimer's disease (AD) with approximately 50% of AD patients carrying at least one APOEε4 allele. Our group identified a protective interaction between APOEε4 with the African-specific A allele of rs10423769, which reduces the AD risk effect of APOEε4 homozygotes by approximately 75%. The protective variant lies 2Mb from APOE in a region of segmental duplications (SD) of chromosome 19 containing a cluster of pregnancy specific beta-1 glycoprotein genes (PSGs) and a long non-coding RNA. Using both short and long read sequencing, we demonstrate that rs10423769_A allele lies within a unique single haplotype inside this region of segmental duplication. We identified the protective haplotype in all African ancestry populations studied, including both West and East Africans, suggesting the variant has an old origin. Long-read sequencing identified both structural and DNA methylation differences between the protective rs10423769_A allele and non-protective haplotypes. An expanded variable number tandem repeat (VNTR) containing multiple MEF2 family transcription factor binding motifs was found associated with the protective haplotype (p-value = 2.9e-10). These findings provide novel insights into the mechanisms of this African-origin protective variant for AD in APOEε4 carriers and supports the importance of including all ancestries in AD research.

Very low levels of ABCA7 in the cerebrum and Alzheimer's disease onset between the ages of 60 and 80 independently of APOE

This cross-sectional study addressed the ABCA7-Alzheimer's disease (AD) association. ABCA7 protein levels were quantified in 3 cerebral regions of brain donors with Braak neurofibrillary tangle (NFT) stages 0-V. Ordinal regression models were implemented to estimate the effect of ABCA7 on stopping in an earlier Braak NFT stage versus progressing to the later stages in 2 prespecified age segments. In the final model, high ABCA7 levels in the parietal cortex increased the odds of remaining cognitively healthy (ie, in stages 0/I) versus experiencing AD onset (ie, progressing to stages II-V) in the 61-80 age segment (OR = 2.87, adj 95% CI = 1.41-7.86, adj P = .007, n = 109), after controlling for APOE and other covariates. No ABCA7-AD association was found in the 81-98 age segment (n = 113). Parietal ABCA7 levels in 61-80-year-old with stages II-V were very low, even significantly lower than in 81-98-year-old with stages II-V. ABCA7 levels in the prefrontal cortex and hippocampus predicted AD onset in the 61-80 age segment after adjustment for APOE. ABCA7 levels were also the lowest in 61-80-year-old with frequent neuritic plaques. Thus, very low ABCA7 levels in the cerebrum are associated with AD onset in the 7th-8th decade of life.

Single-cell atlas of ABCA7 loss-of-function reveals impaired neuronal respiration via choline-dependent lipid imbalances

Loss-of-function (LoF) variants in the lipid transporter ABCA7 significantly increase the risk of Alzheimer's disease (odds ratio ∼2), yet the pathogenic mechanisms and the neural cell types affected by these variants remain largely unknown. Here, we performed single-nuclear RNA sequencing of 36 human post-mortem samples from the prefrontal cortex of 12 ABCA7 LoF carriers and 24 matched non-carrier control individuals. ABCA7 LoF was associated with gene expression changes in all major cell types. Excitatory neurons, which expressed the highest levels of ABCA7, showed transcriptional changes related to lipid metabolism, mitochondrial function, cell cycle-related pathways, and synaptic signaling. ABCA7 LoF-associated transcriptional changes in neurons were similarly perturbed in carriers of the common AD missense variant ABCA7 p.Ala1527Gly (n = 240 controls, 135 carriers), indicating that findings from our study may extend to large portions of the at-risk population. Consistent with ABCA7's function as a lipid exporter, lipidomic analysis of isogenic iPSC-derived neurons (iNs) revealed profound intracellular triglyceride accumulation in ABCA7 LoF, which was accompanied by a relative decrease in phosphatidylcholine abundance. Metabolomic and biochemical analyses of iNs further indicated that ABCA7 LoF was associated with disrupted mitochondrial bioenergetics that suggested impaired lipid breakdown by uncoupled respiration. Treatment of ABCA7 LoF iNs with CDP-choline (a rate-limiting precursor of phosphatidylcholine synthesis) reduced triglyceride accumulation and restored mitochondrial function, indicating that ABCA7 LoF-induced phosphatidylcholine dyshomeostasis may directly disrupt mitochondrial metabolism of lipids. Treatment with CDP-choline also rescued intracellular amyloid β -42 levels in ABCA7 LoF iNs, further suggesting a link between ABCA7 LoF metabolic disruptions in neurons and AD pathology. This study provides a detailed transcriptomic atlas of ABCA7 LoF in the human brain and mechanistically links ABCA7 LoF-induced lipid perturbations to neuronal energy dyshomeostasis. In line with a growing body of evidence, our study highlights the central role of lipid metabolism in the etiology of Alzheimer's disease.

A comprehensive tandem repeat catalog of the human genome

With the increasing availability of long-read sequencing data, high-quality human genome assemblies, and software for fully characterizing tandem repeats, genome-wide genotyping of tandem repeat loci on a population scale becomes more feasible. Such efforts not only expand our knowledge of the tandem repeat landscape in the human genome but also enhance our ability to differentiate pathogenic tandem repeat mutations from benign polymorphisms. To this end, we analyzed 272 genomes assembled using datasets from three public initiatives that employed different long-read sequencing technologies. Here, we report a catalog of over 18 million tandem repeat loci, many of which were previously unannotated. Some of these loci are highly polymorphic, and many of them reside within coding sequences.

Structural characterization of a polymorphic repeat at the CACNA1C schizophrenia locus

Genetic variation within intron 3 of the CACNA1C calcium channel gene is associated with schizophrenia and bipolar disorder, but analysis of the causal variants and their effect is complicated by a nearby variable-number tandem repeat (VNTR). Here, we used 155 long-read genome assemblies from 78 diverse individuals to delineate the structure and population variability of the CACNA1C intron 3 VNTR. We categorized VNTR sequences into 7 Types of structural alleles using sequence differences among repeat units. Only 12 repeat units at the 5' end of the VNTR were shared across most Types, but several Types were related through a series of large and small duplications. The most diverged Types were rare and present only in individuals with African ancestry, but the multiallelic structural polymorphism Variable Region 2 was present across populations at different frequencies, consistent with expansion of the VNTR preceding the emergence of early hominins. VR2 was in complete linkage disequilibrium with fine-mapped schizophrenia variants (SNPs) from genome-wide association studies (GWAS). This risk haplotype was associated with decreased CACNA1C gene expression in brain tissues profiled by the GTEx project. Our work suggests that sequence variation within a human-specific VNTR affects gene expression, and provides a detailed characterization of new alleles at a flagship neuropsychiatric locus.

Brain cholesterol and Alzheimer's disease: challenges and opportunities in probe and drug development

Cholesterol homeostasis is impaired in Alzheimer's disease; however, attempts to modulate brain cholesterol biology have not translated into tangible clinical benefits for patients to date. Several recent milestone developments have substantially improved our understanding of how excess neuronal cholesterol contributes to the pathophysiology of Alzheimer's disease. Indeed, neuronal cholesterol was linked to the formation of amyloid-β and neurofibrillary tangles through molecular pathways that were recently delineated in mechanistic studies. Furthermore, remarkable advances in translational molecular imaging have now made it possible to probe cholesterol metabolism in the living human brain with PET, which is an important prerequisite for future clinical trials that target the brain cholesterol machinery in Alzheimer's disease patients-with the ultimate aim being to develop disease-modifying treatments. This work summarizes current concepts of how the biosynthesis, transport and clearance of brain cholesterol are affected in Alzheimer's disease. Further, current strategies to reverse these alterations by pharmacotherapy are critically discussed in the wake of emerging translational research tools that support the assessment of brain cholesterol biology not only in animal models but also in patients with Alzheimer's disease.

Emerging Role of ABC Transporters in Glia Cells in Health and Diseases of the Central Nervous System

ATP-binding cassette (ABC) transporters play a crucial role for the efflux of a wide range of substrates across different cellular membranes. In the central nervous system (CNS), ABC transporters have recently gathered significant attention due to their pivotal involvement in brain physiology and neurodegenerative disorders, such as Alzheimer's disease (AD). Glial cells are fundamental for normal CNS function and engage with several ABC transporters in different ways. Here, we specifically highlight ABC transporters involved in the maintenance of brain homeostasis and their implications in its metabolic regulation. We also show new aspects related to ABC transporter function found in less recognized diseases, such as Huntington's disease (HD) and experimental autoimmune encephalomyelitis (EAE), as a model for multiple sclerosis (MS). Understanding both their impact on the physiological regulation of the CNS and their roles in brain diseases holds promise for uncovering new therapeutic options. Further investigations and preclinical studies are warranted to elucidate the complex interplay between glial ABC transporters and physiological brain functions, potentially leading to effective therapeutic interventions also for rare CNS disorders.

Exploring the Relationship between CLPTM1L-MS2 Variants and Susceptibility to Bladder Cancer

CLPTM1L (Cleft Lip and Palate Transmembrane Protein 1-Like) has previously been implicated in tumorigenesis and drug resistance in cancer. However, the genetic link between CLPTM1L and bladder cancer remains uncertain. In this study, we investigated the genetic association of variable number of tandem repeats (VNTR; minisatellites, MS) regions within CLPTM1L with bladder cancer. We identified four CLPTM1L-MS regions (MS1~MS4) located in intron regions. To evaluate the VNTR polymorphic alleles, we analyzed 441 cancer-free controls and 181 bladder cancer patients. Our analysis revealed a higher frequency of specific repeat sizes within the MS2 region in bladder cancer cases compared to controls. Notably, 25 and 27 repeats were exclusively present in the bladder cancer group. Moreover, rare alleles within the medium-length repeat range (25-29 repeats) were associated with an elevated bladder cancer risk (odds ratio [OR] = 5.78, 95% confidence interval [CI]: 1.49-22.47, p = 0.004). We confirmed that all MS regions followed Mendelian inheritance, and demonstrated that MS2 alleles increased CLPTM1L promoter activity in the UM-UC3 bladder cancer cells through a luciferase assay. Our findings propose the utility of CLPTM1L-MS regions as DNA typing markers, particularly highlighting the potential of middle-length rare alleles within CLPTM1L-MS2 as predictive markers for bladder cancer risk.

Mutated Toll-like receptor 9 increases Alzheimer's disease risk by compromising innate immunity protection

The development of Alzheimer's disease (AD) involves central and peripheral immune deregulation. Gene identification and studies of AD genetic variants of peripheral immune components may aid understanding of peripheral-central immune crosstalk and facilitate new opportunities for therapeutic intervention. In this study, we have identified in a Flanders-Belgian family a novel variant p.E317D in the Toll-like receptor 9 gene (TLR9), co-segregating with EOAD in an autosomal dominant manner. In human, TLR9 is an essential innate and adaptive immune component predominantly expressed in peripheral immune cells. The p.E317D variant caused 50% reduction in TLR9 activation in the NF-κB luciferase assay suggesting that p.E317D is a loss-of-function mutation. Cytokine profiling of human PBMCs upon TLR9 activation revealed a predominantly anti-inflammatory response in contrast to the inflammatory responses from TLR7/8 activation. The cytokines released upon TLR9 activation suppressed inflammation and promoted phagocytosis of Aβ42 oligomers in human iPSC-derived microglia. Transcriptome analysis identified upregulation of AXL, RUBICON and associated signaling pathways, which may underline the effects of TLR9 signaling-induced cytokines in regulating the inflammatory status and phagocytic property of microglia. Our data suggest a protective role of TLR9 signaling in AD pathogenesis, and we propose that TLR9 loss-of-function may disrupt a peripheral-central immune crosstalk that promotes dampening of inflammation and clearance of toxic protein species, leading to the build-up of neuroinflammation and pathogenic protein aggregates in AD development.

CSF biomarker analysis of ABCA7 mutation carriers suggests altered APP processing and reduced inflammatory response

BACKGROUND

The Alzheimer's disease (AD) risk gene ABCA7 has suggested functions in lipid metabolism and the immune system. Rare premature termination codon (PTC) mutations and an expansion of a variable number of tandem repeats (VNTR) polymorphism in the gene, both likely cause a lower ABCA7 expression and hereby increased risk for AD. However, the exact mechanism of action remains unclear. By studying CSF biomarkers reflecting different types of AD-related pathological processes, we aim to get a better insight in those processes and establish a biomarker profile of mutation carriers.

METHODS

The study population consisted of 229 AD patients for whom CSF was available and ABCA7 sequencing and VNTR genotyping had been performed. This included 28 PTC mutation and 16 pathogenic expansion carriers. CSF levels of Aβ1-42, Aβ1-40, P-tau181, T-tau, sAPPα, sAPPβ, YKL-40, and hFABP were determined using ELISA and Meso Scale Discovery assays. We compared differences in levels of these biomarkers and the Aβ ratio between AD patients with or without an ABCA7 PTC mutation or expansion using linear regression on INT-transformed data with APOE-status, age and sex as covariates.

RESULTS

Carriers of ABCA7 expansion mutations had significantly lower Aβ1-42 levels (P = 0.022) compared with non-carrier patients. The effect of the presence of ABCA7 mutations on CSF levels was especially pronounced in APOE ε4-negative carriers. In addition, VNTR expansion carriers had reduced Aβ1-40 (P = 0.023), sAPPα (P = 0.047), sAPPβ (P = 0.016), and YKL-40 (P = 0.0036) levels.

CONCLUSIONS

Our results are suggestive for an effect on APP processing by repeat expansions given the changes in the amyloid-related CSF biomarkers that were found in carriers. The decrease in YKL-40 levels in expansion carriers moreover suggests that these patients potentially have a reduced inflammatory response to AD damage. Moreover, our findings suggest the existence of a mechanism, independent of lowered expression, affecting neuropathology in expansion carriers.

ABCA7-Associated Clinical Features and Molecular Mechanisms in Alzheimer's Disease

Alzheimer's disease (AD) is the most common type of neurodegenerative disease and its pathogenesis is still unclear. Genetic factors are thought to account for a large proportion of the overall AD phenotypes. ATP-binding cassette transporter A7 (ABCA7) is one of the most important risk gene for AD. Multiple forms of ABCA7 variants significantly increase the risk of AD, such as single-nucleotide polymorphisms, premature termination codon variants, missense variants, variable number tandem repeat, mutations, and alternative splicing. AD patients with ABCA7 variants usually exhibit typical clinical and pathological features of traditional AD with a wide age of onset range. ABCA7 variants can alter ABCA7 protein expression levels and protein structure to affect protein functions such as abnormal lipid metabolism, amyloid precursor protein (APP) processing, and immune cell function. Specifically, ABCA7 deficiency can cause neuronal apoptosis by inducing endoplasmic reticulum stress through the PERK/eIF2α pathway. Second, ABCA7 deficiency can increase Aβ production by upregulating the SREBP2/BACE1 pathway and promoting APP endocytosis. In addition, the ability of microglia to phagocytose and degrade Aβ is destroyed by ABCA7 deficiency, leading to reduced clearance of Aβ. Finally, disturbance of lipid metabolism may also be an important method by which ABCA7 variants influence the incidence rate of AD. In the future, more attention should be given to different ABCA7 variants and ABCA7 targeted therapies for AD.

A landscape of complex tandem repeats within individual human genomes

Markedly expanded tandem repeats (TRs) have been correlated with ~60 diseases. TR diversity has been considered a clue toward understanding missing heritability. However, haplotype-resolved long TRs remain mostly hidden or blacked out because their complex structures (TRs composed of various units and minisatellites containing >10-bp units) make them difficult to determine accurately with existing methods. Here, using a high-precision algorithm to determine complex TR structures from long, accurate reads of PacBio HiFi, an investigation of 270 Japanese control samples yields several genome-wide findings. Approximately 322,000 TRs are difficult to impute from the surrounding single-nucleotide variants. Greater genetic divergence of TR loci is significantly correlated with more events of younger replication slippage. Complex TRs are more abundant than single-unit TRs, and a tendency for complex TRs to consist of <10-bp units and single-unit TRs to be minisatellites is statistically significant at loci with ≥500-bp TRs. Of note, 8909 loci with extended TRs (>100b longer than the mode) contain several known disease-associated TRs and are considered candidates for association with disorders. Overall, complex TRs and minisatellites are found to be abundant and diverse, even in genetically small Japanese populations, yielding insights into the landscape of long TRs.

Repeat polymorphisms underlie top genetic risk loci for glaucoma and colorectal cancer

Many regions in the human genome vary in length among individuals due to variable numbers of tandem repeats (VNTRs). To assess the phenotypic impact of VNTRs genome-wide, we applied a statistical imputation approach to estimate the lengths of 9,561 autosomal VNTR loci in 418,136 unrelated UK Biobank participants and 838 GTEx participants. Association and statistical fine-mapping analyses identified 58 VNTRs that appeared to influence a complex trait in UK Biobank, 18 of which also appeared to modulate expression or splicing of a nearby gene. Non-coding VNTRs at TMCO1 and EIF3H appeared to generate the largest known contributions of common human genetic variation to risk of glaucoma and colorectal cancer, respectively. Each of these two VNTRs associated with a >2-fold range of risk across individuals. These results reveal a substantial and previously unappreciated role of non-coding VNTRs in human health and gene regulation.

vamos: variable-number tandem repeats annotation using efficient motif sets

Roughly 3% of the human genome is composed of variable-number tandem repeats (VNTRs): arrays of motifs at least six bases. These loci are highly polymorphic, yet current approaches that define and merge variants based on alignment breakpoints do not capture their full diversity. Here we present a method vamos: VNTR Annotation using efficient Motif Sets that instead annotates VNTR using repeat composition under different levels of motif diversity. Using vamos we estimate 7.4-16.7 alleles per locus when applied to 74 haplotype-resolved human assemblies, compared to breakpoint-based approaches that estimate 4.0-5.5 alleles per locus.

Evaluating the Contribution of Cell Type-Specific Alternative Splicing to Variation in Lipid Levels

BACKGROUND

Genome-wide association studies have identified hundreds of loci associated with lipid levels. However, the genetic mechanisms underlying most of these loci are not well-understood. Recent work indicates that changes in the abundance of alternatively spliced transcripts contribute to complex trait variation. Consequently, identifying genetic loci that associate with alternative splicing in disease-relevant cell types and determining the degree to which these loci are informative for lipid biology is of broad interest.

METHODS

We analyze gene splicing in 83 sample-matched induced pluripotent stem cell (iPSC) and hepatocyte-like cell lines (n=166), as well as in an independent collection of primary liver tissues (n=96) to perform discovery of splicing quantitative trait loci (sQTLs).

RESULTS

We observe that transcript splicing is highly cell type specific, and the genes that are differentially spliced between iPSCs and hepatocyte-like cells are enriched for metabolism pathway annotations. We identify 1384 hepatocyte-like cell sQTLs and 1455 iPSC sQTLs at a false discovery rate of <5% and find that sQTLs are often shared across cell types. To evaluate the contribution of sQTLs to variation in lipid levels, we conduct colocalization analysis using lipid genome-wide association data. We identify 19 lipid-associated loci that colocalize either with an hepatocyte-like cell expression quantitative trait locus or sQTL. Only 2 loci colocalize with both a sQTL and expression quantitative trait locus, indicating that sQTLs contribute information about genome-wide association studies loci that cannot be obtained by analysis of steady-state gene expression alone.

CONCLUSIONS

These results provide an important foundation for future efforts that use iPSC and iPSC-derived cells to evaluate genetic mechanisms influencing both cardiovascular disease risk and complex traits in general.

Aberrant DNA Methylation, Expression, and Occurrence of Transcript Variants of the ABC Transporter ABCA7 in Breast Cancer

The ABC transporter ABCA7 has been found to be aberrantly expressed in a variety of cancer types, including breast cancer. We searched for specific epigenetic and genetic alterations and alternative splicing variants of ABCA7 in breast cancer and investigated whether these alterations are associated with ABCA7 expression. By analyzing tumor tissues from breast cancer patients, we found CpGs at the exon 5-intron 5 boundary aberrantly methylated in a molecular subtype-specific manner. The detection of altered DNA methylation in tumor-adjacent tissues suggests epigenetic field cancerization. In breast cancer cell lines, DNA methylation levels of CpGs in promoter-exon 1, intron 1, and at the exon 5-intron 5 boundary were not correlated with ABCA7 mRNA levels. By qPCR involving intron-specific and intron-flanking primers, we identified intron-containing ABCA7 mRNA transcripts. The occurrence of intron-containing transcripts was neither molecular subtype-specific nor directly correlated with DNA methylation at the respective exon-intron boundaries. Treatment of breast cancer cell lines MCF-7, BT-474, SK-BR3, and MDA-MB-231 with doxorubicin or paclitaxel for 72 h resulted in altered ABCA7 intron levels. Shotgun proteomics revealed that an increase in intron-containing transcripts was associated with significant dysregulation of splicing factors linked to alternative splicing.

A novel hypervariable variable number tandem repeat in the dopamine transporter gene (SLC6A3)

The dopamine transporter gene, SLC6A3, has received substantial attention in genetic association studies of various phenotypes. Although some variable number tandem repeats (VNTRs) present in SLC6A3 have been tested in genetic association studies, results have not been consistent. VNTRs in SLC6A3 that have not been examined genetically were characterized. The Tandem Repeat Annotation Library was used to characterize the VNTRs of 64 unrelated long-read haplotype-phased SLC6A3 sequences. Sequence similarity of each repeat unit of the five VNTRs is reported, along with the correlations of SNP-SNP, SNP-VNTR, and VNTR-VNTR alleles across the gene. One of these VNTRs is a novel hyper-VNTR (hyVNTR) in intron 8 of SLC6A3, which contains a range of 3.4-133.4 repeat copies and has a consensus sequence length of 38 bp, with 82% G+C content. The 38-base repeat was predicted to form G-quadruplexes in silico and was confirmed by circular dichroism spectroscopy. In addition, this hyVNTR contains multiple putative binding sites for PRDM9, which, in combination with low levels of linkage disequilibrium around the hyVNTR, suggests it might be a recombination hotspot.

Genome-wide identification of tandem repeats associated with splicing variation across 49 tissues in humans

Tandem repeats (TRs) are one of the largest sources of polymorphism, and their length is associated with gene regulation. Although previous studies reported several tandem repeats regulating gene splicing in cis (spl-TRs), no large-scale study has been conducted. In this study, we established a genome-wide catalog of 9537 spl-TRs with a total of 58,290 significant TR-splicing associations across 49 tissues (false discovery rate 5%) by using Genotype-Tissue expression (GTex) Project data. Regression models explaining splicing variation by using spl-TRs and other flanking variants suggest that at least some of the spl-TRs directly modulate splicing. In our catalog, two spl-TRs are known loci for repeat expansion diseases, spinocerebellar ataxia 6 (SCA6) and 12 (SCA12). Splicing alterations by these spl-TRs were compatible with those observed in SCA6 and SCA12. Thus, our comprehensive spl-TR catalog may help elucidate the pathomechanism of genetic diseases.

The role of structural variations in Alzheimer's disease and other neurodegenerative diseases

Dozens of single nucleotide polymorphisms (SNPs) related to Alzheimer's disease (AD) have been discovered by large scale genome-wide association studies (GWASs). However, only a small portion of the genetic component of AD can be explained by SNPs observed from GWAS. Structural variation (SV) can be a major contributor to the missing heritability of AD; while SV in AD remains largely unexplored as the accurate detection of SVs from the widely used array-based and short-read technology are still far from perfect. Here, we briefly summarized the strengths and weaknesses of available SV detection methods. We reviewed the current landscape of SV analysis in AD and SVs that have been found associated with AD. Particularly, the importance of currently less explored SVs, including insertions, inversions, short tandem repeats, and transposable elements in neurodegenerative diseases were highlighted.

Variable number tandem repeats (VNTRs) as modifiers of breast cancer risk in carriers of BRCA1 185delAG

Despite substantial efforts in identifying both rare and common variants affecting disease risk, in the majority of diseases, a large proportion of unexplained genetic risk remains. We propose that variable number tandem repeats (VNTRs) may explain a proportion of the missing genetic risk. Herein, in a pilot study with a retrospective cohort design, we tested whether VNTRs are causal modifiers of breast cancer risk in 347 female carriers of the BRCA1 185delAG pathogenic variant, an important group given their high risk of developing breast cancer. We performed targeted-capture to sequence VNTRs, called genotypes with adVNTR, tested the association of VNTRs and breast cancer risk using Cox regression models, and estimated the effect size using a retrospective likelihood approach. Of 303 VNTRs that passed quality control checks, 4 VNTRs were significantly associated with risk to develop breast cancer at false discovery rate [FDR] < 0.05 and an additional 4 VNTRs had FDR < 0.25. After determining the specific risk alleles, there was a significantly earlier age at diagnosis of breast cancer in carriers of the risk alleles compared to those without the risk alleles for seven of eight VNTRs. One example is a VNTR in exon 2 of LINC01973 with a per-allele hazard ratio of 1.58 (1.07-2.33) and 5.28 (2.79-9.99) for the homozygous risk-allele genotype. Results from this first systematic study of VNTRs demonstrate that VNTRs may explain a proportion of the unexplained genetic risk for breast cancer.

Using coding and non-coding rare variants to target candidate genes in patients with severe tinnitus

Tinnitus is the phantom percept of an internal non-verbal set of noises and tones. It is reported by 15% of the population and it is usually associated with hearing and/or brain disorders. The role of structural variants (SVs) in coding and non-coding regions has not been investigated in patients with severe tinnitus. In this study, we performed whole-genome sequencing in 97 unrelated Swedish individuals with chronic tinnitus (TIGER cohort). Rare single nucleotide variants (SNV), large structural variants (LSV), and copy number variations (CNV) were retrieved to perform a gene enrichment analysis in TIGER and in a subgroup of patients with severe tinnitus (SEVTIN, n = 34), according to the tinnitus handicap inventory (THI) scores. An independent exome sequencing dataset of 147 Swedish tinnitus patients was used as a replication cohort (JAGUAR cohort) and population-specific datasets from Sweden (SweGen) and Non-Finish Europeans (NFE) from gnomAD were used as control groups. SEVTIN patients showed a higher prevalence of hyperacusis, hearing loss, and anxiety when they were compared to individuals in the TIGER cohort. We found an enrichment of rare missense variants in 6 and 8 high-constraint genes in SEVTIN and TIGER cohorts, respectively. Of note, an enrichment of missense variants was found in the CACNA1E gene in both SEVTIN and TIGER. We replicated the burden of missense variants in 9 high-constrained genes in the JAGUAR cohort, including the gene NAV2, when data were compared with NFE. Moreover, LSVs in constrained regions overlapping CACNA1E, NAV2, and TMEM132D genes were observed in TIGER and SEVTIN.

SNP-by-CpG Site Interactions in ABCA7 Are Associated with Cognition in Older African Americans

SNPs in ABCA7 confer the largest genetic risk for Alzheimer's Disease (AD) in African Americans (AA) after APOE ε4. However, the relationship between ABCA7 and cognitive function has not been thoroughly examined. We investigated the effects of five known AD risk SNPs and 72 CpGs in ABCA7, as well as their interactions, on general cognitive function (cognition) in 634 older AA without dementia from Genetic Epidemiology Network of Arteriopathy (GENOA). Using linear mixed models, no SNP or CpG was associated with cognition after multiple testing correction, but five CpGs were nominally associated (p < 0.05). Four SNP-by-CpG interactions were associated with cognition (FDR q < 0.1). Contrast tests show that methylation is associated with cognition in some genotype groups (p < 0.05): a 1% increase at cg00135882 and cg22271697 is associated with a 0.68 SD decrease and 0.14 SD increase in cognition for those with the rs3764647 GG/AG (p = 0.004) and AA (p = 2 × 10-4) genotypes, respectively. In addition, a 1% increase at cg06169110 and cg17316918 is associated with a 0.37 SD decrease (p = 2 × 10-4) and 0.33 SD increase (p = 0.004), respectively, in cognition for those with the rs115550680 GG/AG genotype. While AD risk SNPs in ABCA7 were not associated with cognition in this sample, some have interactions with proximal methylation on cognition.

A polymorphic transcriptional regulatory domain in the amyotrophic lateral sclerosis risk gene CFAP410 correlates with differential isoform expression

We describe the characterisation of a variable number tandem repeat (VNTR) domain within intron 1 of the amyotrophic lateral sclerosis (ALS) risk gene CFAP410 (Cilia and flagella associated protein 410) (previously known as C21orf2), providing insight into how this domain could support differential gene expression and thus be a modulator of ALS progression or risk. We demonstrated the VNTR was functional in a reporter gene assay in the HEK293 cell line, exhibiting both the properties of an activator domain and a transcriptional start site, and that the differential expression was directed by distinct repeat number in the VNTR. These properties embedded in the VNTR demonstrated the potential for this VNTR to modulate CFAP410 expression. We extrapolated these findings in silico by utilisation of tagging SNPs for the two most common VNTR alleles to establish a correlation with endogenous gene expression. Consistent with in vitro data, CFAP410 isoform expression was found to be variable in the brain. Furthermore, although the number of matched controls was low, there was evidence for one specific isoform being correlated with lower expression in those with ALS. To address if the genotype of the VNTR was associated with ALS risk, we characterised the variation of the CFAP410 VNTR in ALS cases and matched controls by PCR analysis of the VNTR length, defining eight alleles of the VNTR. No significant difference was observed between cases and controls, we noted, however, the cohort was unlikely to contain sufficient power to enable any firm conclusion to be drawn from this analysis. This data demonstrated that the VNTR domain has the potential to modulate CFAP410 expression as a regulatory element that could play a role in its tissue-specific and stimulus-inducible regulation that could impact the mechanism by which CFAP410 is involved in ALS.

Pathogenesis, therapeutic strategies and biomarker development based on "omics" analysis related to microglia in Alzheimer's disease

Alzheimer's disease (AD) is the most common neurodegenerative disease and the most common cause of dementia. Among various pathophysiological aspects, microglia are considered to play important roles in the pathogenesis of AD. Genome wide association studies (GWAS) showed that the majority of AD risk genes are highly or exclusively expressed in microglia, underscoring the critical roles of microglia in AD pathogenesis. Recently, omics technologies have greatly advanced our knowledge of microglia biology in AD. Omics approaches, including genomics, epigenomics, transcriptomics, proteomics, and metabolomics/lipidomics, present remarkable opportunities to delineate the underlying mechanisms, discover novel diagnostic biomarkers, monitor disease progression, and shape therapeutic strategies for diseases. In this review, we summarized research based on microglial "omics" analysis in AD, especially the recent research advances in the identification of AD-associated microglial subsets. This review reinforces the important role of microglia in AD and advances our understanding of the mechanism of microglia in AD pathogenesis. Moreover, we proposed the value of microglia-based omics in the development of therapeutic strategies and biomarkers for AD.

Association between ABCA7 gene polymorphisms and Parkinson's disease susceptibility in a northern Chinese Han population

OBJECTIVE

As a typical member of the ABC transporter superfamily, ABCA7 has been shown to play an important role in stalling the pathogenesis of neurodegenerative disorders through maintaining the normal microglial function, regulating cellular responses to inflammation and ER stress, and modulating lipid metabolism. Variants in the ABCA7 locus have been hypothesized to be correlated with the genetic predisposition of several neurodegenerative disorders. The goal of this study was to examine whether there is a link between three specific single nucleotide polymorphisms in the ABCA7 gene, namely, rs3764650, rs4147929, and rs3752246, with the risk of developing Parkinson's disease (PD) in a northern Chinese Han community.

METHODS

In this case-control study, we recruited 821 participants, including 411 patients with sporadic PD and 410 independent, healthy controls. A Polymerase Chain Reaction-Restriction Fragment Length Polymorphism genotyping assay was used to identify polymorphisms of the three selected single nucleotide polymorphisms (rs3764650, rs4147929, and rs3752246) of the ABCA7 gene. Sanger sequencing was further applied to identify the accuracy of the genotyping results. The chi-square test was used to compare the frequencies of alleles and genotypes in patients and controls. Odds ratios and 95% confidence intervals were calculated using logistic regression.

RESULTS

We found significant between-group differences in the alleles (A vs. G, nominal P = 0.014) and dominant models (AA + GA vs. GG, nominal P = 0.015) of rs4147929. Subgroup analysis showed that the frequency of the rs4147929 A allele in male patients with PD was significantly higher than that in male controls (nominal P = 0.036). For the rs3752246 polymorphism, the frequency of the G allele was significantly higher in patients with PD than in controls, and the dominant model fit the data best when considering the nominal P-values (nominal P = 0.019, nominal P = 0.033, respectively). Differences in G allele and genotypes frequencies between patients and controls remained significant in women (nominal P = 0.032 for allele, nominal P = 0.015 for genotype), as well as in individuals aged more than 50 years (nominal P = 0.044, nominal P = 0.020, respectively). No significant differences were observed in allele or genotype frequencies between patients with PD and healthy controls for rs3764650. The frequency of the TCG (rs3764650-rs3752246-rs4147929) haplotype was significantly lower in the PD group than in the healthy control group (odds ratio = 0.772; 95% confidence interval = 0.634-0.940; P = 0.011).

CONCLUSION

The rs4147929 polymorphism was significantly associated with PD susceptibility in the northern Chinese Han population. The A allele of rs4147929 was a risk factor for developing PD. The TCG haplotype presented a protective role in the pathogenesis of PD. Further studies using larger sample sizes, considering different clinical and biochemical parameters such as the cognitive status of subjects at the same time, are warranted to better clarify the effects of these common variants on the pathogenesis and development of PD.

A phenome-wide association study identifies effects of copy-number variation of VNTRs and multicopy genes on multiple human traits

The human genome contains tens of thousands of large tandem repeats and hundreds of genes that show common and highly variable copy-number changes. Due to their large size and repetitive nature, these variable number tandem repeats (VNTRs) and multicopy genes are generally recalcitrant to standard genotyping approaches and, as a result, this class of variation is poorly characterized. However, several recent studies have demonstrated that copy-number variation of VNTRs can modify local gene expression, epigenetics, and human traits, indicating that many have a functional role. Here, using read depth from whole-genome sequencing to profile copy number, we report results of a phenome-wide association study (PheWAS) of VNTRs and multicopy genes in a discovery cohort of ∼35,000 samples, identifying 32 traits associated with copy number of 38 VNTRs and multicopy genes at 1% FDR. We replicated many of these signals in an independent cohort and observed that VNTRs showing trait associations were significantly enriched for expression QTLs with nearby genes, providing strong support for our results. Fine-mapping studies indicated that in the majority (∼90%) of cases, the VNTRs and multicopy genes we identified represent the causal variants underlying the observed associations. Furthermore, several lie in regions where prior SNV-based GWASs have failed to identify any significant associations with these traits. Our study indicates that copy number of VNTRs and multicopy genes contributes to diverse human traits and suggests that complex structural variants potentially explain some of the so-called "missing heritability" of SNV-based GWASs.

The role of ATP-binding cassette subfamily A in the etiology of Alzheimer's disease

Background

Alzheimer’s disease (AD) is the leading cause of dementia, clinically characterized by memory deficits and progressive cognitive decline. Despite decades of research effective therapies are lacking, and a large part of the genetic heritability remains unidentified. ABCA7 and ABCA1, members of the ATP-binding cassette subfamily A (ABCA), were identified as AD risk genes in genome-wide association studies. Nevertheless, genetic and/or functional studies propose a link between AD and two other members of the ABCA subclass, i.e., ABCA2 and ABCA5.

Main body

Changes in expression or dysfunction of these transporters were found to increase amyloid β levels. This might be related to the common role of ABCA transporters in cellular cholesterol homeostasis, for which a prominent role in AD development has been suggested. In this review, we provide a comprehensive overview and discussion on the contribution of the ABCA subfamily to the etiopathogenesis of AD.

Conclusions

A better understanding of the function and identification of disease-associated genetic variants in ABCA transporters can contribute to the development of novel therapeutic strategies for AD.

Rare missense mutations in ABCA7 might increase Alzheimer's disease risk by plasma membrane exclusion

The adenosine triphosphate-binding cassette subfamily A member 7 gene (ABCA7) is associated with Alzheimer's disease (AD) in large genome-wide association studies. Targeted sequencing of ABCA7 suggests a role for rare premature termination codon (PTC) mutations in AD, with haploinsufficiency through nonsense-mediated mRNA decay as a plausible pathogenic mechanism. Since other classes of rare variants in ABCA7 are poorly understood, we investigated the contribution and pathogenicity of rare missense, indel and splice variants in ABCA7 in Belgian AD patient and control cohorts. We identified 8.36% rare variants in the patient cohort versus 6.05% in the control cohort. For 10 missense mutations identified in the Belgian cohort we analyzed the pathogenetic effect on protein localization in vitro using immunocytochemistry. Our results demonstrate that rare ABCA7 missense mutations can contribute to AD by inducing protein mislocalization, resulting in a lack of functional protein at the plasma membrane. In one pedigree, a mislocalization-inducing missense mutation in ABCA7 (p.G1820S) co-segregated with AD in an autosomal dominant inheritance pattern. Brain autopsy of six patient missense mutation carriers showed typical AD neuropathological characteristics including cerebral amyloid angiopathy type 1. Also, among the rare ABCA7 missense mutations, we observed mutations that affect amino acid residues that are conserved in ABCA1 and ABCA4, of which some correspond to established ABCA1 or ABCA4 disease-causing mutations involved in Tangier or Stargardt disease.

trfermikit: a tool to discover VNTR-associated deletions

SUMMARY

We present trfermikit, a software tool designed to detect deletions larger than 50 bp occurring in Variable Number Tandem Repeats using Illumina DNA sequencing reads. In such regions, it achieves a better tradeoff between sensitivity and false discovery than a state-of-the-art structural variation caller, Manta and complements it by recovering a significant number of deletions that Manta missed. trfermikit is based upon the fermikit pipeline, which performs read assembly, maps the assembly to the reference genome and calls variants from the alignment.

AVAILABILITY AND IMPLEMENTATION

https://github.com/petermchale/trfermikit.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Pharmacogenomics of Alzheimer's Disease: Novel Strategies for Drug Utilization and Development

Alzheimer's disease (AD) is a priority health problem in developed countries with a high cost to society. Approximately 20% of direct costs are associated with pharmacological treatment. Over 90% of patients require multifactorial treatments, with risk of adverse drug reactions (ADRs) and drug-drug interactions (DDIs) for the treatment of concomitant diseases such as hypertension (>25%), obesity (>70%), diabetes mellitus type 2 (>25%), hypercholesterolemia (40%), hypertriglyceridemia (20%), metabolic syndrome (20%), hepatobiliary disorder (15%), endocrine/metabolic disorders (>20%), cardiovascular disorder (40%), cerebrovascular disorder (60-90%), neuropsychiatric disorders (60-90%), and cancer (10%).For the past decades, pharmacological studies in search of potential treatments for AD focused on the following categories: neurotransmitter enhancers (11.38%), multitarget drugs (2.45%), anti-amyloid agents (13.30%), anti-tau agents (2.03%), natural products and derivatives (25.58%), novel synthetic drugs (8.13%), novel targets (5.66%), repository drugs (11.77%), anti-inflammatory drugs (1.20%), neuroprotective peptides (1.25%), stem cell therapy (1.85%), nanocarriers/nanotherapeutics (1.52%), and other compounds (<1%).Pharmacogenetic studies have shown that the therapeutic response to drugs in AD is genotype-specific in close association with the gene clusters that constitute the pharmacogenetic machinery (pathogenic, mechanistic, metabolic, transporter, pleiotropic genes) under the regulatory control of epigenetic mechanisms (DNA methylation, histone/chromatin remodeling, microRNA regulation). Most AD patients (>60%) are carriers of over ten pathogenic genes. The genes that most frequently (>50%) accumulate pathogenic variants in the same AD case are A2M (54.38%), ACE (78.94%), BIN1 (57.89%), CLU (63.15%), CPZ (63.15%), LHFPL6 (52.63%), MS4A4E (50.87%), MS4A6A (63.15%), PICALM (54.38%), PRNP (80.7059), and PSEN1 (77.19%). There is also an accumulation of 15 to 26 defective pharmagenes in approximately 85% of AD patients. About 50% of AD patients are carriers of at least 20 mutant pharmagenes, and over 80% are deficient metabolizers for the most common drugs, which are metabolized via the CYP2D6, CYP2C9, CYP2C19, and CYP3A4/5 enzymes.The implementation of pharmacogenetics can help optimize drug development and the limited therapeutic resources available to treat AD, and personalize the use of anti-dementia drugs in combination with other medications for the treatment of concomitant disorders.

Advancing genomic technologies and clinical awareness accelerates discovery of disease-associated tandem repeat sequences

Expansions of gene-specific DNA tandem repeats (TRs), first described in 1991 as a disease-causing mutation in humans, are now known to cause >60 phenotypes, not just disease, and not only in humans. TRs are a common form of genetic variation with biological consequences, observed, so far, in humans, dogs, plants, oysters, and yeast. Repeat diseases show atypical clinical features, genetic anticipation, and multiple and partially penetrant phenotypes among family members. Discovery of disease-causing repeat expansion loci accelerated through technological advances in DNA sequencing and computational analyses. Between 2019 and 2021, 17 new disease-causing TR expansions were reported, totaling 63 TR loci (>69 diseases), with a likelihood of more discoveries, and in more organisms. Recent and historical lessons reveal that properly assessed clinical presentations, coupled with genetic and biological awareness, can guide discovery of disease-causing unstable TRs. We highlight critical but underrecognized aspects of TR mutations. Repeat motifs may not be present in current reference genomes but will be in forthcoming gapless long-read references. Repeat motif size can be a single nucleotide to kilobases/unit. At a given locus, repeat motif sequence purity can vary with consequence. Pathogenic repeats can be "insertions" within nonpathogenic TRs. Expansions, contractions, and somatic length variations of TRs can have clinical/biological consequences. TR instabilities occur in humans and other organisms. TRs can be epigenetically modified and/or chromosomal fragile sites. We discuss the expanding field of disease-associated TR instabilities, highlighting prospects, clinical and genetic clues, tools, and challenges for further discoveries of disease-causing TR instabilities and understanding their biological and pathological impacts-a vista that is about to expand.

Revisiting tandem repeats in psychiatric disorders from perspectives of genetics, physiology, and brain evolution

Genome-wide association studies (GWASs) have revealed substantial genetic components comprised of single nucleotide polymorphisms (SNPs) in the heritable risk of psychiatric disorders. However, genetic risk factors not covered by GWAS also play pivotal roles in these illnesses. Tandem repeats, which are likely functional but frequently overlooked by GWAS, may account for an important proportion in the "missing heritability" of psychiatric disorders. Despite difficulties in characterizing and quantifying tandem repeats in the genome, studies have been carried out in an attempt to describe impact of tandem repeats on gene regulation and human phenotypes. In this review, we have introduced recent research progress regarding the genomic distribution and regulatory mechanisms of tandem repeats. We have also summarized the current knowledge of the genetic architecture and biological underpinnings of psychiatric disorders brought by studies of tandem repeats. These findings suggest that tandem repeats, in candidate psychiatric risk genes or in different levels of linkage disequilibrium (LD) with psychiatric GWAS SNPs and haplotypes, may modulate biological phenotypes related to psychiatric disorders (e.g., cognitive function and brain physiology) through regulating alternative splicing, promoter activity, enhancer activity and so on. In addition, many tandem repeats undergo tight natural selection in the human lineage, and likely exert crucial roles in human brain evolution. Taken together, the putative roles of tandem repeats in the pathogenesis of psychiatric disorders is strongly implicated, and using examples from previous literatures, we wish to call for further attention to tandem repeats in the post-GWAS era of psychiatric disorders.

Sweet Taste Preference: Relationships with Other Tastes, Liking for Sugary Foods and Exploratory Genome-Wide Association Analysis in Subjects with Metabolic Syndrome

Taste perception and its association with nutrition and related diseases (type 2 diabetes, obesity, metabolic syndrome, cardiovascular, etc.) are emerging fields of biomedicine. There is currently great interest in investigating the environmental and genetic factors that influence sweet taste and sugary food preferences for personalized nutrition. Our aims were: (1) to carry out an integrated analysis of the influence of sweet taste preference (both in isolation and in the context of other tastes) on the preference for sugary foods and its modulation by type 2 diabetes status; (2) as well as to explore new genetic factors associated with sweet taste preference. We studied 425 elderly white European subjects with metabolic syndrome and analyzed taste preference, taste perception, sugary-foods liking, biochemical and genetic markers. We found that type 2 diabetic subjects (38%) have a small, but statistically higher preference for sweet taste (p = 0.021) than non-diabetic subjects. No statistically significant differences (p > 0.05) in preferences for the other tastes (bitter, salty, sour or umami) were detected. For taste perception, type 2 diabetic subjects have a slightly lower perception of all tastes (p = 0.026 for the combined “total taste score”), bitter taste being statistically lower (p = 0.023). We also carried out a principal component analysis (PCA), to identify latent variables related to preferences for the five tastes. We identified two factors with eigenvalues >1. Factor 2 was the one with the highest correlation with sweet taste preference. Sweet taste preference was strongly associated with a liking for sugary foods. In the exploratory SNP-based genome-wide association study (GWAS), we identified some SNPs associated with sweet taste preference, both at the suggestive and at the genome-wide level, especially a lead SNP in the PTPRN2 (Protein Tyrosine Phosphatase Receptor Type N2) gene, whose minor allele was associated with a lower sweet taste preference. The PTPRN2 gene was also a top-ranked gene obtained in the gene-based exploratory GWAS analysis. In conclusion, sweet taste preference was strongly associated with sugary food liking in this population. Our exploratory GWAS identified an interesting candidate gene related with sweet taste preference, but more studies in other populations are required for personalized nutrition.

Uncovering the impact of noncoding variants in neurodegenerative brain diseases

Neurodegenerative brain diseases (NBDs) are characterized by cognitive decline and movement impairments caused by neuronal loss in different brain regions. A large fraction of the genetic heritability of NBDs is not explained by the current known mutations. Genome-wide association studies identified novel disease-risk loci, adding to the genetic basis of NBDs. Many of the associated variants reside in noncoding regions with distinct molecular functions. Genetic variation in these regions can alter functions and contribute to disease pathogenesis. Here, we discuss noncoding variants associated with NBDs. Methods for better functional interpretation of noncoding variation will expand our knowledge of the genetic architecture of NBDs and broaden the routes for therapeutic strategies.

Straglr: discovering and genotyping tandem repeat expansions using whole genome long-read sequences

Tandem repeat (TR) expansion is the underlying cause of over 40 neurological disorders. Long-read sequencing offers an exciting avenue over conventional technologies for detecting TR expansions. Here, we present Straglr, a robust software tool for both targeted genotyping and novel expansion detection from long-read alignments. We benchmark Straglr using various simulations, targeted genotyping data of cell lines carrying expansions of known diseases, and whole genome sequencing data with chromosome-scale assembly. Our results suggest that Straglr may be useful for investigating disease-associated TR expansions using long-read sequencing.

Characterizing nucleotide variation and expansion dynamics in human-specific variable number tandem repeats

There are more than 55,000 variable number tandem repeats (VNTRs) in the human genome, notable for both their striking polymorphism and mutability. Despite their role in human evolution and genomic variation, they have yet to be studied collectively and in detail, partially owing to their large size, variability, and predominant location in noncoding regions. Here, we examine 467 VNTRs that are human-specific expansions, unique to one location in the genome, and not associated with retrotransposons. We leverage publicly available long-read genomes, including from the Human Genome Structural Variant Consortium, to ascertain the exact nucleotide composition of these VNTRs and compare their composition of alleles. We then confirm repeat unit composition in more than 3000 short-read samples from the 1000 Genomes Project. Our analysis reveals that these VNTRs contain highly structured repeat motif organization, modified by frequent deletion and duplication events. Although overall VNTR compositions tend to remain similar between 1000 Genomes Project superpopulations, we describe a notable exception with substantial differences in repeat composition (in PCBP3), as well as several VNTRs that are significantly different in length between superpopulations (in ART1, PROP1, DYNC2I1, and LOC102723906). We also observe that most of these VNTRs are expanded in archaic human genomes, yet remain stable in length between single generations. Collectively, our findings indicate that repeat motif variability, repeat composition, and repeat length are all informative modalities to consider when characterizing VNTRs and their contribution to genomic variation.

Brain Regional Differences in Hexanucleotide Repeat Length in X-Linked Dystonia-Parkinsonism Using Nanopore Sequencing

Objective

Our study investigated the presence of regional differences in hexanucleotide repeat number in postmortem brain tissues of 2 patients with X-linked dystonia-parkinsonism (XDP), a combined dystonia-parkinsonism syndrome modified by a (CCCTCT)_n repeat within the causal SINE-VNTR-Alu retrotransposon insertion in the TAF1 gene.

Methods

Genomic DNA was extracted from blood and postmortem brain samples, including the basal ganglia and cortex from both patients and from the cerebellum, midbrain, and pituitary gland from 1 patient. Repeat sizing was performed using fragment analysis, small-pool PCR-based Southern blotting, and Oxford nanopore sequencing.

Results

The basal ganglia (p < 0.001) and cerebellum (p < 0.001) showed higher median repeat numbers and higher degrees of repeat instability compared with blood.

Conclusions

Somatic repeat instability may predominate in brain regions selectively affected in XDP, thereby hinting at its potential role in disease manifestation and modification.

Genome-wide characterization of human minisatellite VNTRs: population-specific alleles and gene expression differences

Variable Number Tandem Repeats (VNTRs) are tandem repeat (TR) loci that vary in copy number across a population. Using our program, VNTRseek, we analyzed human whole genome sequencing datasets from 2770 individuals in order to detect minisatellite VNTRs, i.e., those with pattern sizes ≥7 bp. We detected 35 638 VNTR loci and classified 5676 as commonly polymorphic (i.e. with non-reference alleles occurring in >5% of the population). Commonly polymorphic VNTR loci were found to be enriched in genomic regions with regulatory function, i.e. transcription start sites and enhancers. Investigation of the commonly polymorphic VNTRs in the context of population ancestry revealed that 1096 loci contained population-specific alleles and that those could be used to classify individuals into super-populations with near-perfect accuracy. Search for quantitative trait loci (eQTLs), among the VNTRs proximal to genes, indicated that in 187 genes expression differences correlated with VNTR genotype. We validated our predictions in several ways, including experimentally, through the identification of predicted alleles in long reads, and by comparisons showing consistency between sequencing platforms. This study is the most comprehensive analysis of minisatellite VNTRs in the human population to date.

30 years of repeat expansion disorders: What have we learned and what are the remaining challenges?

Tandem repeats represent one of the most abundant class of variations in human genomes, which are polymorphic by nature and become highly unstable in a length-dependent manner. The expansion of repeat length across generations is a well-established process that results in human disorders mainly affecting the central nervous system. At least 50 disorders associated with expansion loci have been described to date, with half recognized only in the last ten years, as prior methodological difficulties limited their identification. These limitations still apply to the current widely used molecular diagnostic methods (exome or gene panels) and thus result in missed diagnosis detrimental to affected individuals and their families, especially for disorders that are very rare and/or clinically not recognizable. Most of these disorders have been identified through family-driven approaches and many others likely remain to be identified. The recent development of long-read technologies provides a unique opportunity to systematically investigate the contribution of tandem repeats and repeat expansions to the genetic architecture of human disorders. In this review, we summarize the current and most recent knowledge about the genetics of repeat expansion disorders and the diversity of their pathophysiological mechanisms and outline the perspectives of developing personalized treatments in the future.

Pervasive cis effects of variation in copy number of large tandem repeats on local DNA methylation and gene expression

Variable number tandem repeats (VNTRs) are composed of large tandemly repeated motifs, many of which are highly polymorphic in copy number. However, because of their large size and repetitive nature, they remain poorly studied. To investigate the regulatory potential of VNTRs, we used read-depth data from Illumina whole-genome sequencing to perform association analysis between copy number of ∼70,000 VNTRs (motif size ≥ 10 bp) with both gene expression (404 samples in 48 tissues) and DNA methylation (235 samples in peripheral blood), identifying thousands of VNTRs that are associated with local gene expression (eVNTRs) and DNA methylation levels (mVNTRs). Using an independent cohort, we validated 73%-80% of signals observed in the two discovery cohorts, while allelic analysis of VNTR length and CpG methylation in 30 Oxford Nanopore genomes gave additional support for mVNTR loci, thus providing robust evidence to support that these represent genuine associations. Further, conditional analysis indicated that many eVNTRs and mVNTRs act as QTLs independently of other local variation. We also observed strong enrichments of eVNTRs and mVNTRs for regulatory features such as enhancers and promoters. Using the Human Genome Diversity Panel, we define sets of VNTRs that show highly divergent copy numbers among human populations and show that these are enriched for regulatory effects and preferentially associate with genes that have been linked with human phenotypes through GWASs. Our study provides strong evidence supporting functional variation at thousands of VNTRs and defines candidate sets of VNTRs, copy number variation of which potentially plays a role in numerous human phenotypes.

Premature termination codon mutations in ABCA7 contribute to Alzheimer's disease risk in Belgian patients

The ATP-Binding Cassette Subfamily A Member 7 gene (ABCA7) was identified as a risk gene for Alzheimer's disease (AD) in genome-wide association studies of large cohorts of late-onset AD (LOAD) patients. Extended resequencing of the ABCA7 coding regions identified mutations that lead to premature termination codons (PTC) and loss of function of ABCA7. PTC mutations were enriched in LOAD patients and were frequently present in patients with early-onset AD (EOAD). We aimed at assessing the contribution of ABCA7 PTC mutations to AD in the Belgian population by screening the ABCA7 coding region in a Belgian AD cohort of 1376 patients, including LOAD and EOAD patients, and in a Belgian control cohort of 976 individuals. We identified a PTC mutation in 67 AD patients (4.9%) and in 18 control individuals (1.8%) confirming the enrichment of ABCA7 PTC mutations in Belgian AD patients. The patient carriers had a mean onset age of 69.7 ± 9.8 years with a wide onset age range of 42 years (48-90 years). In 77.3% of the families of ABCA7 carriers, there were AD patients present suggestive of a positive family history of disease, but a Mendelian co-segregation of ABCA7 PTC mutations with disease is not clear. Overall, our genetic data predict that PTC mutations in ABCA7 are common in the Belgian population and are present in LOAD and EOAD patients.