Interplay between probe design and test performance: overlap between genomic regions of interest, capture regions and high quality reference calls influence performance of WES-based assays.

Whole exome sequencing (WES)-based assays undergo rigorous validation before being implemented in diagnostic laboratories. This validation process generates experimental evidence that allows laboratories to predict the performance of the intended assay. The NA12878 Genome in a Bottle (GIAB) HapMap reference sample is commonly used for validation in diagnostic laboratories. We investigated what data points should be taken into consideration when validating WES-based assays using the GIAB reference in a diagnostic setting. We delineate specific factors that require special consideration and identify OMIM genes associated with diseases that may 'bypass' validation. Four replicates of the NA12878 sample were sequenced at the CHEO Genetics Diagnostic Laboratory on a NextSeq 500; the data were analyzed using the bcbio_nexgen v1.1.2 pipeline. The hap.py validation engine, Real Time Genomics vcfeval tool, and high confidence (HC) variant calls in HC regions available for the GIAB sample were used to validate the obtained variant calls. The same validation process was then used to evaluate variant calls obtained for the same sample by two other clinical diagnostic laboratories. We showed that variant calls in NA12878 can be confidently measured only in the regions that intersect between the GIAB HC regions and the target regions of exome capture. Of the 4139 (as of October 2019) OMIM genes associated with a phenotype and having a known molecular basis of disease, 84 were fully outside of the GIAB HC regions and many of the remaining OMIM genes were only partially covered by the HC regions. A significant proportion of variants identified in the NA12878 sample outside of the HC regions have unknown (UNK) status due to the absence of HC reference alleles. Verification of such calls is possible either by an alternative truth set or by orthogonal testing. Similarly, many variants outside of exome capture regions, if not accounted for, will be deemed false negatives due to insufficient probe coverage. Our results demonstrate the importance of the intersection between genomic regions of interest, capture regions, and the high confidence regions. If not considered, false and ambiguous variant calls could have a negative impact on diagnostic accuracy of the intended WES-based diagnostic assay and increase the need for confirmatory testing. To enable laboratories to identify 'problematic' regions and optimize validation efforts, we have made our VCF and BED files available in UCSC Genome Browser: NA12878 WES Benchmark. Relevant genes and genome annotations are evolving, we implemented a general purpose algorithm to cross-reference OMIM genes with the genomic regions of interest that can be applied to capture genes/regions outside HC regions (see repository of data material section).

Leveraging the power of new molecular technologies in the clinical setting requires unprecedented awareness of limitations and drawbacks: experience of one diagnostic laboratory.

BACKGROUND: Advances in molecular technologies and in-silico variant prediction tools offer wide-ranging opportunities in diagnostic settings, yet they also present with significant limitations. OBJECTIVE: Here, we contextualise the limitations of next-generation sequencing (NGS), multiplex ligation-dependent probe amplification (MLPA) and in-silico prediction tools routinely used by diagnostic laboratories by reviewing specific experiences from our diagnostic laboratory. METHODS: We investigated discordant annotations and/or incorrect variant 'callings' in exons of 56 genes constituting our cardiomyopathy and connective tissue disorder NGS panels. Discordant variants and segmental duplications (SD) were queried using the National Center for Biotechnology Information (NCBI) Basic Local Alignment Search Tool and the University of California Santa Cruz genome browser, respectively, to identify regions of high homology. Discrepant variant analyses by in-silico models were re-evaluated using updated file entries. RESULTS: We observed a 5% error rate in MYH7 variant 'calling' using MLPA, which resulted from >90% homology of the MYH7 probe-binding site to MYH6. SDs were detected in TTN, PKP2 and MYLK. SDs in MYLK presented the highest risk (15.7%) of incorrect variant 'calling'. The inaccurate 'callings' and discrepant in-silico predictions were resolved following detailed investigation into the source of error. CONCLUSION: Recognising the limitations described here may help avoid incorrect diagnoses and leverage the power of new molecular technologies in diagnostic settings.

Genetic Variation in the Ontario Neurodegenerative Disease Research Initiative.

BACKGROUND/OBJECTIVE: Apolipoprotein E (APOE) E4 is the main genetic risk factor for Alzheimer's disease (AD). Due to the consistent association, there is interest as to whether E4 influences the risk of other neurodegenerative diseases. Further, there is a constant search for other genetic biomarkers contributing to these phenotypes, such as microtubule-associated protein tau (MAPT) haplotypes. Here, participants from the Ontario Neurodegenerative Disease Research Initiative were genotyped to investigate whether the APOE E4 allele or MAPT H1 haplotype are associated with five neurodegenerative diseases: (1) AD and mild cognitive impairment (MCI), (2) amyotrophic lateral sclerosis, (3) frontotemporal dementia (FTD), (4) Parkinson's disease, and (5) vascular cognitive impairment. METHODS: Genotypes were defined for their respective APOE allele and MAPT haplotype calls for each participant, and logistic regression analyses were performed to identify the associations with the presentations of neurodegenerative diseases. RESULTS: Our work confirmed the association of the E4 allele with a dose-dependent increased presentation of AD, and an association between the E4 allele alone and MCI; however, the other four diseases were not associated with E4. Further, the APOE E2 allele was associated with decreased presentation of both AD and MCI. No associations were identified between MAPT haplotype and the neurodegenerative disease cohorts; but following subtyping of the FTD cohort, the H1 haplotype was significantly associated with progressive supranuclear palsy. CONCLUSION: This is the first study to concurrently analyze the association of APOE isoforms and MAPT haplotypes with five neurodegenerative diseases using consistent enrollment criteria and broad phenotypic analysis.

Étude de variance génétique dans le cadre de l'initiative de recherche sur les maladies neurodégénératives en Ontario. Contexte/Objectif : L'apolipoprotéine E4 (ApoE4) constitue le principal facteur de risque génétique de la maladie d'Alzheimer. En raison de cette association systématique, il existe un intérêt certain à savoir dans quelle mesure cette classe d'apolipoprotéines peut influencer le risque d'autres maladies neurodégénératives. En outre, le milieu de la recherche n'a de cesse d'identifier d'autres biomarqueurs génétiques, par exemple les haplotypes H1 de la protéine tau associée aux microtubules, qui contribuent à certains phénotypes, Dans le cadre de cette étude, des participants à l'initiative de recherche sur les maladies neurodégénératives en Ontario ont été « génotypés ¼ afin de déterminer si l'ApoE4 ou l'haplotype H1 mentionné ci-dessus peuvent être associés à cinq maladies neurodégénératives : 1) la maladie d'Alzheimer et d'autres troubles cognitifs légers ; 2) la sclérose latérale amyotrophique ; 3) la démence fronto-temporale ; 4) la maladie de Parkinson ; 5) et finalement les déficits cognitifs d'origine vasculaire. Méthodes : Pour chaque participant, la cartographie des génotypes a été établie en fonction de leur ApoE4 respectif et de la présence d'haplotypes H1 de la protéine tau associée aux microtubules. Des analyses de régression logistique ont été ensuite effectuées dans le but d'identifier de possibles liens avec ces maladies neurodégénératives. Résultats : Nos travaux ont confirmé l'association entre l'ApoE4 et une plus grande occurrence de cas d'Alzheimer, et ce, en tenant compte de l'effet d'une dose de médicament. Ils ont aussi montré une association entre la seule ApoE4 et des troubles cognitifs légers. Cela dit, il convient de préciser que les quatre autres maladies n'ont pas été associées à cet allèle. Plus encore, nous avons trouvé que l'allèle E2 de l'apolipoprotéine était associé à une occurrence plus faible de cas d'Alzheimer et de troubles cognitifs légers. Fait à souligner, aucune association n'a été détectée entre l'haplotype H1 de la protéine tau associée aux microtubules et nos cohortes atteintes de maladies neurodégénératives. Toutefois, à la suite du sous-typage de la cohorte de participants atteints de démence fronto-temporale, il s'est avéré que l'haplotype H1 était associé de façon notable à la paralysie supra-nucléaire progressive. Conclusion : Il s'agit de la première étude à analyser simultanément, au moyen de critères de participation cohérents et d'une analyse phénotypique élargie, les associations entre les isoformes de l'ApoE, l'haplotype H1 de la protéine tau associée aux microtubules et cinq maladies neurodégénératives.

Reinterpretation of sequence variants: one diagnostic laboratory's experience, and the need for standard guidelines.

PurposeThe advent of next-generation sequencing resulted in substantial increases in the number of variants detected, interpreted, and reported by molecular genetics diagnostic laboratories. Recent publications have provided standards for the interpretation of sequence variants, but there are currently no standards regarding reinterpretation of these variants. Recognizing that significant changes in variant classification may occur over time, many genetics diagnostic laboratories have independently developed practices for variant reinterpretation. The purpose of this study is to describe our laboratory approach to variant reinterpretation.MethodsWe surveyed eight genetics diagnostic laboratories in Canada and the United States.ResultsEach laboratory had differing protocols, but most felt that clinically relevant changes to variant classifications should be communicated to ordering providers. Based on results of this survey and our experience, we developed a cost-effective and resource-efficient approach to variant reinterpretation.ConclusionOngoing variant reinterpretation is required to maintain the highest standards for delivering genetics laboratory services. Our approach to variant reinterpretation offers an efficient solution that does not compromise accuracy or timely delivery of genetics laboratory services.

MTHFSD and DDX58 are novel RNA-binding proteins abnormally regulated in amyotrophic lateral sclerosis.

Tar DNA-binding protein 43 (TDP-43) is an RNA-binding protein normally localized to the nucleus of cells, where it elicits functions related to RNA metabolism such as transcriptional regulation and alternative splicing. In amyotrophic lateral sclerosis, TDP-43 is mislocalized from the nucleus to the cytoplasm of diseased motor neurons, forming ubiquitinated inclusions. Although mutations in the gene encoding TDP-43, TARDBP, are found in amyotrophic lateral sclerosis, these are rare. However, TDP-43 pathology is common to over 95% of amyotrophic lateral sclerosis cases, suggesting that abnormalities of TDP-43 play an active role in disease pathogenesis. It is our hypothesis that a loss of TDP-43 from the nucleus of affected motor neurons in amyotrophic lateral sclerosis will lead to changes in RNA processing and expression. Identifying these changes could uncover molecular pathways that underpin motor neuron degeneration. Here we have used translating ribosome affinity purification coupled with microarray analysis to identify the mRNAs being actively translated in motor neurons of mutant TDP-43(A315T) mice compared to age-matched non-transgenic littermates. No significant changes were found at 5 months (presymptomatic) of age, but at 10 months (symptomatic) the translational profile revealed significant changes in genes involved in RNA metabolic process, immune response and cell cycle regulation. Of 28 differentially expressed genes, seven had a ≥ 2-fold change; four were validated by immunofluorescence labelling of motor neurons in TDP-43(A315T) mice, and two of these were confirmed by immunohistochemistry in amyotrophic lateral sclerosis cases. Both of these identified genes, DDX58 and MTHFSD, are RNA-binding proteins, and we show that TDP-43 binds to their respective mRNAs and we identify MTHFSD as a novel component of stress granules. This discovery-based approach has for the first time revealed translational changes in motor neurons of a TDP-43 mouse model, identifying DDX58 and MTHFSD as two TDP-43 targets that are misregulated in amyotrophic lateral sclerosis.

Characterizing familial corticobasal syndrome due to Alzheimer's disease pathology and PSEN1 mutations.

INTRODUCTION: Corticobasal syndrome (CBS) resulting from genetic Alzheimer's disease (AD) has been described only once. Whether familial CBS-AD is a distinct clinical entity with its own imaging signature remains unknown. METHODS: Four individuals with CBS from two families underwent detailed assessment. For two individuals, regional atrophy and hypoperfusion were compared to autopsy-confirmed typical late-onset AD and corticobasal degeneration, as well as genetically proven PSEN1 cases with an amnestic presentation. RESULTS: One family harbored a novel mutation in PSEN1:p.Phe283Leu. MRI demonstrated severe parietal, perirolandic, and temporal atrophy, with relative sparing of frontal and ipsilateral hippocampal regions. Autopsy confirmed pure AD pathology. The other family harbored a known PSEN1 mutation:p.Gly378Val. DISCUSSION: This report confirms familial CBS-AD as a distinct clinical entity, with a parietal-perirolandic-temporal atrophy signature. It illustrates the clinical heterogeneity that can occur despite a shared genetic cause and underscores the need for biomarkers such as amyloid imaging during life.

Hypermethylation of the CpG island near the G4C2 repeat in ALS with a C9orf72 expansion.

The G4C2 repeat expansion in C9orf72 is the most common known cause of amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD). We tested the hypothesis that the repeat expansion causes aberrant CpG methylation near the G4C2 repeat, which could be responsible for the downregulation of gene expression. We investigated the CpG methylation profile by two methods using genomic DNA from the blood of individuals with ALS (37 expansion carriers and 64 noncarriers), normal controls (n = 76), and family members of 7 ALS probands with the expansion. We report that hypermethylation of the CpG island 5' of the G4C2 repeat is associated with the presence of the expansion (p < 0.0001). A higher degree of methylation was significantly correlated with a shorter disease duration (p < 0.01), associated with familial ALS (p = 0.009) and segregated with the expansion in 7 investigated families. Notably, we did not detect methylation for either normal or intermediate alleles (up to 43 repeats), bringing to question the current cutoff of 30 repeats for pathological alleles. Our study raises several important questions for the future investigation of large data sets, such as whether the degree of methylation corresponds to clinical presentation (ALS versus FTLD).

Coding mutations in SORL1 and Alzheimer disease.

OBJECTIVE: Common single nucleotide polymorphisms in the SORL1 gene have been associated with late onset Alzheimer disease (LOAD), but causal variants have not been fully characterized nor has the mechanism been established. The study was undertaken to identify functional SORL1 mutations in patients with LOAD. METHODS: This was a family- and cohort-based genetic association study. Caribbean Hispanics with familial and sporadic LOAD and similarly aged controls were recruited from the United States and the Dominican Republic, and patients with sporadic disease of Northern European origin were recruited from Canada. Prioritized coding variants in SORL1 were detected by targeted resequencing and validated by genotyping in additional family members and unrelated healthy controls. Variants transfected into human embryonic kidney 293 cell lines were tested for Aß40 and Aß42 secretion, and the amount of the amyloid precursor protein (APP) secreted at the cell surface was determined. RESULTS: Seventeen coding exonic variants were significantly associated with disease. Two rare variants (rs117260922-E270K and rs143571823-T947M) with minor allele frequency (MAF) < 1% and 1 common variant (rs2298813-A528T) with MAF = 14.9% segregated within families and were deemed deleterious to the coding protein. Transfected cell lines showed increased Aß40 and Aß42 secretion for the rare variants (E270K and T947M) and increased Aß42 secretion for the common variant (A528T). All mutants increased the amount of APP at the cell surface, although in slightly different ways, thereby failing to direct full-length APP into the retromer-recycling endosome pathway. INTERPRETATION: Common and rare variants in SORL1 elevate the risk of LOAD by directly affecting APP processing, which in turn can result in increased Aß40 and Aß42 secretion.

Rare coding mutations identified by sequencing of Alzheimer disease genome-wide association studies loci.

OBJECTIVE: To detect rare coding variants underlying loci detected by genome-wide association studies (GWAS) of late onset Alzheimer disease (LOAD). METHODS: We conducted targeted sequencing of ABCA7, BIN1, CD2AP, CLU, CR1, EPHA1, MS4A4A/MS4A6A, and PICALM in 3 independent LOAD cohorts: 176 patients from 124 Caribbean Hispanics families, 120 patients and 33 unaffected individuals from the 129 National Institute on Aging LOAD Family Study; and 263 unrelated Canadian individuals of European ancestry (210 sporadic patients and 53 controls). Rare coding variants found in at least 2 data sets were genotyped in independent groups of ancestry-matched controls. Additionally, the Exome Aggregation Consortium was used as a reference data set for population-based allele frequencies. RESULTS: Overall we detected a statistically significant 3.1-fold enrichment of the nonsynonymous mutations in the Caucasian LOAD cases compared with controls (p = 0.002) and no difference in synonymous variants. A stop-gain mutation in ABCA7 (E1679X) and missense mutation in CD2AP (K633R) were highly significant in Caucasian LOAD cases, and mutations in EPHA1 (P460L) and BIN1 (K358R) were significant in Caribbean Hispanic families with LOAD. The EPHA1 variant segregated completely in an extended Caribbean Hispanic family and was also nominally significant in the Caucasians. Additionally, BIN1 (K358R) segregated in 2 of the 6 Caribbean Hispanic families where the mutations were discovered. INTERPRETATION: Targeted sequencing of confirmed GWAS loci revealed an excess burden of deleterious coding mutations in LOAD, with the greatest burden observed in ABCA7 and BIN1. Identifying coding variants in LOAD will facilitate the creation of tractable models for investigation of disease-related mechanisms and potential therapies.

Inbreeding among Caribbean Hispanics from the Dominican Republic and its effects on risk of Alzheimer disease.

BACKGROUND: Inbreeding can be associated with a modification of disease risk due to excess homozygosity of recessive alleles affecting a wide range of phenotypes. We estimated the inbreeding coefficient in Caribbean Hispanics and examined its effects on risk of late-onset Alzheimer disease. METHODS: The inbreeding coefficient was calculated in 3,392 subjects (1,451 late-onset Alzheimer disease patients and 1,941 age-matched healthy controls) of Caribbean Hispanic ancestry using 177,997 nearly independent single-nucleotide polymorphisms from genome-wide array. The inbreeding coefficient was estimated using the excess homozygosity method with and without adjusting for admixture. RESULTS: The average inbreeding coefficient in Caribbean Hispanics without accounting for admixture was F = 0.018 (±0.048), suggesting a mating equivalent to that of second cousins or second cousins once removed. Adjusting for admixture from three parent populations, the average inbreeding coefficient was found to be 0.0034 (±0.019) or close to third-cousin mating. Inbreeding coefficient was a significant predictor of Alzheimer disease when age, sex, and APOE genotype were used as adjusting covariates (P = 0.03). CONCLUSION: The average inbreeding coefficient of this population is significantly higher than that of the general Caucasian populations in North America. The high rate of inbreeding resulting in increased frequency of recessive variants is advantageous for the identification of rare variants associated with late-onset Alzheimer disease.Genet Med 17 8, 639-643.

Association of Dual-Task Gait Cost and White Matter Hyperintensity Burden Poststroke: Results From the ONDRI.

BACKGROUND: Acute change in gait speed while performing a mental task [dual-task gait cost (DTC)], and hyperintensity magnetic resonance imaging signals in white matter are both important disability predictors in older individuals with history of stroke (poststroke). It is still unclear, however, whether DTC is associated with overall hyperintensity volume from specific major brain regions in poststroke. METHODS: This is a cohort study with a total of 123 older (69 ± 7 years of age) participants with history of stroke were included from the Ontario Neurodegenerative Disease Research Initiative. Participants were clinically assessed and had gait performance assessed under single- and dual-task conditions. Structural neuroimaging data were analyzed to measure both, white matter hyperintensity (WMH) and normal appearing volumes. Percentage of WMH volume in frontal, parietal, occipital, and temporal lobes as well as subcortical hyperintensities in basal ganglia + thalamus were the main outcomes. Multivariate models investigated associations between DTC and hyperintensity volumes, adjusted for age, sex, years of education, global cognition, vascular risk factors, APOE4 genotype, residual sensorimotor symptoms from previous stroke and brain volume. RESULTS: There was a significant positive global linear association between DTC and hyperintensity burden (adjusted Wilks' λ = .87, P = .01). Amongst all WMH volumes, hyperintensity burden from basal ganglia + thalamus provided the most significant contribution to the global association (adjusted ß = .008, η2 = .03; P = .04), independently of brain atrophy. CONCLUSIONS: In poststroke, increased DTC may be an indicator of larger white matter damages, specifically in subcortical regions, which can potentially affect the overall cognitive processing and decrease gait automaticity by increasing the cortical control over patients' locomotion.

Axial Impairment Following Deep Brain Stimulation in Parkinson's Disease: A Surgicogenomic Approach.

BACKGROUND: Postoperative outcome following deep brain stimulation (DBS) of the subthalamic nucleus is variable, particularly with respect to axial motor improvement. We hypothesized a genetic underpinning to the response to surgical intervention, termed "surgicogenomics". OBJECTIVE: We aimed to identify genetic variants associated with clinical heterogeneity in DBS outcome of Parkinson's disease (PD) patients that could then be applied clinically to target selection leading to improved surgical outcome. METHODS: Retrospective clinical data was extracted from 150 patient's charts. Each individual was genotyped using the genome-wide NeuroX array tailored to study neurologic diseases. Genetic data were clustered based on surgical outcome assessed by comparing pre- and post-operative scores of levodopa equivalent daily dose and axial impairment at one and five years post-surgery. Allele frequencies were compared between patients with excellent vs. moderate/poor outcomes grouped using a priori defined cut-offs. We analyzed common variants, burden of rare coding variants, and PD polygenic risk score. RESULTS: NeuroX identified 2,917 polymorphic markers at 113 genes mapped to known PD loci. The gene-burden analyses of 202 rare nonsynonymous variants suggested a nominal association of axial impairment with 14 genes (most consistent with CRHR1, IP6K2, and PRSS3). The strongest association with surgical outcome was detected between a reduction in levodopa equivalent daily dose and common variations tagging two linkage disequilibrium blocks with SH3GL2. CONCLUSION: Once validated in independent populations, our findings may be implemented to improve patient selection for DBS in PD.

The Intersection between COVID-19, the Gene Family of ACE2 and Alzheimer's Disease.

We reviewed factors that might influence COVID-19 outcomes (eg, neurological symptoms), including the link to Alzheimer's disease. Since the virus triggers COVID-19 infection through binding to ACE2, we focused on the ACE2 gene family, including ACE. Both ACE2 and ACE are involved in the renin-angiotensin system (RAS). In general, ACE causes inflammation and vasoconstriction, while ACE2 leads to anti-inflammation activity and vasodilation. The disturbed balance between these counter-regulatory pathways could influence susceptibility to COVID-19. Notably, dysregulation of the RAS-equilibrium contributes to Alzheimer's disease. Differences in the incidence and symptoms of COVID-19 in diverse populations could be attributed to variability in the human genome. For example, ACE and ACE2 variations could modify the outcome of COVID-19 in different populations. It would be important to conduct genome-wide studies to detect variants influencing COVID-19 presentation, with a special focus on variants affecting immune-related pathways and expression of RAS-related genes.

ALU transposition induces familial hypertrophic cardiomyopathy.

BACKGROUND: Hypertrophic cardiomyopathy (HCM) is characterized by left ventricular hypertrophy (LVH) in the absence of predisposing cardiovascular conditions. Pathogenic variants in at least 16 cardiac sarcomeric genes have been implicated in HCM, most of which act in a dominant-negative fashion. However loss-of-function (haploinsufficiency) is the most common disease mechanism for pathogenic variants in MYBPC3, suggesting that MYBPC3 complete deletion may play a role in HCM pathogenesis. Here, we investigate MYBPC3 complete deletion as a disease mechanism in HCM by analyzing two unrelated patients with confirmed diagnosis of HCM that tested negative by Sanger sequencing analysis. METHODS: MYBPC3 complete deletion was investigated by Multiplex ligation-dependent probe amplification (MLPA) and microarray analyses. The mechanism of deletion was investigated by interrogating the SINEBase database. RESULTS: Patient-1 was diagnosed with nonobstructive HCM in his mid-40s while undergoing assessment for palpitations, and patient-2 with obstructive HCM in his late-20s while undergoing systolic heart murmur assessment for an unrelated illness. MLPA testing revealed a heterozygous deletion of all MYBPC3 exons in both patients. Subsequent microarray testing confirmed these deletions which extended beyond the 5' and 3' ends of MYBPC3. Genomic assessment suggested that these deletions resulted from Alu/Alu-homologous recombination. CONCLUSION: Our results demonstrate that haploinsufficiency resulting from MYBPC3 complete deletion, potentially mediated by Alu recombination, is an important disease mechanism in cardiomyopathy and emphasizes the importance of copy number variation analysis in patients clinically suspected of HCM.

Adopting High-Resolution Allele Frequencies Substantially Expedites Variant Interpretation in Genetic Diagnostic Laboratories.

A cohort of 1242 individuals tested in a clinical diagnostic laboratory was used to test whether the filtering allele frequencies (FAFs)-based framework, recently recommended for MHY7-associated cardiomyopathy, is extendable to 45 cardiomyopathy genes. Statistical analysis revealed a threshold of 0.00164% for the extreme outlier allele frequencies (AFs), based on the Genome Aggregation Database (exome fraction) total AFs of 138 unique pathogenic and likely pathogenic variants; 135 of them (97.8%) had AFs of <0.004%, the recommended threshold to apply moderate pathogenicity evidence for MYH7-associated cardiomyopathy. Of the 460 cases reported with only variant(s) of unknown clinical significance (VUCSs), 97 (21%) solely had VUCSs with FAFs >0.03%, frequencies above which were estimated herein as strong evidence against pathogenicity. Interestingly, 74.5% (172/231) of the unique VUCSs with FAFs >0.03% had Genome Aggregation Database maximum allele frequencies across all populations AFs >0.1%, deemed herein as stand-alone evidence against pathogenicity. Accordingly, using an FAF threshold of >0.1%, compared with AF >1%, led us to issue considerably more (25.9% versus 41.3%) negative patient reports. Also, 82.7% (N = 629) of the unique classified benign or likely benign variants with AFs <1% had FAFs >0.1%, reinforcing the use of this filtering strategy. Together, these data demonstrate that implementing FAF thresholds may considerably decrease the amount of variant interpretations and significantly reduce the cost of genetic testing for clinical genetic laboratories, without compromising the accuracy of genetic diagnostic services.

Genetic Diagnostic Testing for Inherited Cardiomyopathies: Considerations for Offering Multi-Gene Tests in a Health Care Setting.

Inherited cardiomyopathies (ICs) are a major cause of heart disease. Given their marked clinical and genetic heterogeneity, the content and clinical utility of IC multi-gene panels has been the topic of continuous debate. Our genetics diagnostic laboratory has been providing clinical diagnostic testing for ICs since 2012. We began by testing nine genes and expanded our panel by fivefold in 2015. Here, we describe the implementation of a cost-effective next-generation sequencing (NGS)-based assay for testing of IC genes, including a protocol that minimizes the amount of Sanger sequencing required to confirm variants identified by NGS, which reduces the cost and time of testing. The NGS assay was developed for the simultaneous analysis of 45 IC genes and was assessed for the impact of panel expansion on variant detection, turnaround time, and cost of testing in a cohort of 993 patients. The assay led to a considerable reduction in test cost and turnaround time. However, only a marginal increase was observed in the diagnostic yield, whereas the rate of inconclusive findings increased considerably. These findings suggest that the ongoing evaluation of gene content and monitoring of clinical utility for multi-gene tests are essential to achieve maximum clinical utility of multi-gene tests in a publicly funded health care setting.

An APOE-independent cis-eSNP on chromosome 19q13.32 influences tau levels and late-onset Alzheimer's disease risk.

Although multiple susceptibility loci for late-onset Alzheimer's disease (LOAD) have been identified, a large portion of the genetic risk for this disease remains unexplained. LOAD risk may be associated with single-nucleotide polymorphisms responsible for changes in gene expression (eSNPs). To detect eSNPs associated with LOAD, we integrated data from LOAD genome-wide association studies and expression quantitative trait loci using Sherlock (a Bayesian statistical method). We identified a cis-regulatory eSNP (rs2927438) located on chromosome 19q13.32, for which subsequent analyses confirmed the association with both LOAD risk and the expression level of several nearby genes. Importantly, rs2927438 may represent an APOE-independent LOAD eSNP according to the weak linkage disequilibrium of rs2927438 with the 2 polymorphisms (rs7412 and rs429358) defining the APOE-ε2, -ε3, and -ε4 alleles. Furthermore, rs2927438 does not influence chromatin interaction events at the APOE locus or cis-regulation of APOE expression. Further exploratory analysis revealed that rs2927438 is significantly associated with tau levels in the cerebrospinal fluid. Our findings suggest that rs2927438 may confer APOE-independent risk for LOAD.

Targeted Next-generation Sequencing and Bioinformatics Pipeline to Evaluate Genetic Determinants of Constitutional Disease.

Next-generation sequencing (NGS) is quickly revolutionizing how research into the genetic determinants of constitutional disease is performed. The technique is highly efficient with millions of sequencing reads being produced in a short time span and at relatively low cost. Specifically, targeted NGS is able to focus investigations to genomic regions of particular interest based on the disease of study. Not only does this further reduce costs and increase the speed of the process, but it lessens the computational burden that often accompanies NGS. Although targeted NGS is restricted to certain regions of the genome, preventing identification of potential novel loci of interest, it can be an excellent technique when faced with a phenotypically and genetically heterogeneous disease, for which there are previously known genetic associations. Because of the complex nature of the sequencing technique, it is important to closely adhere to protocols and methodologies in order to achieve sequencing reads of high coverage and quality. Further, once sequencing reads are obtained, a sophisticated bioinformatics workflow is utilized to accurately map reads to a reference genome, to call variants, and to ensure the variants pass quality metrics. Variants must also be annotated and curated based on their clinical significance, which can be standardized by applying the American College of Medical Genetics and Genomics Pathogenicity Guidelines. The methods presented herein will display the steps involved in generating and analyzing NGS data from a targeted sequencing panel, using the ONDRISeq neurodegenerative disease panel as a model, to identify variants that may be of clinical significance.