Genome Sequencing for Cases Unsolved by Exome Sequencing: Identifying a Single-Exon Deletion in TBCK in a Case from 30 Years Ago.

In patients with neurodevelopmental disorders (NDDs), exome sequencing (ES), the diagnostic gold standard, reveals an underlying monogenic condition in only approximately 40% of cases. We report the case of a female patient with profound NDD who died 30 years ago at the age of 3 years and for whom genome sequencing (GS) now identified a single-exon deletion in TBCK previously missed by ExomeDepth, the copy number variation (CNV) detection algorithm in ES.Deoxyribonucleic acid (DNA) was extracted from frozen muscle tissue of the index patient and the parents' blood. Genome data were analyzed for structural variants and single nucleotide variants (SUVs)/indels as part of the Bavarian Genomes consortium project.Biallelic variants in TBCK, which are linked to the autosomal recessive disorder TBCK syndrome, were detected in the affected individual: a novel frameshift variant and a deletion of exon 23, previously established as common but underrecognized pathogenic variant in individuals with TBCK syndrome. While in the foregoing ES analysis, calling algorithms for (SNVs)/indels were able to identify the frameshift variant, ExomeDepth failed to call the intragenic deletion.Our case illustrates the added value of GS for the detection of single-exon deletions for which calling from ES data remains challenging and confirms that the deletion of exon 23 in TBCK may be underdiagnosed in patients with NDDs. Furthermore, it shows the importance of "molecular or genetic autopsy" allowing genetic risk counseling for family members as well as the end of a diagnostic odyssey of 30 years.

Identification of Restless Legs Syndrome Genes by Mutational Load Analysis.

OBJECTIVE: Restless legs syndrome is a frequent neurological disorder with substantial burden on individual well-being and public health. Genetic risk loci have been identified, but the causatives genes at these loci are largely unknown, so that functional investigation and clinical translation of molecular research data are still inhibited. To identify putatively causative genes, we searched for highly significant mutational burden in candidate genes. METHODS: We analyzed 84 candidate genes in 4,649 patients and 4,982 controls by next generation sequencing using molecular inversion probes that targeted mainly coding regions. The burden of low-frequency and rare variants was assessed, and in addition, an algorithm (binomial performance deviation analysis) was established to estimate independently the sequence variation in the probe binding regions from the variation in sequencing depth. RESULTS: Highly significant results (considering the number of genes in the genome) of the conventional burden test and the binomial performance deviation analysis overlapped significantly. Fourteen genes were highly significant by one method and confirmed with Bonferroni-corrected significance by the other to show a differential burden of low-frequency and rare variants in restless legs syndrome. Nine of them (AAGAB, ATP2C1, CNTN4, COL6A6, CRBN, GLO1, NTNG1, STEAP4, VAV3) resided in the vicinity of known restless legs syndrome loci, whereas 5 (BBS7, CADM1, CREB5, NRG3, SUN1) have not previously been associated with restless legs syndrome. Burden test and binomial performance deviation analysis also converged significantly in fine-mapping potentially causative domains within these genes. INTERPRETATION: Differential burden with intragenic low-frequency variants reveals putatively causative genes in restless legs syndrome. ANN NEUROL 2020;87:184-193.

Targeted resequencing and systematic in vivo functional testing identifies rare variants in MEIS1 as significant contributors to restless legs syndrome.

Restless legs syndrome (RLS) is a common neurologic condition characterized by nocturnal dysesthesias and an urge to move, affecting the legs. RLS is a complex trait, for which genome-wide association studies (GWASs) have identified common susceptibility alleles of modest (OR 1.2-1.7) risk at six genomic loci. Among these, variants in MEIS1 have emerged as the largest risk factors for RLS, suggesting that perturbations in this transcription factor might be causally related to RLS susceptibility. To establish this causality, direction of effect, and total genetic burden of MEIS1, we interrogated 188 case subjects and 182 control subjects for rare alleles not captured by previous GWASs, followed by genotyping of ∼3,000 case subjects and 3,000 control subjects, and concluded with systematic functionalization of all discovered variants using a previously established in vivo model of neurogenesis. We observed a significant excess of rare MEIS1 variants in individuals with RLS. Subsequent assessment of all nonsynonymous variants by in vivo complementation revealed an excess of loss-of-function alleles in individuals with RLS. Strikingly, these alleles compromised the function of the canonical MEIS1 splice isoform but were irrelevant to an isoform known to utilize an alternative 3' sequence. Our data link MEIS1 loss of function to the etiopathology of RLS, highlight how combined sequencing and systematic functional annotation of rare variation at GWAS loci can detect risk burden, and offer a plausible explanation for the specificity of phenotypic expressivity of loss-of-function alleles at a locus broadly necessary for neurogenesis and neurodevelopment.

Restless legs syndrome-associated intronic common variant in Meis1 alters enhancer function in the developing telencephalon.

Genome-wide association studies (GWAS) identified the MEIS1 locus for Restless Legs Syndrome (RLS), but causal single nucleotide polymorphisms (SNPs) and their functional relevance remain unknown. This locus contains a large number of highly conserved noncoding regions (HCNRs) potentially functioning as cis-regulatory modules. We analyzed these HCNRs for allele-dependent enhancer activity in zebrafish and mice and found that the risk allele of the lead SNP rs12469063 reduces enhancer activity in the Meis1 expression domain of the murine embryonic ganglionic eminences (GE). CREB1 binds this enhancer and rs12469063 affects its binding in vitro. In addition, MEIS1 target genes suggest a role in the specification of neuronal progenitors in the GE, and heterozygous Meis1-deficient mice exhibit hyperactivity, resembling the RLS phenotype. Thus, in vivo and in vitro analysis of a common SNP with small effect size showed allele-dependent function in the prospective basal ganglia representing the first neurodevelopmental region implicated in RLS.

Phenotypic and genome-wide studies on dicarbonyls: major associations to glomerular filtration rate and gamma-glutamyltransferase activity.

BACKGROUND: The dicarbonyl compounds methylglyoxal (MG), glyoxal (GO) and 3-deoxyglucosone (3-DG) have been linked to various diseases. However, disease-independent phenotypic and genotypic association studies with phenome-wide and genome-wide reach, respectively, have not been provided. METHODS: MG, GO and 3-DG were measured by LC-MS in 1304 serum samples of two populations (KORA, n = 482; BiDirect, n = 822) and assessed for associations with genome-wide SNPs (GWAS) and with phenome-wide traits. Redundancy analysis (RDA) was used to identify major independent trait associations. FINDINGS: Mutual correlations of dicarbonyls were highly significant, being stronger between MG and GO (ρ = 0.6) than between 3-DG and MG or GO (ρ = 0.4). Significant phenotypic results included associations of all dicarbonyls with sex, waist-to-hip ratio, glomerular filtration rate (GFR), gamma-glutamyltransferase (GGT), and hypertension, of MG and GO with age and C-reactive protein, of GO and 3-DG with glucose and antidiabetics, of MG with contraceptives, of GO with ferritin, and of 3-DG with smoking. RDA revealed GFR, GGT and, in case of 3-DG, glucose as major contributors to dicarbonyl variance. GWAS did not identify genome-wide significant loci. SNPs previously associated with glyoxalase activity did not reach nominal significance. When multiple testing was restricted to the lead SNPs of GWASs on the traits selected by RDA, 3-DG was found to be associated (p = 2.3 × 10-5) with rs1741177, an eQTL of NF-κB inhibitor NFKBIA. INTERPRETATION: This large-scale, population-based study has identified numerous associations, with GFR and GGT being of pivotal importance, providing unbiased perspectives on dicarbonyls beyond the current state. FUNDING: Deutsche Forschungsgemeinschaft, Helmholtz Munich, German Centre for Cardiovascular Research (DZHK), German Federal Ministry of Research and Education (BMBF).

Episignature analysis of moderate effects and mosaics.

DNA methylation classifiers ("episignatures") help to determine the pathogenicity of variants of uncertain significance (VUS). However, their sensitivity is limited due to their training on unambiguous cases with strong-effect variants so that the classification of variants with reduced effect size or in mosaic state may fail. Moreover, episignature evaluation of mosaics as a function of their degree of mosaicism has not been developed so far. We improved episignatures with respect to three categories. Applying (i) minimum-redundancy-maximum-relevance feature selection we reduced their length by up to one order of magnitude without loss of accuracy. Performing (ii) repeated re-training of a support vector machine classifier by step-wise inclusion of cases in the training set that reached probability scores larger than 0.5, we increased the sensitivity of the episignature-classifiers by 30%. In the newly diagnosed patients we confirmed the association between DNA methylation aberration and age at onset of KMT2B-deficient dystonia. Moreover, we found evidence for allelic series, including KMT2B-variants with moderate effects and comparatively mild phenotypes such as late-onset focal dystonia. Retrained classifiers also can detect mosaics that previously remained below the 0.5-threshold, as we showed for KMT2D-associated Kabuki syndrome. Conversely, episignature-classifiers are able to revoke erroneous exome calls of mosaicism, as we demonstrated by (iii) comparing presumed mosaic cases with a distribution of artificial in silico-mosaics that represented all the possible variation in degree of mosaicism, variant read sampling and methylation analysis.

ExomeChip-based rare variant association study in restless legs syndrome.

Restless legs syndrome (RLS) is a common sleep-related movement disorder in populations of European descent and disease risk is strongly influenced by genetic factors. Common variants have been assessed extensively in several genome-wide association studies, but the contribution of rarer genetic variation has not been investigated at this scale. We therefore genotyped a case-control set of 9246 individuals for mainly rare and low frequency exonic variants using the Illumina ExomeChip. However, standard single variant and gene-level association tests were negative. This does not preclude a role of rare variants in RLS, but is likely due to the small sample size and the limited selection of rare genetic variation captured on the array. Therefore, exome or whole genome sequencing should be performed rather than increasing the sample size of ExomeChip studies in order to identify rare risk variants for RLS.

Identification of novel risk loci for restless legs syndrome in genome-wide association studies in individuals of European ancestry: a meta-analysis.

BACKGROUND: Restless legs syndrome is a prevalent chronic neurological disorder with potentially severe mental and physical health consequences. Clearer understanding of the underlying pathophysiology is needed to improve treatment options. We did a meta-analysis of genome-wide association studies (GWASs) to identify potential molecular targets. METHODS: In the discovery stage, we combined three GWAS datasets (EU-RLS GENE, INTERVAL, and 23andMe) with diagnosis data collected from 2003 to 2017, in face-to-face interviews or via questionnaires, and involving 15 126 cases and 95 725 controls of European ancestry. We identified common variants by fixed-effect inverse-variance meta-analysis. Significant genome-wide signals (p≤5 × 10-8) were tested for replication in an independent GWAS of 30 770 cases and 286 913 controls, followed by a joint analysis of the discovery and replication stages. We did gene annotation, pathway, and gene-set-enrichment analyses and studied the genetic correlations between restless legs syndrome and traits of interest. FINDINGS: We identified and replicated 13 new risk loci for restless legs syndrome and confirmed the previously identified six risk loci. MEIS1 was confirmed as the strongest genetic risk factor for restless legs syndrome (odds ratio 1·92, 95% CI 1·85-1·99). Gene prioritisation, enrichment, and genetic correlation analyses showed that identified pathways were related to neurodevelopment and highlighted genes linked to axon guidance (associated with SEMA6D), synapse formation (NTNG1), and neuronal specification (HOXB cluster family and MYT1). INTERPRETATION: Identification of new candidate genes and associated pathways will inform future functional research. Advances in understanding of the molecular mechanisms that underlie restless legs syndrome could lead to new treatment options. We focused on common variants; thus, additional studies are needed to dissect the roles of rare and structural variations. FUNDING: Deutsche Forschungsgemeinschaft, Helmholtz Zentrum München-Deutsches Forschungszentrum für Gesundheit und Umwelt, National Research Institutions, NHS Blood and Transplant, National Institute for Health Research, British Heart Foundation, European Commission, European Research Council, National Institutes of Health, National Institute of Neurological Disorders and Stroke, NIH Research Cambridge Biomedical Research Centre, and UK Medical Research Council.

Connecting Anxiety and Genomic Copy Number Variation: A Genome-Wide Analysis in CD-1 Mice.

Genomic copy number variants (CNVs) have been implicated in multiple psychiatric disorders, but not much is known about their influence on anxiety disorders specifically. Using next-generation sequencing (NGS) and two additional array-based genotyping approaches, we detected CNVs in a mouse model consisting of two inbred mouse lines showing high (HAB) and low (LAB) anxiety-related behavior, respectively. An influence of CNVs on gene expression in the central (CeA) and basolateral (BLA) amygdala, paraventricular nucleus (PVN), and cingulate cortex (Cg) was shown by a two-proportion Z-test (p = 1.6 x 10-31), with a positive correlation in the CeA (p = 0.0062), PVN (p = 0.0046) and Cg (p = 0.0114), indicating a contribution of CNVs to the genetic predisposition to trait anxiety in the specific context of HAB/LAB mice. In order to confirm anxiety-relevant CNVs and corresponding genes in a second mouse model, we further examined CD-1 outbred mice. We revealed the distribution of CNVs by genotyping 64 CD 1 individuals using a high-density genotyping array (Jackson Laboratory). 78 genes within those CNVs were identified to show nominally significant association (48 genes), or a statistical trend in their association (30 genes) with the time animals spent on the open arms of the elevated plus-maze (EPM). Fifteen of them were considered promising candidate genes of anxiety-related behavior as we could show a significant overlap (permutation test, p = 0.0051) with genes within HAB/LAB CNVs. Thus, here we provide what is to our knowledge the first extensive catalogue of CNVs in CD-1 mice and potential corresponding candidate genes linked to anxiety-related behavior in mice.

Using the MCF10A/MCF10CA1a Breast Cancer Progression Cell Line Model to Investigate the Effect of Active, Mutant Forms of EGFR in Breast Cancer Development and Treatment Using Gefitinib.

BACKGROUND: Basal-like and triple negative breast cancer (TNBC) share common molecular features, poor prognosis and a propensity for metastasis to the brain. Amplification of epidermal growth factor receptor (EGFR) occurs in ~50% of basal-like breast cancer, and mutations in the epidermal growth factor receptor (EGFR) have been reported in up to ~ 10% of Asian TNBC patients. In non-small cell lung cancer several different mutations in the EGFR tyrosine kinase domain confer sensitivity to receptor tyrosine kinase inhibitors, but the tumourigenic potential of EGFR mutations in breast cells and their potential for targeted therapy is unknown. MATERIALS AND METHODS: Constructs containing wild type, G719S or E746-A750 deletion mutant forms of EGFR were transfected into the MCF10A breast cells and their tumorigenic derivative, MCF10CA1a. The effects of EGFR over-expression and mutation on proliferation, migration, invasion, response to gefitinib, and tumour formation in vivo was investigated. Copy number analysis and whole exome sequencing of the MCF10A and MCF10CA1a cell lines were also performed. RESULTS: Mutant EGFR increased MCF10A and MCF10CA1a proliferation and MCF10A gefitinib sensitivity. The EGFR-E746-A750 deletion increased MCF10CA1a cell migration and invasion, and greatly increased MCF10CA1a xenograft tumour formation and growth. Compared to MCF10A cells, MCF10CA1a cells exhibited large regions of gain on chromosomes 3 and 9, deletion on chromosome 7, and mutations in many genes implicated in cancer. CONCLUSIONS: Mutant EGFR enhances the oncogenic properties of MCF10A cell line, and increases sensitivity to gefitinib. Although the addition of EGFR E746-A750 renders the MCF10CA1a cells more tumourigenic in vivo it is not accompanied by increased gefitinib sensitivity, perhaps due to additional mutations, including the PIK3CA H1047R mutation, that the MCF10CA1a cell line has acquired. Screening TNBC/basal-like breast cancer for EGFR mutations may prove useful for directing therapy but, as in non-small cell lung cancer, accompanying mutations in PIK3CA may confer gefitinib resistance.