Unravelling undiagnosed rare disease cases by HiFi long-read genome sequencing.

Solve-RD is a pan-European rare disease (RD) research program that aims to identify disease-causing genetic variants in previously undiagnosed RD families. We utilised 10-fold coverage HiFi long-read sequencing (LRS) for detecting causative structural variants (SVs), single nucleotide variants (SNVs), insertion-deletions (InDels), and short tandem repeat (STR) expansions in extensively studied RD families without clear molecular diagnoses. Our cohort includes 293 individuals from 114 genetically undiagnosed RD families selected by European Rare Disease Network (ERN) experts. Of these, 21 families were affected by so-called 'unsolvable' syndromes for which genetic causes remain unknown, and 93 families with at least one individual affected by a rare neurological, neuromuscular, or epilepsy disorder without genetic diagnosis despite extensive prior testing. Clinical interpretation and orthogonal validation of variants in known disease genes yielded thirteen novel genetic diagnoses due to de novo and rare inherited SNVs, InDels, SVs, and STR expansions. In an additional four families, we identified a candidate disease-causing SV affecting several genes including an MCF2 / FGF13 fusion and PSMA3 deletion. However, no common genetic cause was identified in any of the 'unsolvable' syndromes. Taken together, we found (likely) disease-causing genetic variants in 13.0% of previously unsolved families and additional candidate disease-causing SVs in another 4.3% of these families. In conclusion, our results demonstrate the added value of HiFi long-read genome sequencing in undiagnosed rare diseases.

Systematic analysis of paralogous regions in 41,755 exomes uncovers clinically relevant variation.

The short lengths of short-read sequencing reads challenge the analysis of paralogous genomic regions in exome and genome sequencing data. Most genetic variants within these homologous regions therefore remain unidentified in standard analyses. Here, we present a method (Chameleolyser) that accurately identifies single nucleotide variants and small insertions/deletions (SNVs/Indels), copy number variants and ectopic gene conversion events in duplicated genomic regions using whole-exome sequencing data. Application to a cohort of 41,755 exome samples yields 20,432 rare homozygous deletions and 2,529,791 rare SNVs/Indels, of which we show that 338,084 are due to gene conversion events. None of the SNVs/Indels are detectable using regular analysis techniques. Validation by high-fidelity long-read sequencing in 20 samples confirms >88% of called variants. Focusing on variation in known disease genes leads to a direct molecular diagnosis in 25 previously undiagnosed patients. Our method can readily be applied to existing exome data.

Evolution of age-related mutation-driven clonal haematopoiesis over 20 years is associated with metabolic dysfunction in obesity.

BACKGROUND: Haematopoietic clones caused by somatic mutations with ≥2% variant allele frequency (VAF) increase with age and are linked to risk of haematological malignancies and cardiovascular disease. Recent observations suggest that smaller clones (VAF<2%) are also associated with adverse outcomes. Our aims were to determine the prevalence of clonal haematopoiesis driven by clones of variable sizes in individuals with obesity treated by usual care or bariatric surgery (a treatment that improves metabolic status), and to examine the expansion of clones in relation to age and metabolic dysregulation over up to 20 years. METHODS: Clonal haematopoiesis-driver mutations (CHDMs) were identified in blood samples from participants of the Swedish Obese Subjects intervention study. Using an ultrasensitive assay, we analysed single-timepoint samples from 1050 individuals treated by usual care and 841 individuals who had undergone bariatric surgery, and multiple-timepoint samples taken over 20 years from a subset (n = 40) of the individuals treated by usual care. FINDINGS: In this explorative study, prevalence of CHDMs was similar in the single-timepoint usual care and bariatric surgery groups (20.6% and 22.5%, respectively, P = 0.330), with VAF ranging from 0.01% to 31.15%. Clone sizes increased with age in individuals with obesity, but not in those who underwent bariatric surgery. In the multiple-timepoint analysis, VAF increased by on average 7% (range -4% to 24%) per year and rate of clone growth was negatively associated with HDL-cholesterol (R = -0.68, 1.74 E-04). INTERPRETATION: Low HDL-C was associated with growth of haematopoietic clones in individuals with obesity treated by usual care. FUNDING: The Swedish Research Council, The Swedish state under an agreement between the Swedish government and the county councils, the ALF (Avtal om Läkarutbildning och Forskning) agreement, The Swedish Heart-Lung Foundation, The Novo Nordisk Foundation, The European Research Council, The Netherlands Organisation for Scientific Research.

Clonal Hematopoiesis Is Associated With Low CD4 Nadir and Increased Residual HIV Transcriptional Activity in Virally Suppressed Individuals With HIV.

Clonal hematopoiesis, a common age-related phenomenon marked by expansion of cells with clonal hematopoiesis driver mutations, has been associated with all-cause mortality, cancer, and cardiovascular disease. People with HIV (PWH) are at risk for non-AIDS-related comorbidities such as atherosclerotic cardiovascular disease and cancer. In a cross-sectional cohort study, we compared clonal hematopoiesis prevalence in PWH on stable antiretroviral therapy with prevalence in a cohort of overweight individuals and a cohort of age- and sex-matched population controls. The prevalence of clonal hematopoiesis adjusted for age was increased and clone size was larger in PWH compared to population controls. Clonal hematopoiesis is associated with low CD4 nadir, increased residual HIV-1 transcriptional activity, and coagulation factors in PWH. Future studies on the effect of clonal hematopoiesis on the HIV reservoir and non-AIDS-related comorbidities are warranted.

Impact of rare and common genetic variation in the interleukin-1 pathway on human cytokine responses.

BACKGROUND: The interleukin (IL)-1 pathway is primarily associated with innate immunological defense and plays a major role in the induction and regulation of inflammation. Both common and rare genetic variation in this pathway underlies various inflammation-mediated diseases, but the role of rare variants relative to common variants in immune response variability in healthy individuals remains unclear. METHODS: We performed molecular inversion probe sequencing on 48 IL-1 pathway-related genes in 463 healthy individuals from the Human Functional Genomics Project. We functionally grouped common and rare variants, over gene, subpathway, and inflammatory levels and performed the Sequence Kernel Association Test to test for association with in vitro stimulation-induced cytokine responses; specifically, IL-1ß and IL-6 cytokine measurements upon stimulations that represent an array of microbial infections: lipopolysaccharide (LPS), phytohaemagglutinin (PHA), Candida albicans (C. albicans), and Staphylococcus aureus (S. aureus). RESULTS: We identified a burden of NCF4 rare variants with PHA-induced IL-6 cytokine and showed that the respective carriers are in the 1% lowest IL-6 producers. Collapsing rare variants in IL-1 subpathway genes produces a bidirectional association with LPS-induced IL-1ß cytokine levels, which is reflected by a significant Spearman correlation. On the inflammatory level, we identified a burden of rare variants in genes encoding for proteins with an anti-inflammatory function with S. aureus-induced IL-6 cytokine. In contrast to these rare variant findings which were based on different types of stimuli, common variant associations were exclusively identified with C. albicans-induced cytokine over various levels of grouping, from the gene, to subpathway, to inflammatory level. CONCLUSIONS: In conclusion, this study shows that functionally grouping common and rare genetic variants enables the elucidation IL-1-mediated biological mechanisms, specifically, for IL-1ß and IL-6 cytokine responses induced by various stimuli. The framework used in this study may allow for the analysis of rare and common genetic variants in a wider variety of (non-immune) complex phenotypes and therefore has the potential to contribute to better understanding of unresolved, complex traits and diseases.

Clinical Validation of Whole Genome Sequencing for Cancer Diagnostics.

Whole genome sequencing (WGS) using fresh-frozen tissue and matched blood samples from cancer patients may become the most complete genetic tumor test. With the increasing availability of small biopsies and the need to screen more number of biomarkers, the use of a single all-inclusive test is preferable over multiple consecutive assays. To meet high-quality diagnostics standards, we optimized and clinically validated WGS sample and data processing procedures, resulting in a technical success rate of 95.6% for fresh-frozen samples with sufficient (≥20%) tumor content. Independent validation of identified biomarkers against commonly used diagnostic assays showed a high sensitivity (recall; 98.5%) and precision (positive predictive value; 97.8%) for detection of somatic single-nucleotide variants and insertions and deletions (across 22 genes), and high concordance for detection of gene amplification (97.0%; EGFR and MET) as well as somatic complete loss (100%; CDKN2A/p16). Gene fusion analysis showed a concordance of 91.3% between DNA-based WGS and an orthogonal RNA-based gene fusion assay. Microsatellite (in)stability assessment showed a sensitivity of 100% with a precision of 94%, and virus detection (human papillomavirus), an accuracy of 100% compared with standard testing. In conclusion, whole genome sequencing has a >95% sensitivity and precision compared with routinely used DNA techniques in diagnostics, and all relevant mutation types can be detected reliably in a single assay.

A Genetics-First Approach Revealed Monogenic Disorders in Patients With ARM and VACTERL Anomalies.

Background: The VATER/VACTERL association (VACTERL) is defined as the non-random occurrence of the following congenital anomalies: Vertebral, Anal, Cardiac, Tracheal-Esophageal, Renal, and Limb anomalies. As no unequivocal candidate gene has been identified yet, patients are diagnosed phenotypically. The aims of this study were to identify patients with monogenic disorders using a genetics-first approach, and to study whether variants in candidate genes are involved in the etiology of VACTERL or the individual features of VACTERL: Anorectal malformation (ARM) or esophageal atresia with or without trachea-esophageal fistula (EA/TEF). Methods: Using molecular inversion probes, a candidate gene panel of 56 genes was sequenced in three patient groups: VACTERL (n = 211), ARM (n = 204), and EA/TEF (n = 95). Loss-of-function (LoF) and additional likely pathogenic missense variants, were prioritized and validated using Sanger sequencing. Validated variants were tested for segregation and patients were clinically re-evaluated. Results: In 7 out of the 510 patients (1.4%), pathogenic or likely pathogenic variants were identified in SALL1, SALL4, and MID1, genes that are associated with Townes-Brocks, Duane-radial-ray, and Opitz-G/BBB syndrome. These syndromes always include ARM or EA/TEF, in combination with at least two other VACTERL features. We did not identify LoF variants in the remaining candidate genes. Conclusions: None of the other candidate genes were identified as novel unequivocal disease genes for VACTERL. However, a genetics-first approach allowed refinement of the clinical diagnosis in seven patients, in whom an alternative molecular-based diagnosis was found with important implications for the counseling of the families.

Presence of Genetic Variants Among Young Men With Severe COVID-19.

Importance: Severe coronavirus disease 2019 (COVID-19) can occur in younger, predominantly male, patients without preexisting medical conditions. Some individuals may have primary immunodeficiencies that predispose to severe infections caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Objective: To explore the presence of genetic variants associated with primary immunodeficiencies among young patients with COVID-19. Design, Setting, and Participants: Case series of pairs of brothers without medical history meeting the selection criteria of young (age <35 years) brother pairs admitted to the intensive care unit (ICU) due to severe COVID-19. Four men from 2 unrelated families were admitted to the ICUs of 4 hospitals in the Netherlands between March 23 and April 12, 2020. The final date of follow-up was May 16, 2020. Available family members were included for genetic variant segregation analysis and as controls for functional experiments. Exposure: Severe COVID-19. Main Outcome and Measures: Results of rapid clinical whole-exome sequencing, performed to identify a potential monogenic cause. Subsequently, basic genetic and immunological tests were performed in primary immune cells isolated from the patients and family members to characterize any immune defects. Results: The 4 male patients had a mean age of 26 years (range, 21-32), with no history of major chronic disease. They were previously well before developing respiratory insufficiency due to severe COVID-19, requiring mechanical ventilation in the ICU. The mean duration of ventilatory support was 10 days (range, 9-11); the mean duration of ICU stay was 13 days (range, 10-16). One patient died. Rapid clinical whole-exome sequencing of the patients and segregation in available family members identified loss-of-function variants of the X-chromosomal TLR7. In members of family 1, a maternally inherited 4-nucleotide deletion was identified (c.2129_2132del; p.[Gln710Argfs*18]); the affected members of family 2 carried a missense variant (c.2383G>T; p.[Val795Phe]). In primary peripheral blood mononuclear cells from the patients, downstream type I interferon (IFN) signaling was transcriptionally downregulated, as measured by significantly decreased mRNA expression of IRF7, IFNB1, and ISG15 on stimulation with the TLR7 agonist imiquimod as compared with family members and controls. The production of IFN-γ, a type II IFN, was decreased in patients in response to stimulation with imiquimod. Conclusions and Relevance: In this case series of 4 young male patients with severe COVID-19, rare putative loss-of-function variants of X-chromosomal TLR7 were identified that were associated with impaired type I and II IFN responses. These preliminary findings provide insights into the pathogenesis of COVID-19.

Rare genetic variants in interleukin-37 link this anti-inflammatory cytokine to the pathogenesis and treatment of gout.

OBJECTIVE: Gout is characterised by severe interleukin (IL)-1-mediated joint inflammation induced by monosodium urate crystals. Since IL-37 is a pivotal anti-inflammatory cytokine suppressing the activity of IL-1, we conducted genetic and functional studies aimed at elucidating the role of IL-37 in the pathogenesis and treatment of gout. METHODS: Variant identification was performed by DNA sequencing of all coding bases of IL37 using molecular inversion probe-based resequencing (discovery cohort: gout n=675, controls n=520) and TaqMan genotyping (validation cohort: gout n=2202, controls n=2295). Predictive modelling of the effects of rare variants on protein structure was followed by in vitro experiments evaluating the impact on protein function. Treatment with recombinant IL-37 was evaluated in vitro and in vivo in a mouse model of gout. RESULTS: We identified four rare variants in IL37 in six of the discovery gout patients; p.(A144P), p.(G174Dfs*16), p.(C181*) and p.(N182S), whereas none emerged in healthy controls (Fisher's exact p-value=0.043). All variants clustered in the functional domain of IL-37 in exon 5 (p-value=5.71×10-5). Predictive modelling and functional studies confirmed loss of anti-inflammatory functions and we substantiated the therapeutic potential of recombinant IL-37 in the treatment of gouty inflammation. Furthermore, the carrier status of p.(N182S)(rs752113534) was associated with increased risk (OR=1.81, p-value=0.031) of developing gout in hyperuricaemic individuals of Polynesian ancestry. CONCLUSION: Here, we provide genetic as well as mechanistic evidence for the role of IL-37 in the pathogenesis of gout, and highlight the therapeutic potential of recombinant IL-37 for the treatment of gouty arthritis.

Exome chip association study excluded the involvement of rare coding variants with large effect sizes in the etiology of anorectal malformations.

INTRODUCTION: Anorectal malformations (ARM) are rare congenital malformations, resulting from disturbed hindgut development. A genetic etiology has been suggested, but evidence for the involvement of specific genes is scarce. We evaluated the contribution of rare and low-frequency coding variants in ARM etiology, assuming a multifactorial model. METHODS: We analyzed 568 Caucasian ARM patients and 1,860 population-based controls using the Illumina HumanExome Beadchip array, which contains >240,000 rare and low-frequency coding variants. GenomeStudio clustering and calling was followed by re-calling of 'no-calls' using zCall for patients and controls simultaneously. Single variant and gene-based analyses were performed to identify statistically significant associations, applying Bonferroni correction. Following an extra quality control step, candidate variants were selected for validation using Sanger sequencing. RESULTS: When we applied a MAF of ≥1.0%, no variants or genes showed statistically significant associations with ARM. Using a MAF cut-off at 0.4%, 13 variants initially reached statistical significance, but had to be discarded upon further inspection: ten variants represented calling errors of the software, while the minor alleles of the remaining three variants were not confirmed by Sanger sequencing. CONCLUSION: Our results show that rare and low-frequency coding variants with large effect sizes, present on the exome chip do not contribute to ARM etiology.

Deletions and loss-of-function variants in TP63 associated with orofacial clefting.

We aimed to identify novel deletions and variants of TP63 associated with orofacial clefting (OFC). Copy number variants were assessed in three OFC families using microarray analysis. Subsequently, we analyzed TP63 in a cohort of 1072 individuals affected with OFC and 706 population-based controls using molecular inversion probes (MIPs). We identified partial deletions of TP63 in individuals from three families affected with OFC. In the OFC cohort, we identified several TP63 variants predicting to cause loss-of-function alleles, including a frameshift variant c.569_576del (p.(Ala190Aspfs*5)) and a nonsense variant c.997C>T (p.(Gln333*)) that introduces a premature stop codon in the DNA-binding domain. In addition, we identified the first missense variants in the oligomerization domain c.1213G>A (p.(Val405Met)), which occurred in individuals with OFC. This variant was shown to abrogate oligomerization of mutant p63 protein into oligomeric complexes, and therefore likely represents a loss-of-function allele rather than a dominant-negative. All of these variants were inherited from an unaffected parent, suggesting reduced penetrance of such loss-of-function alleles. Our data indicate that loss-of-function alleles in TP63 can also give rise to OFC as the main phenotype. We have uncovered the dosage-dependent functions of p63, which were previously rejected.

Somatic variants in autosomal dominant genes are a rare cause of sporadic Alzheimer's disease.

INTRODUCTION: A minority of patients with sporadic early-onset Alzheimer's disease (AD) exhibit de novo germ line mutations in the autosomal dominant genes such as APP, PSEN1, or PSEN2. We hypothesized that negatively screened patients may harbor somatic variants in these genes. METHODS: We applied an ultrasensitive approach based on single-molecule molecular inversion probes followed by deep next generation sequencing of 11 genes to 100 brain and 355 blood samples from 445 sporadic patients with AD (>80% exhibited an early onset, <66 years). RESULTS: We identified and confirmed nine somatic variants (allele fractions: 0.2%-10.8%): two APP, five SORL1, one NCSTN, and one MARK4 variants by independent amplicon-based deep sequencing. DISCUSSION: Two of the SORL1 variant might have contributed to the disease, the two APP variants were interpreted as likely benign and the other variants remained of unknown significance. Somatic variants in the autosomal dominant AD genes may not be a common cause of sporadic AD, including early onset cases.

A genotype-first approach identifies an intellectual disability-overweight syndrome caused by PHIP haploinsufficiency.

Genotype-first combined with reverse phenotyping has shown to be a powerful tool in human genetics, especially in the era of next generation sequencing. This combines the identification of individuals with mutations in the same gene and linking these to consistent (endo)phenotypes to establish disease causality. We have performed a MIP (molecular inversion probe)-based targeted re-sequencing study in 3,275 individuals with intellectual disability (ID) to facilitate a genotype-first approach for 24 genes previously implicated in ID.Combining our data with data from a publicly available database, we confirmed 11 of these 24 genes to be relevant for ID. Amongst these, PHIP was shown to have an enrichment of disruptive mutations in the individuals with ID (5 out of 3,275). Through international collaboration, we identified a total of 23 individuals with PHIP mutations and elucidated the associated phenotype. Remarkably, all 23 individuals had developmental delay/ID and the majority were overweight or obese. Other features comprised behavioral problems (hyperactivity, aggression, features of autism and/or mood disorder) and dysmorphisms (full eyebrows and/or synophrys, upturned nose, large ears and tapering fingers). Interestingly, PHIP encodes two protein-isoforms, PHIP/DCAF14 and NDRP, each involved in neurodevelopmental processes, including E3 ubiquitination and neuronal differentiation. Detailed genotype-phenotype analysis points towards haploinsufficiency of PHIP/DCAF14, and not NDRP, as the underlying cause of the phenotype.Thus, we demonstrated the use of large scale re-sequencing by MIPs, followed by reverse phenotyping, as a constructive approach to verify candidate disease genes and identify novel syndromes, highlighted by PHIP haploinsufficiency causing an ID-overweight syndrome.

Ultra-sensitive Sequencing Identifies High Prevalence of Clonal Hematopoiesis-Associated Mutations throughout Adult Life.

Clonal hematopoiesis results from somatic mutations in hematopoietic stem cells, which give an advantage to mutant cells, driving their clonal expansion and potentially leading to leukemia. The acquisition of clonal hematopoiesis-driver mutations (CHDMs) occurs with normal aging and these mutations have been detected in more than 10% of individuals ≥65 years. We aimed to examine the prevalence and characteristics of CHDMs throughout adult life. We developed a targeted re-sequencing assay combining high-throughput with ultra-high sensitivity based on single-molecule molecular inversion probes (smMIPs). Using smMIPs, we screened more than 100 loci for CHDMs in more than 2,000 blood DNA samples from population controls between 20 and 69 years of age. Loci screened included 40 regions known to drive clonal hematopoiesis when mutated and 64 novel candidate loci. We identified 224 somatic mutations throughout our cohort, of which 216 were coding mutations in known driver genes (DNMT3A, JAK2, GNAS, TET2, and ASXL1), including 196 point mutations and 20 indels. Our assay's improved sensitivity allowed us to detect mutations with variant allele frequencies as low as 0.001. CHDMs were identified in more than 20% of individuals 60 to 69 years of age and in 3% of individuals 20 to 29 years of age, approximately double the previously reported prevalence despite screening a limited set of loci. Our findings support the occurrence of clonal hematopoiesis-associated mutations as a widespread mechanism linked with aging, suggesting that mosaicism as a result of clonal evolution of cells harboring somatic mutations is a universal mechanism occurring at all ages in healthy humans.

Quantification of differential gene expression by multiplexed targeted resequencing of cDNA.

Whole-transcriptome or RNA sequencing (RNA-Seq) is a powerful and versatile tool for functional analysis of different types of RNA molecules, but sample reagent and sequencing cost can be prohibitive for hypothesis-driven studies where the aim is to quantify differential expression of a limited number of genes. Here we present an approach for quantification of differential mRNA expression by targeted resequencing of complementary DNA using single-molecule molecular inversion probes (cDNA-smMIPs) that enable highly multiplexed resequencing of cDNA target regions of ∼100 nucleotides and counting of individual molecules. We show that accurate estimates of differential expression can be obtained from molecule counts for hundreds of smMIPs per reaction and that smMIPs are also suitable for quantification of relative gene expression and allele-specific expression. Compared with low-coverage RNA-Seq and a hybridization-based targeted RNA-Seq method, cDNA-smMIPs are a cost-effective high-throughput tool for hypothesis-driven expression analysis in large numbers of genes (10 to 500) and samples (hundreds to thousands).

Overlapping SETBP1 gain-of-function mutations in Schinzel-Giedion syndrome and hematologic malignancies.

Schinzel-Giedion syndrome (SGS) is a rare developmental disorder characterized by multiple malformations, severe neurological alterations and increased risk of malignancy. SGS is caused by de novo germline mutations clustering to a 12bp hotspot in exon 4 of SETBP1. Mutations in this hotspot disrupt a degron, a signal for the regulation of protein degradation, and lead to the accumulation of SETBP1 protein. Overlapping SETBP1 hotspot mutations have been observed recurrently as somatic events in leukemia. We collected clinical information of 47 SGS patients (including 26 novel cases) with germline SETBP1 mutations and of four individuals with a milder phenotype caused by de novo germline mutations adjacent to the SETBP1 hotspot. Different mutations within and around the SETBP1 hotspot have varying effects on SETBP1 stability and protein levels in vitro and in in silico modeling. Substitutions in SETBP1 residue I871 result in a weak increase in protein levels and mutations affecting this residue are significantly more frequent in SGS than in leukemia. On the other hand, substitutions in residue D868 lead to the largest increase in protein levels. Individuals with germline mutations affecting D868 have enhanced cell proliferation in vitro and higher incidence of cancer compared to patients with other germline SETBP1 mutations. Our findings substantiate that, despite their overlap, somatic SETBP1 mutations driving malignancy are more disruptive to the degron than germline SETBP1 mutations causing SGS. Additionally, this suggests that the functional threshold for the development of cancer driven by the disruption of the SETBP1 degron is higher than for the alteration in prenatal development in SGS. Drawing on previous studies of somatic SETBP1 mutations in leukemia, our results reveal a genotype-phenotype correlation in germline SETBP1 mutations spanning a molecular, cellular and clinical phenotype.

BRCA Testing by Single-Molecule Molecular Inversion Probes.

BACKGROUND: Despite advances in next generation DNA sequencing (NGS), NGS-based single gene tests for diagnostic purposes require improvements in terms of completeness, quality, speed, and cost. Single-molecule molecular inversion probes (smMIPs) are a technology with unrealized potential in the area of clinical genetic testing. In this proof-of-concept study, we selected 2 frequently requested gene tests, those for the breast cancer genes BRCA1 and BRCA2, and developed an automated work flow based on smMIPs. METHODS: The BRCA1 and BRCA2 smMIPs were validated using 166 human genomic DNA samples with known variant status. A generic automated work flow was built to perform smMIP-based enrichment and sequencing for BRCA1, BRCA2, and the checkpoint kinase 2 (CHEK2) c.1100del variant. RESULTS: Pathogenic and benign variants were analyzed in a subset of 152 previously BRCA-genotyped samples, yielding an analytical sensitivity and specificity of 100%. Following automation, blind analysis of 65 in-house samples and 267 Norwegian samples correctly identified all true-positive variants (>3000), with no false positives. Consequent to process optimization, turnaround times were reduced by 60% to currently 10-15 days. Copy number variants were detected with an analytical sensitivity of 100% and an analytical specificity of 88%. CONCLUSIONS: smMIP-based genetic testing enables automated and reliable analysis of the coding sequences of BRCA1 and BRCA2. The use of single-molecule tags, double-tiled targeted enrichment, and capturing and sequencing in duplo, in combination with automated library preparation and data analysis, results in a robust process and reduces routine turnaround times. Furthermore, smMIP-based copy number variation analysis could make independent copy number variation tools like multiplex ligation-dependent probes amplification dispensable.

Chromosomal abnormalities in hepatic cysts point to novel polycystic liver disease genes.

Autosomal dominant polycystic liver disease (ADPLD) is caused by variants in PRKCSH, SEC63, and LRP5, whereas autosomal dominant polycystic kidney disease is caused by variants in PKD1 and PKD2. Liver cyst development in these disorders is explained by somatic loss-of-heterozygosity (LOH) of the wild-type allele in the developing cyst. We hypothesize that we can use this mechanism to identify novel disease genes that reside in LOH regions. In this study, we aim to map abnormal genomic regions using high-density SNP microarrays to find novel PLD genes. We collected 46 cysts from 23 patients with polycystic or sporadic hepatic cysts, and analyzed DNA from those cysts using high-resolution microarray (n=24) or Sanger sequencing (n=22). We here focused on regions of homozygosity on the autosomes (>3.0 Mb) and large CNVs (>1.0 Mb). We found frequent LOH in PRKCSH (22/29) and PKD1/PKD2 (2/3) cysts of patients with known heterozygous germline variants in the respective genes. In the total cohort, 12/23 patients harbored abnormalities outside of familiar areas. In individual ADPLD cases, we identified germline events: a 2q13 complex rearrangement resulting in BUB1 haploinsufficiency, a 47XXX karyotype, chromosome 9q copy-number loss, and LOH on chromosome 3p. The latter region was overlapping with an LOH region identified in two other cysts. Unique germline and somatic abnormalities occur frequently in and outside of known genes underlying cysts. Each liver cyst has a unique genetic makeup. LOH driver gene BUB1 may imply germline causes of genetic instability in PLD.

Novel mutations in LRP6 highlight the role of WNT signaling in tooth agenesis.

PURPOSE: We aimed to identify a novel genetic cause of tooth agenesis (TA) and/or orofacial clefting (OFC) by combining whole-exome sequencing (WES) and targeted resequencing in a large cohort of TA and OFC patients. METHODS: WES was performed in two unrelated patients: one with severe TA and OFC and another with severe TA only. After deleterious mutations were identified in a gene encoding low-density lipoprotein receptor-related protein 6 (LRP6), all its exons were resequenced with molecular inversion probes in 67 patients with TA, 1,072 patients with OFC, and 706 controls. RESULTS: We identified a frameshift (c.4594delG, p.Cys1532fs) and a canonical splice-site mutation (c.3398-2A>C, p.?) in LRP6, respectively, in the patient with TA and OFC and in the patient with severe TA only. The targeted resequencing showed significant enrichment of unique LRP6 variants in TA patients but not in nonsyndromic OFC patients. Of the five variants in patients with TA, two affected the canonical splice site and three were missense variants; all variants segregated with the dominant phenotype, and in one case the missense mutation occurred de novo. CONCLUSION: Mutations in LRP6 cause TA in humans.Genet Med 18 11, 1158-1162.

DVL3 Alleles Resulting in a -1 Frameshift of the Last Exon Mediate Autosomal-Dominant Robinow Syndrome.

Robinow syndrome is a rare congenital disorder characterized by mesomelic limb shortening, genital hypoplasia, and distinctive facial features. Recent reports have identified, in individuals with dominant Robinow syndrome, a specific type of variant characterized by being uniformly located in the penultimate exon of DVL1 and resulting in a -1 frameshift allele with a premature termination codon that escapes nonsense-mediated decay. Here, we studied a cohort of individuals who had been clinically diagnosed with Robinow syndrome but who had not received a molecular diagnosis from variant studies of DVL1, WNT5A, and ROR2. Because of the uniform location of frameshift variants in DVL1-mediated Robinow syndrome and the functional redundancy of DVL1, DVL2, and DVL3, we elected to pursue direct Sanger sequencing of the penultimate exon of DVL1 and its paralogs DVL2 and DVL3 to search for potential disease-associated variants. Remarkably, targeted sequencing identified five unrelated individuals harboring heterozygous, de novo frameshift variants in DVL3, including two splice acceptor mutations and three 1 bp deletions. Similar to the variants observed in DVL1-mediated Robinow syndrome, all variants in DVL3 result in a -1 frameshift, indicating that these highly specific alterations might be a common cause of dominant Robinow syndrome. Here, we review the current knowledge of these peculiar variant alleles in DVL1- and DVL3-mediated Robinow syndrome and further elucidate the phenotypic features present in subjects with DVL1 and DVL3 frameshift mutations.