Next-generation phenotyping integrated in a national framework for patients with ultrarare disorders improves genetic diagnostics and yields new molecular findings.

Individuals with ultrarare disorders pose a structural challenge for healthcare systems since expert clinical knowledge is required to establish diagnoses. In TRANSLATE NAMSE, a 3-year prospective study, we evaluated a novel diagnostic concept based on multidisciplinary expertise in Germany. Here we present the systematic investigation of the phenotypic and molecular genetic data of 1,577 patients who had undergone exome sequencing and were partially analyzed with next-generation phenotyping approaches. Molecular genetic diagnoses were established in 32% of the patients totaling 370 distinct molecular genetic causes, most with prevalence below 1:50,000. During the diagnostic process, 34 novel and 23 candidate genotype-phenotype associations were identified, mainly in individuals with neurodevelopmental disorders. Sequencing data of the subcohort that consented to computer-assisted analysis of their facial images with GestaltMatcher could be prioritized more efficiently compared with approaches based solely on clinical features and molecular scores. Our study demonstrates the synergy of using next-generation sequencing and phenotyping for diagnosing ultrarare diseases in routine healthcare and discovering novel etiologies by multidisciplinary teams.

Identification of de novo variants in nonsyndromic cleft lip with/without cleft palate patients with low polygenic risk scores.

BACKGROUND: Nonsyndromic cleft lip with/without cleft palate (nsCL/P) is a congenital malformation of multifactorial etiology. Research has identified >40 genome-wide significant risk loci, which explain less than 40% of nsCL/P heritability. Studies show that some of the hidden heritability is explained by rare penetrant variants. METHODS: To identify new candidate genes, we searched for highly penetrant de novo variants (DNVs) in 50 nsCL/P patient/parent-trios with a low polygenic risk for the phenotype (discovery). We prioritized DNV-carrying candidate genes from the discovery for resequencing in independent cohorts of 1010 nsCL/P patients of diverse ethnicities and 1574 population-matched controls (replication). Segregation analyses and rare variant association in the replication cohort, in combination with additional data (genome-wide association data, expression, protein-protein-interactions), were used for final prioritization. CONCLUSION: In the discovery step, 60 DNVs were identified in 60 genes, including a variant in the established nsCL/P risk gene CDH1. Re-sequencing of 32 prioritized genes led to the identification of 373 rare, likely pathogenic variants. Finally, MDN1 and PAXIP1 were prioritized as top candidates. Our findings demonstrate that DNV detection, including polygenic risk score analysis, is a powerful tool for identifying nsCL/P candidate genes, which can also be applied to other multifactorial congenital malformations.

A Resilience Related Glial-Neurovascular Network Is Transcriptionally Activated after Chronic Social Defeat in Male Mice.

Upon chronic stress, a fraction of individuals shows stress resilience, which can prevent long-term mental dysfunction. The underlying molecular mechanisms are complex and have not yet been fully understood. In this study, we performed a data-driven behavioural stratification together with single-cell transcriptomics of the hippocampus in a mouse model of chronic social defeat stress. Our work revealed that in a sub-group exhibiting molecular responses upon chronic stress, the dorsal hippocampus is particularly involved in neuroimmune responses, angiogenesis, myelination, and neurogenesis, thereby enabling brain restoration and homeostasis after chronic stress. Based on these molecular insights, we applied rapamycin after the stress as a proof-of-concept pharmacological intervention and were able to substantially increase stress resilience. Our findings serve as a data resource and can open new avenues for further understanding of molecular processes underlying stress response and for targeted interventions supporting resilience.

Multi-omics signatures of alcohol use disorder in the dorsal and ventral striatum.

Alcohol Use Disorder (AUD) is a major contributor to global mortality and morbidity. Postmortem human brain tissue enables the investigation of molecular mechanisms of AUD in the neurocircuitry of addiction. We aimed to identify differentially expressed (DE) genes in the ventral and dorsal striatum between individuals with AUD and controls, and to integrate the results with findings from genome- and epigenome-wide association studies (GWAS/EWAS) to identify functionally relevant molecular mechanisms of AUD. DNA-methylation and gene expression (RNA-seq) data was generated from postmortem brain samples of 48 individuals with AUD and 51 controls from the ventral striatum (VS) and the dorsal striatal regions caudate nucleus (CN) and putamen (PUT). We identified DE genes using DESeq2, performed gene-set enrichment analysis (GSEA), and tested enrichment of DE genes in results of GWASs using MAGMA. Weighted correlation network analysis (WGCNA) was performed for DNA-methylation and gene expression data and gene overlap was tested. Differential gene expression was observed in the dorsal (FDR < 0.05), but not the ventral striatum of AUD cases. In the VS, DE genes at FDR < 0.25 were overrepresented in a recent GWAS of problematic alcohol use. The ARHGEF15 gene was upregulated in all three brain regions. GSEA in CN and VS pointed towards cell-structure associated GO-terms and in PUT towards immune pathways. The WGCNA modules most strongly associated with AUD showed strong enrichment for immune response and inflammation pathways. Our integrated analysis of multi-omics data sets provides further evidence for the importance of immune- and inflammation-related processes in AUD.

Sensitivity and specificity of loss of heterozygosity analysis for the classification of rare germline variants in BRCA1/2: results of the observational AGO-TR1 study (NCT02222883).

Variant-specific loss of heterozygosity (LOH) analyses may be useful to classify BRCA1/2 germline variants of unknown significance (VUS). The sensitivity and specificity of this approach, however, remains unknown. We performed comparative next-generation sequencing analyses of the BRCA1/2 genes using blood-derived and tumour-derived DNA of 488 patients with ovarian cancer enrolled in the observational AGO-TR1 trial (NCT02222883). Overall, 94 pathogenic, 90 benign and 24 VUS were identified in the germline. A significantly increased variant fraction (VF) of a germline variant in the tumour indicates loss of the wild-type allele; a decreased VF indicates loss of the variant allele. We demonstrate that significantly increased VFs predict pathogenicity with high sensitivity (0.84, 95% CI 0.77 to 0.91), poor specificity (0.63, 95% CI 0.53 to 0.73) and poor positive predictive value (PPV; 0.71, 95% CI 0.62 to 0.79). Significantly decreased VFs predict benignity with low sensitivity (0.26, 95% CI 0.17 to 0.35), high specificity (1.0, 95% CI 0.96 to 1.00) and PPV (1.0, 95% CI 0.85 to 1.00). Variant classification based on significantly increased VFs results in an unacceptable proportion of false-positive results. A significantly decreased VF in the tumour may be exploited as a reliable predictor for benignity, with no false-negative result observed. When applying the latter approach, VUS identified in four patients can now be considered benign. Trial registration number NCT02222883.

F8 Inversions at Xq28 Causing Hemophilia A Are Associated With Specific Methylation Changes: Implication for Molecular Epigenetic Diagnosis.

Diverse DNA structural variations (SVs) in human cancers and several other diseases are well documented. For genomic inversions in particular, the disease causing mechanism may not be clear, especially if the inversion border does not cross a coding sequence. Understanding about the molecular processes of these inverted genomic sequences, in a mainly epigenetic context, may provide additional information regarding sequence-specific regulation of gene expression in human diseases. Herein, we study one such inversion hotspot at Xq28, which leads to the disruption of F8 gene and results in hemophilia A phenotype. To determine the epigenetic consequence of this rearrangement, we evaluated DNA methylation levels of 12 CpG rich regions with the coverage of 550 kb by using bisulfite-pyrosequencing and next-generation sequencing (NGS)-based bisulfite re-sequencing enrichment assay. Our results show that this inversion prone area harbors widespread methylation changes at the studied regions. However, only 5/12 regions showed significant methylation changes, specifically in case of intron 1 inversion (two regions), intron 22 inversion (two regions) and one common region in both inversions. Interestingly, these aberrant methylated regions were found to be overlapping with the inversion proximities. In addition, two CpG sites reached 100% sensitivity and specificity to discriminate wild type from intron 22 and intron 1 inversion samples. While we found age to be an influencing factor on methylation levels at some regions, covariate analysis still confirms the differential methylation induced by inversion, regardless of age. The hemophilia A methylation inversion "HAMI" assay provides an advantage over conventional PCR-based methods, which may not detect novel rare genomic rearrangements. Taken together, we showed that genomic inversions in the F8 (Xq28) region are associated with detectable changes in methylation levels and can be used as an epigenetic diagnostic marker.

Deleterious somatic variants in 473 consecutive individuals with ovarian cancer: results of the observational AGO-TR1 study (NCT02222883).

BACKGROUND: For individuals with ovarian cancer (OC), therapy options mainly depend on BRCA1/2 germline status. What is the prevalence of deleterious somatic variants, that is, does genetic tumour testing identify subgroups of individuals who also might benefit from targeted therapy? METHODS: Paired analysis of tumour-derived versus blood-derived DNA to determine the prevalence of deleterious somatic variants in OC predisposition genes (ATM, BRCA1/2, BRIP1, MSH2/6, PALB2, RAD51C/D and TP53) and the PIK3CA and PTEN genes in individuals with OC (AGO-TR1 study, NCT02222883). Results were complemented by BRCA1, PALB2 and RAD51C promoter methylation analyses and stratified by histological subtype; 473 individuals were included. RESULTS: The combined analyses revealed that deleterious germline variants in established OC predisposition genes (all: 125/473, 26.4%; BRCA1/2: 97/473, 20.5%), deleterious somatic variants in established OC predisposition genes excluding TP53 (all: 39/473, 8.2%; BRCA1/2: 30/473, 6.3%) and promoter methylation (all: 67/473, 14.2%; BRCA1: 57/473, 12.1%; RAD51C: 10/473, 2.1%; PALB2: 0/473) were mutually exclusive, with a few exceptions. The same holds true for deleterious somatic PIK3CA and/or PTEN variants (33/473, 7.0%) found to be enriched in endometrioid and clear cell OC (16/35, 45.7%); 84.3 % of the deleterious single-nucleotide/indel germline variants in established OC predisposition genes showed significantly higher variant fractions (VFs) in the tumour-derived versus blood-derived DNA, indicating a loss of the wild-type alleles. CONCLUSION: Tumour sequencing of the BRCA1, BRCA2, PIK3CA and PTEN genes along with BRCA1 and RAD51C promoter methylation analyses identified large subgroups of germline mutation-negative individuals who may be addressed in interventional studies using PARP or PI3K/AKT/mTOR inhibitors. TRIAL REGISTRATION NUMBER: NCT02222883.

Compound heterozygous SZT2 mutations in two siblings with early-onset epilepsy, intellectual disability and macrocephaly.

PURPOSE: Mutations in SZT2 have been previously reported in several cases of early onset epilepsy and intellectual disability. In this study we investigate potential causal mutations in two male siblings affected by early onset epilepsy, intellectual disability and macrocephaly. METHODS: We use family-based whole-exome sequencing to identify candidate variants. RESULTS: We report the identification of two potential causal SZT2 mutations in compound heterozygous state. We observe considerable differences in the clinical phenotype severity of the two affected individuals. The cerebral MRI revealed no abnormalities in the older affected brother, while in the youngest one it revealed a right frontal polymicrogiria. Moreover, while good seizure control was achieved in the older affected individual the younger brother is affected by pharmacoresistant epilepsy, progressive spastic paraplegia, cortical myoclonus and a more severe intellectual disability. We also analyzed the relative location of the reported pathogenic mutations in the SZT2 protein. CONCLUSION: Variable phenotypic expressivity is observed for this condition, while the location and type of mutations in SZT2 also has a potential impact on epilepsy severity. These findings extend our knowledge of epileptogenic conditions related to SZT2 and mTOR signaling.

The severity of human peri-implantitis lesions correlates with the level of submucosal microbial dysbiosis.

AIM: To cross-sectionally analyse the submucosal microbiome of peri-implantitis (PI) lesions at different severity levels. MATERIALS AND METHODS: Microbial signatures of 45 submucosal plaque samples from untreated PI lesions obtained from 30 non-smoking, systemically healthy subjects were assessed by 16s sequencing. Linear mixed models were used to identify taxa with differential abundance by probing depth, after correction for age, gender, and multiple samples per subject. Network analyses were performed to identify groups of taxa with mutual occurrence or exclusion. Subsequently, the effects of peri-implant probing depth on submucosal microbial dysbiosis were calculated using the microbial dysbiosis index. RESULTS: In total, we identified 337 different taxa in the submucosal microbiome of PI. Total abundance of 12 taxa correlated significantly with increasing probing depth; a significant relationship with lower probing depth was found for 16 taxa. Network analysis identified two mutually exclusive complexes associated with shallow pockets and deeper pockets, respectively. Deeper peri-implant pockets were associated with significantly increased dysbiosis. CONCLUSION: Increases in peri-implant pocket depth are associated with substantial changes in the submucosal microbiome and increasing levels of dysbiosis.

Prevalence of deleterious germline variants in risk genes including BRCA1/2 in consecutive ovarian cancer patients (AGO-TR-1).

BACKGROUND: Identification of families at risk for ovarian cancer offers the opportunity to consider prophylactic surgery thus reducing ovarian cancer mortality. So far, identification of potentially affected families in Germany was solely performed via family history and numbers of affected family members with breast or ovarian cancer. However, neither the prevalence of deleterious variants in BRCA1/2 in ovarian cancer in Germany nor the reliability of family history as trigger for genetic counselling has ever been evaluated. METHODS: Prospective counseling and germline testing of consecutive patients with primary diagnosis or with platinum-sensitive relapse of an invasive epithelial ovarian cancer. Testing included 25 candidate and established risk genes. Among these 25 genes, 16 genes (ATM, BRCA1, BRCA2, CDH1, CHEK2, MLH1, MSH2, MSH6, NBN, PMS2, PTEN, PALB2, RAD51C, RAD51D, STK11, TP53) were defined as established cancer risk genes. A positive family history was defined as at least one relative with breast cancer or ovarian cancer or breast cancer in personal history. RESULTS: In total, we analyzed 523 patients: 281 patients with primary diagnosis of ovarian cancer and 242 patients with relapsed disease. Median age at primary diagnosis was 58 years (range 16-93) and 406 patients (77.6%) had a high-grade serous ovarian cancer. In total, 27.9% of the patients showed at least one deleterious variant in all 25 investigated genes and 26.4% in the defined 16 risk genes. Deleterious variants were most prevalent in the BRCA1 (15.5%), BRCA2 (5.5%), RAD51C (2.5%) and PALB2 (1.1%) genes. The prevalence of deleterious variants did not differ significantly between patients at primary diagnosis and relapse. The prevalence of deleterious variants in BRCA1/2 (and in all 16 risk genes) in patients <60 years was 30.2% (33.2%) versus 10.6% (18.9%) in patients ≥60 years. Family history was positive in 43% of all patients. Patients with a positive family history had a prevalence of deleterious variants of 31.6% (36.0%) versus 11.4% (17.6%) and histologic subtype of high grade serous ovarian cancer versus other showed a prevalence of deleterious variants of 23.2% (29.1%) and 10.2% (14.8%), respectively. Testing only for BRCA1/2 would miss in our series more than 5% of the patients with a deleterious variant in established risk genes. CONCLUSIONS: 26.4% of all patients harbor at least one deleterious variant in established risk genes. The threshold of 10% mutation rate which is accepted for reimbursement by health care providers in Germany was observed in all subgroups analyzed and neither age at primary diagnosis nor histo-type or family history sufficiently enough could identify a subgroup not eligible for genetic counselling and testing. Genetic testing should therefore be offered to every patient with invasive epithelial ovarian cancer and limiting testing to BRCA1/2 seems to be not sufficient.

CHD1L: a new candidate gene for congenital anomalies of the kidneys and urinary tract (CAKUT).

BACKGROUND: Recently, we identified a microduplication in chromosomal band 1q21.1 encompassing the CHD1L/ALC1 gene encoding a chromatin-remodelling enzyme in congenital anomalies of the kidneys and urinary tract (CAKUT) patient. METHODS: To explore the role of CHD1L in CAKUT, we screened 85 CAKUT patients for mutations in the CHD1L gene and performed functional analyses of the three heterozygous missense variants detected. In addition, we quantitatively determined CHD1L expression in multiple human fetal and adult tissues and analysed expression of CHD1L protein in human embryonal, adult and hydronephrotic kidney sections. RESULTS: Two of three novel heterozygous missense variants identified in three patients were not found in >400 control chromosomes. All variants lead to amino acid substitutions in or near the CHD1L macro domain, a poly-ADP-ribose (PAR)-binding module interacting with PAR polymerase 1 (PARP1), and showed decreased interaction with PARP1 by pull-down assay of transfected cell lysates. Quantitative messenger RNA analysis demonstrated high CHD1L expression in human fetal kidneys, and levels were four times higher than in adult kidneys. In the human embryo at 7-11 weeks gestation, CHD1L immunolocalized in the early ureteric bud and the S- and comma-shaped bodies, critical stages of kidney development. In normal postnatal sections, CHD1L was expressed in the cytoplasm of tubular cells in all tubule segments. CHD1L expression appeared higher in the hydronephrotic kidney of one patient with a hypofunctional CHD1L variant than in normal kidneys, recapitulating high fetal levels. CONCLUSION: Our data suggest that CHD1L plays a role in kidney development and may be a new candidate gene for CAKUT.

Pathogenic effects of novel mutations in the P-type ATPase ATP13A2 (PARK9) causing Kufor-Rakeb syndrome, a form of early-onset parkinsonism.

Kufor-Rakeb syndrome (KRS) is a rare form of autosomal recessive juvenile or early-onset, levodopa responsive parkinsonism and has been associated with mutations in ATP13A2(also known as PARK9), a lysosomal type 5 P-type ATPase. Recently, we identified novel compound heterozygous mutations, c.3176T>G (p.L1059R) and c.3253delC (p.L1085WfsX1088) in ATP13A2 of two siblings affected with KRS. When overexpressed, wild-type ATP13A2 localized to Lysotracker-positive and LAMP2-positive lysosomes while both truncating and missense mutated ATP13A2 were retained in the endoplasmic reticulum (ER). Both mutant proteins were degraded by the proteasomal but not the lysosomal pathways. In addition, ATP13A2 mRNA with c.3253delC was degraded by nonsense-mediated mRNA decay (NMD), which was protected by cycloheximide treatment. To validate our findings in a biologically relevant setting, we used patient-derived human olfactory neurosphere cultures and fibroblasts and demonstrated persistent ER stress by detecting upregulation of unfolded protein response-related genes in the patient-derived cells. We also confirmed NMD degraded ATP13A2 c.3253delC mRNA in the cells. These findings indicate that these novel ATP13A2 mutations are indeed pathogenic and support the notion that mislocalization of the mutant ATP13A2, resultant ER stress, alterations in the proteasomal pathways and premature degradation of mutant ATP13A2 mRNA contribute to the aetiology of KRS.

Hereditary parkinsonism with dementia is caused by mutations in ATP13A2, encoding a lysosomal type 5 P-type ATPase.

Neurodegenerative disorders such as Parkinson and Alzheimer disease cause motor and cognitive dysfunction and belong to a heterogeneous group of common and disabling disorders. Although the complex molecular pathophysiology of neurodegeneration is largely unknown, major advances have been achieved by elucidating the genetic defects underlying mendelian forms of these diseases. This has led to the discovery of common pathophysiological pathways such as enhanced oxidative stress, protein misfolding and aggregation and dysfunction of the ubiquitin-proteasome system. Here, we describe loss-of-function mutations in a previously uncharacterized, predominantly neuronal P-type ATPase gene, ATP13A2, underlying an autosomal recessive form of early-onset parkinsonism with pyramidal degeneration and dementia (PARK9, Kufor-Rakeb syndrome). Whereas the wild-type protein was located in the lysosome of transiently transfected cells, the unstable truncated mutants were retained in the endoplasmic reticulum and degraded by the proteasome. Our findings link a class of proteins with unknown function and substrate specificity to the protein networks implicated in neurodegeneration and parkinsonism.