Twist exome capture allows for lower average sequence coverage in clinical exome sequencing

in Human Genomics (2023), 17(1), 39

Background Exome and genome sequencing are the predominant techniques in the diagnosis and research of genetic disorders. Sufficient, uniform and reproducible/consistent sequence coverage is a main ... [more ▼]

Background Exome and genome sequencing are the predominant techniques in the diagnosis and research of genetic disorders. Sufficient, uniform and reproducible/consistent sequence coverage is a main determinant for the sensitivity to detect single-nucleotide (SNVs) and copy number variants (CNVs). Here we compared the ability to obtain comprehensive exome coverage for recent exome capture kits and genome sequencing techniques. Results We compared three different widely used enrichment kits (Agilent SureSelect Human All Exon V5, Agilent SureSelect Human All Exon V7 and Twist Bioscience) as well as short-read and long-read WGS. We show that the Twist exome capture significantly improves complete coverage and coverage uniformity across coding regions compared to other exome capture kits. Twist performance is comparable to that of both short- and long-read whole genome sequencing. Additionally, we show that even at a reduced average coverage of 70× there is only minimal loss in sensitivity for SNV and CNV detection. Conclusion We conclude that exome sequencing with Twist represents a significant improvement and could be performed at lower sequence coverage compared to other exome capture techniques. [less ▲]

Accurate long-read sequencing identified GBA variants as a major genetic risk factor in the Luxembourg Parkinson\textquoterights study 2023.03.29.23287880

E-print/Working paper (2023)

Heterozygous variants in the glucocerebrosidase GBA gene are an increasingly recognized risk factor for Parkinson's disease (PD). Due to the pseudogene GBAP1 that shares 96\% sequence homology with the ... [more ▼]

Heterozygous variants in the glucocerebrosidase GBA gene are an increasingly recognized risk factor for Parkinson's disease (PD). Due to the pseudogene GBAP1 that shares 96\% sequence homology with the GBA coding region, accurate variant calling by array-based or short-read sequencing methods remains a major challenge in understanding the genetic landscape of GBA-related PD. We established a novel long-read sequencing technology for assessing the full length of the GBA gene. We used subsequent regression models for genotype-phenotype analyses. We sequenced 752 patients with parkinsonism and 806 healthy controls of the Luxembourg Parkinson's study. All GBA variants identified showed a 100% true positive rate by Sanger validation. We found 12% of unrelated PD patients carrying GBA variants. Three novel variants of unknown significance (VUS) were identified. Using a structure-based approach, we defined a potential risk prediction method for VUS. This study describes the full landscape of GBA-related parkinsonism in Luxembourg, showing a high prevalence of GBA variants as the major genetic risk for PD. Our approach provides an important advancement for highly accurate GBA variant calling, which is essential for providing access to emerging causative therapies for GBA carriers. [less ▲]

Transferability of European-derived cardiometabolic polygenic risk scores in the South Asians and their interplay with family history 2023.03.20.23287470

E-print/Working paper (2023)

Background & Aims We aimed to investigate the effect of polygenic risk scores (PRSs) derived from individuals of European (EUR) ancestry on common diseases among individuals of South Asian (SAS) ancestry ... [more ▼]

Background & Aims We aimed to investigate the effect of polygenic risk scores (PRSs) derived from individuals of European (EUR) ancestry on common diseases among individuals of South Asian (SAS) ancestry in the UK Biobank (UKB). Additionally, we studied the interaction between PRS and family history (FH) in the same population.Methods To calculate the PRS, we used a previously published panel of SNPs derived from the EUR population and applied it to the individuals of SAS ancestry from the UKB study. We applied the PRS using summary statistics from genome-wide association studies (GWAS) for cardiometabolic and lifestyle diseases such as coronary artery disease (CAD), obesity, and type 2 diabetes (T2D). Each PRS was adjusted according to an individual\textquoterights predicted genetic ancestry to derive an adjusted PRS (aPRS). We calculated the percentiles based on aPRS and divided them according to the percentiles into three categories: low, intermediate, and high. Considering the intermediate-aPRS percentile as a reference, we compared the low and high aPRS categories and generated the odds ratio (OR) estimates.Results The risk of developing severe obesity for individuals of SAS ancestry was almost threefold higher for individuals with high aPRS than for those with intermediate aPRS, with an OR of 3.67 (95% CI = 2.47-5.48, P < 0.01). While the risk of severe obesity was lower in the low-aPRS group (OR = 0.19, CI = 0.05\textendash0.52, P < 0.01). Comparable results were found in the EUR data, where the low-PRS group had an OR of 0.26 (95% CI= 0.24-0.3, P < 0.01) and the high-PRS group had an OR of 3.2 (95% CI = 3.1-3.3, P < 0.01). We observed similar results for CAD and T2D. Further, we show that SAS individuals with a familial history of CAD and T2D with high-aPRS exhibit further higher risk to these diseases, thereby implying a greater genetic predisposition to these conditions.Conclusion Our findings suggest that using CAD, obesity, and T2D GWAS summary statistics predominantly from the EUR population have sufficient power to identify SAS individuals with higher genetic risk. With future GWAS recruiting more SAS participants and tailoring the PRSs towards SAS ancestry, we believe that the predictive power of PRS would improve. [less ▲]

Forensic DNA Phenotyping: a review on SNP panels, genotyping techniques, and prediction models

E-print/Working paper (2023)

In the past few years, forensic DNA phenotyping (FDP) has attracted a strong interest in the forensic research. Among the increasing publications, many have focused on testing the available panels to ... [more ▼]

In the past few years, forensic DNA phenotyping (FDP) has attracted a strong interest in the forensic research. Among the increasing publications, many have focused on testing the available panels to infer biogeographical ancestry (BGA) on less represented populations and understanding the genetic mechanisms underlying externally visible characteristics (EVC). However, there are currently no publications that gather all the existing panels limited to FDP and discuss the main technical limitations of the technique. In this review, we performed a bibliographic search in Scopus database of FDP-related literature, which resulted in a total of 48, 43 and 15 panels for BGA, EVC and both BGA and EVC inference, respectively. Here we provide a list of commercial and non-commercial FDP panels and the FDP limitations regarding the lack of harmonization in terms of terminology (i.e., categorization and measurement of traits) and reporting, the lack of genetic knowledge and environment influence to select markers and develop panels, and the debate surrounding the selection of genotyping technologies and prediction models and algorithms. In conclusion, this review aims to be an updated guide and to present an overview of the current FDP-related literature. [less ▲]

Investigation of Shared Genetic Risk Factors Between Parkinson's Disease and Cancers

in Movement Disorders (2023), n/a(n/a),

Abstract Background Epidemiological studies that examined the association between Parkinson's disease (PD) and cancers led to inconsistent results, but they face a number of methodological difficulties ... [more ▼]

Abstract Background Epidemiological studies that examined the association between Parkinson's disease (PD) and cancers led to inconsistent results, but they face a number of methodological difficulties. Objective We used results from genome-wide association studies (GWASs) to study the genetic correlation between PD and different cancers to identify common genetic risk factors. Methods We used individual data for participants of European ancestry from the Courage-PD (Comprehensive Unbiased Risk Factor Assessment for Genetics and Environment in Parkinson's Disease; PD, N = 16,519) and EPITHYR (differentiated thyroid cancer, N = 3527) consortia and summary statistics of GWASs from iPDGC (International Parkinson Disease Genomics Consortium; PD, N = 482,730), Melanoma Meta-Analysis Consortium (MMAC), Breast Cancer Association Consortium (breast cancer), the Prostate Cancer Association Group to Investigate Cancer Associated Alterations in the Genome (prostate cancer), International Lung Cancer Consortium (lung cancer), and Ovarian Cancer Association Consortium (ovarian cancer) (N comprised between 36,017 and 228,951 for cancer GWASs). We estimated the genetic correlation between PD and cancers using linkage disequilibrium score regression. We studied the association between PD and polymorphisms associated with cancers, and vice versa, using cross-phenotypes polygenic risk score (PRS) analyses. Results We confirmed a previously reported positive genetic correlation of PD with melanoma (Gcorr = 0.16 [0.04; 0.28]) and reported an additional significant positive correlation of PD with prostate cancer (Gcorr = 0.11 [0.03; 0.19]). There was a significant inverse association between the PRS for ovarian cancer and PD (odds ratio [OR] = 0.89 [0.84; 0.94]). Conversely, the PRS of PD was positively associated with breast cancer (OR = 1.08 [1.06; 1.10]) and inversely associated with ovarian cancer (OR = 0.95 [0.91; 0.99]). The association between PD and ovarian cancer was mostly driven by rs183211 located in an intron of the NSF gene (17q21.31). Conclusions We show evidence in favor of a contribution of pleiotropic genes to the association between PD and specific cancers. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA. [less ▲]

Embracing Monogenic Parkinson's Disease: The MJFF Global Genetic PD Cohort

in Movement Disorders (2023), 38(2), 286--303

BACKGROUND: As gene-targeted therapies are increasingly being developed for Parkinson's disease (PD), identifying and characterizing carriers of specific genetic pathogenic variants is imperative. Only a ... [more ▼]

BACKGROUND: As gene-targeted therapies are increasingly being developed for Parkinson's disease (PD), identifying and characterizing carriers of specific genetic pathogenic variants is imperative. Only a small fraction of the estimated number of subjects with monogenic PD worldwide are currently represented in the literature and availability of clinical data and clinical trial-ready cohorts is limited. OBJECTIVE: The objectives are to (1) establish an international cohort of affected and unaffected individuals with PD-linked variants; (2) provide harmonized and quality-controlled clinical characterization data for each included individual; and (3) further promote collaboration of researchers in the field of monogenic PD. METHODS: We conducted a worldwide, systematic online survey to collect individual-level data on individuals with PD-linked variants in SNCA, LRRK2 VPS35, PRKN, PINK1, DJ-1, as well as selected pathogenic and risk variants in GBA and corresponding demographic, clinical, and genetic data. All registered cases underwent thorough quality checks, and pathogenicity scoring of the variants and genotype-phenotype relationships were analyzed. RESULTS: We collected 3888 variant carriers for our analyses, reported by 92 centers (42 countries) worldwide. Of the included individuals 3185 had a diagnosis of PD (ie, 1306 LRRK2, 115 SNCA, 23 VPS35 429 PRKN, 75 PINK1, 13 DJ-1, and 1224 GBA) and 703 were unaffected (ie, 328 LRRK2, 32 SNCA, 3 VPS35, 1 PRKN, 1 PINK1, and 338 GBA). In total, we identified 269 different pathogenic variants; 1322 individuals in our cohort (34\%) were indicated as not previously published. CONCLUSIONS: Within the MJFF Global Genetic PD Study Group, we (1) established the largest international cohort of affected and unaffected individuals carrying PD-linked variants; (2) provide harmonized and quality-controlled clinical and genetic data for each included individual; (3) promote collaboration in the field of genetic PD with a view toward clinical and genetic stratification of patients for gene-targeted clinical trials. 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society. [less ▲]

Data-driven historical characterization of epilepsy-associated genes.

in European Journal of Paediatric Neurology (2022)

Many epilepsy-associated genes have been identified over the last three decades, revealing a remarkable molecular heterogeneity with the shared outcome of recurrent seizures. Information about the genetic ... [more ▼]

Many epilepsy-associated genes have been identified over the last three decades, revealing a remarkable molecular heterogeneity with the shared outcome of recurrent seizures. Information about the genetic landscape of epilepsies is scattered throughout the literature and answering the simple question of how many genes are associated with epilepsy is not straightforward. Here, we present a computationally driven analytical review of epilepsy-associated genes using the complete scientific literature in PubMed. Based on our search criteria, we identified a total of 738 epilepsy-associated genes. We further classified these genes into two Tiers. A broad gene list of 738 epilepsy-associated genes (Tier 2) and a narrow gene list composed of 143 epilepsy-associated genes (Tier 1). Our search criteria do not reflect the degree of association. The average yearly number of identified epilepsy-associated genes between 1992 and 2021 was 4.8. However, most of these genes were only identified in the last decade (2010–2019). Ion channels represent the largest class of epilepsy-associated genes. For many of these, both gain- and loss-of-function effects have been associated with epilepsy in recent years. We identify 28 genes frequently reported with heterogenous variant effects which should be considered for variant interpretation. Overall, our study provides an updated and manually curated list of epilepsy-related genes together with additional annotations and classifications reflecting the current genetic landscape of epilepsy. [less ▲]

Alterations of oral microbiota and impact on the gut microbiome in type 1 diabetes mellitus revealed by integrated multi-omic analyses

in Microbiome (2022)

Background: Alterations to the gut microbiome have been linked to multiple chronic diseases. However, the drivers of such changes remain largely unknown. The oral cavity acts as a major route of exposure ... [more ▼]

Background: Alterations to the gut microbiome have been linked to multiple chronic diseases. However, the drivers of such changes remain largely unknown. The oral cavity acts as a major route of exposure to exogenous factors including pathogens, and processes therein may affect the communities in the subsequent compartments of the gastrointestinal tract. Here, we perform strain‑resolved, integrated meta‑genomic, transcriptomic, and proteomic analyses of paired saliva and stool samples collected from 35 individuals from eight families with multiple cases of type 1 diabetes mellitus (T1DM). Results: We identified distinct oral microbiota mostly reflecting competition between streptococcal species. More specifically, we found a decreased abundance of the commensal Streptococcus salivarius in the oral cavity of T1DM individuals, which is linked to its apparent competition with the pathobiont Streptococcus mutans. The decrease in S. salivarius in the oral cavity was also associated with its decrease in the gut as well as higher abundances in facultative anaerobes including Enterobacteria. In addition, we found evidence of gut inflammation in T1DM as reflected in the expression profiles of the Enterobacteria as well as in the human gut proteome. Finally, we were able to follow transmitted strain‑variants from the oral cavity to the gut at the individual omic levels, highlighting not only the transfer, but also the activity of the transmitted taxa along the gastrointestinal tract. Conclusions: Alterations of the oral microbiome in the context of T1DM impact the microbial communities in the lower gut, in particular through the reduction of “mouth‑to‑gut” transfer of Streptococcus salivarius. Our results indicate that the observed oral‑cavity‑driven gut microbiome changes may contribute towards the inflammatory processes involved in T1DM. Through the integration of multi‑omic analyses, we resolve strain‑variant “mouth‑to‑gut” transfer in a disease context. [less ▲]

binny: an automated binning algorithm to recover high-quality genomes from complex metagenomic datasets

in Briefings in Bioinformatics (2022)

The reconstruction of genomes is a critical step in genome-resolved metagenomics and for multi-omic data integration from microbial communities. Here, we present binny, a binning tool that produces high ... [more ▼]

The reconstruction of genomes is a critical step in genome-resolved metagenomics and for multi-omic data integration from microbial communities. Here, we present binny, a binning tool that produces high-quality metagenome-assembled genomes (MAG) from both contiguous and highly fragmented genomes. Based on established metrics, binny outperforms or is highly competitive with commonly used and state-of-the-art binning methods and finds unique genomes that could not be detected by other methods. binny uses k-mer-composition and coverage by metagenomic reads for iterative, nonlinear dimension reduction of genomic signatures as well as subsequent automated contig clustering with cluster assessment using lineage-specific marker gene sets. When compared with seven widely used binning algorithms, binny provides substantial amounts of uniquely identified MAGs and almost always recovers the most near-complete (⁠>95% pure, >90% complete) and high-quality (⁠>90% pure, >70% complete) genomes from simulated datasets from the Critical Assessment of Metagenome Interpretation initiative, as well as substantially more high-quality draft genomes, as defined by the Minimum Information about a Metagenome-Assembled Genome standard, from a real-world benchmark comprised of metagenomes from various environments than any other tested method. [less ▲]

Genotype-phenotype correlations in SCN8A-related disorders reveal prognostic and therapeutic implications

in Brain (2022)

We report detailed functional analyses and genotype-phenotype correlations in 392 individuals carrying disease-causing variants in SCN8A, encoding the voltage-gated Na+ channel NaV1.6, with the aim of ... [more ▼]

We report detailed functional analyses and genotype-phenotype correlations in 392 individuals carrying disease-causing variants in SCN8A, encoding the voltage-gated Na+ channel NaV1.6, with the aim of describing clinical phenotypes related to functional effects. Six different clinical subgroups could be identified: 1) Benign familial infantile epilepsy (BFIE) (n = 15, normal cognition, treatable seizures), 2) intermediate epilepsy (n = 33, mild ID, partially pharmaco-responsive), 3) developmental and epileptic encephalopathy (DEE, n = 177, severe ID, majority pharmaco-resistant), 4) generalized epilepsy (n = 20, mild to moderate ID, frequently with absence seizures), 5) unclassifiable epilepsy (n = 127), and 6) neurodevelopmental disorder without epilepsy (n = 20, mild to moderate ID). Groups 1–3 presented with focal or multifocal seizures (median age of onset: four months) and focal epileptiform discharges, whereas the onset of seizures in group 4 was later (median: 42 months) with generalized epileptiform discharges. We performed functional studies expressing missense variants in ND7/23 neuroblastoma cells and primary neuronal cultures using recombinant tetrodotoxin-insensitive human NaV1.6 channels and whole-cell patch-clamping. Two variants causing DEE showed a strong gain-of-function (GOF, hyperpolarising shift of steady-state activation, strongly increased neuronal firing rate), and one variant causing BFIE or intermediate epilepsy showed a mild GOF (defective fast inactivation, less increased firing). In contrast, all three variants causing generalized epilepsy induced a loss-of-function (LOF, reduced current amplitudes, depolarising shift of steady-state activation, reduced neuronal firing). Including previous studies, functional effects were known for 170 individuals. All 136 individuals carrying a functionally tested GOF variant had either focal (97, groups 1–3), or unclassifiable epilepsy (39), whereas 34 with a LOF variant had either generalized (14), no (11) or unclassifiable (6) epilepsy; only three had DEE. Computational modeling in the GOF group revealed a significant correlation between the severity of the electrophysiological and clinical phenotypes. GOF variant carriers responded significantly better to sodium channel blockers (SCBs) than to other anti-seizure medications, and the same applied for all individuals of groups 1–3.In conclusion, our data reveal clear genotype-phenotype correlations between age at seizure onset, type of epilepsy and gain- or loss-of-function effects of SCN8A variants. Generalized epilepsy with absence seizures is the main epilepsy phenotype of LOF variant carriers and the extent of the electrophysiological dysfunction of the GOF variants is a main determinant of the severity of the clinical phenotype in focal epilepsies. Our pharmacological data indicate that SCBs present a treatment option in SCN8A-related focal epilepsy with onset in the first year of life. [less ▲]